Intra-host Diversity Research Articles

Using phylogenetics to infer HIV-1 transmission direction between known transmission pairs

Inferring the transmission direction between linked individuals living with HIV provides unparalleled power to understand the epidemiology that determines transmission. Phylogenetic ancestral-state reconstruction approaches infer the transmission direction by identifying the individual in whom the most recent common ancestor of the virus populations originated. While these methods vary in accuracy, it is unclear why. To evaluate the performance of phylogenetic ancestral-state reconstruction to determine the transmission direction of HIV-1 infection, we inferred the transmission direction for 112 transmission pairs where transmission direction and detailed additional information were available. We then fit a statistical model to evaluate the extent to which epidemiological, sampling, genetic, and phylogenetic factors influenced the outcome of the inference. Finally, we repeated the analysis under real-life conditions with only routinely available data. We found that whether ancestral-state reconstruction correctly infers the transmission direction depends principally on the phylogeny's topology. For example, under real-life conditions, the probability of identifying the correct transmission direction increases from 32%-when a monophyletic-monophyletic or paraphyletic-polyphyletic tree topology is observed and when the tip closest to the root does not agree with the state at the root-to 93% when a paraphyletic-monophyletic topology is observed and when the tip closest to the root agrees with the root state. Our results suggest that documenting larger differences in relative intrahost diversity increases our confidence in the transmission direction inference of linked pairs for population-level studies of HIV. These findings provide a practical starting point to determine our confidence in transmission direction inference from ancestral-state reconstruction.

In Depth Viral Diversity Analysis in Atypical Neurological and Neonatal Chikungunya Infections in Rio de Janeiro, Brazil.

Chikungunya virus (CHIKV) is an arthropod-borne virus (arbovirus) transmitted by Aedes mosquitoes. The human infection usually manifests as a febrile and incapacitating arthritogenic illness, self-limiting and non-lethal. However, since 2013, CHIKV spreading through the tropics and to the Americas was accompanied by an increasing number of cases of atypical disease presentation, namely severe neuropathies and neonatal infection due to intrapartum vertical transmission. The pathophysiological mechanisms underlying these conditions have not been fully elucidated. However, arbovirus intrahost genetic diversity is thought to be linked to viral pathogenesis. To determine whether particular viral variants could be somehow associated, we analyzed the intrahost genetic diversity of CHIKV in three infected patients with neurological manifestations and three mothers infected during the intrapartum period, as well as their babies following vertical transmission. No statistically supported differences were observed for the genetic variability (nucleotide substitutions/gene length) along the genome between the groups. However, the newborn and cerebrospinal fluid samples (corresponding to virus passed through the placenta and/or the blood–brain barrier (BBB)) presented a different composition of their intrahost mutant ensembles compared to maternal or patient serum samples, even when concurrent. This finding could be consistent with the unidirectional virus transmission through these barriers, and the effect of selective bottlenecks during the transmission event. In addition, a higher proportion of defective variants (insertions/deletions and stop codons) was detected in the CSF and maternal samples and those were mainly distributed within the viral non-structural genes. Since defective viral genomes in RNA viruses are known to contribute to the outcome of acute viral infections and influence disease severity, their role in these atypical cases should be further investigated. Finally, with the in silico approach adopted, we detected no relevant non-conservative mutational pattern that could provide any hint of the pathophysiological mechanisms underlying these atypical cases. The present analysis represents a unique contribution to our understanding of the transmission events in these cases and generates hypotheses regarding underlying mechanisms, that can be explored further.

Clinical, microbiological, and genomic characteristics of clade-III Candida auris colonization and infection in southern California, 2019-2022.

Candida auris is an emerging fungal pathogen causing outbreaks in healthcare facilities. Five distinctive genomic clades exhibit clade-unique characteristics, highlighting the importance of real-time genomic surveillance and incorporating genotypic information to inform infection prevention practices and treatment algorithms. Both active and passive surveillance were used to screen hospitalized patients. C. auris polymerase chain reaction (PCR) assay on inguinal-axillary swabs was performed on high-risk patients upon admission. All clinical yeast isolates were identified to the species level. C. auris isolates were characterized by both phenotypic antifungal susceptibility tests and whole-genome sequencing. From late 2019 to early 2022, we identified 45 patients with C. auris. Most had a tracheostomy or were from a facility with a known outbreak. Moreover, 7 patients (15%) were only identified through passive surveillance. Also, 8 (18%) of the patients had a history of severe COVID-19. The overall mortality was 18%. Invasive C. auris infections were identified in 13 patients (29%), 9 (69%) of whom had bloodstream infections. Patients with invasive infection were more likely to have a central line. All C. auris isolates were resistant to fluconazole but susceptible to echinocandins. Genomic analysis showed that 1 dominant clade-III lineage is circulating in Los Angeles, with very limited intrahost and interhost genetic diversity. We have demonstrated that a robust C. auris surveillance program can be established using both active and passive surveillance, with multidisciplinary efforts involving the microbiology laboratory and the hospital epidemiology team. In Los Angeles County, C. auris strains are highly related and echinocandins should be used for empiric therapy.

Human cytomegalovirus strain diversity and dynamics reveal the donor lung as a major contributor after transplantation.

Mixed human cytomegalovirus (HCMV) strain infections are frequent in lung transplant recipients (LTRs). To date, the influence of the donor (D) and recipient (R) HCMV serostatus on intra-host HCMV strain composition and viral population dynamics after transplantation is only poorly understood. Here, we investigated ten pre-transplant lungs from HCMV-seropositive donors and 163 sequential HCMV-DNA-positive plasma and bronchoalveolar lavage samples from fifty LTRs with multiviremic episodes post-transplantation. The study cohort included D+R+ (38 per cent), D+R− (36 per cent), and D−R+ (26 per cent) patients. All samples were subjected to quantitative genotyping by short amplicon deep sequencing, and twenty-four of them were additionally PacBio long-read sequenced for genotype linkages. We find that D+R+ patients show a significantly elevated intra-host strain diversity compared to D+R− and D−R+ patients (P = 0.0089). Both D+ patient groups display significantly higher viral population dynamics than D− patients (P = 0.0061). Five out of ten pre-transplant donor lungs were HCMV DNA positive, whereof three multiple HCMV strains were detected, indicating that multi-strain transmission via lung transplantation is likely. Using long reads, we show that intra-host haplotypes can share distinctly linked genotypes, which limits overall intra-host diversity in mixed infections. Together, our findings demonstrate donor-derived strains as the main source of increased HCMV strain diversity and dynamics post-transplantation. These results foster strategies to mitigate the potential transmission of the donor strain reservoir to the allograft, such as ex vivo delivery of HCMV-selective immunotoxins prior to transplantation to reduce latent HCMV.

The Second Human Pegivirus, a Non-Pathogenic RNA Virus with Low Prevalence and Minimal Genetic Diversity.

The second human pegivirus (HPgV-2) is a virus discovered in the plasma of a hepatitis C virus (HCV)-infected patient in 2015 belonging to the pegiviruses of the family Flaviviridae. HPgV-2 has been proved to be epidemiologically associated with and structurally similar to HCV but unrelated to HCV disease and non-pathogenic, but its natural history and tissue tropism remain unclear. HPgV-2 is a unique RNA virus sharing the features of HCV and the first human pegivirus (HPgV-1 or GBV-C). Moreover, distinct from most RNA viruses such as HCV, HPgV-1 and human immunodeficiency virus (HIV), HPgV-2 exhibits much lower genomic diversity, with a high global sequence identity ranging from 93.5 to 97.5% and significantly lower intra-host variation than HCV. The mechanisms underlying the conservation of the HPgV-2 genome are not clear but may include efficient innate immune responses, low immune selection pressure and, possibly, the unique features of the viral RNA-dependent RNA polymerase (RdRP). In this review, we summarize the prevalence, pathogenicity and genetic diversity of HPgV-2 and discuss the possible reasons for the uniformity of its genome sequence, which should elucidate the implications of RNA virus fidelity for attenuated viral vaccines.

Clinical sequencing uncovers the genomic characteristics and mutation spectrum of the 2018 African swine fever virus in Guangdong, China.

African swine fever (ASF) outbreak have caused tremendous economic loss to the pig industry in China since its emergence in August 2018. Previous studies revealed that many published sequences are not suitable for detailed analyses due to the lack of data regarding quality parameters and methodology, and outdated annotations. Thus, high-quality genomes of highly pathogenic strains that can be used as references for early Chinese ASF outbreaks are still lacking, and little is known about the features of intra-host variants of ASF virus (ASFV). In this study, a full genome sequencing of clinical samples from the first ASF outbreak in Guangdong in 2018 was performed using MGI (MGI Tech Co., Ltd., Shenzhen, China) and Nanopore sequencing platforms, followed by Sanger sequencing to verify the variations. With 22 sequencing corrections, we obtained a high-quality genome of one of the earliest virulent isolates, GZ201801_2. After proofreading, we improved (add or modify) the annotations of this isolate using the whole genome alignment with Georgia 2007/1. Based on the complete genome sequence, we constructed the methylation profiles of early ASFV strains in China and predicted the potential 5mC and 6mA methylation sites, which are likely involved in metabolism, transcription, and replication. Additionally, the intra-host single nucleotide variant distribution and mutant allele frequency in the clinical samples of early strain were determined for the first time and found a strong preference for A and T substitution mutation, non-synonymous mutations, and mutations that resulted in amino acid substitutions into Lysine. In conclusion, this study provides a high-quality genome sequence, updated genome annotation, methylation profile, and mutation spectrum of early ASFV strains in China, thereby providing a reference basis for further studies on the evolution, transmission, and virulence of ASFV.

Remdesivir-induced emergence of SARS-CoV2 variants in patients with prolonged infection

We here investigate the impact of antiviral treatments such as remdesivir on intra-host genomic diversity and emergence of SARS-CoV2 variants in patients with a prolonged course of infection.Sequencing and variant analysis performed in 112 longitudinal respiratory samples from 14 SARS-CoV2-infected patients with severe disease progression show that major frequency variants do not generally arise during prolonged infection. However, remdesivir treatment can increase intra-host genomic diversity and result in the emergence of novel major variant species harboring fixed mutations. This is particularly evident in a patient with B-cell depletion who rapidly developed mutations in the RNA-dependent RNA polymerase gene following remdesivir treatment.Remdesivir treatment-associated emergence of novel variants is of great interest in light of current treatment guidelines for hospitalized patients suffering from severe SARS-CoV2 disease, as well as the potential use of remdesivir to preventively treat non-hospitalized patients at high risk for severe disease progression.

Improved SARS-CoV-2 sequencing surveillance allows the identification of new variants and signatures in infected patients

BackgroundGenomic surveillance of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the only approach to rapidly monitor and tackle emerging variants of concern (VOC) of the COVID-19 pandemic. Such scrutiny is crucial to limit the spread of VOC that might escape the immune protection conferred by vaccination strategies or previous virus exposure. It is also becoming clear now that efficient genomic surveillance would require monitoring of the host gene expression to identify prognostic biomarkers of treatment efficacy and disease progression. Here we propose an integrative workflow to both generate thousands of SARS-CoV-2 genome sequences per week and analyze host gene expression upon infection.MethodsIn this study we applied an integrated workflow for RNA extracted from nasal swabs to obtain in parallel the full genome of SARS-CoV-2 and transcriptome of host respiratory epithelium. The RNA extracted from each sample was reverse transcribed and the viral genome was specifically enriched through an amplicon-based approach. The very same RNA was then used for patient transcriptome analysis. Samples were collected in the Campania region, Italy, for viral genome sequencing. Patient transcriptome analysis was performed on about 700 samples divided into two cohorts of patients, depending on the viral variant detected (B.1 or delta).ResultsWe sequenced over 20,000 viral genomes since the beginning of the pandemic, producing the highest number of sequences in Italy. We thus reconstructed the pandemic dynamics in the regional territory from March 2020 to December 2021. In addition, we have matured and applied novel proof-of-principle approaches to prioritize possible gain-of-function mutations by leveraging patients’ metadata and isolated patient-specific signatures of SARS-CoV-2 infection. This allowed us to (i) identify three new viral variants that specifically originated in the Campania region, (ii) map SARS-CoV-2 intrahost variability during long-term infections and in one case identify an increase in the number of mutations in the viral genome, and (iii) identify host gene expression signatures correlated with viral load in upper respiratory ways.ConclusionIn conclusion, we have successfully generated an optimized and cost-effective strategy to monitor SARS-CoV-2 genetic variability, without the need of automation. Thus, our approach is suitable for any lab with a benchtop sequencer and a limited budget, allowing an integrated genomic surveillance on premises. Finally, we have also identified a gene expression signature defining SARS-CoV-2 infection in real-world patients’ upper respiratory ways.

ViReMaShiny: an interactive application for analysis of viral recombination data.

MotivationRecombination is an essential driver of virus evolution and adaption, giving rise to new chimeric viruses, structural variants, sub-genomic RNAs and defective RNAs. Next-generation sequencing (NGS) of virus samples, either from experimental or clinical settings, has revealed a complex distribution of recombination events that contributes to intrahost diversity. We and others have previously developed alignment tools to discover and map these diverse recombination events in NGS data. However, there is no standard for data visualization to contextualize events of interest, and downstream analysis often requires bespoke coding.ResultsWe present ViReMaShiny, a web-based application built using the R Shiny framework to allow interactive exploration and point-and-click visualization of viral recombination data provided in BED format generated by computational pipelines such as ViReMa (Viral-Recombination-Mapper).Availability and implementationThe application is hosted at https://routhlab.shinyapps.io/ViReMaShiny/ with associated documentation at https://jayeung12.github.io/. Code is available at https://github.com/routhlab/ViReMaShiny.Supplementary information Supplementary data are available at Bioinformatics online.

SARS-CoV-2 shedding dynamics and transmission in immunosuppressed patients

ABSTRACT Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants of concern have been emerging. However, knowledge of temporal and spatial dynamics of SARS-CoV-2 is limited. This study characterized SARS-CoV-2 evolution in immunosuppressed patients with long-term SARS-CoV-2 shedding for 73–250 days, without specific treatment. We conducted whole-genome sequencing of 27 serial samples, including 26 serial samples collected from various anatomic sites of two patients and the first positive sample from patient 2‘s mother. We analysed the intrahost temporal dynamics and genomic diversity of the viral population within different sample types. Intrahost variants emerging during infection showed diversity between individual hosts. Remarkably, N501Y, P681R, and E484K, key substitutions within spike protein, emerged in vivo during infection and became the dominant population. P681R, which had not yet been detected in the publicly available genome in Korea, appeared within patient 1 during infection. Mutually exclusive substitutions at residues R346 (R346S and R346I) and E484 (E484K and E484A) of spike protein and continuous turnover of these substitutions occurred. Unique genetic changes were observed in urine samples. A household transmission from patient 2 to his mother, at least 38 days after the diagnosis, was characterized. Viruses may differently mutate and adjust to the host selective pressure, which could enable the virus to replicate efficiently for fitness in each host. Intrahost variants could be candidate variants likely to spread to the population eventually. Our findings may provide new insights into the dynamics of SARS-CoV-2 in response to interactions between the virus and host.

Time Series Genomics of Pseudomonas aeruginosa Reveals the Emergence of a Hypermutator Phenotype and Within-Host Evolution in Clinical Inpatients.

ABSTRACTPseudomonas aeruginosa, a common opportunistic pathogen, is one of the leading etiological agents of nosocomial infections. Many previous studies have reported the nosocomial transmission and epidemiology of P. aeruginosa infections. However, longitudinal studies regarding the dynamics of P. aeruginosa colonization and infection in health care settings are limited. We obtained longitudinal samples from aged patients with prolonged intensive care unit (ICU) stays (~4 to 19 months). P. aeruginosa was isolated from 71 samples obtained from seven patients and characterized by whole-genome sequencing. The P. aeruginosa isolates were assigned to 10 clonal complexes, and turnover of main clones was observed in sequential sputum samples from two patients. By comparing intraclonal genomic diversities, we identified two clones that had significantly higher numbers of single nucleotide polymorphisms and variations in homopolymeric sequences than the other clones, indicating a hypermutator phenotype. These hypermutator clones were associated with mutations T147I/G521S and P27L in the MutL protein, and their mutation rates were estimated to be 3.20 × 10−5 and 6.59 × 10−5 per year per nucleotide, respectively. We also identified 24 recurrently mutated genes that exhibited intraclonal diversity in two or more clones. Notably, one recurrent mutation, S698F in FptA, was observed in four clones. These findings suggest that convergent microevolution and adaption of P. aeruginosa occur in long-term ICU patients.IMPORTANCE Pseudomonas aeruginosa is a predominant opportunistic pathogen that causes nosocomial infections. Inappropriate empirical therapy can lead to prolonged hospital stays and increased mortality. In our study of sequential P. aeruginosa isolates from inpatients, high intrahost diversity was observed, including switching of clones and the emergence of a hypermutator phenotype. Recurrently mutated genes also suggested that convergent microevolution and adaption of P. aeruginosa occur in inpatients, and genomic diversity is associated with differences in multiple-drug-resistance profiles. Taken together, our findings highlight the importance of longitudinal surveillance of nosocomial P. aeruginosa clones.

Longitudinal follow-up of HPV16 sequence after cervical infection: Low intrahost variation and no correlation with clinical evolution.

Human papillomavirus (HPV) 16 exhibits different variants that may differ in their carcinogenic risk. To identify some high-risk variants, we sequenced and compared HPV16 whole genomes obtained from a longitudinal cohort of 34 HPV16-infected women who had either spontaneously cleared their infection (clearance group or "C"), or developed cervical high-grade lesions following a viral persistence (group persistence or "P"). Phylogenetic analysis of paired samples obtained at the beginning (C0 or P0) and at the end (C2 or P2) of the follow-up (median intervals between C0-C2 and between P0-P2 were 16 and 36.5 months, respectively) revealed a low genetic variability within the host compared to the genetic interhost diversity. By comparing our HPV16 sequences to a reference sequence, we observed 301 different substitutions, more often transitions (60.9%) than transversions (39.1%), that occurred throughout the viral genome, but with a low frequency in E6 and E7 oncogenes (10 and 9 substitutions), suggesting a high conservation of these genes. Deletions and insertions were mostly observed in intergenic regions of the virus. The only significant substitution found between the subgroups C2 and P2 was observed in the L2 gene (L330F), with an unclear biological relevance. Our results suggest a low longitudinal intrahost evolution of HPV16 sequences and no correlation between genetic variations and clinical evolution.

Intra-host genetic population diversity: Role in emergence and persistence of drug resistance among Mycobacterium tuberculosis complex minor variants.

Drug resistant tuberculosis (DR-TB) is an important public health issue in different parts of the world. Mycobacterium tuberculosis complex variants (MTBC vars) preferentially infect certain hosts, limiting their distribution to different ecosystems. However, MTBC vars can infect other hosts beyond their preferred target potentially contributing to persistence of drug resistance (DR) in other niches. Here, we performed a comprehensive intra-host genetic analysis for the identification of DR-related mutations among all MTBC minor vars whole genome sequences (8,095 strains) publicly available worldwide. High confidence drug-resistance mutations in katG (isoniazid), rpsL (streptomycin), pncA (pyrazinamide), rpoB (rifampicin) and gyrA (fluoroquinolones) genes were identified among intrahost minor sub-populations in 197 different strains (2.43%) belonging to vars africanum, bovis, caprae, microti, orygis and pinnipedii. In addition, a three-dimensional structure modeling analysis to assess the role of novel mutations was also performed. Our findings highlight the importance of detecting discrete intra-host populations carrying DR mutations.

Low Intrahost and Interhost Genetic Diversity of Carnivore Protoparvovirus 1 in Domestic Cats during a Feline Panleukopenia Outbreak.

Feline panleukopenia (FPL), a highly contagious and frequently fatal disease of cats, is caused by Feline parvovirus (FPV) and Canine parvovirus (CPV). We characterised the diversity of these Carnivore protoparvovirus 1 variants in 18 faecal samples collected from domestic cats with FPL during an outbreak, using targeted parvoviral DNA metagenomics to a mean depth of >10,000 × coverage per site. All samples comprised FPV alone. Compared with the reference FPV genome, isolated in 1967, 44 mutations were detected. Ten of these were nonsynonymous, including 9 in nonstructural genes and one in VP1/VP2 (Val232Ile), which was the only one to exhibit interhost diversity, being present in five sequences. There were five other polymorphic nucleotide positions, all with synonymous mutations. Intrahost diversity at all polymorphic positions was low, with subconsensus variant frequencies (SVF) of <1% except for two positions (2108 and 3208) in two samples with SVF of 1.1–1.3%. Intrahost nucleotide diversity was measured across the whole genome (0.7–1.5%) and for each gene and was highest in the NS2 gene of four samples (1.2–1.9%). Overall, intrahost viral genetic diversity was limited and most mutations observed were synonymous, indicative of a low background mutation rate and strong selective constraints.

Probe Capture Enrichment Methods for HIV and HCV Genome Sequencing and Drug Resistance Genotyping.

Human immunodeficiency virus (HIV) infections remain a significant public health concern worldwide. Over the years, sophisticated sequencing technologies such as next-generation sequencing (NGS) have emerged and been utilized to monitor the spread of HIV drug resistance (HIVDR), identify HIV drug resistance mutations, and characterize transmission dynamics. Similar applications also apply to the Hepatitis C virus (HCV), another bloodborne viral pathogen with significant intra-host genetic diversity. Several advantages to using NGS over conventional Sanger sequencing include increased data throughput, scalability, cost-effectiveness when batched sample testing is performed, and sensitivity for quantitative detection of minority resistant variants. However, NGS alone may fail to detect genomes from pathogens present in low copy numbers. As with all sequencing platforms, the primary determinant in achieving quality sequencing data is the quality and quantity of the initial template input. Samples containing degraded RNA/DNA and/or low copy number have been a consistent sequencing challenge. To overcome this limitation probe capture enrichment is a method that has recently been employed to target, enrich, and sequence the genome of a pathogen present in low copies, and for compromised specimens that contain poor quality nucleic acids. It involves the hybridization of sequence-specific DNA or RNA probes to a target sequence, which is followed by an enrichment step via PCR to increase the number of copies of the targeted sequences after which the samples are subjected to NGS procedures. This method has been performed on pathogens such as bacteria, fungus, and viruses and allows for the sequencing of complete genomes, with high coverage. Post NGS, data analysis can be performed through various bioinformatics pipelines which can provide information on genetic diversity, genotype, virulence, and drug resistance. This article reviews how probe capture enrichment helps to increase the likelihood of sequencing HIV and HCV samples that contain low viral loads and/or are compromised.

Cover Image, Volume 94, Number 8, August 2022

Front Cover Caption: The cover image is based on the Research Article Intrahost variation and evolutionary dynamics of adenoviruses correlate to neutrophils in infected patients by Shuting Huo et al., https://doi.org/10.1002/jmv.27744.

Whole genome sequencing of clinical specimens reveals the genomic diversity of porcine reproductive and respiratory syndrome viruses emerging in China.

The Porcine reproductive and respiratory syndrome virus (PRRSV), a single (+) RNA virus, is characterized by high genome variability and constant evolution. Owing to increasingly complex mutations, there is a growing difficulty in accessing the whole genome. Additionally, there is limited knowledge on PRRSV intra-host nucleotide variants, which may reflect the complex viral-host dynamics. Here, we performed next-generation sequencing on four clinical lung tissues to reveal the genomic diversity and highlight virus-host interactions. The complete genomes of the HN0713 and GDYJ1224 strains shared 90.7% and 91.3% homology with the lineage 1 strain NADC30, respectively, while the GDGZ0408 and GDHY0425 strains shared 92.0% and 91.6% homology with the JXA1 strain, respectively. Recombination analysis showed that the ORF5-7 genes of the GDGZ0408 and GDHY0425 strains, whose complete genomes belong to lineage 8.7 based on the phylogenetic tree, are both recombined with lineage 3 strains. Furthermore, nsp3, nsp10-12, ORF2 genes and a part of the 3'-UTR of the GDHY0425 strain were provided by the lineage 5.1 strain. Two lineage 1 strains (GDYJ1224 and HN0713) were produced by a recombination of lineages 8.7 and 1. Additionally, the lineage 3 strain was associated with the recombinant HN0713 strain. We determined the intra-host single nucleotide variant frequencies and found more than 200 sites at a frequency of>1% in all samples. GDGZ0408 with parts of the nsp9 and nsp10 genes of HP-PRRSV lineage 8.7 presented more genetically diverse populations than others, indicating that lineage 8.7 might drive robust intra-host single nucleotide variants (iSNVs). Moreover, in the iSNV pools, nsp2 and ORF2a presented the highest mutation dynamic. Overall, this study provided evidence for the alarmingly increasing recombination and ever-changing evolutionary dynamics of PRRSV, and revealed the potential causes of vaccine escape, providing a novel insight into the nucleotide variant population in clinical samples.

SARS-CoV-2 Within-Host and in vitro Genomic Variability and Sub-Genomic RNA Levels Indicate Differences in Viral Expression Between Clinical Cohorts and in vitro Culture.

BackgroundLow frequency intrahost single nucleotide variants (iSNVs) of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) have been increasingly recognised as predictive indicators of positive selection. Particularly as growing numbers of SARS-CoV-2 variants of interest (VOI) and concern (VOC) emerge. However, the dynamics of subgenomic RNA (sgRNA) expression and its impact on genomic diversity and infection outcome remain poorly understood. This study aims to investigate and quantify iSNVs and sgRNA expression in single and longitudinally sampled cohorts over the course of mild and severe SARS-CoV-2 infection, benchmarked against an in vitro infection model.MethodsTwo clinical cohorts of SARS-CoV-2 positive cases in New South Wales, Australia collected between March 2020 and August 2021 were sequenced. Longitudinal samples from cases hospitalised due to SARS-CoV-2 infection (severe) (n = 16) were analysed and compared with cases that presented with SARS-CoV-2 symptoms but were not hospitalised (mild) (n = 23). SARS-CoV-2 genomic diversity profiles were also examined from daily sampling of culture experiments for three SARS-CoV-2 variants (Lineage A, B.1.351, and B.1.617.2) cultured in VeroE6 C1008 cells (n = 33).ResultsIntrahost single nucleotide variants were detected in 83% (19/23) of the mild cohort cases and 100% (16/16) of the severe cohort cases. SNP profiles remained relatively fixed over time, with an average of 1.66 SNPs gained or lost, and an average of 4.2 and 5.9 low frequency variants per patient were detected in severe and mild infection, respectively. sgRNA was detected in 100% (25/25) of the mild genomes and 92% (24/26) of the severe genomes. Total sgRNA expressed across all genes in the mild cohort was significantly higher than that of the severe cohort. Significantly higher expression levels were detected in the spike and the nucleocapsid genes. There was significantly less sgRNA detected in the culture dilutions than the clinical cohorts.Discussion and ConclusionThe positions and frequencies of iSNVs in the severe and mild infection cohorts were dynamic overtime, highlighting the importance of continual monitoring, particularly during community outbreaks where multiple SARS-CoV-2 variants may co-circulate. sgRNA levels can vary across patients and the overall level of sgRNA reads compared to genomic RNA can be less than 1%. The relative contribution of sgRNA to the severity of illness warrants further investigation given the level of variation between genomes. Further monitoring of sgRNAs will improve the understanding of SARS-CoV-2 evolution and the effectiveness of therapeutic and public health containment measures during the pandemic.

Tumors and Cytomegalovirus: An Intimate Interplay.

Human cytomegalovirus (HCMV) is a herpesvirus that alternates lytic and latent infection, infecting between 40 and 95% of the population worldwide, usually without symptoms. During its lytic cycle, HCMV can result in fever, asthenia, and, in some cases, can lead to severe symptoms such as hepatitis, pneumonitis, meningitis, retinitis, and severe cytomegalovirus disease, especially in immunocompromised individuals. Usually, the host immune response keeps the virus in a latent stage, although HCMV can reactivate in an inflammatory context, which could result in sequential lytic/latent viral cycles during the lifetime and thereby participate in the HCMV genomic diversity in humans and the high level of HCMV intrahost genomic variability. The oncomodulatory role of HCMV has been reported, where the virus will favor the development and spread of cancerous cells. Recently, an oncogenic role of HCMV has been highlighted in which the virus will directly transform primary cells and might therefore be defined as the eighth human oncovirus. In light of these new findings, it is critical to understand the role of the immune landscape, including the tumor microenvironment present in HCMV-harboring tumors. Finally, the oncomodulatory/oncogenic potential of HCMV could lead to the development of novel adapted therapeutic approaches against HCMV, especially since immunotherapy has revolutionized cancer therapeutic strategies and new therapeutic approaches are actively needed, particularly to fight tumors of poor prognosis.

Target Enrichment Metagenomics Reveals Human Pegivirus-1 in Pediatric Hematopoietic Stem Cell Transplantation Recipients.

Human pegivirus-1 (HPgV-1) is a lymphotropic human virus, typically considered nonpathogenic, but its infection can sometimes cause persistent viremia both in immunocompetent and immunosuppressed individuals. In a viral discovery research program in hematopoietic stem cell transplant (HSCT) pediatric patients, HPgV-1 was detected in 3 out of 14 patients (21.4%) using a target enrichment next-generation sequencing method, and the presence of the viruses was confirmed by agent-specific qRT-PCR assays. For the first time in this patient cohort, complete genomes of HPgV-1 were acquired and characterized. Phylogenetic analyses indicated that two patients had HPgV-1 genotype 2 and one had HPgV-1 genotype 3. Intra-host genomic variations were described and discussed. Our results highlight the necessity to screen HSCT patients and blood and stem cell donors to reduce the potential risk of HPgV-1 transmission.