Evolution in health and medicine Sackler colloquium: The comparative genomics of viral emergence

Natural selection differences detected in key protein domains between non-pathogenic and pathogenic feline coronavirus phenotypes

Feline coronaviruses (FCoVs) commonly cause mild enteric infections in felines worldwide (termed feline enteric coronavirus [FECV]), with around 12 per cent developing into deadly feline infectious peritonitis (FIP; feline infectious peritonitis virus [FIPV]). Genomic differences between FECV and FIPV have been reported, yet the putative genotypic basis of the highly pathogenic phenotype remains unclear. Here, we used state-of-the-art molecular evolutionary genetic statistical techniques to identify and compare differences in natural selection pressure between FECV and FIPV sequences, as well as to identify FIPV- and FECV-specific signals of positive selection. We analyzed full-length FCoV protein coding genes thought to contain mutations associated with FIPV (Spike, ORF3abc, and ORF7ab). We identified two sites exhibiting differences in natural selection pressure between FECV and FIPV: one within the S1/S2 furin cleavage site (FCS) and the other within the fusion domain of Spike. We also found fifteen sites subject to positive selection associated with FIPV within Spike, eleven of which have not previously been suggested as possibly relevant to FIP development. These sites fall within Spike protein subdomains that participate in host cell receptor interaction, immune evasion, tropism shifts, host cellular entry, and viral escape. There were fourteen sites (twelve novel sites) within Spike under positive selection associated with the FECV phenotype, almost exclusively within the S1/S2 FCS and adjacent to C domain, along with a signal of relaxed selection in FIPV relative to FECV, suggesting that furin cleavage functionality may not be needed for FIPV. Positive selection inferred in ORF7b was associated with the FECV phenotype and included twenty-four positively selected sites, while ORF7b had signals of relaxed selection in FIPV. We found evidence of positive selection in ORF3c in FCoV-wide analyses, but no specific association with the FIPV or FECV phenotype. We hypothesize that some combination of mutations in FECV may contribute to FIP development, and that it is unlikely to be one singular 'switch' mutational event. This work expands our understanding of the complexities of FIP development and provides insights into how evolutionary forces may alter pathogenesis in coronavirus genomes.

Combination Therapy for the Treatment of Shingles with an Immunostimulatory Vaccine Virus and Acyclovir

Practically the entire global population is infected by herpesviruses that establish lifelong latency and can be reactivated. Alpha-herpesviruses, herpes simplex viruses 1 and 2 (HSV-1/HSV-2) and varicella zoster virus (VZV), establish latency in sensory neurons and then reactivate to infect epithelial cells in the mucosa or skin, resulting in a vesicular rash. Licensed antivirals inhibit virus replication, but do not affect latency. On reactivation, VZV causes herpes zoster, also known as shingles. The 76-year-old first author of this paper published an autobiography of his own severe herpes zoster ophthalmicus (HZO) infection with orbital edema, which is considered an emergency condition. Acyclovir (ACV) treatment was complemented with an immunostimulatory viral therapy, which resolved most symptoms within a few days. The orally administered live-attenuated infectious bursal disease vaccine virus (IBDV) delivers its double-stranded RNA (dsRNA) cargo to host cells and activates the natural antiviral interferon (IFN) gene defense system from within the host cells. IBDV has already been demonstrated to be safe and effective against five different families of viruses, hepatitis A virus (HAV), hepatitis B and C virus (HBV/HCV), severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), and varicella zoster virus (VZV). Here we propose a short phase I/II trial in elderly shingles patients who will be assigned to receive either ACV monotherapy or ACV combined with R903/78, an attenuated immunostimulatory IBDV strain. The primary endpoints will be safety, but the efficacy of the combination therapy against the ACV monotherapy also will be assessed.

Natural selection differences detected in key protein domains between non-pathogenic and pathogenic Feline Coronavirus phenotypes

Feline Coronaviruses (FCoVs) commonly cause mild enteric infections in felines worldwide (termed Feline Enteric Coronavirus [FECV]), with around 12% developing into deadly Feline Infectious Peritonitis (FIP; Feline Infectious Peritonitis Virus [FIPV]). Genomic differences between FECV and FIPV have been reported, yet the putative genotypic basis of the highly pathogenic phenotype remains unclear. Here, we used state-of-the-art molecular evolutionary genetic statistical techniques to identify and compare differences in natural selection pressure between FECV and FIPV sequences, as well as to identify FIPV and FECV specific signals of positive selection. We analyzed full length FCoV protein coding genes thought to contain mutations associated with FIPV (Spike, ORF3abc, and ORF7ab). We identified two sites exhibiting differences in natural selection pressure between FECV and FIPV: one within the S1/S2 furin cleavage site, and the other within the fusion domain of Spike. We also found 15 sites subject to positive selection associated with FIPV within Spike, 11 of which have not previously been suggested as possibly relevant to FIP development. These sites fall within Spike protein subdomains that participate in host cell receptor interaction, immune evasion, tropism shifts, host cellular entry, and viral escape. There were 14 sites (12 novel) within Spike under positive selection associated with the FECV phenotype, almost exclusively within the S1/S2 furin cleavage site and adjacent C domain, along with a signal of relaxed selection in FIPV relative to FECV, suggesting that furin cleavage functionality may not be needed for FIPV. Positive selection inferred in ORF7b was associated with the FECV phenotype, and included 24 positively selected sites, while ORF7b had signals of relaxed selection in FIPV. We found evidence of positive selection in ORF3c in FCoV wide analyses, but no specific association with the FIPV or FECV phenotype. We hypothesize that some combination of mutations in FECV may contribute to FIP development, and that is unlikely to be one singular "switch" mutational event. This work expands our understanding of the complexities of FIP development and provides insights into how evolutionary forces may alter pathogenesis in coronavirus genomes.

Non-uniform aspects of the SARS-CoV-2 intraspecies evolution reopen question of its origin

Several hypotheses have been presented on the origin of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) from its identification as the agent causing the current coronavirus disease 19 (COVID-19) pandemic. So far, no solid evidence has been found to support any hypothesis on the origin of this virus, and the issue continue to resurface over and over again. Here we have unfolded a pattern of distribution of several mutations in the SARS-CoV-2 proteins in 24 geo-locations across different continents. The results showed an evenly uneven distribution of the unique protein variants, distinct mutations, unique frequency of common conserved residues, and mutational residues across these 24 geo-locations. Furthermore, ample mutations were identified in the evolutionarily conserved invariant regions in the SARS-CoV-2 proteins across almost all geo-locations studied. This pattern of mutations potentially breaches the law of evolutionary conserved functional units of the beta-coronavirus genus. These mutations may lead to several novel SARS-CoV-2 variants with a high degree of transmissibility and virulence. A thorough investigation on the origin and characteristics of SARS-CoV-2 needs to be conducted in the interest of science and for the preparation of meeting the challenges of potential future pandemics.

Geographic Range Overlap Rather than Phylogenetic Distance Explains Rabies Virus Transmission among Closely Related Bat Species

The cross-species transmission (CST) of pathogens can have dramatic consequences, as highlighted by recent disease emergence events affecting human, animal and plant health. Understanding the ecological and evolutionary factors that increase the likelihood of disease agents infecting and establishing in a novel host is therefore an important research area. Previous work across different pathogens, including rabies virus (RABV), found that increased evolutionary distance between hosts reduces the frequency of cross-species transmission and of permanent host shifts. However, whether this effect of host relatedness still holds for transmission among recently diverged hosts is not well understood. We aimed to ask if high host relatedness can still increase the probability of a host shift between more recently diverged hosts, and the importance of this effect relative to ecological predictors. We first addressed this question by quantifying the CST frequency of RABV between North American bat species within the genus Myotis, using a multi-decade data set containing 128 nucleoprotein (N) RABV sequences from ten host species. We compared RABV CST frequency within Myotis to the rates of CST between nine genera of North American bat species. We then examined whether host relatedness or host range overlap better explains the frequency of CST seen between Myotis species. We found that at the within genus scale, host range overlap, rather than host relatedness best explains the frequency of CST events. Moreover, we found evidence of CST occurring among a higher proportion of species, and CST more frequently resulting in sustained transmission in the novel host in the Myotis dataset compared to the multi-genus dataset. Our results suggest that among recently diverged species, the ability to infect a novel host is no longer restricted by physiological barriers but instead is limited by physical contact. Our results improve predictions of where future CST events for RABV might occur and clarify the relationship between host divergence and pathogen emergence.

Developing an empirical model for spillover and emergence: Orsay virus host range in Caenorhabditis

A lack of tractable experimental systems in which to test hypotheses about the ecological and evolutionary drivers of disease spillover and emergence has limited our understanding of these processes. Here we introduce a promising system: Caenorhabditis hosts and Orsay virus, a positive-sense single-stranded RNA virus that naturally infects C. elegans. We assayed species across the Caenorhabditis tree and found Orsay virus susceptibility in 21 of 84 wild strains belonging to 14 of 44 species. Confirming patterns documented in other systems, we detected effects of host phylogeny on susceptibility. We then tested whether susceptible strains were capable of transmitting Orsay virus by transplanting exposed hosts and determining whether they transmitted infection to conspecifics during serial passage. We found no evidence of transmission in 10 strains (virus undetectable after passaging in all replicates), evidence of low-level transmission in 5 strains (virus lost between passage 1 and 5 in at least one replicate) and evidence of sustained transmission in 6 strains (including all three experimental C. elegans strains) in at least one replicate. Transmission was strongly associated with viral amplification in exposed populations. Variation in Orsay virus susceptibility and transmission among Caenorhabditis strains suggests that the system could be powerful for studying spillover and emergence.

Discovery of two novel Torque Teno viruses in Callithrix penicillata provides insights on Anelloviridae diversification dynamics

The development of high-throughput sequencing (HTS) technologies and metagenomics protocols deeply impacted the discovery of viral diversity. Moreover, the characterization of novel viruses in the Neotropical primates (NP) is central for the comprehension of viral evolution dynamics in those hosts, due to their evolutionary proximity to Old World primates, including humans. In the present work, novel anelloviruses were detected and characterized through HTS protocols in the NP Callithrix penicillata, the common black-tufted marmoset. De novo assembly of generated sequences was carried out, and a total of 15 contigs were identified with complete Anelloviridae ORF1 gene, two of them including a flanking GC-rich region, confirming the presence of two whole novel genomes of ~3 kb. The identified viruses were monophyletic within the Epsilontorquevirus genus, a lineage harboring previously reported anelloviruses infecting hosts from the Cebidae family. The genetic divergence found in the new viruses characterized two novel species, named Epsilontorquevirus callithrichensis I and II. The phylogenetic pattern inferred for the Epsilontorquevirus genus was consistent with the topology of their host species tree, echoing a virus-host diversification model observed in other viral groups. This study expands the host span of Anelloviridae and provides insights into their diversification dynamics, highlighting the importance of sampling animal viral genomes to obtain a clearer depiction of their long-term evolutionary processes.

An Orally Administered Nonpathogenic Attenuated Vaccine Virus Can Be Used to Control SARS-CoV-2 Infection: A Complementary Plan B to COVID-19 Vaccination

Background

Coronavirus disease 2019 (COVID-19) vaccination has substantially altered the course of the pandemic, saving tens of millions of lives globally. The problem is that despite such spectacular results, vaccination alone will not be able to control the COVID-19 pandemic because of the rapid evolution of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) even in vaccinated human populations. Therefore, the development of a post-infection, broad-based, orally administered antiviral therapy that would complement vaccination efforts is urgently needed.

Methodology

The so-called viral superinfection therapy (SIT) administers a nonpathogenic attenuated double-stranded RNA (dsRNA) vaccine virus drug candidate, the infectious bursal disease virus serotype R903/78 (IBDV-R903/78) that activates the interferon (IFN) genes, which are the natural, antiviral defense system of host cells.

Results

Here we present two cases of properly vaccinated (with BNT162b2-Pfizer) and booster-dosed COVID-19 patients with vaccine breakthrough infections whose disease duration was shortened to a few days by oral SIT.

Conclusions

SIT has already been demonstrated to be safe and effective against five different families of viruses, hepatitis A virus, hepatitis B virus, hepatitis C virus, SARS-CoV-2, and herpes zoster virus. The R903/78 drug candidate is simple to manufacture and easy to administer in an outpatient setting. The expected cost of SIT will be affordable even in resource-limited countries.

Evolution and host adaptability of plant RNA viruses: Research insights on compositional biases

During recent decades, many new emerging or re-emerging RNA viruses have been found in plants through the development of deep-sequencing technology and big data analysis. These findings largely changed our understanding of the origin, evolution and host range of plant RNA viruses. There is evidence that their genetic composition originates from viruses, and host populations play a key role in the evolution and host adaptability of plant RNA viruses. In this mini-review, we describe the state of our understanding of the evolution of plant RNA viruses in view of compositional biases and explore how they adapt to the host. It appears that adenine rich (A-rich) coding sequences, low CpG and UpA dinucleotide frequencies and lower codon usage patterns were found in the vast majority of plant RNA viruses. The codon usage pattern of plant RNA viruses was influenced by both natural selection and mutation pressure, and natural selection mostly from hosts was the dominant factor. The codon adaptation analyses support that plant RNA viruses probably evolved a dynamic balance between codon adaptation and deoptimization to maintain efficient replication cycles in multiple hosts with various codon usage patterns. In the future, additional combinations of computational and experimental analyses of the nucleotide composition and codon usage of plant RNA viruses should be addressed.

Circulation and Evolution of SARS-CoV-2 in India: Let the Data Speak

The COVID-19 pandemic is a global challenge that impacted 200+ countries. India ranks in the second and third positions in terms of number of reported cases and deaths. Being a populous country with densely packed cities, SARS-CoV-2 spread exponentially. India sequenced ≈0.14% isolates from confirmed cases for pandemic surveillance and contributed ≈1.58% of complete genomes sequenced globally. This study was designed to map the circulating lineage diversity and to understand the evolution of SARS-CoV-2 in India using comparative genomics and population genetics approaches. Despite varied sequencing coverage across Indian States and Union Territories, isolates belonging to variants of concern (VoC) and variants of interest (VoI) circulated, persisted, and diversified during the first seventeen months of the pandemic. Delta and Kappa lineages emerged in India and spread globally. The phylogenetic tree shows lineage-wise monophyletic clusters of VoCs/VoIs and diversified tree topologies for non-VoC/VoI lineages designated as 'Others' in this study. Evolutionary dynamics analyses substantiate a lack of spatio-temporal clustering, which is indicative of multiple global and local introductions. Sites under positive selection and significant variations in spike protein corroborate with the constellation of mutations to be monitored for VoC/VoI as well as substitutions that are characteristic of functions with implications in virus-host interactions, differential glycosylation, immune evasion, and escape from neutralization.

Microglial Implications in SARS-CoV-2 Infection and COVID-19: Lessons From Viral RNA Neurotropism and Possible Relevance to Parkinson's Disease

Since December 2019, humankind has been experiencing a ravaging severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) outbreak, the second coronavirus pandemic in a decade after the Middle East respiratory syndrome coronavirus (MERS-CoV) disease in 2012. Infection with SARS-CoV-2 results in Coronavirus disease 2019 (COVID-19), which is responsible for over 3.1 million deaths worldwide. With the emergence of a second and a third wave of infection across the globe, and the rising record of multiple reinfections and relapses, SARS-CoV-2 infection shows no sign of abating. In addition, it is now evident that SARS-CoV-2 infection presents with neurological symptoms that include early hyposmia, ischemic stroke, meningitis, delirium and falls, even after viral clearance. This may suggest chronic or permanent changes to the neurons, glial cells, and/or brain vasculature in response to SARS-CoV-2 infection or COVID-19. Within the central nervous system (CNS), microglia act as the central housekeepers against altered homeostatic states, including during viral neurotropic infections. In this review, we highlight microglial responses to viral neuroinfections, especially those with a similar genetic composition and route of entry as SARS-CoV-2. As the primary sensor of viral infection in the CNS, we describe the pathogenic and neuroinvasive mechanisms of RNA viruses and SARS-CoV-2 vis-à-vis the microglial means of viral recognition. Responses of microglia which may culminate in viral clearance or immunopathology are also covered. Lastly, we further discuss the implication of SARS-CoV-2 CNS invasion on microglial plasticity and associated long-term neurodegeneration. As such, this review provides insight into some of the mechanisms by which microglia could contribute to the pathophysiology of post-COVID-19 neurological sequelae and disorders, including Parkinson's disease, which could be pervasive in the coming years given the growing numbers of infected and re-infected individuals globally.

Interplay between H1N1 influenza a virus infection, extracellular and intracellular respiratory tract pH, and host responses in a mouse model

During influenza A virus (IAV) entry, the hemagglutinin (HA) protein is triggered by endosomal low pH to undergo irreversible structural changes that mediate membrane fusion. HA proteins from different isolates vary in the pH at which they become activated in endosomes or become irreversible inactivated if exposed to extracellular acid. Little is known about extracellular pH in the upper respiratory tracts of mammals, how pH may shift during IAV infection, and its impact on replication of viruses that vary in HA activation pH. Here, we inoculated DBA/2J mice intranasally with A/TN/1-560/2009 (H1N1) (activation pH 5.5) or a mutant containing the destabilizing mutation HA1-Y17H (pH 6.0). We measured the kinetics of extracellular pH during infection using an optical pH-sensitive microsensor probe placed in the naris, nasal sinus, soft palate, and trachea. We also measured intracellular pH of single-cell suspensions of live, primary lung epithelial cells with various wavelength pH-sensitive dyes localized to cell membranes, cytosol, endosomes, secretory vesicles, microtubules, and lysosomes. Infection with either virus decreased extracellular pH and increased intracellular pH. Peak host immune responses were observed at 2 days post infection (DPI) and peak pH changes at 5 DPI. Extracellular and intracellular pH returned to baseline by 7 DPI in mice infected with HA1-Y17H and was restored later in wildtype-infected. Overall, IAV infection altered respiratory tract pH, which in turn modulated replication efficiency. This suggests a virus-host pH feedback loop that may select for IAV strains containing HA proteins of optimal pH stability, which may be approximately pH 5.5 in mice but may differ in other species.

SARS-CoV-2 exhibits intra-host genomic plasticity and low-frequency polymorphic quasispecies

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

The phylogenetic range of bacterial and viral pathogens of vertebrates

Many major human pathogens are multihost pathogens, able to infect other vertebrate species. Describing the general patterns of host-pathogen associations across pathogen taxa is therefore important to understand risk factors for human disease emergence. However, there is a lack of comprehensive curated databases for this purpose, with most previous efforts focusing on viruses. Here, we report the largest manually compiled host-pathogen association database, covering 2,595 bacteria and viruses infecting 2,656 vertebrate hosts. We also build a tree for host species using nine mitochondrial genes, giving a quantitative measure of the phylogenetic similarity of hosts. We find that the majority of bacteria and viruses are specialists infecting only a single host species, with bacteria having a significantly higher proportion of specialists compared to viruses. Conversely, multihost viruses have a more restricted host range than multihost bacteria. We perform multiple analyses of factors associated with pathogen richness per host species and the pathogen traits associated with greater host range and zoonotic potential. We show that factors previously identified as important for zoonotic potential in viruses-such as phylogenetic range, research effort, and being vector-borne-are also predictive in bacteria. We find that the fraction of pathogens shared between two hosts decreases with the phylogenetic distance between them. Our results suggest that host phylogenetic similarity is the primary factor for host-switching in pathogens.

Porcine circovirus 2 (PCV-2) genotype update and proposal of a new genotyping methodology

Porcine circovirus 2 (PCV-2) is one of the most widespread viral infections of swine, causing a remarkable economic impact because of direct losses and indirect costs for its control. As other ssDNA viruses, PCV-2 is characterized by a high evolutionary rate, leading to the emergence of a plethora of variants with different biological and epidemiological features. Over time, several attempts have been made to organize PCV-2 genetic heterogeneity in recognized genotypes. This categorization has clearly simplified the epidemiological investigations, allowing to identify different spatial and temporal patterns among genotypes. Additionally, variable virulence and vaccine effectiveness have also been hypothesized. However, the rapid increase in sequencing activity, coupled with the per se high viral variability, has challenged the previously established nomenclature, leading to the definition of several study-specific genotypes and hindering the capability of performing comparable epidemiological studies.

Based on these premises, an updated classification scheme is herein reported. Recognizing the impossibility of defining a clear inter-cluster p-distance cut-off, the present study proposes a phylogeny-grounded genotype definition based on three criteria: maximum intra-genotype p-distance of 13% (calculated on the ORF2 gene), bootstrap support at the corresponding internal node higher than 70% and at least 15 available sequences. This scheme allowed defining 8 genotypes (PCV-2a to PCV-2h), which six of those had been previously proposed. To minimize the inconvenience of implementing a new classification, the most common names already adopted have been maintained when possible. The analysis of sequence-associated metadata highlighted a highly unbalanced sequencing activity in terms of geographical, host and temporal distribution. The PCV-2 molecular epidemiology scenario appears therefore characterized by a severe bias that could lead to spurious associations between genetic and epidemiological/biological viral features. While the suggested classification can establish a “common language” for future studies, further efforts should be paid to achieve a more homogeneous and informative representation of the PCV-2 global scenario.

Transmissibility of emerging viral zoonoses

Effective public health research and preparedness requires an accurate understanding of which virus species possess or are at risk of developing human transmissibility. Unfortunately, our ability to identify these viruses is limited by gaps in disease surveillance and an incomplete understanding of the process of viral adaptation. By fitting boosted regression trees to data on 224 human viruses and their associated traits, we developed a model that predicts the human transmission ability of zoonotic viruses with over 84% accuracy. This model identifies several viruses that may have an undocumented capacity for transmission between humans. Viral traits that predicted human transmissibility included infection of nonhuman primates, the absence of a lipid envelope, and detection in the human nervous system and respiratory tract. This predictive model can be used to prioritize high-risk viruses for future research and surveillance, and could inform an integrated early warning system for emerging infectious diseases.

Emergency Services of Viral RNAs: Repair and Remodeling

Reproduction of RNA viruses is typically error-prone due to the infidelity of their replicative machinery and the usual lack of proofreading mechanisms. The error rates may be close to those that kill the virus. Consequently, populations of RNA viruses are represented by heterogeneous sets of genomes with various levels of fitness. This is especially consequential when viruses encounter various bottlenecks and new infections are initiated by a single or few deviating genomes. Nevertheless, RNA viruses are able to maintain their identity by conservation of major functional elements. This conservatism stems from genetic robustness or mutational tolerance, which is largely due to the functional degeneracy of many protein and RNA elements as well as to negative selection. Another relevant mechanism is the capacity to restore fitness after genetic damages, also based on replicative infidelity. Conversely, error-prone replication is a major tool that ensures viral evolvability. The potential for changes in debilitated genomes is much higher in small populations, because in the absence of stronger competitors low-fit genomes have a choice of various trajectories to wander along fitness landscapes. Thus, low-fit populations are inherently unstable, and it may be said that to run ahead it is useful to stumble. In this report, focusing on picornaviruses and also considering data from other RNA viruses, we review the biological relevance and mechanisms of various alterations of viral RNA genomes as well as pathways and mechanisms of rehabilitation after loss of fitness. The relationships among mutational robustness, resilience, and evolvability of viral RNA genomes are discussed.

Live Attenuated Influenza Vaccines engineered to express the nucleoprotein of a recent isolate stimulate human influenza CD8+ T cells more relevant to current infections

Live attenuated influenza vaccines (LAIV) induce CD8⁺ T lymphocyte responses that play an important role in killing virus-infected cells. Despite the relative conservation of internal influenza A proteins, the epitopes recognized by T cells can undergo drift under immune pressure. The internal proteins of Russian LAIVs are derived from the master donor virus A/Leningrad/134/17/57 (Len/17) isolated 60 years ago and as such, some CD8⁺ T cell epitopes may vary between the vaccine and circulating wild-type strains. To partially overcome this issue, the nucleoprotein (NP) gene of wild-type virus can be incorporated into LAIV reassortant virus, along with the HA and NA genes. The present study compares the human CD8+ T cell memory responses to H3N2 LAIVs with the Len/17 or the wild-type NP using an in vitro model.

Public health concerns of highly pathogenic avian influenza H5N1 endemicity in Africa

Highly pathogenic avian influenza virus (HPAIV) H5N1 was first officially reported in Africa in 2006; thereafter this virus has spread rapidly from Nigeria to 11 other African countries. This study was aimed at utilizing data from confirmed laboratory reports to carry out a qualitative evaluation of the factors responsible for HPAI H5N1 persistence in Africa and the public health implications; and to suggest appropriate control measures. Relevant publications were sought from data banks and repositories of FAO, OIE, WHO, and Google scholars. Substantiated data on HPAI H5N1 outbreaks in poultry in Africa and in humans across the world were mined. HPAI H5N1 affects poultry and human populations, with Egypt having highest human cases (346) globally. Nigeria had a reinfection from 2014 to 2015, with outbreaks in Côte d'Ivoire, Ghana, Niger, Nigeria, and Burkina Faso throughout 2016 unabated. The persistence of this virus in Africa is attributed to the survivability of HPAIV, ability to evolve other subtypes through genetic reassortment, poor biosecurity compliance at the live bird markets and poultry farms, husbandry methods and multispecies livestock farming, poultry vaccinations, and continuous shedding of HPAIV, transboundary transmission of HPAIV through poultry trades; and transcontinental migratory birds. There is, therefore, the need for African nations to realistically reassess their status, through regular surveillance and be transparent with HPAI H5N1 outbreak data. Also, it is important to have an understanding of HPAIV migration dynamics which will be helpful in epidemiological modeling, disease prevention, control and eradication measures.

Viral recombination blurs taxonomic lines: examination of single-stranded DNA viruses in a wastewater treatment plant

Understanding the structure and dynamics of microbial communities, especially those of economic concern, is of paramount importance to maintaining healthy and efficient microbial communities at agricultural sites and large industrial cultures, including bioprocessors. Wastewater treatment plants are large bioprocessors which receive water from multiple sources, becoming reservoirs for the collection of many viral families that infect a broad range of hosts. To examine this complex collection of viruses, full-length genomes of circular ssDNA viruses were isolated from a wastewater treatment facility using a combination of sucrose-gradient size selection and rolling-circle amplification and sequenced on an Illumina MiSeq. Single-stranded DNA viruses are among the least understood groups of microbial pathogens due to genomic biases and culturing difficulties, particularly compared to the larger, more often studied dsDNA viruses. However, the group contains several notable well-studied examples, including agricultural pathogens which infect both livestock and crops (Circoviridae and Geminiviridae), and model organisms for genetics and evolution studies (Microviridae). Examination of the collected viral DNA provided evidence for 83 unique genotypic groupings, which were genetically dissimilar to known viral types and exhibited broad diversity within the community. Furthermore, although these genomes express similarities to known viral families, such as Circoviridae, Geminiviridae, and Microviridae, many are so divergent that they may represent new taxonomic groups. This study demonstrated the efficacy of the protocol for separating bacteria and large viruses from the sought after ssDNA viruses and the ability to use this protocol to obtain an in-depth analysis of the diversity within this group.

Positive Selection Drives Rapid Evolution of the meq Oncogene of Marek's Disease Virus

Marek’s disease (MD), caused by Marek’s disease virus (MDV), a poultry-borne alphaherpesvirus, is a devastating disease of poultry causing an estimated annual loss of one billion dollars to poultry producers, worldwide. Despite decades of control through vaccination, MDV field strains continue to emerge having increased virulence. The evolutionary mechanism driving the emergence of this continuum of strains to increased MDV virulence, however, remains largely enigmatic. Increase in MDV virulence has been associated with specific amino acid changes within the C-terminus domain of Mareks’s EcoRI-Q (meq)-encoded oncoprotein. In this study, we sought to determine whether the meq gene has evolved adaptively and whether past vaccination efforts have had any significant effect on the reduction or increase of MDV diversity over time. Our analysis suggests that meq is estimated to be evolving at a much faster rate than most dsDNA viruses, and is comparable with the evolutionary rate of RNA viruses. Interestingly, most of the polymorphisms in meq gene appear to have evolved under positive selection and the time of divergence at the meq locus coincides with the period during which the poultry industry had undergone transitions in management practices including the introduction and widespread use of live attenuated vaccines. Our study has revealed that the decades-long use of vaccines did not reduce MDV diversity, but rather had a stimulating effect on the emergence of field strains with increased genetic diversity until the early 2000s. During the years 2004–2005, there was an abrupt decline in the genetic diversity of field isolates followed by a recovery from this bottleneck in the year 2010. Collectively, these data suggest that vaccination seems to not have had any effect on MDV eradication, but rather had a stimulating effect on MDV emergence through adaptation.

The transcriptomics of an experimentally evolved plant-virus interaction

Models of plant-virus interaction assume that the ability of a virus to infect a host genotype depends on the matching between virulence and resistance genes. Recently, we evolved tobacco etch potyvirus (TEV) lineages on different ecotypes of Arabidopsis thaliana, and found that some ecotypes selected for specialist viruses whereas others selected for generalists. Here we sought to evaluate the transcriptomic basis of such relationships. We have characterized the transcriptomic responses of five ecotypes infected with the ancestral and evolved viruses. Genes and functional categories differentially expressed by plants infected with local TEV isolates were identified, showing heterogeneous responses among ecotypes, although significant parallelism existed among lineages evolved in the same ecotype. Although genes involved in immune responses were altered upon infection, other functional groups were also pervasively over-represented, suggesting that plant resistance genes were not the only drivers of viral adaptation. Finally, the transcriptomic consequences of infection with the generalist and specialist lineages were compared. Whilst the generalist induced very similar perturbations in the transcriptomes of the different ecotypes, the perturbations induced by the specialist were divergent. Plant defense mechanisms were activated when the infecting virus was specialist but they were down-regulated when infecting with generalist.

Virological factors that increase the transmissibility of emerging human viruses

The early detection of pathogens with epidemic potential is of major importance to public health. Most emerging infections result in dead-end "spillover" events in which a pathogen is transmitted from an animal reservoir to a human but is unable to achieve the sustained human-to-human transmission necessary for a full-blown epidemic. It is therefore critical to determine why only some virus infections are efficiently transmitted among humans whereas others are not. We sought to determine which biological features best characterized those viruses that have achieved sustained human transmission. Accordingly, we compiled a database of 203 RNA and DNA human viruses and used an information theoretic approach to assess which of a set of key biological variables were the best predictors of human-to-human transmission. The variables analyzed were as follows: taxonomic classification; genome length, type, and segmentation; the presence or absence of an outer envelope; recombination frequency; duration of infection; host mortality; and whether or not a virus exhibits vector-borne transmission. This comparative analysis revealed multiple strong associations. In particular, we determined that viruses with low host mortality, that establish long-term chronic infections, and that are nonsegmented, nonenveloped, and, most importantly, not transmitted by vectors were more likely to be transmissible among humans. In contrast, variables including genome length, genome type, and recombination frequency had little predictive power. In sum, we have identified multiple biological features that seemingly determine the likelihood of interhuman viral transmissibility, in turn enabling general predictions of whether viruses of a particular type will successfully emerge in human populations.

Virus classification in 60-dimensional protein space

Due to vast sequence divergence among different viral groups, sequence alignment is not directly applicable to genome-wide comparative analysis of viruses. More and more attention has been paid to alignment-free methods for whole genome comparison and phylogenetic tree reconstruction. Among alignment-free methods, the recently proposed "Natural Vector (NV) representation" has successfully been used to study the phylogeny of multi-segmented viruses based on a 12-dimensional genome space derived from the nucleotide sequence structure. But the preference of proteomes over genomes for the determination of viral phylogeny was not deeply investigated. As the translated products of genes, proteins directly form the shape of viral structure and are vital for all metabolic pathways. In this study, using the NV representation of a protein sequence along with the Hausdorff distance suitable to compare point sets, we construct a 60-dimensional protein space to analyze the evolutionary relationships of 4021 viruses by whole-proteomes in the current NCBI Reference Sequence Database (RefSeq). We also take advantage of the previously developed natural graphical representation to recover viral phylogeny. Our results demonstrate that the proposed method is efficient and accurate for classifying viruses. The accuracy rates of our predictions such as for Baltimore II viruses are as high as 95.9% for family labels, 95.7% for subfamily labels and 96.5% for genus labels. Finally, we discover that proteomes lead to better viral classification when reliable protein sequences are abundant. In other cases, the accuracy rates using proteomes are still comparable to that of genomes.

Utilization of multiple "omics" studies in microbial pathogeny for microbiology insights

In the present day, bioinformatics becomes the modern science with several advantages. Several new "omics" sciences have been introduced for a few years and those sciences can be applied in biomedical work. Here, the author will summarize and discuss on important applications of omics studies in microbiology focusing on microbial pathogeny. It can be seen that genomics and proteinomics can be well used in this area of biomedical studies.

Full genomic analysis of new variant rabbit hemorrhagic disease virus revealed multiple recombination events

Rabbit hemorrhagic disease virus (RHDV), a Lagovirus of the family Caliciviridae, causes rabbit hemorrhagic disease (RHD) in the European rabbit (Oryctolagus cuniculus). The disease was first documented in 1984 in China and rapidly spread worldwide. In 2010, a new RHDV variant emerged, tentatively classified as 'RHDVb'. RHDVb is characterized by affecting vaccinated rabbits and those <2 months old, and is genetically distinct (~20 %) from older strains. To determine the evolution of RHDV, including the new variant, we generated 28 full-genome sequences from samples collected between 1994 and 2014. Phylogenetic analysis of the gene encoding the major capsid protein, VP60, indicated that all viruses sampled from 2012 to 2014 were RHDVb. Multiple recombination events were detected in the more recent RHDVb genomes, with a single major breakpoint located in the 5' region of VP60. This breakpoint divides the genome into two regions: one that encodes the non-structural proteins and another that encodes the major and minor structural proteins, VP60 and VP10, respectively. Additional phylogenetic analysis of each region revealed two types of recombinants with distinct genomic backgrounds. Recombinants always include the structural proteins of RHDVb, with non-structural proteins from non-pathogenic lagoviruses or from pathogenic genogroup 1 strains. Our results show that in contrast to the evolutionary history of older RHDV strains, recombination plays an important role in generating diversity in the newly emerged RHDVb.

Global comparison of multiple-segmented viruses in 12-dimensional genome space

We have recently developed a computational approach in a vector space for genome-based virus classification. This approach, called the "Natural Vector (NV) representation", which is an alignment-free method, allows us to classify single-segmented viruses with high speed and accuracy. For multiple-segmented viruses, typically phylogenetic trees of each segment are reconstructed for discovering viral phylogeny. Consensus tree methods may be used to combine the phylogenetic trees based on different segments. However, consensus tree methods were not developed for instances where the viruses have different numbers of segments or where their segments do not match well. We propose a novel approach for comparing multiple-segmented viruses globally, even in cases where viruses contain different numbers of segments. Using our method, each virus is represented by a set of vectors in R(12). The Hausdorff distance is then used to compare different sets of vectors. Phylogenetic trees can be reconstructed based on this distance. The proposed method is used for predicting classification labels of viruses with n-segments (n ⩾ 1). The correctness rates of our predictions based on cross-validation are as high as 96.5%, 95.4%, 99.7%, and 95.6% for Baltimore class, family, subfamily, and genus, respectively, which are comparable to the rates for single-segmented viruses only. Our method is not affected by the number or order of segments. We also demonstrate that the natural graphical representation based on the Hausdorff distance is more reasonable than the consensus tree for a recent public health threat, the influenza A (H7N9) viruses.

Increased virulence of rabbit haemorrhagic disease virus associated with genetic resistance in wild Australian rabbits (Oryctolagus cuniculus)

The release of myxoma virus (MYXV) and Rabbit Haemorrhagic Disease Virus (RHDV) in Australia with the aim of controlling overabundant rabbits has provided a unique opportunity to study the initial spread and establishment of emerging pathogens, as well as their co-evolution with their mammalian hosts. In contrast to MYXV, which attenuated shortly after its introduction, rapid attenuation of RHDV has not been observed. By studying the change in virulence of recent field isolates at a single field site we show, for the first time, that RHDV virulence has increased through time, likely because of selection to overcome developing genetic resistance in Australian wild rabbits. High virulence also appears to be favoured as rabbit carcasses, rather than diseased animals, are the likely source of mechanical insect transmission. These findings not only help elucidate the co-evolutionary interaction between rabbits and RHDV, but reveal some of the key factors shaping virulence evolution.

Arboviruses: variations on an ancient theme

ABSTRACT 

Arboviruses utilize different strategies to complete their transmission cycle between vertebrate and invertebrate hosts. Most possess an RNA genome coupled with an RNA polymerase lacking proofreading activity and generate large populations of genetically distinct variants, permitting rapid adaptation to environmental changes. With mutation rates of between 10- 6 and 10-4 substitutions per nucleotide, arboviral genomes rapidly acquire mutations that can lead to viral emergence. Arboviruses can be described in seven families, four of which have medical importance: Togaviridae, Flaviviridae, Bunyaviridae and Reoviridae. The Togaviridae and Flaviviridae both have ssRNA genomes, while the Bunyaviridae and Reoviridae possess segmented RNA genomes. Recent epidemics caused by these arboviruses have been associated with specific mutations leading to enhanced host ranges, vector shifts and virulence.

Viral genome phylogeny based on Lempel-Ziv complexity and Hausdorff distance

In this paper, we develop a novel method to study the viral genome phylogeny. We apply Lempel-Ziv complexity to define the distance between two nucleic acid sequences. Then, based on this distance we use the Hausdorff distance (HD) and a modified Hausdorff distance (MHD) to make the phylogenetic analysis for multi-segmented viral genomes. The results show the MHD can provide more accurate phylogenetic relationship. Our method can have global comparison of all multi-segmented genomes simultaneously, that is, we treat the multi-segmented viral genome as an entirety to make the comparative analysis. Our method is not affected by the number or order of segments, and each segment can make contribution for the phylogeny of whole genomes. We have analyzed several groups of real multi-segmented genomes from different viral families. The results show that our method will provide a new powerful tool for studying the classification of viral genomes and their phylogenetic relationships.

Role of genomic and proteomic tools in the study of host-virus interactions and virus evolution

Viruses have short replication cycles and produce genomic variants within a host, a process that seems to adapt to their specific host and also enable them to infect new hosts. The recent emergence of viral genomic variants from the circulating pool within the host population and re-emergence of the old ones are posing serious threat to agriculture, animal husbandry and humanity as a whole. This review assesses the potential role of genomic and proteomic tools that can monitor not only the course of infection and pathogenesis, but also predict the pandemic or zoonotic epidemic potential of a virus in a previously exposed or immunologically naive biological population.

Multiple host transfers, but only one successful lineage in a continent-spanning emergent pathogen

Emergence of a new disease in a novel host is thought to be a rare outcome following frequent pathogen transfers between host species. However, few opportunities exist to examine whether disease emergence stems from a single successful pathogen transfer, and whether this successful lineage represents only one of several pathogen transfers between hosts. We examined the successful host transfer and subsequent evolution of the bacterial pathogen Mycoplasma gallisepticum, an emergent pathogen of house finches (Haemorhous (formerly Carpodacus) mexicanus). Our principal goals were to assess whether host transfer has been a repeated event between the original poultry hosts and house finches, whether only a single host transfer was ultimately responsible for the emergence of M. gallisepticum in these finches, and whether the spread of the pathogen from east to west across North America has resulted in spatial structuring in the pathogen. Using a phylogeny of M. gallisepticum based on 107 isolates from domestic poultry, house finches and other songbirds, we infer that the bacterium has repeatedly jumped between these two groups of hosts but with only a single lineage of M. gallisepticum persisting and evolving in house finches; bacterial evolution has produced monophyletic eastern and western North American subclades.

Real time classification of viruses in 12 dimensions

The International Committee on Taxonomy of Viruses authorizes and organizes the taxonomic classification of viruses. Thus far, the detailed classifications for all viruses are neither complete nor free from dispute. For example, the current missing label rates in GenBank are 12.1% for family label and 30.0% for genus label. Using the proposed Natural Vector representation, all 2,044 single-segment referenced viral genomes in GenBank can be embedded in [Formula: see text]. Unlike other approaches, this allows us to determine phylogenetic relations for all viruses at any level (e.g., Baltimore class, family, subfamily, genus, and species) in real time. Additionally, the proposed graphical representation for virus phylogeny provides a visualization of the distribution of viruses in [Formula: see text]. Unlike the commonly used tree visualization methods which suffer from uniqueness and existence problems, our representation always exists and is unique. This approach is successfully used to predict and correct viral classification information, as well as to identify viral origins; e.g. a recent public health threat, the West Nile virus, is closer to the Japanese encephalitis antigenic complex based on our visualization. Based on cross-validation results, the accuracy rates of our predictions are as high as 98.2% for Baltimore class labels, 96.6% for family labels, 99.7% for subfamily labels and 97.2% for genus labels.

Molecular epidemiology of circulating highly pathogenic avian influenza (H5N1) virus in chickens, in Bangladesh, 2007-2010

Bangladesh has been severely hit by highly pathogenic avian influenza H5N1 (HPAI-H5N1). However, little is known about the genetic diversity and the evolution of the circulating viruses in Bangladesh. In the present study, we analyzed the hemagglutinin gene of 30 Bangladeshi chicken isolates from 2007 through 2010. We analyzed the polybasic amino acid sequence motif of the cleavage site and amino acid substitution pattern. Phylogenetic history was reconstructed using neighbor-joining and Bayesian time-scaled methods. In addition, we used Mantel correlation tests to analyze the relation between genetic relatedness and spatial and temporal distances. Neighbor-joining phylogeography revealed that virus circulating in Bangladesh from 2007 through 2010 belonged to clade 2.2. The results suggest that clade 2.2 viruses are firmly entrenched and have probably become endemic in Bangladesh. We detected several amino acid substitutions, but they are not indicative of adaptation toward human infection. The Mantel correlation test confirmed significant correlation between genetic distances and temporal distances between the viruses. The Bayesian tree shows that isolates from waves 3 and 4 derived from a subgroup of isolates from the previous waves grouping with a high posterior probability (pp=1.0). This indicates the possibility of formation of local subclades. One surprising finding of spatio-temporal analysis was that genetically identical virus caused independent outbreaks over a distance of more than 200 km and within 14 days of each other. This might indicate long distance dispersal through vectors such as migratory birds and vehicles, and challenges the effectiveness of movement restriction around 10 km radius of an outbreak. The study indicates possible endemicity of the clade 2.2 HPAI-H5N1 virus in Bangladesh. Furthermore, the formation of a subclade capable of transmission to humans cannot be ruled out. The findings of this study might provide valuable information for future surveillance, prevention and control programme.

Implications of altered replication fidelity on the evolution and pathogenesis of coronaviruses

RNA virus evolution results from viral replication fidelity and mutational robustness in combination with selection. Recent studies have confirmed the impact of increased fidelity on RNA virus replication and pathogenesis; however, the impact of decreased fidelity is less defined. Coronaviruses have the largest RNA genomes, and encode an exoribonuclease activity that is required for high-fidelity replication. Genetically stable exoribonuclease mutants will allow direct testing of viral mutational tolerance to RNA mutagens and other selective pressures. Recent studies support the hypothesis that coronavirus replication fidelity may result from a multi-protein complex, suggesting multiple pathways to disrupt or alter virus fidelity and diversity, and attenuate pathogenesis.

Evolutionary history and phylodynamics of influenza A and B neuraminidase (NA) genes inferred from large-scale sequence analyses

Background

Influenza neuraminidase (NA) is an important surface glycoprotein and plays a vital role in viral replication and drug development. The NA is found in influenza A and B viruses, with nine subtypes classified in influenza A. The complete knowledge of influenza NA evolutionary history and phylodynamics, although critical for the prevention and control of influenza epidemics and pandemics, remains lacking.

Methodology/Principal findings

Evolutionary and phylogenetic analyses of influenza NA sequences using Maximum Likelihood and Bayesian MCMC methods demonstrated that the divergence of influenza viruses into types A and B occurred earlier than the divergence of influenza A NA subtypes. Twenty-three lineages were identified within influenza A, two lineages were classified within influenza B, and most lineages were specific to host, subtype or geographical location. Interestingly, evolutionary rates vary not only among lineages but also among branches within lineages. The estimated tMRCAs of influenza lineages suggest that the viruses of different lineages emerge several months or even years before their initial detection. The d_N/d_S ratios ranged from 0.062 to 0.313 for influenza A lineages, and 0.257 to 0.259 for influenza B lineages. Structural analyses revealed that all positively selected sites are at the surface of the NA protein, with a number of sites found to be important for host antibody and drug binding.

Conclusions/Significance

The divergence into influenza type A and B from a putative ancestral NA was followed by the divergence of type A into nine NA subtypes, of which 23 lineages subsequently diverged. This study provides a better understanding of influenza NA lineages and their evolutionary dynamics, which may facilitate early detection of newly emerging influenza viruses and thus improve influenza surveillance.

Does history repeat itself? Wavelets and the phylodynamics of influenza A

Unprecedented global surveillance of viruses will result in massive sequence data sets that require new statistical methods. These data sets press the limits of Bayesian phylogenetics as the high-dimensional parameters that comprise a phylogenetic tree increase the already sizable computational burden of these techniques. This burden often results in partitioning the data set, for example, by gene, and inferring the evolutionary dynamics of each partition independently, a compromise that results in stratified analyses that depend only on data within a given partition. However, parameter estimates inferred from these stratified models are likely strongly correlated, considering they rely on data from a single data set. To overcome this shortfall, we exploit the existing Monte Carlo realizations from stratified Bayesian analyses to efficiently estimate a nonparametric hierarchical wavelet-based model and learn about the time-varying parameters of effective population size that reflect levels of genetic diversity across all partitions simultaneously. Our methods are applied to complete genome influenza A sequences that span 13 years. We find that broad peaks and trends, as opposed to seasonal spikes, in the effective population size history distinguish individual segments from the complete genome. We also address hypotheses regarding intersegment dynamics within a formal statistical framework that accounts for correlation between segment-specific parameters.

Influenza emergence in the face of evolutionary constraints

Different influenza subtypes can evolve at very different rates, but the causes are not well understood. In this paper, we explore whether differences in transmissibility between subtypes can play a role if there are fitness constraints on antigenic evolution. We investigate the problem using a mathematical model that separates the interaction of strains through cross-immunity from the process of emergence for new antigenic variants. Evolutionary constraints are also included with antigenic mutation incurring a fitness cost. We show that the transmissibility of a strain can become disproportionately important in dictating the rate of antigenic drift: strains that spread only slightly more easily can have a much higher rate of emergence. Further, we see that the effect continues when vaccination is considered; a small increase in the rate of transmission can make it much harder to control the frequency at which new strains emerge. Our results not only highlight the importance of considering both transmission and fitness constraints when modelling influenza evolution, but may also help in understanding the differences between the emergence of H1N1 and H3N2 subtypes.

Evolutionary dynamics of influenza A nucleoprotein (NP) lineages revealed by large-scale sequence analyses

Influenza A viral nucleoprotein (NP) plays a critical role in virus replication and host adaptation, however, the underlying molecular evolutionary dynamics of NP lineages are less well-understood. In this study, large-scale analyses of 5094 NP nucleotide sequences revealed eight distinct evolutionary lineages, including three host-specific lineages (human, classical swine and equine), two cross-host lineages (Eurasian avian-like swine and swine-origin human pandemic H1N1 2009) and three geographically isolated avian lineages (Eurasian, North American and Oceanian). The average nucleotide substitution rate of the NP lineages was estimated to be 2.4 × 10(-3) substitutions per site per year, with the highest value observed in pandemic H1N1 2009 (3.4 × 10(-3)) and the lowest in equine (0.9 × 10(-3)). The estimated time of most recent common ancestor (TMRCA) for each lineage demonstrated that the earliest human lineage was derived around 1906, and the latest pandemic H1N1 2009 lineage dated back to December 17, 2008. A marked time gap was found between the times when the viruses emerged and were first sampled, suggesting the crucial role for long-term surveillance of newly emerging viruses. The selection analyses showed that human lineage had six positive selection sites, whereas pandemic H1N1 2009, classical swine, Eurasian avian and Eurasian swine had only one or two sites. Protein structure analyses revealed several positive selection sites located in epitope regions or host adaptation regions, indicating strong adaptation to host immune system pressures in influenza viruses. Along with previous studies, this study provides new insights into the evolutionary dynamics of influenza A NP lineages. Further lineage analyses of other gene segments will allow better understanding of influenza A virus evolution and assist in the improvement of global influenza surveillance.

Host shifts and molecular evolution of H7 avian influenza virus hemagglutinin

Evolutionary consequences of host shifts represent a challenge to identify the mechanisms involved in the emergence of influenza A (IA) viruses. In this study we focused on the evolutionary history of H7 IA virus in wild and domestic birds, with a particular emphasis on host shifts consequences on the molecular evolution of the hemagglutinin (HA) gene. Based on a dataset of 414 HA nucleotide sequences, we performed an extensive phylogeographic analysis in order to identify the overall genetic structure of H7 IA viruses. We then identified host shift events and investigated viral population dynamics in wild and domestic birds, independently. Finally, we estimated changes in nucleotide substitution rates and tested for positive selection in the HA gene. A strong association between the geographic origin and the genetic structure was observed, with four main clades including viruses isolated in North America, South America, Australia and Eurasia-Africa. We identified ten potential events of virus introduction from wild to domestic birds, but little evidence for spillover of viruses from poultry to wild waterbirds. Several sites involved in host specificity (addition of a glycosylation site in the receptor binding domain) and virulence (insertion of amino acids in the cleavage site) were found to be positively selected in HA nucleotide sequences, in genetically unrelated lineages, suggesting parallel evolution for the HA gene of IA viruses in domestic birds. These results highlight that evolutionary consequences of bird host shifts would need to be further studied to understand the ecological and molecular mechanisms involved in the emergence of domestic bird-adapted viruses.

Genetic diversity in the 3' terminal 4.7-kb region of grapevine leafroll-associated virus 3

Grapevine leafroll-associated virus 3 (GLRaV-3; Ampelovirus, Closteroviridae), associated with grapevine leafroll disease, is an important pathogen found across all major grape-growing regions of the world. The genetic diversity of GLRaV-3 in Napa Valley, CA, was studied by sequencing 4.7 kb in the 3' terminal region of 50 isolates obtained from Vitis vinifera 'Merlot'. GLRaV-3 isolates were subdivided into four distinct phylogenetic clades. No evidence of positive selection was observed in the data set, although neutral selection (ratio of nonsynonymous to synonymous substitution rates = 1.1) was observed in one open reading frame (ORF 11, p4). Additionally, the four clades had variable degrees of overall nucleotide diversity. Moreover, no geographical structure among isolates was observed, and isolates belonging to different phylogenetic clades were found in distinct vineyards, with one exception. Considered with the evidence of purifying selection (i.e., against deleterious mutations), these data indicate that the population of GLRaV-3 in Napa Valley is not expanding and its effective population size is not increasing. Furthermore, research on the biological characterization of GLRaV-3 strains might provide valuable insights on the biology of this species that may have epidemiological relevance.

Implications of global and regional patterns of highly pathogenic avian influenza virus H5N1 clades for risk management

This paper analyses the publicly available data on the distribution and evolution of highly pathogenic avian influenza virus (HPAIV) H5N1 clades, whilst acknowledging the biases resulting from the non-random selection of isolates for gene sequencing. The data indicate molecular heterogeneity in the global distribution of HPAIV H5N1, in particular in different parts of East and Southeast Asia. Analysis of the temporal pattern of haemagglutinin clade data shows a progression from clade 0 (the 'dominant' clade between 1996 and 2002) to clade 1 (2003-2005) and then to clade 2.3.4 (2005 onwards). This process continuously produces variants, depending on the frequency of virus multiplication in the host population, which is influenced by geographical variation in poultry density, poultry production systems and also HPAI risk management measures such as vaccination. Increased multilateral collaboration needs to focus on developing enhanced disease surveillance and control targeted at evolutionary 'hotspots'.

Influenza genome diversity and evolution

The influenza viruses contain highly variable genomes and are able to infect a wide range of host species. Large-scale sequencing projects have collected abundant influenza sequence data for assessing influenza genome diversity and evolution. This work reviews current influenza sequence databases characteristics and statistics, as well as recent studies utilizing these databases to unravel influenza virus diversity and evolution. Also discussed are the newest deep sequencing methods and their applications to influenza virus research.

PhEVER: a database for the global exploration of virus-host evolutionary relationships

Fast viral adaptation and the implication of this rapid evolution in the emergence of several new infectious diseases have turned this issue into a major challenge for various research domains. Indeed, viruses are involved in the development of a wide range of pathologies and understanding how viruses and host cells interact in the context of adaptation remains an open question. In order to provide insights into the complex interactions between viruses and their host organisms and namely in the acquisition of novel functions through exchanges of genetic material, we developed the PhEVER database. This database aims at providing accurate evolutionary and phylogenetic information to analyse the nature of virus-virus and virus-host lateral gene transfers. PhEVER (http://pbil.univ-lyon1.fr/databases/phever) is a unique database of homologous families both (i) between sequences from different viruses and (ii) between viral sequences and sequences from cellular organisms. PhEVER integrates extensive data from up-to-date completely sequenced genomes (2426 non-redundant viral genomes, 1007 non-redundant prokaryotic genomes, 43 eukaryotic genomes ranging from plants to vertebrates) and offers a clustering of proteins into homologous families containing at least one viral sequences, as well as alignments and phylogenies for each of these families. Public access to PhEVER is available through its webpage and through all dedicated ACNUC retrieval systems.

Microevolution of canine influenza virus in shelters and its molecular epidemiology in the United States

Canine influenza virus (CIV) emerged around 2000 when an equine influenza virus (EIV) was transmitted to dogs in Florida. After 2003, the canine virus was carried by infected greyhounds to various parts of the United States and then became established in several large animal shelters, where it has continued to circulate. To better understand the evolution of CIV since its emergence, and particularly its microevolution in spatially restricted populations, we examined multiple gene segments of CIV from dogs resident in two large animal shelters in New York City during the period 2006 to 2009. In particular, we focused on viruses circulating in the two shelters in 2008 and 2009, which we found shared a common ancestor. While viruses in each shelter were generally monophyletic, we observed some gene flow between them. These shelter sequences were compared to earlier CIV isolates. The shelter viruses differed in 1 to 6 amino acids in each gene segment compared to viruses isolated in Florida between 2003 and 2005 and in Colorado in 2006 and 2008. A comparison of the sequences of equine and canine viruses revealed amino acid replacements that distinguished the viruses from the two hosts, but no clear evidence of positive selection indicative of host adaptation was detected, suggesting that any host range adaptation in CIV occurred early in the emergence of this virus or even before it transferred to dogs.

Influenza virus evolution, host adaptation, and pandemic formation

Newly emerging or "re-emerging" viral diseases continue to pose significant global public health threats. Prototypic are influenza viruses that are major causes of human respiratory infections and mortality. Influenza viruses can cause zoonotic infections and adapt to humans, leading to sustained transmission and emergence of novel viruses. Mechanisms by which viruses evolve in one host, cause zoonotic infection, and adapt to a new host species remain unelucidated. Here, we review the evolution of influenza A viruses in their reservoir hosts and discuss genetic changes associated with introduction of novel viruses into humans, leading to pandemics and the establishment of seasonal viruses.

Microevolution of Helicobacter pylori during prolonged infection of single hosts and within families

Our understanding of basic evolutionary processes in bacteria is still very limited. For example, multiple recent dating estimates are based on a universal inter-species molecular clock rate, but that rate was calibrated using estimates of geological dates that are no longer accepted. We therefore estimated the short-term rates of mutation and recombination in Helicobacter pylori by sequencing an average of 39,300 bp in 78 gene fragments from 97 isolates. These isolates included 34 pairs of sequential samples, which were sampled at intervals of 0.25 to 10.2 years. They also included single isolates from 29 individuals (average age: 45 years) from 10 families. The accumulation of sequence diversity increased with time of separation in a clock-like manner in the sequential isolates. We used Approximate Bayesian Computation to estimate the rates of mutation, recombination, mean length of recombination tracts, and average diversity in those tracts. The estimates indicate that the short-term mutation rate is 1.4×10⁻⁶ (serial isolates) to 4.5×10⁻⁶ (family isolates) per nucleotide per year and that three times as many substitutions are introduced by recombination as by mutation. The long-term mutation rate over millennia is 5–17-fold lower, partly due to the removal of non-synonymous mutations due to purifying selection. Comparisons with the recent literature show that short-term mutation rates vary dramatically in different bacterial species and can span a range of several orders of magnitude.

Mutation rates in bacteria have generally been considered to be much slower than in viruses. This is partly because estimates of long-term mutation rates for the evolution of distinct species have been inappropriately used for dating divergence within species. Furthermore, the most commonly used long-term mutation rate is based on geological dates that are no longer accepted. In addition, only few short-term mutation rates have been calculated within bacterial species, and these differ with the species by several orders of magnitude. Here, we provide robust estimates for short-term mutation and recombination rates within Helicobacter pylori, a bacterium that commonly infects the human gastric mucosa, based on serial isolates from long-term infections and on differences between isolates from multiple family members. These short-term mutation rates are 5–17-fold faster than long-term mutation rates in H. pylori that have been calibrated by parallel ancient migrations of humans. Short-term mutation rates in bacteria, including those for H. pylori, can be quite fast, partially overlapping with those for viruses. Future calculations of ages of bacterial species will need to account for dramatic differences in mutation rate between species and for dramatic differences between short- and long-term mutation rates.