A maximum likelihood method for detecting directional evolution in protein sequences and its application to influenza A virus

Characterization of Conserved Evolution in H7N9 Avian Influenza Virus Prior Mass Vaccination

Vaccination is crucial for the prevention and mitigation of avian influenza infections in China. The inactivated H7N9 vaccine, when administered to poultry, significantly lowers the risk of infection among both poultry and humans, while also markedly decreasing the prevalence of H7N9 detections. Highly pathogenic (HP) H7N9 viruses occasionally appear, whereas their low pathogenicity (LP) counterparts have been scarcely detected since 2018. However, these contributing factors remain poorly understood. We conducted an exploratory investigation of the mechanics via the application of comprehensive bioinformatic approaches. We delineated the Yangtze River Delta (YRD) H7N9 lineage into 5 clades (YRD-A to E). Our findings highlight the emergence and peak occurrence of the LP H7N9-containing YRD-E clade during the 5th epidemic wave in China's primary poultry farming areas. A more effective control of LP H7N9 through vaccination was observed compared to that of its HP H7N9 counterpart. YRD-E exhibited a tardy evolutionary trajectory, denoted by the conservation of its genetic and antigenic variation. Our analysis of YRD-E revealed only minimal amino acid substitutions along its phylogenetic tree and a few selective sweep mutations since 2016. In terms of epidemic fitness, the YRD-E was measured to be lower than that of the HP variants. Collectively, these findings underscore the conserved evolutionary patterns distinguishing the YRD-E. Given the conservation presented in its evolutionary patterns, the YRD-E LP H7N9 is hypothesized to be associated with a reduction following the mass vaccination in a relatively short period owing to its lower probability of antigenic variation that might affect vaccine efficiency.

Epitopes in the HA and NA of H5 and H7 avian influenza viruses that are important for antigenic drift

Avian influenza viruses evolve antigenically to evade host immunity. Two influenza A virus surface glycoproteins, the haemagglutinin and neuraminidase, are the major targets of host immunity and undergo antigenic drift in response to host pre-existing humoral and cellular immune responses. Specific sites have been identified as important epitopes in prominent subtypes such as H5 and H7, which are of animal and public health significance due to their panzootic and pandemic potential. The haemagglutinin is the immunodominant immunogen, it has been extensively studied, and the antigenic reactivity is closely monitored to ensure candidate vaccine viruses are protective. More recently, the neuraminidase has received increasing attention for its role as a protective immunogen. The neuraminidase is expressed at a lower abundance than the haemagglutinin on the virus surface but does elicit a robust antibody response. This review aims to compile the current information on haemagglutinin and neuraminidase epitopes and immune escape mutants of H5 and H7 highly pathogenic avian influenza viruses. Understanding the evolution of immune escape mutants and the location of epitopes is critical for identification of vaccine strains and development of broadly reactive vaccines that can be utilized in humans and animals.

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.

Next-generation development and application of codon model in evolution

To date, numerous nucleotide, amino acid, and codon substitution models have been developed to estimate the evolutionary history of any sequence/organism in a more comprehensive way. Out of these three, the codon substitution model is the most powerful. These models have been utilized extensively to detect selective pressure on a protein, codon usage bias, ancestral reconstruction and phylogenetic reconstruction. However, due to more computational demanding, in comparison to nucleotide and amino acid substitution models, only a few studies have employed the codon substitution model to understand the heterogeneity of the evolutionary process in a genome-scale analysis. Hence, there is always a question of how to develop more robust but less computationally demanding codon substitution models to get more accurate results. In this review article, the authors attempted to understand the basis of the development of different types of codon-substitution models and how this information can be utilized to develop more robust but less computationally demanding codon substitution models. The codon substitution model enables to detect selection regime under which any gene or gene region is evolving, codon usage bias in any organism or tissue-specific region and phylogenetic relationship between different lineages more accurately than nucleotide and amino acid substitution models. Thus, in the near future, these codon models can be utilized in the field of conservation, breeding and medicine.

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.

RASCL: Rapid Assessment of Selection in CLades through molecular sequence analysis

An important unmet need revealed by the COVID-19 pandemic is the near-real-time identification of potentially fitness-altering mutations within rapidly growing SARS-CoV-2 lineages. Although powerful molecular sequence analysis methods are available to detect and characterize patterns of natural selection within modestly sized gene-sequence datasets, the computational complexity of these methods and their sensitivity to sequencing errors render them effectively inapplicable in large-scale genomic surveillance contexts. Motivated by the need to analyze new lineage evolution in near-real time using large numbers of genomes, we developed the Rapid Assessment of Selection within CLades (RASCL) pipeline. RASCL applies state of the art phylogenetic comparative methods to evaluate selective processes acting at individual codon sites and across whole genes. RASCL is scalable and produces automatically updated regular lineage-specific selection analysis reports: even for lineages that include tens or hundreds of thousands of sampled genome sequences. Key to this performance is (i) generation of automatically subsampled high quality datasets of gene/ORF sequences drawn from a selected "query" viral lineage; (ii) contextualization of these query sequences in codon alignments that include high-quality "background" sequences representative of global SARS-CoV-2 diversity; and (iii) the extensive parallelization of a suite of computationally intensive selection analysis tests. Within hours of being deployed to analyze a novel rapidly growing lineage of interest, RASCL will begin yielding JavaScript Object Notation (JSON)-formatted reports that can be either imported into third-party analysis software or explored in standard web-browsers using the premade RASCL interactive data visualization dashboard. By enabling the rapid detection of genome sites evolving under different selective regimes, RASCL is well-suited for near-real-time monitoring of the population-level selective processes that will likely underlie the emergence of future variants of concern in measurably evolving pathogens with extensive genomic surveillance.

Thermostability engineering of an inulin fructotransferase for the biosynthesis of difructose anhydride I

The thermostability of enzymes is an essential factor that performs a vital role during practical applications. Inulin fructotransferases can efficiently convert inulin into bio-functional difructose anhydrides (DFAs). The present study aimed to improve the thermostability of a previously reported inulin fructotransferase, SpIFTase, and apply it to the biosynthesis of DFA I. In silico rational design was used to predict mutation sites, based on sequential and structural information. Two triple-site mutants, Q69L/Q234L/K310G and E201I/Q234L/K310G, were characterized and exhibited enhanced thermostability with approximately 5 °C higher in melting temperature (T_m), respectively, and a 45-fold longer half-life (t_1/2) at 70 °C, compared to that of SpIFTase. Molecular dynamic simulations and elaborate structural analysis suggested that the combinations of hydrophobic interaction, electrostatic potential distribution, and decreased flexibility via stabilization of loops and α-helix improved the thermostability of SpIFTase. Additionally, the promising mutants exhibited great potential to the industrial production of DFA I.

Deep phylogenetic-based clustering analysis uncovers new and shared mutations in SARS-CoV-2 variants as a result of directional and convergent evolution

Nearly two decades after the last epidemic caused by a severe acute respiratory syndrome coronavirus (SARS-CoV), newly emerged SARS-CoV-2 quickly spread in 2020 and precipitated an ongoing global public health crisis. Both the continuous accumulation of point mutations, owed to the naturally imposed genomic plasticity of SARS-CoV-2 evolutionary processes, as well as viral spread over time, allow this RNA virus to gain new genetic identities, spawn novel variants and enhance its potential for immune evasion. Here, through an in-depth phylogenetic clustering analysis of upwards of 200,000 whole-genome sequences, we reveal the presence of previously unreported and hitherto unidentified mutations and recombination breakpoints in Variants of Concern (VOC) and Variants of Interest (VOI) from Brazil, India (Beta, Eta and Kappa) and the USA (Beta, Eta and Lambda). Additionally, we identify sites with shared mutations under directional evolution in the SARS-CoV-2 Spike-encoding protein of VOC and VOI, tracing a heretofore-undescribed correlation with viral spread in South America, India and the USA. Our evidence-based analysis provides well-supported evidence of similar pathways of evolution for such mutations in all SARS-CoV-2 variants and sub-lineages. This raises two pivotal points: (i) the co-circulation of variants and sub-lineages in close evolutionary environments, which sheds light onto their trajectories into convergent and directional evolution, and (ii) a linear perspective into the prospective vaccine efficacy against different SARS-CoV-2 strains.

Rapid epidemic expansion of the SARS-CoV-2 Omicron variant in southern Africa

The SARS-CoV-2 epidemic in southern Africa has been characterized by three distinct waves. The first was associated with a mix of SARS-CoV-2 lineages, while the second and third waves were driven by the Beta (B.1.351) and Delta (B.1.617.2) variants, respectively1-3. In November 2021, genomic surveillance teams in South Africa and Botswana detected a new SARS-CoV-2 variant associated with a rapid resurgence of infections in Gauteng province, South Africa. Within three days of the first genome being uploaded, it was designated a variant of concern (Omicron, B.1.1.529) by the World Health Organization and, within three weeks, had been identified in 87 countries. The Omicron variant is exceptional for carrying over 30 mutations in the spike glycoprotein, which are predicted to influence antibody neutralization and spike function4. Here we describe the genomic profile and early transmission dynamics of Omicron, highlighting the rapid spread in regions with high levels of population immunity.

RASCL: Rapid Assessment Of SARS-CoV-2 Clades Through Molecular Sequence Analysis

An important component of efforts to manage the ongoing COVID19 pandemic is the R apid A ssessment of how natural selection contributes to the emergence and proliferation of potentially dangerous S ARS-CoV-2 lineages and CL ades (RASCL). The RASCL pipeline enables continuous comparative phylogenetics-based selection analyses of rapidly growing clade-focused genome surveillance datasets, such as those produced following the initial detection of potentially dangerous variants. From such datasets RASCL automatically generates down-sampled codon alignments of individual genes/ORFs containing contextualizing background reference sequences, analyzes these with a battery of selection tests, and outputs results as both machine readable JSON files, and interactive notebook-based visualizations.

AVAILABILITY

RASCL is available from a dedicated repository at https://github.com/veg/RASCL and as a Galaxy workflow https://usegalaxy.eu/u/hyphy/w/rascl . Existing clade/variant analysis results are available here: https://observablehq.com/@aglucaci/rascl .

CONTACT

Dr. Sergei L Kosakovsky Pond ( spond@temple.edu ).

SUPPLEMENTARY INFORMATION

N/A.

[The molecular evolution of Dabie bandavirus (Phenuiviridae: Bandavirus: Dabie bandavirus), the agent of severe fever with thrombocytopenia syndrome]

Since the Dabie bandavirus (DBV; former SFTS virus, SFTSV) was identified, the epidemics of severe fever with thrombocytopenic syndrome (SFTS) caused by this virus have occurred in several countries in East Asia. The rapid increase in incidence indicates that this infectious agent has a pandemic potential and poses an imminent global public health threat.The analysis of molecular evolution of SFTS agent that includes its variants isolated in China, Japan and South Korea was performed in this review. The evolution rate of DBV and the estimated dates of existence of the common ancestor were ascertained, and the possibility of reassortation was demonstrated.The evolutionary rates of DBV genome segments were estimated to be 2.28 × 10-4 nucleotides/site/year for S-segment, 2.42 × 10-4 for M-segment, and 1.19 × 10-4 for L-segment. The positions of positive selection were detected in the viral genome.Phylogenetic analyses showed that virus may be divided into two clades, containing six different genotypes. The structures of phylogenetic trees for S-, M- and L-segments showed that all genotypes originate from the common ancestor.Data of sequence analysis suggest that DBV use several mechanisms to maintain the high level of its genetic diversity. Understanding the phylogenetic factors that determine the virus transmission is important for assessing the epidemiological characteristics of the disease and predicting its possible outbreaks.

Allelic diversity at BoLA DRB3 locus and association with predisposition to clinical mastitis in indicus and crossbred cattle

Bovine lymphocyte antigen (BoLA) DRB3 locus in healthy and mastitis affected cattle has been genotyped by a polymerase chain reaction and restriction fragment length polymorphisms (PCR-RLFP) using RsaI restriction enzyme, followed by sequencing. In 130 farm animals, 25 BoLA DRB3 alleles have been detected by PCR-RFLP. Three distinct allelic patterns significantly associated with mastitis in Karan Fries crossbred and Sahiwal indicus cattle have been identified, whereas, four other allelic patterns were significantly high in frequency among healthy animals. Sequencing of RFLP genotypes revealed 25 and 47 alleles among healthy Sahiwal and Karan Fries, respectively, while 17 and 38 patterns observed in mastitis affected Sahiwal and Karan Fries animals, respectively. From Tajima's D-test of neutrality, it was concluded that alleles associated with mastitis were expanding in the population, whereas those of healthy were under contraction. Phylogenetic analysis carried out to delineate the evolutionary relationship of the farm and field animals at DRB3 locus, differentiating allelic patterns into six different clusters. Among the phylogenetic lineages, five patterns DRB3*028:01, DRB3*011:03, DRB3*031:01, DRB3*001:01 and DRB3*043:01, were previously reported, whereas one novel allelic variant was observed in indicus and crossbred cattle. This information will help in further exploring the association between BoLA-DRB3 genetic diversity and disease resistance in distinct cattle breeds, important in designing breeding strategies for increasing the distribution of favorable alleles.

Characterization of changes in the hemagglutinin that accompanied the emergence of H3N2/1968 pandemic influenza viruses

The hemagglutinin (HA) of A/H3N2 pandemic influenza viruses (IAVs) of 1968 differed from its inferred avian precursor by eight amino acid substitutions. To determine their phenotypic effects, we studied recombinant variants of A/Hong Kong/1/1968 virus containing either human-type or avian-type amino acids in the corresponding positions of HA. The precursor HA displayed receptor binding profile and high conformational stability typical for duck IAVs. Substitutions Q226L and G228S, in addition to their known effects on receptor specificity and replication, marginally decreased HA stability. Substitutions R62I, D63N, D81N and N193S reduced HA binding avidity. Substitutions R62I, D81N and A144G promoted viral replication in human airway epithelial cultures. Analysis of HA sequences revealed that substitutions D63N and D81N accompanied by the addition of N-glycans represent common markers of avian H3 HA adaptation to mammals. Our results advance understanding of genotypic and phenotypic changes in IAV HA required for avian-to-human adaptation and pandemic emergence.

Impact of viral features, host jumps and phylogeography on the rapid evolution of Aleutian mink disease virus (AMDV)

Aleutian mink disease virus (AMDV) is one the most relevant pathogens of domestic mink, where it can cause significant economic losses, and wild species, which are considered a threat to mink farms. Despite their relevance, many aspects of the origin, evolution, and geographic and host spreading patterns of AMDV have never been investigated on a global scale using a comprehensive biostatistical approach. The present study, benefitting from a large dataset of sequences collected worldwide and several phylodynamic-based approaches, demonstrates the ancient origin of AMDV and its broad, unconstrained circulation from the initial intercontinental spread to the massive among-country circulation, especially within Europe, combined with local persistence and evolution. Clear expansion of the viral population size occurred over time until more effective control measures started to be applied. The role of frequent changes in epidemiological niches, including different hosts, in driving the high nucleotide and amino acid evolutionary rates was also explored by comparing the strengths of selective pressures acting on different populations. The obtained results suggest that the viral passage among locations and between wild and domesticated animals poses a double threat to farm profitability and animal welfare and health, which is particularly relevant for endangered species. Therefore, further efforts must be made to limit viral circulation and to refine our knowledge of factors enhancing AMDV spread, particularly at the wild-domestic interface.

The virus-host interface: Molecular interactions of Alphacoronavirus-1 variants from wild and domestic hosts with mammalian aminopeptidase N

The Alphacoronavirus‐1 species include viruses that infect numerous mammalian species. To better understand the wide host range of these viruses, better knowledge on the molecular determinants of virus–host cell entry mechanisms in wildlife hosts is essential. We investigated Alphacoronavirus‐1 infection in carnivores using long‐term data on Serengeti spotted hyenas (Crocuta crocuta) and molecular analyses guided by the tertiary structure of the viral spike (S) attachment protein's interface with the host receptor aminopeptidase N (APN). We sequenced the complete 3′‐end region of the genome of nine variants from wild African carnivores, plus the APN gene of 15 wild carnivore species. Our results revealed two outbreaks of Alphacoronavirus‐1 infection in spotted hyenas associated with genetically distinct canine coronavirus type II (CCoVII) variants. Within the receptor binding domain (RBD) of the S gene the residues that directly bind to the APN receptor were conserved in all variants studied, even those infecting phylogenetically diverse host taxa. We identified a variable region within RBD located next to a region that directly interacts with the APN receptor. Two residues within this variable region were under positive selection in hyena variants, indicating that both sites were associated with adaptation of CCoVII to spotted hyena APN. Analysis of APN sequences revealed that most residues that interact with the S protein are conserved in wild carnivores, whereas some adjacent residues are highly variable. Of the variable residues, four that are critical for virus–host binding were under positive selection and may modulate the efficiency of virus attachment to carnivore APN.

Relative Efficiencies of Simple and Complex Substitution Models in Estimating Divergence Times in Phylogenomics

The conventional wisdom in molecular evolution is to apply parameter-rich models of nucleotide and amino acid substitutions for estimating divergence times. However, the actual extent of the difference between time estimates produced by highly complex models compared with those from simple models is yet to be quantified for contemporary data sets that frequently contain sequences from many species and genes. In a reanalysis of many large multispecies alignments from diverse groups of taxa, we found that the use of the simplest models can produce divergence time estimates and credibility intervals similar to those obtained from the complex models applied in the original studies. This result is surprising because the use of simple models underestimates sequence divergence for all the data sets analyzed. We found three fundamental reasons for the observed robustness of time estimates to model complexity in many practical data sets. First, the estimates of branch lengths and node-to-tip distances under the simplest model show an approximately linear relationship with those produced by using the most complex models applied on data sets with many sequences. Second, relaxed clock methods automatically adjust rates on branches that experience considerable underestimation of sequence divergences, resulting in time estimates that are similar to those from complex models. And, third, the inclusion of even a few good calibrations in an analysis can reduce the difference in time estimates from simple and complex models. The robustness of time estimates to model complexity in these empirical data analyses is encouraging, because all phylogenomics studies use statistical models that are oversimplified descriptions of actual evolutionary substitution processes.

Identifying Potentially Beneficial Genetic Mutations Associated with Monophyletic Selective Sweep and a Proof-of-Concept Study with Viral Genetic Data

Genetic mutations play a central role in evolution. For a significantly beneficial mutation, a one-time mutation event suffices for the species to prosper and predominate through the process called “monophyletic selective sweep.” However, existing methods that rely on counting the number of mutation events to detect selection are unable to find such a mutation in selective sweep. We here introduce a method to detect mutations at the single amino acid/nucleotide level that could be responsible for monophyletic selective sweep evolution. The method identifies a genetic signature associated with selective sweep using the population genetic test statistic Tajima’s D. We applied the algorithm to ebolavirus, influenza A virus, and severe acute respiratory syndrome coronavirus 2 to identify known biologically significant mutations and unrecognized mutations associated with potential selective sweep. The method can detect beneficial mutations, possibly leading to discovery of previously unknown biological functions and mechanisms related to those mutations.

IMPORTANCE In biology, research on evolution is important to understand the significance of genetic mutation. When there is a significantly beneficial mutation, a population of species with the mutation prospers and predominates, in a process called “selective sweep.” However, there are few methods that can find such a mutation causing selective sweep from genetic data. We here introduce a novel method to detect such mutations. Applying the method to the genomes of ebolavirus, influenza viruses, and the novel coronavirus, we detected known biologically significant mutations and identified mutations the importance of which is previously unrecognized. The method can deepen our understanding of molecular and evolutionary biology.

Parallel evolution in the emergence of highly pathogenic avian influenza A viruses

Parallel molecular evolution and adaptation are important phenomena commonly observed in viruses. Here, we exploit parallel molecular evolution to understand virulence evolution in avian influenza viruses (AIV). Highly-pathogenic AIVs evolve independently from low-pathogenic ancestors via acquisition of polybasic cleavage sites. Why some AIV lineages but not others evolve in this way is unknown. We hypothesise that the parallel emergence of highly-pathogenic AIV may be facilitated by permissive or compensatory mutations occurring across the viral genome. We combine phylogenetic, statistical and structural approaches to discover parallel mutations in AIV genomes associated with the highly-pathogenic phenotype. Parallel mutations were screened using a statistical test of mutation-phenotype association and further evaluated in the contexts of positive selection and protein structure. Our resulting mutational panel may help to reveal new links between virulence evolution and other traits, and raises the possibility of predicting aspects of AIV evolution.

Positive Selection of ORF1ab, ORF3a, and ORF8 Genes Drives the Early Evolutionary Trends of SARS-CoV-2 During the 2020 COVID-19 Pandemic

In this study, we analyzed full-length SARS-CoV-2 genomes from multiple countries to determine early trends in the evolutionary dynamics of the novel COVID-19 pandemic. Results indicated SARS-CoV-2 evolved early into at least three phylogenetic groups, characterized by positive selection at specific residues of the accessory proteins ORF3a and ORF8. Also, we are reporting potential relevant sites under positive selection at specific sites of non-structural proteins nsp6 and helicase. Our analysis of co-evolution showed evidence of epistatic interactions among sites in the genome that may be important in the generation of variants adapted to humans. These observations might impact not only public health but also suggest that more studies are needed to understand the genetic mechanisms that may affect the development of therapeutic and preventive tools, like antivirals and vaccines. Collectively, our results highlight the identification of ongoing selection even in a scenario of conserved sequences collected over the first 3 months of this pandemic.

Limited Cross-Protection Provided by Prior Infection Contributes to High Prevalence of Influenza D Viruses in Cattle

Since its detection in swine, influenza D virus (IDV) has been shown to be present in multiple animal hosts, and bovines have been identified as its natural reservoir. However, it remains unclear how IDVs emerge, evolve, spread, and maintain in bovine populations. Through multiple years of virological and serological surveillance in a single order-buyer cattle facility in Mississippi, we showed consistently high seroprevalence of IDVs in cattle and recovered a total of 32 IDV isolates from both healthy and sick animals, including those with antibodies against IDV. Genomic analyses of these isolates along with those isolated from other areas showed that active genetic reassortment occurred in IDV and that five reassortants were identified in the Mississippian facility. Two antigenic groups were identified through antigenic cartography analyses for these 32 isolates and representative IDVs from other areas. Remarkably, existing antibodies could not protect cattle from experimental reinfection with IDV. Additional phenotypic analyses demonstrated variations in growth dynamics and pathogenesis in mice between viruses independent of genomic constellation. In summary, this study suggests that, in addition to epidemiological factors, the ineffectiveness of preexisting immunity and cocirculation of a diverse viral genetic pool could facilitate its high prevalence in animal populations.IMPORTANCE Influenza D viruses (IDVs) are panzootic in multiple animal hosts, but the underlying mechanism is unclear. Through multiple years of surveillance in the same order-buyer cattle facility, 32 IDV isolates were recovered from both healthy and sick animals, including those with evident antibodies against IDV. Active reassortment occurred in the cattle within this facility and in those across other areas, and multiple reassortants cocirculated in animals. These isolates are shown with a large extent of phenotypic diversity in replication efficiency and pathogenesis but little in antigenic properties. Animal experiments demonstrated that existing antibodies could not protect cattle from experimental reinfection with IDV. This study suggests that, in addition to epidemiological factors, limited protection from preexisting immunity against IDVs in cattle herds and cocirculation of a diverse viral genetic pool likely facilitate the high prevalence of IDVs in animal populations.

Characterization of the FMDV-serotype-O isolates collected during 1962 and 1997 discloses new topotypes, CEY-1 and WCSA-1, and six new lineages

The genetic diversity of the FMD viruses collected from the outbreaks during the second half of the 20th Century in Sri Lanka was assessed in the present study. We sequenced the VP1 genomic region of the samples collected during FMDV epidemics caused by serotype O in Sri Lanka during 1962 and 1997. For comparison, we sequenced the VP1 of the related viral isolates collected from other Asian countries. We analyzed the VP1 sequences of the viral strains using the UPGMA method with uncorrected pairwise distances. Nucleotide divergence (ND) thresholds of 15%-20% and 5%-<15% were used to differentiate topotypes and lineages, respectively. We calibrated the divergence times and lineage-specific substitution rates using Bayesian-skyline models. Based on the ND estimations and phylogenetic relationships, we identified and named two new topotypes [CEYLON 1 (CEY-1) and WEST, CENTRAL AND SOUTH ASIA 1 (WCSA-1)] and six new lineages (Syr-62, Srl-77, Tur-69, May-78, Tai-87 and Bur-77) of serotype O. We believe that the novel topotypes and lineages named may have disappeared although they have similar substitution rates for epizootic outbreaks. Because the amino acid selection analysis revealed that the two topotypes and six lineages identified were under purifying selection during the outbreaks.

Genetic and Biochemical Mechanisms for Bacterial Lipid A Modifiers Associated with Polymyxin Resistance

Polymyxins are a group of detergent-like antimicrobial peptides that are the ultimate line of defense against carbapenem-resistant pathogens in clinical settings. Polymyxin resistance primarily originates from structural remodeling of lipid A anchored on bacterial surfaces. We integrate genetic, structural, and biochemical aspects of three major types of lipid A modifiers that have been shown to confer intrinsic colistin resistance. Namely, we highlight ArnT, a glycosyltransferase, EptA, a phosphoethanolamine transferase, and the AlmEFG tripartite system, which is restricted to EI Tor biotype of Vibrio cholerae O1. We also discuss the growing family of mobile colistin resistance (MCR) enzymes, each of which is analogous to EptA, and which pose great challenges to global public health.

The evolution and phylodynamics of serotype A and SAT2 foot-and-mouth disease viruses in endemic regions of Africa

Foot-and-mouth disease (FMD) is a major livestock disease with direct clinical impacts as well as indirect trade implications. Control through vaccination and stamping-out has successfully reduced or eradicated the disease from Europe and large parts of South America. However, sub-Saharan Africa remains endemically affected with 5/7 serotypes currently known to be circulating across the continent. This has significant implications both locally for livestock production and poverty reduction but also globally as it represents a major reservoir of viruses, which could spark new epidemics in disease free countries or vaccination zones. This paper describes the phylodynamics of serotypes A and SAT2 in Africa including recent isolates from Cameroon in Central Africa. We estimated the most recent common ancestor for serotype A was an East African virus from the 1930s (median 1937; HPD 1922-1950) compared to SAT2 which has a much older common ancestor from the early 1700s (median 1709; HPD 1502-1814). Detailed analysis of the different clades shows clearly that different clades are evolving and diffusing across the landscape at different rates with both serotypes having a particularly recent clade that is evolving and spreading more rapidly than other clades within their serotype. However, the lack of detailed sequence data available for Africa seriously limits our understanding of FMD epidemiology across the continent. A comprehensive view of the evolutionary history and dynamics of FMD viruses is essential to understand many basic epidemiological aspects of FMD in Africa such as the scale of persistence and the role of wildlife and thus the opportunities and scale at which vaccination and other controls could be applied. Finally we ask endemic countries to join the OIE/FAO supported regional networks and take advantage of new cheap technologies being rolled out to collect isolates and submit them to the World Reference Laboratory.

Evolutionary and phylodynamic analyses of Dengue virus serotype I in Guangdong Province, China, between 1985 and 2015

To determine the evolutionary and phylodynamic history of DENV-1 in Guangdong, the strains detected between 1985 and 2015 were determined with phylogenetic and Bayesian analyses of the E gene. Three DENV-1 genotypes (I, V, and VI) were circulating in Guangdong, and genotype I was detected most frequently. The evolutionary rate of DENV-1 was estimated to be 1.03 × 10^-3 nucleotide substitutions/site/year. The most recent ancestor of the viruses existed approximately 141 years ago. The observed epidemiological dynamics correlated with similar fluctuations in diversity, and the epidemiological dynamics of DENV-1 transmission reflect dramatic changes in the viral population sizes. Two recombination events were identified in those strains. The selection pressures were estimated and revealed an abundance of negatively selected sites but few positively selected sites. These data improve our understanding of the evolution and molecular epidemiology of DENV-1 and provide insights that will facilitate the surveillance and control of DENV-1.

Sweep Dynamics (SD) plots: Computational identification of selective sweeps to monitor the adaptation of influenza A viruses

Monitoring changes in influenza A virus genomes is crucial to understand its rapid evolution and adaptation to changing conditions e.g. establishment within novel host species. Selective sweeps represent a rapid mode of adaptation and are typically observed in human influenza A viruses. We describe Sweep Dynamics (SD) plots, a computational method combining phylogenetic algorithms with statistical techniques to characterize the molecular adaptation of rapidly evolving viruses from longitudinal sequence data. SD plots facilitate the identification of selective sweeps, the time periods in which these occurred and associated changes providing a selective advantage to the virus. We studied the past genome-wide adaptation of the 2009 pandemic H1N1 influenza A (pH1N1) and seasonal H3N2 influenza A (sH3N2) viruses. The pH1N1 influenza virus showed simultaneous amino acid changes in various proteins, particularly in seasons of high pH1N1 activity. Partially, these changes resulted in functional alterations facilitating sustained human-to-human transmission. In the evolution of sH3N2 influenza viruses, we detected changes characterizing vaccine strains, which were occasionally revealed in selective sweeps one season prior to the WHO recommendation. Taken together, SD plots allow monitoring and characterizing the adaptive evolution of influenza A viruses by identifying selective sweeps and their associated signatures.

Population Diversity and Collective Interactions during Influenza Virus Infection

Influenza A virus (IAV) continues to pose an enormous and unpredictable global public health threat, largely due to the continual evolution of escape from preexisting immunity and the potential for zoonotic emergence. Understanding how the unique genetic makeup and structure of IAV populations influences their transmission and evolution is essential for developing more-effective vaccines, therapeutics, and surveillance capabilities. Owing to their mutation-prone replicase and unique genome organization, IAV populations exhibit enormous amounts of diversity both in terms of sequence and functional gene content. Here, I review what is currently known about the genetic and genomic diversity present within IAV populations and how this diversity may shape the replicative and evolutionary dynamics of these viruses.

The Evolutionary Pathway to Virulence of an RNA Virus

Paralytic polio once afflicted almost half a million children each year. The attenuated oral polio vaccine (OPV) has enabled world-wide vaccination efforts, which resulted in nearly complete control of the disease. However, poliovirus eradication is hampered globally by epidemics of vaccine-derived polio. Here, we describe a combined theoretical and experimental strategy that describes the molecular events leading from OPV to virulent strains. We discover that similar evolutionary events occur in most epidemics. The mutations and the evolutionary trajectories driving these epidemics are replicated using a simple cell-based experimental setup where the rate of evolution is intentionally accelerated. Furthermore, mutations accumulating during epidemics increase the replication fitness of the virus in cell culture and increase virulence in an animal model. Our study uncovers the evolutionary strategies by which vaccine strains become pathogenic and provides a powerful framework for rational design of safer vaccine strains and for forecasting virulence of viruses. VIDEO ABSTRACT.

Antigenic Drift Defines a New D4 Subgenotype of Measles Virus

The measles virus hemagglutinin (MeV-H) protein is the main target of protective neutralizing antibodies. Using a panel of monoclonal antibodies (MAbs) that recognize known major antigenic sites in MeV-H, we identified a D4 genotype variant that escapes neutralization by MAbs targeting the neutralizing epitope (NE) antigenic site. By site-directed mutagenesis, L249P was identified as the critical mutation disrupting the NE in this genotype D4 variant. Forty-two available D4 genotype gene sequences were subsequently analyzed and divided into 2 groups according to the presence or absence of the L249P MeV-H mutation. Further analysis of the MeV-N gene sequences of these 2 groups confirmed that they represent clearly definable, sequence-divergent D4 subgenotypes, which we named subgenotypes D4.1 and D4.2. The subgenotype D4.1 MeVs were isolated predominantly in Kenya and Ethiopia, whereas the MAb-resistant subgenotype D4.2 MeVs were isolated predominantly in France and Great Britain, countries with higher vaccine coverage rates. Interestingly, D4.2 subgenotype viruses showed a trend toward diminished susceptibility to neutralization by human sera pooled from approximately 60 to 80 North American donors. Escape from MAb neutralization may be a powerful epidemiological surveillance tool to monitor the evolution of new MeV subgenotypes.IMPORTANCE Measles virus is a paradigmatic RNA virus, as the antigenic composition of the vaccination has not needed to be updated since its discovery. The vaccine confers protection by inducing neutralizing antibodies that interfere with the function of the hemagglutinin protein. Viral strains are indistinguishable serologically, although characteristic nucleotide sequences differentiate 24 genotypes. In this work, we describe a distant evolutionary branch within genotype D4. Designated subgenotype D4.2, this virus is distinguishable by neutralization with vaccine-induced monoclonal antibodies that target the neutralizing epitope (NE). The subgenotype D4.2 viruses have a higher predominance in countries with intermediary levels of vaccine coverage. Our studies demonstrate that subgenotype D4.2 lacks epitopes associated with half of the known antigenic sites, which significantly impacts our understanding of measles virus evolution.

Plant Actin-Depolymerizing Factors Possess Opposing Biochemical Properties Arising from Key Amino Acid Changes throughout Evolution

Functional divergence in paralogs is an important genetic source of evolutionary innovation. Actin-depolymerizing factors (ADFs) are among the most important actin binding proteins and are involved in generating and remodeling actin cytoskeletal architecture via their conserved F-actin severing or depolymerizing activity. In plants, ADFs coevolved with actin, but their biochemical properties are diverse. Unfortunately, the biochemical function of most plant ADFs and the potential mechanisms of their functional divergence remain unclear. Here, in vitro biochemical analyses demonstrated that all 11 ADF genes in Arabidopsis thaliana exhibit opposing biochemical properties. Subclass III ADFs evolved F-actin bundling (B-type) function from conserved F-actin depolymerizing (D-type) function, and subclass I ADFs have enhanced D-type function. By tracking historical mutation sites on ancestral proteins, several fundamental amino acid residues affecting the biochemical functions of these proteins were identified in Arabidopsis and various plants, suggesting that the biochemical divergence of ADFs has been conserved during the evolution of angiosperm plants. Importantly, N-terminal extensions on subclass III ADFs that arose from intron-sliding events are indispensable for the alteration of D-type to B-type function. We conclude that the evolution of these N-terminal extensions and several conserved mutations produced the diverse biochemical functions of plant ADFs from a putative ancestor.

Finding Direction in the Search for Selection

Tests for positive selection have mostly been developed to look for diversifying selection where change away from the current amino acid is often favorable. However, in many cases we are interested in directional selection where there is a shift toward specific amino acids, resulting in increased fitness in the species. Recently, a few methods have been developed to detect and characterize directional selection on a molecular level. Using the results of evolutionary simulations as well as HIV drug resistance data as models of directional selection, we compare two such methods with each other, as well as against a standard method for detecting diversifying selection. We find that the method to detect diversifying selection also detects directional selection under certain conditions. One method developed for detecting directional selection is powerful and accurate for a wide range of conditions, while the other can generate an excessive number of false positives.

Dynamic Convergent Evolution Drives the Passage Adaptation across 48 Years' History of H3N2 Influenza Evolution

Influenza viruses are often propagated in a diverse set of culturing media and additional substitutions known as passage adaptation can cause extra evolution in the target strain, leading to ineffective vaccines. Using 25,482 H3N2 HA1 sequences curated from Global Initiative on Sharing All Influenza Data and National Center for Biotechnology Information databases, we found that passage adaptation is a very dynamic process that changes over time and evolves in a seesaw like pattern. After crossing the species boundary from bird to human in 1968, the influenza H3N2 virus evolves to be better adapted to the human environment and passaging them in embryonated eggs (i.e., an avian environment) leads to increasingly stronger positive selection. On the contrary, passage adaptation to the mammalian cell lines changes from positive selection to negative selection. Using two statistical tests, we identified 19 codon positions around the receptor binding domain strongly contributing to passage adaptation in the embryonated egg. These sites show strong convergent evolution and overlap extensively with positively selected sites identified in humans, suggesting that passage adaptation can confound many of the earlier studies on influenza evolution. Interestingly, passage adaptation in recent years seems to target a few codon positions in antigenic surface epitopes, which makes it difficult to produce antigenically unaltered vaccines using embryonic eggs. Our study outlines another interesting scenario whereby both convergent and adaptive evolution are working in synchrony driving viral adaptation. Future studies from sequence analysis to vaccine production need to take careful consideration of passage adaptation.

Genomic Analysis of the Emergence, Evolution, and Spread of Human Respiratory RNA Viruses

The emergence and reemergence of rapidly evolving RNA viruses-particularly those responsible for respiratory diseases, such as influenza viruses and coronaviruses-pose a significant threat to global health, including the potential of major pandemics. Importantly, recent advances in high-throughput genome sequencing enable researchers to reveal the genomic diversity of these viral pathogens at much lower cost and with much greater precision than they could before. In particular, the genome sequence data generated allow inferences to be made on the molecular basis of viral emergence, evolution, and spread in human populations in real time. In this review, we introduce recent computational methods that analyze viral genomic data, particularly in combination with metadata such as sampling time, geographic location, and virulence. We then outline the insights these analyses have provided into the fundamental patterns and processes of evolution and emergence in human respiratory RNA viruses, as well as the major challenges in such genomic analyses.

A new look at the LTR retrotransposon content of the chicken genome

Background

LTR retrotransposons contribute approximately 10 % of the mammalian genome, but it has been previously reported that there is a deficit of these elements in the chicken relative to both mammals and other birds. A novel LTR retrotransposon classification pipeline, LocaTR, was developed and subsequently utilised to re-examine the chicken LTR retrotransposon annotation, and determine if the proposed chicken deficit is biologically accurate or simply a technical artefact.

Results

Using LocaTR 3.01 % of the chicken galGal4 genome assembly was annotated as LTR retrotransposon-derived elements (nearly double the previous annotation), including 1,073 that were structurally intact. Element distribution is significantly correlated with chromosome size and is non-random within each chromosome. Elements are significantly depleted within coding regions and enriched in gene sparse areas of the genome. Over 40 % of intact elements are found in clusters, unrelated by age or genera, generally in poorly recombining regions. The transcription of most LTR retrotransposons were suppressed or incomplete, but individual domain and full length retroviral transcripts were produced in some cases, although mostly with regularly interspersed stop codons in all reading frames. Furthermore, RNAseq data from 23 diverse tissues enabled greater characterisation of the co-opted endogenous retrovirus Ovex1. This gene was shown to be expressed ubiquitously but at variable levels across different tissues. LTR retrotransposon content was found to be very variable across the avian lineage and did not correlate with either genome size or phylogenetic position. However, the extent of previous, species-specific LTR retrotransposon annotation appears to be a confounding factor.

Conclusions

Use of the novel LocaTR pipeline has nearly doubled the annotated LTR retrotransposon content of the chicken genome compared to previous estimates. Further analysis has described element distribution, clustering patterns and degree of expression in a variety of adult tissues, as well as in three embryonic stages. This study also enabled better characterisation of the co-opted gamma retroviral envelope gene Ovex1. Additionally, this work suggests that there is no deficit of LTR retrotransposons within the Galliformes relative to other birds, or to mammalian genomes when scaled for the three-fold difference in genome size.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3043-1) contains supplementary material, which is available to authorized users.

Specialization of a polyphenism switch gene following serial duplications in Pristionchus nematodes

Polyphenism is an extreme manifestation of developmental plasticity, requiring distinct developmental programs and the addition of a switch mechanism. Because the genetic basis of polyphenism switches has only begun to be understood, how their mechanisms arise is unclear. In the nematode Pristionchus pacificus, which has a mouthpart polyphenism specialized for alternative diets, a gene (eud-1) executing the polyphenism switch was recently identified as the product of lineage-specific duplications. Here, we infer the role of gene duplications in producing a switch gene. Using reverse genetics and population genetic analyses, we examine evidence for competing scenarios of degeneration and complementation, neutral evolution, and functional specialization. Of the daughter genes, eud-1 alone has assumed switch-like regulation of the mouth polyphenism. Measurements of life-history traits in single, double, and triple sulfatase mutants did not, given a benign environment, identify alternative or complementary roles for eud-1 paralogs. Although possible roles are still unknown, selection analyses of the sister species and 104 natural isolates of P. pacificus detected purifying selection on the genes, suggesting their functionality by their fixation and evolutionary maintenance. Our approach shows the tractability of reverse genetics in a nontraditional model system to study evolution by gene duplication.

Evolutionary Dynamics of Floral Homeotic Transcription Factor Protein-Protein Interactions

Protein–protein interactions (PPIs) have widely acknowledged roles in the regulation of development, but few studies have addressed the timing and mechanism of shifting PPIs over evolutionary history. The B-class MADS-box transcription factors, PISTILLATA (PI) and APETALA3 (AP3) are key regulators of floral development. PI-like (PI^L) and AP3-like (AP3^L) proteins from a number of plants, including Arabidopsis thaliana (Arabidopsis) and the grass Zea mays (maize), bind DNA as obligate heterodimers. However, a PI^L protein from the grass relative Joinvillea can bind DNA as a homodimer. To ascertain whether Joinvillea PI^L homodimerization is an anomaly or indicative of broader trends, we characterized PI^L dimerization across the Poales and uncovered unexpected evolutionary lability. Both obligate B-class heterodimerization and PI^L homodimerization have evolved multiple times in the order, by distinct molecular mechanisms. For example, obligate B-class heterodimerization in maize evolved very recently from PI^L homodimerization. A single amino acid change, fixed during domestication, is sufficient to toggle one maize PI^L protein between homodimerization and obligate heterodimerization. We detected a signature of positive selection acting on residues preferentially clustered in predicted sites of contact between MADS-box monomers and dimers, and in motifs that mediate MADS PPI specificity in Arabidopsis. Changing one positively selected residue can alter PI^L dimerization activity. Furthermore, ectopic expression of a Joinvillea PI^L homodimer in Arabidopsis can homeotically transform sepals into petals. Our results provide a window into the evolutionary remodeling of PPIs, and show that novel interactions have the potential to alter plant form in a context-dependent manner.

Determination of antigenicity-altering patches on the major surface protein of human influenza A/H3N2 viruses

Human influenza viruses are rapidly evolving RNA viruses that cause short-term respiratory infections with substantial morbidity and mortality in annual epidemics. Uncovering the general principles of viral coevolution with human hosts is important for pathogen surveillance and vaccine design. Protein regions are an appropriate model for the interactions between two macromolecules, but the currently used epitope definition for the major antigen of influenza viruses, namely hemagglutinin, is very broad. Here, we combined genetic, evolutionary, antigenic, and structural information to determine the most relevant regions of the hemagglutinin of human influenza A/H3N2 viruses for interaction with human immunoglobulins. We estimated the antigenic weights of amino acid changes at individual sites from hemagglutination inhibition data using antigenic tree inference followed by spatial clustering of antigenicity-altering protein sites on the protein structure. This approach determined six relevant areas (patches) for antigenic variation that had a key role in the past antigenic evolution of the viruses. Previous transitions between successive predominating antigenic types of H3N2 viruses always included amino acid changes in either the first or second antigenic patch. Interestingly, there was only partial overlap between the antigenic patches and the patches under strong positive selection. Therefore, besides alterations of antigenicity, other interactions with the host may shape the evolution of human influenza A/H3N2 viruses.

Molecular evolution of fever, thrombocytopenia and leukocytopenia virus (FTLSV) based on whole-genome sequences

FTLSV is a novel bunyavirus that was discovered in 2007 in the Henan province of China and has reported case fatality rates of up to 30%. Despite the high case fatality rate, knowledge of the evolution and molecular epidemiology of FTLSV is limited. In this study, detailed phylogenetic analyses were performed on whole-genome sequences to examine the virus's evolutionary rates, estimate dates of common ancestry, and determine the population dynamics and selection pressure for FTLSV. The evolutionary rates of FTLSV were estimated to be 2.28×10(-4), 2.42×10(-4) and 1.19×10(-4) nucleotide substitutions/site/year for the S, M and L segments, respectively. The most recent ancestor of the viruses existed approximately 182-294 years ago. Evidence of RNA segment reassortment was found in FTLSV. A Bayesian skyline plot showed that after a period of genetic stability following high variability, the FTLSV population appeared to have contracted it. Selection pressures were estimated and revealed an abundance of negatively selected sites and sparse positively selected sites. These data will be valuable in understanding the evolution and molecular epidemiology of FTLSV, eventually helping to determine mechanisms of emergence and pathogenicity and the level of the virus's threat to public health.

Modelling the role of immunity in reversion of viral antigenic sites

Antigenic sites in viral pathogens exhibit distinctive evolutionary dynamics due to their role in evading recognition by host immunity. Antigenic selection is known to drive higher rates of non-synonymous substitution; less well understood is why differences are observed between viruses in their propensity to mutate to a novel or previously encountered amino acid. Here, we present a model to explain patterns of antigenic reversion and forward substitution in terms of the epidemiological and molecular processes of the viral population. We develop an analytical three-strain model and extend the analysis to a multi-site model to predict characteristics of observed sequence samples. Our model provides insight into how the balance between selection to escape immunity and to maintain viability is affected by the rate of mutational input. We also show that while low probabilities of reversion may be due to either a low cost of immune escape or slowly decaying host immunity, these two scenarios can be differentiated by the frequency patterns at antigenic sites. Comparison between frequency patterns of human influenza A (H3N2) and human RSV-A suggests that the increased rates of antigenic reversion in RSV-A is due to faster decaying immunity and not higher costs of escape.

Phylodynamic analysis of HIV sub-epidemics in Mochudi, Botswana

Southern Africa continues to be the epicenter of the HIV/AIDS epidemic. This HIV-1 subtype C epidemic has a predominantly heterosexual mode of virus transmission and high (>15%) HIV prevalence among adults. The epidemiological dynamics of the HIV-1C epidemic in southern Africa are still poorly understood. Here, we aim at a better understanding of HIV transmission dynamics by analyzing HIV-1 subtype C sequences from Mochudi, a peri-urban village in Botswana. HIV-1C env gene sequences (gp120 V1C5) were obtained through enhanced household-based HIV testing and counseling in Mochudi. More than 1200 sequences were generated and phylogenetically distinct sub-epidemics within Mochudi identified. The Bayesian birth-death skyline plot was used to estimate the effective reproductive number, R, and the timing of virus transmission, to classify sub-epidemics as "acute" (those with recent viral transmissions) or "historic" (those without recent viral transmissions). We identified two of the 15 sub-epidemics as "acute." The median estimates of R among the clusters ranged from 0.72 to 1.77. The majority of HIV lineages, 11 out of 15 clusters with 5+ members, appear to have been introduced to Mochudi between 1996 and 2002. The median peak duration of viral transmissions was 7.1 years (range 2.9-9.7 years). The median life span of identified HIV sub-epidemics, i.e., the time between the inferred epidemic origin and its most recent sample, was 13.1 years (range 10.2-22.1 years). Most viral transmissions within the sub-epidemics occurred between 1997 and 2007. The time period during which infected people are infectious appears to have decreased since the introduction of the national ART program in Botswana. Real-time HIV genotyping and breaking down local HIV epidemics into phylogenetically distinct sub-epidemics may help to reveal the structure and dynamics of HIV transmission networks in communities, and aid in the design of targeted interventions for members of the acute sub-epidemics that likely fuel local HIV/AIDS epidemics.

The Molecular Epidemiology and Evolution of Murray Valley Encephalitis Virus: Recent Emergence of Distinct Sub-lineages of the Dominant Genotype 1

Background

Recent increased activity of the mosquito-borne Murray Valley encephalitis virus (MVEV) in Australia has renewed concerns regarding its potential to spread and cause disease.

Methodology/Principal Findings

To better understand the genetic relationships between earlier and more recent circulating strains, patterns of virus movement, as well as the molecular basis of MVEV evolution, complete pre-membrane (prM) and Envelope (Env) genes were sequenced from sixty-six MVEV strains from different regions of the Australasian region, isolated over a sixty year period (1951–2011). Phylogenetic analyses indicated that, of the four recognized genotypes, only G1 and G2 are contemporary. G1 viruses were dominant over the sampling period and found across the known geographic range of MVEV. Two distinct sub-lineages of G1 were observed (1A and 1B). Although G1B strains have been isolated from across mainland Australia, Australian G1A strains have not been detected outside northwest Australia. Similarly, G2 is comprised of only Western Australian isolates from mosquitoes, suggesting G1B and G2 viruses have geographic or ecological restrictions. No evidence of recombination was found and a single amino acid substitution in the Env protein (S332G) was found to be under positive selection, while several others were found to be under directional evolution. Evolutionary analyses indicated that extant genotypes of MVEV began to diverge from a common ancestor approximately 200 years ago. G2 was the first genotype to diverge, followed by G3 and G4, and finally G1, from which subtypes G1A and G1B diverged between 1964 and 1994.

Conclusions/Significance

The results of this study provides new insights into the genetic diversity and evolution of MVEV. The demonstration of co-circulation of all contemporary genetic lineages of MVEV in northwestern Australia, supports the contention that this region is the enzootic focus for this virus.

Murray Valley encephalitis virus is the most significant cause of mosquito-borne encephalitis in humans in Australia, and can also cause neurological disease in horses. This study reports an expanded phylogenetic study of this virus and the first molecular evolutionary analysis. Of the four recognized genotypes of Murray Valley encephalitis virus, only two were found to be actively circulating (genotypes 1 and 2), and genotype 1 was dominant. Distinct genetic sub-lineages within genotype 1 were found to have recently emerged. Molecular clock analysis indicated that genotype 2 viruses are the oldest genetic lineage while genotype 1 viruses are the most recent to diverge. The co-circulation of distinct genetic lineages of this virus in northwestern Australia, comprising the oldest and youngest lineages, supports previous findings that MVEV circulates endemically in this region.

Widespread convergence in toxin resistance by predictable molecular evolution

The question about whether evolution is unpredictable and stochastic or intermittently constrained along predictable pathways is the subject of a fundamental debate in biology, in which understanding convergent evolution plays a central role. At the molecular level, documented examples of convergence are rare and limited to occurring within specific taxonomic groups. Here we provide evidence of constrained convergent molecular evolution across the metazoan tree of life. We show that resistance to toxic cardiac glycosides produced by plants and bufonid toads is mediated by similar molecular changes to the sodium-potassium-pump (Na(+)/K(+)-ATPase) in insects, amphibians, reptiles, and mammals. In toad-feeding reptiles, resistance is conferred by two point mutations that have evolved convergently on four occasions, whereas evidence of a molecular reversal back to the susceptible state in varanid lizards migrating to toad-free areas suggests that toxin resistance is maladaptive in the absence of selection. Importantly, resistance in all taxa is mediated by replacements of 2 of the 12 amino acids comprising the Na(+)/K(+)-ATPase H1-H2 extracellular domain that constitutes a core part of the cardiac glycoside binding site. We provide mechanistic insight into the basis of resistance by showing that these alterations perturb the interaction between the cardiac glycoside bufalin and the Na(+)/K(+)-ATPase. Thus, similar selection pressures have resulted in convergent evolution of the same molecular solution across the breadth of the animal kingdom, demonstrating how a scarcity of possible solutions to a selective challenge can lead to highly predictable evolutionary responses.

Positive Selection in CD8+ T-Cell Epitopes of Influenza Virus Nucleoprotein Revealed by a Comparative Analysis of Human and Swine Viral Lineages

Numerous experimental studies have demonstrated that CD8(+) T cells contribute to immunity against influenza by limiting viral replication. It is therefore surprising that rigorous statistical tests have failed to find evidence of positive selection in the epitopes targeted by CD8(+) T cells. Here we use a novel computational approach to test for selection in CD8(+) T-cell epitopes. We define all epitopes in the nucleoprotein (NP) and matrix protein (M1) with experimentally identified human CD8(+) T-cell responses and then compare the evolution of these epitopes in parallel lineages of human and swine influenza viruses that have been diverging since roughly 1918. We find a significant enrichment of substitutions that alter human CD8(+) T-cell epitopes in NP of human versus swine influenza virus, consistent with the idea that these epitopes are under positive selection. Furthermore, we show that epitope-altering substitutions in human influenza virus NP are enriched on the trunk versus the branches of the phylogenetic tree, indicating that viruses that acquire these mutations have a selective advantage. However, even in human influenza virus NP, sites in T-cell epitopes evolve more slowly than do nonepitope sites, presumably because these epitopes are under stronger inherent functional constraint. Overall, our work demonstrates that there is clear selection from CD8(+) T cells in human influenza virus NP and illustrates how comparative analyses of viral lineages from different hosts can identify positive selection that is otherwise obscured by strong functional constraint.

IMPORTANCE

There is a strong interest in correlates of anti-influenza immunity that are protective against diverse virus strains. CD8(+) T cells provide such broad immunity, since they target conserved viral proteins. An important question is whether T-cell immunity is sufficiently strong to drive influenza virus evolution. Although many studies have shown that T cells limit viral replication in animal models and are associated with decreased symptoms in humans, no studies have proven with statistical significance that influenza virus evolves under positive selection to escape T cells. Here we use comparisons of human and swine influenza viruses to rigorously demonstrate that human influenza virus evolves under pressure to fix mutations in the nucleoprotein that promote escape from T cells. We further show that viruses with these mutations have a selective advantage since they are preferentially located on the "trunk" of the phylogenetic tree. Overall, our results show that CD8(+) T cells targeting nucleoprotein play an important role in shaping influenza virus evolution.

Evolution and Coevolution of PRC2 Genes in Vertebrates and Mammals

Recruited by noncoding RNAs (ncRNAs) to specific genomic sites, polycomb repressive complexes 2 (PRC2) modify chromatin states in nearly all eukaryotes. The limited ncRNAs in Drosophila but abundant in mammals should have made PRC2 genes evolved significantly in Deuterostomia to adapt to the much increased ncRNAs. This study analyzes the evolution and coevolution of seven PRC2 genes in 29 Deuterostomia. These genes, previously assumed highly conserved, are found to have obtained multiple insertions in vertebrates and mammals and undergone significant positive selections in marsupials and prosimians, indicating adaptions to substantially increased lncRNAs (long noncoding RNAs) in mammals and in primates. Some insertions occur notably in homologous sequences of human nonsense-mediated decay (NMD) transcripts. Moreover, positive selections and signals of convergent evolution imply the independent increase of lncRNAs in mammals and in primates. Coevolutionary analysis reveals that patterns of interaction between PRC2 proteins have also much evolved from vertebrates to mammals, indicating adaptation at the protein complex level. The potential functions of mammalian-specific insertions and NMD transcripts deserve further experimental examination.

Molecular Epidemiology and Phylogenetic Analysis of Human Adenovirus Caused an Outbreak in Taiwan during 2011

An outbreak of adenovirus has been surveyed in Taiwan in 2011. To better understand the evolution and epidemiology of adenovirus in Taiwan, full-length sequence of hexon and fiber coapsid protein was analyzed using series of phylogenetic and dynamic evolution tools. Six different serotypes were identified in this outbreak and the species B was predominant (HAdV-3, 71.50%; HAdV-7, 15.46%). The most frequent diagnosis was acute tonsillitis (54.59%) and bronchitis (47.83%). Phylogenetic analysis revealed that hexon protein gene sequences were highly conserved for HAdV-3 and HAdV-7 circulation in Taiwan. However, comparison of restriction fragment length polymorphism (RFLP) analysis and phylogenetic trees of fiber gene in HAdV-7 clearly indicated that the predominant genotype in Taiwan has shifted from 7b to 7d. Several positive selection sites were observed in hexon protein. The estimated nucleotide substitution rates of hexon protein of HAdV-3 and HAdV-7 were 0.234×10^-3 substitutions/site/year (95% HPD: 0.387~0.095×10^-3) and 1.107×10^-3 (95% HPD: 0. 541~1.604) respectively; those of the fiber protein of HAdV-3 and HAdV-7 were 1.085×10^-3 (95% HPD: 1.767~0.486) and 0.132×10^-3 (95% HPD: 0.283~0.014) respectively. Phylodynamic analysis by Bayesian skyline plot (BSP) suggested that using individual gene to evaluate the effective population size might possibly cause miscalculation. In summary, the virus evolution is ongoing, and continuous surveillance of this virus evolution will contribute to the control of the epidemic.

Evolutionary analysis of Tembusu virus: evidence for the emergence of a dominant genotype

We recently identified Tembusu virus (TMUV) as a causative agent of duck infectious disease, which has spread in China since 2010. A recent study has indicated a potential case of human infection by TMUV, highlighting the need for further study of TMUV, especially its evolution. Because no information exists regarding the evolution of TMUV, we conducted comprehensive phylogenetic analyses using the largest collection of complete open reading frame (ORF) sequences of TMUV. Our results indicated that two lineages of viruses were associated with the 2010 outbreak in China, and lineage II, in particular sublineage II-c, has arisen as the dominant lineage currently circulating. We inferred that the most recent common ancestor (MRCA) of this TMUV was emerged around 1996. Evidence of natural recombination was also detected in TMUV. Molecular adaptation analyses revealed that strong negative selection shaped the evolution of TMUV, while a number of codons subjected to positive pressure were also identified. Our study, for the first time, illustrated the evolutionary history and character of TMUV and will be helpful for vaccine and diagnostic development.

Analysis of the evolution and structure of a complex intrahost viral population in chronic hepatitis C virus mapped by ultradeep pyrosequencing

Hepatitis C virus (HCV) causes chronic infection in up to 50% to 80% of infected individuals. Hypervariable region 1 (HVR1) variability is frequently studied to gain an insight into the mechanisms of HCV adaptation during chronic infection, but the changes to and persistence of HCV subpopulations during intrahost evolution are poorly understood. In this study, we used ultradeep pyrosequencing (UDPS) to map the viral heterogeneity of a single patient over 9.6 years of chronic HCV genotype 4a infection. Informed error correction of the raw UDPS data was performed using a temporally matched clonal data set. The resultant data set reported the detection of low-frequency recombinants throughout the study period, implying that recombination is an active mechanism through which HCV can explore novel sequence space. The data indicate that polyvirus infection of hepatocytes has occurred but that the fitness quotients of recombinant daughter virions are too low for the daughter virions to compete against the parental genomes. The subpopulations of parental genomes contributing to the recombination events highlighted a dynamic virome where subpopulations of variants are in competition. In addition, we provide direct evidence that demonstrates the growth of subdominant populations to dominance in the absence of a detectable humoral response.

IMPORTANCE

Analysis of ultradeep pyrosequencing data sets derived from virus amplicons frequently relies on software tools that are not optimized for amplicon analysis, assume random incorporation of sequencing errors, and are focused on achieving higher specificity at the expense of sensitivity. Such analysis is further complicated by the presence of hypervariable regions. In this study, we made use of a temporally matched reference sequence data set to inform error correction algorithms. Using this methodology, we were able to (i) detect multiple instances of hepatitis C virus intrasubtype recombination at the E1/E2 junction (a phenomenon rarely reported in the literature) and (ii) interrogate the longitudinal quasispecies complexity of the virome. Parallel to the UDPS, isolation of IgG-bound virions was found to coincide with the collapse of specific viral subpopulations.

Computational prediction of vaccine strains for human influenza A (H3N2) viruses

Human influenza A viruses are rapidly evolving pathogens that cause substantial morbidity and mortality in seasonal epidemics around the globe. To ensure continued protection, the strains used for the production of the seasonal influenza vaccine have to be regularly updated, which involves data collection and analysis by numerous experts worldwide. Computer-guided analysis is becoming increasingly important in this problem due to the vast amounts of generated data. We here describe a computational method for selecting a suitable strain for production of the human influenza A virus vaccine. It interprets available antigenic and genomic sequence data based on measures of antigenic novelty and rate of propagation of the viral strains throughout the population. For viral isolates sampled between 2002 and 2007, we used this method to predict the antigenic evolution of the H3N2 viruses in retrospective testing scenarios. When seasons were scored as true or false predictions, our method returned six true positives, three false negatives, eight true negatives, and one false positive, or 78% accuracy overall. In comparison to the recommendations by the WHO, we identified the correct antigenic variant once at the same time and twice one season ahead. Even though it cannot be ruled out that practical reasons such as lack of a sufficiently well-growing candidate strain may in some cases have prevented recommendation of the best-matching strain by the WHO, our computational decision procedure allows quantitative interpretation of the growing amounts of data and may help to match the vaccine better to predominating strains in seasonal influenza epidemics. Importance: Human influenza A viruses continuously change antigenically to circumvent the immune protection evoked by vaccination or previously circulating viral strains. To maintain vaccine protection and thereby reduce the mortality and morbidity caused by infections, regular updates of the vaccine strains are required. We have developed a data-driven framework for vaccine strain prediction which facilitates the computational analysis of genetic and antigenic data and does not rely on explicit evolutionary models. Our computational decision procedure generated good matches of the vaccine strain to the circulating predominant strain for most seasons and could be used to support the expert-guided prediction made by the WHO; it thus may allow an increase in vaccine efficacy.

Evolution of the hemagglutinin gene of H3N8 canine influenza virus in dogs

With the widespread use of a recently developed canine influenza virus (CIV) H3N8 vaccine, continual molecular evaluation of circulating CIVs is necessary for monitoring antigenic drift. The aim of this project was to further describe the genetic evolution of CIV, as well as determine any genetic variation within potential antigenic regions that might result in antigenic drift. To this end, the hemagglutinin gene of 19 CIV isolates from dogs residing in Colorado, New York, and South Carolina humane shelters was sequenced and compared to CIV strains isolated during 2003–2012. Phylogenetic analysis suggests that CIV might be diverging into two geographically distinct lineages. Using a mixed-effects model for evolution and single likelihood ancestor counting methods, several amino acid sites were found to be undergoing selection pressure. Additionally, a total of six amino acid changes were observed in two possible antigenic sites for CIVs isolated from Colorado and New York humane shelters between 2009 and 2011. As CIV isolates might be diverging into geographically distinct lineages, further experiments are warranted to determine the extent of antigenic drift occurring within circulating CIV.

Developmental pathway for potent V1V2-directed HIV-neutralizing antibodies

Antibodies capable of neutralizing HIV-1 often target variable regions 1 and 2 (V1V2) of the HIV-1 envelope, but the mechanism of their elicitation has been unclear. Here we define the developmental pathway by which such antibodies are generated and acquire the requisite molecular characteristics for neutralization. Twelve somatically related neutralizing antibodies (CAP256-VRC26.01-12) were isolated from donor CAP256 (from the Centre for the AIDS Programme of Research in South Africa (CAPRISA)); each antibody contained the protruding tyrosine-sulphated, anionic antigen-binding loop (complementarity-determining region (CDR) H3) characteristic of this category of antibodies. Their unmutated ancestor emerged between weeks 30-38 post-infection with a 35-residue CDR H3, and neutralized the virus that superinfected this individual 15 weeks after initial infection. Improved neutralization breadth and potency occurred by week 59 with modest affinity maturation, and was preceded by extensive diversification of the virus population. HIV-1 V1V2-directed neutralizing antibodies can thus develop relatively rapidly through initial selection of B cells with a long CDR H3, and limited subsequent somatic hypermutation. These data provide important insights relevant to HIV-1 vaccine development.