Increasing evolution, prevalence, and outbreaks for rift valley fever virus in the process of breaking geographical barriers

BACKGROUND

Rift valley fever (RVF) is listed as one of prioritized diseases by WHO. This study aims to describe RVF virus' landscape distribution globally, and to insight dynamics change of its evolution, prevalence, and outbreaks in the process of breaking geographical barriers.

METHODS

A systematic literature review and meta-analyses was conducted to estimate RVF prevalence by hosts using a random-effect model. Molecular clock-based phylogenetic analyses were performed to estimate RVF virus nucleotide substitution rates using nucleotide sequences in NCBI database. RVF virus prevalence, nucleotide substitution rates, and outbreaks were compared before and after breaking geographical barriers twice, respectively.

RESULTS

RVF virus was reported from 26 kinds of hosts covering 48 countries from 1930 to 2022. Since RVF broke geographical barriers, (1) nucleotide substitution rates significantly increased after firstly spreading out of Africa in 2000, (2) prevalence in humans significantly increased from 1.92 % (95 % CI: 0.86-3.25 %) to 3.03 % (95 % CI: 2.09-4.12 %) after it broke Sahara Desert geographical barriers in 1977, and to 5.24 % (95 % CI: 3.81-6.82 %) after 2000, (3) RVF outbreaks in humans and the number of wildlife hosts presented increasing trends. RVF virus spillover may exist between bats and humans, and accelerate viral substitution rates in humans. During outbreaks, the RVF virus substitution rates accelerated in humans. 60.00 % RVF outbreaks occurred 0-2 months after floods and (or) heavy rainfall.

CONCLUSION

RVF has the increasing risk to cause pandemics, and global collaboration on "One Health" is needed to prevent potential pandemics.

Genomic Evolution and Recombination Dynamics of Human Adenovirus D Species: Insights from Comprehensive Bioinformatic Analysis

Human adenoviruses (HAdVs) can infect various epithelial mucosal cells, ultimately causing different symptoms in infected organ systems. With more than 110 types classified into seven species (A-G), HAdV-D species possess the highest number of viruses and are the fastest proliferating. The emergence of new adenovirus types and increased diversity are driven by homologous recombination (HR) between viral genes, primarily in structural elements such as the penton base, hexon and fiber proteins, and the E1 and E3 regions. A comprehensive analysis of the HAdV genome provides valuable insights into the evolution of human adenoviruses and identifies genes that display high variation across the entire genome to determine recombination patterns. Hypervariable regions within genetic sequences correlate with functional characteristics, thus allowing for adaptation to new environments and hosts. Proteotyping of newly emerging and already established adenoviruses allows for prediction of the characteristics of novel viruses. HAdV-D species evolved in a direction that increased diversity through gene recombination. Bioinformatics analysis across the genome, particularly in highly variable regions, allows for the verification or re-evaluation of recombination patterns in both newly introduced and pre-existing viruses, ultimately aiding in tracing various biological traits such as virus tropism and pathogenesis. Our research does not only assist in predicting the emergence of new adenoviruses but also offers critical guidance in regard to identifying potential regulatory factors of homologous recombination hotspots.

Successful treatment of SARS-CoV-2 in an immunocompromised patient with persistent infection for 245 days: A case report

Background

Immunocompromised patients receiving B-cell-depleting therapies are at increased risk of persistent SARS-CoV-2 infection, with many experiencing fatal outcomes. We report a successful outcome in a patient with rheumatoid arthritis (RA) on rituximab diagnosed with COVID-19 in July 2020 with persistent infection for over 245 days.

Results

The patient received numerous treatment courses for persistent COVID-19 infection, including remdesivir, baricitinib, immunoglobulin and high doses of corticosteroids followed by a prolonged taper due to persistent respiratory symptoms and cryptogenic organizing pneumonia. Her clinical course was complicated by Pseudomonas aeruginosa sinusitis with secondary bacteremia, and cytomegalovirus (CMV) viremia and pneumonitis. SARS-CoV-2 positive RNA samples were extracted from two nasopharyngeal swabs and sequenced using targeted amplicon Next-Generation Sequencing which were analyzed for virus evolution over time. Viral sequencing indicated lineage B.1.585.3 SARS-CoV-2 accumulated Spike protein mutations associated with immune evasion and resistance to therapeutics. Upon slowly decreasing the patient's steroids, she had resolution of her symptoms and had a negative nasopharyngeal SARS-CoV-2 PCR and serum CMV PCR in March 2021.

Conclusion

A patient with RA on B-cell depleting therapy developed persistent SARS-CoV-2 infection allowing for virus evolution and had numerous complications, including viral and bacterial co-infections with opportunistic pathogens. Despite intra-host evolution with a more immune evasive SARS-CoV-2 lineage, it was cleared after 245 days with reconstitution of the patient's immune system.

Whole genome sequencing reveals insights into hepatitis E virus genome diversity, and virus compartmentalization in chronic hepatitis E

BACKGROUND

Hepatitis E virus (HEV) is a leading cause of acute hepatitis and can cause chronic infections in immunocompromised patients. Although HEV infections can be treated with ribavirin, antiviral efficacy is hampered by resistance mutations, normally detected by virus sequencing.

OBJECTIVES

High-throughput sequencing (HTS) allows for cost-effective complete viral genome sequencing. This enables the discovery and delineation of new subtypes, and revised the recognition of quasispecies and putative resistance mutations. However, HTS is challenged by factors including low viral load, sample degradation, high host background, and high viral diversity.

STUDY DESIGN

We apply complete genome sequencing strategies for HEV, including a targeted enrichment approach. These approaches were used to investigate sequence diversity in HEV RNA-positive animal and human samples and intra-host diversity in a chronically infected patient.

RESULTS

Here, we describe the identification of potential novel subtypes in a blood donation (genotype 3) and in an ancient livestock sample (genotype 7). In a chronically infected patient, we successfully investigated intra-host virus diversity, including the presence of ribavirin resistance mutations. Furthermore, we found convincing evidence for HEV compartmentalization, including the central nervous system, in this patient.

CONCLUSIONS

Targeted enrichment of viral sequences enables the generation of complete genome sequences from a variety of difficult sample materials. Moreover, it enables the generation of greater sequence coverage allowing more advanced analyses. This is key for a better understanding of virus diversity. Investigation of existing ribavirin resistance, in the context of minorities or compartmentalization, may be critical in treatment strategies of HEV patients.

Molecular insights into the SARS-CoV-2 Omicron variant from Bangladesh suggest diverse and continuous evolution

The study analyzed the molecular dynamics of the circulating SARS-CoV-2 Omicron variant from its identification in November 2021 to January 2023. The SARS-CoV-2 sequences from Bangladesh revealed three distinct waves of the Omicron variant. More than 50 sub-lineages of Omicron variant were introduced into the country, with the majority belonging to the major lineages of BA.1-like (24.91%), BA.2-like (43.35%), BA.5-like (5.76%), XBB (10.47%), and "Others and Unassigned" (18.64%). Furthermore, the relative frequencies over time revealed that Omicron lineages existed for a short period of time before being replaced by other sub-lineages. Many potential mutations were found in the receptor binding domain of the Spike protein including G339D/H, S371 L/F, K417 N, T478K, E484A, Q493R, Q498R, and N501Y. In conclusion, the SARS-CoV-2 Omicron variant from Bangladesh showed diverse genetic features and continuous evolution. Therefore, the choice of vaccine and monitoring of hospitalized patients is important alongside genetic characterization of the circulating SARS-CoV-2.

Phylodynamics of dengue virus 2 in Nicaragua leading up to the 2019 epidemic reveals a role for lineage turnover

BACKGROUND

Dengue is a mosquito-borne viral disease posing a significant threat to public health. Dengue virus (DENV) evolution is often characterized by lineage turnover, which, along with ecological and immunological factors, has been linked to changes in dengue phenotype affecting epidemic dynamics. Utilizing epidemiologic and virologic data from long-term population-based studies (the Nicaraguan Pediatric Dengue Cohort Study and Nicaraguan Dengue Hospital-based Study), we describe a lineage turnover of DENV serotype 2 (DENV-2) prior to a large dengue epidemic in 2019. Prior to this epidemic, Nicaragua had experienced relatively low levels of DENV transmission from 2014 to 2019, a period dominated by chikungunya in 2014/15 and Zika in 2016.

RESULTS

Our phylogenetic analyses confirmed that all Nicaraguan DENV-2 isolates from 2018 to 2019 formed their own clade within the Nicaraguan lineage of the Asian/American genotype. The emergence of the new DENV-2 lineage reflects a replacement of the formerly dominant clade presiding from 2005 to 2009, a lineage turnover marked by several shared derived amino acid substitutions throughout the genome. To elucidate evolutionary drivers of lineage turnover, we performed selection pressure analysis and reconstructed the demographic history of DENV-2. We found evidence of adaptive evolution by natural selection at the codon level as well as in branch formation.

CONCLUSIONS

The timing of its emergence, along with a statistical signal of adaptive evolution and distinctive amino acid substitutions, the latest in the NS5 gene, suggest that this lineage may have increased fitness relative to the prior dominant DENV-2 strains. This may have contributed to the intensity of the 2019 DENV-2 epidemic, in addition to previously identified immunological factors associated with pre-existing Zika virus immunity.

Dengue virus serotype 2 genotype III evolution during the 2019 outbreak in Mato Grosso, Midwestern Brazil

DENV-2 was the main responsible for a 70% increase in dengue incidence in Brazil during 2019. That year, our metagenomic study by Illumina NextSeq on serum samples from acute febrile patients (n = 92) with suspected arbovirus infection, sampled in 22 cities of the state of Mato Grosso (MT), in the middle west of Brazil, revealed eight complete genomes and two near-complete sequences of DENV-2 genotype III, one Human parvovirus B19 genotype I (5,391 nt) and one Coxsackievirus A6 lineage D (4,514 nt). These DENV-2 sequences share the aminoacidic identities of BR4 lineage on E protein domains I, II and III, and were included in a clade with sequences of the same lineage circulating in the southeast of Brazil in the same year. Nevertheless, 11/34 non-synonymous mutations are unique to three strains inthis study, distributed in the E (n = 6), NS3 (n = 2) and NS5 (n = 3) proteins. Other 14 aa changes on C (n = 1), E (n = 3), NS1 (n = 2), NS2A (n = 1) and NS5 (n = 7) were first reported in a genotype III lineage, having been already reported only in other DENV-2 genotypes. All 10 sequences have mutations in the NS5 protein (14 different aa changes). Nine E protein aa changes found in two sequences, six of which are unique, are in the ectodomain; where the E:M272T change is on the hinge of the E protein at domain II, in a region critical for the anchoring to the host cell receptor. The NS5:G81R mutation, in the methyltransferase domain, was found in one strain of this study. Altogether, these data points to an important evolution of DENV-2 genotype III lineage BR4 during this outbreak in 2019 in MT. Genomic surveillance is essential to detect virus etiology and evolution, possibly related to immune evasion and viral fitness changes leading to future novel outbreaks.

Genomic evolution of bovine viral diarrhea virus based on complete genome and individual gene analyses

Bovine viral diarrhea virus (BVDV) genome consists of a single-stranded, positive-sense RNA with high genetic diversity. In the last years, significant progress has been achieved in BVDV knowledge evolution through phylodynamic analysis based on the partial 5'UTR sequences, whereas a few studies have used other genes or the complete coding sequence (CDS). However, no research has evaluated and compared BVDV evolutionary history based on the complete genome (CG), CDS, and individual genes. In this study, phylodynamic analyses were carried out with BVDV-1 (Pestivirus A) and BVDV-2 (Pestivirus B) CG sequences available on the GenBank database and each genomic region: CDS, UTRs, and individual genes. In comparison to the CG, the estimations for both BVDV species varied according to the dataset used, pointing out the importance of considering the analyzed genomic region when concluding. This study may provide new insight into BVDV evolution history while highlighting the need to increase the available BVDV CG sequences to perform more comprehensive phylodynamic studies in the future.

A classification algorithm based on dynamic ensemble selection to predict mutational patterns of the envelope protein in HIV-infected patients

BACKGROUND

Therapeutics against the envelope (Env) proteins of human immunodeficiency virus type 1 (HIV-1) effectively reduce viral loads in patients. However, due to mutations, new therapy-resistant Env variants frequently emerge. The sites of mutations on Env that appear in each patient are considered random and unpredictable. Here we developed an algorithm to estimate for each patient the mutational state of each position based on the mutational state of adjacent positions on the three-dimensional structure of the protein.

METHODS

We developed a dynamic ensemble selection algorithm designated k-best classifiers. It identifies the best classifiers within the neighborhood of a new observation and applies them to predict the variability state of each observation. To evaluate the algorithm, we applied amino acid sequences of Envs from 300 HIV-1-infected individuals (at least six sequences per patient). For each patient, amino acid variability values at all Env positions were mapped onto the three-dimensional structure of the protein. Then, the variability state of each position was estimated by the variability at adjacent positions of the protein.

RESULTS

The proposed algorithm showed higher performance than the base learner and a panel of classification algorithms. The mutational state of positions in the high-mannose patch and CD4-binding site of Env, which are targeted by multiple therapeutics, was predicted well. Importantly, the algorithm outperformed other classification techniques for predicting the variability state at multi-position footprints of therapeutics on Env.

CONCLUSIONS

The proposed algorithm applies a dynamic classifier-scoring approach that increases its performance relative to other classification methods. Better understanding of the spatiotemporal patterns of variability across Env may lead to new treatment strategies that are tailored to the unique mutational patterns of each patient. More generally, we propose the algorithm as a new high-performance dynamic ensemble selection technique.

Jawless vertebrates do not escape retrovirus infection

On occasion, retroviruses infect the genome of germline cell, forming endogenous retroviruses (ERVs), which provide molecular fossils for studying the deep evolution of retroviruses. While ERVs have been extensively characterized in the genomes of jawed vertebrates, much remains contentious and unexplored about the diversity and evolution of ERVs within jawless vertebrates. Here, we report the discovery of a novel ERV lineage, designated EbuERVs, in the genome of a hagfish Eptatretus burgeri. Phylogenetic analyses show that EbuERVs pertain to epsilon-retroviruses and might have derived from cross-species transmission from jawed vertebrates. EbuERVs are estimated to have invaded in the hagfish genome at least tens of millions of years ago. Evolutionary dynamics analyses indicate that EbuERVs might have experienced one proliferation peak and have been not active in transposition anymore. However, some EbuERVs can transcribe in embryo and might serve as lncRNA. Overall, these findings expand the distribution of retroviruses from jawed vertebrates to jawless vertebrates.

Genome-based survey of the SARS-CoV-2 BF.7 variant from Asia

The SARS-CoV-2 BF.7 variant represents one of the most recent subvariant under monitoring. At the beginning of the 2023 if caused several concerns especially in Asia because of a resurge in COVID-19 cases. Here we perform a genome-based integrative approach on SARS-CoV-2 BF.7 in order to shed light on this emerging lineage and produce some consideration on its real dangerousness. Both genetic and structural data suggest that this new variant currently does not show evidence of an high expansion capability. It is very common in Asia, but it appears less virulent than other Omicron variants as proved by its relatively low evolutionary rate (5.62 x 10^-4 subs/sites/years). The last plateau has been reached around December 14, 2022 and then the genetic variability, and thus the viral population size, no longer increased. As already seen for several previous variants, the features that may be theoretically related to advantages are due to genetic drift that allows to the virus a constant adaptability to the host, but is not strictly connected to a fitness advantage. These results have further pointed that the genome-based monitoring must continue uninterruptedly in order to be prepared and well documented on the real situation. This article is protected by copyright. All rights reserved.

Differential Neural Networks Prediction Using Slow and Fast Hybrid Learning: Application to Prognosis of Infectionsand Deaths of COVID-19 Dynamics

This essay discusses a potential method for predicting the behavior of various physical processes and uses the COVID-19 outbreak to demonstrate its applicability. This study assumes that the current data set reflects the output of a dynamic system that is governed by a nonlinear ordinary differential equation. This dynamic system may be described by a Differential Neural Network (DNN) with time-varying weights matrix parameters. A new hybrid learning scheme based on the decomposition of the signal to be predicted. The decomposition considers the slow and fast components of the signal which is more natural to signals such as the ones corresponding to the number of infected and deceased patients who suffered of COVID 2019 sickness. The paper results demonstrate the recommended method offers competitive performance (70 days of COVID prediction) in comparison to similar studies.

The emergence of an imported variant of dengue virus serotype 2 in the Jazan region, southwestern Saudi Arabia

BACKGROUND

Dengue virus (DENV) infection is a global economic and public health concern, particularly in tropical and subtropical countries where it is endemic. Saudi Arabia has seen an increase in DENV infections, especially in the western and southwestern regions. This study aims to investigate the genetic variants of DENV-2 that were circulating during a serious outbreak in Jazan region in 2019.

METHODS

A total of 482 serum samples collected during 2019 from Jazan region were tested with reverse transcription-polymerase chain reaction (RT-PCR) to detect and classify DENV; positive samples underwent sequencing and bioinformatics analyses.

RESULTS

Out of 294 positive samples, type-specific RT-PCR identified 58.8% as DENV-2 but could not identify 41.2%. Based on sequencing and bioinformatics analyses, the samples tested PCR positive in the first round but PCR negative in the second round were found to be imported genetic variant of DENV-2. The identified DENV-2 imported variant showed similarities to DENV-2 sequences reported in Malaysia, Singapore, Korea and China. The results revealed the imported genetic variant of DENV-2 was circulating in Jazan region that was highly prevalent and it was likely a major factor in this outbreak.

CONCLUSIONS

The emergence of imported DENV variants is a serious challenge for the dengue fever surveillance and control programmes in endemic areas. Therefore, further investigations and continuous surveillance of existing and new viral strains in the region are warranted.

The evolution, characterization and phylogeography of avian influenza H9N2 viruses from India

The low pathogenic avian influenza H9N2 virus is a significant zoonotic agent and contributes genes to highly pathogenic avian influenza (HPAI) viruses. H9N2 viruses are prevalent in India with a reported human case. We elucidate the spatio-temporal origins of the H9N2 viruses from India. A total of 30H9N2 viruses were isolated from poultry and environmental specimens (years 2015-2020). Genome sequences of H9N2 viruses (2003-2020) from India were analyzed, revealing several substitutions. We found five reassortant genotypes. The HA, NA and PB2 genes belonged to the Middle-Eastern B sublineage; NP and M to the classical G1 lineage; PB1, PA and NS showed resemblance to genes from either HPAI-H7N3/H5N1 viruses. Molecular clock and phylogeography revealed that the introduction of all the genes to India took place around the year 2000. This is the first report of the genesis and evolution of the H9N2 viruses from India, and highlights the need for surveillance.

Mycoviruses as a part of the global virome: Diversity, evolutionary links and lifestyle

Knowledge of mycovirus diversity, evolution, horizontal gene transfer and shared ancestry with viruses infecting distantly related hosts, such as plants and arthropods, has increased vastly during the last few years due to advances in the high throughput sequencing methodologies. This also has enabled the discovery of novel mycoviruses with previously unknown genome types, mainly new positive and negative single-stranded RNA mycoviruses ((+) ssRNA and (-) ssRNA) and single-stranded DNA mycoviruses (ssDNA), and has increased our knowledge of double-stranded RNA mycoviruses (dsRNA), which in the past were thought to be the most common viruses infecting fungi. Fungi and oomycetes (Stramenopila) share similar lifestyles and also have similar viromes. Hypothesis about the origin and cross-kingdom transmission events of viruses have been raised and are supported by phylogenetic analysis and by the discovery of natural exchange of viruses between different hosts during virus-fungus coinfection in planta. In this review we make a compilation of the current information on the genome organization, diversity and taxonomy of mycoviruses, discussing their possible origins. Our focus is in recent findings suggesting the expansion of the host range of many viral taxa previously considered to be exclusively fungal, but we also address factors affecting virus transmissibility and coexistence in single fungal or oomycete isolates, as well as the development of synthetic mycoviruses and their use in investigating mycovirus replication cycles and pathogenicity.

The Role of Extensive Recombination in the Evolution of Geminiviruses

Mutation, recombination and pseudo-recombination are the major forces driving the evolution of viruses by the generation of variants upon which natural selection, genetic drift and gene flow can act to shape the genetic structure of viral populations. Recombination between related virus genomes co-infecting the same cell usually occurs via template swapping during the replication process and produces a chimeric genome. The family Geminiviridae shows the highest evolutionary success among plant virus families, and the common presence of recombination signatures in their genomes reveals a key role in their evolution. This review describes the general characteristics of members of the family Geminiviridae and associated DNA satellites, as well as the extensive occurrence of recombination at all taxonomic levels, from strain to family. The review also presents an overview of the recombination patterns observed in nature that provide some clues regarding the mechanisms involved in the generation and emergence of recombinant genomes. Moreover, the results of experimental evolution studies that support some of the conclusions obtained in descriptive or in silico works are summarized. Finally, the review uses a number of case studies to illustrate those recombination events with evolutionary and pathological implications as well as recombination events in which DNA satellites are involved.

Virus Evolution on Fitness Landscapes

The landscape paradigm is revisited in the light of evolution in simple systems. A brief overview of different classes of fitness landscapes is followed by a more detailed discussion of the RNA model, which is currently the only evolutionary model that allows for a comprehensive molecular analysis of a fitness landscape. Neutral networks of genotypes are indispensable for the success of evolution. Important insights into the evolutionary mechanism are gained by considering the topology of sequence and shape spaces. The dynamic concept of molecular quasispecies is viewed in the light of the landscape paradigm. The distribution of fitness values in state space is mirrored by the population structures of mutant distributions. Two classes of thresholds for replication error or mutations are important: (i) the-conventional-genotypic error threshold, which separates ordered replication from random drift on neutral networks, and (ii) a phenotypic error threshold above which the molecular phenotype is lost. Empirical landscapes are reviewed and finally, the implications of the landscape concept for virus evolution are discussed.

Ribovirus classification by a polymerase barcode sequence

RNA viruses encoding a polymerase gene (riboviruses) dominate the known eukaryotic virome. High-throughput sequencing is revealing a wealth of new riboviruses known only from sequence, precluding classification by traditional taxonomic methods. Sequence classification is often based on polymerase sequences, but standardised methods to support this approach are currently lacking. To address this need, we describe the polymerase palmprint, a segment of the palm sub-domain robustly delineated by well-conserved catalytic motifs. We present an algorithm, Palmscan, which identifies palmprints in nucleotide and amino acid sequences; PALMdb, a collection of palmprints derived from public sequence databases; and palmID, a public website implementing palmprint identification, search, and annotation. Together, these methods demonstrate a proof-of-concept workflow for high-throughput characterisation of RNA viruses, paving the path for the continued rapid growth in RNA virus discovery anticipated in the coming decade.

Regional mutations in CHIKV-ECSA genomes and detection of other viruses in the serum of acute febrile patients by a metagenomic approach in Mato Grosso, Central-Western Brazil, 2018

Mato Grosso (MT) State is part of central western Brazil and has a tropical permissive environment that favors arbovirus outbreaks. A metagenomic approach was used to identify viral genomes in seven pools of serum from patients (n=65) with acute febrile disease. Seven chikungunya virus (CHIKV) genomes were determined, showing four amino acid changes found only in CHIKV genomes obtained in MT since 2018: nsP2:T31I, nsP3: A388V, E3:T201I and E3:H57R, in addition to other mutations in E1, nsP2 and nsP4. Six parvovirus B19 (B19V) genotype I genomes (4771-5131 nt) showed four aa alterations (NS1:N473D, R579Q; VP1:I716T; and 11 kDa:V44A) compared to most similar B19V from the USA. Coinfection between CHIKV and B19V was evidenced in 22/65 (33.8%) patients by RT‒PCR and PCR, respectively. Other viruses found in these pools include human pegivirus C, torque teno virus 3, an unclassified TTV and torque teno mini virus. Metagenomics represents a useful approach to detect viruses in the serum of acute febrile patients suspected of arbovirus disease.

The origins and molecular evolution of SARS-CoV-2 lineage B.1.1.7 in the UK

The first SARS-CoV-2 variant of concern (VOC) to be designated was lineage B.1.1.7, later labelled by the World Health Organization as Alpha. Originating in early autumn but discovered in December 2020, it spread rapidly and caused large waves of infections worldwide. The Alpha variant is notable for being defined by a long ancestral phylogenetic branch with an increased evolutionary rate, along which only two sequences have been sampled. Alpha genomes comprise a well-supported monophyletic clade within which the evolutionary rate is typical of SARS-CoV-2. The Alpha epidemic continued to grow despite the continued restrictions on social mixing across the UK and the imposition of new restrictions, in particular, the English national lockdown in November 2020. While these interventions succeeded in reducing the absolute number of cases, the impact of these non-pharmaceutical interventions was predominantly to drive the decline of the SARS-CoV-2 lineages that preceded Alpha. We investigate the only two sampled sequences that fall on the branch ancestral to Alpha. We find that one is likely to be a true intermediate sequence, providing information about the order of mutational events that led to Alpha. We explore alternate hypotheses that can explain how Alpha acquired a large number of mutations yet remained largely unobserved in a region of high genomic surveillance: an under-sampled geographical location, a non-human animal population, or a chronically infected individual. We conclude that the latter provides the best explanation of the observed behaviour and dynamics of the variant, although the individual need not be immunocompromised, as persistently infected immunocompetent hosts also display a higher within-host rate of evolution. Finally, we compare the ancestral branches and mutation profiles of other VOCs and find that Delta appears to be an outlier both in terms of the genomic locations of its defining mutations and a lack of the rapid evolutionary rate on its ancestral branch. As new variants, such as Omicron, continue to evolve (potentially through similar mechanisms), it remains important to investigate the origins of other variants to identify ways to potentially disrupt their evolution and emergence.

Viral metagenomics unveiled extensive communications of viruses within giant pandas and their associated organisms in the same ecosystem

Cross-species transmission events were commonplace, with numerous cases of host-switching during the viral evolutionary history, but relatively little evidence for onward transmission in different species living in the same ecosystem. For understanding the communications of viruses in giant pandas (Ailuropoda melanoleuca) and their associated organisms, based on a large size of samples (N = 2305) collected between 2015 and 2020 from giant panda (N = 776) and other four giant panda-associated organisms in the same ecosystem, red pandas (N = 700), stray cats (N = 32), wild rats (N = 42), and mosquitoes (N = 755), viromics was used for the virus identification and subsequent virus traceability. The results showed that a feline panleukopenia virus (FPV) was found in giant pandas with clinical signs of vomiting and mild diarrhea. Meanwhile, the same FPV strain was also prevalent in the healthy red panda (Ailurus fulgens) population. From the viromes of the five different organisms, 250 virus genomes were determined. Our data revealed that besides FPV, other putative pathogenic viruses, such as red panda amdoparvoviruses (RPAVs) and Getah viruses (GETVs) were responsible for previous disease or death of some red pandas. We also demonstrated that a number of viruses were involved in potential interspecies jumping events between giant pandas and their associated species. Collectively, our results shed light on the genetic diversity and relationship of diverse viral pathogens in 'Giant pandas-Associated animals-Arthropods' and report some cases of possible viral host-switching among these host species living in the same ecosystem.

On the intrinsic nature of viral pathogenesis: The assumption of a Darwinian paradigm to describe COVID-19 pandemic

Our hypothesis about evolution of the COVID-19 pandemic foresees an inverse relation between infectivity (R0) and lethality (L) of SARS-CoV-2. The above parameters are driven by a continuing mutation process granting the virus a clear survival advantage over virulence. For interpreting this relation we adopted a simple equation, R0 × L ≈ k, by which R0 and L depend upon a constant k, that corresponds to an intrinsic property of the viral species involved. The hypothesis was verified by following changes of the R0 and L terms of the formula in the different variants of SARS-CoV-2 that progressively appeared. A further validation came when the equation was applied to pandemic and epidemic influenza type A viruses, Ebola virus and measles virus. We believe this equation that considers virus biology in Darwinian terms could be extremely useful to better face infectious viral threats and validate virus-host molecular interactions relevant to viral pathogenesis.

A simple model for how the risk of pandemics from different virus families depends on viral and human traits

Different virus families, like influenza or corona viruses, exhibit characteristic traits such as typical modes of transmission and replication as well as specific animal reservoirs in which each family of viruses circulate. These traits of genetically related groups of viruses influence how easily an animal virus can adapt to infect humans, how well novel human variants can spread in the population, and the risk of causing a global pandemic. Relating the traits of virus families to their risk of causing future pandemics, and identification of the key time scales within which public health interventions can control the spread of a new virus that could cause a pandemic, are obviously significant. We address these issues using a minimal model whose parameters are related to characteristic traits of different virus families. A key trait of viruses that "spillover" from animal reservoirs to infect humans is their ability to propagate infection through the human population (fitness). We find that the risk of pandemics emerging from virus families characterized by a wide distribution of the fitness of spillover strains is much higher than if such strains were characterized by narrow fitness distributions around the same mean. The dependences of the risk of a pandemic on various model parameters exhibit inflection points. We find that these inflection points define informative thresholds. For example, the inflection point in variation of pandemic risk with time after the spillover represents a threshold time beyond which global interventions would likely be too late to prevent a pandemic.

Multiple novel filamentous phages detected in the cloacal swab samples of birds using viral metagenomics approach

Members of the family Inoviridae (inoviruses) are characterized by their unique filamentous morphology and infection cycle. The viral genome of inovirus is able to integrate into the host genome and continuously releases virions without lysing the host, establishing chronic infection. A large number of inoviruses have been obtained from microbial genomes and metagenomes recently, but putative novel inoviruses remaining to be identified. Here, using viral metagenomics, we identified four novel inoviruses from cloacal swab samples of wild and breeding birds. The circular genome of those four inoviruses are 6732 to 7709 nt in length with 51.4% to 56.5% GC content and encodes 9 to 13 open reading frames, respectively. The zonula occludens toxin gene implicated in the virulence of pathogenic host bacteria were identified in all four inoviruses and shared the highest amino acid sequences identity (< 37.3%) to other reference strains belonging to different genera of the family Inoviridae and among themselves. Phylogenetic analysis indicated that all the four inoviruses were genetically far away from other strains belonging to the family Inoviridae and formed an independent clade. According to the genetic distance-based criteria, all the four inoviruses identified in the present study respectively belong to four novel putative genera in the family Inoviridae.

Molecular archeology of human viruses

The evolution of human-virus associations is usually reconstructed from contemporary patterns of genomic diversity. An intriguing, though still rarely implemented, alternative is to search for the genetic material of viruses in archeological and medical archive specimens to document evolution as it happened. In this chapter, we present lessons from ancient DNA research and incorporate insights from virology to explore the potential range of applications and likely limitations of archeovirological approaches. We also highlight the numerous questions archeovirology will hopefully allow us to tackle in the near future, and the main expected roadblocks to these avenues of research.

Deciphering transmission dynamics and spillover of avian influenza viruses from avian species to swine populations globally

Genome sequences of eleven avian influenza virus (AIV) subtypes have been reported in swine populations from seven countries until August 2020. To unravel the transmission dynamics and spillover events of AIVs from avian reservoirs to swine, full-length hemagglutinin (HA) sequences of AIV subtypes (n = 11) reported from various avian species and swine were retrieved from the ‘Influenza Research Database’. Phylogenetic analysis identified closely related avian and swine AIV sequences suggesting potential spillover events from multiple domestic and wild avian species, including chicken, duck, pigeon, goose, quail, and aquatic birds to swine. Furthermore, N-linked glycosylation analysis of these closely related AIV sequences supported the possibility of multiple spillover events of highly pathogenic H5N1 and low pathogenic H9N2 viruses from various avian species to swine. The principal coordinate analysis further validated these findings for H5N1 and H9N2 viruses; however, spillover events of the other nine AIV subtypes were limited. Interestingly, the presence of potential mammalian adaptation markers, particularly in some of the swine H5N1, H7N9, and H9N2 viruses, suggested that these viruses may have already adapted in swine. The occurrence and circulation of these AIVs in swine, especially the H5N1 and H9N2 viruses with numerous spillover events from the avian reservoirs to swine, pose a significant threat in terms of their reassortment with endemic swine viruses or circulating human influenza viruses within the swine which may facilitate the emergence of a novel influenza virus strain with pandemic potential.

N-terminal domain mutations of the spike protein are structurally implicated in epitope recognition in emerging SARS-CoV-2 strains

During the past two years, the world has been ravaged by a global pandemic caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Acquired mutations in the SARS-CoV-2 genome affecting virus infectivity and/or immunogenicity have led to a number of novel strains with higher transmissibility compared to the original Wuhan strain. Mutations in the receptor binding domain (RBD) of the SARS-CoV-2 spike protein have been extensively studied in this context. However, mutations and deletions within the N-terminal domain (NTD) located adjacent to the RBD are less studied. Many of these are found within certain β sheet-linking loops, which are surprisingly long in SARS-CoV-2 in comparison to SARS-CoV and other related β coronaviruses. Here, we perform a structural and epidemiological study of novel strains carrying mutations and deletions within these loops. We identify short and long-distance interactions that stabilize the NTD loops and form a critical epitope that is essential for the recognition by a wide variety of neutralizing antibodies from convalescent plasma. Among the different mutations/deletions found in these loops, Ala 67 and Asp 80 mutations as well as His 69/Val 70 and Tyr 144 deletions have been identified in different fast-spreading strains. Similarly, deletions in amino acids 241-243 and 246-252 have been found to affect the network of NTD loops in strains with high transmissibility. Our structural findings provide insight regarding the role of these mutations/deletions in altering the epitope structure and thus affecting the immunoreactivity of the NTD region of spike protein.

Early dynamic changes of quasispecies in the reverse transcriptase region of hepatitis B virus in telbivudine treatment

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.

Ancient Gene Capture and Recent Gene Loss Shape the Evolution of Orthopoxvirus-Host Interaction Genes

The survival of viruses depends on their ability to resist host defenses and, of all animal virus families, the poxviruses have the most antidefense genes. Orthopoxviruses (ORPV), a genus within the subfamily Chordopoxvirinae, infect diverse mammals and include one of the most devastating human pathogens, the now eradicated smallpox virus. ORPV encode ∼200 genes, of which roughly half are directly involved in virus genome replication and expression as well as virion morphogenesis. The remaining ∼100 "accessory" genes are responsible for virus-host interactions, particularly counter-defense of innate immunity. Complete sequences are currently available for several hundred ORPV genomes isolated from a variety of mammalian hosts, providing a rich resource for comparative genomics and reconstruction of ORPV evolution. To identify the provenance and evolutionary trends of the ORPV accessory genes, we constructed clusters including the orthologs of these genes from all chordopoxviruses. Most of the accessory genes were captured in three major waves early in chordopoxvirus evolution, prior to the divergence of ORPV and the sister genus Centapoxvirus from their common ancestor. The capture of these genes from the host was followed by extensive gene duplication, yielding several paralogous gene families. In addition, nine genes were gained during the evolution of ORPV themselves. In contrast, nearly every accessory gene was lost, some on multiple, independent occasions in numerous lineages of ORPV, so that no ORPV retains them all. A variety of functional interactions could be inferred from examination of pairs of ORPV accessory genes that were either often or rarely lost concurrently. IMPORTANCE Orthopoxviruses (ORPV) include smallpox (variola) virus, one of the most devastating human pathogens, and vaccinia virus, comprising the vaccine used for smallpox eradication. Among roughly 200 ORPV genes, about half are essential for genome replication and expression as well as virion morphogenesis, whereas the remaining half consists of accessory genes counteracting the host immune response. We reannotated the accessory genes of ORPV, predicting the functions of uncharacterized genes, and reconstructed the history of their gain and loss during the evolution of ORPV. Most of the accessory genes were acquired in three major waves antedating the origin of ORPV from chordopoxviruses. The evolution of ORPV themselves was dominated by gene loss, with numerous genes lost at the base of each major group of ORPV. Examination of pairs of ORPV accessory genes that were either often or rarely lost concurrently during ORPV evolution allows prediction of different types of functional interactions.

Prediction of two novel overlapping ORFs in the genome of SARS-CoV-2

Six candidate overlapping genes have been detected in SARS-CoV-2, yet current methods struggle to detect overlapping genes that recently originated. However, such genes might encode proteins beneficial to the virus, and provide a model system to understand gene birth. To complement existing detection methods, I first demonstrated that selection pressure to avoid stop codons in alternative reading frames is a driving force in the origin and retention of overlapping genes. I then built a detection method, CodScr, based on this selection pressure. Finally, I combined CodScr with methods that detect other properties of overlapping genes, such as a biased nucleotide and amino acid composition. I detected two novel ORFs (ORF-Sh and ORF-Mh), overlapping the spike and membrane genes respectively, which are under selection pressure and may be beneficial to SARS-CoV-2. ORF-Sh and ORF-Mh are present, as ORF uninterrupted by stop codons, in 100% and 95% of the SARS-CoV-2 genomes, respectively.

Phylogenetic and evolutionary analysis of foot-and-mouth disease virus A/ASIA/Sea-97 lineage

Foot-and-mouth disease virus (FMDV) A/ASIA/Sea-97 is a predominant lineage in Southeast Asia and East Asia. However, Sea-97 lineage has not been well studied since its first outbreak in Thailand in 1997. Thus, we conducted phylogenetic and evolutionary analysis of Sea-97 using 224 VP1 sequences of FMDV A/ASIA during 1960 and 2018. Phylogenetic analysis revealed that Sea-97 lineage can be classified into five groups (G1–G5). After the emergence of G2 from G1, the genetic diversity of Sea-97 increased sharply, causing divergence into G3, G4 and G5. During this evolutionary process, Sea-97 lineage, which was initially found only in some countries in Southeast Asia, gradually spread to East Asia. The evolution rate of this lineage was estimated to be 1.2 × 10^–2 substitutions/site/year and there were many differences in amino acid residues compared to vaccine strain. Substitutions at antigenically important sites may affect the efficacy of the vaccine, suggesting the need for appropriate vaccine strains. Our results could provide meaningful information to understand comprehensive characteristic of Sea-97 lineage.

A holistic perspective on herpes simplex virus (HSV) ecology and evolution

Herpes simplex viruses (HSV) cause chronic infection in humans that are characterized by periodic episodes of mucosal shedding and ulcerative disease. HSV causes millions of infections world-wide, with lifelong bouts of viral reactivation from latency in neuronal ganglia. Infected individuals experience different levels of disease severity and frequency of reactivation. There are two distantly related HSV species, with HSV-1 infections historically found most often in the oral niche and HSV-2 infections in the genital niche. Over the last two decades, HSV-1 has emerged as the leading cause of first-episode genital herpes in multiple countries. While HSV-1 has the highest level of genetic diversity among human alpha-herpesviruses, it is not yet known how quickly the HSV-1 viral population in a human host adapts over time, or if there are population bottlenecks associated with viral reactivation and/or transmission. It is also unknown how the ecological environments in which HSV infections occur influence their evolutionary trajectory, or that of co-occurring viruses and microbes. In this review, we explore how HSV accrues genetic diversity within each new infection, and yet maintains its ability to successfully infect most of the human population. A holistic examination of the ecological context of natural human infections can expand our awareness of how HSV adapts as it moves within and between human hosts, and reveal the complexity of these lifelong human-virus interactions. These insights may in turn suggest new areas of exploration for other chronic pathogens that successfully evolve and persist among their hosts.

The Evolution of Severe Acute Respiratory Syndrome Coronavirus-2 during Pandemic and Adaptation to the Host

Severe Acute Respiratory Syndrome Coronavirus-2 is a zoonotic virus with a possible origin in bats and potential transmission to humans through an intermediate host. When zoonotic viruses jump to a new host, they undergo both mutational and natural selective pressures that result in non-synonymous and synonymous adaptive changes, necessary for efficient replication and rapid spread of diseases in new host species. The nucleotide composition and codon usage pattern of SARS-CoV-2 indicate the presence of a highly conserved, gene-specific codon usage bias. The codon usage pattern of SARS-CoV-2 is mostly antagonistic to human and bat codon usage. SARS-CoV-2 codon usage bias is mainly shaped by the natural selection, while mutational pressure plays a minor role. The time-series analysis of SARS-CoV-2 genome indicates that the virus is slowly evolving. Virus isolates from later stages of the outbreak have more biased codon usage and nucleotide composition than virus isolates from early stages of the outbreak.

Experimental evolution of cowpea mild mottle virus reveals recombination-driven reduction in virulence accompanied by increases in diversity and viral fitness

Major themes in pathogen evolution are emergence, evolution of virulence, host adaptation and the processes that underlie them. RNA viruses are of particular interest due to their rapid evolution. The in vivo molecular evolution of an RNA plant virus was demonstrated here using a necrotic isolate of cowpea mild mottle virus (CPMMV) and a susceptible soybean genotype submitted to serial inoculations. We show that the virus lost the capacity to cause necrosis after six passages through the host plant. When a severe bottleneck was imposed, virulence reduction occurred in the second passage. The change to milder symptoms had fitness benefits for the virus (higher RNA accumulation) and for its vector, the whitefly Bemisia tabaci. Genetic polymorphisms were highest in ORF1 (viral replicase) and were independent of the symptom pattern. Recombination was a major contributor to this diversity - even with the strong genetic bottleneck, recombination events and hot spots were detected within ORF1. Virulence reduction was associated with different sites in ORF1 associated to recombination events in both experiments. Overall, the results demonstrate that the reduction in virulence was a consequence of the emergence of new variants, driven by recombination. Besides providing details of the evolutionary mechanisms behind a reduction in virulence and its effect under viral and vector fitness, we propose that this recombination-driven switch in virulence allows the pathogen to rapidly adapt to a new host and, potentially, switch back.

Genomic diversity in a population of Spodoptera frugiperda nucleopolyhedrovirus

Spodoptera frugiperda multiple nucleopolyhedrovirus (SfMNPV) represents a strong candidate to develop environmental-friendly pesticides against the fall armyworm (Spodoptera frugiperda), a widespread pest that poses a severe threat to different crops around the world. To date, SfMNPV genomic diversity of different isolates has been mainly studied by means of restriction pattern analyses and by sequencing of the egt region. Here, the genomic diversity present inside an isolate of SfMNPV was explored using high-throughput sequencing for the first time. We identified 704 intrahost single nucleotide variants, from which 184 are nonsynonymous mutations distributed among 82 different coding sequences. We detected several structural variants affecting SfMNPV genome, including two previously reported deletions inside the egt region. A comparative analysis between polymorphisms present in different SfMNPV isolates and our intraisolate diversity data suggests that coding regions with higher genetic diversity are associated with oral infectivity or unknown functions. In this context, through molecular evolution studies we provide evidence of diversifying selection acting on sf29, a putative collagenase which could contribute to the oral infectivity of SfMNPV. Overall, our results contribute to deepen our understanding of the coevolution between SfMNPV and the fall armyworm and will be useful to improve the applicability of this virus as a biological control agent.

Genetic variation and evolution of foot-and-mouth disease virus serotype A in relation to vaccine matching

Foot-and-mouth disease (FMD) is a severe, highly contagious viral disease that affects a wide variety of domestic and wild cloven-hoofed animals. FMD vaccines can play a vital role in disease control and are very widely used globally each year. However, due to the diversity of FMDV, the choice of FMD vaccine is still a huge challenge. In this study, 45 FMDV/A isolates were phylogenetically categorized into three topotypes: ASIA (n = 31), AFRICA (n = 10), and EURO-SA (n = 4). Three sera collected from vaccinated cattle with FMDV A22/IRQ/24/64, A/IRN/05, and A/ARG/01 were used to evaluate their antigenic relationship (r₁) with the field isolates. The IRQ/24/64 serum demonstrated a 39% (17/44) match (r₁ ≥ 0.3) to the field isolates, whereas IRN/05 serum and ARG/01serum showed an 18% (8/44) and a 2% (1/44) match (r₁ ≥ 0.3) to the field isolates, respectively. The A22/IRQ/24/64 matched with isolates mainly from topotype ASIA, with limited cross-topotype match with isolates from topotypes AFRICA and EURO-SA. However, the A/IRN/05 did not show a cross-topotype match with topotype AFRICA isolates and A/ARG/01 failed to match any isolates from topotypes ASIA and AFRICA. After analyzing the amino acid variation of the known antigenic sites of 45 strains of FMDV/A, it was found that together antigenic sites 1 and 3 contributed about 71% of the amino acid changes to the vaccine evaluated. Based on the capsid sequences, the FMDV/A evolved unequally among topotypes. The topotypes of ASIA and AFRICA evolves faster than that of EURO-SA. The FMDV/A continues to show a high level of genetic diversity driven by a high substitution rate, purifying selection, and positive selection concentrated on antigenic sites or near antigenic sites. The current research shows the challenges of the FMDV/A vaccine selection and emphasizes the importance of continuous monitoring of antigenic evolution for the selection of effective vaccines.

Host-directed editing of the SARS-CoV-2 genome

The extensive sequence data generated from SARS-CoV-2 during the 2020 pandemic has facilitated the study of viral genome evolution over a brief period of time. This has highlighted instances of directional mutation pressures exerted on the SARS-CoV-2 genome from host antiviral defense systems. In this brief review we describe three such human defense mechanisms, the apolipoprotein B mRNA editing catalytic polypeptide-like proteins (APOBEC), adenosine deaminase acting on RNA proteins (ADAR), and reactive oxygen species (ROS), and discuss their potential implications on SARS-CoV-2 evolution.

Potential APOBEC-mediated RNA editing of the genomes of SARS-CoV-2 and other coronaviruses and its impact on their longer term evolution

Members of the APOBEC family of cytidine deaminases show antiviral activities in mammalian cells through lethal editing in the genomes of small DNA viruses, herpesviruses and retroviruses, and potentially those of RNA viruses such as coronaviruses. Consistent with the latter, APOBEC-like directional C→U transitions of genomic plus-strand RNA are greatly overrepresented in SARS-CoV-2 genome sequences of variants emerging during the COVID-19 pandemic. A C→U mutational process may leave evolutionary imprints on coronavirus genomes, including extensive homoplasy from editing and reversion at targeted sites and the occurrence of driven amino acid sequence changes in viral proteins. If sustained over longer periods, this process may account for the previously reported marked global depletion of C and excess of U bases in human seasonal coronavirus genomes. This review synthesizes the current knowledge on APOBEC evolution and function and the evidence of their role in APOBEC-mediated genome editing of SARS-CoV-2 and other coronaviruses.

-

SARS-CoV-2 sequence variants contain an overabundance of C- > U transitions

-

C- > U transitions are the hallmark of the activity of APOBEC cytosine deaminases

-

Further work is needed to determine APOBEC's role in coronavirus evolution

A recent origin of Orf3a from M protein across the coronavirus lineage arising by sharp divergence

The genome of SARS-CoV-2, the coronavirus responsible for the Covid-19 pandemic, encodes a number of accessory genes. The longest accessory gene, Orf3a, plays important roles in the virus lifecycle indicated by experimental findings, known polymorphisms, its evolutionary trajectory and a distinct three-dimensional fold. Here we show that supervised, sensitive database searches with Orf3a detect weak, yet significant and highly specific similarities to the M proteins of coronaviruses. The similarity profiles can be used to derive low-resolution three-dimensional models for M proteins based on Orf3a as a structural template. The models also explain the emergence of Orf3a from M proteins and suggest a recent origin across the coronavirus lineage, enunciated by its restricted phylogenetic distribution. This study provides evidence for the common origin of M and Orf3a families and proposes for the first time a working model for the structure of the universally distributed M proteins in coronaviruses, consistent with the properties of both protein families.

A novel nyavirus lacking matrix and glycoprotein genes from Argas japonicus ticks

Viruses are highly diverse and are the sole agents that can infect organisms in all domains of life. Viruses are defined as capsid-encoding organisms as opposed to ribosome-encoding cellular organisms. However, recent advances in virology indicate the existence of unique viruses that do not meet this basic definition, such as capsidless viruses. During virome analysis of the soft tick Argas japonicus, we identified virus-like sequences closely related to the members of genus Nyavirus (family Nyamiviridae). Further analysis revealed sequences derived from a novel nyavirus that lacks two structural protein genes, matrix (M) and glycoprotein (G). This unique nyavirus is tentatively named Sekira virus (SEKRV). To our knowledge, this is the first study to report a nyavirus deficient in M and G genes in nature. The mechanism of infection, replication, and persistence of SEKRV remain unknown, yet this finding provides new insight into virus evolution and the diverse way of viral life in nature.

Evolution and genetic diversity of SARS-CoV-2 in Africa using whole genome sequences

Background

The ongoing SARS-CoV-2 pandemic was introduced into Africa on 14th February 2020 and has rapidly spread across the continent causing a severe public health crisis and mortality. We investigated the genetic diversity and evolution of this virus during the early outbreak months, between 14th February to 24th April 2020, using whole genome sequences.

Methods

We performed recombination analysis against closely related CoV strains, Bayesian time scaled phylogeny, and investigation of spike protein amino acid mutations.

Results

Recombination signals were observed between the Afr-SARS-CoV-2 sequences and reference sequences within the RdRPs and S genes. The evolutionary rate of the Afr-SARS-CoV-2 was 4.133 × 10⁻⁴ Highest Posterior Density (HPD 4.132 × 10⁻⁴ to 4.134 × 10⁻⁴) substitutions/site/year. The time to most recent common ancestor (TMRCA) of the African strains was December 7th 2019, (95% HPD November 12th 2019-December 29th 2019). The Afr-SARCoV-2 sequences diversified into two lineages A and B, with B being more diverse with multiple sub-lineages confirmed by both maximum clade credibility (MCC) tree and PANGOLIN software. There was a high prevalence of the D614G spike protein amino acid mutation 59/69 (82.61%) among the African strains.

Conclusion

This study has revealed a rapidly diversifying viral population with the G614G spike protein variant dominatinge advocate for up scaling NGS sequencing platforms across Africa to enhance surveillance and aid control effort of SARS-CoV-2 in Africa.

Intra-host evolution of the ssDNA virus tomato severe rugose virus (ToSRV)

To evaluate and quantify the evolutionary dynamics of the bipartite begomovirus tomato severe rugose virus (ToSRV) in a cultivated and a non-cultivated host, plants of tomato and Nicandra physaloides were biolistically inoculated with an infectious clone and systemically infected leaves were sampled at 30, 75 and 120 days after inoculation. Total DNA was extracted and sequenced in the Illumina HiSeq 2000 platform. The datasets were trimmed with the quality score limit set to 0.01, and the assembly was performed using the infectious clone sequence as reference. SNPs were filtered using a minimum p-value of 0.001 and the sum frequencies were used to calculate the deviation from the original clone sequence. Nucleotide substitution rates were calculated for the two DNA components in both hosts: 1.73 × 10^-3 and 3.07 × 10^-4 sub/site/year for the DNA-A and DNA-B, respectively, in N. physaloides, and 8.05 × 10^-4 and 7.02 × 10^-5 sub/site/year the for DNA-A and DNA-B, respectively, in tomato. These values are in the same range of those estimated for viruses with single-stranded RNA genomes and for other begomoviruses. Strikingly, the number of substitutions decreased over time, suggesting the presence of bottlenecks during systemic infection. Determination of Shannon's entropy indicated different patterns of variation in the DNA-A and the DNA-B, suggesting distinct evolutionary forces acting upon each component.

HBV evolution and genetic variability: Impact on prevention, treatment and development of antivirals

Hepatitis B virus (HBV) poses a major global health burden with 260 million people being chronically infected and 890,000 dying annually from complications in the course of the infection. HBV is a small enveloped virus with a reverse-transcribed DNA genome that infects hepatocytes and can cause acute and chronic infections of the liver. HBV is endemic in humans and apes representing the prototype member of the viral family Hepadnaviridae and can be divided into 10 genotypes. Hepadnaviruses have been found in all vertebrate classes and constitute an ancient viral family that descended from non-enveloped progenitors more than 360 million years ago. The de novo emergence of the envelope protein gene was accompanied with the liver-tropism and resulted in a tight virus-host association. The oldest HBV genomes so far have been isolated from human remains of the Bronze Age and the Neolithic (~7000 years before present). Despite the remarkable stability of the hepadnaviral genome over geological eras, HBV is able to rapidly evolve within an infected individual under pressure of the immune response or during antiviral treatment. Treatment with currently available antivirals blocking intracellular replication of HBV allows controlling of high viremia and improving liver health during long-term therapy of patients with chronic hepatitis B (CHB), but they are not sufficient to cure the disease. New therapy options that cover all HBV genotypes and emerging viral variants will have to be developed soon. In addition to the antiviral treatment of chronically infected patients, continued efforts to expand the global coverage of the currently available HBV vaccine will be one of the key factors for controlling the rising global spread of HBV. Certain improvements of the vaccine (e.g. inclusion of PreS domains) could counteract known problems such as low or no responsiveness of certain risk groups and waning anti-HBs titers leading to occult infections, especially with HBV genotypes E or F. But even with an optimal vaccine and a cure for hepatitis B, global eradication of HBV would be difficult to achieve because of an existing viral reservoir in primates and bats carrying closely related hepadnaviruses with zoonotic potential.

Selection and antigenic characterization of immune-escape mutants of H7N2 low pathogenic avian influenza virus using homologous polyclonal sera

Understanding the dynamics of the selection of influenza A immune escape variants by serum antibody is critical for designing effective vaccination programs for animals, especially poultry where large populations have a short generation time and may be vaccinated with high frequency. In this report, immune-escape mutants of A/turkey/New York/4450/1994 H7N2 low pathogenic avian influenza virus, were selected by serially passaging the virus in the presence of continuously increasing concentrations of homologous chicken polyclonal sera. Amino acid mutations were identified by sequencing the parental hemagglutinin (HA) gene and every 10 passages by both Sanger and deep sequencing, and the antigenic distance of the mutants to the parent strain was determined. Progressively, a total of five amino acid mutations were observed over the course of 30 passages. Based on their absence from the parental virus with deep sequencing, the mutations appear to have developed de novo. The antigenic distance between the selected mutants and the parent strain increased as the number of amino acid mutations accumulated and the concentration of antibodies had to be periodically increased to maintain the same reduction in virus titer during selection. This selection system demonstrates how H7 avian influenza viruses behave under selection with homologous sera, and provides a glimpse of their evolutionary dynamics, which can be applied to developing vaccination programs that maximize the effectiveness of a vaccine over time.

Structure and assembly of double-stranded RNA mycoviruses

Mycoviruses are a diverse group that includes ssRNA, dsRNA, and ssDNA viruses, with or without a protein capsid, as well as with a complex envelope. Most mycoviruses are transmitted by cytoplasmic interchange and are thought to lack an extracellular phase in their infection cycle. Structural analysis has focused on dsRNA mycoviruses, which usually package their genome in a 120-subunit T=1 icosahedral capsid, with a capsid protein (CP) dimer as the asymmetric unit. The atomic structure is available for four dsRNA mycovirus from different families: Saccharomyces cerevisiae virus L-A (ScV-L-A), Penicillium chrysogenum virus (PcV), Penicillium stoloniferum virus F (PsV-F), and Rosellinia necatrix quadrivirus 1 (RnQV1). Their capsids show structural variations of the same framework, with asymmetric or symmetric CP dimers respectively for ScV-L-A and PsV-F, dimers of similar domains of a single CP for PcV, or of two different proteins for RnQV1. The CP dimer is the building block, and assembly proceeds through dimers of dimers or pentamers of dimers, in which the genome is packed as ssRNA by interaction with CP and/or viral polymerase. These capsids remain structurally undisturbed throughout the viral cycle. The T=1 capsid participates in RNA synthesis, organizing the viral polymerase (1-2 copies) and a single loosely packaged genome segment. It also acts as a molecular sieve, to allow the passage of viral transcripts and nucleotides, but to prevent triggering of host defense mechanisms. Due to the close mycovirus-host relationship, CP evolved to allocate peptide insertions with enzyme activity, as reflected in a rough outer capsid surface.

Immuno-epidemiology and pathophysiology of coronavirus disease 2019 (COVID-19)

Occasional zoonotic viral attacks on immunologically naive populations result in massive death tolls that are capable of threatening human survival. Currently, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the infectious agent that causes coronavirus disease (COVID-19), has spread from its epicenter in Wuhan China to all parts of the globe. Real-time mapping of new infections across the globe has revealed that variable transmission patterns and pathogenicity are associated with differences in SARS-CoV-2 lineages, clades, and strains. Thus, we reviewed how changes in the SARS-CoV-2 genome and its structural architecture affect viral replication, immune evasion, and transmission within different human populations. We also looked at which immune dominant regions of SARS-CoV-2 and other coronaviruses are recognized by Major Histocompatibility Complex (MHC)/Human Leukocyte Antigens (HLA) genes and how this could impact on subsequent disease pathogenesis. Efforts were also placed on understanding immunological changes that occur when exposed individuals either remain asymptomatic or fail to control the virus and later develop systemic complications. Published autopsy studies that reveal alterations in the lung immune microenvironment, morphological, and pathological changes are also explored within the context of the review. Understanding the true correlates of protection and determining how constant virus evolution impacts on host-pathogen interactions could help identify which populations are at high risk and later inform future vaccine and therapeutic interventions.

Multiple approaches for massively parallel sequencing of SARS-CoV-2 genomes directly from clinical samples

BACKGROUND

COVID-19 (coronavirus disease 2019) has caused a major epidemic worldwide; however, much is yet to be known about the epidemiology and evolution of the virus partly due to the scarcity of full-length SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2) genomes reported. One reason is that the challenges underneath sequencing SARS-CoV-2 directly from clinical samples have not been completely tackled, i.e., sequencing samples with low viral load often results in insufficient viral reads for analyses.

METHODS

We applied a novel multiplex PCR amplicon (amplicon)-based and hybrid capture (capture)-based sequencing, as well as ultra-high-throughput metatranscriptomic (meta) sequencing in retrieving complete genomes, inter-individual and intra-individual variations of SARS-CoV-2 from serials dilutions of a cultured isolate, and eight clinical samples covering a range of sample types and viral loads. We also examined and compared the sensitivity, accuracy, and other characteristics of these approaches in a comprehensive manner.

RESULTS

We demonstrated that both amplicon and capture methods efficiently enriched SARS-CoV-2 content from clinical samples, while the enrichment efficiency of amplicon outran that of capture in more challenging samples. We found that capture was not as accurate as meta and amplicon in identifying between-sample variations, whereas amplicon method was not as accurate as the other two in investigating within-sample variations, suggesting amplicon sequencing was not suitable for studying virus-host interactions and viral transmission that heavily rely on intra-host dynamics. We illustrated that meta uncovered rich genetic information in the clinical samples besides SARS-CoV-2, providing references for clinical diagnostics and therapeutics. Taken all factors above and cost-effectiveness into consideration, we proposed guidance for how to choose sequencing strategy for SARS-CoV-2 under different situations.

CONCLUSIONS

This is, to the best of our knowledge, the first work systematically investigating inter- and intra-individual variations of SARS-CoV-2 using amplicon- and capture-based whole-genome sequencing, as well as the first comparative study among multiple approaches. Our work offers practical solutions for genome sequencing and analyses of SARS-CoV-2 and other emerging viruses.

Protein Evolution in the Flaviviruses

Proteins are commonly used as molecular targets against pathogens such as viruses and bacteria. However, pathogens can evolve rapidly permitting their populations to increase in protein diversity over time and thus escape to the activity of a molecular therapy. Subsequently, in order to design more durable and robust therapies as well as to understand viral evolution in a host and subsequent transmission, it is central to understand the evolution of pathogen proteins. This understanding can enable the detection of protein regions that can be potential targets for therapies and predict the emergence of molecular resistance against therapies. In this direction, two articles published recently in the Journal of Molecular Evolution investigated the evolution of proteomes of diverse flaviviruses, including Zika virus, Dengue virus and West Nile virus. Here I discuss the importance of considering the evolution of viral proteins, with the use of as realistic as possible models and methods that mimic protein evolution, to improve the design of antiviral therapies.

A synthesis of virus-vector associations reveals important deficiencies in studies on host and vector manipulation by plant viruses

Plant viruses face many challenges in agricultural environments. Although crop fields appear to be abundant resources for these pathogens, it may be difficult for viruses to "escape" from crop environments prior to host senescence or harvesting. One way for viruses to increase the odds of persisting outside of agricultural fields across seasons is by evolving traits that increase transmission opportunities between crops and wild plant communities. There is accumulating evidence that some viruses can achieve this by manipulating crop plant phenotypes in ways that enhance transmission by vectors. Putative manipulations occur through alteration of plant cues (color, size, texture, foliar volatiles, in-leaf metabolites, defenses, and leaf cuticles) that mediate vector orientation, feeding, and dispersal behaviors. Virus effects on host phenotypes are not uniform but appear to exhibit convergence depending on virus traits underlying transmission, particularly the duration of probing and feeding required to acquire and inoculate distinct types of plant viruses. This shared congruence in manipulation strategies and mechanisms across divergent virus lineages suggests that such effects may be adaptive. To discern if this is the case, researchers must consider molecular and environmental constraints on virus evolution, including those imposed by insect vectors from organismal to landscape scales. In this review, we synthesize applied research on vector-borne virus transmission in laboratory and field settings to identify the main factors determining transmission opportunities for plant viruses, and thus, selection pressure to evolve manipulative traits. We then examine these outputs in the context of studies reporting putative instances of plant virus manipulation. Our synthesis reveals important disconnects between virus manipulation studies and actual selection pressures imposed by vectors in real-world contexts.

New insights into the evolutionary features of viral overlapping genes by discriminant analysis

Overlapping genes originate by a mechanism of overprinting, in which nucleotide substitutions in a pre-existing frame induce the expression of a de novo protein from an alternative frame. In this study, I assembled a dataset of 319 viral overlapping genes, which included 82 overlaps whose expression is experimentally known and the respective 237 homologs. Principal component analysis revealed that overlapping genes have a common pattern of nucleotide and amino acid composition. Discriminant analysis separated overlapping from non-overlapping genes with an accuracy of 97%. When applied to overlapping genes with known genealogy, it separated ancestral from de novo frames with an accuracy close to 100%. This high discriminant power was crucial to computationally design variants of de novo viral proteins known to possess selective anticancer toxicity (apoptin) or protection against neurodegeneration (X protein), as well as to detect two new potential overlapping genes in the genome of the new coronavirus SARS-CoV-2.