More effective purifying selection on RNA viruses than in DNA viruses

Comparative genomics and evolutionary analysis of dengue virus strains circulating in Pakistan

Dengue fever virus (DENV) poses a significant public health risk in tropical and subtropical regions across the world. Although the dengue fever virus (DENV) exhibits significant genetic diversity and has the potential to evolve, there is a lack of comprehensive research on the comparative genomics and evolutionary dynamics of the virus in Pakistan. Phylogenetic analysis demonstrated the circulation of all four dengue virus serotypes (DENV-1, - 2, - 3, and - 4) with prevalent genotypes III and V within DENV-1, cosmopolitan genotype within DENV-2, genotype III within DENV-3, and genotype I within DENV-4 during 2006-2014. Based on the complete envelope region, genome-wide residue signature and genetic diversity indicate that there is a high level of genetic diversity among DENV-1 strains, while DENV-3 strains exhibit the least genetic diversity. Comparative analysis of all four DENV serotypes revealed that certain codons in DENV-2 and -4 were subject to strong purifying selection, while a few codon sites in the envelope region showed evidence of positive selection. These findings provided valuable insights into the comparative genomics and evolutionary pattern of DENV strains reported from Pakistan. Whether those characteristics conferred a fitness advantage to DENV-1 genotypes within a specific geography and time interval warrants further investigations. The findings of the current study will contribute to tracking disease dynamics, understanding virus transmission and evolution, and formulating effective disease control strategies.

Genome-informed investigation of the molecular evolution and genetic reassortment of severe fever with thrombocytopenia syndrome virus

BACKGROUND

Severe fever with thrombocytopenia syndrome virus (SFTSV) is a viral pathogen causing significant clinical signs from mild fever with thrombocytopenia to severe hemorrhages. World Health Organization has paid special attention to the dramatic increase in human SFTS cases in China, Japan, and South Korea since the 2010s. The present study investigated the molecular evolution and genetic reassortment of SFTSVs using complete genomic sequences.

METHODS/PRINCIPAL FINDING

We collected the complete genome sequences of SFTSVs globally isolated until 2019 (L segment, n = 307; M segment, n = 326; and S segment, n = 564) and evaluated the evolutionary profiles of SFTSVs based on phylogenetic and molecular selection pressure analyses. By employing a time-scaled Bayesian inference method, we found the geographical heterogeneity of dominant SFTSV genotypes in China, Japan, and South Korea around several centuries before and locally spread by tick-born spillover with infrequent long-distance transmission. Purifying selection predominated the molecular evolution of SFTSVs with limited gene reassortment and fixed substitution, but almost all three gene segments appeared to harbor at least one amino acid residue under positive selection. Specifically, the nonstructural protein and glycoprotein (Gn/Gc) genes were preferential selective targets, and the Gn region retained the highest number of positively selected residues.

CONCLUSION/SIGNIFICANCE

Here, the large-scale genomic analyses of SFTSVs improved prior knowledge of how this virus emerged and evolved in China, Japan, and South Korea. Our results highlight the importance of SFTSV surveillance in both human and non-human reservoirs at the molecular level to fight against fatal human infection with the virus.

The Innovative Informatics Approaches of High-Throughput Technologies in Livestock: Spearheading the Sustainability and Resiliency of Agrigenomics Research

For more than a decade, next-generation sequencing (NGS) has been emerging as the mainstay of agrigenomics research. High-throughput technologies have made it feasible to facilitate research at the scale and cost required for using this data in livestock research. Scale frameworks of sequencing for agricultural and livestock improvement, management, and conservation are partly attributable to innovative informatics methodologies and advancements in sequencing practices. Genome-wide sequence-based investigations are often conducted worldwide, and several databases have been created to discover the connections between worldwide scientific accomplishments. Such studies are beginning to provide revolutionary insights into a new era of genomic prediction and selection capabilities of various domesticated livestock species. In this concise review, we provide selected examples of the current state of sequencing methods, many of which are already being used in animal genomic studies, and summarize the state of the positive attributes of genome-based research for cattle (Bos taurus), sheep (Ovis aries), pigs (Sus scrofa domesticus), horses (Equus caballus), chickens (Gallus gallus domesticus), and ducks (Anas platyrhyncos). This review also emphasizes the advantageous features of sequencing technologies in monitoring and detecting infectious zoonotic diseases. In the coming years, the continued advancement of sequencing technologies in livestock agrigenomics will significantly influence the sustained momentum toward regulatory approaches that encourage innovation to ensure continued access to a safe, abundant, and affordable food supplies for future generations.

Biochemical Characterization of Emerging SARS-CoV-2 Nsp15 Endoribonuclease Variants

Global sequencing efforts from the ongoing COVID-19 pandemic, caused by the novel coronavirus SARS-CoV-2, continue to provide insight into the evolution of the viral genome. Coronaviruses encode 16 nonstructural proteins, within the first two-thirds of their genome, that facilitate viral replication and transcription as well as evasion of the host immune response. However, many of these viral proteins remain understudied. Nsp15 is a uridine-specific endoribonuclease conserved across all coronaviruses. The nuclease activity of Nsp15 helps the virus evade triggering an innate immune response. Understanding how Nsp15 has changed over the course of the pandemic, and how mutations affect its RNA processing function, will provide insight into the evolution of an oligomerization-dependent endoribonuclease and inform drug design. In combination with previous structural data, bioinformatics analyses of 1.9 + million SARS-CoV-2 sequences revealed mutations across Nsp15's three structured domains (N-terminal, Middle, EndoU). Selected Nsp15 variants were characterized biochemically and compared to wild type Nsp15. We found that mutations to important catalytic residues decreased cleavage activity but increased the hexamer/monomer ratio of the recombinant protein. Many of the highly prevalent variants we analyzed led to decreased nuclease activity as well as an increase in the inactive, monomeric form. Overall, our work establishes how Nsp15 variants seen in patient samples affect nuclease activity and oligomerization, providing insight into the effect of these variants in vivo.

High genetic stability of co-circulating human adenovirus type 31 lineages over 59 years

Type 31 of human adenovirus species A (HAdV-A31) is a significant pathogen primarily associated with diarrhoea in children but also with life-threatening disseminated disease in allogeneic haematopoietic stem cell transplant (HSCT) recipients. Nosocomial outbreaks of HAdV-A31 have been frequently described. However, the evolution of HAdV-A31 has not been studied in detail. The evolution of other HAdV types is driven either by intertypic recombination, where different types exchange genome regions, or by immune escape selection of neutralisation determinants. Complete genomic HAdV-A31 sequences from sixty diagnostic specimens of the past 18 years (2003-21) were generated, including fourteen specimens of a presumed outbreak on two HSCT wards. Additionally, twenty-three complete genomes from GenBank were added to our phylogenetic analysis as well as in silico generated and previously published restriction fragment polymorphism (RFLP) data. Phylogenetic analysis of eighty-three genomes indicated that HAdV-A31 evolved slowly with six lineages co-circulating. The two major lineages were lineage 1, which included the prototype from 1962 and nine recent isolates, and lineage 2, which split into four sublineages and included most isolates from 2003 to 2021. The average nucleotide identity within lineages was high (99.8 per cent) and identity between lineages was 98.7 and 99.2 per cent. RFLP data allowed the construction of a lower-resolution phylogeny with two additional putative lineages. Surprisingly, regions of higher diversity separating lineages were found in gene regions coding for non-structural and minor capsid proteins. Intertypic recombinations were not observed, but the phylogeny of lineage 3 was compatible with an interlineage recombination event in the fibre gene. Applying the phylogenetic analysis to the presumed nosocomial outbreak excluded two suspected transmission events and separated it into two different, simultaneous outbreaks caused by different sublineages of lineage 2. However, due to the high nucleotide identity within HAdV-A31 lineages, the proof of infection chains remains debatable. This in-depth study on the molecular phylogeny of HAdV-A31 highlights the high genetic stability of co-circulating HAdV-A31 lineages over almost six decades. It also supports the epidemiological hypothesis that HAdV-A31 circulates as an etiological agent of a childhood disease infecting immunologically naive patients without strong positive selection of immune escape variants and recombinants.

Effects of Boswellia species on viral infections with particular attention to SARS-CoV-2

The emergence of pathogenic viruses is a worldwide frequent cause of diseases and, therefore, the design of treatments for viral infections stands as a significant research topic. Despite many efforts, the production of vaccines is faced with many obstacles and the high rate of viral resistance caused a severe reduction in the efficacy of antiviral drugs. However, the attempt of developing novel natural drugs, as well as the exertion of medicinal plants, may be an applicable solution for the treatment of viral diseases. Boswellia species exhibited a wide range of pharmacological activities in various conditions such as bronchial asthma, rheumatism, and Crohn’s illness. Additionally, pharmacological studies reported the observance of practical antiviral activities from different parts of this substance, especially the oleo-gum-resin. Therefore, this work provided an overview on the antiviral properties of Boswellia species and their potential therapeutic effects in the field of COVID-19 pandemic.

M2-2 gene as a new alternative molecular marker for phylogenetic, phylodynamic, and evolutionary studies of hRSV

The human Respiratory Syncytial Virus (hRSV) is the main causative agent of acute respiratory infections (ARI), such as pneumonia and bronchiolitis. One of the factors that lead to success in viral replication is the interaction of the M2-2 protein with the ribosomal complex. This interaction is responsible for the phase change of viral activity, acting as an inhibitor or inducer of viral replication, according to the concentration of mRNA. Based on the importance of M2-2 gene and protein have to viral physiology, we performed here evaluations of genetic diversity, phylogenetic reconstructions, phylodynamics, and selection test. Our results suggested an alternative way of classifying this virus in clades A and B, based on a new phylogenetic marker, the M2-2 gene. Therefore, our study is the first one to investigate the dynamics of the evolutionary diversification process of hRSV from the perspective of the M2-2 viral gene. In our study was also identified that the M2-2 gene is under the effect of purifying selection originated by population genetic bottlenecks. Therefore, the M2-2 gene demonstrated an interesting potential to be applied in evolutionary studies involving hRSV, recovering phylogenetic signals and traits of natural selection under the evolution of this virus.

Biochemical Characterization of Emerging SARS-CoV-2 Nsp15 Endoribonuclease Variants

Global sequencing efforts from the ongoing COVID-19 pandemic, caused by the novel coronavirus SARS-CoV-2, continue to provide insight into the evolution of the viral genome. Coronaviruses encode 16 nonstructural proteins, within the first two-thirds of their genome, that facilitate viral replication and transcription as well as evasion of the host immune response. However, many of these viral proteins remain understudied. Nsp15 is a uridine-specific endoribonuclease conserved across all coronaviruses. The nuclease activity of Nsp15 helps the virus evade triggering an innate immune response. Understanding how Nsp15 has changed over the course of the pandemic, and how mutations affect its RNA processing function, will provide insight into the evolution of an oligomerization-dependent endoribonuclease and inform drug design. In combination with previous structural data, bioinformatics analyses of 1.9+ million SARS-CoV-2 sequences revealed mutations across Nsp15’s three structured domains (N-terminal, Middle, EndoU). Selected Nsp15 variants were characterized biochemically and compared to wild type Nsp15. We found that mutations to important catalytic residues decreased cleavage activity but increased the hexamer/monomer ratio of the recombinant protein. Many of the highly prevalent variants we analyzed led to decreased nuclease activity as well as an increase in the inactive, monomeric form. Overall, our work establishes how Nsp15 variants seen in patient samples affect nuclease activity and oligomerization, providing insight into the effect of these variants in vivo .

Genetic diversity and evolution of goose astrovirus in the east of China

Goose astrovirus (GAstV), an agent of fatal visceral gout in goslings, has been widely circulating in eastern China since 2017, but little is known about its genetic diversity and systematic evolution. In this study, we isolated and sequenced nine nearly full-length GAstV genomes and conducted comprehensive genetic diversity and evolutionary analysis and compared them with other reported GAstV sequences. Our results indicated that two genotypic species of GAstV were circulating in China, and GAstV-2 subgenotype II-c had arisen as the dominant genotype in Shandong province and across the whole country. Multiple alignment of GAstV amino acid sequences revealed several characteristic mutations in GAstV-2 II-c strains, as well as additional residues in the nine new isolates which varied over time. Phylogenetic analysis of three open reading frames demonstrated different evolutionary histories. Evidence of natural recombination was also detected in GAstV, with most of the recombination occurring in the GAstV-2 II-c subgenotype. Molecular adaptation analyses revealed that the evolution of GAstV was shaped by strong negative selection, although a number of amino acids, which potentially affect host infection and cell entry, were subjected to positive pressure. Overall, these findings improve our understanding of the epidemiology and evolution of GAstV and may help in the development of vaccines and diagnostics. This article is protected by copyright. All rights reserved.

Banana bunchy top virus genetic diversity in Pakistan and association of diversity with recombination in its genomes

Banana Bunchy top virus (BBTV) is a multipartite circular single strand DNA virus that belongs to genus Babuvirus and family Nanoviridae. It causes significant crop losses worldwide and also in Pakistan. BBTV is present in Pakistan since 1988 however, till now only few (about twenty only) sequence of genomic components have been reported from the country. To have insights into current genetic diversity in Pakistan fifty-seven genomic components including five complete genomes (comprises of DNA-R, -U3, -S, -M, -C and -N components) were sequenced in this study. The genetic diversity analysis of populations from Pakistan showed that DNA-R is highly conserved followed by DNA-N, whereas DNA-U3 is highly diverse with the most diverse Common Region Stem-loop (CR-SL) in BBTV genome, a functional region, which previously been reported to have undergone recombination in Pakistani population. A Maximum Likelihood (ML) phylogenetic analysis of entire genomes of isolates by using sequence of all the components concatenated together with the reported genomes around the world revealed deeper insights about the origin of the disease in Pakistan. A comparison of the genetic diversity of Pakistani and entire BBTV populations around the world indicates that there exists a correlation between genetic diversity and recombination. Population genetics analysis indicated that the degree of selection pressure differs depending on the area and genomic component. A detailed analysis of recombination across various components and functional regions suggested that recombination is closely associated with the functional parts of BBTV genome showing high genetic diversity. Both genetic diversity and recombination analyses suggest that the CR-SL is a recombination hotspot in all BBTV genomes and among the six components DNA-U3 is the only recombined component that has extensively undergone inter and intragenomic recombination. Diversity analysis of recombinant regions results on average one and half fold increase and, in some cases up to four-fold increase due to recombination. These results suggest that recombination is significantly contributing to the genetic diversity of BBTV populations around the world.

Purifying Selection Determines the Short-Term Time Dependency of Evolutionary Rates in SARS-CoV-2 and pH1N1 Influenza

High-throughput sequencing enables rapid genome sequencing during infectious disease outbreaks and provides an opportunity to quantify the evolutionary dynamics of pathogens in near real-time. One difficulty of undertaking evolutionary analyses over short timescales is the dependency of the inferred evolutionary parameters on the timespan of observation. Crucially, there are an increasing number of molecular clock analyses using external evolutionary rate priors to infer evolutionary parameters. However, it is not clear which rate prior is appropriate for a given time window of observation due to the time-dependent nature of evolutionary rate estimates. Here, we characterize the molecular evolutionary dynamics of SARS-CoV-2 and 2009 pandemic H1N1 (pH1N1) influenza during the first 12 months of their respective pandemics. We use Bayesian phylogenetic methods to estimate the dates of emergence, evolutionary rates, and growth rates of SARS-CoV-2 and pH1N1 over time and investigate how varying sampling window and data set sizes affect the accuracy of parameter estimation. We further use a generalized McDonald-Kreitman test to estimate the number of segregating nonneutral sites over time. We find that the inferred evolutionary parameters for both pandemics are time dependent, and that the inferred rates of SARS-CoV-2 and pH1N1 decline by ∼50% and ∼100%, respectively, over the course of 1 year. After at least 4 months since the start of sequence sampling, inferred growth rates and emergence dates remain relatively stable and can be inferred reliably using a logistic growth coalescent model. We show that the time dependency of the mean substitution rate is due to elevated substitution rates at terminal branches which are 2-4 times higher than those of internal branches for both viruses. The elevated rate at terminal branches is strongly correlated with an increasing number of segregating nonneutral sites, demonstrating the role of purifying selection in generating the time dependency of evolutionary parameters during pandemics.

Chikungunya virus produced by a persistently infected mosquito cell line comprises a shorter genome and is non-infectious to mammalian cells

Although RNA viruses have high mutation rates, host cells and organisms work as selective environments, maintaining the viability of virus populations by eliminating deleterious genotypes. In serial passages of RNA viruses in a single cell line, most of these selective bottlenecks are absent, with no virus circulation and replication in different tissues or host alternation. In this work, Aedes aegypti Aag-2 cells were accidentally infected with Chikungunya virus (CHIKV) and Mayaro virus (MAYV). After numerous passages to achieve infection persistency, the infectivity of these viruses was evaluated in Ae. albopictus C6/36 cells, African green monkey Vero cells and primary-cultured human fibroblasts. While these CHIKV and MAYV isolates were still infectious to mosquito cells, they lost their ability to infect mammalian cells. After genome sequencing, it was observed that CHIKV accumulated many nonsynonymous mutations and a significant deletion in the coding sequence of the hypervariable domain in the nsP3 gene. Since MAYV showed very low titres, it was not sequenced successfully. Persistently infected Aag-2 cells also accumulated high loads of short and recombinant CHIKV RNAs, which seemed to have been originated from virus-derived DNAs. In conclusion, the genome of this CHIKV isolate could guide mutagenesis strategies for the production of attenuated or non-infectious (to mammals) CHIKV vaccine candidates. Our results also reinforce that a paradox is expected during passages of cells persistently infected by RNA viruses: more loosening for the development of more diverse virus genotypes and more pressure for virus specialization to this constant cellular environment.

High Mutation Frequency and Significant Population Differentiation in Papaya Ringspot Virus-W Isolates

A total of 101 papaya ringspot virus-W (PRSV-W) isolates were collected from five different cucurbit hosts in six counties of Oklahoma during the 2016–2018 growing seasons. The coat protein (CP) coding region of these isolates was amplified by reverse transcription-polymerase chain reaction, and 370 clones (3–5 clones/isolate) were sequenced. Phylogenetic analysis revealed three phylogroups while host, location, and collection time of isolates had minimal impact on grouping pattern. When CP gene sequences of these isolates were compared with sequences of published PRSV isolates (both P and W strains), they clustered into four phylogroups based on geographical location. Oklahoman PRSV-W isolates formed one of the four distinct major phylogroups. The permutation-based tests, including Ks, Ks *, Z *, Snn, and neutrality tests, indicated significant genetic differentiation and polymorphisms among PRSV-W populations in Oklahoma. The selection analysis confirmed that the CP gene is undergoing purifying selection. The mutation frequencies among all PRSV-W isolates were within the range of 1 × 10⁻³. The substitution mutations in 370 clones of PRSV-W isolates showed a high proportion of transition mutations, which gave rise to higher GC content. The N-terminal region of the CP gene mostly contained the variable sites with numerous mutational hotspots, while the core region was highly conserved.

Molecular Analysis of Full-Length VP2 of Canine Parvovirus Reveals Antigenic Drift in CPV-2b and CPV-2c Variants in Central Chile

Canine parvovirus (CPV) is a major pathogen in canines, with a high mortality rate in unvaccinated puppies. CPV is traditionally classified into three antigenic variants (CPV-2a, CPV-2b and CPV-2c) based on the amino acid sequence of the VP2 protein. Currently, various mutations are described in the receptor-binding area or in the regions of greatest antigenicity of the VP2 protein, giving rise to new viral variants that are capable of immunological escape, affecting the protective immunity of traditional vaccines. In the present study, a molecular characterization of the VP2 gene was performed, which included phylogenetic analysis, amino acid characterization and determination of selection pressures. Blood samples were initially collected from canine patients with clinical signs of gastrointestinal infection, of which 69 were positive for CPV as measured by means of PCR and 18 samples were selected for the amplification of the complete VP2 gene. The analysis revealed a higher rate of CPV-2c-positive patients compared to CPV-2b. Furthermore, the amino acid characterization of VP2 indicated mutations in the regions of highest antigenicity previously described in the literature (CPV-2b: 297 and 324; CPV-2c: 440), as well as others not previously documented (CPV-2b: 514; CPV-2c: 188, 322, 379, 427 and 463). Our analysis of selection pressure showed that the VP2 gene is under negative selection. However, positive selection point sites were identified, both in CPV-2c (324, 426 and 440) and CPV-2b (297 and 324), at sites that have been associated with evasion of the immune response via antigenic drift, which possibly has implications for the protective immunity generated by traditional vaccines.

The interplay of SARS-CoV-2 evolution and constraints imposed by the structure and functionality of its proteins

The unprecedented pace of the sequencing of the SARS-CoV-2 virus genomes provides us with unique information about the genetic changes in a single pathogen during ongoing pandemic. By the analysis of close to 200,000 genomes we show that the patterns of the SARS-CoV-2 virus mutations along its genome are closely correlated with the structural and functional features of the encoded proteins. Requirements of foldability of proteins' 3D structures and the conservation of their key functional regions, such as protein-protein interaction interfaces, are the dominant factors driving evolutionary selection in protein-coding genes. At the same time, avoidance of the host immunity leads to the abundance of mutations in other regions, resulting in high variability of the missense mutation rate along the genome. "Unexplained" peaks and valleys in the mutation rate provide hints on function for yet uncharacterized genomic regions and specific protein structural and functional features they code for. Some of these observations have immediate practical implications for the selection of target regions for PCR-based COVID-19 tests and for evaluating the risk of mutations in epitopes targeted by specific antibodies and vaccine design strategies.

Genomic and immunogenic changes of Piscine novirhabdovirus (Viral Hemorrhagic Septicemia Virus) over its evolutionary history in the Laurentian Great Lakes

A unique and highly virulent subgenogroup (-IVb) of Piscine novirhabdovirus, also known as Viral Hemorrhagic Septicemia Virus (VHSV), suddenly appeared in the Laurentian Great Lakes, causing large mortality outbreaks in 2005 and 2006, and affecting >32 freshwater fish species. Periods of apparent dormancy have punctuated smaller and more geographically-restricted outbreaks in 2007, 2008, and 2017. In this study, we conduct the largest whole genome sequencing analysis of VHSV-IVb to date, evaluating its evolutionary changes from 48 isolates in relation to immunogenicity in cell culture. Our investigation compares genomic and genetic variation, selection, and rates of sequence changes in VHSV-IVb, in relation to other VHSV genogroups (VHSV-I, VHSV-II, VHSV-III, and VHSV-IVa) and with other Novirhabdoviruses. Results show that the VHSV-IVb isolates we sequenced contain 253 SNPs (2.3% of the total 11,158 nucleotides) across their entire genomes, with 85 (33.6%) of them being non-synonymous. The most substitutions occurred in the non-coding region (NCDS; 4.3%), followed by the Nv- (3.8%), and M- (2.8%) genes. Proportionally more M-gene substitutions encoded amino acid changes (52.9%), followed by the Nv- (50.0%), G- (48.6%), N- (35.7%) and L- (23.1%) genes. Among VHSV genogroups and subgenogroups, VHSV-IVa from the northeastern Pacific Ocean has shown the fastest substitution rate (2.01x10^-3), followed by VHSV-IVb (6.64x10^-5) and by the VHSV-I, -II and-III genogroups from Europe (4.09x10^-5). A 2016 gizzard shad (Dorosoma cepedianum) from Lake Erie possessed the most divergent VHSV-IVb sequence. The in vitro immunogenicity analysis of that sample displayed reduced virulence (as did the other samples from 2016), in comparison to the original VHSV-IVb isolate (which had been traced back to 2003, as an origin date). The 2016 isolates that we tested induced milder impacts on fish host cell innate antiviral responses, suggesting altered phenotypic effects. In conclusion, our overall findings indicate that VHSV-IVb has undergone continued sequence change and a trend to lower virulence over its evolutionary history (2003 through present-day), which may facilitate its long-term persistence in fish host populations.

Mutational pressure by host APOBEC3s more strongly affects genes expressed early in the lytic phase of herpes simplex virus-1 (HSV-1) and human polyomavirus (HPyV) infection

Herpes-Simplex Virus 1 (HSV-1) infects most humans when they are young, sometimes with fatal consequences. Gene expression occurs in a temporal order upon lytic HSV-1 infection: immediate early (IE) genes are expressed, then early (E) genes, followed by late (L) genes. During this infection cycle, the HSV-1 genome has the potential for exposure to APOBEC3 (A3) proteins, a family of cytidine deaminases that cause C>U mutations on single-stranded DNA (ssDNA), often resulting in a C>T transition. We developed a computational model for the mutational pressure of A3 on the lytic cycle of HSV-1 to determine which viral kinetic gene class is most vulnerable to A3 mutations. Using in silico stochastic methods, we simulated the infectious cycle under varying intensities of A3 mutational pressure. We found that the IE and E genes are more vulnerable to A3 than L genes. We validated this model by analyzing the A3 evolutionary footprints in 25 HSV-1 isolates. We find that IE and E genes have evolved to underrepresent A3 hotspot motifs more so than L genes, consistent with greater selection pressure on IE and E genes. We extend this model to two-step infections, such as those of polyomavirus, and find that the same pattern holds for over 25 human Polyomavirus (HPyVs) genomes. Genes expressed earlier during infection are more vulnerable to mutations than those expressed later.

Variable routes to genomic and host adaptation among coronaviruses

Natural selection operating on the genomes of viral pathogens in different host species strongly contributes to adaptation facilitating host colonization. Here, we analyse, quantify and compare viral adaptation in genomic sequence data derived from seven zoonotic events in the Coronaviridae family among primary, intermediate and human hosts. Rates of nonsynonymous (d_N) and synonymous (d_S) changes on specific amino acid positions were quantified for each open reading frame (ORF). Purifying selection accounted for 77% of all sites under selection. Diversifying selection was most frequently observed in viruses infecting the primary hosts of each virus and predominantly occurred in the orf1ab genomic region. Within all four intermediate hosts, diversifying selection on the spike gene was observed either solitarily or in combination with orf1ab and other genes. Consistent with previous evidence, pervasive diversifying selection on coronavirus spike genes corroborates the role this protein plays in host cellular entry, adaptation to new hosts and evasion of host cellular immune responses. Structural modelling of spike proteins identified a significantly higher proportion of sites for SARS‐CoV‐2 under positive selection in close proximity to sites of glycosylation relative to the other coronaviruses. Among human coronaviruses, there was a significant inverse correlation between the number of sites under positive selection and the estimated years since the virus was introduced into the human population. Abundant diversifying selection observed in SARS‐CoV‐2 suggests the virus remains in the adaptive phase of the host switch, typical of recent host switches. A mechanistic understanding of where, when and how genomic adaptation occurs in coronaviruses following a host shift is crucial for vaccine design, public health responses and predicting future pandemics.

Natural Selection Plays an Important Role in Shaping the Codon Usage of Structural Genes of the Viruses Belonging to the Coronaviridae Family

Viruses belonging to the Coronaviridae family have a single-stranded positive-sense RNA with a poly-A tail. The genome has a length of ~29.9 kbps, which encodes for genes that are essential for cell survival and replication. Different evolutionary constraints constantly influence the codon usage bias (CUB) of different genes. A virus optimizes its codon usage to fit the host environment on which it savors. This study is a comprehensive analysis of the CUB for the different genes encoded by viruses of the Coronaviridae family. Different methods including relative synonymous codon usage (RSCU), an Effective number of codons (ENc), parity plot 2, and Neutrality plot, were adopted to analyze the factors responsible for the genetic evolution of the Coronaviridae family. Base composition and RSCU analyses demonstrated the presence of A-ended and U-ended codons being preferred in the 3rd codon position and are suggestive of mutational selection. The lesser ENc value for the spike ‘S’ gene suggests a higher bias in the codon usage of this gene compared to the other structural genes. Parity plot 2 and neutrality plot analyses demonstrate the role and the extent of mutational and natural selection towards the codon usage pattern. It was observed that the structural genes of the Coronaviridae family analyzed in this study were at the least under 84% influence of natural selection, implying a major role of natural selection in shaping the codon usage.

SARS-CoV-2, an evolutionary perspective of interaction with human ACE2 reveals undiscovered amino acids necessary for complex stability

The emergence of SARS-CoV-2 has resulted in nearly 1,280,000 infections and 73,000 deaths globally so far. This novel virus acquired the ability to infect human cells using the SARS-CoV cell receptor hACE2. Because of this, it is essential to improve our understanding of the evolutionary dynamics surrounding the SARS-CoV-2 hACE2 interaction. One way theory predicts selection pressures should shape viral evolution is to enhance binding with host cells. We first assessed evolutionary dynamics in select betacoronavirus spike protein genes to predict whether these genomic regions are under directional or purifying selection between divergent viral lineages, at various scales of relatedness. With this analysis, we determine a region inside the receptor-binding domain with putative sites under positive selection interspersed among highly conserved sites, which are implicated in structural stability of the viral spike protein and its union with human receptor ACE2. Next, to gain further insights into factors associated with recognition of the human host receptor, we performed modeling studies of five different betacoronaviruses and their potential binding to hACE2. Modeling results indicate that interfering with the salt bridges at hot spot 353 could be an effective strategy for inhibiting binding, and hence for the prevention of SARS-CoV-2 infections. We also propose that a glycine residue at the receptor-binding domain of the spike glycoprotein can have a critical role in permitting bat SARS-related coronaviruses to infect human cells.

Evolutionary trajectory of fish Piscine novirhabdovirus (=Viral Hemorrhagic Septicemia Virus) across its Laurentian Great Lakes history: Spatial and temporal diversification

Piscine novirhabdovirus = Viral Hemorrhagic Septicemia Virus (VHSV) first appeared in the Laurentian Great Lakes with large outbreaks from 2005 to 2006, as a new and novel RNA rhabdovirus subgenogroup (IVb) that killed >30 fish species. Interlude periods punctuated smaller more localized outbreaks in 2007, 2010, and 2017, although some fishes tested positive in the intervals. There have not been reports of outbreaks or positives from 2018, 2019, or 2020. Here, we employ a combined population genetics and phylogenetic approach to evaluate spatial and temporal evolutionary trajectory on its G-gene sequence variation, in comparison with whole-genome sequences (11,083 bp) from a subset of 44 individual isolates (including 40 newly sequenced ones). Our results show that IVb (N = 184 individual fish isolates) diversified into 36 G-gene haplotypes from 2003 to 2017, stemming from two originals ("a" and "b"). G-gene haplotypes "a" and "b" differed by just one synonymous single-nucleotide polymorphism (SNP) substitution, remained the most abundant until 2011, then disappeared. Group "a" descendants (14 haplotypes) remained most prevalent in the Upper and Central Great Lakes, with eight (51%) having nonsynonymous substitutions. Group "b" descendants primarily have occurred in the Lower Great Lakes, including 22 haplotypes, of which 15 (68%) contained nonsynonymous changes. Evolutionary patterns of the whole-genome sequences (which had 34 haplotypes among 44 isolates) appear congruent with those from the G-gene. Virus populations significantly diverged among the Upper, Central, and Lower Great Lakes, diversifying over time. Spatial divergence was apparent in the overall patterns of nucleotide substitutions, while amino acid changes increased temporally. VHSV-IVb thus significantly differentiated across its less than two decades in the Great Lakes, accompanied by declining outbreaks and virulence. Continuing diversification likely allowed the virus to persist at low levels in resident fish populations, and may facilitate its potential for further and future spread to new habitats and nonacclimated hosts.

The interplay of SARS-CoV-2 evolution and constraints imposed by the structure and functionality of its proteins

Fast evolution of the SARS-CoV-2 virus provides us with unique information about the patterns of genetic changes in a single pathogen in the timescale of months. This data is used extensively to track the phylodynamic of the pandemic's spread and its split into distinct clades. Here we show that the patterns of SARS-CoV-2 virus mutations along its genome are closely correlated with the structural features of the coded proteins. We show that the foldability of proteins' 3D structures and conservation of their functions are the universal factors driving evolutionary selection in protein-coding genes. Insights from the analysis of mutation distribution in the context of the SARS-CoV-2 proteins' structures and functions have practical implications including evaluating potential antigen epitopes or selection of primers for PCR-based COVID-19 tests.

Many human RNA viruses show extraordinarily stringent selective constraints on protein evolution

How negative selection, positive selection, and population size contribute to the large variation in nucleotide substitution rates among RNA viruses remains unclear. Here, we studied the ratios of nonsynonymous-to-synonymous substitution rates (d _N/d _S) in protein-coding genes of human RNA and DNA viruses and mammals. Among the 21 RNA viruses studied, 18 showed a genome-average d _N/d _S from 0.01 to 0.10, indicating that over 90% of nonsynonymous mutations are eliminated by negative selection. Only HIV-1 showed a d _N/d _S (0.31) higher than that (0.22) in mammalian genes. By comparing the d _N/d _S values among genes in the same genome and among species or strains, we found that both positive selection and population size play significant roles in the d _N/d _S variation among genes and species. Indeed, even in flaviviruses and picornaviruses, which showed the lowest ratios among the 21 species studied, positive selection appears to have contributed significantly to d _N/d _S We found the view that positive selection occurs much more frequently in influenza A subtype H3N2 than subtype H1N1 holds only for the hemagglutinin and neuraminidase genes, but not for other genes. Moreover, we found no support for the view that vector-borne RNA viruses have lower d _N/d _S ratios than non-vector-borne viruses. In addition, we found a correlation between d _N and d _S, implying a correlation between d _N and the mutation rate. Interestingly, only 2 of the 8 DNA viruses studied showed a d _N/d _S < 0.10, while 4 showed a d _N/d _S > 0.22. These observations increase our understanding of the mechanisms of RNA virus evolution.

Human metapneumovirus prevalence and patterns of subgroup persistence identified through surveillance of pediatric pneumonia hospital admissions in coastal Kenya, 2007-2016

BACKGROUND

Human metapneumovirus (HMPV) is an important respiratory pathogen that causes seasonal epidemics of acute respiratory illness and contributes significantly to childhood pneumonia. Current knowledge and understanding on its patterns of spread, prevalence and persistence in communities in low resource settings is limited.

METHODS

We present findings of a molecular-epidemiological analysis of nasal samples from children < 5 years of age admitted with syndromic pneumonia between 2007 and 2016 to Kilifi County Hospital, coastal Kenya. HMPV infection was detected using real-time RT-PCR and positives sequenced in the fusion (F) and attachment (G) genes followed by phylogenetic analysis. The association between disease severity and HMPV subgroup was assessed using Fisher's exact test.

RESULTS

Over 10 years, 274/6756 (4.1%) samples screened were HMPV positive. Annual prevalence fluctuated between years ranging 1.2 to 8.7% and lowest in the recent years (2014-2016). HMPV detections were most frequent between October of one year to April of the following year. Genotyping was successful for 205/274 (74.8%) positives revealing clades A2b (41.0%) and A2c (10.7%), and subgroups B1 (23.4%) and B2 (24.9%). The dominance patterns were: clade A2b between 2007 and 11, subgroup B1 between 2012 and 14, and clade A2c in more recent epidemics. Subgroup B2 viruses were present in all the years. Temporal phylogenetic clustering within the subgroups for both local and global sequence data was seen. Subgroups occurring in each epidemic season were comprised of multiple variants. Pneumonia severity did not vary by subgroup (p = 0.264). In both the F and G gene, the sequenced regions were found to be predominantly under purifying selection.

CONCLUSION

Subgroup patterns from this rural African setting temporally map with global strain distribution, suggesting a well-mixed global virus transmission pool of HMPV. Persistence in the local community is characterized by repeated introductions of HMPV variants from the global pool. The factors underlying the declining prevalence of HMPV in this population should be investigated.

Cosmopolitan A1 lineage of dengue virus serotype 2 is circulating in Pakistan: A study from 2017 dengue viral outbreak

Dengue viral infection has become a challenge in tropical and subtropical countries where dengue virus is endemic. Its epidemics are occurring at higher rates amid its circulation throughout the year. Since the first documented outbreak in Pakistan in 1994, this region has reported many sporadic cases and epidemics. There is availability of small scale demographic and epidemiological studies on dengue viral infection in Pakistan. The year 2017 witnessed a huge dengue outbreak in Peshawar city of Pakistan with 69 deaths and 24 807 laboratory-confirmed cases. We suspect that the circulation of a different lineage or genotype could be responsible for the enhanced number of infected patients in Pakistan's 2017 outbreak since previous studies have already described this phenomenon in other countries. For this, we collected 1447 suspected blood samples and their epidemiological data. After serotyping through polymerase chain reaction nine samples of Dengue virus2 (DENV2) were randomly selected and were subjected to Sanger's sequencing for genotyping analysis. The mean distance, genetic diversity, and phylogenetic analysis were carried out using K2 model. The phylogenetic analysis split Pakistani isolates into two lineages, the sequences from 2017 outbreak in Peshawar grouped within A1 lineage of cosmopolitan genotype (IV) of DENV2. The difference in distance, genetic diversity, and amino acids composition strongly back the results that the new lineage is circulating in the region. This is significant as Pakistan is struggling to control dengue epidemics which have caused much loss in both monetary and health sectors.

Assessing the Diversity of Endogenous Viruses Throughout Ant Genomes

Endogenous viral elements (EVEs) can play a significant role in the evolution of their hosts and have been identified in animals, plants, and fungi. Additionally, EVEs potentially provide an important snapshot of the evolutionary frequency of viral infection. The purpose of this study is to take a comparative host-centered approach to EVE discovery in ant genomes to better understand the relationship of EVEs to their ant hosts. Using a comprehensive bioinformatic pipeline, we screened all nineteen published ant genomes for EVEs. Once the EVEs were identified, we assessed their phylogenetic relationships to other closely related exogenous viruses. A diverse group of EVEs were discovered in all screened ant host genomes and in many cases are similar to previously identified exogenous viruses. EVEs similar to ssRNA viral proteins are the most common viral lineage throughout the ant hosts, which is potentially due to more chronic infection or more effective endogenization of certain ssRNA viruses in ants. In addition, both EVEs similar to viral glycoproteins and retrovirus-derived proteins are also abundant throughout ant genomes, suggesting their tendency to endogenize. Several of these newly discovered EVEs are found to be potentially functional within the genome. The discovery and analysis of EVEs is essential in beginning to understand viral–ant interactions over evolutionary time.

Effects of life history and ecology on virus evolutionary potential

The life history traits of viruses pose many consequences for viral population structure. In turn, population structure may influence the evolutionary trajectory of a virus. Here we review factors that affect the evolutionary potential of viruses, including rates of mutation and recombination, bottlenecks, selection pressure, and ecological factors such as the requirement for hosts and vectors. Mutation, while supplying a pool of raw genetic material, also results in the generation of numerous unfit mutants. The infection of multiple host species may expand a virus' ecological niche, although it may come at a cost to genetic diversity. Vector-borne viruses often experience a diminished frequency of positive selection and exhibit little diversity, and resistance against vector-borne viruses may thus be more durable than against non-vectored viruses. Evidence indicates that adaptation to a vector is more evolutionarily difficult than adaptation to a host. Overall, a better understanding of how various factors influence viral dynamics in both plant and animal pathosystems will lead to more effective anti-viral treatments and countermeasures.

Dynamic interactions between RNA viruses and human hosts unravelled by a decade of next generation sequencing

BACKGROUND

Next generation sequencing (NGS) methods have significantly contributed to a paradigm shift in genomic research for nearly a decade now. These methods have been useful in studying the dynamic interactions between RNA viruses and human hosts.

SCOPE OF THE REVIEW

In this review, we summarise and discuss key applications of NGS in studying the host - pathogen interactions in RNA viral infections of humans with examples.

MAJOR CONCLUSIONS

Use of NGS to study globally relevant RNA viral infections have revolutionized our understanding of the within host and between host evolution of these viruses. These methods have also been useful in clinical decision-making and in guiding biomedical research on vaccine design.

GENERAL SIGNIFICANCE

NGS has been instrumental in viral genomic studies in resolving within-host viral genomic variants and the distribution of nucleotide polymorphisms along the full-length of viral genomes in a high throughput, cost effective manner. In the future, novel advances such as long read, single molecule sequencing of viral genomes and simultaneous sequencing of host and pathogens may become the standard of practice in research and clinical settings. This will also bring on new challenges in big data analysis.

Genetic Diversity and the Spatio-Temporal Analyses of Hantaviruses in Shandong Province, China

Hemorrhagic fever with renal syndrome (HFRS) is a serious public health problem in Shandong Province, China. We conducted an epizootiologic investigation and phylogeographic and phylodynamic analyses to infer the phylogenetic relationships of hantaviruses in space and time, and gain further insights into their evolutionary dynamics in Shandong Province. Our data indicated that the Seoul virus (SEOV) is distributed throughout Shandong, whereas Hantaan virus (HTNV) co-circulates with SEOV in the eastern and southern areas of Shandong. Their distribution showed strong geographic clustering. In addition, our analyses indicated multiple evolutionary paths, long-distance transmission, and demographic expansion events for SEOV in some areas. Selection pressure analyses revealed that negative selection on hantaviruses acted as the principal evolutionary force, whereas a little evidence of positive selection exists. We found that several positively selected sites were located within major functional regions and indicated the importance of these residues for adaptive evolution of hantaviruses.

The diversity, evolution and epidemiology of plant viruses: A phylogenetic view

During the past four decades, the scientific community has seen an exponential advance in the number, sophistication, and quality of molecular techniques and bioinformatics tools for the genetic characterization of plant virus populations. Predating these advances, the field of Phylogenetics has significantly contributed to understand important aspects of plant virus evolution. This review aims at summarizing the impact of Phylogenetics in the current knowledge on three major aspects of plant virus evolution that have benefited from the development of phylogenetic inference: (1) The identification and classification of plant virus diversity. (2) The mechanisms and forces shaping the evolution of plant virus populations. (3) The understanding of the interaction between plant virus evolution, epidemiology and ecology. The work discussed here highlights the important role of phylogenetic approaches in the study of the dynamics of plant virus populations.

Molecular evolution of porcine epidemic diarrhea virus and porcine deltacoronavirus strains in Central China

Porcine epidemic diarrhea virus (PEDV) and porcine deltacoronavirus (PDCoV) are epizootic swine viruses. To detect and study the evolution of PEDV and PDCoV in central China (Shanxi, Henan, Hubei province), 70 clinical intestinal and fecal samples from piglets with severe watery diarrhea during August 2015 and June 2016 were collected, tested and analyzed. PEDV was more frequently detected by PCR than PDCoV. Phylogenetic analysis of S genes showed that the 10 PEDV strains from this study clustered into G2a (n = 7) and G2b (n = 3) groups. Additionally, the three G2b strains (PEDV S2△) contained the same specific 3 nt deletion in S2 as other reference strains in G2b. Interestingly, complete genome analysis indicated that CH/hubei/2016 was closer to the US INDEL strain and G2a group. CH/hubei/2016 had one recombination event in S2 gene which may have resulted from AH2012-12 (from G2b group) and CH-ZMDZY-11 (from G2a group). Furthermore, 10 purifying selection sites in S gene indicated an adaptive evolution of PEDV in central China swine herds. These results suggested that Pandemic G2a and G2b are predominant PEDV genotype circulating in central China. In addition, the deletion and recombination identified in S gene suggested PEDV strains of central exhibited an evolutionary variety. However, whether these changes affect the pathogenicity and antigenicity of wild PEDV is unknown and is worth for further investigation.

•

PEDV (84.2%) infection could be more commonly detected than PDCoV (2.9%) in central China.

•

A specific 3 nt-deletion in S2 gene was firstly reported in PEDV strains of central China.

•

The further analyses provided evidence of the relationship between PEDV S2△ and previous PEDV stains (3-deletion in S2).

Longitudinal Surveillance of Porcine Rotavirus B Strains from the United States and Canada and In Silico Identification of Antigenically Important Sites

Rotavirus B (RVB) is an important swine pathogen, but control and prevention strategies are limited without an available vaccine. To develop a subunit RVB vaccine with maximal effect, we characterized the amino acid sequence variability and predicted antigenicity of RVB viral protein 7 (VP7), a major neutralizing antibody target, from clinically infected pigs in the United States and Canada. We identified genotype-specific antigenic sites that may be antibody neutralization targets. While some antigenic sites had high amino acid functional group diversity, nine antigenic sites were completely conserved. Analysis of nucleotide substitution rates at amino acid sites (dN/dS) suggested that negative selection appeared to be playing a larger role in the evolution of the identified antigenic sites when compared to positive selection, and was identified in six of the nine conserved antigenic sites. These results identified important characteristics of RVB VP7 variability and evolution and suggest antigenic residues on RVB VP7 that are negatively selected and highly conserved may be good candidate regions to include in a subunit vaccine design due to their tendency to remain stable.

Baculovirus Molecular Evolution via Gene Turnover and Recurrent Positive Selection of Key Genes

Hosts and viruses are locked in an evolutionary arms race. Hosts are constantly evolving to suppress virulence and replication, while viruses, which are reliant on host machinery for survival and reproduction, develop counterstrategies to escape this immune defense. Viruses must also adapt to novel conditions while establishing themselves in a host species. Both processes provide strong selection for viral adaptation. Understanding adaptive evolution in insect viruses can help us to better understand adaptive evolution in general and is important due to the use of these viruses as biocontrol agents and for protecting ecologically or economically important species from outbreaks. Here we examine the molecular evolution of baculoviruses and nudiviruses, a group of insect-infecting viruses with key roles in biocontrol. We looked for signatures of selection between genomes of baculoviruses infecting a range of species and within a population of baculoviruses. Both analyses found only a few strong signatures of positive selection, primarily in replication- and transcription-associated genes and several structural protein genes. In both analyses, we detected a conserved complex of genes, including the helicase gene, showing consistently high levels of adaptive evolution, suggesting that they may be key in antagonistic coevolution to escape host suppression. These genes are integral to the baculovirus life cycle and may be good focal genes for developing baculoviruses as effective biocontrol agents or for targeting baculoviruses infecting ecologically relevant species. Recombination and complex genomes make evolution in these double-stranded DNA viruses more efficient than that in smaller RNA viruses with error-prone replication, as seen via signatures of selection in specific genes within a population of baculoviruses.

IMPORTANCE Most viral evolutionary studies focus on RNA viruses. While these viruses cause many human and animal diseases, such studies leave us with a lesser understanding of how DNA viruses adapt to hosts and how the host responds to these pathogens. In this paper, we focus on the evolution of baculoviruses, a group of insect-infecting DNA viruses, many of which have been used in biocontrol. We find that most of the genome is under purifying selection, with only a few key genes evolving adaptively. Our results provide a glimpse into how DNA viruses differ from RNA viruses in their evolutionary dynamics and identify genes that are key to DNA virus adaptation, improving our understanding of how this group of pathogens evolves.

Constraints on the Genetic and Antigenic Variability of Measles Virus

Antigenic drift and genetic variation are significantly constrained in measles virus (MeV). Genetic stability of MeV is exceptionally high, both in the lab and in the field, and few regions of the genome allow for rapid genetic change. The regions of the genome that are more tolerant of mutations (i.e., the untranslated regions and certain domains within the N, C, V, P, and M proteins) indicate genetic plasticity or structural flexibility in the encoded proteins. Our analysis reveals that strong constraints in the envelope proteins (F and H) allow for a single serotype despite known antigenic differences among its 24 genotypes. This review describes some of the many variables that limit the evolutionary rate of MeV. The high genomic stability of MeV appears to be a shared property of the Paramyxovirinae, suggesting a common mechanism that biologically restricts the rate of mutation.

Genome-Wide Analysis of Human Metapneumovirus Evolution

Human metapneumovirus (HMPV) has been described as an important etiologic agent of upper and lower respiratory tract infections, especially in young children and the elderly. Most of school-aged children might be introduced to HMPVs, and exacerbation with other viral or bacterial super-infection is common. However, our understanding of the molecular evolution of HMPVs remains limited. To address the comprehensive evolutionary dynamics of HMPVs, we report a genome-wide analysis of the eight genes (N, P, M, F, M2, SH, G, and L) using 103 complete genome sequences. Phylogenetic reconstruction revealed that the eight genes from one HMPV strain grouped into the same genetic group among the five distinct lineages (A1, A2a, A2b, B1, and B2). A few exceptions of phylogenetic incongruence might suggest past recombination events, and we detected possible recombination breakpoints in the F, SH, and G coding regions. The five genetic lineages of HMPVs shared quite remote common ancestors ranging more than 220 to 470 years of age with the most recent origins for the A2b sublineage. Purifying selection was common, but most protein genes except the F and M2-2 coding regions also appeared to experience episodic diversifying selection. Taken together, these suggest that the five lineages of HMPVs maintain their individual evolutionary dynamics and that recombination and selection forces might work on shaping the genetic diversity of HMPVs.

Genetic Diversity and Positive Selection Analysis of Classical Swine Fever Virus Envelope Protein Gene E2 in East China under C-Strain Vaccination

Classical swine fever virus (CSFV) causes an economically important and highly contagious disease of pigs worldwide. C-strain vaccination is one of the most effective ways to contain this disease. Since 2014, sporadic CSF outbreaks have been occurring in some C-strain vaccinated provinces of China. To decipher the disease etiology, 25 CSFV E2 genes from 169 clinical samples were cloned and sequenced. Phylogenetic analyses revealed that all 25 isolates belonged to subgenotype 2.1. Twenty-three of the 25 isolates were clustered in a newly defined subgenotype, 2.1d, and shared some consistent molecular characteristics. To determine whether the complete E2 gene was under positive selection pressure, we used a site-by-site analysis to identify specific codons that underwent evolutionary selection, and seven positively selected codons were found. Three positively selected sites (amino acids 17, 34, and 72) were identified in antigenicity-relevant domains B/C of the amino-terminal half of the E2 protein. In addition, another positively selected site (amino acid 200) exhibited a polarity change from hydrophilic to hydrophobic, which may change the antigenicity and virulence of CSFV. The results indicate that the circulating CSFV strains in Shandong province were mostly clustered in subgenotype 2.1d. Moreover, the identification of these positively selected sites could help to reveal molecular determinants of virulence or pathogenesis, and to clarify the driving force of CSFV evolution in East China.

Analysis of Soybean mosaic virus genetic diversity in Iran allows the characterization of a new mutation resulting in overcoming Rsv4-resistance

The genetic variation and population structure of Soybean mosaic virus (SMV) in Iran was analysed through the characterization of a set of isolates collected in the soybean-growing provinces of Iran. The partial nucleotide sequence of these isolates showed a single, undifferentiated population with low genetic diversity, highly differentiated from other SMV world populations. These traits are compatible with a population bottleneck associated with the recent introduction of SMV in Iran. Phylogenetic analyses suggest that SMV was introduced into Iran from East Asia, with at least three introduction events. The limited genetic diversification of SMV in Iran may be explained by strong negative selection in most viral genes eliminating the majority of mutations, together with recombination purging deleterious mutations. The pathogenicity of Iranian SMV isolates was typified on a set of soybean differential lines either susceptible or carrying different resistance genes or alleles to SMV. Two pathotypes were distinguished according to the ability to overcome Rsv4 resistance in line V94-5152. Amino acid sequence comparisons of virulent and avirulent isolates on V94-5152 (Rsv4), plus site-directed mutagenesis in a biologically active cDNA clone, identified mutation S1053N in the P3 protein as the determinant for virulence on V94-5152. Codon 1053 was shown to be under positive selection, and S1053N-determined Rsv4-virulence occurred in isolates with different genealogies. The V94-5152 (Rsv4)-virulence determinant in Iranian isolates maps into a different amino acid position in the P3 protein than those previously reported, indicating different evolutionary pathways towards resistance breaking that might be conditioned by sequence context.

First dating of a recombination event in mammalian tick-borne flaviviruses

The mammalian tick-borne flavivirus group (MTBFG) contains viruses associated with important human and animal diseases such as encephalitis and hemorrhagic fever. In contrast to mosquito-borne flaviviruses where recombination events are frequent, the evolutionary dynamic within the MTBFG was believed to be essentially clonal. This assumption was challenged with the recent report of several homologous recombinations within the Tick-borne encephalitis virus (TBEV). We performed a thorough analysis of publicly available genomes in this group and found no compelling evidence for the previously identified recombinations. However, our results show for the first time that demonstrable recombination (i.e., with large statistical support and strong phylogenetic evidences) has occurred in the MTBFG, more specifically within the Louping ill virus lineage. Putative parents, recombinant strains and breakpoints were further tested for statistical significance using phylogenetic methods. We investigated the time of divergence between the recombinant and parental strains in a Bayesian framework. The recombination was estimated to have occurred during a window of 282 to 76 years before the present. By unravelling the temporal setting of the event, we adduce hypotheses about the ecological conditions that could account for the observed recombination.

Ongoing purifying selection on intergenic spacers in group A streptococcus

Bacterial intergenic spacers are non-coding genomic regions enriched with cis-regulatory elements for gene expression. A population genetics approach was used to investigate the evolutionary force shaping the genetic diversity of intergenic spacers among 13 genomes of group A streptococcus (GAS). Analysis of 590 genes and their linked 5' intergenic spacers showed reduced nucleotide diversity in spacers compared to synonymous nucleotide diversity in protein-coding regions, suggestive of past purifying selection on spacers. Certain spacers showed elevated nucleotide diversity indicative of past homologous recombination with divergent genotypes. In addition, analysis of the difference between mean nucleotide difference and number of segregating sites showed evidence of an excess of rare variants both at nonsynonymous sites in genes and at sites in spacers, which is evidence that there are numerous slightly deleterious variants in GAS populations with potential effects on both protein sequences and gene expression.

Mapping the landscape of host-pathogen coevolution: HLA class I binding and its relationship with evolutionary conservation in human and viral proteins

The high diversity of HLA binding preferences has been driven by the sequence diversity of short segments of relevant pathogenic proteins presented by HLA molecules to the immune system. To identify possible commonalities in HLA binding preferences, we quantify these using a novel measure termed "targeting efficiency," which captures the correlation between HLA-peptide binding affinities and the conservation of the targeted proteomic regions. Analysis of targeting efficiencies for 95 HLA class I alleles over thousands of human proteins and 52 human viruses indicates that HLA molecules preferentially target conserved regions in these proteomes, although the arboviral Flaviviridae are a notable exception where nonconserved regions are preferentially targeted by most alleles. HLA-A alleles and several HLA-B alleles that have maintained close sequence identity with chimpanzee homologues target conserved human proteins and DNA viruses such as Herpesviridae and Adenoviridae most efficiently, while all HLA-B alleles studied efficiently target RNA viruses. These patterns of host and pathogen specialization are both consistent with coevolutionary selection and functionally relevant in specific cases; for example, preferential HLA targeting of conserved proteomic regions is associated with improved outcomes in HIV infection and with protection against dengue hemorrhagic fever. Efficiency analysis provides a novel perspective on the coevolutionary relationship between HLA class I molecular diversity, self-derived peptides that shape T-cell immunity through ontogeny, and the broad range of viruses that subsequently engage with the adaptive immune response.

The evolutionary biology of poxviruses

The poxviruses (family Poxviridae) are a family of double-stranded viruses including several species that infect humans and their domestic animals, most notably Variola virus (VARV), the causative agent of smallpox. The evolutionary biology of these viruses poses numerous questions, for which we have only partial answers at present. Here we review evidence regarding the origin of poxviruses, the frequency of host transfer in poxvirus history, horizontal transfer of host genes to poxviruses, and the population processes accounting for patterns of nucleotide sequence polymorphism.

Small effective population sizes and rare nonsynonymous variants in potyviruses

Analysis of nucleotide sequence polymorphism in complete genomes of 12 species of potyviruses (single-stranded, positive-sense RNA viruses, family Potyviridae) revealed evidence that long-term effective population sizes of these viruses are on the order of 10(4). Comparison of nucleotide diversity in non-coding regions and at synonymous and nonsynonymous sites in coding regions showed that purifying selection has acted to eliminate numerous deleterious mutations both at nonsynonymous sites and in non-coding regions. The ratio of nonsynonymous to synonymous polymorphic sites increased as a function of the number of genomes sampled, whereas mean gene diversity at nonsynonymous polymorphic sites decreased with increasing sample size at a substantially faster rate than does mean gene diversity at synonymous polymorphic sites. Very similar relationships were observed both in available genomic sequences of 12 potyvirus species and in subsets created by randomly sampling from among 98 TuMV genomes. Taken together, these observations imply that a greater proportion of nonsynonymous than of synonymous variants are relatively rare as the result of ongoing purifying selection, and thus many nonsynonymous variants are unlikely to be discovered without extensive sampling.

Detecting natural selection in RNA virus populations using sequence summary statistics

At present, most analyses that aim to detect the action of natural selection upon viral gene sequences use phylogenetic estimates of the ratio of silent to replacement mutations. Such methods, however, are impractical to compute on large data sets comprising hundreds of complete viral genomes, which are becoming increasingly common due to advances in genome sequencing technology. Here we investigate the statistical performance of computationally efficient tests that are based on sequence summary statistics, and explore their applicability to RNA virus data sets in two ways. Firstly, we perform extensive simulations in order to measure the type I error of two well-known summary statistic methods - Tajima's D and the McDonald-Kreitman test - under a range of virus-like mutational and demographic scenarios. Secondly, we apply these methods to a compilation of approximately 100 RNA virus alignments that represent natural RNA virus populations. In addition, we develop and introduce a new implementation of the McDonald-Kreitman test and show that it greatly improves the test's statistical reliability on typical viral data sets. Our results suggest that variants of the McDonald-Kreitman test could prove useful in the analysis of very large sets of highly diverse viral genetic data.

Relaxation of purifying selection on the SAD lineage of live attenuated oral vaccines for rabies virus

Analysis of patterns of nucleotide sequence diversity in wild-type rabies virus (RABV) genomes and in the SAD live attenuated oral vaccine lineage was used to test for the relaxation of purifying selection in the latter and provide evidence regarding the genomic regions where such relaxation of selection occurs. The wild-type sequences showed evidence of strong past and ongoing purifying selection both on nonsynonymous sites in coding regions and on non-coding regions, particularly the start, end and 5' UTR regions. SAD vaccine sequences showed a relaxation of purifying selection at nonsynonymous sites in coding regions, resulting a substantial number of amino acid sequence polymorphisms at sites that were invariant in the wild-type sequences. Moreover, SAD vaccine sequences showed high levels of mutation accumulation in the non-coding regions that were most conserved in the wild-type sequences. Understanding the biological effects of the unique mutations accumulated in the vaccine lineage is important because of their potential effects on antigenicity and effectiveness of the vaccine.

More radical amino acid replacements in primates than in rodents: support for the evolutionary role of effective population size

We examined the pattern of nucleotide substitution in 4933 conserved single-copy orthologous protein-coding genes of human, rhesus, mouse, and rat. Consistent with previous studies, the median ratio of the number of nonsynonymous substitutions per nonsynonymous site (d(N)) to the number of synonymous substitutions per synonymous site (d(S)) was significantly higher in the comparison between the two primates than in the comparison between the two rodents. This pattern was particularly strong in the case of genes expressed in the immune system, but also occurred in other genes, including a set of highly conserved genes involved in the regulation of transcription. Both synonymous and nonsynonymous differences occurred independently in the same codons in the primates and in the rodents to a greater extent than expected by chance, but the extent of the deviation from random expectation was much greater in the case of nonsynonymous differences. Parallel amino acid replacements occurred at the same sites in the primates and rodents far more frequently than expected by chance, but tended to involve very conservative amino acid changes. Divergent amino acid changes involved more chemically different amino acids than parallel changes, and divergent amino acid replacements between the primates were significantly more radical than those between the rodents. These results are most easily explained on the hypothesis that the evolution of these genes has been shaped largely by purifying selection, which has been less effective in primates than in rodents, presumably as a consequence of lower long-term effective population sizes in the former.

Relaxation of purifying selection on live attenuated vaccine strains of the family Paramyxoviridae

In wild-type sequences of three paramyxoviruses (measles virus, mumps virus, and Newcastle disease virus), nucleotide diversity at both non-coding sites and at nonsynonymous sites in coding regions was significantly reduced in comparison to that at synonymous sites. Likewise, both the mean and variance of gene diversity at nonsynonymous polymorphic sites were reduced in comparison to non-coding and synonymous sites. Neither of these patterns, which reflect the action of purifying selection against deleterious mutations at nonsynonymous and non-coding sites, were seen in the case of live attenuated vaccine strains, implying that purifying selection has been substantially relaxed on the latter, potentially affecting their biological properties, including antigenicity and vaccine effectiveness. Since the accumulation of mutations increases as a function of the number of generations of replication, these findings highlight the utility of minimizing the number of generations between the original vaccine master seed and the strains used in vaccination, along with periodic monitoring of the extent of sequence evolution.

Synonymous and nonsynonymous polymorphisms versus divergences in bacterial genomes

Comparison of the ratio of nonsynonymous to synonymous polymorphisms within species with the ratio of nonsynonymous to synonymous substitutions between species has been widely used as a supposed indicator of positive Darwinian selection, with the ratio of these 2 ratios being designated as a neutrality index (NI). Comparison of genome-wide polymorphism within 12 species of bacteria with divergence from an outgroup species showed substantial differences in NI among taxa. A low level of nonsynonymous polymorphism at a locus was the best predictor of NI < 1, rather than a high level of nonsynonymous substitution between species. Moreover, genes with NI < 1 showed a strong tendency toward the occurrence of rare nonsynonymous polymorphisms, as expected under the action of ongoing purifying selection. Thus, our results are more consistent with the hypothesis that a high relative rate of between-species nonsynonymous substitution reflects mainly the action of purifying selection within species to eliminate slightly deleterious mutations rather than positive selection between species. This conclusion is consistent with previous results highlighting an important role of slightly deleterious variants in bacterial evolution and suggests caution in the use of the McDonald-Kreitman test and related statistics as tests of positive selection.

Distinctive pattern of sequence polymorphism in the NS3 protein of hepatitis C virus type 1b reflects conflicting evolutionary pressures

Analysis of complete polyprotein-encoding sequences of hepatitis C virus genotype 1b (HCV-1b) showed evidence not only of past purifying selection but also of abundant slightly deleterious non-synonymous variants subject to ongoing purifying selection. The NS3 protein (with protease and NTPase/helicase activity) revealed less evidence of purifying selection acting on the cytotoxic T cells (CTL) epitopes than did the other proteins, whereas outside the CTL epitopes NS3 was more conserved than the other proteins. Moreover, NS3 showed a high incidence of forward-and-backward or parallel non-synonymous changes in CTL epitopes, as measured by the consistency index across the phylogeny of HCV-1b genomes computed at non-singleton non-synonymous polymorphic sites. This result implies that certain non-synonymous mutations have recurred frequently throughout the phylogeny in the codons encoding the epitopes in NS3. This pattern is most easily explained by the frequent re-occurrence of the same set of escape mutations in CTL epitopes of NS3, which are selectively favoured within hosts expressing the presenting class I major histocompatibility complex molecule, but are subject to purifying selection at the population level. The fact that this pattern is most strikingly observed in the case of NS3 suggests that the evolutionary conflict between immune escape and functional constraint on the protein is more acute in the case of NS3 than any of the other proteins of HCV-1b.