FUBAR: a fast, unconstrained bayesian approximation for inferring selection

A Comparison of One-Rate and Two-Rate Inference Frameworks for Site-Specific dN/dS Estimation

Two broad paradigms exist for inferring [Formula: see text] the ratio of nonsynonymous to synonymous substitution rates, from coding sequences: (i) a one-rate approach, where [Formula: see text] is represented with a single parameter, or (ii) a two-rate approach, where [Formula: see text] and [Formula: see text] are estimated separately. The performances of these two approaches have been well studied in the specific context of proper model specification, i.e., when the inference model matches the simulation model. By contrast, the relative performances of one-rate vs. two-rate parameterizations when applied to data generated according to a different mechanism remain unclear. Here, we compare the relative merits of one-rate and two-rate approaches in the specific context of model misspecification by simulating alignments with mutation-selection models rather than with [Formula: see text]-based models. We find that one-rate frameworks generally infer more accurate [Formula: see text] point estimates, even when [Formula: see text] varies among sites. In other words, modeling [Formula: see text] variation may substantially reduce accuracy of [Formula: see text] point estimates. These results appear to depend on the selective constraint operating at a given site. For sites under strong purifying selection ([Formula: see text]), one-rate and two-rate models show comparable performances. However, one-rate models significantly outperform two-rate models for sites under moderate-to-weak purifying selection. We attribute this distinction to the fact that, for these more quickly evolving sites, a given substitution is more likely to be nonsynonymous than synonymous. The data will therefore be relatively enriched for nonsynonymous changes, and modeling [Formula: see text] contributes excessive noise to [Formula: see text] estimates. We additionally find that high levels of divergence among sequences, rather than the number of sequences in the alignment, are more critical for obtaining precise point estimates.

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

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

Nonstructural Proteins Are Preferential Positive Selection Targets in Zika Virus and Related Flaviviruses

The Flavivirus genus comprises several human pathogens such as dengue virus (DENV), Japanese encephalitis virus (JEV), and Zika virus (ZIKV). Although ZIKV usually causes mild symptoms, growing evidence is linking it to congenital birth defects and to increased risk of Guillain-Barré syndrome. ZIKV encodes a polyprotein that is processed to produce three structural and seven nonstructural (NS) proteins. We investigated the evolution of the viral polyprotein in ZIKV and in related flaviviruses (DENV, Spondweni virus, and Kedougou virus). After accounting for saturation issues, alignment uncertainties, and recombination, we found evidence of episodic positive selection on the branch that separates DENV from the other flaviviruses. NS1 emerged as the major selection target, and selected sites were located in immune epitopes or in functionally important protein regions. Three of these sites are located in an NS1 region that interacts with structural proteins and is essential for virion biogenesis. Analysis of the more recent evolutionary history of ZIKV lineages indicated that positive selection acted on NS5 and NS4B, this latter representing the preferential target. All selected sites were located in the N-terminal portion of NS4B, which inhibits interferon response. One of the positively selected sites (26M/I/T/V) in ZIKV also represents a selection target in sylvatic DENV2 isolates, and a nearby residue evolves adaptively in JEV. Two additional positively selected sites are within a protein region that interacts with host (e.g. STING) and viral (i.e. NS1, NS4A) proteins. Notably, mutations in the NS4B region of other flaviviruses modulate neurovirulence and/or neuroinvasiveness. These results suggest that the positively selected sites we identified modulate viral replication and contribute to immune evasion. These sites should be prioritized in future experimental studies. However, analyses herein detected no selective events associated to the spread of the Asian/American ZIKV lineage.

Comparative genomic evidence for duplication of TLR1 subfamily and miiuy croaker TLR1 perceives LPS stimulation via MyD88 and TIRAP

Being indispensable pattern recognition receptors in innate immune responses in host protection, Toll-like receptors (TLRs) play an important role in pathogen recognition. Fish TLRs exhibit high variety and distinct features, although little is known about their function on ligand recognition and signaling pathway in fish. This paper reports the evolutionary spectrum of the TLR1 subfamily (referred to as TLR1, TLR6, and TLR10) as determined using the comparative genomic approach. We hypothesized that the TLR1 subfamily underwent two rounds of gene duplication events; the first duplication occurred prior to the divergence of amphibians, and the second one occurred prior to the divergence of eutherians. To further study the function of fish TLR1, we identified miiuy croaker (Miichthys miiuy) TLR1 (mmiTLR1) and determined its potential ability to perceive Vibrio anguillarum and lipopolysaccharide stimulation. Data further suggested that mmiTLR1 is dependent on TIRAP and MyD88 for signal transmission. In addition, immunocytochemistry showed the speculative interaction between MyD88 and mmiTLR1 TIR domain. Overall, we systematically and comprehensively analyzed evolution of TLR1 subfamily and the function of mmiTLR1, which will provide the basis for future scientific research on fish TLRs.

The Evolutionary Histories of Antiretroviral Proteins SERINC3 and SERINC5 Do Not Support an Evolutionary Arms Race in Primates

Molecular evolutionary arms races between viruses and their hosts are important drivers of adaptation. These Red Queen dynamics have been frequently observed in primate retroviruses and their antagonists, host restriction factor genes, such as APOBEC3F/G, TRIM5-α, SAMHD1, and BST-2. Host restriction factors have experienced some of the most intense and pervasive adaptive evolution documented in primates. Recently, two novel host factors, SERINC3 and SERINC5, were identified as the targets of HIV-1 Nef, a protein crucial for the optimal infectivity of virus particles. Here, we compared the evolutionary fingerprints of SERINC3 and SERINC5 to those of other primate restriction factors and to a set of other genes with diverse functions. SERINC genes evolved in a manner distinct from the canonical arms race dynamics seen in the other restriction factors. Despite their antiviral activity against HIV-1 and other retroviruses, SERINC3 and SERINC5 have a relatively uneventful evolutionary history in primates.

IMPORTANCE

Restriction factors are host proteins that block viral infection and replication. Many viruses, like HIV-1 and related retroviruses, evolved accessory proteins to counteract these restriction factors. The importance of these interactions is evidenced by the intense adaptive selection pressures that dominate the evolutionary histories of both the host and viral genes involved in this so-called arms race. The dynamics of these arms races can point to mechanisms by which these viral infections can be prevented. Two human genes, SERINC3 and SERINC5, were recently identified as targets of an HIV-1 accessory protein important for viral infectivity. Unexpectedly, we found that these SERINC genes, unlike other host restriction factor genes, show no evidence of a recent evolutionary arms race with viral pathogens.

Molecular Evolution and Phylogeography of Co-circulating IHNV and VHSV in Italy

Infectious haematopoietic necrosis virus (IHNV) and viral haemorrhagic septicaemia virus (VHSV) are the most important viral pathogens impacting rainbow trout farming. These viruses are persistent in Italy, where they are responsible for severe disease outbreaks (epizootics) that affect the profitability of the trout industry. Despite the importance of IHNV and VHSV, little is known about their evolution at a local scale, although this is likely to be important for virus eradication and control. To address this issue we performed a detailed molecular evolutionary and epidemiological analysis of IHNV and VHSV in trout farms from northern Italy. Full-length glycoprotein gene sequences of a selection of VHSV (n = 108) and IHNV (n = 89) strains were obtained. This revealed that Italian VHSV strains belong to sublineages Ia1 and Ia2 of genotype Ia and are distributed into 7 genetic clusters. In contrast, all Italian IHNV isolates fell within genogroup E, for which only a single genetic cluster was identified. More striking was that IHNV has evolved more rapidly than VHSV (mean rates of 11 and 7.3 × 10⁻⁴ nucleotide substitutions per site, per year, respectively), indicating that these viruses exhibit fundamentally different evolutionary dynamics. The time to the most recent common ancestor of both IHNV and VHSV was consistent with the first reports of these pathogens in Italy. By combining sequence data with epidemiological information it was possible to identify different patterns of virus spread among trout farms, in which adjacent facilities can be infected by either genetically similar or different viruses, and farms located in different water catchments can be infected by identical strains. Overall, these findings highlight the importance of combining molecular and epidemiological information to identify the determinants of IHN and VHS spread, and to provide data that is central to future surveillance strategies and possibly control.

Characterisation and expression analysis of UBC9 and UBS27 genes in developing gonads of cicindelids (Coleoptera: Cicindelidae)

Ubiquitin and small ubiquitin-like modifiers (SUMO) are post-translational modifiers essential in a variety of cellular processes, including gametogenesis. SUMO-conjugating enzyme (UBC9) and the ubiquitin ribosomal fusion protein UBS27 have been characterised in several model species. However, their expression in coleopteran remains unstudied. In this study, UBC9 and UBS27 genes have been characterised in the tiger beetle Cicindela campestris for the first time. Bioinformatic analysis showed that the Cc-UBC9 gene encoded a 159 amino acid protein with a predicted molecular weight of 18.18kDa, and the Cc-UBS27 gene encoded a 156 amino acid protein with a predicted molecular weight of 17.71kDa. Selection analyses carried out in several cicindelid species revealed that both genes were affected by purifying selection. Real time quantitative PCR analysis demonstrated that Cc-UBC9 and Cc-UBS27 were expressed in different tissues. The highest expression on both genes was found in the ovary and testis, and there were differential expression levels between immature and mature stages of testis development. The expression patterns of Cc-UBC9 and Cc-UBS27 suggest that these genes play important roles in gametogenesis in C. campestris. This information is relevant to better understand the reproductive process in cicindelids and the function of ubiquitin and small ubiquitin-related modifier genes in the Coleoptera.

Adaptive evolution of virus-sensing toll-like receptor 8 in bats

Recently, bats have gained attention as potential reservoir hosts for emerging zoonotic single-stranded (ssRNA) viruses that may prove fatal for humans and other mammals. It has been hypothesized that some features of their innate immune system may enable bats to trigger an efficient early immune response. Toll-like receptors (TLRs) represent a first line defense within the innate immune system and lie directly at the host–pathogen interface in targeting specific microbe-molecular patterns. However, the direction and strength of selection acting on TLRs are largely unknown for bats. Here, we studied the selection on viral ssRNA sensing TLR8 based on sequence data of 21 bat species. The major part (63 %) of the TLR8 gene evolved under purifying selection, likely due to functional constraints. We also found evidence for persistent positive selection acting on specific amino acid sites (7 %), especially when compared to viral TLR evolution of other mammals. All of these putatively positively selected codons were located in the ligand-binding ectodomain, some coincidenced or were in close proximity to functional sites, as suggested by the crystallographic structure of the human TLR8. This might contribute to the inter-species variation in the ability to recognize molecular patterns of viruses. TLR8 evolution within bats revealed that branches leading to ancestral and recent lineages evolved under episodic positive selection, indicating selective selection pressures in restricted bat lineages. Altogether, we found that the TLR8 displays extensive sequence variation within bats and that unique features separate them from humans and other mammals.

Environmental variability and the evolution of the glucocorticoid receptor (Nr3c1) in African starlings

One of the primary ways that organisms cope with environmental change is through regulation of the hypothalamo-pituitary-adrenal (HPA) axis, the neuroendocrine system that controls reactions to stress. Variation in genes regulating the HPA axis - particularly the glucocorticoid receptor - may facilitate adaptation to changing climatic conditions by altering expression. Here we examine signatures of selection on the glucocorticoid receptor gene (Nr3c1) in African starlings that inhabit a range of environments, including those with variable climatic conditions. To investigate potential adaptive mechanisms underlying the vertebrate stress response, we sequence the Nr3c1 gene in 27 species of African starlings. Although we find some evidence of positive selection, substitution rate is negatively correlated with variance in precipitation. This suggests climatic cycling in sub-Saharan Africa may have resulted in lower substitution rates to maintain a successful coping strategy. When environmental conditions fluctuate rapidly, variation in the strength of purifying selection can explain evolutionary rate variation.

Analyzing endocrine system conservation and evolution

Analyzing variation in rates of evolution can provide important insights into the factors that constrain trait evolution, as well as those that promote diversification. Metazoan endocrine systems exhibit apparent variation in evolutionary rates of their constituent components at multiple levels, yet relatively few studies have quantified these patterns and analyzed them in a phylogenetic context. This may be in part due to historical and current data limitations for many endocrine components and taxonomic groups. However, recent technological advancements such as high-throughput sequencing provide the opportunity to collect large-scale comparative data sets for even non-model species. Such ventures will produce a fertile data landscape for evolutionary analyses of nucleic acid and amino acid based endocrine components. Here I summarize evolutionary rate analyses that can be applied to categorical and continuous endocrine traits, and also those for nucleic acid and protein-based components. I emphasize analyses that could be used to test whether other variables (e.g., ecology, ontogenetic timing of expression, etc.) are related to patterns of rate variation and endocrine component diversification. The application of phylogenetic-based rate analyses to comparative endocrine data will greatly enhance our understanding of the factors that have shaped endocrine system evolution.

Human Management of a Wild Plant Modulates the Evolutionary Dynamics of a Gene Determining Recessive Resistance to Virus Infection

This work analyses the genetic variation and evolutionary patterns of recessive resistance loci involved in matching-allele (MA) host-pathogen interactions, focusing on the pvr2 resistance gene to potyviruses of the wild pepper Capsicum annuum glabriusculum (chiltepin). Chiltepin grows in a variety of wild habitats in Mexico, and its cultivation in home gardens started about 25 years ago. Potyvirus infection of Capsicum plants requires the physical interaction of the viral VPg with the pvr2 product, the translation initiation factor eIF4E1. Mutations impairing this interaction result in resistance, according to the MA model. The diversity of pvr2/eIF4E1 in wild and cultivated chiltepin populations from six biogeographical provinces in Mexico was analysed in 109 full-length coding sequences from 97 plants. Eleven alleles were found, and their interaction with potyvirus VPg in yeast-two-hybrid assays, plus infection assays of plants, identified six resistance alleles. Mapping resistance mutations on a pvr2/eIF4E1 model structure showed that most were around the cap-binding pocket and strongly altered its surface electrostatic potential, suggesting resistance-associated costs due to functional constraints. The pvr2/eIF4E1 phylogeny established that susceptibility was ancestral and resistance was derived. The spatial structure of pvr2/eIF4E1 diversity differed from that of neutral markers, but no evidence of selection for resistance was found in wild populations. In contrast, the resistance alleles were much more frequent, and positive selection stronger, in cultivated chiltepin populations, where diversification of pvr2/eIF4E1 was higher. This analysis of the genetic variation of a recessive resistance gene involved in MA host-pathogen interactions in populations of a wild plant show that evolutionary patterns differ according to the plant habitat, wild or cultivated. It also demonstrates that human management of the plant population has profound effects on the diversity and the evolution of the resistance gene, resulting in the selection of resistance alleles.

Genetic Characterization of Northwestern Colombian Chikungunya Virus Strains from the 2014-2015 Epidemic

Chikungunya fever, an acute and often chronic arthralgic disease caused by the mosquito-borne alphavirus, chikungunya virus (CHIKV), spread into the Americas in late 2013. Since then it has caused epidemics in nearly all New World countries, the second largest being Colombia with over 450,000 suspected cases beginning in September, 2014, and focused in Bolivar Department in the north. We examined 32 human sera from suspected cases, including diverse age groups and both genders, and sequenced the CHIKV envelope glycoprotein genes, known determinants of vector host range. As expected for Asian lineage CHIKV strains, these isolates lacked known Aedes albopictus-adaptive mutations. All the Colombian strains were closely related to those from the Virgin Islands, Saint Lucia, Mexico, Puerto Rico, and Brazil, consistent with a single, point-source introduction from the southeast Asia/Pacific region. Two substitutions in the E2 and E1 envelope glycoprotein genes were found in the Colombian strains, especially E1-K211E involving a residue shown previously to affect epistatically the penetrance of the E1-A226V A. albopictus-adaptive substitution. We also identified two amino acid substitutions unique to all American CHIKV sequences: E2-V368A and 6K-L20M. Only one codon, 6K-47, had a high nonsynonymous substitution rate suggesting positive selection.

Rapid Sequencing of Complete env Genes from Primary HIV-1 Samples

The ability to study rapidly evolving viral populations has been constrained by the read length of next-generation sequencing approaches and the sampling depth of single-genome amplification methods. Here, we develop and characterize a method using Pacific Biosciences' Single Molecule, Real-Time (SMRT®) sequencing technology to sequence multiple, intact full-length human immunodeficiency virus-1 env genes amplified from viral RNA populations circulating in blood, and provide computational tools for analyzing and visualizing these data.

Identification of 48 full-length MHC-DAB functional alleles in miiuy croaker and evidence for positive selection

Major histocompatibility complex (MHC) molecules play a vital role in the immune response and are a highly polymorphic gene superfamily in vertebrates. As the molecular marker associated with polymorphism and disease susceptibility/resistance, the polymorphism of MHC genes has been investigated in many tetrapods and teleosts. Most studies were focused on the polymorphism of the second exon, which encodes the peptide-binding region (PBR) in the α1- or β1-domain, but few studies have examined the full-length coding region. To comprehensive investigate the polymorphism of MHC gene, we identified 48 full-length miiuy croaker (Miichthys miiuy) MHC class IIB (Mimi-DAB) functional alleles from 26 miiuy croaker individuals. All of the alleles encode 34 amino acid sequences, and a high level of polymorphism was detected in Mimi-DAB alleles. The rate of non-synonymous substitutions (dN) occurred at a significantly higher frequency than that of synonymous substitutions (dS) in the PBR, and this result suggests that balancing selection maintains polymorphisms at the Mimi-DAB locus. Phylogenetic analysis based on the full-length and exon 2 sequences of Mimi-DAB alleles both showed that the Mimi-DAB alleles were clustered into two major groups. A total of 19 positive selected sites were identified on the Mimi-DAB alleles after testing for positive selection, and 14 sites were predicted to be associated with antigen-binding sites, which suggests that most of selected sites are significant for disease resistance. The polymorphism of Mimi-DAB alleles provides an important resource for analyzing the association between the polymorphism of MHC gene and disease susceptibility/resistance, and for researching the molecular selective breeding of miiuy croaker with enhanced disease resistance.

Evolutionary inactivation of a sialidase in group B Streptococcus

Group B Streptococcus (GBS) is a leading cause of bacterial sepsis and meningitis in newborns. GBS possesses a protein with homology to the pneumococcal virulence factor, NanA, which has neuraminidase (sialidase) activity and promotes blood-brain barrier penetration. However, phylogenetic sequence and enzymatic analyses indicate the GBS NanA ortholog has lost sialidase function – and for this distinction we designate the gene and encoded protein nonA/NonA. Here we analyze NonA function in GBS pathogenesis, and through heterologous expression of active pneumococcal NanA in GBS, potential costs of maintaining sialidase function. GBS wild-type and ΔnonA strains lack sialidase activity, but forced expression of pneumococcal NanA in GBS induced degradation of the terminal sialic acid on its exopolysaccharide capsule. Deletion of nonA did not change GBS-whole blood survival or brain microvascular cell invasion. However, forced expression of pneumococcal NanA in GBS removed terminal sialic acid residues from the bacterial capsule, restricting bacterial proliferation in human blood and in vivo upon mouse infection. GBS expressing pneumococcal NanA had increased invasion of human brain microvascular endothelial cells. Thus, we hypothesize that nonA lost enzyme activity allowing the preservation of an effective survival factor, the sialylated exopolysaccharide capsule.

Out of the blue: adaptive visual pigment evolution accompanies Amazon invasion

Incursions of marine water into South America during the Miocene prompted colonization of freshwater habitats by ancestrally marine species and present a unique opportunity to study the molecular evolution of adaptations to varying environments. Freshwater and marine environments are distinct in both spectra and average intensities of available light. Here, we investigate the molecular evolution of rhodopsin, the photosensitive pigment in the eye that activates in response to light, in a clade of South American freshwater anchovies derived from a marine ancestral lineage. Using likelihood-based comparative sequence analyses, we found evidence for positive selection in the rhodopsin of freshwater anchovy lineages at sites known to be important for aspects of rhodopsin function such as spectral tuning. No evidence was found for positive selection in marine lineages, nor in three other genes not involved in vision. Our results suggest that an increased rate of rhodopsin evolution was driven by diversification into freshwater habitats, thereby constituting a rare example of molecular evolution mirroring large-scale palaeogeographic events.

Phylodynamic analysis of porcine circovirus type 2: Methodological approach and datasets

Since its first description, PCV2 has emerged as one of the most economically relevant diseases for the swine industry. Despite the introduction of vaccines effective in controlling clinical syndromes, PCV2 spread was not prevented and some potential evidences of vaccine immuno escape have recently been reported ("Complete genome sequence of a novel porcine circovirus type 2b variant present in cases of vaccine failures in the United States" (Xiao and Halbur, 2012) [1], "Genetic and antigenic characterization of a newly emerging porcine circovirus type 2b mutant first isolated in cases of vaccine failure in Korea" (Seo et al., 2014) [2]). In this article, we used a collection of PCV2 full genomes, provided in the present manuscript, and several phylogentic, phylodynamic and bioinformatic methods to investigate different aspects of PCV2 epidemiology, history and evolution (more thoroughly described in "PHYLODYNAMIC ANALYSIS of PORCINE CIRCOVIRUS TYPE 2 REVEALS GLOBAL WAVES of EMERGING GENOTYPES and the CIRCULATION of RECOMBINANT FORMS"[3]). The methodological approaches used to consistently detect recombiantion events and estimate population dymanics and spreading patterns of rapidly evolving ssDNA viruses are herein reported. Programs used are described and original scripts have been provided. Ensembled databases used are also made available. These consist of a broad collection of complete genome sequences (i.e. 843 sequences; 63 complete genomes of PCV2a, 310 of PCV2b, 4 of PCV2c, 217 of PCV2d, 64 of CRF01, 140 of CRF02 and 45 of CRF03.), divided in differnt ORF (i.e. ORF1, ORF2 and intergenic regions), of PCV2 genotypes and major Circulating Recombinat Forms (CRF) properly annotated with respective collection data and country. Globally, all of these data can be used as a starting point for further studies and for classification purpose.

Evolution of Neuroadaptation in the Periphery and Purifying Selection in the Brain Contribute to Compartmentalization of Simian Immunodeficiency Virus (SIV) in the Brains of Rhesus Macaques with SIV-Associated Encephalitis

The emergence of a distinct subpopulation of human or simian immunodeficiency virus (HIV/SIV) sequences within the brain (compartmentalization) during infection is hypothesized to be linked to AIDS-related central nervous system (CNS) neuropathology. However, the exact evolutionary mechanism responsible for HIV/SIV brain compartmentalization has not been thoroughly investigated. Using extensive viral sampling from several different peripheral tissues and cell types and from three distinct regions within the brain from two well-characterized rhesus macaque models of the neurological complications of HIV infection (neuroAIDS), we have been able to perform in-depth evolutionary analyses that have been unattainable in HIV-infected subjects. The results indicate that, despite multiple introductions of virus into the brain over the course of infection, brain sequence compartmentalization in macaques with SIV-associated CNS neuropathology likely results from late viral entry of virus that has acquired through evolution in the periphery sufficient adaptation for the distinct microenvironment of the CNS.

IMPORTANCE

HIV-associated neurocognitive disorders remain prevalent among HIV type 1-infected individuals, whereas our understanding of the critical components of disease pathogenesis, such as virus evolution and adaptation, remains limited. Building upon earlier findings of specific viral subpopulations in the brain, we present novel yet fundamental results concerning the evolutionary patterns driving this phenomenon in two well-characterized animal models of neuroAIDS and provide insight into the timing of entry of virus into the brain and selective pressure associated with viral adaptation to this particular microenvironment. Such knowledge is invaluable for therapeutic strategies designed to slow or even prevent neurocognitive impairment associated with AIDS.

Grapevine red blotch-associated virus is Present in Free-Living Vitis spp. Proximal to Cultivated Grapevines

Red blotch is an emerging disease of grapevine associated with grapevine red blotch-associated virus (GRBaV). The virus spreads with infected planting stocks but no vector of epidemiological significance has been conclusively identified. A vineyard block of red-blotch-affected Vitis vinifera 'Cabernet franc' clone 214 was observed in California, with a clustering of infected, symptomatic vines focused along one edge of the field proximal to a riparian habitat with free-living Vitis spp. No genetic heterogeneity was observed in a 587-nucleotide region of the GRBaV genome in a population of 44 Cabernet franc clone 214 isolates. By contrast, genetic differences were observed in isolates from other cultivars and clones growing in adjacent blocks. GRBaV was confirmed infecting four free-living vines, two of which were shown to be V. californica × V. vinifera hybrids. The genomes of three free-living GRBaV vine isolates and seven from V. vinifera cultivars were compared; free-living vine isolates were shown to be more similar to each other and a 'Merlot' isolate than to the other cultivated vine isolates. The finding that GRBaV is present in free-living Vitis spp. indicates the virus can be spread by natural (nonhuman-mediated) means, and we hypothesize that in-field spread of GRBaV is occurring.

Bioinformatics analysis of the structural and evolutionary characteristics for toll-like receptor 15

Toll-like receptors (TLRs) play important role in the innate immune system. TLR15 is reported to have a unique role in defense against pathogens, but its structural and evolution characterizations are still poorly understood. In this study, we identified 57 completed TLR15 genes from avian and reptilian genomes. TLR15 clustered into an individual clade and was closely related to family 1 on the phylogenetic tree. Unlike the TLRs in family 1 with the broken asparagine ladders in the middle, TLR15 ectodomain had an intact asparagine ladder that is critical to maintain the overall shape of ectodomain. The conservation analysis found that TLR15 ectodomain had a highly evolutionarily conserved region on the convex surface of LRR11 module, which is probably involved in TLR15 activation process. Furthermore, the protein-protein docking analysis indicated that TLR15 TIR domains have the potential to form homodimers, the predicted interaction interface of TIR dimer was formed mainly by residues from the BB-loops and αC-helixes. Although TLR15 mainly underwent purifying selection, we detected 27 sites under positive selection for TLR15, 24 of which are located on its ectodomain. Our observations suggest the structural features of TLR15 which may be relevant to its function, but which requires further experimental validation.

Evolution of CCL11: genetic characterization in lagomorphs and evidence of positive and purifying selection in mammals

The interactions between chemokines and their receptors are crucial for differentiation and activation of inflammatory cells. CC chemokine ligand 11 (CCL11) binds to CCR3 and to CCR5 that in leporids underwent gene conversion with CCR2. Here, we genetically characterized CCL11 in lagomorphs (leporids and pikas). All lagomorphs have a potentially functional CCL11, and the Pygmy rabbit has a mutation in the stop codon that leads to a longer protein. Other mammals also have mutations at the stop codon that result in proteins with different lengths. By employing maximum likelihood methods, we observed that, in mammals, CCL11 exhibits both signatures of purifying and positive selection. Signatures of purifying selection were detected in sites important for receptor binding and activation. Of the three sites detected as under positive selection, two were located close to the stop codon. Our results suggest that CCL11 is functional in all lagomorphs, and that the signatures of purifying and positive selection in mammalian CCL11 probably reflect the protein's biological roles.

Genetic characterization of Enterovirus 71 strains circulating in Vietnam in 2012

BACKGROUND

Enterovirus 71 subgenogroup C4 caused the largest outbreak of Hand, Foot and Mouth Disease (HFMD) in Vietnam during 2011-2012, resulting in over 200,000 hospitalisations and 207 fatalities.

METHODS

A total of 1917 samples with adequate volume for RT-PCR analysis were collected from patients hospitalised with HFMD throughout Vietnam and 637 were positive for EV71. VP1 gene (n=87) and complete genome (n=9) sequencing was performed. Maximum-likelihood phylogenetic analysis was performed to characterise the B5, C4 and C5 strains detected.

RESULTS

Sequence analyses revealed that the dominant subgenogroup associated with the 2012 outbreak was C4, with B5 and C5 strains representing a small proportion of these cases.

CONCLUSIONS

Numerous countries in the region including Malaysia, Taiwan and China have a large influence on strain diversity in Vietnam and understanding the transmission of EV71 throughout Southeast Asia is vital to inform preventative public health measures and vaccine development efforts.

Quantifying evolutionary constraints on B-cell affinity maturation

The antibody repertoire of each individual is continuously updated by the evolutionary process of B-cell receptor (BCR) mutation and selection. It has recently become possible to gain detailed information concerning this process through high-throughput sequencing. Here, we develop modern statistical molecular evolution methods for the analysis of B-cell sequence data, and then apply them to a very deep short-read dataset of BCRs. We find that the substitution process is conserved across individuals but varies significantly across gene segments. We investigate selection on BCRs using a novel method that side-steps the difficulties encountered by previous work in differentiating between selection and motif-driven mutation; this is done through stochastic mapping and empirical Bayes estimators that compare the evolution of in-frame and out-of-frame rearrangements. We use this new method to derive a per-residue map of selection, which provides a more nuanced view of the constraints on framework and variable regions.

Conservation and Innovation of APOBEC3A Restriction Functions during Primate Evolution

LINE-1 (long interspersed element-1) retroelements are the only active autonomous endogenous retroelements in human genomes. Their retrotransposition activity has created close to 50% of the current human genome. Due to the apparent costs of this proliferation, host genomes have evolved multiple mechanisms to curb LINE-1 retrotransposition. Here, we investigate the evolution and function of the LINE-1 restriction factor APOBEC3A, a member of the APOBEC3 cytidine deaminase gene family. We find that APOBEC3A genes have evolved rapidly under diversifying selection in primates, suggesting changes in APOBEC3A have been recurrently selected in a host-pathogen "arms race." Nonetheless, in contrast to previous reports, we find that the LINE-1 restriction activity of APOBEC3A proteins has been strictly conserved throughout simian primate evolution in spite of its pervasive diversifying selection. Based on these results, we conclude that LINE-1s have not driven the rapid evolution of APOBEC3A in primates. In contrast to this conserved LINE-1 restriction, we find that a subset of primate APOBEC3A genes have enhanced antiviral restriction. We trace this gain of antiviral restriction in APOBEC3A to the common ancestor of a subset of Old World monkeys. Thus, APOBEC3A has not only maintained its LINE-1 restriction ability, but also evolved a gain of antiviral specificity against other pathogens. Our findings suggest that while APOBEC3A has evolved to restrict additional pathogens, only those adaptive amino acid changes that leave LINE-1 restriction unperturbed have been tolerated.

Molecular characterization and prevalence of two capulaviruses: Alfalfa leaf curl virus from France and Euphorbia caput-medusae latent virus from South Africa

Little is known about the prevalence, diversity, evolutionary processes, genomic structures and population dynamics of viruses in the divergent geminivirus lineage known as the capulaviruses. We determined and analyzed full genome sequences of 13 Euphorbia caput-medusae latent virus (EcmLV) and 26 Alfalfa leaf curl virus (ALCV) isolates, and partial genome sequences of 23 EcmLV and 37 ALCV isolates. While EcmLV was asymptomatic in uncultivated southern African Euphorbia caput-medusae, severe alfalfa disease symptoms were associated with ALCV in southern France. The prevalence of both viruses exceeded 10% in their respective hosts. Besides using patterns of detectable negative selection to identify ORFs that are probably functionally expressed, we show that ALCV and EcmLV both display evidence of inter-species recombination and biologically functional genomic secondary structures. Finally, we show that whereas the EcmLV populations likely experience restricted geographical dispersion, ALCV is probably freely moving across the French Mediterranean region.

Low Major Histocompatibility Complex Class II Variation in the Endangered Indo-Pacific Humpback Dolphin (Sousa chinensis): Inferences About the Role of Balancing Selection

It has been widely reported that the major histocompatibility complex (MHC) is under balancing selection due to its immune function across terrestrial and aquatic mammals. The comprehensive studies at MHC and other neutral loci could give us a synthetic evaluation about the major force determining genetic diversity of species. Previously, a low level of genetic diversity has been reported among the Indo-Pacific humpback dolphin (Sousa chinensis) in the Pearl River Estuary (PRE) using both mitochondrial marker and microsatellite loci. Here, the expression and sequence polymorphism of 2 MHC class II genes (DQB and DRB) in 32 S. chinensis from PRE collected between 2003 and 2011 were investigated. High ratios of non-synonymous to synonymous substitution rates, codon-based selection analysis, and trans-species polymorphism (TSP) support the hypothesis that balancing selection acted on S. chinensis MHC sequences. However, only 2 haplotypes were detected at either DQB or DRB loci. Moreover, the lack of deviation from the Hardy-Weinberg expectation at DRB locus combined with the relatively low heterozygosity at both DQB locus and microsatellite loci suggested that balancing selection might not be sufficient, which further suggested that genetic drift associated with historical bottlenecks was not mitigated by balancing selection in terms of the loss of MHC and neutral variation in S. chinensis. The combined results highlighted the importance of maintaining the genetic diversity of the endangered S. chinensis.

Ebolavirus evolves in human to minimize the detection by immune cells by accumulating adaptive mutations

The current outbreak of Zaire ebolavirus (EBOV) lasted longer than the previous outbreaks and there is as yet no proven treatment or vaccine available. Understanding host immune pressure and associated EBOV immune evasion that drive the evolution of EBOV is vital for diagnosis as well as designing a highly effective vaccine. The aim of this study was to deduce adaptive selection pressure acting on each amino acid sites of EBOV responsible for the recent 2014 outbreak. Multiple statistical methods employed in the study include SLAC, FEL, REL, IFEL, FUBAR and MEME. Results show that a total of 11 amino acid sites from sGP and ssGP, and 14 sites from NP, VP40, VP24 and L proteins were inferred as positively and negatively selected, respectively. Overall, the function of 11 out of 25 amino acid sites under selection pressure exactly found to be involved in T cell and B-cell epitopes. We identified that the EBOV had evolved through purifying selection pressure, which is a predictor that is known to aid the virus to adapt better to the human host and subsequently reduce the efficiency of existing immunity. Furthermore, computational RNA structure prediction showed that the three synonymous nucleotide mutations in NP gene altered the RNA secondary structure and optimal base-pairing energy, implicating a possible effect on genome replication. Here, we have provided evidence that the EBOV strains involved in the recent 2014 outbreak have evolved to minimize the detection by T and B cells by accumulating adaptive mutations to increase the survival fitness.

Pleistocene divergence across a mountain range and the influence of selection on mitogenome evolution in threatened Australian freshwater cod species

Climatic differences across a taxon's range may be associated with specific bioenergetic demands and may result in genetics-based metabolic adaptation, particularly in aquatic ectothermic organisms that rely on heat exchange with the environment to regulate key physiological processes. Extending down the east coast of Australia, the Great Dividing Range (GDR) has a strong influence on climate and the evolutionary history of freshwater fish species. Despite the GDR acting as a strong contemporary barrier to fish movement, many species, and species with shared ancestries, are found on both sides of the GDR, indicative of historical dispersal events. We sequenced complete mitogenomes from the four extant species of the freshwater cod genus Maccullochella, two of which occur on the semi-arid, inland side of the GDR, and two on the mesic coastal side. We constructed a dated phylogeny and explored the relative influences of purifying and positive selection in the evolution of mitogenome divergence among species. Results supported mid- to late-Pleistocene divergence of Maccullochella across the GDR (220-710 thousand years ago), bringing forward previously reported dates. Against a background of pervasive purifying selection, we detected potentially functionally relevant fixed amino acid differences across the GDR. Although many amino acid differences between inland and coastal species may have become fixed under relaxed purifying selection in coastal environments rather than positive selection, there was evidence of episodic positive selection acting on specific codons in the Mary River coastal lineage, which has consistently experienced the warmest and least extreme climate in the genus.

Visualization of Host-Polerovirus Interaction Topologies Using Protein Interaction Reporter Technology

Demonstrating direct interactions between host and virus proteins during infection is a major goal and challenge for the field of virology. Most protein interactions are not binary or easily amenable to structural determination. Using infectious preparations of a polerovirus (Potato leafroll virus [PLRV]) and protein interaction reporter (PIR), a revolutionary technology that couples a mass spectrometric-cleavable chemical cross-linker with high-resolution mass spectrometry, we provide the first report of a host-pathogen protein interaction network that includes data-derived, topological features for every cross-linked site that was identified. We show that PLRV virions have hot spots of protein interaction and multifunctional surface topologies, revealing how these plant viruses maximize their use of binding interfaces. Modeling data, guided by cross-linking constraints, suggest asymmetric packing of the major capsid protein in the virion, which supports previous epitope mapping studies. Protein interaction topologies are conserved with other species in the Luteoviridae and with unrelated viruses in the Herpesviridae and Adenoviridae. Functional analysis of three PLRV-interacting host proteins in planta using a reverse-genetics approach revealed a complex, molecular tug-of-war between host and virus. Structural mimicry and diversifying selection-hallmarks of host-pathogen interactions-were identified within host and viral binding interfaces predicted by our models. These results illuminate the functional diversity of the PLRV-host protein interaction network and demonstrate the usefulness of PIR technology for precision mapping of functional host-pathogen protein interaction topologies.

IMPORTANCE

The exterior shape of a plant virus and its interacting host and insect vector proteins determine whether a virus will be transmitted by an insect or infect a specific host. Gaining this information is difficult and requires years of experimentation. We used protein interaction reporter (PIR) technology to illustrate how viruses exploit host proteins during plant infection. PIR technology enabled our team to precisely describe the sites of functional virus-virus, virus-host, and host-host protein interactions using a mass spectrometry analysis that takes just a few hours. Applications of PIR technology in host-pathogen interactions will enable researchers studying recalcitrant pathogens, such as animal pathogens where host proteins are incorporated directly into the infectious agents, to investigate how proteins interact during infection and transmission as well as develop new tools for interdiction and therapy.

Population-Level Immune-Mediated Adaptation in HIV-1 Polymerase during the North American Epidemic

Human leukocyte antigen (HLA) class I-associated polymorphisms in HIV-1 that persist upon transmission to HLA-mismatched hosts may spread in the population as the epidemic progresses. Transmission of HIV-1 sequences containing such adaptations may undermine cellular immune responses to the incoming virus in future hosts. Building upon previous work, we investigated the extent of HLA-associated polymorphism accumulation in HIV-1 polymerase (Pol) through comparative analysis of linked HIV-1/HLA class I genotypes sampled during historic (1979 to 1989; n = 338) and modern (2001 to 2011; n = 278) eras from across North America (Vancouver, BC, Canada; Boston, MA; New York, NY; and San Francisco, CA). Phylogenies inferred from historic and modern HIV-1 Pol sequences were star-like in shape, with an inferred most recent common ancestor (epidemic founder virus) sequence nearly identical to the modern North American subtype B consensus sequence. Nevertheless, modern HIV-1 Pol sequences exhibited roughly 2-fold-higher patristic (tip-to-tip) genetic distances than historic sequences, with HLA pressures likely driving ongoing diversification. Moreover, the frequencies of published HLA-associated polymorphisms in individuals lacking the selecting HLA class I allele was on average ∼2.5-fold higher in the modern than in the historic era, supporting their spread in circulation, though some remained stable in frequency during this time. Notably, polymorphisms restricted by protective HLA alleles appear to be spreading to a greater relative extent than others, though these increases are generally of modest absolute magnitude. However, despite evidence of polymorphism spread, North American hosts generally remain at relatively low risk of acquiring an HIV-1 polymerase sequence substantially preadapted to their HLA profiles, even in the present era.

IMPORTANCE

HLA class I-restricted cytotoxic T-lymphocyte (CTL) escape mutations in HIV-1 that persist upon transmission may accumulate in circulation over time, potentially undermining host antiviral immunity to the transmitted viral strain. We studied >600 experimentally collected HIV-1 polymerase sequences linked to host HLA information dating back to 1979, along with phylogenetically reconstructed HIV-1 sequences dating back to the virus' introduction into North America. Overall, our results support the gradual spread of many-though not all-HIV-1 polymerase immune escape mutations in circulation over time. This is consistent with recent observations from other global regions, though the extent of polymorphism accumulation in North America appears to be lower than in populations with high seroprevalence, older epidemics, and/or limited HLA diversity. Importantly, the risk of acquiring an HIV-1 polymerase sequence at transmission that is substantially preadapted to one's HLA profile remains relatively low in North America, even in the present era.

Geographically driven adaptation of chilli veinal mottle virus revealed by genetic diversity analysis of the coat protein gene

Chilli veinal mottle virus (ChiVMV) is an important plant pathogen with a wide host range. The genetic structure of ChiVMV was investigated by analyzing the coat protein (CP) genes of 87 ChiVMV isolates from seven Asian regions. Pairwise F ST values between ChiVMV populations ranged from 0.108 to 0.681, indicating a significant spatial structure for this pathogen. In phylogeny-trait association analysis, the viral isolates from the same region tended to group together, showing a distinct geographic feature. These results suggest that geographic driven adaptation may be an important determinant of the genetic diversity of ChiVMV.

The role of defensive ecological interactions in the evolution of conotoxins

Venoms comprise of complex mixtures of peptides evolved for predation and defensive purposes. Remarkably, some carnivorous cone snails can inject two distinct venoms in response to predatory or defensive stimuli, providing a unique opportunity to study separately how different ecological pressures contribute to toxin diversification. Here, we report the extraordinary defensive strategy of the Rhizoconus subgenus of cone snails. The defensive venom from this worm-hunting subgenus is unusually simple, almost exclusively composed of αD-conotoxins instead of the ubiquitous αA-conotoxins found in the more complex defensive venom of mollusc- and fish-hunting cone snails. A similarly compartmentalized venom gland as those observed in the other dietary groups facilitates the deployment of this defensive venom. Transcriptomic analysis of a Conus vexillum venom gland revealed the αD-conotoxins as the major transcripts, with lower amounts of 15 known and four new conotoxin superfamilies also detected with likely roles in prey capture. Our phylogenetic and molecular evolution analysis of the αD-conotoxins from five subgenera of cone snails suggests they evolved episodically as part of a defensive strategy in the Rhizoconus subgenus. Thus, our results demonstrate an important role for defence in the evolution of conotoxins.

Causes of evolutionary rate variation among protein sites

It has long been recognized that certain sites within a protein, such as sites in the protein core or catalytic residues in enzymes, are evolutionarily more conserved than other sites. However, our understanding of rate variation among sites remains surprisingly limited. Recent progress to address this includes the development of a wide array of reliable methods to estimate site-specific substitution rates from sequence alignments. In addition, several molecular traits have been identified that correlate with site-specific mutation rates, and novel mechanistic biophysical models have been proposed to explain the observed correlations. Nonetheless, current models explain, at best, approximately 60% of the observed variance, highlighting the limitations of current methods and models and the need for new research directions.

Comparative sequence analyses of rhodopsin and RPE65 reveal patterns of selective constraint across hereditary retinal disease mutations

Retinitis pigmentosa (RP) comprises several heritable diseases that involve photoreceptor, and ultimately retinal, degeneration. Currently, mutations in over 50 genes have known links to RP. Despite advances in clinical characterization, molecular characterization of RP remains challenging due to the heterogeneous nature of causal genes, mutations, and clinical phenotypes. In this study, we compiled large datasets of two important visual genes associated with RP: rhodopsin, which initiates the phototransduction cascade, and the retinoid isomerase RPE65, which regenerates the visual cycle. We used a comparative evolutionary approach to investigate the relationship between interspecific sequence variation and pathogenic mutations that lead to degenerative retinal disease. Using codon-based likelihood methods, we estimated evolutionary rates (dN/dS) across both genes in a phylogenetic context to investigate differences between pathogenic and nonpathogenic amino acid sites. In both genes, disease-associated sites showed significantly lower evolutionary rates compared to nondisease sites, and were more likely to occur in functionally critical areas of the proteins. The nature of the dataset (e.g., vertebrate or mammalian sequences), as well as selection of pathogenic sites, affected the differences observed between pathogenic and nonpathogenic sites. Our results illustrate that these methods can serve as an intermediate step in understanding protein structure and function in a clinical context, particularly in predicting the relative pathogenicity (i.e., functional impact) of point mutations and their downstream phenotypic effects. Extensions of this approach may also contribute to current methods for predicting the deleterious effects of candidate mutations and to the identification of protein regions under strong constraint where we expect pathogenic mutations to occur.

Testis-enriched heat shock protein A2 (HSPA2): Adaptive advantages of the birds with internal testes over the mammals with testicular descent

The molecular chaperone heat shock protein A2 (HSPA2), a member of the 70 kDa heat shock protein (HSP70) family, plays an important role in spermatogenesis and male fertility. Although HSPA2 is evolutionarily highly conserved across the metazoan lineages, the observation of striking differences in temperature-sensitive expressions, testicular physiology, spermatogenesis, as well as its role in male fertility indicates that avian and mammalian HSPA2 may exhibit distinct evolutionary trajectory. The present study reports that while mammalian HSPA2 is constrained by intense purifying selection, avian HSPA2 has been subjected to positive selection. The majority of the positively selected amino acid residues fall on the α-helix and β-sheets of the peptide-binding domain located at the carboxyl-terminal region of the avian HSPA2. The detection of positively selected sites at the helix and β-sheets, which are less tolerant to molecular adaptation, indicates an important functional consequence and contribution to the structural and functional diversification of the avian HSPA2. Collectively, avian HSPA2 may have an adaptive advantage over the mammals in response to heat stress, and therefore, mammals with testicular descent may be at a greater risk in the event of scrotal temperature rise.

Driving forces behind the evolution of the Aleutian mink disease parvovirus in the context of intensive farming

Aleutian mink disease virus (AMDV) causes plasmacytosis, an immune complex-associated syndrome that affects wild and farmed mink. The virus can also infect other small mammals (e.g., ferrets, skunks, ermines, and raccoons), but the disease in these hosts has been studied less. In 2007, a mink plasmacytosis outbreak began on the Island of Newfoundland, and the virus has been endemic in farms since then. In this study, we evaluated the molecular epidemiology of AMDV in farmed and wild animals of Newfoundland since before the beginning of the outbreak and investigated the epidemic in a global context by studying AMDV worldwide, thereby examining its diffusion and phylogeography. Furthermore, AMDV evolution was examined in the context of intensive farming, where host population dynamics strongly influence viral evolution. Partial NS1 sequences and several complete genomes were obtained from Newfoundland viruses and analyzed along with numerous sequences from other locations worldwide that were either obtained as part of this study or from public databases. We observed very high viral diversity within Newfoundland and within single farms, where high rates of co-infection, recombinant viruses and polymorphisms were observed within single infected individuals. Worldwide, we documented a partial geographic distribution of strains, where viruses from different countries co-exist within clades but form country-specific subclades. Finally, we observed the occurrence of recombination and the predominance of negative selection pressure on AMDV proteins. A surprisingly low number of immunoepitopic sites were under diversifying pressure, possibly because AMDV gains no benefit by escaping the immune response as viral entry into target cells is mediated through interactions with antibodies, which therefore contribute to cell infection. In conclusion, the high prevalence of AMDV in farms facilitates the establishment of co-infections that can favor the occurrence of recombination and enhance viral diversity. Viruses are then exchanged between different farms and countries and can be introduced into the wild, with the rapidly evolving viruses producing many parallel lineages.

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.

Evolution and Spread of Ebola Virus in Liberia, 2014-2015

The 2013-present Western African Ebola virus disease (EVD) outbreak is the largest ever recorded with >28,000 reported cases. Ebola virus (EBOV) genome sequencing has played an important role throughout this outbreak; however, relatively few sequences have been determined from patients in Liberia, the second worst-affected country. Here, we report 140 EBOV genome sequences from the second wave of the Liberian outbreak and analyze them in combination with 782 previously published sequences from throughout the Western African outbreak. While multiple early introductions of EBOV to Liberia are evident, the majority of Liberian EVD cases are consistent with a single introduction, followed by spread and diversification within the country. Movement of the virus within Liberia was widespread, and reintroductions from Liberia served as an important source for the continuation of the already ongoing EVD outbreak in Guinea. Overall, little evidence was found for incremental adaptation of EBOV to the human host.

Episodic Diversifying Selection Shaped the Genomes of Gibbon Ape Leukemia Virus and Related Gammaretroviruses

Gibbon ape leukemia viruses (GALVs) are part of a larger group of pathogenic gammaretroviruses present across phylogenetically diverse host species of Australasian mammals. Despite the biomedical utility of GALVs as viral vectors and in cancer gene therapy, full genome sequences have not been determined for all of the five identified GALV strains, nor has a comprehensive evolutionary analysis been performed. We therefore generated complete genomic sequences for each GALV strain using hybridization capture and high-throughput sequencing. The four strains of GALV isolated from gibbons formed a monophyletic clade that was closely related to the woolly monkey virus (WMV), which is a GALV strain that likely originated in a gibbon host. The GALV-WMV clade in turn formed a sister group to the koala retroviruses (KoRVs). Genomic signatures of episodic diversifying selection were detected among the gammaretroviruses with concentration in the env gene across the GALV strains that were particularly oncogenic and KoRV strains that were potentially exogenous, likely reflecting their adaptation to the host immune system. In vitro studies involving vectors chimeric between GALV and KoRV-B established that variable regions A and B of the surface unit of the envelope determine which receptor is used by a viral strain to enter host cells.

IMPORTANCE

The gibbon ape leukemia viruses (GALVs) are among the most medically relevant retroviruses due to their use as viral vectors for gene transfer and in cancer gene therapy. Despite their importance, full genome sequences have not been determined for the majority of primate isolates, nor has comprehensive evolutionary analysis been performed, despite evidence that the viruses are facing complex selective pressures associated with cross-species transmission. Using hybridization capture and high-throughput sequencing, we report here the full genome sequences of all the GALV strains and demonstrate that diversifying selection is acting on them, particularly in the envelope gene in functionally important domains, suggesting that host immune pressure is shaping GALV evolution.

Comparative analysis of complete mitochondrial genomes suggests that relaxed purifying selection is driving high nonsynonymous evolutionary rate of the NADH2 gene in whitefish (Coregonus ssp.)

Several studies have recently reported evidence for positive selection acting on the mitochondrial genome (mitogenome), emphasizing its potential role in adaptive divergence and speciation. In this study we searched 107 full mitogenomes of recently diverged species and lineages of whitefish (Coregonus ssp.) for signals of positive selection. These salmonids show several distinct morphological and ecological differences that may be associated with energetics and therefore potentially positive selection at the mitogenome level. We found that purifying selection and genetic drift were the predominant evolutionary forces acting on the analyzed mitogenomes. However, the NADH dehydrogenase 2 gene (ND2) showed a highly elevated dN/dS ratio compared to the other mitochondrial genes, which was significantly higher in whitefish compared to other salmonids. We therefore further examined nonsynonymous evolution in ND2 by (i) mapping amino acid changes to a protein model structure which showed that they were located away from key functional residues of the protein, (ii) locating them in the sequences of other species of fish (Salmonidae, Anguillidae, Scombridae and Percidae) only to find pronounced overlap of nonsynonymous regions. We thus conclude that relaxed purifying selection is driving the evolution of ND2 by affecting mostly regions that have lower functional relevance.

Molecular Evolution of the TET Gene Family in Mammals

Ten-eleven translocation (TET) proteins, a family of Fe²⁺- and 2-oxoglutarate-dependent dioxygenases, are involved in DNA demethylation. They also help regulate various cellular functions. Three TET paralogs have been identified (TET1, TET2, and TET3) in humans. This study focuses on the evolution of mammalian TET genes. Distinct patterns in TET1 and TET2 vs. TET3 were revealed by codon-based tests of positive selection. Results indicate that TET1 and TET2 genes have experienced positive selection more frequently than TET3 gene, and that the majority of codon sites evolved under strong negative selection. These findings imply that the selective pressure on TET3 may have been relaxed in several lineages during the course of evolution. Our analysis of convergent amino acid substitutions also supports the different evolutionary dynamics among TET gene subfamily members. All of the five amino acid sites that are inferred to have evolved under positive selection in the catalytic domain of TET2 are localized at the protein’s outer surface. The adaptive changes of these positively selected amino acid sites could be associated with dynamic interactions between other TET-interacting proteins, and positive selection thus appears to shift the regulatory scheme of TET enzyme function.

Genome-Wide Analysis of Evolutionary Markers of Human Influenza A(H1N1)pdm09 and A(H3N2) Viruses May Guide Selection of Vaccine Strain Candidates

Here we analyzed whole-genome sequences of 3,969 influenza A(H1N1)pdm09 and 4,774 A(H3N2) strains that circulated during 2009–2015 in the world. The analysis revealed changes at 481 and 533 amino acid sites in proteins of influenza A(H1N1)pdm09 and A(H3N2) strains, respectively. Many of these changes were introduced as a result of random drift. However, there were 61 and 68 changes that were present in relatively large number of A(H1N1)pdm09 and A(H3N2) strains, respectively, that circulated during relatively long time. We named these amino acid substitutions evolutionary markers, as they seemed to contain valuable information regarding the viral evolution. Interestingly, influenza A(H1N1)pdm09 and A(H3N2) viruses acquired non-overlapping sets of evolutionary markers. We next analyzed these characteristic markers in vaccine strains recommended by the World Health Organization for the past five years. Our analysis revealed that vaccine strains carried only few evolutionary markers at antigenic sites of viral hemagglutinin (HA) and neuraminidase (NA). The absence of these markers at antigenic sites could affect the recognition of HA and NA by human antibodies generated in response to vaccinations. This could, in part, explain moderate efficacy of influenza vaccines during 2009–2014. Finally, we identified influenza A(H1N1)pdm09 and A(H3N2) strains, which contain all the evolutionary markers of influenza A strains circulated in 2015, and which could be used as vaccine candidates for the 2015/2016 season. Thus, genome-wide analysis of evolutionary markers of influenza A(H1N1)pdm09 and A(H3N2) viruses may guide selection of vaccine strain candidates.

Deep venomics of the Pseudonaja genus reveals inter- and intra-specific variation

Australian elapid venom remains an under-investigated resource of novel bioactive peptides. In this study, the venom gland transcriptomes and proteomes of the Australian western brown snakes, Pseudonaja aspidorhyncha and Pseudonaja nuchalis, were compared to Pseudonaja textilis. A deep venomics strategy incorporating high throughput 454 pyrosequencing gave a total of 200,911 raw reads for the three venoms. Subsequent annotation identified 5716 transcripts from 20 different toxin families with inter-specific variation between species observed in eight of the less abundant families. Integration of each venom proteome with the corresponding annotated reads identified 65 isoforms from six toxin families; high sequence coverage highlighted subtle differences between sequences and intra and inter-specific variation between species. High quality MS/MS data identified unusual glycoforms with natriuretic peptides from P. aspidorhyncha and P. nuchaliscontaining O-linked trisaccharides with high homology to the glycosylated region of TNPc. Molecular evolutionary assessments indicated the accelerated evolution of all toxin families with the exception of both natriuretic peptides and P. aspidorhyncha PLA2s that were found to be evolutionarily constrained under purifying selection pressures. This study has revealed a wide range of novel peptide sequences from six bioactive peptide families and highlights the subtle differences between toxins in these closely related species.

BIOLOGICAL SIGNIFICANCE

Mining Australia's vastly untapped source of toxins from its venomous creatures has been significantly advanced by employing deep venomics methodology. Technological advances in transcriptome analysis using next generation sequencing platforms and proteome analysis by highly sensitive tandem mass spectrometry allowed a more comprehensive interrogation of three underinvestigated brown snake (Pseudonaja) venoms uncovering many novel peptide sequences that are unique to these closely related species. This generic strategy will provide invaluable information when applied to other venomous snakes for a deeper understanding of venom composition, envenomation, venom evolution, as well as identifying research tools and drug leads.