Viruses as quasispecies: biological implications

The molecular epidemiology of a dengue virus outbreak in Taiwan: population wide versus infrapopulation mutation analysis

Dengue virus (DENV) causes approximately 390 million dengue infections worldwide every year. There were 22,777 reported DENV infections in Tainan, Taiwan in 2015. In this study, we sequenced the C-prM-E genes from 45 DENV 2015 strains, and phylogenetic analysis based on C-prM-E genes revealed that all strains were classified as DENV serotype 2 Cosmopolitan genotype. Sequence analysis comparing different DENV-2 genotypes and Cosmopolitan DENV-2 sequences prior to 2015 showed a clade replacement event in the DENV-2 Cosmopolitan genotype. Additionally, a major substitution C-A314G (K73R) was found in the capsid region which may have contributed to the clade replacement event. Reverse genetics virus rgC-A314G (K73R) showed slower replication in BHK-21 and C6/36 cells compared to wildtype virus, as well as a decrease in NS1 production in BHK-21-infected cells. After a series of passaging, the C-A314G (K73R) mutation reverted to wildtype and was thus considered to be unstable. Next generation sequencing (NGS) of three sera collected from a single DENV2-infected patient at 1-, 2-, and 5-days post-admission was employed to examine the genetic diversity over-time and mutations that may work in conjunction with C-A314G (K73R). Results showed that the number of haplotypes decreased with time in the DENV-infected patient. On the fifth day after admission, two new haplotypes emerged, and a single non-synonymous NS4A-L115I mutation was identified. Therefore, we have identified a persistent mutation C-A314G (K73R) in all of the DENV-2 isolates, and during the course of an infection, a single new non-synonymous mutation in the NS4A region appears in the virus population within a single host. The C-A314G (K73R) thus may have played a role in the DENV-2 2015 outbreak while the NS4A-L115I may be advantageous during DENV infection within the host.

Cell-cell communication and initial population composition shape the structure of potato spindle tuber viroid quasispecies

RNA viruses and viroids replicate with high mutation rates, forming quasispecies, population of variants centered around dominant sequences. The mechanisms governing quasispecies remain unclear. Plasmodesmata regulate viroid movement and were hypothesized to impact viroid quasispecies. Here, we sequenced the progeny of potato spindle tuber viroid intermediate (PSTVd-I) strain from mature guard cells lacking plasmodesmal connections and from in vitro-cultivated mesophyll cell protoplasts from systemic leaves of early-infected tomato (Solanum lycopersicum) plants. Remarkably, more variants accumulated in guard cells compared to whole leaves. Similarly, after extended cell culture, we observed more variants in cultivated mesophyll protoplasts. Coinfection and single-cell sequencing experiments demonstrated that the same plant cell can be infected multiple times by the same or different PSTVd sequences. To study the impact of initial population composition on PSTVd-I quasispecies, we conducted coinfections with PSTVd-I and variants. Two inoculum ratios (10:1 or 1:10) established quasispecies with or without PSTVd-I as the master sequence. In the absence of the master sequence, the percentage of novel variants initially increased. Moreover, a 1:1 PSTVd-I/variant RNA ratio resulted in PSTVd-I dominating (>50%), while the variants reached 20%. After PSTVd-I-only infection, the variants reached around 10%, while after variant-only infection, the variants were significantly more than 10%. These results emphasize the role of cell-to-cell communication and initial population composition in shaping PSTVd quasispecies.

The dynamic variation position and predominant quasispecies of hepatitis B virus: Novel predictors of early hepatocarcinoma

To find the predictors of early HCC based on the dynamic changes of HBV quasispecies, this study utilizing the second-generation sequencing (NGS) and high-order multiplex droplet digital PCR (ddPCR) technology to examine the HBV quasispecies in serum of total 247 subjects recruited from high-incidence area of HCC. In the discovery stage, 15 non-synonymous Single Nucleotide Polymorphisms (SNPs) with higher variant proportion in HCC case group were founded (all P<0.05). Furthermore, the variant proportions in some of these SNPs were observed changing regularly within 5 years before the onset of HCC, and 5 of them located in HBX, 2 in HBS and 2 in HBC. The HBV predominant quasispecies and their consensus sequences were identified by genetic evolution analysis, in which the high HBS and HBC quasispecies heterogeneity were found associated with the forming of multifarious quasispecies clones, and the HBX gene had the highest proportion of predominant quasispecies (46.7 % in HBX vs 12.7 % and 13.8 % in HBS and HBC respectively) with the key variations (G1512A, A1630G, T1753C/G/A, A1762T and G1764A) determined. In the validation stage, we confirmed that the combined double mutations of G1512A+A1630G, A1762T+G1764A, and the combined triple mutations of T1753C/G/A + A1762T+G1764A, all expressed higher in early HCC cases when comparing with control group (all P<0.05). We also demonstrated the advantages of ddPCR using in multi-variations detection in large-sample for early HCC surveillance and screening. So we think that the dynamic of key HBV variation positions and their different combinations determined by quasispecies anlysis in this study can act as the novel predictors of early hepatocarcinoma and suitable to popularize and apply in HCC screening.

The central nervous system is a potential reservoir and possible origin of drug resistance in hepatitis B infection

Background

The significance of hepatitis B virus (HBV) in cerebrospinal fluid (CSF) is unclear.

Methods

Synchronous serum and CSF samples were collected from 13 patients. HBV DNA, full-length genome, quasispecies, phylogenetic tree, compartmentalization and mutation of the reverse transcriptase (RT) region were performed based on PCR and sequencing methods.

Results

HBV DNA was detected in the CSF of 3 antiviral-naïve individuals and 1 individual after successful antiviral therapy. Complete full-length HBV genomes were isolated from the CSF of 5 individuals, including 2 with undetectable serum HBV DNA. Ten individuals exhibited distinct CSF-serum quasispecies, 8 harbored independent CSF-serum genetic compartmentalization and phylogenetic trees, and 5 lamivudine/entecavir-associated resistance mutations only in the CSF. The frequencies of rtL180M and rtM204I/V mutations in both serum and CSF were higher in HIV-HBV-coinfected individuals than in the HBV-monoinfected ones (serum: rtL180M: 3.9% vs. 0, P = 0.004; rtM204I/V: 21.3% vs. 0, P < 0.001; CSF: rtL180M: 7.6% vs. 0, P = 0.026; rtM204I/V 7.6% vs. 1.6%, P = 0.097).

Conclusion

CSF is a potential HBV reservoir, and HBV in CSF harbors distinct evolution and mutation characteristics from those in serum. HIV infection increases the possibility of HBV rtL180M and rtM204I/V mutations in both serum and CSF.

Tautomerism and Rotamerism of Favipiravir and Halogenated Analogues in Solution and in the Solid State

Favipiravir is an important selective antiviral against RNA-based viruses, and currently, it is being repurposed as a potential drug for the treatment of COVID-19. This type of chemical system presents different carboxamide-rotameric and hydroxyl-tautomeric states, which could be essential for interpreting its selective antiviral activity. Herein, the tautomeric 3-hydroxypyrazine/3-pyrazinone pair of favipiravir and its 6-substituted analogues, 6-Cl, 6-Br, 6-I, and 6-H, were fully investigated in solution and in the solid state through ultraviolet-visible, 1H nuclear magnetic resonance, infrared spectroscopy, and X-ray diffraction techniques. Also, a study of the gas phase was performed using density functional theory calculations. In general, the keto-enol balance in these 3-hydroxy-2-pyrazinecarboxamides is finely modulated by external and internal electrical variations via changes in solvent polarity or by replacement of substituents at position 6. The enol tautomer was prevalent in an apolar environment, whereas an increase in the level of the keto tautomer was favored by an increase in solvent polarity and, even moreso, with a strong hydrogen-donor solvent. Keto tautomerization was favored either in solution or in the solid state with a decrease in 6-substituent electronegativity as follows: H ≫ I ≈ Br > Cl ≥ F. Specific rotameric states based on carboxamide, "cisoide" and "transoide", were identified for the enol and keto tautomer, respectively; their rotamerism is dependent on the tautomerism and not the aggregation state.

The genetic variability and evolution of red-spotted grouper nervous necrosis virus quasispecies can be associated with its virulence

Nervous necrosis virus, NNV, is a neurotropic virus that causes viral nervous necrosis disease in a wide range of fish species, including European sea bass (Dicentrarchus labrax). NNV has a bisegmented (+) ssRNA genome consisting of RNA1, which encodes the RNA polymerase, and RNA2, encoding the capsid protein. The most prevalent NNV species in sea bass is red-spotted grouper nervous necrosis virus (RGNNV), causing high mortality in larvae and juveniles. Reverse genetics studies have associated amino acid 270 of the RGNNV capsid protein with RGNNV virulence in sea bass. NNV infection generates quasispecies and reassortants able to adapt to various selective pressures, such as host immune response or switching between host species. To better understand the variability of RGNNV populations and their association with RGNNV virulence, sea bass specimens were infected with two RGNNV recombinant viruses, a wild-type, rDl956, highly virulent to sea bass, and a single-mutant virus, Mut270Dl965, less virulent to this host. Both viral genome segments were quantified in brain by RT-qPCR, and genetic variability of whole-genome quasispecies was studied by Next Generation Sequencing (NGS). Copies of RNA1 and RNA2 in brains of fish infected with the low virulent virus were 1,000-fold lower than those in brains of fish infected with the virulent virus. In addition, differences between the two experimental groups in the Ts/Tv ratio, recombination frequency and genetic heterogeneity of the mutant spectra in the RNA2 segment were found. These results show that the entire quasispecies of a bisegmented RNA virus changes as a consequence of a single point mutation in the consensus sequence of one of its segments. Sea bream (Sparus aurata) is an asymptomatic carrier for RGNNV, thus rDl965 is considered a low-virulence isolate in this species. To assess whether the quasispecies characteristics of rDl965 were conserved in another host showing different susceptibility, juvenile sea bream were infected with rDl965 and analyzed as above described. Interestingly, both viral load and genetic variability of rDl965 in seabream were similar to those of Mut270Dl965 in sea bass. This result suggests that the genetic variability and evolution of RGNNV mutant spectra may be associated with its virulence.

RNA structure-altering mutations underlying positive selection on Spike protein reveal novel putative signatures to trace crossing host-species barriers in Betacoronavirus

Similar to other RNA viruses, the emergence of Betacoronavirus relies on cross-species viral transmission, which requires careful health surveillance monitoring of protein-coding information as well as genome-wide analysis. Although the evolutionary jump from natural reservoirs to humans may be mainly traced-back by studying the effect that hotspot mutations have on viral proteins, it is largely unexplored if other impacts might emerge on the structured RNA genome of Betacoronavirus. In this survey, the protein-coding and viral genome architecture were simultaneously studied to uncover novel insights into cross-species horizontal transmission events. We analysed 1,252,952 viral genomes of SARS-CoV, MERS-CoV, and SARS-CoV-2 distributed across the world in bats, intermediate animals, and humans to build a new landscape of changes in the RNA viral genome. Phylogenetic analyses suggest that bat viruses are the most closely related to the time of most recent common ancestor of Betacoronavirus, and missense mutations in viral proteins, mainly in the S protein S1 subunit: SARS-CoV (G > T; A577S); MERS-CoV (C > T; S746R and C > T; N762A); and SARS-CoV-2 (A > G; D614G) appear to have driven viral diversification. We also found that codon sites under positive selection on S protein overlap with non-compensatory mutations that disrupt secondary RNA structures in the RNA genome complement. These findings provide pivotal factors that might be underlying the eventual jumping the species barrier from bats to intermediate hosts. Lastly, we discovered that nearly half of the Betacoronavirus genomes carry highly conserved RNA structures, and more than 90% of these RNA structures show negative selection signals, suggesting essential functions in the biology of Betacoronavirus that have not been investigated to date. Further research is needed on negatively selected RNA structures to scan for emerging functions like the potential of coding virus-derived small RNAs and to develop new candidate antiviral therapeutic strategies.

Twenty years of evolution and diversification of digitaria streak virus in Digitaria setigera

Within the family Geminiviridae, the emergence of new species results from their high mutation and recombination rates. In this study, we report the variability and evolution of digitaria streak virus (DSV), a mastrevirus isolated in 1986 from the grass Digitaria setigera in an island of the Vanuatu archipelago. Viral DNA of DSV samples was amplified from D. setigera specimens, derived from the naturally infected original plant, which were propagated in different laboratories in France and Italy for more than 20 years. From the consensus sequences, the nucleotide substitution rate was estimated for the period between a sample and the original sequence published in 1987, as well as for the period between samples. In addition, the intra-host genetic complexity and diversity of 8 DSV populations with a total of 165 sequenced haplotypes was characterized. The evolutionary rate of DSV was estimated to be between 1.13 × 10⁻⁴ and 9.87 × 10⁻⁴ substitutions/site/year, within the ranges observed in other single-stranded DNA viruses and RNA viruses. Bioinformatic analyses revealed high variability and heterogeneity in DSV populations, which confirmed that mutant spectra are continuously generated and are organized as quasispecies. The analysis of polymorphisms revealed nucleotide substitution biases in viral genomes towards deamination and oxidation of single-stranded DNA. The differences in variability in each of the genomic regions reflected a dynamic and modular evolution in the mutant spectra that was not reflected in the consensus sequences. Strikingly, the most variable region of the DSV genome, encoding the movement protein, showed rapid fixation of the mutations in the consensus sequence and a concomitant dN/dS ratio of 6.130, which suggests strong positive selection in this region. Phylogenetic analyses revealed a possible divergence in three genetic lineages from the original Vanuatu DSV isolate.

SARS-CoV-2 Variant Identification Using a Genome Tiling Array and Genotyping Probes

With over three million deaths worldwide attributed to the respiratory disease COVID-19 caused by the novel coronavirus SARS-CoV-2, it is essential that continued efforts be made to track the evolution and spread of the virus globally. We previously presented a rapid and cost-effective method to sequence the entire SARS-CoV-2 genome with 95% coverage and 99.9% accuracy. This method is advantageous for identifying and tracking variants in the SARS-CoV-2 genome when compared to traditional short read sequencing methods which can be time consuming and costly. Herein we present the addition of genotyping probes to our DNA chip which target known SARS-CoV-2 variants. The incorporation of the genotyping probe sets along with the advent of a moving average filter have improved our sequencing coverage and accuracy of the SARS-CoV-2 genome.

Multiplex PCR-Based Nanopore Sequencing and Epidemiological Surveillance of Hantaan orthohantavirus in Apodemus agrarius, Republic of Korea

Whole-genome sequencing of infectious agents enables the identification and characterization of emerging viruses. The MinION device is a portable sequencer that allows real-time sequencing in fields or hospitals. Hantaan orthohantavirus (Hantaan virus, HTNV), harbored by Apodemus agrarius, causes hemorrhagic fever with renal syndrome (HFRS) and poses a critical public health threat worldwide. In this study, we aimed to evaluate the feasibility of using nanopore sequencing for whole-genome sequencing of HTNV from samples having different viral copy numbers. Amplicon-based next-generation sequencing was performed in A. agrarius lung tissues collected from the Republic of Korea. Genomic sequences of HTNV were analyzed based on the viral RNA copy numbers. Amplicon-based nanopore sequencing provided nearly full-length genomic sequences of HTNV and showed sufficient read depth for phylogenetic analysis after 8 h of sequencing. The average identity of the HTNV genome sequences for the nanopore sequencer compared to those of generated from Illumina MiSeq revealed 99.8% (L and M segments) and 99.7% (S segment) identities, respectively. This study highlights the potential of the portable nanopore sequencer for rapid generation of accurate genomic sequences of HTNV for quicker decision making in point-of-care testing of HFRS patients during a hantavirus outbreak.

High Prevalence of Preexisting HBV Polymerase Mutations in Pregnant Women Does Not Limit the Antiviral Therapy Efficacy

Background

HBV-resistant mutants in treatment-naïve patients may lead to antiviral treatment failure. It is not clear if HBV mutants are present in pregnant women and about the influence of the preexisting mutants on the short-term antiviral therapy during pregnancy.

Method

We enrolled 73 pregnant women with high HBV DNA load and telbivudine (TBV) treatment during pregnancy in this retrospective study. The UDPS was used to detect the HBV mutations before and after the TBV treatment.

Results

Before TBV treatment, the complexity of HBV quasispecies of all subjects was 0.40 ± 0.09; 41.1% (30/73) and 53.4% (39/73) subjects had rtM204I/V and rtN236 T/A detected, respectively; and 9.6% (7/73) patients had more than 20% frequency mutation of rtM204I/V, which was also similar with high frequency of rtN236 T/A mutation (41.1% vs. 53.4%, P=0.136; frequencies >20%: 9.6% vs. 5.5%, P=0.347). After TBV treatment, 71.2% (52/73) subjects had HBV DNA load ≥ 10³ IU/mL at delivery. Among them, 75.0% of patients with rtM204I positive had HBV DNA load ≥10³ IU/mL at delivery, which was comparable with the subjects without rtM204I (75.0% vs. 70.8%, P=0.710). No changes were found in the frequencies and the complexity of HBV quasispecies of rtM204I mutation after the TVB treatment.

Conclusion

The prevalence of preexisting drug-resistant mutations among pregnant women was high using UPDS. However, the preexisting HBV mutation had limited influence on the efficacy of short-term TBV treatment, and TBV treatment during late pregnancy seemed not to increase the risk of emerging HBV-resistant mutants.

Population Disequilibrium as Promoter of Adaptive Explorations in Hepatitis C Virus

Replication of RNA viruses is characterized by exploration of sequence space which facilitates their adaptation to changing environments. It is generally accepted that such exploration takes place mainly in response to positive selection, and that further diversification is boosted by modifications of virus population size, particularly bottleneck events. Our recent results with hepatitis C virus (HCV) have shown that the expansion in sequence space of a viral clone continues despite prolonged replication in a stable cell culture environment. Diagnosis of the expansion was based on the quantification of diversity indices, the occurrence of intra-population mutational waves (variations in mutant frequencies), and greater individual residue variations in mutant spectra than those anticipated from sequence alignments in data banks. In the present report, we review our previous results, and show additionally that mutational waves in amplicons from the NS5A-NS5B-coding region are equally prominent during HCV passage in the absence or presence of the mutagenic nucleotide analogues favipiravir or ribavirin. In addition, by extending our previous analysis to amplicons of the NS3- and NS5A-coding region, we provide further evidence of the incongruence between amino acid conservation scores in mutant spectra from infected patients and in the Los Alamos National Laboratory HCV data banks. We hypothesize that these observations have as a common origin a permanent state of HCV population disequilibrium even upon extensive viral replication in the absence of external selective constraints or changes in population size. Such a persistent disequilibrium—revealed by the changing composition of the mutant spectrum—may facilitate finding alternative mutational pathways for HCV antiviral resistance. The possible significance of our model for other genetically variable viruses is discussed.

Divergent Mutational Landscapes of Consensus and Minority Genotypes of West Nile Virus Demonstrate Host and Gene-Specific Evolutionary Pressures

Our current understanding of the natural evolution of RNA viruses comes largely from consensus level genetic analyses which ignore the diverse mutant swarms that comprise within-host viral populations. The breadth and composition of viral mutant swarms impact viral fitness and adaptation, and the capacity for swarm plasticity is likely to be particularly important for arthropod-borne viruses (arboviruses) that cycle between taxonomically divergent hosts. Despite this, characterization of the relationship between the selective pressures and genetic signatures of the mutant swarm and consensus sequences is lacking. To clarify this, we analyzed previously generated whole genome, deep-sequencing data from 548 West Nile virus samples isolated from avian tissues or mosquitoes in New York State from 1999-2018. Both consensus level (interhost) and minority level (intrahost) nucleotide and amino acid sequences were analyzed, and diversity at each position was calculated across the genome in order to assess the relationship between minority and consensus sequences for individual genes and hosts. Our results indicate that consensus sequences are an inept representation of the overall genetic diversity. Unique host and gene-specific signatures and selective pressures were identified. These data demonstrate that an accurate and comprehensive understanding of arbovirus evolution and adaptation within and between hosts requires consideration of minority genotypes.

Structures and Functions of Viral 5' Non-Coding Genomic RNA Domain-I in Group-B Enterovirus Infections

Group-B enteroviruses (EV-B) are ubiquitous naked single-stranded positive RNA viral pathogens that are responsible for common acute or persistent human infections. Their genome is composed in the 5′ end by a non-coding region, which is crucial for the initiation of the viral replication and translation processes. RNA domain-I secondary structures can interact with viral or cellular proteins to form viral ribonucleoprotein (RNP) complexes regulating viral genomic replication, whereas RNA domains-II to -VII (internal ribosome entry site, IRES) are known to interact with cellular ribosomal subunits to initiate the viral translation process. Natural 5′ terminally deleted viral forms lacking some genomic RNA domain-I secondary structures have been described in EV-B induced murine or human infections. Recent in vitro studies have evidenced that the loss of some viral RNP complexes in the RNA domain-I can modulate the viral replication and infectivity levels in EV-B infections. Moreover, the disruption of secondary structures of RNA domain-I could impair viral RNA sensing by RIG-I (Retinoic acid inducible gene I) or MDA5 (melanoma differentiation-associated protein 5) receptors, a way to overcome antiviral innate immune response. Overall, natural 5′ terminally deleted viral genomes resulting in the loss of various structures in the RNA domain-I could be major key players of host–cell interactions driving the development of acute or persistent EV-B infections.

Within-Host Viral Diversity: A Window into Viral Evolution

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.

An integrated software for virus community sequencing data analysis

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.

Engineering the Live-Attenuated Polio Vaccine to Prevent Reversion to Virulence

The live-attenuated oral poliovirus vaccine (OPV or Sabin vaccine) replicates in gut-associated tissues, eliciting mucosa and systemic immunity. OPV protects from disease and limits poliovirus spread. Accordingly, vaccination with OPV is the primary strategy used to end the circulation of all polioviruses. However, the ability of OPV to regain replication fitness and establish new epidemics represents a significant risk of polio re-emergence should immunization cease. Here, we report the development of a poliovirus type 2 vaccine strain (nOPV2) that is genetically more stable and less likely to regain virulence than the original Sabin2 strain. We introduced modifications within at the 5' untranslated region of the Sabin2 genome to stabilize attenuation determinants, 2C coding region to prevent recombination, and 3D polymerase to limit viral adaptability. Prior work established that nOPV2 is immunogenic in preclinical and clinical studies, and thus may enable complete poliovirus eradication.

Identification of highly conserved, serotype-specific dengue virus sequences: implications for vaccine design

Background

The sequence diversity of dengue virus (DENV) is one of the challenges in developing an effective vaccine against the virus. Highly conserved, serotype-specific (HCSS), immune-relevant DENV sequences are attractive candidates for vaccine design, and represent an alternative to the approach of selecting pan-DENV conserved sequences. The former aims to limit the number of possible cross-reactive epitope variants in the population, while the latter aims to limit the cross-reactivity between the serotypes to favour a serotype-specific response. Herein, we performed a large-scale systematic study to map and characterise HCSS sequences in the DENV proteome.

Methods

All reported DENV protein sequence data for each serotype was retrieved from the NCBI Entrez Protein (nr) Database (txid: 12637). The downloaded sequences were then separated according to the individual serotype proteins by use of BLASTp search, and subsequently removed for duplicates and co-aligned across the serotypes. Shannon’s entropy and mutual information (MI) analyses, by use of AVANA, were performed to measure the diversity within and between the serotype proteins to identify HCSS nonamers. The sequences were evaluated for the presence of promiscuous T-cell epitopes by use of NetCTLpan 1.1 and NetMHCIIpan 3.2 server for human leukocyte antigen (HLA) class I and class II supertypes, respectively. The predicted epitopes were matched to reported epitopes in the Immune Epitope Database.

Results

A total of 2321 nonamers met the HCSS selection criteria of entropy < 0.25 and MI > 0.8. Concatenating these resulted in a total of 337 HCSS sequences. DENV4 had the most number of HCSS nonamers; NS5, NS3 and E proteins had among the highest, with none in the C and only one in prM. The HCSS sequences were immune-relevant; 87 HCSS sequences were both reported T-cell epitopes/ligands in human and predicted epitopes, supporting the accuracy of the predictions. A number of the HCSS clustered as immunological hotspots and exhibited putative promiscuity beyond a single HLA supertype. The HCSS sequences represented, on average, ~ 40% of the proteome length for each serotype; more than double of pan-DENV sequences (conserved across the four serotypes), and thus offer a larger choice of sequences for vaccine target selection. HCSS sequences of a given serotype showed significant amino acid difference to all the variants of the other serotypes, supporting the notion of serotype-specificity.

Conclusion

This work provides a catalogue of HCSS sequences in the DENV proteome, as candidates for vaccine target selection. The methodology described herein provides a framework for similar application to other pathogens.

Kinetics of hepatitis B surface antigen quasispecies during REP 2139-Ca therapy in HBeAg-positive chronic HBV infection

Chronic HBV infection results in various clinical manifestations due to different levels of immune response. In recent years, hepatitis B treatment has improved by long-term administration of nucleos(t)ide analogues (NUCs) and peg-interferon. Nucleic acid polymers (NAPs; REP 2139-Ca and REP 2139-Mg) are new antiviral drugs that block the assembly of subviral particles, thus preventing the release of HBsAg and allowing its clearance and restoration of functional control of infection when combined with various immunotherapies. In the REP 102 study (NCT02646189), 9 of 12 patients showed substantial reduction of HBsAg and seroconversion to anti-HBs in response to REP 2139-Ca, whereas 3 of 12 patients did not show responses (>1 log reduction of HBsAg and HBV DNA from baseline). We characterized the dynamic changes of HBV quasispecies (QS) within the major hydrophilic region (MHR) of the 'pre-S/S' open reading frame including the 'a' determinant in responders and nonresponders of the REP 102 study and four untreated matched controls. HBV QS complexity at baseline varied slightly between responders and nonresponders (P = .28). However, these responders showed significant decline in viral complexity (P = .001) as REP 2139-Ca therapy progressed but no significant change in complexity was observed among the nonresponders (P = .99). The MHR mutations were more frequently observed in responders than in nonresponders and matched controls. No mutations were observed in 'a' determinant of major QS population which may interfere with the detection of HBsAg by diagnostic assays. No specific mutations were found within the MHR which could explain patients' poor HBsAg response during REP 2139-Ca therapy.

The ReFRAME library as a comprehensive drug repurposing library to identify mammarenavirus inhibitors

Several mammarenaviruses, chiefly Lassa virus (LASV) in Western Africa and Junín virus (JUNV) in the Argentine Pampas, cause severe disease in humans and pose important public health problems in their endemic regions. Moreover, mounting evidence indicates that the worldwide-distributed mammarenavirus lymphocytic choriomeningitis virus (LCMV) is a neglected human pathogen of clinical significance. The lack of licensed mammarenavirus vaccines and partial efficacy of current anti-mammarenavirus therapy limited to an off-label use of the nucleoside analog ribavirin underscore an unmet need for novel therapeutics to combat human pathogenic mammarenavirus infections. This task can be facilitated by the implementation of "drug repurposing" strategies to reduce the time and resources required to advance identified antiviral drug candidates into the clinic. We screened a drug repurposing library of 11,968 compounds (Repurposing, Focused Rescue and Accelerated Medchem [ReFRAME]) and identified several potent inhibitors of LCMV multiplication that had also strong anti-viral activity against LASV and JUNV. Our findings indicate that enzymes of the rate-limiting steps of pyrimidine and purine biosynthesis, the pro-viral MCL1 apoptosis regulator, BCL2 family member protein and the mitochondrial electron transport complex III, play critical roles in the completion of the mammarenavirus life cycle, suggesting they represent potential druggable targets to counter human pathogenic mammarenavirus infections.

Characterization of EIAV env Quasispecies during Long-Term Passage In Vitro: Gradual Loss of Pathogenicity

As the only widely used live lentiviral vaccine, the equine infectious anima virus (EIAV) attenuated vaccine was developed by in vitro passaging of a virulent strain for 121 generations. In our previous study, we observed that the attenuated vaccine was gradually selected under increased environmental pressure at the population level (termed a quasispecies). To further elucidate the potential correlation between viral quasispecies evolution and pathogenesis, a systematic study was performed by sequencing env using several methods. Some key mutations were identified within Env, and we observed that increased percentages of these mutations were accompanied by an increased passage number and attenuated virulence. Phylogenetic analysis revealed that env mutations related to the loss of virulence might have occurred evolutionarily. Among these mutations, deletion of amino acid 236 in the V4 region of Env resulted in the loss of one N-glycosylation site that was crucial for virulence. Notably, the 236-deleted sequence represented a "vaccine-specific" mutation that was also found in wild EIAV_LN40 strains based on single genome amplification (SGA) analysis. Therefore, our results suggest that the EIAV attenuated vaccine may originate from a branch of quasispecies of EIAV_LN40. Generally, the presented results may increase our understanding of the attenuation mechanism of the EIAV vaccine and provide more information about the evolution of other lentiviruses.

Subtype-Specific Prevalence of Hepatitis C Virus NS5A Resistance Associated Substitutions in Mainland China

Resistance associated substitutions (RASs) can reduce the efficacy of direct-acting antiviral agents (DAAs) targeting hepatitis C virus (HCV) and lead to treatment failure. Clinical data of HCV NS5A RASs prevalence are limited in China and need to be investigated. A total of 878 unique patient samples with different genotypes (GT) (1b: n = 489, 2a: n = 203, 3a: n = 60, 3b: n = 78, 6a: n = 48) were collected from around mainland China by KingMed Laboratory and analyzed for NS5A RASs distribution by Sanger sequencing. Phylogeographic analyses based on NS5A domain 1 sequences indicated circulation of both locally and nationally epidemic strains. Relatively high frequency of Y93H (14.1%) was only detected in GT1b but not in other subtypes. High frequency of L31M was found in both GT2a (95.6%) and GT3b (98.7%) sequences. Due to the overlapping incidence of A30K, 96% of GT3b isolates had NS5A RASs combination A30K + L31M, which confers high levels of resistance to most NS5A inhibitors. No RASs were detected in GT6a strains. Meanwhile, baseline NS5A RASs fingerprints were also evaluated in 185 DAA treatment-naive GT1b patients with next generation sequencing method. Patients presenting with Y93H had statistically higher entropy of HCV NS5A sequences. Taken together, subtype-specific distribution patterns of NS5A RASs were observed. GT1b patients with higher HCV complexity tend to have a greater chance of Y93H presence, while GT3b patients are naturally resistant to current NS5A inhibitors and their treatment may pose a challenge to real-world DAA application.

Analysis of HBV X gene quasispecies characteristics by next-generation sequencing and cloning-based sequencing and its association with hepatocellular carcinoma progression

OBJECTIVES

This study aimed to describe the differences between next-generation sequencing (NGS) and cloning-based sequencing (CBS) in HBX quasispecies research and primitively investigate the relationship between the dominant HBX quasispecies and hepatocellular carcinoma (HCC).

METHODS

A total of 12 serum samples were collected. Serum hepatitis B virus (HBV) DNA was extracted, and the HBV X-region (HBX) was amplified by nested polymerase chain reaction (PCR). The PCR products were simultaneously tested with NGS and CBS to detect quasispecies of the HBX.

RESULTS

A total of 9348 eligible quasispecies sequences were obtained by NGS, which were much larger than the 98 of that by CBS. By the phylogenetic tree, the dominant quasispecies sequence of each sample could be found, although they had several nucleotides differences between the dominant quasispecies sequences found by CBS and NGS. By comparing the quasispecies heterogeneity, it was found that the quasispecies complexity value of HBV X-region obtained by NGS was higher than CBS (P < 0.05). The diversity values, including d, dS, dN, an d d N/ dS obtained by NGS were lower than by CBS (all of P < 0.01). The relativity of Spearman(rs) in d, dS, and dN were statistically significant (rs_ d = 0.865, P = 0.001; rs_ dS = 0.722, P = 0.014; and rs_ dN = 0.738, P = 0.011, respectively). There were 21 different bases between the HBX quasispecies of case A and control B.

CONCLUSION

The results of this can be used as guidance when researchers plan to choose a suitable method to study quasispecies, especially the HBV X gene quasispecies. Some high-risk mutations of HBX quasispecies were also found in this study and their relationship with HCC need deeper exploration.

Viral quasispecies of hepatitis B virus in patients with YMDD mutation and lamivudine resistance may not predict the efficacy of lamivudine/adefovir rescue therapy

The association between hepatitis B virus (HBV) quasispecies (QS) and the efficacy of nucleos(t)ide analog therapy is currently not well defined, particularly in the case of lamivudine (LAM)/adefovir (ADV) combination rescue therapy for patients with chronic HBV infection (CHB) presenting with LAM resistance. In the present study, 16 CHB patients with the rtM204I/V mutation in the tyrosine-methionine-aspartate-aspartate motif of the C domain of the polymerase gene who switched to LAM/ADV treatment due to LAM resistance were assessed. HBV DNA was isolated from these patients and the reverse transcriptase (RT) region was sequenced. The QS heterogeneity and distribution was analyzed, the mutation sites were recorded and the phylogenetic trees were constructed. The results indicated that QS heterogeneity and distribution in the RT and S regions were not significantly different between responders (RS) and non-RS (NRS) at baseline (P>0.05), except for the higher frequency of a dominant strain in the RT region at the nucleotide level in the RS group (P=0.039). In addition, in NRS, no significant difference in QS heterogeneity or distribution in these regions was identified at six months vs. the baseline. Furthermore, although in the non-responder group the frequency of the LAM resistance-associated mutations (rtM204V/I) decreased at 6 months compared with the baseline, it did not disappear in any of the patients after six months of treatment. Analysis of individual patients did not indicate any consistent selection of specific HBV mutants during LAM/ADV rescue therapy. In conclusion, the baseline HBV QS within the RT and S regions may not be a valid predictor of the response to LAM/ADV rescue treatment in CHB patients with LAM resistance.

Stochastic Modeling and Simulation of Viral Evolution

RNA viruses comprise vast populations of closely related, but highly genetically diverse, entities known as quasispecies. Understanding the mechanisms by which this extreme diversity is generated and maintained is fundamental when approaching viral persistence and pathobiology in infected hosts. In this paper, we access quasispecies theory through a mathematical model based on the theory of multitype branching processes, to better understand the roles of mechanisms resulting in viral diversity, persistence and extinction. We accomplish this understanding by a combination of computational simulations and the theoretical analysis of the model. In order to perform the simulations, we have implemented the mathematical model into a computational platform capable of running simulations and presenting the results in a graphical format in real time. Among other things, we show that the establishment of virus populations may display four distinct regimes from its introduction into new hosts until achieving equilibrium or undergoing extinction. Also, we were able to simulate different fitness distributions representing distinct environments within a host which could either be favorable or hostile to the viral success. We addressed the most used mechanisms for explaining the extinction of RNA virus populations called lethal mutagenesis and mutational meltdown. We were able to demonstrate a correspondence between these two mechanisms implying the existence of a unifying principle leading to the extinction of RNA viruses.

Clinical implication and viral mutation in basal core promoter/pre-core of hepatitis B virus C/D recombinant

BACKGROUND AND AIM

Hepatitis B virus (HBV) C/D recombinant is predominant in Tibet in Western China. Although the geographical and ethnic distributions of the C/D recombinant have been described, the clinical implication and the characteristics of viral mutation in the basal core promoter (BCP)/pre-core (PC) region remain unclear.

METHODS

A total of 174 chronic HBV carriers, including 115 with chronic hepatitis B, 45 with liver cirrhosis, and 14 with hepatocellular carcinoma, were enrolled. Using next-generation sequencing, the S and BCP/PC genes were determined and analyzed.

RESULTS

Genotypes B, C2, D, and C/D recombinant were detected in 1.1% (2/174), 19.5% (34/174), 0.6% (1/174) and 78.7% (137/174) of the patients, respectively. The clinical parameters and viral mutation frequency in the BCP/PC region were compared between C2- and C/D recombinant-infected patients. The distribution of C2 and C/D did not differ by disease status or liver function. Significantly higher levels of HBV DNA (6.7 ± 1.6 vs. 5.9 ± 1.5, p = 0.014), HBeAg (263.5 vs. 20.0, p = 0.013) and A1762T/G1764A double-mutations (81.0 vs. 61.8%, p = 0.018), but a lower frequency of G1896A stop mutation (33.6 vs. 76.5%, p < 0.001) was observed in patients with the C/D recombinant than in patients with genotype C2. The clonal frequencies of A1762T, G1764A, G1896A and A1846T were lower in patients with C/D than C2.

CONCLUSION

The C/D recombinant has different mutation pattern in the BCP/PC region compared with genotype C2. The lower clonal frequencies of BCP/PC mutations may explain the higher levels of HBV DNA and HBeAg in C/D-infected patients.

Differential Shape of Geminivirus Mutant Spectra Across Cultivated and Wild Hosts With Invariant Viral Consensus Sequences

Geminiviruses (family Geminiviridae) possess single-stranded circular DNA genomes that are replicated by cellular polymerases in plant host cell nuclei. In their hosts, geminivirus populations behave as ensembles of mutant and recombinant genomes, known as viral quasispecies. This favors the emergence of new geminiviruses with altered host range, facilitating new or more severe diseases or overcoming resistance traits. In warm and temperate areas several whitefly-transmitted geminiviruses of the genus Begomovirus cause the tomato yellow leaf curl disease (TYLCD) with significant economic consequences. TYLCD is frequently controlled in commercial tomatoes by using the dominant Ty-1 resistance gene. Over a 45 day period we have studied the diversification of three begomoviruses causing TYLCD: tomato yellow leaf curl virus (TYLCV), tomato yellow leaf curl Sardinia virus (TYLCSV) and tomato yellow leaf curl Malaga virus (TYLCMaV, a natural recombinant between TYLCV and TYLCSV). Viral quasispecies resulting from inoculation of geminivirus infectious clones were examined in plants of susceptible tomato (ty-1/ty-1), heterozygous resistant tomato (Ty-1/ty-1), common bean, and the wild reservoir Solanum nigrum. Differences in virus fitness across hosts were observed while viral consensus sequences remained invariant. However, the complexity and heterogeneity of the quasispecies were high, especially in common bean and the wild host. Interestingly, the presence or absence of the Ty-1 allele in tomato did not lead to differences in begomovirus mutant spectra. However, the fitness decrease of TYLCSV and TYLCV in tomato at 45 dpi might be related to an increase in CP (Coat protein) mutation frequency. In Solanum nigrum the recombinant TYLCMaV, which showed lower fitness than TYLCSV, at 45 dpi actively explored Rep (Replication associated protein) ORF but not the overlapping C4. Our results underline the importance of begomovirus mutant spectra during infections. This is especially relevant in the wild reservoir of the viruses, which has the potential to maintain highly diverse mutant spectra without modifying their consensus sequences.

Mining a Kröhnke Pyridine Library for Anti-Arenavirus Activity

Several arenaviruses cause hemorrhagic fever (HF) disease in humans and represent important public health problems in their endemic regions. In addition, evidence indicates that the worldwide-distributed prototypic arenavirus lymphocytic choriomeningitis virus is a neglected human pathogen of clinical significance. There are no licensed arenavirus vaccines, and current antiarenavirus therapy is limited to an off-label use of ribavirin that is only partially effective. Therefore, there is an unmet need for novel therapeutics to combat human pathogenic arenaviruses, a task that will be facilitated by the identification of compounds with antiarenaviral activity that could serve as probes to identify arenavirus-host interactions suitable for targeting, as well as lead compounds to develop future antiarenaviral drugs. Screening of a combinatorial library of Krönhke pyridines identified compound KP-146 [(5-(5-(2,3-dihydrobenzo[ b][1,4] dioxin-6-yl)-4'-methoxy-[1,1'-biphenyl]-3-yl)thiophene-2-carboxamide] as having strong anti-lymphocytic choriomeningitis virus (LCMV) activity in cultured cells. KP-146 did not inhibit LCMV cell entry but rather interfered with the activity of the LCMV ribonucleoprotein (vRNP) responsible for directing virus RNA replication and gene transcription, as well as with the budding process mediated by the LCMV matrix Z protein. LCMV variants with increased resistance to KP-146 did not emerge after serial passages in the presence of KP-146. Our findings support the consideration of Kröhnke pyridine scaffold as a valuable source to identify compounds that could serve as tools to dissect arenavirus-host interactions, as well as lead candidate structures to develop antiarenaviral drugs.

Characterization of Equine Infectious Anemia Virus Long Terminal Repeat Quasispecies In Vitro and In Vivo

The equine infectious anemia virus (EIAV) attenuated vaccine was developed by long-term passaging of a field-isolated virulent strain in cross-species hosts, followed by successive cultivation in cells in vitro. To explore the molecular mechanism underlying the evolution of the EIAV attenuated vaccine, a systematic study focusing on long-terminal-repeat (LTR) variation in numerous virus strains ranging from virulent EIAV to attenuated EIAV was performed over time both in vitro and in vivo. Two hypervariable regions were identified within the U3 region in the enhancer region (EHR) and the negative regulatory element (NRE) and within the R region in the transcription start site (TSS) and the Tat-activating region (TAR). Among these sites, variation in the U3 region resulted in the formation of additional transcription factor binding sites; this variation of the in vitro-adapted strains was consistent with the loss of pathogenicity. Notably, the same LTR variation pattern was observed both in vitro and in vivo. Generally, the LTR variation in both the attenuated virus and the virulent strain fluctuated over time in vivo. Interestingly, the attenuated-virus-specific LTR variation was also detected in horses infected with the virulent strain, supporting the hypothesis that the evolution of an attenuated virus might have involved branching from EIAV quasispecies. This hypothesis was verified by phylogenetic analysis. The present systematic study examining the molecular evolution of attenuated EIAV from EIAV quasispecies may provide an informative model reflecting the evolution of similar lentiviruses.

IMPORTANCE The attenuated EIAV vaccine was the first lentiviral vaccine used to successfully control for equine infectious anemia in China. This vaccine provides an important reference for studying the relationship between EIAV gene variation and changes in biological characteristics. Importantly, the vaccine provides a model for the investigation of lentiviral quasispecies evolution. This study followed the “natural” development of the attenuated EIAV vaccine by use of a systematic analysis of LTR evolution in vitro and in vivo. The results revealed that the increase in LTR variation with passaging was accompanied by a decrease in virulence, which indicated that LTR variability might parallel the attenuation of virulence. Interestingly, the attenuated-virus-specific LTR variation was also detected in virulent-strain-infected horses, a finding consistent with those of previous investigations of gp90 and S2 evolution. Therefore, we present a hypothesis that the evolution of the attenuated virus may involve branching from EIAV quasispecies present in vivo.

Lethal mutagenesis of an RNA plant virus via lethal defection

Lethal mutagenesis is an antiviral therapy that relies on increasing the viral mutation rate with mutagenic nucleoside or base analogues. Currently, the molecular mechanisms that lead to virus extinction through enhanced mutagenesis are not fully understood. Increasing experimental evidence supports the lethal defection model of lethal mutagenesis of RNA viruses, where replication-competent-defectors drive infective virus towards extinction. Here, we address lethal mutagenesis in vivo using 5-fluorouracil (5-FU) during the establishment of tobacco mosaic virus (TMV) systemic infections in N. tabacum. The results show that 5-FU decreased the infectivity of TMV without affecting its viral load. Analysis of molecular clones spanning two genomic regions showed an increase of the FU-related base transitions A → G and U → C. Although the mutation frequency or the number of mutations per molecule did not increase, the complexity of the mutant spectra and the distribution of the mutations were altered. Overall, our results suggest that 5-FU antiviral effect on TMV is associated with the perturbation of the mutation-selection balance in the genomic region of the RNA-dependent RNA polymerase (RdRp). Our work supports the lethal defection model for lethal mutagenesis in vivo in a plant RNA virus and opens the way to study lethal mutagens in plant-virus systems.

The pH stability of foot-and-mouth disease virus

ᅟ

This review summarized the molecular determinants of the acid stability of FMDV in order to explore the uncoating mechanism of FMDV and improve the acid stability of vaccines.

Background

The foot-and-mouth disease virus (FMDV) capsid is highly acid labile and tends to dissociate into pentameric subunits at acidic condition to release viral RNA for initiating virus replication. However, the acid stability of virus capsid is greatly required for the maintenance of intact virion during the process of virus culture and vaccine production. The conflict between the acid lability in vivo and acid stability in vitro of FMDV capsid promotes the selection of a series of amino acid substitutions which can confer resistance to acid-induced FMDV inactivation. In order to explore the uncoating activity of FMDV and enhance the acid stability of vaccines, we summarized the available works about the pH stability of FMDV.

Main body of the abstract

In this review, we analyzed the intrinsic reasons for the acid instability of FMDV from the structural and functional aspects. We also listed all substitutions obtained by different research methods and showed them in the partial capsid of FMDV. We found that a quadrangle region in the viral capsid was the place where a great many pH-sensitive residues were distributed. As the uncoating event of FMDV is dependent on the pH-sensitive amino acid residues in the capsid, this most pH-sensitive position indicates a potential candidate location for RNA delivery triggered by the acid-induced coat disassociation.

Short conclusion

This review provided an overview of the pH stability of FMDV. The study of pH stability of FMDV not only contributes to the exploration of molecule and mechanism information for FMDV uncoating, but also enlightens the development of FMDV vaccines, including the traditionally inactivated vaccines and the new VLP (virus-like particle) vaccines.

A Selective Bottleneck Shapes the Evolutionary Mutant Spectra of Enterovirus A71 during Viral Dissemination in Humans

RNA viruses accumulate mutations to rapidly adapt to environmental changes. Enterovirus A71 (EV-A71) causes various clinical manifestations with occasional severe neurological complications. However, the mechanism by which EV-A71 evolves within the human body is unclear. Utilizing deep sequencing and haplotype analyses of viruses from various tissues of an autopsy patient, we sought to define the evolutionary pathway by which enterovirus A71 evolves fitness for invading the central nervous system in humans. Broad mutant spectra with divergent mutations were observed at the initial infection sites in the respiratory and digestive systems. After viral invasion, we identified a haplotype switch and dominant haplotype, with glycine at VP1 residue 31 (VP1-31G) in viral particles disseminated into the integumentary and central nervous systems. In vitro viral growth and fitness analyses indicated that VP1-31G conferred growth and a fitness advantage in human neuronal cells, whereas VP1-31D conferred enhanced replication in human colorectal cells. A higher proportion of VP1-31G was also found among fatal cases, suggesting that it may facilitate central nervous system infection in humans. Our data provide the first glimpse of EV-A71 quasispecies from oral tissues to the central nervous system within humans, showing broad implications for the surveillance and pathogenesis of this reemerging viral pathogen.IMPORTANCE EV-A71 continues to be a worldwide burden to public health. Although EV-A71 is the major etiological agent of hand, foot, and mouth disease, it can also cause neurological pulmonary edema, encephalitis, and even death, especially in children. Understanding selection processes enabling dissemination and accurately estimating EV-A71 diversity during invasion in humans are critical for applications in viral pathogenesis and vaccine studies. Here, we define a selection bottleneck appearing in respiratory and digestive tissues. Glycine substitution at VP1 residue 31 helps viruses break through the bottleneck and invade the central nervous system. This substitution is also advantageous for replication in neuronal cells in vitro Considering that fatal cases contain enhanced glycine substitution at VP1-31, we suggest that the increased prevalence of VP1-31G may alter viral tropism and aid central nervous system invasion. Our findings provide new insights into a dynamic mutant spectral switch active during acute viral infection with emerging viral pathogens.

A comparative study on the characterization of hepatitis B virus quasispecies by clone-based sequencing and third-generation sequencing

Hepatitis B virus (HBV) has a high mutation rate due to the extremely high replication rate and the proofreading deficiency during reverse transcription. The generated variants with genetic heterogeneity are described as viral quasispecies (QS). Clone-based sequencing (CBS) is thought to be the ‘gold standard’ for assessing QS complexity and diversity of HBV, but an important issue about CBS is cost-effectiveness and laborious. In this study, we investigated the utility of the third-generation sequencing (TGS) DNA sequencing to characterize genetic heterogeneity of HBV QS and assessed the possible contribution of TGS technology in HBV QS studies. Parallel experiments including 3 control samples, which consisted of HBV full gene genotype B and genotype C plasmids, and 10 patients samples were performed by using CBS and TGS to analyze HBV whole-genome QS. Characterization of QS heterogeneity was conducted by using comprehensive statistical analysis. The results showed that TGS had a high consistency with CBS when measuring the complexity and diversity of QS. In addition, to detect rare variants, there were strong advantages conferred by TGS. In summary, TGS was considered to be practicable in HBV QS studies and it might have a relevant role in the clinical management of HBV infection in the future.

Sequence-specific detection of different strains of LCMV in a single sample using tentacle probes

BACKGROUND

Virus infections often result in quasispecies of viral strains that can have dramatic impacts on disease outcomes. However, sequencing of viruses to determine strain composition is time consuming and often cost-prohibitive. Rapid, cost-effective methods are needed for accurate measurement of virus diversity to understand virus evolution and can be useful for experimental systems.

METHODS

We have developed a novel molecular method for sequence-specific detection of RNA virus genetic variants called Tentacle Probes. The probes are modified molecular beacons that have dramatically improved false positive rates and specificity in routine qPCR. To validate this approach, we have designed Tentacle Probes for two different strains of Lymphocytic Choriomeningitis Virus (LCMV) that differ by only 3 nucleotide substitutions, the parental Armstrong and the more virulent Clone-13 strain. One of these mutations is a missense mutation in the receptor protein GP1 that leads to the Armstrong strain to cause an acute infection and Clone-13 to cause a chronic infection instead. The probes were designed using thermodynamic calculations for hybridization between target or non-target sequences and the probe.

RESULTS

Using this approach, we were able to distinguish these two strains of LCMV individually by a single nucleotide mutation. The assay showed high reproducibility among different concentrations of viral cDNA, as well as high specificity and sensitivity, especially for the Clone-13 Tentacle Probe. Furthermore, in virus mixing experiments we were able to detect less than 10% of Clone-13 cDNA diluted in Armstrong cDNA.

CONCLUSIONS

Thus, we have developed a fast, cost-effective approach for identifying Clone-13 strain in a mix of other LCMV strains.

Molecular and pathobiological characterization of 61 Potato mop-top virus full-length cDNAs reveals great variability of the virus in the centre of potato domestication, novel genotypes and evidence for recombination

The evolutionary divergence of Potato mop-top virus (PMTV), a tri-partite, single-stranded RNA virus, is exceptionally low, based on the analysis of sequences obtained from isolates from Europe, Asia and North America. In general, RNA viruses exist as dynamic populations of closely related and recombinant genomes that are subjected to continuous genetic variation. The reason behind the low genetic variation of PMTV remains unclear. The question remains as to whether the low variability is a shared property of all PMTV isolates or is a result of the limited number of isolates characterized so far. We hypothesized that higher divergence of the virus might exist in the Andean regions of South America, the centre of potato domestication. Here, we report high variability of PMTV isolates collected from 12 fields in three locations in the Andean region of Peru. To evaluate PMTV genetic variation in Peru, we generated full-length cDNA clones, which allowed reliable comparative molecular and pathobiological characterization of individual isolates. We found significant divergence of the CP-RT and 8K sequences. The 8K cistron, which encodes a viral suppressor of RNA silencing, was found to be under diversifying selection. Phylogenetic analysis determined that, based on the CP-RT sequence, all PMTV isolates could be categorized into three separate lineages (clades). Moreover, we found evidence for recombination between two clades. Using infectious cDNA clones of the representatives of these two clades, as well as reassortants for the RNA-CP genomic component, we determined the pathobiological differences between the lineages, which we coined as S (for severe) and M (for mild) types. Interestingly, all isolates characterized previously (from Europe, Asia and North America) fall into the S-type clade, whereas most of the Peruvian isolates belong to the M-type. Taken together, our results support the notion of the single introduction of PMTV from the centre of potato origin to Europe, and subsequent spread of the S-type into Asia and USA. This is also supported by the suggested novel classification of isolates based on genetic constellations.

Strain-dependent mutational effects for Pepino mosaic virus in a natural host

Background

Pepino mosaic virus (PepMV) is an emerging plant pathogen that infects tomatoes worldwide. Understanding the factors that influence its evolutionary success is essential for developing new control strategies that may be more robust against the evolution of new viral strains. One of these evolutionary factors is the distribution of mutational fitness effect (DMFE), that is, the fraction of mutations that are lethal, deleterious, neutral, and beneficial on a given viral strain and host species. The goal of this study was to characterize the DMFE of introduced nonsynonymous mutations on a mild isolate of PepMV from the Chilean 2 strain (PepMV-P22). Additionally, we also explored whether the fitness effect of a given mutation depends on the gene where it appears or on epistatic interactions with the genetic background. To address this latter possibility, a subset of mutations were also introduced in a mild isolate of the European strain (PepMV-P11) and the fitness of the resulting clones measured.

Results

A collection of 25 PepMV clones each containing a single nucleotide nonsynonymous substitution was created by site-directed mutagenesis and the fitness of each mutant was determined. PepMV-P22 genome showed a high degree of robustness against point mutations, with 80% of mutations being either neutral or even beneficial and only 20% being deleterious or lethal. We found that the effect of mutations strongly depended on the gene in which they were introduced. Mutations with the largest average beneficial effects were those affecting the RdRp gene, in contrast to mutations affecting TGB1 and CP genes, for which the average effects were deleterious. Moreover, significant epistatic interactions were observed between nonsynonymous mutations and the genetic background, meaning that the effect of a given nucleotide substitution on a particular genomic context cannot be predicted by knowing its effect in a different one.

Conclusions

Our results indicated that PepMV genome has a surprisingly high robustness against mutations. We also found that fitness consequences of a given mutation differ between the two strains analyzed. This discovery suggests that the strength of selection, and thus the rates of evolution, vary among PepMV strains.

Lack of impact of pre-existing T97A HIV-1 integrase mutation on integrase strand transfer inhibitor resistance and treatment outcome

T97A is an HIV-1 integrase polymorphism associated with integrase strand transfer inhibitor (INSTI) resistance. Using pooled data from 16 clinical studies, we investigated the prevalence of T97A (pre-existing and emergent) and its impact on INSTI susceptibility and treatment response in INSTI-naive patients who enrolled on elvitegravir (EVG)- or raltegravir (RAL)-based regimens. Prior to INSTI-based therapy, primary INSTI resistance-associated mutations (RAMs) were absent and T97A pre-existed infrequently (1.4%; 47 of 3367 integrase sequences); most often among non-B (5.3%) than B (0.9%) HIV-1 subtypes. During INSTI-based therapy, few patients experienced virologic failure with emergent INSTI RAMs (3%; 122 of 3881 patients), among whom T97A emerged infrequently in the presence (n = 6) or absence (n = 8) of primary INSTI RAMs. A comparison between pre-existing and emergent T97A patient populations (i.e., in the absence of primary INSTI RAMs) showed no significant differences in EVG or RAL susceptibility in vitro. Furthermore, among all T97A-containing viruses tested, only 38-44% exhibited reduced susceptibility to EVG and/or RAL (all of low magnitude; <11-fold), while all maintained susceptibility to dolutegravir. Of the patients with pre-existing T97A, 17 had available clinical follow-up: 16 achieved virologic suppression and 1 maintained T97A and INSTI sensitivity without further resistance development. Overall, T97A is an infrequent integrase polymorphism that is enriched among non-B HIV-1 subtypes and can confer low-level reduced susceptibility to EVG and/or RAL. However, detection of T97A does not affect response to INSTI-based therapy with EVG or RAL. These results suggest a very low risk of initiating INSTI-based therapy in patients with pre-existing T97A.

Characteristics of CpG Islands and their quasispecies of full-length hepatitis B virus genomes from patients at different phases of infection

BACKGROUND

CpG islands in hepatitis B virus (HBV) genome are potential targets for methylation mediated gene silencing, and may be involved in the pathogenesis of HBV infection. To date, their characteristics in HBV quasispecies (QS) remain largely unknown. The purpose of this study was to investigate the characteristics of CpG islands in HBV QS.

METHODS

Forty patients diagnosed as acute hepatitis B (AHB, n = 10), immune-tolerant HBV carriers (IT, n = 9), chronic hepatitis B (CHB, n = 11), or acute on chronic liver failure (ACLF, n = 10), were enrolled in this case-control study. A total of 599 clones were isolated, and full-length HBV genomes were sequenced.

RESULTS

CpG island II (CGII) in AHB group was shorter in length and its QS heterogeneity was lower than that in the chronic infection group. Among the chronic infection subgroups, CGII and CpG island III (CGIII) in IT group were longer and their heterogeneity was lower compared to CHB and ACLF groups. Length of CGII correlated with HBV DNA levels positively while the complexity and diversity of CGII correlated with HBV DNA levels negatively. Moreover, CGII and CGIII were shorter in genotype B than those in genotype C, while QS complexity and diversity of either CGII or CGIII had no significant difference between genotype B and C.

CONCLUSIONS

Overall, our results suggest that the distribution, length and QS heterogeneity of CpG islands in full-length HBV genome differ across clinical phases of infection, of which the mechanism warrants further study.

DNA from Dust: Comparative Genomics of Large DNA Viruses in Field Surveillance Samples

The intensification of the poultry industry over the last 60 years facilitated the evolution of increased virulence and vaccine breaks in Marek's disease virus (MDV-1). Full-genome sequences are essential for understanding why and how this evolution occurred, but what is known about genome-wide variation in MDV comes from laboratory culture. To rectify this, we developed methods for obtaining high-quality genome sequences directly from field samples without the need for sequence-based enrichment strategies prior to sequencing. We applied this to the first characterization of MDV-1 genomes from the field, without prior culture. These viruses were collected from vaccinated hosts that acquired naturally circulating field strains of MDV-1, in the absence of a disease outbreak. This reflects the current issue afflicting the poultry industry, where virulent field strains continue to circulate despite vaccination and can remain undetected due to the lack of overt disease symptoms. We found that viral genomes from adjacent field sites had high levels of overall DNA identity, and despite strong evidence of purifying selection, had coding variations in proteins associated with virulence and manipulation of host immunity. Our methods empower ecological field surveillance, make it possible to determine the basis of viral virulence and vaccine breaks, and can be used to obtain full genomes from clinical samples of other large DNA viruses, known and unknown. IMPORTANCE Despite both clinical and laboratory data that show increased virulence in field isolates of MDV-1 over the last half century, we do not yet understand the genetic basis of its pathogenicity. Our knowledge of genome-wide variation between strains of this virus comes exclusively from isolates that have been cultured in the laboratory. MDV-1 isolates tend to lose virulence during repeated cycles of replication in the laboratory, raising concerns about the ability of cultured isolates to accurately reflect virus in the field. The ability to directly sequence and compare field isolates of this virus is critical to understanding the genetic basis of rising virulence in the wild. Our approaches remove the prior requirement for cell culture and allow direct measurement of viral genomic variation within and between hosts, over time, and during adaptation to changing conditions.

An eco-epidemiological study of Morbilli-related paramyxovirus infection in Madagascar bats reveals host-switching as the dominant macro-evolutionary mechanism

An eco-epidemiological investigation was carried out on Madagascar bat communities to better understand the evolutionary mechanisms and environmental factors that affect virus transmission among bat species in closely related members of the genus Morbillivirus, currently referred to as Unclassified Morbilli-related paramyxoviruses (UMRVs). A total of 947 bats were investigated originating from 52 capture sites (22 caves, 18 buildings, and 12 outdoor sites) distributed over different bioclimatic zones of the island. Using RT-PCR targeting the L-polymerase gene of the Paramyxoviridae family, we found that 10.5% of sampled bats were infected, representing six out of seven families and 15 out of 31 species analyzed. Univariate analysis indicates that both abiotic and biotic factors may promote viral infection. Using generalized linear modeling of UMRV infection overlaid on biotic and abiotic variables, we demonstrate that sympatric occurrence of bats is a major factor for virus transmission. Phylogenetic analyses revealed that all paramyxoviruses infecting Malagasy bats are UMRVs and showed little host specificity. Analyses using the maximum parsimony reconciliation tool CoRe-PA, indicate that host-switching, rather than co-speciation, is the dominant macro-evolutionary mechanism of UMRVs among Malagasy bats.

Dynamics of Genotypic Mutations of the Hepatitis B Virus Associated With Long-Term Entecavir Treatment Determined With Ultradeep Pyrosequencing: A Retrospective Observational Study

The aim of the study is to explore the evolution of genotypic mutations within the reverse transcriptase region in partial virological responders (PVRs) receiving long-term entecavir (ETV) treatment. A total of 32 patients were classified as completely virological responders (CVRs) (n = 12) or PVRs (n = 20). Five partial responders were hepatitis B virus (HBV)-DNA positive after long-term therapy, which lasted for >3 years. A total of 71 serum samples from these 32 patients were assayed by ultra-deep pyrosequencing (UDPS): 32 samples were from all patients at baseline, and 39 were from PVRs with sequential inter-treatment. Approximately 84,708 sequences were generated per sample. At baseline, the quasispecies heterogeneity did not significantly differ between the 2 groups. The frequencies of substitutions indicating pre-existence of nucleos(t)ide analog resistant (NAr) mutants ranged from 0.10% to 6.70%, which did not statistically differ between groups either. However, the substitutions associated with the NAr mutants were significantly different from those associated with the non-NAr mutants in 13 patients; 6 of these patients were PVRs and the others were CVRs. Five patients were HBV DNA positive after regular ETV monotherapy for >3 years, and 4 of these patients underwent mild NAr substitution fluctuations (<20%). One patient developed virological breakthrough while bearing single, double, and triple (rtL180 M, rtM204 V, rtS202G) substitutions. In addition to the common substitutions, unknown amino acid substitutions, such as rtL145 M/S, rtF151Y/L, rtR153Q, rtI224 V, rtN248H, rtS223A, rtS256C, need to be further verified. NAr substitutions are observed at frequencies of 0.10% to 6.7% before therapy. Long-term ETV therapy generally results in virological responses, as long as the proportion of resistance mutations remains at a relatively low level. Genotypic resistance to ETV is detected in all PVRs receiving long-term ETV therapy.

Mixed group of Rhizobiales microbes in lung and blood of a patient with fatal pulmonary illness

We examined the microbial composition in the diseased lung and early-phase microbial cultures from the blood of a patient with a rapidly progressing fatal pulmonary illness. Although no microbes could be isolated from such cultures during the initial study, the HTS-microbiome study revealed the presence of a unique mixture of alphaproteobacteria, composed mainly of different families of Rhizobiales microbes. Microbial 16S rDNA sequences matching closely to Afipia cberi were identified mainly in the patient's diseased lung tissue, but only rarely in the early-phase blood cultures. Conversely, the high abundance of sequences found in early-phase blood cultures of different broth media matched closely with those of the families Methylobacteriaceae, Phyllobacteriaceae and Sphingomonadaceae. The two species that successfully adapted to grow in a laboratory culture system were A. cberi and Mesorhizobium hominis, which eventually were isolated from a previously cryopreserved blood culture of SP4 broth. Many other species, including members of the Bradyrhizobiaceae and Phyllobacteriaceae families, and all members of the Methylobacteriaceae and Sphingomonadaceae families identified by HTS remained non-cultivated. We developed specific PCR primers and FISH probes, which detected the target Rhizobiales microbes in former blood cultures and autopsy lung tissues. It is unclear what role these Rhizobiales microbes might have played in the patient's complex disease process. However, the above mentioned assays should help in rapidly detecting and identifying these previously unrecognized Rhizobiales microbes in patients.

Deep Sequencing for Evaluation of Genetic Stability of Influenza A/California/07/2009 (H1N1) Vaccine Viruses

Virus growth during influenza vaccine manufacture can lead to mutations that alter antigenic properties of the virus, and thus may affect protective potency of the vaccine. Different reassortants of pandemic "swine" H1N1 influenza A vaccine (121XP, X-179A and X-181) viruses as well as wild type A/California/07/2009(H1N1) and A/PR/8/34 strains were propagated in embryonated eggs and used for DNA/RNA Illumina HiSeq and MiSeq sequencing. The RNA sequences of these viruses published in NCBI were used as references for alignment of the sequencing reads generated in this study. Consensus sequences of these viruses differed from the NCBI-deposited sequences at several nucleotides. 121XP stock derived by reverse genetics was more heterogeneous than X-179A and X-181 stocks prepared by conventional reassortant technology. Passaged 121XP virus contained four non-synonymous mutations in the HA gene. One of these mutations (Lys₂₂₆Glu) was located in the Ca antigenic site of HA (present in 18% of the population). Two non-synonymous mutations were present in HA of viruses derived from X-179A: Pro₃₁₄Gln (18%) and Asn₁₄₆Asp (78%). The latter mutation located in the Sa antigenic site was also detected at a low level (11%) in the wild-type A/California/07/2009(H1N1) virus, and was present as a complete substitution in X-181 viruses derived from X-179A virus. In the passaged X-181 viruses, two mutations emerged in HA: a silent mutation A₁₃₉₈G (31%) in one batch and G₇₅₆T (Glu₂₅₂Asp, 47%) in another batch. The latter mutation was located in the conservative region of the antigenic site Ca. The protocol for RNA sequencing was found to be robust, reproducible, and suitable for monitoring genetic consistency of influenza vaccine seed stocks.

Antigenic cooperation among intrahost HCV variants organized into a complex network of cross-immunoreactivity

Hepatitis C virus (HCV) has the propensity to cause chronic infection. Continuous immune escape has been proposed as a mechanism of intrahost viral evolution contributing to HCV persistence. Although the pronounced genetic diversity of intrahost HCV populations supports this hypothesis, recent observations of long-term persistence of individual HCV variants, negative selection increase, and complex dynamics of viral subpopulations during infection as well as broad cross-immunoreactivity (CR) among variants are inconsistent with the immune-escape hypothesis. Here, we present a mathematical model of intrahost viral population dynamics under the condition of a complex CR network (CRN) of viral variants and examine the contribution of CR to establishing persistent HCV infection. The model suggests a mechanism of viral adaptation by antigenic cooperation (AC), with immune responses against one variant protecting other variants. AC reduces the capacity of the host's immune system to neutralize certain viral variants. CRN structure determines specific roles for each viral variant in host adaptation, with variants eliciting broad-CR antibodies facilitating persistence of other variants immunoreacting with these antibodies. The proposed mechanism is supported by empirical observations of intrahost HCV evolution. Interference with AC is a potential strategy for interruption and prevention of chronic HCV infection.

Genetic Variability of Hepatitis C Virus (HCV) 5' Untranslated Region in HIV/HCV Coinfected Patients Treated with Pegylated Interferon and Ribavirin

Association between hepatitis C virus (HCV) quasispecies and treatment outcome among patients with chronic hepatitis C has been the subject of many studies. However, these studies focused mainly on viral variable regions (E1 and E2) and usually did not include human immunodeficiency virus (HIV)-positive patients. The aim of the present study was to analyze heterogeneity of the 5' untranslated region (5'UTR) in HCV/HIV coinfected patients treated with interferon and ribavirin. The HCV 5'UTR was amplified from serum and peripheral blood mononuclear cells (PBMC) samples in 37 HCV/HIV coinfected patients treated for chronic hepatitis C. Samples were collected right before treatment, and at 2, 4, 6, 8, 12, 20, 24, 36, 44, 48, 60, and 72 weeks. Heterogeneity of the 5'UTR was analyzed by single strand conformational polymorphism (SSCP), cloning and sequencing. Sustained virological response (SVR) was achieved in 46% of analyzed HCV/HIV co-infected patients. Stable SSCP band pattern was observed in 22 patients (62.9%) and SVR rate among these patients was 23%. Decline in the number of bands and/or shift in band positions were found in 6 patients (17.1%), 5 (83%) of whom achieved SVR (p=0.009). A novel viral genotype was identified in all but one of these patients. In 5 of these 6 patients a new genotype was dominant. 5'UTR heterogeneity may correlate with interferon and ribavirin treatment outcome. In the analyzed group of HCV/HIV coinfected patients, viral quasispecies stability during treatment favored viral persistence, whereas decrease in the number of variants and/or emergence of new variants was associated with SVR. Among injection drug users (IDU) patients, a new genotype may become dominant during treatment, probably due to the presence of mixed infections with various strains, which have different susceptibility to treatment.

Deep sequencing and phylogenetic analysis of variants resistant to interferon-based protease inhibitor therapy in chronic hepatitis induced by genotype 1b hepatitis C virus

Because of recent advances in deep sequencing technology, detailed analysis of hepatitis C virus (HCV) quasispecies and their dynamic changes in response to direct antiviral agents (DAAs) became possible, although the role of quasispecies is not fully understood. In this study, to clarify the evolution of viral quasispecies and the origin of drug-resistant mutations induced by interferon (IFN)-based protease inhibitor therapy, the nonstructural-3 (NS3) region of genotype 1b HCV in 34 chronic hepatitis patients treated with telaprevir (TVR)/pegylated interferon (PEG-IFN)/ribavirin (RBV) was subjected to a deep sequencing study coupled with phylogenetic analysis. Twenty-six patients (76.5%) achieved a sustained viral response (SVR), while 8 patients did not (non-SVR; 23.5%). When the complexity of the quasispecies was expressed as the mutation frequency or Shannon entropy value, a significant decrease in the IFNL3 (rs8099917) TT group and a marginal decrease in the SVR group were found soon (12 h) after the introduction of treatment, whereas there was no decrease in the non-SVR group and no significant decrease in mutation frequency in the IFNL3 TG/GG group. In the analysis of viral quasispecies composition in non-SVR patients, major populations greatly changed, accompanied by the appearance of resistance, and the compositions were unlikely to return to the pretreatment composition even after the end of therapy. Clinically TVR-resistant variants were observed in 5 non-SVR patients (5/8, 62.5%), all of which were suspected to have acquired resistance by mutations through phylogenetic analysis. In conclusion, results of the study have important implications for treatment response and outcome in interferon-based protease inhibitor therapy.

IMPORTANCE

In the host, hepatitis C virus (HCV) consists of a variety of populations (quasispecies), and it is supposed that dynamic changes in quasispecies are closely related to pathogenesis, although this is poorly understood. In this study, recently developed deep sequencing technology was introduced, and changes in quasispecies associated with telaprevir (TVR)/pegylated interferon (PEG-IFN)/ribavirin (RBV) triple therapy and their clinical significance were investigated extensively by phylogenetic tree analysis. Through this study, the associations among treatment response, changes in viral quasispecies complexity in the early stage of treatment, changes in the quasispecies composition, and origin of TVR-resistant variant HCV were elucidated.

Similar prevalence of low-abundance drug-resistant variants in treatment-naive patients with genotype 1a and 1b hepatitis C virus infections as determined by ultradeep pyrosequencing

Background and Objectives

Hepatitis C virus (HCV) variants that confer resistance to direct-acting-antiviral agents (DAA) have been detected by standard sequencing technology in genotype (G) 1 viruses from DAA-naive patients. It has recently been shown that virological response rates are higher and breakthrough rates are lower in G1b infected patients than in G1a infected patients treated with certain classes of HCV DAAs. It is not known whether this corresponds to a difference in the composition of G1a and G1b HCV quasispecies in regards to the proportion of naturally occurring DAA-resistant variants before treatment.

Methods

We used ultradeep pyrosequencing to determine the prevalence of low-abundance (<25% of the sequence reads) DAA-resistant variants in 191 NS3 and 116 NS5B isolates from 208 DAA-naive G1-infected patients.

Results

A total of 3.5 million high-quality reads of ≥200 nucleotides were generated. The median coverage depth was 4150x and 4470x per NS3 and NS5B amplicon, respectively. Both G1a and G1b populations showed Shannon entropy distributions, with no difference between G1a and G1b in NS3 or NS5B region at the nucleotide level. A higher number of substitutions that confer resistance to protease inhibitors were observed in G1a isolates (mainly at amino acid 80 of the NS3 region). The prevalence of amino acid substitutions that confer resistance to NS5B non-nucleoside inhibitors was similar in G1a and G1b isolates. The NS5B S282T variant, which confers resistance to the polymerase inhibitors mericitabine and sofosbuvir, was not detected in any sample.

Conclusion

The quasispecies genetic diversity and prevalence of DAA-resistant variants was similar in G1a and G1b isolates and in both NS3 and NS5B regions, suggesting that this is not a determinant for the higher level of DAA resistance observed across G1a HCV infected patients upon treatment.

Tautomerism provides a molecular explanation for the mutagenic properties of the anti-HIV nucleoside 5-aza-5,6-dihydro-2'-deoxycytidine

Viral lethal mutagenesis is a strategy whereby the innate immune system or mutagenic pool nucleotides increase the error rate of viral replication above the error catastrophe limit. Lethal mutagenesis has been proposed as a mechanism for several antiviral compounds, including the drug candidate 5-aza-5,6-dihydro-2'-deoxycytidine (KP1212), which causes A-to-G and G-to-A mutations in the HIV genome, both in tissue culture and in HIV positive patients undergoing KP1212 monotherapy. This work explored the molecular mechanism(s) underlying the mutagenicity of KP1212, and specifically whether tautomerism, a previously proposed hypothesis, could explain the biological consequences of this nucleoside analog. Establishing tautomerism of nucleic acid bases under physiological conditions has been challenging because of the lack of sensitive methods. This study investigated tautomerism using an array of spectroscopic, theoretical, and chemical biology approaches. Variable temperature NMR and 2D infrared spectroscopic methods demonstrated that KP1212 existed as a broad ensemble of interconverting tautomers, among which enolic forms dominated. The mutagenic properties of KP1212 were determined empirically by in vitro and in vivo replication of a single-stranded vector containing a single KP1212. It was found that KP1212 paired with both A (10%) and G (90%), which is in accord with clinical observations. Moreover, this mutation frequency is sufficient for pushing a viral population over its error catastrophe limit, as observed before in cell culture studies. Finally, a model is proposed that correlates the mutagenicity of KP1212 with its tautomeric distribution in solution.

RNA virus evolution at variable error rate

ABSTRACT: 

RNA viruses replicate their genomes with very high error rates, which leads to the generation of a large genetic diversity that makes them highly adaptable to most environmental pressures, including antiviral drugs and immune responses. However, since most mutations are deleterious, an excess of errors can be very negative for RNA viruses, entailing that error rates must be finely regulated. Currently, the manipulation of the error rate is emerging as a promising antiviral therapy that could minimize the problem of virus adaptation to classical treatments. This review provides a detailed analysis of the different outcomes that can result from the variation of the error rate in RNA viruses, on the basis of the more relevant findings obtained in experimental studies.

Standards for sequencing viral genomes in the era of high-throughput sequencing

Thanks to high-throughput sequencing technologies, genome sequencing has become a common component in nearly all aspects of viral research; thus, we are experiencing an explosion in both the number of available genome sequences and the number of institutions producing such data. However, there are currently no common standards used to convey the quality, and therefore utility, of these various genome sequences. Here, we propose five "standard" categories that encompass all stages of viral genome finishing, and we define them using simple criteria that are agnostic to the technology used for sequencing. We also provide genome finishing recommendations for various downstream applications, keeping in mind the cost-benefit trade-offs associated with different levels of finishing. Our goal is to define a common vocabulary that will allow comparison of genome quality across different research groups, sequencing platforms, and assembly techniques.