Search

We study a non-equilibrium dynamical model: a marked continuous contact model in $d$-dimensional space, $d \ge 1$. In contrast with the continuous contact model in a critical regime, see \cite{KKP}, \cite{KPZ}, the model under consideration is in the subcritical regime and it contains an additional spontaneous spatially homogeneous birth from an external source. We prove that this system has an invariant measure. We prove also that the process starting from any initial distribution converges to this invariant measure., arXiv admin note: substantial text overlap with arXiv:1601.07841

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., 34 pages, 10 figures, 4 tables, 1 appendix

Introduction: Genotype VI of avian avulavirus 1 (AAvV-1) has pigeons and doves as its reservoir and is often termed pigeon paramyxovirus type-1 (PPMV-1). The pathogenesis of PPMV-1 infections in poultry is largely obscure. It is known that PPMV-1 requires a series of passages in chickens before it becomes adapted to gallinaceous poultry. Material and Methods: Changes in the genome of PPMV-1 were analysed after serial passages in specific pathogen free (SPF) chickens, using high-throughput sequencing. Additionally, histopathological lesions induced by PPMV-1 in experimentally inoculated pigeons, chickens, and turkeys were evaluated. Results: Following six passages of PPMV-1 in chickens, 10 nonsynonymous substitutions were found including one (in the NP protein) which dominated the genetic pool of viral quasispecies. Histopathological changes induced by the post-passage PPMV-1 strain were more prominent than changes wrought by the pre-passaged PPMV-1 strain and the lesions were most intense in pigeons followed by chickens and turkeys. Conclusion: PPMV-1 is highly adapted to pigeons and passaging through chickens results in the acquisition of novel amino acids in the polymerase complex, which may alter the pathogenic potential of the virus.

Geoghegan and Holmes describe the history of evolutionary ideas in the study of viruses, showing that two different approaches to studying virus evolution—the comparative and the experimental—were both established in seminal papers published in the late... RNA viruses are diverse, abundant, and rapidly evolving. Genetic data have been generated from virus populations since the late 1970s and used to understand their evolution, emergence, and spread, culminating in the generation and analysis of many thousands of viral genome sequences. Despite this wealth of data, evolutionary genetics has played a surprisingly small role in our understanding of virus evolution. Instead, studies of RNA virus evolution have been dominated by two very different perspectives, the experimental and the comparative, that have largely been conducted independently and sometimes antagonistically. Here, we review the insights that these two approaches have provided over the last 40 years. We show that experimental approaches using in vitro and in vivo laboratory models are largely focused on short-term intrahost evolutionary mechanisms, and may not always be relevant to natural systems. In contrast, the comparative approach relies on the phylogenetic analysis of natural virus populations, usually considering data collected over multiple cycles of virus–host transmission, but is divorced from the causative evolutionary processes. To truly understand RNA virus evolution it is necessary to meld experimental and comparative approaches within a single evolutionary genetic framework, and to link viral evolution at the intrahost scale with that which occurs over both epidemiological and geological timescales. We suggest that the impetus for this new synthesis may come from methodological advances in next-generation sequencing and metagenomics.

Accurate and sensitive approaches to detect HIV-1 drug resistance mutations (DRMs) are indispensable for the paradigm of ‘treatment as prevention’. While HIV-1 proviral DNA allows sensitive high throughput sequencing (HTS)-based DRM detection, its applicability is limited by presence of defective genomes. This study demonstrates application of quasispecies reconstruction algorithms (QRAs) to improve DRM detection sensitivity from proviral DNA. A robust benchmarking of 5 QRAs was performed with 2 distinct experimental control-datasets including a stringent, novel control: DCPM, simulating in-vivo variant distribution (0.08%-86.5%). Selected QRA was further evaluated for its ability to differentiate DRMs from hypermutated sequences using an in-silico control. PredictHaplo outperformed all others in terms of precision and was selected for further analysis. Near full-genome HTS was performed on proviral DNA from 20 HIV-1C infected individuals, at different stages of ART, from Mumbai, India. DRM detection was performed through residue-wise variation analysis and implementation of QRAs. Both analyses were highly concordant for DRM frequencies >10% (spearman r=0.91, p

RNA viruses demonstrate a vast range of variants, called quasispecies, due to error-prone replication by viral RNA-dependent RNA polymerase. Although live attenuated vaccines are effective in preventing RNA virus infection, there is a risk of reversal to virulence after their administration. To test the hypothesis that high-fidelity viral polymerase reduces the diversity of influenza virus quasispecies, resulting in inhibition of reversal of the attenuated phenotype, we first screened for a high-fidelity viral polymerase using serial virus passages under selection with a guanosine analog ribavirin. Consequently, we identified a Leu66-to-Val single amino acid mutation in polymerase basic protein 1 (PB1). The high-fidelity phenotype of PB1-L66V was confirmed using next-generation sequencing analysis and biochemical assays with the purified influenza viral polymerase. As expected, PB1-L66V showed at least two-times-lower mutation rates and decreased misincorporation rates, compared to the wild type (WT). Therefore, we next generated an attenuated PB1-L66V virus with a temperature-sensitive (ts) phenotype based on FluMist, a live attenuated influenza vaccine (LAIV) that can restrict virus propagation by ts mutations, and examined the genetic stability of the attenuated PB1-L66V virus using serial virus passages. The PB1-L66V mutation prevented reversion of the ts phenotype to the WT phenotype, suggesting that the high-fidelity viral polymerase could contribute to generating an LAIV with high genetic stability, which would not revert to the pathogenic virus. IMPORTANCE The LAIV currently in use is prescribed for actively immunizing individuals aged 2 to 49 years. However, it is not approved for infants and elderly individuals, who actually need it the most, because it might prolong virus propagation and cause an apparent infection in these individuals, due to their weak immune systems. Recently, reversion of the ts phenotype of the LAIV strain currently in use to a pathogenic virus was demonstrated in cultured cells. Thus, the generation of mutations associated with enhanced virulence in LAIV should be considered. In this study, we isolated a novel influenza virus strain with a Leu66-to-Val single amino acid mutation in PB1 that displayed a significantly higher fidelity than the WT. We generated a novel LAIV candidate strain harboring this mutation. This strain showed higher genetic stability and no ts phenotype reversion. Thus, our high-fidelity strain might be useful for the development of a safer LAIV.

Background Telbivudine (LdT) - a synthetic thymidine β -L-nucleoside analogue (NA) - is an effective inhibitor for hepatitis B virus (HBV) replication. The quasispecies spectra in the reverse transcriptase (RT) region of the HBV genome and their dynamic changes associated with LdT treatment remains largely unknown. Methods We prospectively recruited a total of 21 treatment-naive patients with chronic HBV infection and collected sequential serum samples at five time points (baseline, weeks 1, 3, 12, and 24 after LdT treatment). The HBV RT region was amplified and shotgun-sequenced by the Ion Torrent Personal Genome Machine (PGM)® system. We reconstructed full-length haplotypes of the RT region using an integrated bioinformatics framework, including de novo contig assembly and full-length haplotype reconstruction. In addition, we investigated the quasispecies' dynamic changes and evolution history and characterized potential NAs resistant mutations over the treatment course. Results Viral quasispecies differed obviously between patients with complete (n = 8) and incomplete/no response (n = 13) at 12 weeks after LdT treatment. A reduced d N / d S ratio in quasispecies demonstrated a selective constraint resulting from antiviral therapy. The temporal clustering of sequential quasispecies showed different patterns along with a 24-week observation, although its statistic did not differ significantly. Several patients harboring pre-existing resistant mutations showed different clinical responses, while NAs resistant mutations were rare within a short-term treatment. Conclusion A complete profile of quasispecies reconstructed from in-depth shotgun sequencing may has important implications for enhancing clinical decision in adjusting antiviral therapy timely.

Introduction Hepatitis B virus (HBV) causes a complex and persistent infection with a major impact on patients health. Viral-genome sequencing can provide valuable information for characterizing virus genotype, infection dynamics and drug and vaccine resistance. Areas covered This article reviews the current literature to describe the next-generation sequencing progress that facilitated a more comprehensive study of HBV quasispecies in diagnosis and clinical monitoring. Expert opinion HBV variability plays a key role in liver disease progression and treatment efficacy. Second-generation sequencing improved the sensitivity for detecting and quantifying mutations, mixed genotypes and viral recombination. Third-generation sequencing enables the analysis of the entire HBV genome, although the high error rate limits its use in clinical practice.

The most insidious pandemic of modern life does not arise from an infectious agent but rather from malnutrition. With its global incidence tripling over the past three decades, obesity is a major public health concern. Obesity’s rising prevalence has also illuminated its impact on communicable diseases. Following the 2009 H1N1 influenza A virus pandemic, obesity was identified as a risk factor for increased disease severity and mortality in infected individuals. Obesity causes a chronic state of meta-inflammation with systemic implications for immunity, including delayed antiviral responses to influenza virus infection, poor recovery, and impaired immunological memory. However, the majority of past work neglected the innate responses, especially the impact on the epithelium. To better understand the innate response to influenza infection, obese animal models and human respiratory epithelial cells were inoculated with influenza viruses to quantify the innate responses. In obese hosts, antiviral effector production, including type I and type III interferons, was consistently blunted across all model systems and linked to a highly oxidative, pro-inflammatory environment. This interferon-deficient state altered the host microenvironment in which influenza replicated, spread, and caused disease. Due to the inherent mutability and segmented nature of influenza virus, this altered microenvironment also had implications for the genetic drift and shift of the viral population. Compared to infection of lean hosts, minor variants rapidly emerged within the viral population in obese hosts which exhibited increased viral replication, enhanced virulence, and antiviral resistance. The delayed action of interferon also permitted greater rates of reassortment. Together, this suggested obesity permitted the emergence of novel, and potentially pathogenic, viral variants. Influenza viruses pose constant seasonal and pandemic threats with high morbidity and mortality suffered by overweight and obese patients. If current trends continue, nearly half of the worldwide population will be obese by 2050. This population will have growing impacts on both non-communicable and communicable diseases and may affect global evolutionary trends of influenza virus.

In large, genetically diverse populations (e.g., RNA viruses), neutral evolution is no less directional than adaptive evolution., Even if they have no impact on phenotype, neutral mutations are not equivalent in the eyes of evolution: A robust neutral variant—one which remains functional after further mutations—is more likely to spread in a large, diverse population than a fragile one. Quasispecies theory shows that the equilibrium frequency of a genotype is proportional to its eigenvector centrality in the neutral network. This paper explores the link between the selection for mutational robustness and the navigability of neutral networks. I show that sequences of neutral mutations follow a “maximal entropy random walk,” a canonical Markov chain on graphs with nonlocal, nondiffusive dynamics. I revisit M. Smith’s word-game model of evolution in this light, finding that the likelihood of certain sequences of substitutions can decrease with the population size. These counterintuitive results underscore the fertility of the interface between evolutionary dynamics, information theory, and physics.

The live attenuated yellow fever (YF) vaccine was developed in the 1930s. Currently, the 17D and 17DD attenuated substrains are used for vaccine production. The 17D strain is used for vaccine production by several countries, while the 17DD strain is used exclusively in Brazil. The cell passages carried out through the seed-lot system of vaccine production influence the presence of quasispecies causing changes in the stability and immunogenicity of attenuated genotypes by increasing attenuation or virulence. Using next-generation sequencing, we carried out genomic characterization and genetic diversity analysis between vaccine lots of the Brazilian YF vaccine, produced by BioManguinhos–Fiocruz, and used during 11 years of vaccination in Brazil. We present 20 assembled and annotated genomes from the Brazilian 17DD vaccine strain, eight single nucleotide polymorphisms and the quasispecies spectrum reconstruction for the 17DD vaccine, through a pipeline here introduced. The V2IDA pipeline provided a relationship between low genetic diversity, maintained through the seed lot system, and the confirmation of genetic stability of lots of the Brazilian vaccine against YF. Our study sets precedents for use of V2IDA in genetic diversity analysis and in silico stability investigation of attenuated viral vaccines, facilitating genetic surveillance during the vaccine production process.

The entire H5N1 highly pathogenic avian influenza viral genomes were identified in the frozen autopsy specimens: the trachea, lung, colon, and intestinal feces from a patient who died of the disease in 2006. Phylogenetic analysis of the viral genomes showed that these viruses belonged to clade 1 and were the reassortants generated from the reassortment of the viruses within the same clade. The sequencing data from the autopsy specimens revealed at least 8 quasispecies of the H5N1 viruses across all 4 specimen types. These sequences were compared to those derived from the virus isolates grown in Madin Darby canine kidney (MDCK) cells. The virus isolates from the trachea, lung, and fecal specimens showed 27 nucleotide substitutions, leading to the changes of 18 amino acid residues. However, there was no change in the amino acid residues that determined the viral virulence. The changes were more commonly observed in the lung, particularly in the HA and NA genes. Our study suggested that the adaptation changes for the viral fitness to survive in a new host species (MDCK cells) might involve many genes, for example, the amino acid substitution 177G or 177W adjacent to the receptor-binding residues in the HA1 globular head and the substitution M315I in PB2. However, a mutation changes near the receptor binding domain may play an important role in determining the cell tropism and is needed to be further explored.

Some of the newly emerging corona viral variants show high numbers of mutations. This is unexpected for a virus with a low mutation rate due to an inherent proof-reading system. Could such a variant arise under very special conditions occurring in a host where the virus replicates and mutates in a rather unlimited fashion, such as in immune compromised patients? The virus was shown to replicate in an immunosuppressed cancer patient for more than 105 days and might be a source of new variants. These patients are asymptomatic and the virus may therefore escape detection and attention and be high-risk. Similarly, HIV-infected individuals may be immunocompromised and support coronavirus replication with increased mutation rates. The patients may promote “within-host evolution”. Some of the viruses present in such a highly mutagenic swarm or quasispecies within one patient may become founders and cause a pandemic by further “between-host evolution”. B.1.1.7 with 23 mutations may be such a case. Immunosuppressed patients can be identified and treated by the synthetic antibody cocktails as passive immunization and kept under control. Immunosuppressed patients can be easily identified and supervised by healthcare workers—once they become aware of the risk—to avoid new variants with pandemic potential.

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., The work at CBMSO was supported by grants SAF2014-52400-R from Ministerio de Economía y Competitividad (MINECO), SAF2017-87846-R and BFU2017-91384-EXP from Ministerio de Ciencia, Innovación y Universidades (MCIU), PI18/00210 from Instituto de Salud Carlos III, S2013/ABI-2906 (PLATESA from Comunidad de Madrid/FEDER), and S2018/BAA-4370 (PLATESA2 from Comunidad de Madrid/FEDER). C.P. is supported by the Miguel Servet program of the Instituto de Salud Carlos III (CPII19/00001), cofinanced by the European Regional Development Fund (ERDF). CIBERehd (Centro de Investigación en Red de Enfermedades Hepáticas y Digestivas) is funded by Instituto de Salud Carlos III. Institutional grants from the Fundación Ramón Areces and Banco Santander to the CBMSO are also acknowledged. The team at CBMSO belongs to the Global Virus Network (GVN). The work in Barcelona was supported by Instituto de Salud Carlos III, cofinanced by the European Regional Development Fund (ERDF) Grant No. PI19/00301 and by the Centro para el Desarrollo Tecnológico Industrial (CDTI) from the MICIU, Grant No. IDI-20200297. Work at CAB was supported by MINECO grant BIO2016-79618R and PID2019-104903RB-I00 (funded by the EU under the FEDER program) and by the Spanish State research agency (AEI) through project number MDM-2017-0737 Unidad de Excelencia “María de Maeztu”-Centro de Astrobiología (CSIC-INTA). C.G.-C. is supported by predoctoral contract PRE2018-083422 from MCIU. B.M.-G. is supported by predoctoral contract PFIS FI19/00119 from Instituto de Salud Carlos III (Ministerio de Sanidad y Consumo), cofinanced by Fondo Social Europeo (FSE).

Hepatitis C is an inflammatory liver disease caused by the single-stranded RNA (ssRNA) hepatitis C virus (HCV). The genetic diversity of the virus and quasispecies produced during replication have resulted in viral resistance to direct-acting antivirals (DAAs) as well as impediments in vaccine development. The recent adaptation of CRISPR-Cas as an alternative antiviral approach has demonstrated degradation of viral nucleic acids in eukaryotes. In particular, the CRISPR-effector Cas13 enzyme has been shown to target ssRNA viruses effectively. In this work, we have employed Cas13a to knockdown HCV in mammalian cells. Using a computational screen, we identified several potential Cas13a target sites within highly conserved regions of the HCV internal ribosomal entry site (IRES). Our results demonstrate significant inhibition of HCV replication as well as translation in huh-7.5 cells with minimal effects on cell viability. These findings were validated using a multi-modality approach involving qRT-PCR, luciferase assay, and MTT cell viability assay. In conclusion, the CRISPR-Cas13a system efficiently targets HCV in vitro, suggesting its potential as a programmable therapeutic antiviral strategy.

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 ≥ 103 IU/mL at delivery. Among them, 75.0% of patients with rtM204I positive had HBV DNA load ≥103 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.

Currently, direct-acting antiviral drugs (DAAs) are widely used as therapeutic methods for hepatitis C virus (HCV)-positive patients, however, patients may experience treatment failure, and the dynamic changes of HCV genomes in these patients are unknown. In this study, three real-world DAAs cohorts were enrolled to observe clinical efficacy. In addition, serum samples from treatment failure patients at baseline and relapse were used to analyze changes of the HCV genomes at near full-length genome level, including resistance-associated variants (RAVs), viral quasispecies diversity and selection analysis. Next-generation sequencing was used as the detection method. The overall sustained virological response at 12 w after the end of treatment was achieved in 91.9% (57/62) of HCV patients, and 3 paired samples obtained from relapsed patients. All the 3 patients harbored baseline NS5A RAVs, the frequency of NS5A RAVs increased in 2 patients and a new NS5A RAV emerged in 1 patient at relapse, and almost all the viral strains existed with NS5A RAVs at relapse. The results of the viral quasispecies diversity analysis revealed that viral quasispecies diversity decreased at relapse compared to baseline, and the results of selection analysis indicated that the virus population experienced a bottleneck phenomenon, recent selective sweep and population expansion or was under purification selection after DAAs treatment. This study indicated that the clinical efficacy was excellent in real-world DAAs cohorts, and the viral strains existed at relapse were selective by DAAs therapy.

Coronaviruses (CoVs) possess the largest and most complex RNA genome (up to 32 kb) that encodes for 16 non-structural proteins regulating RNA synthesis and modification. Coronaviruses are known to infect a wide range of mammalian and avian species causing remarkably diverse disease syndromes. Variable tissue tropism and the ability to easily cross interspecies barriers are the well-known characteristics of certain CoVs. The 21st century epidemics of severe acute respiratory CoV (SARS-CoV), Middle East respiratory CoV and the ongoing SARS-CoV-2 pandemic further highlight these characteristics and emphasize the relevance of CoVs to the global public health. Bovine CoVs (BCoVs) are betacoronaviruses associated with neonatal calf diarrhea, and with winter dysentery and shipping fever in older cattle. Of interest, no distinct genetic or antigenic markers have been identified in BCoVs associated with these distinct clinical syndromes. In contrast, like other CoVs, BCoVs exist as quasispecies. Besides cattle, BCoVs and bovine-like CoVs were identified in various domestic and wild ruminant species (water buffalo, sheep, goat, dromedary camel, llama, alpaca, deer, wild cattle, antelopes, giraffes, and wild goats), dogs and humans. Surprisingly, bovine-like CoVs also cannot be reliably distinguished from BCoVs using comparative genomics. Additionally, there are historical examples of zoonotic transmission of BCoVs. This article will discuss BCoV pathogenesis, epidemiology, interspecies transmission, immune responses, vaccines, and diagnostics.

Norovirus (NoV) is the leading cause of viral-related diarrhea in cancer patients, in whom it can be chronic, contributing to decreased quality of life, interruption of cancer care, malnutrition, and altered mucosal barrier function. Immunosuppressed cancer patients shed NoV for longer periods of time than immunocompetent hosts, favoring quasispecies development and emergence of novel NoV variants. While nucleic acid amplification tests (NAATs) for NoV diagnosis have revolutionized our understanding of NoV burden of disease, not all NAATs provide information on viral load or infecting genotype. There is currently no effective antiviral or vaccine for chronic NoV infections. Screening for inhibitors of NoV replication in intestinal organoid culture models and creation of NoV-specific adoptive T cells are promising new strategies to develop treatments for chronic NoV in immunosuppressed patients. Herein we summarize data on the epidemiology, clinical manifestations, diagnostic challenges, and treatment of NoV infection in patients with cancer.

Ebola virus remains a global threat to public health and biosecurity, yet we still know relatively little about its pathogenesis and the complications that arise following recovery. With nearly 20,000 survivors from the 2013–2016 West African outbreak, as well as over 1,000 survivors of the recent outbreak in the DRC, we must consider the consequences of virus persistence and recrudescent disease, even if they are rare., Ebola virus (EBOV) is responsible for numerous devastating outbreaks throughout Africa, including the 2013–2016 West African outbreak as well as the two recent outbreaks in the Democratic Republic of the Congo (DRC), one of which is ongoing. Although EBOV disease (EVD) has typically been considered a highly lethal acute infection, increasing evidence suggests that the virus can persist in certain immune-privileged sites and occasionally lead to EVD recrudescence. Little is understood about the processes that contribute to EBOV persistence and recrudescence, in part because of the rarity of these phenomena but also because of the absence of an animal model that recapitulates them. Here, we describe a case of EBOV persistence associated with atypical EVD in a nonhuman primate (NHP) following inoculation with EBOV and treatment with an experimental monoclonal antibody cocktail. Although this animal exhibited only mild signs of acute EVD, it developed severe disease 2 weeks later and succumbed shortly thereafter. Viremia was undetectable at the time of death, despite abundant levels of viral RNA in most tissues, each of which appeared to harbor a distinct viral quasispecies. Remarkably, sequence analysis identified a single mutation in glycoprotein (GP) that not only resisted antibody-mediated neutralization but also increased viral growth kinetics and virulence. Overall, this report represents the most thoroughly characterized case of atypical EVD in an NHP described thus far, and it provides valuable insight into factors that may contribute to EBOV persistence and recrudescent disease.

Viral quasispecies dynamics between pre- and post-nucleos(t)ide analog (NA) therapy remains unclear. This study aimed to investigate the HBV quasispecies evolution and its relationship with durability of off-therapy responses in HBeAg-negative chronic hepatitis B (CHB) patients who stopped NA therapy. Fifty-four HBeAg-negative CHB patients who stopped NAs, including 19 virological controllers (VC) who maintained serum HBV DNA

Host immune response and viral factors are involved in disease progression in patients with chronic hepatitis B virus (HBV) infection. However, the relationship between HBV quasispecies and liver fibrosis progression remains unclear. In this study, 447 patients with chronic HBV infection, including 239 with chronic hepatitis B (CHB), 104 with liver cirrhosis (LC) and 104 with hepatocellular carcinoma (HCC) were enrolled. The 239 CHB patients were divided into groups F1, F2, and F3 according to liver fibrosis score. Four fragments of the HBV genome were determined and analyzed using next-generation sequencing. Specific mutations, such as A1762T, G1764A and G1896A, in the BCP/PC region were more common in patients with advanced liver disease and formed the majority of the viral quasispecies pool in patients with LC and HCC. The viral complexity and diversity increased as the fibrosis progressed, especially in patients with CHB who were comparable in age but at different stages of fibrosis. Patients with early-stage fibrosis experienced higher purifying selection pressure in the four sequenced regions, whereas different protein-coding region experienced different negative selection with disease progression. HBV quasispecies diversity may increase fibrosis progression in CHB patients with aging under immune selection.

All life forms on earth ultimately descended from a primordial population dubbed the last universal common ancestor or LUCA via Darwinian evolution. Extant living systems share two salient functional features, a metabolism extracting and transforming energy required for survival, and an evolvable, informational polymer – the genome – conferring heredity. Genome replication invariably generates essential and ubiquitous genetic parasites. Here we model the energetic, replicative conditions of LUCA-like organisms and their parasites, as well as adaptive problem solving of host-parasite pairs. We show using the Lotka-Volterra equations that three host-parasite pairs – individually a unit of a host and a parasite that is itself parasitized – are sufficient for robust and stable homeostasis, forming a life cycle. This catalytic life cycle efficiently captures, channels and transforms energy, enabling dynamic host survival and adaptation. We propose a Malthusian fitness model for an original quasispecies evolving through a host-parasite life cycle.

It is estimated that up to 80% of infections caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) are asymptomatic and asymptomatic patients can still effectively transmit the virus and cause disease. While much of the effort has been placed on decoding single nucleotide variation in SARS-CoV-2 genomes, considerably less is known about their transcript variation and any correlation with clinical severity in human hosts, as defined here by the presence or absence of symptoms. To assess viral genomic signatures of disease severity, we conducted a systematic characterization of SARS-CoV-2 transcripts and genetic variants in 81 clinical specimens collected from symptomatic and asymptomatic individuals using multi-scale transcriptomic analyses including amplicon-seq, short-read metatranscriptome and long-read Iso-seq. Here we show a highly coordinated and consistent pattern of sgRNA expression from individuals with robust SARS-CoV-2 symptomatic infection and their expression is significantly repressed in the asymptomatic infections. We also observe widespread inter- and intra-patient variants in viral RNAs, known as quasispecies frequently found in many RNA viruses. We identify unique sets of deletions preferentially found primarily in symptomatic individuals, with many likely to confer changes in SARS-CoV-2 virulence and host responses. Moreover, these frequently occurring structural variants in SARS-CoV-2 genomes serve as a mechanism to further induce SARS-CoV-2 proteome complexity. Our results indicate that differential sgRNA expression and structural mutational burden are highly correlated with the clinical severity of SARS-CoV-2 infection. Longitudinally monitoring sgRNA expression and structural diversity could further guide treatment responses, testing strategies, and vaccine development. Wong and Ngan et al. characterize the expression and structural variation of SARS-CoV-2 subgenomic RNAs (sgRNAs) in diagnostic specimens from symptomatic and asymptomatic COVID-19 patients. In a series of genomic and transcriptomic analyses, the authors observe reduced sgRNA expression and a distinct set of structural deletions in asymptomatic infections compared with symptomatic infections. Infection with SARS-CoV-2 can lead to symptoms of different severity, with some individuals remaining entirely asymptomatic but still responsible for much of the viral transmission. Here, we sought to identify markers of the severity of symptoms of SARS-CoV-2 infection, as defined by the presence or absence of symptoms. We used a combination of methods to study SARS-CoV-2 genes and their readouts, known as transcripts, in clinical samples from people with symptomatic and asymptomatic infections. We demonstrate that transcripts responsible for making viral proteins are seen at lower levels in asymptomatic infections. We also identify structural changes in the viral genes and transcripts that potentially influence the host’s response to the virus. Our study defines potential markers of symptom severity that may ultimately guide risk mitigation and testing strategies.

International audience; Since the identification of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) as the etiological agent of the current COVID-19 pandemic, a rapid and massive effort has been made to obtain the genomic sequences of this virus to monitor (in near real time) the phylodynamic and diversity of this new pathogen. However, less attention has been given to the assessment of intra-host diversity. RNA viruses such as SARS-CoV-2 inhabit the host as a population of variants called quasispecies. We studied the quasispecies diversity in four of the main SARS-CoV-2 genes (ORF1a, ORF1b, S and N genes), using a dataset consisting of 210 next-generation sequencing (NGS) samples collected between January and early April of 2020 in the State of Victoria, Australia. We found evidence of quasispecies diversity in 68% of the samples, 76% of which was nonsynonymous variants with a higher density in the spike (S) glycoprotein and ORF1a genes. About one-third of the nonsynonymous intra-host variants were shared among the samples, suggesting host-to-host transmission. Quasispecies diversity changed over time. Phylogenetic analysis showed that some of the intra-host single-nucleotide variants (iSNVs) were restricted to specific lineages, highlighting their potential importance in the epidemiology of this virus. A greater effort must be made to determine the magnitude of the genetic bottleneck during transmission and the epidemiological and/or evolutionary factors that may play a role in the changes in the diversity of quasispecies over time.

Hypomyelination is a rare consequence of in utero bovine viral diarrhoea virus (BVDV) infection. We describe a BVDV outbreak in a naïve Holstein dairy herd in northern Italy, with an unusually high prevalence of calves with neurological signs, generalised tremors and ataxia. Histological analysis showed that hypomyelination was the predominant lesion and that the most typical BVDV neuropathological findings (e.g. cerebellar hypoplasia) were absent. Virological and molecular analyses showed that non-cytopathic BVDV genotype 1b was associated with the calves’ neurological signs and excluded other viruses responsible for congenital infection or neurological disorders. Whole-genome sequencing of BVDVs from the brain of a calf with neurological signs and the whole blood of a persistently infected herd-mate with no such sign showed >99.7 % sequence identity. Analysis of the quasispecies distribution revealed the greatest variation rates in regions coding for the structural proteins E1 and E2. Variation was slightly greater in the brain- than in the blood-derived sequence and occurred at different sites, suggesting the occurrence of distinct evolutionary processes in the two persistently infected calves. Molecular characterisation of BVDV genomes from five other calves with neurological signs from the same farm confirmed that the E1 and E2 regions were the most variable. Several factors, including genetic variability and host factors, appear to have contributed to the observed unique BVDV disease phenotype, characterised by hypomyelination and neurological signs.

A married couple who developed hepatitis was referred to our hospital. The husband was diagnosed with acute hepatitis C virus (HCV) infection based on HCV RNA positivity and seroconversion to HCV antibody. The wife was also diagnosed with HCV-related hepatitis; however, she could not be confirmed to have acute hepatitis due to the lack of information on her HCV negativity just before this event. The HCV strains recovered from the couple were genotype 2b and shared 100% identities within the 5’-untranslated region-core region sequence (655 nucleotides/nt) and non-structural (NS)5B region sequence (502 nt). The amplified hypervariable region 1 (HVR-1) sequence indicated that all 10 clones from the wife shared 100% identity and were identical to 3 of 10 heterogeneous clones (separable into 4 groups) from the husband. The husband had a history of intravenous drug use. These results suggested that one of four quasispecies 2b HCV strains was transmitted from the husband to the wife, with the husband being the infectious source for acute HCV infection in the wife, most likely via sexual intercourse. A sequence analysis of the HCV genomes and the further comparison of the HVR-1 amino acid sequence variability may be useful for defining the infectious source of HCV, especially in couples or cluster cases.

The ongoing COVID-19 pandemic, caused by SARS-CoV-2, constitutes a tremendous global health issue. Continuous monitoring of the virus has become a cornerstone to make rational decisions on implementing societal and sanitary measures to curtail the virus spread. Additionally, emerging SARS-CoV-2 variants have increased the need for genomic surveillance to detect particular strains because of their potentially increased transmissibility, pathogenicity and immune escape. Targeted SARS-CoV-2 sequencing of diagnostic and wastewater samples has been explored as an epidemiological surveillance method for the competent authorities. Currently, only the consensus genome sequence of the most abundant strain is taken into consideration for analysis, but multiple variant strains are now circulating in the population. Consequently, in diagnostic samples, potential co-infection(s) by several different variants can occur or quasispecies can develop during an infection in an individual. In wastewater samples, multiple variant strains will often be simultaneously present. Currently, quality criteria are mainly available for constructing the consensus genome sequence, and some guidelines exist for the detection of co-infections and quasispecies in diagnostic samples. The performance of detection and quantification of low-frequency variants using whole genome sequencing (WGS) of SARS-CoV-2 remains largely unknown. Here, we evaluated the detection and quantification of mutations present at low abundances using the mutations defining the SARS-CoV-2 lineage B.1.1.7 (alpha variant) as a case study. Real sequencing data were in silico modified by introducing mutations of interest into raw wild-type sequencing data, or by mixing wild-type and mutant raw sequencing data, to construct mixed samples subjected to WGS using a tiling amplicon-based targeted metagenomics approach and Illumina sequencing. As anticipated, higher variation and lower sensitivity were observed at lower coverages and allelic frequencies. We found that detection of all low-frequency variants at an abundance of 10, 5, 3, and 1%, requires at least a sequencing coverage of 250, 500, 1500, and 10,000×, respectively. Although increasing variability of estimated allelic frequencies at decreasing coverages and lower allelic frequencies was observed, its impact on reliable quantification was limited. This study provides a highly sensitive low-frequency variant detection approach, which is publicly available at https://galaxy.sciensano.be, and specific recommendations for minimum sequencing coverages to detect clade-defining mutations at certain allelic frequencies. This approach will be useful to detect and quantify low-frequency variants in both diagnostic (e.g., co-infections and quasispecies) and wastewater [e.g., multiple variants of concern (VOCs)] samples.

Successful therapeutics and vaccines for coronavirus disease 2019 (COVID-19) have harnessed the immune response to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Evidence that SARS-CoV-2 exists as locally evolving variants suggests that immunological differences may impact the effectiveness of antibody-based treatments such as convalescent plasma and vaccines. Considering that near-sourced convalescent plasma likely reflects the antigenic composition of local viral strains, we hypothesize that convalescent plasma has a higher efficacy, as defined by death within 30 days of transfusion, when the convalescent plasma donor and treated patient were in close geographic proximity. Results of a series of modeling techniques applied to approximately 28,000 patients from the Expanded Access to Convalescent Plasma program (ClinicalTrials.gov number: NCT04338360) support this hypothesis. This work has implications for the interpretation of clinical studies, the ability to develop effective COVID-19 treatments, and, potentially, for the effectiveness of COVID-19 vaccines as additional locally-evolving variants continue to emerge., Regional differences in SARS-CoV-2 variants may affect treatment outcome. Here, the authors show that near-sourced convalescent plasma has higher efficacy, as defined by death within 30 days of transfusion, than plasma sourced more than 150 miles away.

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 antiviral activities at later stages of infection. Here, we characterized the IFN-I sensitivity of 500 clonally derived HIV-1 isolates from the plasma and CD4+ T cells of 26 individuals sampled longitudinally after transmission or after antiretroviral therapy (ART) and analytical treatment interruption. We determined the concentration of IFNα2 and IFNβ that reduced viral replication in vitro by 50% (IC50) and found consistent changes in the sensitivity of HIV-1 to IFN-I inhibition both across individuals and over time. Resistance of HIV-1 isolates to IFN-I was uniformly high during acute infection, decreased in all individuals in the first year after infection, was reacquired concomitant with CD4+ T cell loss, and remained elevated in individuals with accelerated disease. HIV-1 isolates obtained by viral outgrowth during suppressive ART were relatively IFN-I sensitive, resembling viruses circulating just before ART initiation. However, viruses that rebounded after treatment interruption displayed the highest degree of IFNα2 and IFNβ resistance observed at any time during the infection course. These findings indicate a dynamic interplay between host innate responses and the evolving HIV-1 quasispecies, with the relative contribution of IFN-I to HIV-1 control affected by both ART and analytical treatment interruption. Although elevated at transmission, host innate pressures are the highest during viral rebound, limiting the viruses that successfully become reactivated from latency to those that are IFN-I resistant.

Mathematical modeling studies are useful for an improved understanding of the dynamics of human immunodeficiency virus type 1 (HIV-1) infection dynamics. Mutations are considered to facilitate the escape of HIV-1 from the immune system and anti-HIV-1 drugs. Thus, quasispecies, populations of closely related virus variants, are important in HIV-1 infection. Recently, a multiscale model incorporating both intracellular viral genome replication and extracellular viral infection was developed to study viral dynamics. A multiscale model has also been developed for analyzing HIV-1 infection. However, the multiscale model could not reveal high quasispecies diversity. Herein, I present a new quasi-steady-state multiscale model of HIV-1 infection that incorporates a quasi-steady state among the virion, uninfected cell, and virion–uninfected cell complex similar to that among the enzyme, substrate, and enzyme–substrate complex in enzyme kinetics. This new model can reveal stable quasispecies better than the previously described multiscale model. The model and the strategy to calculate the model might be useful in the field of demographics where multiscale models are used. For HIV-1 infection, the model might be useful for analyzing and predicting the breakthrough of drug-resistant mutants in the future.

SummaryBackgroundDetection of incident hepatitis C virus (HCV) infections is crucial for identification of outbreaks and development of public health interventions. However, there is no single diagnostic assay for distinguishing recent and persistent HCV infections. HCV exists in each infected host as a heterogeneous population of genomic variants, whose evolutionary dynamics remain incompletely understood. Genetic analysis of such viral populations can be applied to the detection of incident HCV infections and used to understand intra-host viral evolution.MethodsWe studied intra-host HCV populations sampled using next-generation sequencing from 98 recently and 256 persistently infected individuals. Genetic structure of the populations was evaluated using 245,878 viral sequences from these individuals and a set of selected parameters measuring their diversity, topological structure, complexity, strength of selection, epistasis, evolutionary dynamics, and physico-chemical properties.FindingsDistributions of the viral population parameters differ significantly between recent and persistent infections. A general increase in viral genetic diversity from recent to persistent infections is frequently accompanied by decline in genomic complexity and increase in structuredness of the HCV population, likely reflecting a high level of intra-host adaptation at later stages of infection. Using these findings, we developed a Machine Learning classifier for the infection staging, which yielded a detection accuracy of 95.22%, thus providing a higher accuracy than other genomic-based models.InterpretationThe detection of a strong association between several HCV genetic factors and stages of infection suggests that intra-host HCV population develops in a complex but regular and predictable manner in the course of infection. The proposed models may serve as a foundation of cyber-molecular assays for staging infection, that could potentially complement and/or substitute standard laboratory assays.FundingAZ and PS were supported by NIH grant 1R01EB025022. PIB was supported by GSU MBD fellowship.

Summary Background & Aims We aimed to investigate how viral quasispecies of the HBV whole genome evolves and diversifies in response to HBeAg seroconversion and viral control utilising next-generation sequencing (NGS). Methods Fifty HBeAg-positive chronic hepatitis B patients, including 18 treatment-naïve and 32 interferon (IFN)-treated individuals, were recruited. Serial HBV whole genomes in serum were analysed by NGS to determine sequence characteristics and viral quasispecies. Results HBV quasispecies diversity, measured by nucleotide diversity, was negatively correlated with viral load and hepatitis activity. Spontaneous HBeAg seroconverters exhibited significantly greater viral quasispecies diversity than treatment-naïve non-seroconverters from >1 year before seroconversion (0.0112 vs. 0.0060, p 1 year after seroconversion (0.0103 vs. 0.0068, p, Graphical abstract, Highlights • Deep sequencing of whole HBV genome uncovers the quasispecies changes in chronic hepatitis B patients. • The nucleotide diversity of HBV negatively correlates with viraemia during HBeAg loss/seroconversion. • Viral quasispecies diversity is greater in spontaneous HBeAg seroconverters before and after seroconversion than in treatment-naïve non-seroconverters. • Responders to IFN have greater viral quasispecies diversity than non-responders at 24 weeks after treatment. • The genome positions of non-synonymous intra-host single nucleotide variants (iSNVs) of HBV tend to be located at possible T cell epitopes.

Hepatitis A virus (HAV) is a waterborne pathogen of public health importance. In South Africa (SA), unique HAV subgenotype IB strains have been detected in surface and wastewater samples, as well as on fresh produce at the point of retail. However, due to the use of molecular-based assays, the infectivity of the detected strains was unknown. Considering the potential shift of HAV endemicity from high to intermediate, which could increase the risk of severe symptomatic disease, this study investigated the identity of HAV strains detected before and after viability treatment of selected wastewater discharge samples. For one year, 118 samples consisting of sewage, treated wastewater discharge and downstream dam water were collected from five wastewater treatment plants (WWTP 1, 2, 3, 4 and 5). Unique HAV IB strains were detected in samples from all five WWTPs, with 11 of these strains carrying amino acid mutations at the immunodominant and neutralisation epitopes. A quasispecies dynamic of HAV has also been detected in sewage samples. The subsequent application of viability PCR revealed that potentially infectious HAV strains were discharged from WWTP 1, 2, 4 and 5 into the dam. Therefore, there is a potential risk of HAV exposure to communities using water sources downstream the WWTPs.

In diverse viral infections, the production of excess viral particles containing only viral glycoproteins (subviral particles or SVP) is commonly observed and is a commonly evolved mechanism for immune evasion. In hepatitis B virus (HBV) infection, spherical particles contain the hepatitis B surface antigen, outnumber infectious virus 10 000-100 000 to 1, and have diverse inhibitory effects on the innate and adaptive immune response, playing a major role in the chronic nature of HBV infection. The current goal of therapies in development for HBV infection is a clinical outcome called functional cure, which signals a persistent and effective immune control of the infection. Although removal of spherical SVP (and the HBsAg they carry) is an important milestone in achieving functional cure, this outcome is rarely achieved with current therapies due to distinct mechanisms for assembly, secretion, and persistence of SVP, which are poorly targeted by direct acting antivirals or immunotherapies. In this Review, the current understanding of the distinct mechanisms involved in the production and persistence of spherical SVP in chronic HBV infection and their immunoinhibitory activity will be reviewed as well as current therapies in development with the goal of clearing spherical SVP and achieving functional cure.

Molecular investigation of primary HIV infections (PHI) is crucial to describe current dynamics of HIV transmission. Aim of the study was to investigate HIV transmission clusters (TC) in PHI referred during the years 2013&ndash, 2020 to the National Institute for Infectious Diseases in Rome (INMI), that is the Lazio regional AIDS reference centre, and factors possibly associated with inclusion in TC. These were identified by phylogenetic analysis, based on population sequencing of pol, a more in depth analysis was performed on TC of B subtype, using ultra-deep sequencing (UDS) of env. Of 270 patients diagnosed with PHI during the study period, 229 were enrolled (median follow-up 168 (IQR 96&ndash, 232) weeks). Median age: 39 (IQR 32&ndash, 48) years, 94.8% males, 86.5% Italians, 83.4% MSM, 56.8% carrying HIV-1 subtype B. Of them, 92.6% started early treatment within a median of 4 (IQR 2&ndash, 7) days after diagnosis, median time to sustained suppression was 20 (IQR 8&ndash, 32) weeks. Twenty TC (median size 3, range 2&ndash, 9 individuals), including 68 patients, were identified. A diagnosis prior to 2015 was the unique factor associated with inclusion in a TC. Added value of UDS was the identification of shared quasispecies components in transmission pairs within TC.

Summary In this study, we analyzed norovirus (NoV) evolution in sequential samples of six chronically infected patients. The capsid gene was amplified from stool samples, and deep sequencing was performed. The role of amino acid flexibility in structural changes and ligand binding was studied with molecular dynamics (MD) simulations. Concentrations of capsid-specific antibodies increased in sequential sera. Capsid sequences accumulated mutations during chronic infection, particularly in the surface-exposed antigenic epitopes A, D, and E. The number of quasispecies increased in infections lasting for >1 month. Interestingly, high genetic complexity and distances were followed by ongoing NoV replication, whereas lower genetic complexity and distances preceded cure. MD simulation revealed that surface-exposed amino acid substitutions of the P2 domain caused fluctuation of blockade epitopes. In conclusion, the capsid protein accumulates numerous mutations during chronic infection; however, only those on the protein surface change the protein structure substantially and may lead to immune escape., Graphical abstract, Highlights • Noroviruses accumulated numerous mutations during chronic infection • After 1 month, intraindividual norovirus populations became genetically diverse • High genetic complexity and genetic distances were associated with prolonged infections, Molecular modeling; Virology

Although previous studies have analyzed cross-level CRF01_AE viral genomic data in populations, less is known about intrapatient viral evolutionary dynamics during antiretroviral therapy (ART) failure. We longitudinally sampled plasma and peripheral blood mononuclear cells (PBMC) at different time points from one human immunodeficiency virus type 1 infected patient. The evolution of viral quasispecies was inferred from viral phylogenies. Before treatment, no drug-resistant mutations were found in this patient's plasma, and all viruses had C-C chemokine receptor type 5 (CCR5) tropism. Two months after treatment, the majority of the virus population in plasma and PBMC were drug resistant and X4-tropic. By 5 months after treatment, the viral load increased significantly, and viruses reversed tropism from X4 to R5 in plasma and PBMC. During treatment failure, the effective population of the pol DNA reservoir in PBMC remained stable, whereas the env DNA reservoir increased. The effective population of the R5 tropism virus increased more rapidly than that of the X4 tropism virus. The ratio of non-synonymous to synonymous substitutions in the env gene of R5 tropism virus (0.43) was lower than X4 tropism (0.52). However, four env positive selection sites were identified in R5 tropism viruses (HXB2: 364, 398, 399, and 400) but none were identified in X4 tropism viruses. Our data demonstrated the different intrapatient evolutionary dynamics patterns of env and pol genes in an individual who experienced periods of ART failure. Our findings also suggest the importance of the R5 tropism virus in the DNA reservoir during ART failure.

Prolonged infection of SARS-CoV-2 represents a challenge to the development of effective public health policies to control the COVID-19 pandemic. The reason why some people have persistent infection and how the virus survives for so long are still not fully understood. For this reason, we aimed to investigate the intra-host evolution of SARS-CoV-2 during persistent infection. Thirty-three patients who remained RT-PCR positive in the nasopharynx for at least 16 days were included in this study. Complete SARS-CoV-2 sequences were obtained for each patient at two time points. Phylogenetic, populational, and computational analysis of viral sequences confirmed persistent infection with evidence for a transmission cluster in health care professionals that shared the same workplace. A high number of missense variants targeting crucial structural and non-structural proteins such as Spike and Helicase was found. Interestingly, longitudinal acquisition of substitutions in Spike protein mapped many SARS-CoV-2 predicted T cell epitopes. Furthermore, the mutational profiles observed were suggestive of RNA editing enzyme activities, indicating innate immune mechanisms of the host cell. Viral quasispecies analysis corroborates persistent infection mainly by increasing richness and nucleotide diversity over time. Altogether, our findings highlight a dynamic and complex landscape of host and pathogen interaction during persistent infection suggesting that the host’s innate immunity shapes the increase of intra-host diversity with possible implications for therapeutic strategies and public health decisions during the COVID-19 pandemic.

Introduction: Hand, foot and mouth disease caused by enterovirus A71 (EV-A71) is more frequently associated with neurological complications and deaths compared to other enteroviruses. Areas covered: In this review, the authors put together current understanding of the neuropathogenesis of EV-A71 based on various clinical, human and animal model studies. The authors discuss the important advancements in virus entry, virus dissemination and neuroinvasion. The authors highlight the role of host immune system, host genetic factors, viral quasispecies and heparan sulfate in EV-A71 neuropathogenesis.Expert opinion: Comparison of EV-A71 with EV-D68 and PV shows similarity in primary target sites and dissemination to the central nervous system. More research is needed to understand cellular tropisms, persistence of EV-A71 and other possible invasion routes. EV-A71 infection has varied clinical manifestations which may be attributed to multiple receptors usage. Future development of antivirals and vaccines should target neurotropic enteroviruses. Repurposing drug and immunomodulators used in combination could reduce the severity of EV-A71 infection. Only a few drugs have been tested in clinical trials, and in the absence of antiviral and vaccines (in most parts of the world), active virus surveillance, good hand hygiene and physical distancing should be advocated. A better understanding of EV-A71 neuropathogenesis is critical for antiviral and multivalent vaccines development.

Background The genome of SARS-CoV-2 is susceptible to mutations during viral replication due to the errors generated by RNA-dependent RNA polymerases. These mutations enable the SARS-CoV-2 to evolve into new strains. Viral quasispecies emerge from de novo mutations that occur in individual patients. In combination, these sets of viral mutations provide distinct genetic fingerprints that reveal the patterns of transmission and have utility in contact tracing. Methods Leveraging thousands of sequenced SARS-CoV-2 genomes, we performed a viral pangenome analysis to identify conserved genomic sequences. We used a rapid and highly efficient computational approach that relies on k-mers, short tracts of sequence, instead of conventional sequence alignment. Using this method, we annotated viral mutation signatures that were associated with specific strains. Based on these highly conserved viral sequences, we developed a rapid and highly scalable targeted sequencing assay to identify mutations, detect quasispecies variants, and identify mutation signatures from patients. These results were compared to the pangenome genetic fingerprints. Results We built a k-mer index for thousands of SARS-CoV-2 genomes and identified conserved genomics regions and landscape of mutations across thousands of virus genomes. We delineated mutation profiles spanning common genetic fingerprints (the combination of mutations in a viral assembly) and a combination of mutations that appear in only a small number of patients. We developed a targeted sequencing assay by selecting primers from the conserved viral genome regions to flank frequent mutations. Using a cohort of 100 SARS-CoV-2 clinical samples, we identified genetic fingerprints consisting of strain-specific mutations seen across populations and de novo quasispecies mutations localized to individual infections. We compared the mutation profiles of viral samples undergoing analysis with the features of the pangenome. Conclusions We conducted an analysis for viral mutation profiles that provide the basis of genetic fingerprints. Our study linked pangenome analysis with targeted deep sequenced SARS-CoV-2 clinical samples. We identified quasispecies mutations occurring within individual patients and determined their general prevalence when compared to over 70,000 other strains. Analysis of these genetic fingerprints may provide a way of conducting molecular contact tracing.

“Cross-protection”, a nearly 100 years-old virological term, is suggested to be changed to “close protection”. Evidence for the need of such change has accumulated over the past six decades from the laboratory experiments and field tests conducted by plant pathologists and plant virologists working with different plant viruses, and, in particular, from research on Citrus tristeza virus (CTV). A direct confirmation of such close protection came with the finding that “pre-immunization” of citrus plants with the variants of the T36 strain of CTV but not with variants of other virus strains was providing protection against a fluorescent protein-tagged T36-based recombinant virus variant. Under natural conditions close protection is functional and is closely associated both with the conservation of the CTV genome sequence and prevention of superinfection by closely similar isolates. It is suggested that the mechanism is primarily directed to prevent the danger of virus population collapse that could be expected to result through quasispecies divergence of large RNA genomes of the CTV variants continuously replicating within long-living and highly voluminous fruit trees. This review article provides an overview of the CTV cross-protection research, along with a discussion of the phenomenon in the context of the CTV biology and genetics.

The 2014&ndash, 2016 Zika virus (ZIKV) epidemic in the Americas resulted in large deposits of next-generation sequencing data from clinical samples. This resource was mined to identify emerging mutations and trends in mutations as the outbreak progressed over time. Information on transmission dynamics, prevalence, and persistence of intra-host mutants, and the position of a mutation on a protein were then used to prioritize 544 reported mutations based on their ability to impact ZIKV phenotype. Using this criteria, six mutants (representing naturally occurring mutations) were generated as synthetic infectious clones using a 2015 Puerto Rican epidemic strain PRVABC59 as the parental backbone. The phenotypes of these naturally occurring variants were examined using both cell culture and murine model systems. Mutants had distinct phenotypes, including changes in replication rate, embryo death, and decreased head size. In particular, a NS2B mutant previously detected during in vivo studies in rhesus macaques was found to cause lethal infections in adult mice, abortions in pregnant females, and increased viral genome copies in both brain tissue and blood of female mice. Additionally, mutants with changes in the region of NS3 that interfaces with NS5 during replication displayed reduced replication in the blood of adult mice. This analytical pathway, integrating both bioinformatic and wet lab experiments, provides a foundation for understanding how naturally occurring single mutations affect disease outcome and can be used to predict the of severity of future ZIKV outbreaks. To determine if naturally occurring individual mutations in the Zika virus epidemic genotype affect viral virulence or replication rate in vitro or in vivo, we generated an infectious clone representing the epidemic genotype of stain Puerto Rico, 2015. Using this clone, six mutants were created by changing nucleotides in the genome to cause one to two amino acid substitutions in the encoded proteins. The six mutants we generated represent mutations that differentiated the early epidemic genotype from genotypes that were either ancestral or that occurred later in the epidemic. We assayed each mutant for changes in growth rate, and for virulence in adult mice and pregnant mice. Three of the mutants caused catastrophic embryo effects including increased embryonic death or significant decrease in head diameter. Three other mutants that had mutations in a genome region associated with replication resulted in changes in in vitro and in vivo replication rates. These results illustrate the potential impact of individual mutations in viral phenotype.