Intra-host Diversity Research Articles

Intrahost variation and evolutionary dynamics of adenoviruses correlate to neutrophils in infected patients.

With deep sequencing of virus genomes within the hosts, intrahost single nucleotide variations (iSNVs) have been used for analyses of virus genome variation and evolution, which is indicated to correlate with viral pathogenesis and disease severity. Little is known about the features of iSNVs among DNA viruses. We performed the epidemiological and laboratory investigation of one outbreak of adenovirus. The whole genomes of viruses in both original oral swabs and cell-cultured virus isolates were deeply sequenced. We identified 737 iSNVs in the viral genomes sequenced from original samples and 46 viral iSNVs in cell-cultured isolates, with 33 iSNVs shared by original samples and cultured isolates. Meanwhile, we found these 33 iSNVs were shared by different patients, among which, three hot spot areas 6367-6401, 9213-9247, and 10 584-10 606 within the functional genes of the adenovirus genome were found. Notably, the substitution rates of iSNVs were closely correlated with the clinical and immune indicators of the patients. Especially a positive correlation to neutrophils was found, indicating a predictable biomarker of iSNV dynamics. Our findings demonstrated the neutrophil-correlated dynamic evolution features of the iSNVs within adenoviruses, which indicates a virus-host interaction during human infection of a DNA virus.

CD8+ T-Cell Exhaustion Phenotype in Chronic Hepatitis C Virus Infection Is Associated With Epitope Sequence Variation.

Background and AimsDuring chronic hepatitis C virus (HCV) infection, CD8+ T-cells become functionally exhausted, undergoing progressive phenotypic changes, i.e., overexpression of “inhibitory” molecules such as PD-1 (programmed cell death protein 1) and/or Tim-3 (T-cell immunoglobulin and mucin domain-containing molecule-3). The extreme intrahost genetic diversity of HCV is a major mechanism of immune system evasion, facilitating epitope escape. The aim of the present study was to determine whether T-cell exhaustion phenotype in chronic HCV infection is related to the sequence repertoire of NS3 viral immunodominant epitopes.MethodsThe study population was ninety prospective patients with chronic HCV genotype 1b infection. Populations of peripheral blood CD8+ T-cells expressing PD-1/Tim-3 were assessed by multiparametric flow cytometry, including HCV-specific T-cells after magnetic-based enrichment using MHC-pentamer. Autologous epitope sequences were inferred from next-generation sequencing. The correction of sequencing errors and genetic variants reconstruction was performed using Quasirecomb.ResultsThere was an interplay between the analyzed epitopes sequences and exhaustion phenotype of CD8+ T-cells. A predominance of NS31406 epitope sequence, representing neither prototype KLSGLGLNAV nor cross-reactive variants (KLSSLGLNAV, KLSGLGINAV or KLSALGLNAV), was associated with higher percentage of HCV-specific CD8+PD-1+Tim-3+ T-cells, P=0.0102. Variability (at least two variants) of NS31406 epitope sequence was associated with increased frequencies of global CD8+PD-1+Tim-3+ T-cells (P=0.0197) and lower frequencies of CD8+PD-1−Tim-3− T-cells (P=0.0079). In contrast, infection with NS31073 dominant variant epitope (other than prototype CVNGVCWTV) was associated with lower frequency of global CD8+PD-1+Tim-3+ T-cells (P=0.0054).ConclusionsOur results indicate that PD-1/Tim-3 receptor expression is largely determined by viral epitope sequence and is evident for both HCV-specific and global CD8+ T-cells, pointing to the importance of evaluating autologous viral epitope sequences in the investigation of CD8+ T-cell exhaustion in HCV infection.

Exploring a prolonged enterovirus C104 infection in a severely ill patient using nanopore sequencing

Chronic enterovirus infections can cause significant morbidity, particularly in immunocompromised patients. This study describes a fatal case associated with a chronic untypeable enterovirus infection in an immunocompromised patient admitted to a Dutch university hospital over nine months. We aimed to identify the enterovirus genotype responsible for the infection and to determine potential evolutionary changes. Long-read sequencing was performed using viral targeted sequence capture on four respiratory and one faecal sample. Phylogenetic analysis was performed using a maximum likelihood method, along with a root-to-tip regression and time-scaled phylogenetic analysis to estimate evolutionary changes between sample dates. Intra-host variant detection, using a Fixed Ploidy algorithm, and selection pressure, using a Fixed Effect Likelihood and a Mixed Effects Model of Evolution, were also used to explore the patient samples. Near-complete genomes of enterovirus C104 (EV-C104) were recovered in all respiratory samples but not in the faecal sample. The recovered genomes clustered with a recently reported EV-C104 from Belgium in August 2018. Phylodynamic analysis including ten available EV-C104 genomes, along with the patient sequences, estimated the most recent common ancestor to occur in the middle of 2005 with an overall estimated evolution rate of 2.97 × 10−3 substitutions per year. Although positive selection pressure was identified in the EV-C104 reference sequences, the genomes recovered from the patient samples alone showed an overall negative selection pressure in multiple codon sites along the genome. A chronic infection resulting in respiratory failure from a relatively rare enterovirus was observed in a transplant recipient. We observed an increase in single-nucleotide variations between sample dates from a rapidly declining patient, suggesting mutations are weakly deleterious and have not been purged during selection. This is further supported by the persistence of EV-C104 in the patient, despite the clearance of other viral infections. Next-generation sequencing with viral enrichment could be used to detect and characterise challenging samples when conventional workflows are insufficient.

Rescuing low frequency variants within intra-host viral populations directly from Oxford Nanopore sequencing data

Infectious disease monitoring on Oxford Nanopore Technologies (ONT) platforms offers rapid turnaround times and low cost. Tracking low frequency intra-host variants provides important insights with respect to elucidating within-host viral population dynamics and transmission. However, given the higher error rate of ONT, accurate identification of intra-host variants with low allele frequencies remains an open challenge with no viable computational solutions available. In response to this need, we present Variabel, a novel approach and first method designed for rescuing low frequency intra-host variants from ONT data alone. We evaluate Variabel on both synthetic data (SARS-CoV-2) and patient derived datasets (Ebola virus, norovirus, SARS-CoV-2); our results show that Variabel can accurately identify low frequency variants below 0.5 allele frequency, outperforming existing state-of-the-art ONT variant callers for this task. Variabel is open-source and available for download at: www.gitlab.com/treangenlab/variabel.

Low Pathogenicity H7N3 Avian Influenza Viruses Have Higher Within-Host Genetic Diversity Than a Closely Related High Pathogenicity H7N3 Virus in Infected Turkeys and Chickens.

Within-host viral diversity offers a view into the early stages of viral evolution occurring after a virus infects a host. In recent years, advances in deep sequencing have allowed for routine identification of low-frequency variants, which are important sources of viral genetic diversity and can potentially emerge as a major virus population under certain conditions. We examined within-host viral diversity in turkeys and chickens experimentally infected with closely related H7N3 avian influenza viruses (AIVs), specifically one high pathogenicity AIV (HPAIV) and two low pathogenicity AIV (LPAIVs) with different neuraminidase protein stalk lengths. Consistent with the high mutation rates of AIVs, an abundance of intra-host single nucleotide variants (iSNVs) at low frequencies of 2–10% was observed in all samples collected. Furthermore, a small number of common iSNVs were observed between turkeys and chickens, and between directly inoculated and contact-exposed birds. Notably, the LPAIVs have significantly higher iSNV diversities and frequencies of nonsynonymous changes than the HPAIV in both turkeys and chickens. These findings highlight the dynamics of AIV populations within hosts and the potential impact of genetic changes, including mutations in the hemagglutinin gene that confers the high pathogenicity pathotype, on AIV virus populations and evolution.

Multiorgan and Vascular Tropism of SARS-CoV-2

Although the respiratory tract is the main target of SARS-CoV-2, other tissues and organs are permissive to the infection. In this report, we investigated this wide-spectrum tropism by studying the SARS-CoV-2 genetic intra-host variability in multiple tissues. The virological and histological investigation of multiple specimens from a post-mortem COVID-19 patient was performed. SARS-CoV-2 genome was detected in several tissues, including the lower respiratory system, cardio-vascular biopsies, stomach, pancreas, adrenal gland, mediastinal ganglion and testicles. Subgenomic RNA transcripts were also detected, in favor of an active viral replication, especially in testicles. Ultra-deep sequencing allowed us to highlight several SARS-CoV-2 mutations according to tissue distribution. More specifically, mutations of the spike protein, i.e., V341A (18.3%), E654 (44%) and H655R (30.8%), were detected in the inferior vena cava. SARS-CoV-2 variability can contribute to heterogeneous distributions of viral quasispecies, which may affect the COVID-19 pathogeny.

Circulation, genomic characteristics, and evolutionary dynamics of class I Newcastle disease virus in China

ABSTRACT Newcastle disease caused by Newcastle disease virus (NDV) is one of the most serious threats to chickens and has two clinical forms, typical and atypical, caused by velogenic and lentogenic strains, respectively. To control the epidemic, many vaccines against velogenic class II NDVs have been introduced worldwide, but this has led to accelerated mutation of class II viruses under immune pressure and, on the other hand, to non-vaccine targeting class I NDVs becoming the dominant population in poultry. In this context, this study provided the first large-scale genomic epidemiological and quasispecies dynamic analysis of class I NDVs in China, and found that class I viruses that first appeared in East and South China have spread to central China and become the dominant class with an average evolutionary rate of 1.797 × 10−3. In addition, single nucleotide polymorphism and intra-host single nucleotide variation analyses show that HN and P genes have high mutation rates and may act as front-runners for NDV to expand their host range and enhance their virulence. This study also found that the class I NDV population has accumulated a number of mutations under positive selection and that six isolates with shortened C-terminal extensions of the HN protein are evolving toward increased virulence. These results not only enrich the research resources but also help us to better understand the dynamic evolution and mutational trends of NDV at the genomic level, which is crucial for monitoring, early warning, and controlling the outbreak of Newcastle disease.

Incorporating Within-Host Diversity in Phylogenetic Analyses for Detecting Clusters of New HIV Diagnoses.

BackgroundPhylogenetic analyses of HIV sequences are used to detect clusters and inform public health interventions. Conventional approaches summarize within-host HIV diversity with a single consensus sequence per host of the pol gene, obtained from Sanger or next-generation sequencing (NGS). There is growing recognition that this approach discards potentially important information about within-host sequence variation, which can impact phylogenetic inference. However, whether alternative summary methods that incorporate intra-host variation impact phylogenetic inference of transmission network features is unknown.MethodsWe introduce profile sampling, a method to incorporate within-host NGS sequence diversity into phylogenetic HIV cluster inference. We compare this approach to Sanger- and NGS-derived pol and near-whole-genome consensus sequences and evaluate its potential benefits in identifying molecular clusters among all newly-HIV-diagnosed individuals over six months at the largest HIV center in Rhode Island.ResultsProfile sampling cluster inference demonstrated that within-host viral diversity impacts phylogenetic inference across individuals, and that consensus sequence approaches can obscure both magnitude and effect of these impacts. Clustering differed between Sanger- and NGS-derived consensus and profile sampling sequences, and across gene regions.DiscussionProfile sampling can incorporate within-host HIV diversity captured by NGS into phylogenetic analyses. This additional information can improve robustness of cluster detection.

Assessment of Virological Contributions to COVID-19 Outcomes in a Longitudinal Cohort of Hospitalized Adults.

BackgroundWhile several demographic and clinical correlates of coronavirus disease 2019 (COVID-19) outcome have been identified, their relationship to virological and immunological parameters remains poorly defined.MethodsTo address this, we performed longitudinal collection of nasopharyngeal swabs and blood samples from a cohort of 58 hospitalized adults with COVID-19. Samples were assessed for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) viral load, viral genotype, viral diversity, and antibody titer. Demographic and clinical information, including patient blood tests and several composite measures of disease severity, was extracted from electronic health records.ResultsSeveral factors, including male sex, higher age, higher body mass index, higher 4C Mortality score, and elevated lactate dehydrogenase levels, were associated with intensive care unit admission. Of all measured parameters, only the retrospectively calculated median Deterioration Index score was significantly associated with death. While quantitative polymerase chain reaction cycle threshold (Ct) values and genotype of SARS-CoV-2 were not significantly associated with outcome, Ct value did correlate positively with C-reactive protein levels and negatively with D-dimer, lymphocyte count, and antibody titer. Intrahost viral genetic diversity remained constant through the disease course and resulted in changes in viral genotype in some participants.ConclusionsUltimately, these results suggest that worse outcomes are driven by immune dysfunction rather than by viral load and that SARS-CoV-2 evolution in hospital settings is relatively constant over time.

Limited genomic reconstruction of SARS-CoV-2 transmission history within local epidemiological clusters.

A detailed understanding of how and when severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) transmission occurs is crucial for designing effective prevention measures. Other than contact tracing, genome sequencing provides information to help infer who infected whom. However, the effectiveness of the genomic approach in this context depends on both (high enough) mutation and (low enough) transmission rates. Today, the level of resolution that we can obtain when describing SARS-CoV-2 outbreaks using just genomic information alone remains unclear. In order to answer this question, we sequenced forty-nine SARS-CoV-2 patient samples from ten local clusters in NW Spain for which partial epidemiological information was available and inferred transmission history using genomic variants. Importantly, we obtained high-quality genomic data, sequencing each sample twice and using unique barcodes to exclude cross-sample contamination. Phylogenetic and cluster analyses showed that consensus genomes were generally sufficient to discriminate among independent transmission clusters. However, levels of intrahost variation were low, which prevented in most cases the unambiguous identification of direct transmission events. After filtering out recurrent variants across clusters, the genomic data were generally compatible with the epidemiological information but did not support specific transmission events over possible alternatives. We estimated the effective transmission bottleneck size to be one to two viral particles for sample pairs whose donor–recipient relationship was likely. Our analyses suggest that intrahost genomic variation in SARS-CoV-2 might be generally limited and that homoplasy and recurrent errors complicate identifying shared intrahost variants. Reliable reconstruction of direct SARS-CoV-2 transmission based solely on genomic data seems hindered by a slow mutation rate, potential convergent events, and technical artifacts. Detailed contact tracing seems essential in most cases to study SARS-CoV-2 transmission at high resolution.

Persistent SARS-CoV-2 Infection in a Patient With Non-hodgkin Lymphoma: Intra-Host Genomic Diversity Analysis

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has caused a pandemic, threatening global public health. Several cases of persistent infection have been described, but there are few reports that compared the genetic variability among samples collected from the patient during infection. In the current study, we reported a viral genetic analysis of a diabetic male patient with Non-Hodgkin Lymphoma affected by persistent SARS-CoV-2 infection. We sequenced the patient-derived viral isolated both from oro/nasopharyngeal swab and VeroE6 cell line, collected from the same patient at different points of the infection. Due to the insufficient material of the second swab received, in order to obtain a complete coverage of the viral genome, it was convenient to perform a virus isolation after cell culture. Both genomes belonged to Pangolin Lineage B.1, Nextstrain clade 20A and GISAID clade G. The mutation spectrum predicted for the two viral genomes reveal three additionally mutations in the sequence of second sample when compared with mutations set identified in the first sample. Our findings show the evolution of the intra-host variability during the course of a long-lasting infection.

Genomic diversity and antimicrobial resistance of Prevotella species isolated from chronic lung disease airways.

Cystic fibrosis (CF) and chronic obstructive pulmonary disease (COPD) are characterized by increasingly frequent acute pulmonary exacerbations that reduce life quality and length. Human airways are home to a rich polymicrobial environment, which includes members of the obligately anaerobic genus Prevotella . Despite their commonness, surprisingly little is known about the prevalence, role, genomic diversity and antimicrobial resistance (AMR) potential of Prevotella species and strains in healthy and diseased airways. Here, we used comparative genomics to develop a real-time PCR assay to permit rapid Prevotella species identification and quantification from cultures and clinical specimens. Assay specificity was validated across a panel of Prevotella and non- Prevotella species, followed by PCR screening of CF and COPD respiratory-derived cultures. Next, 35 PCR-positive isolates were subjected to whole-genome sequencing. Of eight identified Prevotella species, P. histicola , P. melaninogenica , P. nanceiensis , P. salivae and P. denticola overlapped between participant cohorts. Phylogenomic analysis revealed considerable interhost but limited intrahost diversity, suggesting patient-specific lineages in the lower airways, probably from oral cavity aspirations. Correlation of phenotypic AMR profiles with AMR genes identified excellent correlation between tetQ presence and decreased doxycycline susceptibility, and ermF presence and decreased azithromycin susceptibility and clindamycin resistance. AMR rates were higher in the CF isolates, reflecting greater antibiotic use in this cohort. All tested Prevotella isolates were tobramycin-resistant, providing a potential selection method to improve Prevotella culture retrieval rates. Our addition of 35 airway-derived Prevotella genomes to public databases will enhance ongoing efforts to unravel the role of this diverse and enigmatic genus in both diseased and healthy lungs.

Intra-host analysis of hepaciviral glycoprotein evolution reveals signatures associated with viral persistence and clearance.

Even 30 years after the discovery of the hepatitis C virus (HCV) in humans there is still no vaccine available. Reasons for this include the high mutation rate of HCV, which allows the virus to escape immune recognition and the absence of an immunocompetent animal model for vaccine development. Phylogenetically distinct hepaciviruses (genus Hepacivirus, family Flaviviridae) have been isolated from diverse species, each with a narrow host range: the equine hepacivirus (EqHV) is the closest known relative of HCV. In this study, we used amplicon-based deep-sequencing to investigate the viral intra-host population composition of the genomic regions encoding the surface glycoproteins E1 and E2. Patterns of E1E2 substitutional evolution were compared in longitudinally sampled EqHV-positive sera of naturally and experimentally infected horses and HCV-positive patients. Intra-host virus diversity was higher in chronically than in acutely infected horses, a pattern which was similar in the HCV-infected patients. However, overall glycoprotein variability was higher in HCV compared to EqHV. Additionally, selection pressure in HCV populations was higher, especially within the N-terminal region of E2, corresponding to the hypervariable region 1 (HVR1) in HCV. An alignment of glycoprotein sequences from diverse hepaciviruses identified the HVR1 as a unique characteristic of HCV: hepaciviruses from non-human species lack this region. Together, these data indicate that EqHV infection of horses could represent a powerful surrogate animal model to gain insights into hepaciviral evolution and HCVs HVR1-mediated immune evasion strategy.

Viral infection and transmission in a large, well-traced outbreak caused by the SARS-CoV-2 Delta variant

The SARS-CoV-2 Delta variant has spread rapidly worldwide. To provide data on its virological profile, we here report the first local transmission of Delta in mainland China. All 167 infections could be traced back to the first index case. Daily sequential PCR testing of quarantined individuals indicated that the viral loads of Delta infections, when they first become PCR-positive, were on average ~1000 times greater compared to lineage A/B infections during the first epidemic wave in China in early 2020, suggesting potentially faster viral replication and greater infectiousness of Delta during early infection. The estimated transmission bottleneck size of the Delta variant was generally narrow, with 1-3 virions in 29 donor-recipient transmission pairs. However, the transmission of minor iSNVs resulted in at least 3 of the 34 substitutions that were identified in the outbreak, highlighting the contribution of intra-host variants to population-level viral diversity during rapid spread.

ViralFlow: A Versatile Automated Workflow for SARS-CoV-2 Genome Assembly, Lineage Assignment, Mutations and Intrahost Variant Detection.

The COVID-19 pandemic is driven by Severe Acute Respiratory Syndrome coronavirus 2 (SARS-CoV-2) that emerged in 2019 and quickly spread worldwide. Genomic surveillance has become the gold standard methodology used to monitor and study this fast-spreading virus and its constantly emerging lineages. The current deluge of SARS-CoV-2 genomic data generated worldwide has put additional pressure on the urgent need for streamlined bioinformatics workflows. Here, we describe a workflow developed by our group to process and analyze large-scale SARS-CoV-2 Illumina amplicon sequencing data. This workflow automates all steps of SARS-CoV-2 reference-based genomic analysis: data processing, genome assembly, PANGO lineage assignment, mutation analysis and the screening of intrahost variants. The pipeline is capable of processing a batch of around 100 samples in less than half an hour on a personal laptop or in less than five minutes on a server with 50 threads. The workflow presented here is available through Docker or Singularity images, allowing for implementation on laptops for small-scale analyses or on high processing capacity servers or clusters. Moreover, the low requirements for memory and CPU cores and the standardized results provided by ViralFlow highlight it as a versatile tool for SARS-CoV-2 genomic analysis.

Spatio-temporal dynamics of intra-host variability in SARS-CoV-2 genomes.

During the course of the COVID-19 pandemic, large-scale genome sequencing of SARS-CoV-2 has been useful in tracking its spread and in identifying variants of concern (VOC). Viral and host factors could contribute to variability within a host that can be captured in next-generation sequencing reads as intra-host single nucleotide variations (iSNVs). Analysing 1347 samples collected till June 2020, we recorded 16 410 iSNV sites throughout the SARS-CoV-2 genome. We found ∼42% of the iSNV sites to be reported as SNVs by 30 September 2020 in consensus sequences submitted to GISAID, which increased to ∼80% by 30th June 2021. Following this, analysis of another set of 1774 samples sequenced in India between November 2020 and May 2021 revealed that majority of the Delta (B.1.617.2) and Kappa (B.1.617.1) lineage-defining variations appeared as iSNVs before getting fixed in the population. Besides, mutations in RdRp as well as RNA-editing by APOBEC and ADAR deaminases seem to contribute to the differential prevalence of iSNVs in hosts. We also observe hyper-variability at functionally critical residues in Spike protein that could alter the antigenicity and may contribute to immune escape. Thus, tracking and functional annotation of iSNVs in ongoing genome surveillance programs could be important for early identification of potential variants of concern and actionable interventions.

Structural insights in cell-type specific evolution of intra-host diversity by SARS-CoV-2

As the global burden of SARS-CoV-2 infections escalates, so does the evolution of viral variants with increased transmissibility and pathology. In addition to this entrenched diversity, RNA viruses can also display genetic diversity within single infected hosts with co-existing viral variants evolving differently in distinct cell types. The BriSΔ variant, originally identified as a viral subpopulation from SARS-CoV-2 isolate hCoV-19/England/02/2020, comprises in the spike an eight amino-acid deletion encompassing a furin recognition motif and S1/S2 cleavage site. We elucidate the structure, function and molecular dynamics of this spike providing mechanistic insight into how the deletion correlates to viral cell tropism, ACE2 receptor binding and infectivity of this SARS-CoV-2 variant. Our results reveal long-range allosteric communication between functional domains that differ in the wild-type and the deletion variant and support a view of SARS-CoV-2 probing multiple evolutionary trajectories in distinct cell types within the same infected host.

Long range PCR-based deep sequencing for haplotype determination in mixed HCMV infections

BackgroundShort read sequencing has been used extensively to decipher the genome diversity of human cytomegalovirus (HCMV) strains, but falls short to reveal individual genomes in mixed HCMV strain populations. Novel third-generation sequencing platforms offer an extended read length and promise to resolve how distant polymorphic sites along individual genomes are linked. In the present study, we established a long amplicon PacBio sequencing workflow to identify the absolute and relative quantities of unique HCMV haplotypes spanning over multiple hypervariable sites in mixtures. Initial validation of this approach was performed with defined HCMV DNA templates derived from cell-culture enriched viruses and was further tested for its suitability on patient samples carrying mixed HCMV infections.ResultsTotal substitution and indel error rate of mapped reads ranged from 0.17 to 0.43% depending on the stringency of quality trimming. Artificial HCMV DNA mixtures were correctly determined down to 1% abundance of the minor DNA source when the total HCMV DNA input was 4 × 104 copies/ml. PCR products of up to 7.7 kb and a GC content < 55% were efficiently generated when DNA was directly isolated from patient samples. In a single sample, up to three distinct haplotypes were identified showing varying relative frequencies. Alignments of distinct haplotype sequences within patient samples showed uneven distribution of sequence diversity, interspersed by long identical stretches. Moreover, diversity estimation at single polymorphic regions as assessed by short amplicon sequencing may markedly underestimate the overall diversity of mixed haplotype populations.ConclusionsQuantitative haplotype determination by long amplicon sequencing provides a novel approach for HCMV strain characterisation in mixed infected samples which can be scaled up to cover the majority of the genome by multi-amplicon panels. This will substantially improve our understanding of intra-host HCMV strain diversity and its dynamic behaviour.

Two-step fitness selection for intra-host variations in SARS-CoV-2.

Spontaneous mutations introduce uncertainty into coronavirus disease 2019 (COVID-19) control procedures and vaccine development. Here, we perform a spatiotemporal analysis on intra-host single-nucleotide variants (iSNVs) in 402 clinical samples from 170 affected individuals, which reveals an increase in genetic diversity over time after symptom onset in individuals. Nonsynonymous mutations are overrepresented in the pool of iSNVs but underrepresented at the single-nucleotide polymorphism (SNP) level, suggesting a two-step fitness selection process: a large number of nonsynonymous substitutions are generated in the host (positive selection), and these substitutions tend to be unfixed as SNPs in the population (negative selection). Dynamic iSNV changes in subpopulations with different gender, age, illness severity, and viral shedding time displayed a varied fitness selection process among populations. Our study highlights that iSNVs provide a mutational pool shaping the rapid global evolution of the virus.

Long-read sequencing reveals the evolutionary drivers of intra-host diversity across natural RNA mycovirus infections.

Intra-host dynamics are a core component of virus evolution but most intra-host data come from a narrow range of hosts or experimental infections. Gaining broader information on the intra-host diversity and dynamics of naturally occurring virus infections is essential to our understanding of evolution across the virosphere. Here we used PacBio long-read HiFi sequencing to characterize the intra-host populations of natural infections of the RNA mycovirus Cryphonectria hypovirus 1 (CHV1). CHV1 is a biocontrol agent for the chestnut blight fungus (Cryphonectria parasitica), which co-invaded Europe alongside the fungus. We characterized the mutational and haplotypic intra-host virus diversity of thirty-eight natural CHV1 infections spread across four locations in Croatia and Switzerland. Intra-host CHV1 diversity values were shaped by purifying selection and accumulation of mutations over time as well as epistatic interactions within the host genome at defense loci. Geographical landscape features impacted CHV1 inter-host relationships through restricting dispersal and causing founder effects. Interestingly, a small number of intra-host viral haplotypes showed high sequence similarity across large geographical distances unlikely to be linked by dispersal.