Search

Soybean (Glycine max L.) is one of the most important crop plants in the Republic of Korea. Here, we conducted a soybean virome study. We harvested a total of 172 soybean leaf samples showing disease symptoms from major soybean-growing regions in the Republic of Korea. Individual samples were examined for virus infection by RT-PCR. Moreover, we generated eight libraries representing eight provinces by pooling samples and four libraries from single samples. RNA-seq followed by bioinformatics analyses revealed 10 different RNA viruses infecting soybean. The proportion of viral reads in each transcriptome ranged from 0.2 to 31.7%. Coinfection of different viruses in soybean plants was very common. There was a single dominant virus in each province, and this geographical difference might be related to the soybean seeds that transmit viruses. In this study, 32 viral genome sequences were assembled and successfully used to analyze the phylogenetic relationships and quasispecies nature of the identified RNA viruses. Moreover, RT-PCR with newly developed primers confirmed infection of the identified viruses in each library. Taken together, our soybean virome study provides a comprehensive overview of viruses infecting soybean in eight geographical regions in the Republic of Korea and four single soybean plants in detail.

Highlights • SARS-CoV-2 exhibits intra-host small- and large-scale genomic variability. • SNVs are collocalized with probes and primers used in molecular diagnostic assays. • SARS-CoV-2 Spike (S) gene host a potential recombination hot-spot., In December 2019, an outbreak of atypical pneumonia (Coronavirus disease 2019 -COVID-19) associated with a novel coronavirus (SARS-CoV-2) was reported in Wuhan city, Hubei province, China. The outbreak was traced to a seafood wholesale market and human to human transmission was confirmed. The rapid spread and the death toll of the new epidemic warrants immediate intervention. The intra-host genomic variability of SARS-CoV-2 plays a pivotal role in the development of effective antiviral agents and vaccines, as well as in the design of accurate diagnostics. We analyzed NGS data derived from clinical samples of three Chinese patients infected with SARS-CoV-2, in order to identify small- and large-scale intra-host variations in the viral genome. We identified tens of low- or higher- frequency single nucleotide variations (SNVs) with variable density across the viral genome, affecting 7 out of 10 protein-coding viral genes. The majority of these SNVs (72/104) corresponded to missense changes. The annotation of the identified SNVs but also of all currently circulating strain variations revealed colocalization of intra-host as well as strain specific SNVs with primers and probes currently used in molecular diagnostics assays. Moreover, we de-novo assembled the viral genome, in order to isolate and validate intra-host structural variations and recombination breakpoints. The bioinformatics analysis disclosed genomic rearrangements over poly-A / poly-U regions located in ORF1ab and spike (S) gene, including a potential recombination hot-spot within S gene. Our results highlight the intra-host genomic diversity and plasticity of SARS-CoV-2, pointing out genomic regions that are prone to alterations. The isolated SNVs and genomic rearrangements reflect the intra-patient capacity of the polymorphic quasispecies, which may arise rapidly during the outbreak, allowing immunological escape of the virus, offering resistance to anti-viral drugs and affecting the sensitivity of the molecular diagnostics assays.

Coxsackievirus A8 (CV-A8) is one of the pathogens associated with hand, foot and mouth disease (HFMD) and herpangina (HA), occasionally leading to severe neurological disorders such as acute flaccid paralysis (AFP). Only one study aimed at CV-A8 has been published to date, and only 12 whole-genome sequences are publicly available. In this study, complete genome sequences from 11 CV-A8 strains isolated from HFMD patients in extensive regions from China between 2013 and 2018 were determined, and all sequences from GenBank were retrieved. A phylogenetic analysis based on a total of 34 complete VP1 sequences of CV-A8 revealed five genotypes: A, B, C, D and E. The newly emerging genotype E presented a highly phylogenetic divergence compared with the other genotypes and was composed of the majority of the strains sequenced in this study. Markov chain Monte Carlo (MCMC) analysis revealed that genotype E has been evolving for nearly a century and somehow arose in approximately 2010. The Bayesian skyline plot showed that the population size of CV-A8 has experienced three dynamic fluctuations since 2001. Amino acid residues of VP1100N, 103Y, 240T and 241V, which were embedded in the potential capsid loops of genotype E, might enhance genotype E adaption to the human hosts. The CV-A8 whole genomes displayed significant intra-genotypic genetic diversity in the non-capsid region, and a total of six recombinant lineages were detected. The Chinese viruses from genotype E might have emerged recently from recombining with European CV-A6 strains. CV-A8 is a less important HFMD pathogen, and the capsid gene diversity and non-capsid recombination variety observed in CV-A8 strains indicated that the constant generation of deleterious genomes and a constant selection pressure against these deleterious mutations is still ongoing within CV-A8 quasispecies. It is possible that CV-A8 could become an important pathogen in the HFMD spectrum in the future. Further surveillance of CV-A8 is greatly needed.

Classical swine fever virus (CSFV) contains a specific motif within the E2 glycoprotein that differs between strains of different virulence. In the highly virulent CSFV strain Koslov, this motif comprises residues S763/L764 in the polyprotein. However, L763/P764 represent the predominant alleles in published CSFV genomes. In this study, changes were introduced into the CSFV strain Koslov (here called vKos_SL) to generate modified CSFVs with substitutions at residues 763 and/or 764 (vKos_LL, vKos_SP, and vKos_LP). The properties of these mutant viruses, in comparison to those of vKos_SL, were determined in pigs. Each of the viruses was virulent and induced typical clinical signs of CSF, but the vKos_LP strain produced them significantly earlier. Full-length CSFV cDNA amplicons (12.3 kb) derived from sera of infected pigs were deep sequenced and cloned to reveal the individual haplotypes that contributed to the single-nucleotide polymorphism (SNP) profiles observed in the virus population. The SNP profiles for vKos_SL and vKos_LL displayed low-level heterogeneity across the entire genome, whereas vKos_SP and vKos_LP displayed limited diversity with a few high-frequency SNPs. This indicated that vKos_SL and vKos_LL exhibited a higher level of fitness in the host and more stability at the consensus level, whereas several consensus changes were observed in the vKos_SP and vKos_LP sequences, pointing to adaptation. For each virus, only a subset of the variants present within the virus inoculums were maintained in the infected pigs. No clear tissue-dependent quasispecies differentiation occurred within inoculated pigs; however, clear evidence for transmission bottlenecks to contact animals was observed, with subsequent loss of sequence diversity. IMPORTANCE The surface-exposed E2 protein of classical swine fever virus is required for its interaction with host cells. A short motif within this protein varies between strains of different virulence. The importance of two particular amino acid residues in determining the properties of a highly virulent strain of the virus has been analyzed. Each of the different viruses tested proved highly virulent, but one of them produced earlier, but not more severe, disease. By analyzing the virus genomes present within infected pigs, it was found that the viruses which replicated within inoculated animals were only a subset of those within the virus inoculum. Furthermore, following contact transmission, it was shown that a very restricted set of viruses had transferred between animals. There were no significant differences in the virus populations present in various tissues of the infected animals. These results indicate mechanisms of virus population change during transmission between animals.

Senecavirus A (SVA), also known as Seneca Valley virus, belongs to the genus Senecavirus in the family Picornaviridae. In this study, a China SVA isolate (CH-LX-01-2016) was rescued from its cDNA clone, and then identified by RT-PCR, indirect immunofluorescence assay and mass spectrometry. The rescued SVA could separately induce typical plaque formations and cytopathic effects in cell monolayers. In order to uncover its evolutionary dynamics, the SVA was subjected to eighty serial passages in vitro. Its progenies per ten passages were analyzed by next-generation sequencing (NGS). The NGS analyses showed that neither sequence-deleting nor -inserting phenotype was detectable in eight progenies, within which a total of forty-one intra-host single-nucleotide variations (SNVs) arose with passaging. Almost all SNVs were identified as the single-nucleotide polymorphism with mixture of two nucleotides. SNVs led to eighteen nonsynonymous mutations, out of which sixteen could directly reflect their own frequencies of amino acid mutation, due to only one SNV occurring in their individual codons. Compared with its parental virus without passaging, the passage-80 SVA progeny had formed a viral quasispecies, as evidenced by a total of twenty-eight SNVs identified in it.

Numerous factors have been identified to influence susceptibility to SARS-CoV-2 infection and disease severity. Cancer patients are more prone to clinically evolve to more severe COVID-19 conditions, but the determinants of such a more severe outcome remain largely unknown. We have determined the full-length SARS-CoV-2 genomic sequences of cancer patients and healthcare workers (HCW; non-cancer controls) by deep sequencing and investigated the within-host viral quasispecies of each infection, quantifying intrahost genetic diversity. Naso- and oropharyngeal SARS-CoV-2+ swabs from 57 cancer patients and 14 healthcare workers (HCW) from the Brazilian Cancer Institute were collected in April–May 2020. Complete genome amplification using ARTIC network V3 multiplex primers was performed followed by next-generation sequencing. Assemblies were conducted in Geneious R11, where consensus sequences were extracted and intrahost single nucleotide variants (iSNVs) were identified. Maximum likelihood phylogenetic analysis was performed using PhyMLv.3.0 and lineages were classified using Pangolin and CoV-GLUE. Phylogenetic analysis showed that all but one strain belonged to clade B1.1. Four genetically linked mutations known as the globally dominant SARS-CoV-2 haplotype (C241T, C3037T, C14408T and A23403G) were found in the majority of consensus sequences. SNV signatures of previously characterized Brazilian genomes were also observed in most samples. Another 85 SNVs were found at a lower frequency (1.4-19.7%). Cancer patients displayed a significantly higher intrahost viral genetic diversity compared to HCW (p = 0.009). Intrahost genetic diversity in cancer patients was independent of SARS-CoV-2 Ct values, and was not associated with disease severity, use of corticosteroids, or use of antivirals, characteristics that could influence viral diversity. Such a feature may explain, at least in part, the more adverse outcomes to which cancer/COVID-19 patients experience.Author SummaryCancer patients are more prone to clinically evolve to more severe COVID-19 conditions, but the determinants of such a more severe outcome remain largely unknown. In this study, phylogenetic and variation analysis of SARS-CoV-2 genomes from cancer patients and non-cancer healthcare workers at the Brazilian National Cancer Institute were characterized by deep sequencing. Viral genomes showed signatures characteristic of Brazilian viruses, consistent with the hypothesis of local, community transmission rather than virus importation from abroad. Despite most genomes in patients and healthcare workers belonging to the same lineage, intrahost variability was higher in cancer patients when compared to non-cancer counterparts. The intrahost genomic diversity analysis presented in our study highlights the relaxed evolution of SARS-CoV-2 in a vulnerable population of cancer patients. The high number of minor variations can result in the selection of immune escape variants, resistance to potential drugs, and/or increased pathogenicity. The impact of this higher intrahost variability over time warrants further investigation.

The viruses of the family Coronaviridae are ubiquitous in nature due to their existence in a wide spectrum of mammals and avian species. The coronaviruses, as RNA viruses, exist as quasispecies because of their high rate of mutations. This review elaborates on the pathogenesis and the developed vaccines of most of the ubiquitous coronavirus' diseases, mainly bovine, dromedary camel, porcine, feline, canine, and avian coronaviruses. The review emphasizes the significant setbacks in the full exploitation of most of the pathogenesis of the coronavirus' diseases, raising the prospect of effective vaccines for these diseases. The therapeutical trials for the treatment of SARS-CoV2 and the setbacks of these trials are also addressed. The review draws attention to the lessons accumulated from the large number of studies of the pathogenesis of animals and birds' coronaviruses and their vaccines, particularly the bovine, feline, and avian coronaviruses. The lessons drawn from the studies will have an immense influence on how the human coronaviruses pathogenesis and vaccine development will proceed. In addition, the extensive efforts to designate suitable animal models to study the lately emerged human coronaviruses are one of the invaluable contributions carried out by veterinarian scientists. Finally, factors and determinants that contribute to the possibility of emerging new coronavirus zoonotic disease are elaborated on and a call goes out to urge transdisciplinary collaboration in the implementation of the "One Health" concept.

One of the features distinguishing SARS-CoV-2 from its more pathogenic counterpart SARS-CoV is the presence of premature stop codons in its ORF3b gene. Here, we show that SARS-CoV-2 ORF3b is a potent interferon antagonist, suppressing the induction of type I interferon more efficiently than its SARS-CoV ortholog. Phylogenetic analyses and functional assays reveal that SARS-CoV-2-related viruses from bats and pangolins also encode truncated ORF3b gene products with strong anti-interferon activity. Furthermore, analyses of approximately 17,000 SARS-CoV-2 sequences identify a natural variant in which a longer ORF3b reading frame was reconstituted. This variant was isolated from two patients with severe disease and further increased the ability of ORF3b to suppress interferon induction. Thus, our findings not only help to explain the poor interferon response in COVID-19 patients but also describe the emergence of natural SARS-CoV-2 quasispecies with an extended ORF3b gene that may potentially affect COVID-19 pathogenesis., Graphical Abstract, Highlights • ORF3b proteins of SARS-CoV-2 and related animal viruses are IFN antagonists • SARS-CoV-2 ORF3b suppresses IFN more efficiently than its SARS-CoV ortholog • The anti-IFN activity of ORF3b depends on the length of its C terminus • An ORF3b with increased IFN antagonism was isolated from two severe COVID-19 cases, COVID-19 pathogenesis is characterized by impaired IFN responses. Konno et al. identify ORF3b proteins of SARS-CoV-2 and related animal viruses as IFN antagonists. Their anti-IFN activity depends on the C-terminal length, and a natural ORF3b variant with increased IFN-suppressive activity was isolated from two severe COVID-19 cases.

Type 1 interferons (IFN-I) are potent innate antiviral effectors that constrain HIV-1 transmission. However, harnessing these cytokines for HIV-1 cure strategies has been hampered by an incomplete understanding of their anti-viral activities at later stages of infection. Here, we characterized the IFN-I sensitivity of 500 clonally-derived HIV-1 isolates from plasma and CD4+ T cells of 26 individuals sampled longitudinally following transmission and/or after antiretroviral therapy (ART) and analytical treatment interruption (ATI). Determining the concentration of IFNα2 and IFNβ that reduced HIV-1 replication by 50% (IC50), we found remarkably consistent changes in the sensitivity of viruses to IFN-I inhibition, both across individuals and over time. IFN-I resistance was uniformly high during acute infection, decreased in all subjects in the first year post-infection, was reacquired concomitant with CD4+ T cell loss, and remained elevated in subjects with accelerated disease. Isolates obtained by viral outgrowth during suppressive ART were relatively IFN-I sensitive, resembling viruses circulating just prior to ART initiation. However, viruses that rebounded following treatment interruption displayed the highest levels of IFNα2 and IFNβ resistance observed at any time during the infection course. These findings indicate a dynamic interplay between host innate immune responses and the evolving HIV-1 quasispecies, with the relative contribution of IFN-I to HIV-1 control impacted by both ART and ATI. Although elevated at transmission, IFN-mediated pressures are the highest during viral rebound, limiting the viruses that successfully reactivate from latency.One Sentence SummaryHIV-1 resistance to IFN-I is highest during acute infection and following analytic treatment interruption, indicating a dynamic interplay between host innate immunity and virus biology.

A novel coronavirus disease (COVID-19) caused by SARS-CoV-2 has been pandemic worldwide. The genetic dynamics of quasispecies afford RNA viruses a great fitness on cell tropism and host range. However, no quasispecies data of SARS-CoV-2 have been reported yet. To explore quasispecies haplotypes and its transmission characteristics, we carried out single-molecule real-time (SMRT) sequencing of the full-length of SARS-CoV-2 spike gene within 14 RNA samples from 2 infection clusters, covering first-to third-generation infected-patients. We observed a special quasispecies structure of SARS-CoV-2 (modeled as ‘One-King’): one dominant haplotype (mean abundance ~70.15%) followed by numerous minor haplotypes (mean abundance < 0.10%). We not only discovered a novel dominant haplotype of F1040but also realized that minor quasispecies were also worthy of attention. Notably, some minor haplotypes (like F1040and currently pandemic one G614) could potentially reveal adaptive and converse into the dominant one. However, minor haplotypes exhibited a high transmission bottleneck (~6% could be stably transmitted), and the new adaptive/dominant haplotypes were likely originated from genetic variations within a host rather than transmission. The evolutionary rate was estimated as 2.68-3.86 × 10−3per site per year, which was larger than the estimation at consensus genome level. The ‘One-King’ model and conversion event expanded our understanding of the genetic dynamics of SARS-CoV-2, and explained the incomprehensible phenomenon at the consensus genome level, such as limited cumulative mutations and low evolutionary rate. Moreover, our findings suggested the epidemic strains may be multi-host origin and future traceability would face huge difficulties.

Hepatitis C virus (HCV) infection remains a major worldwide health problem despite development of highly effective direct-acting antivirals. HCV rapidly evolves upon acute infection and generates multiple viral variants (quasispecies), leading to immune evasion and persistent viral infection. Identification of epitopes of broadly neutralizing anti-HCV antibodies (nAbs) is critical to guide HCV vaccine development. In this study, we developed a new reverse genetics system for HCV infection based on trans-complementation of viral structural proteins. The HCV genome (JFH1 strain) lacking the structural protein-coding sequence can be efficiently rescued by ectopic expression of core-E1-E2-p7-NS2 (core-NS2) or core-E1-E2-p7 (core-p7) in trans, leading to production of single-round infectious virions designated HCVΔS. JFH1-based HCVΔS can be also rescued by expressing core-NS2 of other HCV genotypes, rendering it an efficient tool to display the structural proteins of HCV strains of interests. Furthermore, we successfully rescued HCVΔS with structural proteins from clinical isolates. Multiple viral structural proteins with different sensitivities to nAbs were identified from a same patient serum, demonstrating the genetic diversity of HCV quasispecies in vivo. Interestingly, the structural protein-coding sequences of highly divergent viral quasispecies from the same patient can be clustered based on their hypervariable region 1 (HVR1) in viral envelope protein E2, which critically dictates the sensitivity to neutralizing antibodies. In summary, we developed a novel reverse genetics system that efficiently displays viral structural proteins from HCV clinical isolates, and analysis of quasispecies from the same patient using this system demonstrated that E2 HVR1 is the major determinant of viral evolution in vivo. IMPORTANCE A cell culture model that can recapitulate the diversity of HCV quasispecies in patients is important for analysis of neutralizing epitopes and HCV vaccine development. In this study, we developed a new reverse genetics system for HCV infection based on trans-complementation of viral structural proteins (HCVΔS). This system can be used to display structural proteins of HCV strains of multiple genotypes as well as clinical isolates. By using this system, we showed that multiple different HCV structural proteins from a same patient were displayed on HCVΔS. Interestingly, these variant structural proteins within the same patient can be classified according to the sequence of HVR1in E2, which dictates viral sensitivity to nAbs and viral evolution in vivo. Our work provided a new tool to study highly divergent HCV quasispecies and shed light on underlying mechanisms driving HCV evolution.

Background Hepatitis B virus (HBV) infection is one of the main leading causes of hepatocellular carcinoma (HCC) worldwide. However, it remains uncertain how the reverse-transcriptase (rt) gene contributes to HCC progression. Methods We enrolled a total of 307 patients with chronic hepatitis B (CHB) and 237 with HBV-related HCC from 13 medical centers. Sequence features comprised multidimensional attributes of rt nucleic acid and rt/s amino acid sequences. Machine-learning models were used to establish HCC predictive algorithms. Model performances were tested in the training and independent validation cohorts using receiver operating characteristic curves and calibration plots. Results A random forest (RF) model based on combined metrics (10 features) demonstrated the best predictive performances in both cross and independent validation (AUC, 0.96; accuracy, 0.90), irrespective of HBV genotypes and sequencing depth. Moreover, HCC risk scores for individuals obtained from the RF model (AUC, 0.966; 95% confidence interval, .922–.989) outperformed α-fetoprotein (0.713; .632–.784) in distinguishing between patients with HCC and those with CHB. Conclusions Our study provides evidence for the first time that HBV rt sequences contain vital HBV quasispecies features in predicting HCC. Integrating deep sequencing with feature extraction and machine-learning models benefits the longitudinal surveillance of CHB and HCC risk assessment.

Objectives Human immunodeficiency virus (HIV) quasispecies diversity presents a large barrier to the eradication of HIV. The aim of this study was to investigate intrahost HIV quasispecies diversity and evolutionary patterns underpinning the mechanisms of viral pathogenesis during antiretroviral therapy (ART). Methods Forty-five participants with HIV-1 infection were enrolled in a follow-up cohort for >84 months in 2004, and received a lamivudine-based first-line ART regimen. Blood samples were collected every 6 months to measure viral load and CD4+ cell count. Ultra-deep sequencing and phylogenetic analysis were used to characterize the dynamics governing quasispecies diversity of HIV-1 circulating between plasma RNA and cellular DNA of participants with treatment failure (TF, n = 20) or virologic suppression (VS, n = 25). Results Analysis of the distribution of intrahost single-nucleotide variations (iSNVs) and their mutated allele frequencies revealed that approximately 65% of the quasispecies co-occurred in plasma HIV RNA and cellular DNA either before or after ART. The number and frequency of iSNVs are more representative of intrahost HIV diversity, and have better generalizability than phylogenetic inference by measurement of phylogenetic associations. Furthermore, drug-resistance-associated mutations (DRAMs) accumulated to high levels, dramatically increasing the DRAM-to-total-mutation ratio for TF patients. Linear regression analysis revealed that emergent mutations accumulated faster in TF patients compared with VS patients, at a rate of 0.02 mutations/day/kb. Conclusions Based on iSNV analysis, the results demonstrate the dynamics of intrahost HIV quasispecies diversity in patients on ART, and provide a novel insight into the persistence of HIV and development of DRAMs.

Superspreading events shape the COVID-19 pandemic. Here we provide a national-scale analysis of SARS-CoV-2 outbreaks in Austria, a country that played a major role for virus transmission across Europe and beyond. Capitalizing on a national epidemiological surveillance system, we performed deep whole-genome sequencing of virus isolates from 576 samples to cover major Austrian SARS-CoV-2 clusters. Our data chart a map of early viral spreading in Europe, including the path from low-frequency mutations to fixation. Detailed epidemiological surveys enabled us to calculate the effective SARS-CoV-2 population bottlenecks during transmission and unveil time-resolved intra-patient viral quasispecies dynamics. This study demonstrates the power of integrating deep viral genome sequencing and epidemiological data to better understand how SARS-CoV-2 spreads through populations.Graphical Abstract

Despite the high virological response rates achieved with current directly-acting antiviral agents (DAAs) against hepatitis C virus (HCV), around 2% to 5% of treated patients do not achieve a sustained viral response. Identification of amino acid substitutions associated with treatment failure requires analytical designs, such as subtype-specific ultra-deep sequencing (UDS) methods for HCV characterization and patient management. Using this procedure, we have identified six highly represented amino acid substitutions (HRSs) in NS5A and NS5B of HCV from 220 patients who failed therapy, which are not bona fide resistance-associated substitutions (RAS). They were present frequently in basal and post-treatment virus of patients who failed therapy to different DAA-based therapies. Contrary to several RAS, HRSs belong to the acceptable subset of substitutions according to the PAM250 replacement matrix. Coherently, their mutant frequency, measured by the number of deep sequencing reads within the HCV quasispecies that encode the relevant substitutions, ranged between 90% and 100% in most cases. Also, they have limited predicted disruptive effects on the three-dimensional structures of the proteins harboring them. Possible mechanisms of HRS origin and dominance, as well as their potential predictive value of treatment response are discussed.

Next-generation sequencing (NGS) has revolutionized the scale and depth of biomedical sciences. Because of its unique ability for the detection of sub-clonal variants within genetically diverse populations, NGS has been successfully applied to analyze and quantify the exceptionally-high diversity within viral quasispecies, and many low-frequency drug- or vaccine-resistant mutations of therapeutic importance have been discovered. Although many works have intensively discussed the latest NGS approaches and applications in general, none of them has focused on applying NGS in viral quasispecies studies, mostly due to the limited ability of current NGS technologies to accurately detect and quantify rare viral variants. Here, we summarize several error-correction strategies that have been developed to enhance the detection accuracy of minority variants. We also discuss critical considerations for preparing a sequencing library from viral RNAs and for analyzing NGS data to unravel the mutational landscape.

Infection with human cytomegalovirus (HCMV) can cause severe complications in immunocompromised individuals and congenitally infected children. Characterizing heterogeneous viral populations and their evolution by high-throughput sequencing of clinical specimens requires the accurate assembly of individual strains or sequence variants and suitable variant calling methods. However, the performance of most methods has not been assessed for populations composed of low divergent viral strains with large genomes, such as HCMV. In an extensive benchmarking study, we evaluated 15 assemblers and six variant callers on ten lab-generated benchmark data sets created with two different library preparation protocols, to identify best practices and challenges for analyzing such data.Most assemblers, especially metaSPAdes and IVA, performed well across a range of metrics in recovering abundant strains. However, only one, Savage, recovered low abundant strains and in a highly fragmented manner. Two variant callers, LoFreq and VarScan2, excelled across all strain abundances. Both shared a large fraction of false positive (FP) variant calls, which were strongly enriched in T to G changes in a “G.G” context. The magnitude of this context-dependent systematic error is linked to the experimental protocol. We provide all benchmarking data, results and the entire benchmarking workflow named QuasiModo, Quasispecies Metric determination on omics, under the GNU General Public License v3.0 (https://github.com/hzi-bifo/Quasimodo), to enable full reproducibility and further benchmarking on these and other data.

Human astroviruses (HAstVs) were first identified in 1975 and can be classified into three clades: classic HAstVs (HAstV 1-8), MLB (MLB1-3) and VA (VA1-5), with MLB and VA were newly identified. Recombination and a high mutation rate make HAstV as one of the rapidly evolving infectious agents. This study reported a novel identified recombinant human astrovirus (Y/1-CHN) and its long existence in two immunocompromised patients with diarrhea following allogeneic hematopoietic stem cell transplantation (allo-HSCT). The identified Yu/1-CHN genome contains 6801 base pairs encoding three open reading frames, with ORF1a best hit to the HAstV1 (Pune strain, 97 % nucleotide identity), while ORF1b and ORF2 best hit to HAstV-5 (DL30 strain, 99 % nucleotide identity). Possible recombination breakpoint was predicted to be located in the boundary of ORF1a and ORF1b. Different quasispecies were found in the host, and the dN/dS ratios of the S and P domains were determined to be 1.189 and 1.444, respectively, suggesting a positive selection existed. Fecal samples collected in different clinical phases from the two patients were all positive for Yu/1-CHN, suggesting a long existence of the virus in the host. It was indicated that immunocompromised patients may a reservoir for astrovirus, their excreta should be monitored even after discharge from hospital.

Serum hepatitis B virus (HBV) mutations can predict hepatocellular carcinoma (HCC) occurrence. We aimed to clarify if HBV evolves synchronously in the sera, adjacent liver and tumors and predict HCC prognosis equally. A total of 203 HBV-positive HCC patients with radical hepatectomy in Shanghai, China, during 2011–15 were enrolled in this prospective study. Quasispecies complexity (QC) in HBV core promoter region was assessed using clone-based sequencing. We performed RNA sequencing on tumors and paired adjacent tissues of another 15 HCC patients and analyzed it with three public data sets containing 127 samples. HBV QC was positively correlated to APOBEC3s’ expression level (r = 0.28, P < 0.001), higher in the adjacent tissues than in the tumors (P = 6.50e-3), and higher in early tumors than in advanced tumors (P = 0.039). The evolutionary distance between the sera-derived HBV strains and the tumor-derived ones was significantly longer than that between the sera-derived ones and the adjacent tissue-derived ones (P < 0.001). Multivariate Cox regression analyses indicated that high HBV QC in the sera predicted an unfavorable overall survival (P = 0.002) and recurrence-free survival (RFS; P = 0.004) in HCC, whereas, in the tumors, it predicted a favorable RFS (P < 0.001). APOBECs-related HBV mutations, including G1764A, were more frequent in the sera than in the adjacent tissues. High-frequent A1762T/G1764A in the sera predicted an unfavorable RFS (P < 0.001), whereas, in the tumors, it predicted a favorable RFS (P = 0.035). In conclusion, HBV evolves more advanced in the sera than in the tumors. HBV QC and A1762T/G1764A in the sera and tumors have contrary prognostic effects in HCC.

The evolutionary dynamics of a virus can differ within hosts and across populations. Studies of within-host evolution provide an important link between experimental studies of virus evolution and large-scale phylodynamic analyses. They can determine the extent to which global processes are recapitulated on local scales and how accurately experimental infections model natural ones. They may also inform epidemiologic models of disease spread and reveal how host-level dynamics contribute to a virus's evolution at a larger scale. Over the last decade, advances in viral sequencing have enabled detailed studies of viral genetic diversity within hosts. I review how within-host diversity is sampled, measured, and expressed, and how comparative studies of viral diversity can be leveraged to elucidate a virus's evolutionary dynamics. These concepts are illustrated with detailed reviews of recent research on the within-host evolution of influenza virus, dengue virus, and cytomegalovirus.

Next-generation sequencing (NGS) based virome analyses of mRNA and sRNA have recently become a routine approach for reliable detection of plant viruses and viroids. In the present study we identified the viral/viroidal spectrum of several Indian grapevine cultivars and reconstructed their whole genomes using the publically available mRNAome and sRNAome datasets. Twenty three viruses and viroids (including two variants of grapevine leafroll associated virus 4) were identified from two tissues (fruit peels and young leaves) of three cultivars among which nine unique grapevine viruses and viroids were identified for the first time in India. Irrespective of the assemblers and tissues used, the mRNA based approach identified more acellular pathogens than the sRNA based approach across cultivars. Further, the mRNAome was on par with the whole transcriptome in viral identification. Through de novo assembly of transcriptomes followed by mapping against reference genome, we reconstructed 19 complete/near complete genomes of identified viruses and viroids. The reconstructed viral genomes included four larger RNA genomes (>13 kb), a DNA genome (RG grapevine geminivirus A), a divergent genome (RG grapevine virus B) and a genome for which no reference is available (RG grapevine virus L). A large number of SNPs detected in this study ascertained the quasispecies nature of viruses. Detection of three recombination events and phylogenetic analyses using reconstructed genomes suggested the possible introduction of viruses and viroids into India from several continents through the planting material. The whole genome sequences generated in this study can serve as a resource for reliable indexing of grapevine viruses and viroids in quarantine stations and certification programs.

The SARS-CoV-2 spike (S) protein, the viral mediator for binding and entry into the host cell, has sparked great interest as a target for vaccine development and treatments with neutralizing antibodies. Initial data suggest that the virus has low mutation rates, but its large genome could facilitate recombination, insertions, and deletions, as has been described in other coronaviruses. Here, we deep-sequenced the complete SARS-CoV-2 S gene from 18 patients (10 with mild and 8 with severe COVID-19), and found that the virus accumulates deletions upstream and very close to the S1/S2 cleavage site, generating a frameshift with appearance of a stop codon. These deletions were found in a small percentage of the viral quasispecies (2.2%) in samples from all the mild and only half the severe COVID-19 patients. Our results suggest that the virus may generate free S1 protein released to the circulation. We propose that natural selection has favored a “Don’t burn down the house” strategy, in which free S1 protein may compete with viral particles for the ACE2 receptor, thus reducing the severity of the infection and tissue damage without losing transmission capability.

In the wake of the current SARS-CoV-2 pandemic devastating the world, it is imperative to elucidate the comparative genomics of geographically-diverse strains of this novel coronavirus to gain insights into its microevolution, pathogenesis and control. Here we explore the molecular nature, genome-wide frequency, and gene-wise distribution of mutations in three distinct datasets encompassing 68 SARS-CoV-2 RNA-genomes altogether. While phylogenomic analysis revealed parallelism between the evolutionary paths charted by distinct quasispecies clusters of the virus, occurrence of mutations across genomes was found to be non-random. Whereas deletion mutations are extremely scarce and insertions totally absent, of all the instances of single nucleotide substitution detected, the overwhelming majority were transition mutations with cytidine to uridine being the most prevalent type. Propensity of this transition could be attributed to hydrolytic deamination mediated by ultra-violet irradiation or bisulfite reagent, both of which find wide usage as sterilizer/disinfectant. Transversions, albeit few and predominated by the guanosine to uridine form, were found concentrated in loci encoding the structural proteins of the virus, so might confer versatile tissue-colonization potentials. Mutation frequency of the three distinct genome-sets ranged narrowly between 0.07-1.08 × 10-4 nucleotides positions mutated per nucleotide aligned. Gene-wise mapping of the global mutations illuminated the highly conserved nature of the genes encoding the non-structural proteins Nsp7, Nsp8 (two essential cofactors of the viral RNA-dependent RNA-polymerase) and Nsp9 (Nsp8-interacting single-strand RNA-binding protein), plus the envelope protein E (involved in SARS-CoV-2 assembly, budding and pathogenesis). These mutation-free genomic loci and/or their protein products could be potent targets for future drug designing/targeting.

Viral diversity is an important feature of hepatitis C virus (HCV) infection and an important predictor of disease progression and treatment response. HIV/HCV co-infection is associated with enhanced HCV replication, increased fibrosis, and the development of liver disease. HIV also increases quasispecies diversity of HCV structural genes, although limited data are available regarding the impact of HIV on non-structural genes of HCV, particularly in the absence of direct-acting therapies. The genetic diversity and presence of drug resistance mutations within the RNA-dependent RNA polymerase (NS5B) gene were examined in 3 groups of women with HCV genotype 1a infection, including those with HCV mono-infection, antiretroviral (ART)-naive women with HIV/HCV co-infection and CD4 cell count

One of the features distinguishing SARS-CoV-2 from its more pathogenic counterpart SARS-CoV is the presence of premature stop codons in itsORF3bgene. Here, we show that SARS-CoV-2ORF3bis a potent interferon antagonist, suppressing the induction of type I interferon more efficiently than its SARS-CoV ortholog. Phylogenetic analyses and functional assays revealed that SARS-CoV-2-related viruses from bats and pangolins also encode truncatedORF3bgene products with strong anti-interferon activity. Furthermore, analyses of more than 15,000 SARS-CoV-2 sequences identified a natural variant, in which a longerORF3breading frame was reconstituted. This variant was isolated from two patients with severe disease and further increased the ability of ORF3b to suppress interferon induction. Thus, our findings not only help to explain the poor interferon response in COVID-19 patients, but also describe a possibility of the emergence of natural SARS-CoV-2 quasispecies with extendedORF3bthat may exacerbate COVID-19 symptoms.HighlightsORF3b of SARS-CoV-2 and related bat and pangolin viruses is a potent IFN antagonistSARS-CoV-2 ORF3b suppresses IFN induction more efficiently than SARS-CoV orthologThe anti-IFN activity of ORF3b depends on the length of its C-terminusAn ORF3b with increased IFN antagonism was isolated from two severe COVID-19 cases

Background A virus community is the spectrum of viral strains populating an infected host, which plays a key role in pathogenesis and therapy response in viral infectious diseases. However automatic and dedicated pipeline for interpreting virus community sequencing data has not been developed yet. Results We developed Quasispecies Analysis Package (QAP), an integrated software platform to address the problems associated with making biological interpretations from massive viral population sequencing data. QAP provides quantitative insight into virus ecology by first introducing the definition “virus OTU” and supports a wide range of viral community analyses and results visualizations. Various forms of QAP were developed in consideration of broader users, including a command line, a graphical user interface and a web server. Utilities of QAP were thoroughly evaluated with high-throughput sequencing data from hepatitis B virus, hepatitis C virus, influenza virus and human immunodeficiency virus, and the results showed highly accurate viral quasispecies characteristics related to biological phenotypes. Conclusions QAP provides a complete solution for virus community high throughput sequencing data analysis, and it would facilitate the easy analysis of virus quasispecies in clinical applications.

HIV exists as multiple genotypes in a single infected individual referred to as a quasispecies. Here we reproduced a quasispecies by moderate selective pressure using an HIV reverse transcriptase inhibitor. The drug resistant genotype never completely supplanted the drug sensitive genotype, which stabilized at about 20 percent of viral sequences. Single-cell sequencing showed that resistant genotype frequency plateaued when cells were co-infected with sensitive and resistant genotypes, suggesting a sharing of viral proteins in co-infected cells (complementation) which masks genotypic differences. To test if complementation can confer phenotypic drug resistance, we co-transfected fluorescently labelled molecular clones of sensitive and resistant HIV and observed that genotypically sensitive virus from co-transfected cells was drug resistant. Resistant virus preferentially infected cells in tandem with drug sensitive HIV, explaining how co-infections of sensitive and resistant genotypes were initiated. Modelling this effect, we observed that a stable quasispecies could form at the experimental multiplicities of infection for the drug resistant and drug sensitive virus, showing that complementation can lead to a quasispecies in an HIV evolution experiment.

Background Notwithstanding the efforts of direct-acting antivirals (DAAs) for the treatment of chronically infected hepatitis C virus (HCV) patients, concerns exist regarding the emergence of resistance-associated substitutions (RAS) related to therapy failure. Sanger sequencing is still the reference technique used for the detection of RAS and it detects viral variants present up to 15%, meaning that minority variants are undetectable, using this technique. To date, many studies are focused on the analysis of the impact of HCV low variants using next-generation sequencing (NGS) techniques, but the importance of these minority variants is still debated, and importantly, a common data analysis method is still not defined. Methods Serum samples from four patients failing DAAs therapy were collected at baseline and failure, and amplification of NS3, NS5A and NS5B genes was performed on each sample. The genes amplified were sequenced using Sanger and NGS Illumina sequencing and the data generated were analyzed with different approaches. Three different NGS data analysis methods, two homemade in silico pipeline and one commercially available certified user-friendly software, were used to detect low-level variants. Results The NGS approach allowed to infer also very-low level virus variants. Moreover, data processing allowed to generate high accuracy data which results in reduction in the error rates for each single sequence polymorphism. The results improved the detection of low-level viral variants in the HCV quasispecies of the analyzed patients, and in one patient a low-level RAS related to treatment failure was identified. Importantly, the results obtained from only two out of the three data analysis strategies were in complete agreement in terms of both detection and frequency of RAS. Conclusions These results highlight the need to find a robust NGS data analysis method to standardize NGS results for a better comprehension of the clinical role of low-level HCV variants. Based on the extreme importance of data analysis approaches for wet-data interpretation, a detailed description of the used pipelines and further standardization of the in silico analysis could allow increasing diagnostic laboratory networking to unleash true potentials of NGS.

The trimeric envelope glycoprotein (gp120/gp41)3 of human immunodeficiency virus-1 (HIV-1) mediates viral and host cell membrane fusion, initiated by binding of viral envelope gp120 protein to the CD4 receptor on host immune cells. Functional env genes from infected individuals have been widely used as templates for vaccine design, for setting up viral neutralization assays and to study the viral evolution and pathogenesis. Traditional topoisomerase or T4 DNA polymerase mediated approaches for cloning single genome amplified (SGA) env genes are labor-intensive, cost-ineffective with low-throughput, thereby enabling functional analysis of only a limited number of env genes from the diverse circulating quasispecies in infected individuals. Herein, we report an efficient, easy to optimize and high-throughput approach for cloning diverse HIV-1 env genes. Multiple env/rev gene cassettes, derived from infected infants, were subjected to SGA using Phusion polymerase and utilized as megaprimers in overlap extension PCR mediated cloning (OEC), circumventing the requirement for novel enzymes. Furthermore, utilization of Phusion polymerase for both the amplification of env/rev cassettes and OEC allows convenient monitoring and optimization, thereby providing much greater flexibility and versatility for analysis of env genes from HIV-1 infected individuals.

Crimean-Congo hemorrhagic fever virus (CCHFV) is a tick-borne virus that can cause a hemorrhagic fever in humans, with a case fatality rate of up to 40%. Cases of CCHFV have been reported in Africa, Asia, and southern Europe; and recently, due to the expanding range of its vector, autochthonous cases have been reported in Spain. Although it was discovered over 70 years ago, our understanding of the pathogenesis of this virus remains limited. We used RNA-Seq in two human liver cell lines (HepG2 and Huh7) infected with CCHFV (strain IbAr10200), to examine kinetic changes in host expression and viral replication simultaneously at 1 and 3 days post infection. Through this, numerous host pathways were identified that were modulated by the virus including: antiviral response and endothelial cell leakage. Notably, the genes encoding DDX60, a cytosolic component of the RIG-I signalling pathway and OAS2 were both shown to be dysregulated. Interestingly, PTPRR was induced in Huh7 cells but not HepG2 cells. This has been associated with the TLR9 signalling cascade, and polymorphisms in TLR9 have been associated with poor outcomes in patients. Additionally, we performed whole-genome sequencing on CCHFV to assess viral diversity over time, and its relationship to the host response. As a result, we have demonstrated that through next-generation mRNA deep-sequencing it is possible to not only examine mRNA gene expression, but also to examine viral quasispecies and typing of the infecting strain. This demonstrates a proof-of-principle that CCHFV specimens can be analyzed to identify both the virus and host biomarkers that may have implications for prognosis., Author summary Crimean-Congo hemorrhagic fever virus (CCHFV) is an understudied tick-borne virus that can cause a wide spectrum of disease, ranging from moderate to severe, and can be fatal. Cases have been reported in Asia, Africa and Europe, but the range of the vector continues to expand. Currently, our understanding of the host innate immune response to the virus remains limited. Our aim was to use RNA-seq, a type of next generation sequencing technology, to characterize the host immune response in liver cells as well as sequence the genome of the virus. Results identified numerous genes and pathways that were altered and served as proof of principle that the viral identification and evolution could be investigated from the same sample simultaneously. This study highlights the potential for this technique for both characterization of the virus as well as identification of host biomarkers that could be potential predictors of patient outcome. Additionally it may provide important information about pathogenesis prior to performing animal infection studies.

Bovine viral diarrhoea (BVD) is an important disease of cattle with significant impacts on animal health and welfare. The wide host range of the causative pestiviruses may lead to formation of virus reservoirs in other ruminant or wildlife species, presenting a concern for the long-term success of BVD eradication campaigns. It is likely that the quasispecies nature of these RNA viruses contributes to their interspecies transmission by providing genetic plasticity. Understanding the spectrum of sequence variants present in BVD persistently infected (PI) animals is therefore essential for studies of virus transmission. To analyse quasispecies diversity without amplification bias, we extracted viral RNA from serum of a PI cow, and from cell culture fluid after three passages of the same virus in culture, to produce cDNA without amplification. Sequencing of this material using Illumina 250bp paired-read technology produced full-length virus consensus sequences from both sources and demonstrated the quasispecies diversity of this Pestivirus A type 1a field strain within serum and after culture. We report the distribution and diversity of over 800 single nucleotide polymorphisms and provide evidence for a loss of diversity after only three passages in cell culture, implying that cultured viruses cannot be used to understand quasispecies diversity and may not provide reliable molecular markers for source tracing or transmission studies.Additionally, both serum and cultured viruses could be sequenced as a set of 25 overlapping PCR amplicons that demonstrated the same consensus sequences and the presence of many of the same quasispecies variants. The observation that aspects of the quasispecies structure revealed by massively parallel sequencing are also detected after PCR and Sanger sequencing suggests that this approach may be useful for small or difficult to analyse samples.Impact statementBovine viral diarrhoea viruses are globally important cattle pathogens, which impact performance due to acute infection and BVD-induced immunosuppression. Eradication of BVD in cattle is widely pursued but is hampered by the production of persistently infected (PI) calves – the offspring of cows infected in early pregnancy – which shed virus constantly and drive BVD spread. Genetic variation in BVD viruses is an important feature of their biology, allowing them to adapt to changing conditions and to infect different hosts. Inaccurate virus replication produces a population of viruses with slightly different sequences, a quasispecies, some of which may grow better in other hosts or in culture. Analysing virus sequence variation may help us understand how the virus evolves within and between its hosts. In this paper we show that a BVD virus strain loses quasispecies diversity quickly when cultured and that these changes can be detected even in small diagnostic samples, implying that cultured viruses do not perfectly represent the field strains they were isolated from and therefore may not provide reliable molecular markers for source tracing or transmission studies.Data SummaryPestivirus A genome sequences used in this article are as follows:Sequence data associated with this manuscript has been submitted to the European Nucleotide Archive (www.ebi.ac.uk/ena/) with accession numbers as follows:Consensus genome sequences:MRI103 serum NGS: LR699799MRI103 culture NGS: LR699800MRI103 serum Sanger: LR699801MRI103 culture P3 Sanger: LR699802MRI103 culture P5 Sanger: LR699803NGS raw dataSerum dataset: ERR3624580Culture dataset: ERR3624581

Objective Pretreatment HIV-drug resistance (PDR, HIVDR) to non-nucleoside reverse transcriptase inhibitors (NNRTIs) is increasing globally. NNRTIs continue to be used as first-line antiretroviral therapy (ART) in some communities due to the cost of dolutegravir-based ART or dolutegravir-associated adverse events. A simplified version of the oligonucleotide ligation assay (OLA) - 'OLA-Simple' - is a low-cost, near point-of-care assay that provides ready-to-use lyophilized reagents and reports HIVDR mutations as colored lines on lateral flow strips. Our objective was to design and validate OLA-Simple for a Mexican cohort. Design OLA-Simple probes to detect K65R, K103N/S, Y181C, M184V, and G190A were optimized for HIV Mexican sequences. Sixty clinical plasma specimens were analyzed by OLA-Simple by technicians blinded to Illumina-MiSeq sequences, and HIVDR results were compared. Methods Plasma RNA was tested using OLA-Simple kits. OLA-Simple lateral flow strips were read by in-house software, and were classified as mutant or wild-type at each codon. The comparison of results by OLA-Simple and Miseq was used to generate receiver-operating characteristic curves. Results OLA-Simple PCR amplified 59 of 60 specimens and successfully genotyped 287 of 295 codons, with eight of 295 (2.7%) indeterminate results. Compared to MiSeq, OLA-Simple gave five of 295 (1.7%) false-positive and four of 295 (1.4%) false-negative results. Excluding indeterminate results, OLA-Simple classified mutant with an accuracy of 97.4 and 98.8% when using thresholds at 10 and 25% mutant within an individual's HIV quasispecies, respectively. Conclusions Compared to MiSeq, OLA-Simple detected HIVDR with high sensitivity and accuracy. OLA-Simple could expand access to affordable and rapid HIVDR testing to guide appropriate ART choices in populations using NNRTI-based ART.

We present the first genetic characterization of a COVID-19 cluster in Europe using metagenomic next-generation sequencing (mNGS). Despite low viral loads, the mNGS workflow used herein allowed to characterize the whole genome sequences of SARS-CoV2 isolated from an asymptomatic patient, in 2 clinical samples collected 1 day apart. Comparison of these sequences suggests viral evolution with development of quasispecies. In addition, the present workflow identified a new deletion in nsp2 (Asp268Del) which was found in all 3 samples originating from this cluster as well as in 37 other viruses collected in England and in Netherlands, suggesting the spread of this deletion in Europe. The impact of Asp268Del on SARS-CoV-2 transmission and pathogenicity, as well as on PCR performances and anti-viral strategy should be rapidly evaluated in further studies.

H9N2 avian influenza virus has spread worldwide, and vaccination with an inactivated virus is currently the major prevention method in China. To further understand the effect of the selection pressure from antibodies on the evolution of H9N2 avian influenza virus, F/98 (A/Chicken/Shanghai/F/98), which is the vaccine representative of H9N2 avian influenza virus in East China, was used for serial passaging for 20 generations in chickens with and without vaccination. After plaque purification from trachea and lung tissues, 390 quasispecies were obtained. The second-generation quasispecies under the selection pressure of vaccine antibodies had undergone 100% antigen variation, while after passaging to the fifth generation, only 30–40% of the quasispecies displayed antigen variation when there was no selection pressure of vaccine antibodies, implying that the selection pressure of vaccine antibodies promotes the antigen variation of F/98. We found for the first time that there were three mutation hotspots in the HA genes of the quasispecies under the selection pressure of vaccine antibodies, which were K131R, A168T, and N201D. Moreover, under the selection pressure of vaccine antibodies, 10 amino acids (67–76) of the NA protein of all quasispecies were deleted, and PB2 of the quasispecies had undergone a high-frequency R355K mutation. However, without selection pressure of vaccine antibodies, NP had undergone two high-frequency mutations, namely, V186I and L466I, and a high-frequency mutation of L77I appeared in the NS gene. This result shows that the vaccine antibody selection pressure could control and regulate gene variation of the F/98 virus. Compared to that of the parental virus F/98, the EID50 of the twentieth passaged virus under the selection pressure of vaccine antibodies did not change, while the EID50 of the twentieth passaged virus without selection pressure of vaccine antibodies was significantly enhanced by 794 times. Furthermore, the twentieth passaged virus with selection pressure from vaccine antibodies lost its lethal ability in embryonated chicken eggs, whereas the EID50 of the twentieth passaged virus without selection pressure of vaccine antibodies increased to 6.3 times that of the F/98 strain. All the above results show that the selection pressure of vaccine antibodies promotes the antigen variation of H9N2 avian influenza virus and plays a role in regulating and controlling gene mutation of H9N2 avian influenza virus.

Currently, 50% of the adult population worldwide is overweight or obese. In these studies, we demonstrate that obesity not only enhances the severity of influenza infection but also impacts viral diversity. The altered microenvironment associated with obesity supports a more diverse viral quasispecies and affords the emergence of potentially pathogenic variants capable of inducing greater disease severity in lean hosts. This is likely due to the impaired interferon response, which is seen in both obese mice and obesity-derived human bronchial epithelial cells, suggesting that obesity, aside from its impact on influenza virus pathogenesis, permits the stochastic accumulation of potentially pathogenic viral variants, raising concerns about its public health impact as the prevalence of obesity continues to rise., Obesity is associated with increased disease severity, elevated viral titers in exhaled breath, and significantly prolonged viral shed during influenza A virus infection. Due to the mutable nature of RNA viruses, we questioned whether obesity could also influence influenza virus population diversity. Here, we show that minor variants rapidly emerge in obese mice. The variants exhibit increased viral replication, resulting in enhanced virulence in wild-type mice. The increased diversity of the viral population correlated with decreased type I interferon responses, and treatment of obese mice with recombinant interferon reduced viral diversity, suggesting that the delayed antiviral response exhibited in obesity permits the emergence of a more virulent influenza virus population. This is not unique to obese mice. Obesity-derived normal human bronchial epithelial (NHBE) cells also showed decreased interferon responses and increased viral replication, suggesting that viral diversity also was impacted in this increasing population.

While direct-acting antivirals (DAAs) for hepatitis C virus (HCV) have dramatically progressed, patients still suffer from treatment failures. For the radical eradication of HCV, a deeper understanding of multiple resistance-associated substitutions (RASs) at the single-clone level is essential. To understand HCV quasispecies and their dynamics during DAA treatment, we applied single-molecule real-time (SMRT) deep sequencing on sera from 12 patients with genotype-1b HCV infections with DAA treatment failures, both pre- and post-treatment. We identified >3.2 kbp sequences between NS3 and NS5A genes of 187,539 clones in total, classifying into haplotype codes based on the linkage of seven RAS loci. The number of haplotype codes during the treatment, per sample, significantly decreased from 14.67 ± 9.12 to 6.58 ± 7.1, while the number of nonsynonymous codons on the seven RAS loci, per clone, significantly increased from 1.50 ± 0.92 to 3.64 ± 0.75. In five cases, the minority multi-drug resistant haplotypes at pre-treatment were identical to the major haplotypes at relapse. Moreover, various structural variations (SVs) were detected and their dynamics analysed. These results suggest that SMRT deep sequencing is useful for detecting minority haplotypes and SVs, and to evaluate the dynamics of viral genomes at the single-clone level.

Oncolytic virotherapy is an emerging treatment option for numerous cancers, with several virus families currently being evaluated in clinical trials. More specifically, vaccine-strain measles virus has arisen as a promising candidate for the treatment of different tumour types in several early clinical trials. Replicating viruses, and especially RNA viruses without proofreading polymerases, can rapidly adapt to varying environments by selecting quasispecies with advantageous genetic mutations. Subsequently, these genetic alterations could potentially weaken the safety profile of virotherapy. In this study, we demonstrate that, following an extended period of virus replication in producer or cancer cell lines, the quasispecies consensus sequence of vaccine strain-derived measles virus accrues a remarkably small number of mutations throughout the nonsegmented negative-stranded RNA genome. Interestingly, we detected a nonrandom distribution of genetic alterations within the genome, with an overall decreasing frequency of mutations from the 3' genome start to its 5' end. Comparing the serially passaged viruses to the parental virus on producer cells, we found that the acquired consensus mutations did not drastically change viral replication kinetics or cytolytic potency. Collectively, our data corroborate the genomic stability and excellent safety profile of oncolytic measles virus, thus supporting its continued development and clinical translation as a promising viro-immunotherapeutic.

International audience; Recombinant strains of hepatitis E virus (HEV) with insertions of human genomic fragments or HEV sequence duplications in the sequence encoding the polyproline region (PPR) were previously described in chronically infected patients. Such genomic rearrangements confer a replicative advantage in vitro but little is known about their frequency, location, or origin. As the sequences of only a few virus genomes are available, we analyzed the complete genomes of 114 HEV genotype 3 strains from immunocompromised (n = 85) and immunocompetent (n = 29) patients using the single molecular real-time sequencing method to determine the frequency, location, and origin of inserted genomic fragments, plus the proportions of variants with genomic rearrangements in each virus quasispecies. We also examined the amino acid compositions and post-translational modifications conferred by these rearrangements by comparing them to sequences without human gene insertions or HEV gene duplications. We found genomic rearrangements in 7/114 (6.1%) complete genome sequences (4 HEV-3f, 1 HEV-3e, 1 HEV-3 h, and 1 HEV-3chi-new), all from immunocompromised patients, and 3/7 were found at the acute phase of infection. Six of the seven patients harbored virus-host recombinant variants, including one patient with two different recombinant variants. We also detected three recombinant variants with genome duplications of the PPR or PPR + X domains in a single patient. All the genomic rearrangements (seven human fragment insertions of varying origins and three HEV genome duplications) occurred in the PPR. The sequences with genomic rearrangements had specific characteristics: increased net load (p < 0.001) and more ubiquitination (p < 0.001), phosphorylation (p < 0.001), and acetylation (p < 0.001) sites. The human fragment insertions and HEV genome duplications had slightly different characteristics. We believe this is the first description of HEV strains with genomic rearrangements in patients at the acute phase of infection; perhaps these strains are directly transmitted. Clearly, genomic rearrangements produce a greater net load with duplications and insertions having different features. Further studies are needed to clarify the mechanisms by which such modifications influence HEV replication.

To gain insight into the impact of mutations on the viability of the hepatitis C virus (HCV) genome, we created a set of full-genome mutant libraries, differing from the parent sequence as well as each other, by using a random mutagenesis approach; the proportion of mutations increased across these libraries with declining template amount or dATP concentration. The replication efficiencies of full-genome mutant libraries ranged between 71 and 329 focus-forming units (FFU) per 10(5) Huh7.5 cells. Mutant libraries with low proportions of mutations demonstrated low replication capabilities, whereas those with high proportions of mutations had their replication capabilities restored. Hepatoma cells transfected with selected mutant libraries, with low (4 mutations per 10,000 bp copied), moderate (33 mutations), and high (66 mutations) proportions of mutations, and their progeny were subjected to serial passage. Predominant virus variants (mutants) from these mutant libraries (Mutant(l), Mutant(m), and Mutant(h), respectively) were evaluated for changes in growth kinetics and particle-to-FFU unit ratio, virus protein expression, and modulation of host cell protein synthesis. Mutant(m) and Mutant(l) variants produced >3.0-log-higher extracellular progeny per ml than the parent, and Mutant(h) produced progeny at a rate 1.0-log lower. More than 80% of the mutations were in a nonstructural part of the mutant genomes, the majority were nonsynonymous, and a moderate to large proportion were in the conserved regions. Our results suggest that the HCV genome has the ability to overcome lethal/deleterious mutations because of the high reproduction rate but highly selects for random, beneficial mutations. IMPORTANCE Hepatitis C virus (HCV) in vivo displays high genetic heterogeneity, which is partly due to the high reproduction and random substitutions during error-prone genome replication. It is difficult to introduce random substitutions in vitro because of limitations in inducing mutagenesis from the 5′ end to the 3′ end of the genome. Our study has overcome this limitation. We synthesized full-length genomes with few to several random mutations in the background of an HCV clone that can recapitulate all steps of the life cycle. Our study provides evidence of the capability of the HCV genome to overcome deleterious mutations and remain viable. Mutants that emerged from the libraries had diverse phenotype profiles compared to the parent, and putative adaptive mutations mapped to segments of the conserved nonstructural genome. We demonstrate the potential utility of our system for the study of sequence variation that ensures the survival and adaptation of HCV.

Darwinian principles guide the evolution of viruses. Several interconnected mathematical equations describe general evolutionary dynamics . One of them, the quasispecies equation emphasizes the occurrence of mutation in the course of genome replication and, therefore, it is adequate as a theoretical framework to understand the behavior of viruses displaying error-prone replication. Viral quasispecies are dynamic collections of related variants (termed mutant spectra, clouds or swarms) subjected to a continuous process of genetic variation, competition and selection in favor of those genomes that are more fit in the environment where replication takes place. Selection and random drift shape virus evolution that can be viewed as an exploration of an immense sequence space that results in the replacement of some genome subpopulations by others. Internal interactions among components of a mutant spectrum can be of complementation or interference, and they may modulate the behavior of the ensemble. Viruses exploit the dual potential offered by individual and group selection to

Experimental evolution permits exploring the effect of controlled environmental variables in virus evolution. Several designs in cell culture and in vivo have established basic concepts that can assist in the interpretation of evolutionary events in the field. Important information has come from cytolytic and persistent infections in cell culture that have unveiled the power of virus-cell coevolution in virus and cell diversification. Equally informative are comparisons of the response of viral populations when subjected to different passage regimes. In particular, plaque-to-plaque transfers in cell culture have revealed unusual genotypes and phenotypes that populate minority layers of viral quasispecies. Some of these viruses display properties that contradict features established in virology textbooks. Several hypotheses and principles of population genetics have found experimental confirmation in experimental designs with viruses. The possibilities of using experimental evolution to understand virus behavior are still largely unexploited.

International audience; Prion replication results from the autocatalytic templated assisted conversion of the host-encoded prion protein PrPC into misfolded, polydisperse PrPSc conformers. Structurally distinct PrPSc conformers can give rise to multiple prion strains. Within and between prion strains, the biological activity (replicative efficacy and specific infectivity) of PrPSc assemblies is size dependent and thus reflects an intrinsic structural heterogeneity. The contribution of such PrPSc heterogeneity across species prion adaptation, which is believed to be based on fit adjustment between PrPSc template(s) and host PrPC, has not been explored. To define the structural-to-fitness PrPSc landscape, we measured the relative capacity of size-fractionated PrPSc assemblies from different prion strains to cross mounting species barriers in transgenic mice expressing foreign PrPC. In the absence of a transmission barrier, the relative efficacy of the isolated PrPSc assemblies to induce the disease is like the efficacy observed in the homotypic context. However, in the presence of a transmission barrier, size fractionation overtly delays and even abrogates prion pathogenesis in both the brain and spleen tissues, independently of the infectivity load of the isolated assemblies. Altering by serial dilution PrPSc assembly content of non-fractionated inocula aberrantly reduces their specific infectivity, solely in the presence of a transmission barrier. This suggests that synergy between structurally distinct PrPSc assemblies in the inoculum is requested for crossing the species barrier. Our data support a mechanism whereby overcoming prion species barrier requires complementation between structurally distinct PrPSc assemblies. This work provides key insight into the “quasispecies” concept applied to prions, which would not necessarily rely on prion substrains as constituent but on structural PrPSc heterogeneity within prion population.

Newcastle disease virus (NDV) infections are well known to harbour quasispecies, due to the error-prone nature of the RNA polymerase. Quasispecies variants in the fusion cleavage site of the virus are known to significantly change its virulence. However, little is known about the genomic patterns of diversity and selection in NDV viral swarms. We analyse deep sequencing data from in vitro and in vivo NDV infections to uncover the genomic patterns of diversity and the signatures of selection within NDV swarms. Variants in viruses from in vitro samples are mostly localised in non-coding regions and 3&prime, and 5&prime, untranslated regions (3&prime, UTRs or 5&prime, UTRs), while in vivo samples contain an order of magnitude more variants. We find different patterns of genomic divergence and diversity among NDV genotypes, as well as differences in the genomic distribution of intra-host variants among in vitro and in vivo infections of the same strain. The frequency spectrum shows clear signatures of intra-host purifying selection in vivo on the matrix protein (M) coding gene and positive or diversifying selection on nucleocapsid (NP) and haemagglutinin-neuraminidase (HN). The comparison between within-host polymorphisms and phylogenetic divergence reveals complex patterns of selective pressure on the NDV genome at between- and within-host level. The M sequence is strongly constrained both between and within hosts, fusion protein (F) coding gene is under intra-host positive selection, and NP and HN show contrasting patterns: HN RNA sequence is positively selected between hosts while its protein sequence is positively selected within hosts, and NP is under intra-host positive selection at the RNA level and negative selection at the protein level.

Consider a mathematical model of evolutionary adaptation of fitness landscape and mutation matrix as a reaction to population changes. As a basis, we use an open quasispecies model, which is modified to include explicit death flow. We assume that evolutionary parameters of mutation and selection processes vary in a way to maximize the mean fitness of the system. From this standpoint, Fisher's theorem of natural selection is being rethought and discussed. Another assumption is that system dynamics has two significant timescales. According to our central hypothesis, major evolutionary transitions happen in the steady-state of the corresponding dynamical system, so the evolutionary time is much slower than the one of internal dynamics. For the specific cases of quasispecies systems, we show how our premises form the fitness landscape adaptation process.

The SARS-CoV-2 spike (S) protein, the viral mediator for binding and entry into the host cell, has sparked great interest as a target for vaccine development and treatments with neutralizing antibodies. Initial data suggest that the virus has low mutation rates, but its large genome could facilitate recombination, insertions, and deletions, as has been described in other coronaviruses. Here, we deep-sequenced the complete SARS-CoV-2 S gene from 18 patients (10 with mild and 8 with severe COVID-19), and found that the virus accumulates deletions upstream and very close to the S1/S2 cleavage site (PRRAR/S), generating a frameshift with appearance of a stop codon. These deletions were found in a small percentage of the viral quasispecies (2.2%) in samples from all the mild and only half the severe COVID-19 patients. Our results suggest that the virus may generate free S1 protein released to the circulation. We suggest that natural selection has favoured a “Don’t burn down the house” strategy, in which free S1 protein may compete with viral particles for the ACE2 receptor, thus reducing the severity of the infection and tissue damage without losing transmission capability.

The main biological threat to the western honeybee (Apis mellifera) is the parasitic mite Varroa destructor, largely because it vectors lethal epidemics of honeybee viruses that, in the absence of this mite, are relatively innocuous. The severe pathology is a direct consequence of excessive virus titres caused by this novel transmission route. However, little is known about how the virus adapts genetically during transmission and whether this influences the pathology. Here, we show that upon injection into honeybee pupae, the deformed wing virus type-A (DWV-A) quasispecies undergoes a rapid, extensive expansion of its sequence space, followed by strong negative selection towards a uniform, common shape by the time the pupae have completed their development, with no difference between symptomatic and asymptomatic adults in either DWV titre or genetic composition. This suggests that the physiological and molecular environment during pupal development has a strong, conservative influence on shaping the DWV-A quasispecies in emerging adults. There was furthermore no evidence of any progressive adaptation of the DWV-A quasispecies to serial intra-abdominal injection, simulating mite transmission, despite the generation of ample variation immediately following each transmission, suggesting that the virus either had already adapted to transmission by injection, or was unaffected by it. Many thanks to the Uppsala Genome Centre for sequencing services, to Kaspar Roth at the Institute of Bee Health (IBH) for technical support, to Pilar de la Rua for analytical advice and to the four excellent anonymous reviewers for constructive comments on the manuscript. Bioinformatic analyses were performed at the Uppsala Multidisciplinary Center for Advanced Computational Science (UPPMAX). This project was funded by grant C09.0048 from the EU-COST FA0803 action, grant SNF 31003A_141006/1 from the Swiss National Science Foundation and the Vinetum Foundation. Open access funding provided by Swedish University of Agricultural Sciences. info:eu-repo/semantics/publishedVersion

Avian coronavirus (AvCoV) infects a range of tissues in chickens and several other avian species. Although the viruscan be isolated in chicken embryos, only a few strains of the 6 genotypes/33 lineages can grow in cell lines, with theBeaudette strain (GI-1 lineage) being the most used for in vitro studies. Considering the differences between celllines and chicken embryos as habitats for AvCoV, this study aimed to assess the diversity of the genes coding forthe nonstructural protein 3 (nsp3) and spike envelope protein (S) after serial passages in BHK-21 and Vero cells.After 14 passages of an embryo-adapted Beaudette strain, the virus loads fluctuated in both cell lines, with the highestloads being 8.72 log genome copies/µL for Vero and 6.36 log genome copies/µL for BHK-21 cells. No polymorphismswere found for nsp3; regarding S, not only aa substitutions (Vero: 8th passage A150S, and 14th S150A; BHK-21: 4thS53F, 8th F53Y, and 8th S95R), but also minor variants could be detected on chromatograms with fluctuating intensities.As the regions of these aa substitutions are within the receptor-binding domain of S, it can be speculated that differencesin cell receptors between Vero and BHK-21 cells and the speed of cell death led to the selection of different dominantstrains, while the stability of nsp3 supports its function as a protease involved in AvCoV replication. In conclusion,AvCoV quasispecies evolution is influenced by the biological model under consideration, and a gradual transition isseen for minor and major variants. O Coronavírus aviário AvCoV infecta uma variedade de tecidos de galinhas e de outras espécies aviárias. Apesarde este vírus poder ser isolado em ovos embrionados de galinha, apenas alguns dos 6 genótipos / 33 linhagenspodem crescer em cultivo celular, sendo a cepa Beuadette (linhagem GI-11) a mais utilizada para estudos in vitro.Considerando as diferentes linhagens celulares e ovos embrionados como habitats para o AvCoV, este estudo tevepor objetivo estudar a diversidade de genes que codificam para a proteína não-estrutural 3 (nsp3) e espícula (S)após passagens seriadas em células BHK-21 e VERO. Após 14 passagens, de uma amostra Beuadette adaptada aovos embrionados, os títulos virais variaram em ambas as células, com os maiores títulos sendo de 8,72 log cópiasgenômicas/µL para Vero e 6,36 cópias genômicas/µL para BHK-21. Nenhum polimorfismo foi encontrando paransp3. Considerando a proteína S, não somente foram encontradas substituições de aminoácidos (Vero: 8apassagemA150S e 14apassagem S150A; BHK-21: 4apassagem S53F, 8apassagem F53Y e S95R), mas também, variantessubconsensuais foram detectadas pelos cromatogramas com intensidades flutuantes. Uma vez que as regiões destesaa se encontram no domínio de ligação de receptor de S, pode-se especular que diferenças em receptores celularesentre Vero e BHK-21, além da velocidade da morte celular, levaram à seleção de diferentes cepas dominantes,enquanto que a estabilidade de nsp3 concorda com sua função como protease com papel na replicação de AvCoV.Como conclusão, a evolução de quase-espécies de AvCoV é influenciada pelo modelo biológico sob consideraçãoe uma transição gradual é vista para variantes dominantes e subdominantes.&nbsp