Detecting individual sites subject to episodic diversifying selection

Evolution of bovine ephemeral fever virus in the Australian episystem

Bovine ephemeral fever virus (BEFV) is an arthropod-borne rhabdovirus that causes a debilitating disease of cattle in Africa, Asia, and Australia; however, its global geodynamics are poorly understood. An evolutionary analysis of G gene (envelope glycoprotein) ectodomain sequences of 97 BEFV isolates collected from Australia during 1956 to 2012 revealed that all have a single common ancestor and are phylogenetically distinct from BEFV sampled in other geographical regions. The age of the Australian clade is estimated to be between 56 and 65 years, suggesting that BEFV has entered the continent on few occasions since it was first reported in 1936 and that the 1955-1956 epizootic was the source of all currently circulating viruses. Notably, the Australian clade has evolved as a single genetic lineage across the continent and at a high evolutionary rate of ∼10(-3) nucleotide substitutions/site/year. Screening of 66 isolates using monoclonal antibodies indicated that neutralizing antigenic sites G1, G2, and G4 have been relatively stable, although variations in site G3a/b defined four antigenic subtypes. A shift in an epitope at site G3a, which occurred in the mid-1970s, was strongly associated with a K218R substitution. Similarly, a shift at site G3b was associated primarily with substitutions at residues 215, 220, and 223, which map to the tip of the spike on the prefusion form of the G protein. Finally, we propose that positive selection on residue 215 was due to cross-reacting neutralizing antibody to Kimberley virus (KIMV). This is the first study of the evolution of BEFV in Australia, showing that the virus has entered the continent only once during the past 50 to 60 years, it is evolving at a relatively constant rate as a single genetic lineage, and although the virus is relatively stable antigenically, mutations have resulted in four antigenic subtypes. Furthermore, the study shows that the evolution of BEFV in Australia appears to be driven, at least in part, by cross-reactive antibodies to KIMV which has a similar distribution and ecology but has not been associated with disease. As BEFV and KIMV are each known to be present in Africa and Asia, this interaction may occur on a broader geographic scale.

Identification and diversity of killer cell Ig-like receptors in Aotus vociferans, a New World monkey

Previous BAC clone analysis of the Platyrrhini owl monkey KIRs have shown an unusual genetic structure in some loci. Therefore, cDNAs encoding KIR molecules from eleven Aotus vociferans monkeys were characterized here; ten putative KIR loci were found, some of which encoded atypical proteins such as KIR4DL and transcripts predicted to encode a D0+D1 configuration (AOTVOKIR2DL1*01v1) which appear to be unique in the Aotus genus. Furthermore, alternative splicing was found as a likely mechanism for producing activator receptors in A. vociferans species. KIR proteins from New World monkeys may be split into three new lineages according to domain by domain phylogenetic analysis. Although the A. vociferans KIR family displayed a high divergence among paralogous genes, individual loci were limited in their genetic polymorphism. Selection analysis showed that both constrained and rapid evolution may operate within the AvKIR family. The frequent alternative splicing (as a likely mechanism generating activator receptors), the presence of KIR4DL and KIR2DL1 (D0+D1) molecules and other data reported here suggest that the KIR family in Aotus has had a rapid evolution, independent from its Catarrhini counterparts.

Atractaspis aterrima toxins: the first insight into the molecular evolution of venom in side-stabbers

Although snake venoms have been the subject of intense research, primarily because of their tremendous potential as a bioresource for design and development of therapeutic compounds, some specific groups of snakes, such as the genus Atractaspis, have been completely neglected. To date only limited number of toxins, such as sarafotoxins have been well characterized from this lineage. In order to investigate the molecular diversity of venom from Atractaspis aterrima—the slender burrowing asp, we utilized a high-throughput transcriptomic approach completed with an original bioinformatics analysis pipeline. Surprisingly, we found that Sarafotoxins do not constitute the major ingredient of the transcriptomic cocktail; rather a large number of previously well-characterized snake venom-components were identified. Notably, we recovered a large diversity of three-finger toxins (3FTxs), which were found to have evolved under the significant influence of positive selection. From the normalized and non-normalized transcriptome libraries, we were able to evaluate the relative abundance of the different toxin groups, uncover rare transcripts, and gain new insight into the transcriptomic machinery. In addition to previously characterized toxin families, we were able to detect numerous highly-transcribed compounds that possess all the key features of venom-components and may constitute new classes of toxins.

Phylogenetic diversity of GB virus C at the antigenic site of E2 protein

Amino acids at position 267-298 in E2 protein of GB virus C (GBV-C) were recognized as the antigenic site, and peptides within the region were previously reported to have inhibitory effect on HIV entry. The effect of sequence variability between different types of GBV-C on the antigenic region of the E2 protein was studied by using phylogenetic analysis. Eighty-one unique sequences encompassing this region derived from all seven GBV-C genotypes were compared to each other in this study. The results showed that GBV-C E2 antigenic nucleotide sites belonging to genotype 3 clustered together regardless of synonymous or nonsynonmous sites in the region, whereas, GBV-C E2 antigenic nucleotide sites belonging to the other 6 genotypes clustered together regardless of genotypes. Despite the fact that GBV-C genotypes might confer different degree of 'protection' against HIV, the lack of clustering as a unique group based on the amino acid differences in the antigenic site among the six genotypes suggested some other genomic regions or secondary structure of E2 protein might have played a crucial role in determining the variable protection effect of GBV-C on HIV infection.

Performance of standard and stochastic branch-site models for detecting positive selection among coding sequences

The branch-site model is a widely popular approach that accommodates for the lineage- and the site-specific heterogeneity of natural selection regimes among coding sequences. This model relies on prior knowledge of the (foreground) lineage(s) evolving under positive selection at some sites. Unfortunately, such prior information is not always available in practice. A more recent technique (Guindon S, Rodrigo A, Dyer K, Huelsenbeck J. 2004. Modeling the site-specific variation of selection patterns along lineages. Proc Natl Acad Sci USA 101:12957-12962) alleviates this issue by explicitly modeling the variability of selection patterns using a stochastic process. However, the performance of this approach for deciding whether a set of homologous sequences evolved under positive selection at some point has not been assessed yet. This study compares the sensitivity and specificity of tests for positive selection derived from both the standard and the stochastic approaches using extensive simulations. We show that the two methods have low proportions of type I errors, that is, they tend to be conservative when testing the null hypothesis of no positive selection if sequences truly evolve under neutral or negative selection regimes. Also, the standard approach is more powerful than the stochastic one when the prior knowledge on foreground lineages is correct. When this prior is incorrect, however, the stochastic approach outperforms the standard model in a broad range of conditions. Additional comparisons also suggest that the stochastic branch-site method compares favorably with the recently proposed mixed-effects model of evolution of Murrell et al. (Murrell B, Wertheim JO, Moola S, Weighill T, Scheffler K, Pond SLK. 2012. Detecting individual sites subject to episodic diversifying selection. PLoS Genet. 8:e1002764). Altogether, our results show that the standard branch-site model is well suited to confirmatory analyses, whereas the stochastic approach should be preferred over the standard or the mixed-effects ones for exploratory studies.

The genesis and source of the H7N9 influenza viruses causing human infections in China

A novel H7N9 influenza A virus first detected in March 2013 has since caused more than 130 human infections in China, resulting in 40 deaths. Preliminary analyses suggest that the virus is a reassortant of H7, N9 and H9N2 avian influenza viruses, and carries some amino acids associated with mammalian receptor binding, raising concerns of a new pandemic. However, neither the source populations of the H7N9 outbreak lineage nor the conditions for its genesis are fully known. Using a combination of active surveillance, screening of virus archives, and evolutionary analyses, here we show that H7 viruses probably transferred from domestic duck to chicken populations in China on at least two independent occasions. We show that the H7 viruses subsequently reassorted with enzootic H9N2 viruses to generate the H7N9 outbreak lineage, and a related previously unrecognized H7N7 lineage. The H7N9 outbreak lineage has spread over a large geographic region and is prevalent in chickens at live poultry markets, which are thought to be the immediate source of human infections. Whether the H7N9 outbreak lineage has, or will, become enzootic in China and neighbouring regions requires further investigation. The discovery here of a related H7N7 influenza virus in chickens that has the ability to infect mammals experimentally, suggests that H7 viruses may pose threats beyond the current outbreak. The continuing prevalence of H7 viruses in poultry could lead to the generation of highly pathogenic variants and further sporadic human infections, with a continued risk of the virus acquiring human-to-human transmissibility.

Dynamic changes in viral population structure and compartmentalization during chronic hepatitis C virus infection in children

Classic phylogenetic and modern population-based clustering methods were used to analyze hepatitis C virus (HCV) evolution in plasma and to assess viral compartmentalization within peripheral blood mononuclear cells (PBMCs) in 6 children during 3.2-9.6yr of follow-up. Population structure analysis of cloned amplicons encompassing hypervariable region 1 led to the distinction of two evolutionary patterns, one highly divergent and another one genetically homogeneous. Viral adaptability was reflected by co-evolution of viral communities switching rapidly from one to another in the context of divergence and stability associated with highly homogeneous communities which were replaced by new ones after long periods. Additionally, viral compartmentalization of HCV in PBMCs was statistically demonstrated, suggesting their role as a pool of genetic variability. Our results support the idea of a community-based structure of HCV viral populations during chronic infection and highlight a role of the PBMC compartment in the persistence of such structure.

Phylogenetic analyses of vector mosquito basic helix-loop-helix transcription factors

Basic helix-loop-helix (bHLH) transcription factors play critical roles in the regulation of a wide range of developmental processes in higher organisms and have been identified in more than 20 organisms. Mosquitoes are important vectors of certain human diseases. In this study, Aedes aegypti, Anopheles gambiae str. PEST and Culex quinquefasciatus genomes were found to encode 55, 55 and 57 bHLH genes, respectively. Further phylogenetic analyses and OrthoDB and Kyoto encyclopedia of genes and genomes orthology database searches led us to define orthology for all the identified mosquito bHLHs successfully. This provides useful information with which to update annotations to 40 Ae. aegypti, 55 An. gambiae and 38 C. quinquefasciatus bHLH genes in VectorBase. The mosquito lineage has more bHLH genes in the Atonal, neurogenin (Ngn) and Hes-related with YRPW motif (Hey) families than do other insect species, suggesting that mosquitoes have evolved to be more sensitive to vibration, light and chemicals. Mosquito bHLH genes generally have higher evolutionary rates than other insect species. However, no pervasive positive selection occurred in the evolution of insect bHLH genes. Only episodic positive selection was found to affect evolution of bHLH genes in 11 families. Besides, coding regions of several Ae. aegypti bHLH motifs have unusually long introns in which multiple copies of transposable elements have been identified. These data provide a solid basis for further studies on structures and functions of bHLH proteins in the regulation of mosquito development and for prevention and control of mosquito-mediated human diseases.

Congruence between mcy based genetic type and microcystin composition within the populations of toxic Microcystis in a Plateau Lake, China

Genetic diversity and differential microcystin contributions within the populations of toxic Microcystis in freshwater ecosystems have not been fully discussed. To illustrate this issue, we sequenced clone libraries of the first adenylation domain of mcyB (mcyBA1) from Microcystis populations in a plateau lake. Phylogenetic analysis revealed two divergent groups of McyBA1 with two closely related subgroups within Group I. While neutral evolution was indicated on the whole McyBA1 domain, 1-5 recombination breakpoints and several codons under positive or negative selections were found. Significant seasonal changes of McyBA1 diversity were observed. Group I composed the major proportion of the McyBA1 pool throughout the growing season in the whole lake. Group IA and Group II denominated Microcystis strains isolated from this lake were characterized by preference production of microcystin-RR (62-85%) and microcystin-LR (> 98%) respectively. We detected the intracellular microcystins in lake water and microcystin-RR was a main variant (mostly > 50%). In summary, McyBA1 subgroups were dominant within the population of toxic Microcystis and contributed the predominance of microcystin-RR in the lake. The differences of substrate preference in microcystin biosynthesis among groups were caused by neutral evolution and homologous recombination.

Contrasted evolutionary histories of two Toll-like receptors (Tlr4 and Tlr7) in wild rodents (MURINAE)

BACKGROUND

In vertebrates, it has been repeatedly demonstrated that genes encoding proteins involved in pathogen-recognition by adaptive immunity (e.g. MHC) are subject to intensive diversifying selection. On the other hand, the role and the type of selection processes shaping the evolution of innate-immunity genes are currently far less clear. In this study we analysed the natural variation and the evolutionary processes acting on two genes involved in the innate-immunity recognition of Microbe-Associated Molecular Patterns (MAMPs).

RESULTS

We sequenced genes encoding Toll-like receptor 4 (Tlr4) and 7 (Tlr7), two of the key bacterial- and viral-sensing receptors of innate immunity, across 23 species within the subfamily Murinae. Although we have shown that the phylogeny of both Tlr genes is largely congruent with the phylogeny of rodents based on a comparably sized non-immune sequence dataset, we also identified several potentially important discrepancies. The sequence analyses revealed that major parts of both Tlrs are evolving under strong purifying selection, likely due to functional constraints. Yet, also several signatures of positive selection have been found in both genes, with more intense signal in the bacterial-sensing Tlr4 than in the viral-sensing Tlr7. 92% and 100% of sites evolving under positive selection in Tlr4 and Tlr7, respectively, were located in the extracellular domain. Directly in the Ligand-Binding Region (LBR) of TLR4 we identified two rapidly evolving amino acid residues and one site under positive selection, all three likely involved in species-specific recognition of lipopolysaccharide of gram-negative bacteria. In contrast, all putative sites of LBRTLR7 involved in the detection of viral nucleic acids were highly conserved across rodents. Interspecific differences in the predicted 3D-structure of the LBR of both Tlrs were not related to phylogenetic history, while analyses of protein charges clearly discriminated Rattini and Murini clades.

CONCLUSIONS

In consequence of the constraints given by the receptor protein function purifying selection has been a dominant force in evolution of Tlrs. Nevertheless, our results show that episodic diversifying parasite-mediated selection has shaped the present species-specific variability in rodent Tlrs. The intensity of diversifying selection was higher in Tlr4 than in Tlr7, presumably due to structural properties of their ligands.

Patterns of evolutionary selection pressure in the immune signaling protein TRAF3IP2 in mammals

TRAF3 interacting protein 2 (TRAF3IP2) is important for immune responses to pathogens, inflammatory signals and autoimmunity in mammals. In the present study, we collected 19 mammalian TRAF3IP2 sequences and investigated the various types of selection pressure acting on them. Maximum likelihood estimations of nonsynonymous (dN) to synonymous (dS) substitution (dN/dS) ratios for the aligned coding sequences indicated that, as a whole, TRAF3IP2 has been subject to purifying selection. However, the N-terminus of the protein has been subject to higher selection pressure than the C-terminal domain. While eight amino acid residues within the N-terminus appear to have evolved under positive selection, no evidence for such selection was found in the C-terminus. The positively selected residues, which fall outside the currently known functional sites within TRAF3IP2, may have novel functions. The different selection pressures acting on the N- and C-terminal regions are consistent with their protein structures: the C-terminal structure is an ordered structure, whereas the N-terminus is disordered. Taken together with the results of previous studies, it is plausible that positive selection on the N-terminus of TRAF3IP2 may have occurred by competitive coevolution between mammalian hosts and viruses.

Selectively maintained paleoviruses in Holarctic water fleas reveal an ancient origin for phleboviruses

The ecological model, Daphnia pulex (Cladocera: Daphniidae), is broadly distributed in Holarctic freshwater habitats and has been the subject of multidisciplinary study for over half a century, but never has a natural RNA virus infection been reported in daphnids. Here we report on a group of paleoviruses related to RNA dependent RNA polymerase in the genome of D. pulex. Phylogenetic analysis suggests that these paleoviruses are derived from a viral lineage within the genus Phlebovirus. Comparison of the genomic sequences flanking individual paleoviruses reveal that some are orthologous viral insertions having been present in the common ancestor of the D. pulex species complex, which is millions of years old. Still, we detected some sites that have the signature of purifying selection. In contrast, other paleoviruses in this group seem to be unique to specific host lineages and even contain undisrupted open reading frames, suggesting either more recent acquisition, or selective maintenance.

Pseudoplusia includens single nucleopolyhedrovirus: genetic diversity, phylogeny and hypervariability of the pif-2 gene

The soybean looper (Pseudoplusia includens Walker, 1857) has become a major pest of soybean crops in Brazil. In order to determine the genetic diversity and phylogeny of variants of Pseudoplusia includens single nucleopolyhedrovirus (PsinSNPV-IA to -IG), partial sequences of the genes lef-8, lef-9, pif-2, phr and polh were obtained following degenerate PCR and phylogenetic trees constructed using maximum parsimony and Bayesian methods. The aligned sequences showed polymorphisms among the isolates, where the pif-2 gene was by far the most variable and is predicted to be under positive selection. Furthermore, some of the pif-2 DNA sequence mutations are predicted to result in significant amino acid substitutions, possibly leading to changes in oral infectivity of this baculovirus. Cladistic analysis revealed two closely related monophyletic groups, one containing PsinNPV isolates IB, IC and ID and another containing isolates IA, IE, IF and IG. The phylogeny of PsinSNPV in relation to 56 other baculoviruses was also determined from the concatenated partial LEF-8, LEF-9, PIF-2 and POLH/GRAN deduced amino acid sequences, using maximum-parsimony and Bayesian methods. This analysis clearly places PsinSNPV with the Group II Alphabaculovirus, where PsinSNPV is most closely related to Chrysodeixis chalcites NPV and Trichoplusia ni SNPV.

Evolution and heterogeneity of multiple serotypes of Dengue virus in Pakistan, 2006-2011

Background

Even though dengue has been recognized as one of the major public health threats in Pakistan, the understanding of its molecular epidemiology is still limited. The genotypic diversity of Dengue virus (DENV) serotypes involved in dengue outbreaks since 2005 in Pakistan is not well studied. Here, we investigated the origin, diversity, genetic relationships and geographic distribution of DENV to understand virus evolution during the recent expansion of dengue in Pakistan.

Methods

The study included 200 sera obtained from dengue-suspected patients from 2006 to 2011. DENV infection was confirmed in 94 (47%) sera by a polymerase chain reaction assay. These included 36 (38.3%) DENV-2, 57 DENV-3 (60.6%) and 1 DENV-4 (1.1%) cases. Sequences of 13 whole genomes (6 DENV-2, 6 DENV-3 and 1 DENV-4) and 49 envelope genes (26 DENV-2, 22 DENV-3 and 1 DENV-4) were analysed to determine the origin, phylogeny, diversity and selection pressure during virus evolution.

Results

DENV-2, DENV-3 and DENV-4 in Pakistan from 2006 to 2011 shared 98.5-99.6% nucleotide and 99.3-99.9% amino acid similarity with those circulated in the Indian subcontinent during the last decade. Nevertheless, Pakistan DENV-2 and DENV-3 strains formed distinct clades characterized by amino acid signatures of NS2A-I116T + NS5-K861R and NS3-K590R + NS5-S895L respectively. Each clade consisted of a heterogenous virus population that circulated in Southern (2006–2009) and Northern Pakistan (2011).

Conclusions

DENV-2, DENV-3 and DENV-4 that circulated during 2006–2011 are likely to have first introduced via the southern route of Pakistan. Both DENV-2 and DENV-3 have undergone in-situ evolution to generate heterogenous populations, possibly driven by sustained local DENV transmission during 2006–2011 periods. While both DENV-2 and DENV-3 continued to circulate in Southern Pakistan until 2009, DENV-2 has spread in a Northern direction to establish in Punjab Province, which experienced a massive dengue outbreak in 2011.

Computational analyses of an evolutionary arms race between mammalian immunity mediated by immunoglobulin A and its subversion by bacterial pathogens

IgA is the predominant immunoglobulin isotype in mucosal tissues and external secretions, playing important roles both in defense against pathogens and in maintenance of commensal microbiota. Considering the complexity of its interactions with the surrounding environment, IgA is a likely target for diversifying or positive selection. To investigate this possibility, the action of natural selection on IgA was examined in depth with six different methods: CODEML from the PAML package and the SLAC, FEL, REL, MEME and FUBAR methods implemented in the Datamonkey webserver. In considering just primate IgA, these analyses show that diversifying selection targeted five positions of the Cα1 and Cα2 domains of IgA. Extending the analysis to include other mammals identified 18 positively selected sites: ten in Cα1, five in Cα2 and three in Cα3. All but one of these positions display variation in polarity and charge. Their structural locations suggest they indirectly influence the conformation of sites on IgA that are critical for interaction with host IgA receptors and also with proteins produced by mucosal pathogens that prevent their elimination by IgA-mediated effector mechanisms. Demonstrating the plasticity of IgA in the evolution of different groups of mammals, only two of the eighteen selected positions in all mammals are included in the five selected positions in primates. That IgA residues subject to positive selection impact sites targeted both by host receptors and subversive pathogen ligands highlights the evolutionary arms race playing out between mammals and pathogens, and further emphasizes the importance of IgA in protection against mucosal pathogens.

Rare variants in hypermutable genes underlie common morphology and growth traits in wild Saccharomyces paradoxus

Understanding the molecular basis of common traits is a primary challenge of modern genetics. One model holds that rare mutations in many genetic backgrounds may often phenocopy one another, together explaining the prevalence of the resulting trait in the population. For the vast majority of phenotypes, the role of rare variants and the evolutionary forces that underlie them are unknown. In this work, we use a population of Saccharomyces paradoxus yeast as a model system for the study of common trait variation. We observed an unusual, flocculation and invasive-growth phenotype in one-third of S. paradoxus strains, which were otherwise unrelated. In crosses with each strain in turn, these morphologies segregated as a recessive Mendelian phenotype, mapping either to IRA1 or to IRA2, yeast homologs of the hypermutable human neurofibromatosis gene NF1. The causal IRA1 and IRA2 haplotypes were of distinct evolutionary origin and, in addition to their morphological effects, associated with hundreds of stress-resistance and growth traits, both beneficial and disadvantageous, across S. paradoxus. Single-gene molecular genetic analyses confirmed variant IRA1 and IRA2 haplotypes as causal for these growth characteristics, many of which were independent of morphology. Our data make clear that common growth and morphology traits in yeast result from a suite of variants in master regulators, which function as a mutation-driven switch between phenotypic states.

The murine model for Hantaan virus-induced lethal disease shows two distinct paths in viral evolutionary trajectory with and without ribavirin treatment

In vitro, ribavirin acts as a lethal mutagen in Hantaan virus (HTNV)-infected Vero E6 cells, resulting in an increased mutation load and viral population extinction. In this study, we asked whether ribavirin treatment in the lethal, suckling mouse model of HTNV infection would act similarly. The HTNV genomic RNA (vRNA) copy number and infectious virus were measured in lungs of untreated and ribavirin-treated mice. In untreated, HTNV-infected mice, the vRNA copy number increased for 10 days postinfection (dpi) and thereafter remained constant through 26 dpi. Surprisingly, in ribavirin-treated, HTNV-infected mice, vRNA levels were similar to those in untreated mice between 10 and 26 dpi. Infectious virus levels, however, were different: in ribavirin-treated mice, the amount of infectious HTNV was significantly decreased relative to that in untreated mice, suggesting that ribavirin reduced the specific infectivity of the virus (amount of infectious virus produced per vRNA copy). Mutational analysis revealed a ribavirin-associated elevation in mutation frequency in HTNV vRNA similar to that previously reported in vitro. Codon-based analyses of rates of nonsynonymous (dN) and synonymous (dS) substitutions in the S segment revealed a positive selection for codons within the HTNV N protein gene in the ribavirin-treated vRNA population. In contrast, the vRNA population in untreated, HTNV-infected mice showed a lower level of diversity, reflecting purifying selection for the wild-type genome. In summary, these experiments show two different evolutionary paths that Hantavirus may take during infection in a lethal murine model of disease, as well as the importance of the in vivo host environment in the evolution of the virus, which was not apparent in our prior in vitro model system.

Emerging novel porcine parvoviruses in Europe: origin, evolution, phylodynamics and phylogeography

To elucidate the spatiotemporal phylodynamics, dispersion and evolutionary processes underlying the emergence of novel porcine parvovirus 2 (PPV2), PPV3 and PPV4 species, we analysed all available complete capsid genes, together with ours, obtained in Europe. Bayesian phylogeography indicates that Romania (PPV2 and PPV4) and Croatia (PPV3) are the most likely ancestral areas from which PPVs have subsequently spread to other European countries and regions. The timescale of our reconstruction supported a relatively recent history of the currently circulating novel PPV species (1920s to 1980s) in the domestic or sylvatic host. While PPV2 strains exhibited a large genetic exchange characterized by significant recombination and gene flow between distinct regions and hosts, PPV3 and PPV4 showed a diversification reflected by the accumulation of geographically structured polymorphisms. The RNA-like evolutionary rates detected inter- and intrahost recombination and the positive selection sites provided evidence that the PPV2-4 capsid gene plays a prominent role in host adaptation.

Bimodal evolution of the killer cell Ig-like receptor (KIR) family in New World primates

The immunoglobulin-like receptor (KIR) gene family in New World primates (Platyrrhini) has been characterized only in the owl monkey (Aotus sp.). To gain a better understanding of the KIR system in Platyrrhini, we analyzed a KIR haplotype in Ateles geoffroyi, and sequenced KIR complementary DNAs (cDNAs) from other three Atelidae species, Ateles hybridus, Ateles belzebuth, and Lagothrix lagotricha. Atelidae expressed a variable set of activating and inhibitory KIRs that diversified independently from their Catarrhini counterparts. They had a unique mechanism to generate activating receptors from inhibitory ones, involving a single nucleotide deletion in exon 7 and a change in the donor splice site of intron 7. The A. geoffroyi haplotype contained at least six gene models including a pseudogene, two coding inhibitory receptors, and three coding activating receptors. The centromeric region was in a tail-to-tail orientation with respect to the telomeric region. The owl monkey KIR haplotype shared this organization, and in phylogenetic trees, the centromeric genes clustered together with those of A. geoffroyi, whereas their telomeric genes clustered independently. KIR cDNAs from the other Atelidae species conformed to this pattern. Signatures of positive selection were found in residues predicted to interact with the major histocompatibility complex. Such signatures, however, primarily explained variability between paralogous genes but not between alleles in a locus. Atelidae, therefore, has expanded the KIR family in a bimodal fashion, where an inverted centromeric region has remained relatively conserved and the telomeric region has diversified by a rapid process of gene duplication and divergence, likely favored by positive selection for ligand binding.

Evolution of three parent genes and their retrogene copies in Drosophila species

Retrogenes form a class of gene duplicate lacking the regulatory sequences found outside of the mRNA-coding regions of the parent gene. It is not clear how a retrogene's lack of parental regulatory sequences affects the evolution of the gene pair. To explore the evolution of parent genes and retrogenes, we investigated three such gene pairs in the family Drosophilidae; in Drosophila melanogaster, these gene pairs are CG8331 and CG4960, CG17734 and CG11825, and Sep2 and Sep5. We investigated the embryonic expression patterns of these gene pairs across multiple Drosophila species. Expression patterns of the parent genes and their single copy orthologs are relatively conserved across species, whether or not a species has a retrogene copy, although there is some variation in CG8331 and CG17734. In contrast, expression patterns of the retrogene orthologs have diversified. We used the genome sequences of 20 Drosophila species to investigate coding sequence evolution. The coding sequences of the three gene pairs appear to be evolving predominantly under negative selection; however, the parent genes and retrogenes show some distinct differences in amino acid sequence. Therefore, in general, retrogene expression patterns and coding sequences are distinct compared to their parents and, in some cases, retrogene expression patterns diversify.

Suppressors of RNAi from plant viruses are subject to episodic positive selection

Viral suppressors of RNAi (VSRs) are proteins that actively inhibit the antiviral RNA interference (RNAi) immune response, providing an immune evasion route for viruses. It has been hypothesized that VSRs are engaged in a molecular ‘arms race’ with RNAi pathway genes. Two lines of evidence support this. First, VSRs from plant viruses display high sequence diversity, and are frequently gained and lost over evolutionary time scales. Second, Drosophila antiviral RNAi genes show high rates of adaptive evolution. Here, we investigate whether VSRs diversify faster than other genes and, if so, whether this is a result of positive selection, as might be expected in an arms race. By analysis of 12 plant RNA viruses, we show that the relative rate of protein evolution is higher for VSRs than for other genes, but that this is not attributable to pervasive positive selection. We argue that, because evolutionary time scales are extremely different for viruses and eukaryotes, it is improbable that viral adaptation (as measured by the ratio of non-synonymous to synonymous change) will be dominated by one-to-one coevolution with eukaryotes. Instead, for plant virus VSRs, we find strong evidence of episodic selection—diversifying selection that acts on a subset of lineages—which might be attributable to frequent shifts between different host genotypes or species.

Complete mitochondrial DNA sequences of the threadfin cichlid (Petrochromis trewavasae) and the blunthead cichlid (Tropheus moorii) and patterns of mitochondrial genome evolution in cichlid fishes

The cichlid fishes of the East African Great Lakes represent a model especially suited to study adaptive radiation and speciation. With several African cichlid genome projects being in progress, a promising set of closely related genomes is emerging, which is expected to serve as a valuable data base to solve questions on genotype-phenotype relations. The mitochondrial (mt) genomes presented here are the first results of the assembly and annotation process for two closely related but eco-morphologically highly distinct Lake Tanganyika cichlids, Petrochromis trewavasae and Tropheus moorii. The genomic sequences comprise 16,588 bp (P. trewavasae) and 16,590 bp (T. moorii), and exhibit the typical mitochondrial structure, with 13 protein-coding genes, 2 rRNA genes, 22 tRNA genes, and a non-coding control region. Analyses confirmed that the two species are very closely related with an overall sequence similarity of 96%. We analyzed the newly generated sequences in the phylogenetic context of 21 published labroid fish mitochondrial genomes. Consistent with other vertebrates, the D-loop region was found to evolve faster than protein-coding genes, which in turn are followed by the rRNAs; the tRNAs vary greatly in the rate of sequence evolution, but on average evolve the slowest. Within the group of coding genes, ND6 evolves most rapidly. Codon usage is similar among examined cichlid tribes and labroid families; although a slight shift in usage patterns down the gene tree could be observed. Despite having a clearly different nucleotide composition, ND6 showed a similar codon usage. C-terminal ends of Cox1 exhibit variations, where the varying number of amino acids is related to the structure of the obtained phylogenetic tree. This variation may be of functional relevance for Cox1 synthesis.

Evolution and divergence of the mammalian SAMD9/SAMD9L gene family

BACKGROUND

The physiological functions of the human Sterile Alpha Motif Domain-containing 9 (SAMD9) gene and its chromosomally adjacent paralogue, SAMD9-like (SAMD9L), currently remain unknown. However, the direct links between the deleterious mutations or deletions in these two genes and several human disorders, such as inherited inflammatory calcified tumors and acute myeloid leukemia, suggest their biological importance. SAMD9 and SAMD9L have also recently been shown to play key roles in the innate immune responses to stimuli such as viral infection. We were particularly interested in understanding the mammalian evolutionary history of these two genes. The phylogeny of SAMD9 and SAMD9L genes was reconstructed using the Maximum Likelihood method. Furthermore, six different methods were applied to detect SAMD9 and SAMD9L codons under selective pressure: the site-specific model M8 implemented in the codeml program in PAML software and five methods available on the Datamonkey web server, including the Single Likelihood Ancestor Counting method, the Fixed Effect Likelihood method, the Random Effect Likelihood method, the Mixed Effects Model of Evolution method and the Fast Unbiased Bayesian AppRoximation method. Additionally, the house mouse (Mus musculus) genome has lost the SAMD9 gene, while keeping SAMD9L intact, prompting us to investigate whether this loss is a unique event during evolution.

RESULTS

Our evolutionary analyses suggest that SAMD9 and SAMD9L arose through an ancestral gene duplication event after the divergence of Marsupialia from Placentalia. Additionally, selection analyses demonstrated that both genes have been subjected to positive evolutionary selection. The absence of either SAMD9 or SAMD9L genes from some mammalian species supports a partial functional redundancy between the two genes.

CONCLUSIONS

To the best of our knowledge, this work is the first study on the evolutionary history of mammalian SAMD9 and SAMD9L genes. We conclude that evolutionary selective pressure has acted on both of these two genes since their divergence, suggesting their importance in multiple cellular processes, such as the immune responses to viral pathogens.

Evolution by selection, recombination, and gene duplication in MHC class I genes of two Rhacophoridae species

BACKGROUND

Comparison of major histocompatibility complex (MHC) genes across vertebrate species can reveal molecular mechanisms underlying the evolution of adaptive immunity-related proteins. As the first terrestrial tetrapods, amphibians deserve special attention because of their exposure to probably increased spectrum of microorganisms compared with ancestral aquatic fishes. Knowledge regarding the evolutionary patterns and mechanisms associated with amphibian MHC genes remains limited. The goal of the present study was to isolate MHC class I genes from two Rhacophoridae species (Rhacophorus omeimontis and Polypedates megacephalus) and examine their evolution.

RESULTS

We identified 27 MHC class I alleles spanning the region from exon 2 to 4 in 38 tree frogs. The available evidence suggests that these 27 sequences all belong to classical MHC class I (MHC Ia) genes. Although several anuran species only display one MHC class Ia locus, at least two or three loci were observed in P. megacephalus and R. omeimontis, indicating that the number of MHC class Ia loci varies among anuran species. Recombination events, which mainly involve the entire exons, played an important role in shaping the genetic diversity of the 27 MHC class Ia alleles. In addition, signals of positive selection were found in Rhacophoridae MHC class Ia genes. Amino acid sites strongly suggested by program to be under positive selection basically accorded with the putative antigen binding sites deduced from crystal structure of human HLA. Phylogenetic relationships among MHC class I alleles revealed the presence of trans-species polymorphisms.

CONCLUSIONS

In the two Rhacophoridae species (1) there are two or three MHC class Ia loci; (2) recombination mainly occurs between the entire exons of MHC class Ia genes; (3) balancing selection, gene duplication and recombination all contribute to the diversity of MHC class Ia genes. These findings broaden our knowledge on the evolution of amphibian MHC systems.

Dracula's children: molecular evolution of vampire bat venom

While vampire bat oral secretions have been the subject of intense research, efforts have concentrated only on two components: DSPA (Desmodus rotundus salivary plasminogen activator) and Draculin. The molecular evolutionary history of DSPA has been elucidated, while conversely draculin has long been known from only a very small fragment and thus even the basic protein class was not even established. Despite the fact that vampire bat venom has a multitude of effects unaccounted by the documented bioactivities of DSPA and draculin, efforts have not been made to establish what other bioactive proteins are secreted by their submaxillary gland. In addition, it has remained unclear whether the anatomically distinct anterior and posterior lobes of the submaxillary gland are evolving on separate gene expression trajectories or if they remain under the shared genetic control. Using a combined proteomic and transcriptomic approach, we show that identical proteins are simultaneously expressed in both lobes. In addition to recovering the known structural classes of DSPA, we recovered a novel DSPA isoform as well as obtained a very large sequence stretch of draculin and thus established that it is a mutated version of the lactotransferrin scaffold. This study reveals a much more complex secretion profile than previously recognised. In addition to obtaining novel versions of scaffolds convergently recruited into other venoms (allergen-like, CRiSP, kallikrein, Kunitz, lysozyme), we also documented novel expression of small peptides related to calcitonin, PACAP, and statherin. Other overexpressed protein types included BPI-fold, lacritin, and secretoglobin. Further, we investigate the molecular evolution of various vampire bat venom-components and highlight the dominant role of positive selection in the evolution of these proteins. Conspicuously many of the proteins identified in the proteome were found to be homologous to proteins with known activities affecting vasodilation and platelet aggregation. We show that vampire bat venom proteins possibly evade host immune response by the mutation of the surface chemistry through focal mutagenesis under the guidance of positive Darwinian selection. These results not only contribute to the body of knowledge regarding haematophagous venoms but also provide a rich resource for novel lead compounds for use in drug design and development.

BIOLOGICAL SIGNIFICANCE

These results have direct implications in understanding the molecular evolutionary history of vampire bat venom. The unusual peptides discovered reinforce the value of studying such neglected taxon for biodiscovery.

Genotype patterns of contemporary reassorted H3N2 virus in US swine

To understand the evolution of swine-origin H3N2v influenza viruses that have infected 320 humans in the USA since August 2011, we performed a phylogenetic analysis at a whole genome scale of North American swine influenza viruses (n  =  200). All viral isolates evolved from the prototypical North American H3 cluster 4 (c4), with evidence for further diversification into subclusters. At least ten distinct reassorted H3N2/pandemic H1N1 (rH3N2p) genotypes were identified in swine. Genotype 1 (G1) was most frequently detected in swine and all human H3N2v viruses clustered within a single G1 clade. These data suggest that the genetic requirements for transmission to humans may be restricted to a specific subset of swine viruses. Mutations at putative antigenic sites as well as reduced serological cross-reactivity among the H3 subclusters suggest antigenic drift of these contemporary viruses.

Evolutionary dynamics of West Nile virus in the United States, 1999-2011: phylogeny, selection pressure and evolutionary time-scale analysis

West Nile virus (WNV), an arbovirus maintained in a bird-mosquito enzootic cycle, can infect other vertebrates including humans. WNV was first reported in the US in 1999 where, to date, three genotypes belonging to WNV lineage I have been described (NY99, WN02, SW/WN03). We report here the WNV sequences obtained from two birds, one mosquito, and 29 selected human samples acquired during the US epidemics from 2006–2011 and our examination of the evolutionary dynamics in the open-reading frame of WNV isolates reported from 1999–2011. Maximum-likelihood and Bayesian methods were used to perform the phylogenetic analyses and selection pressure analyses were conducted with the HyPhy package. Phylogenetic analysis identified human WNV isolates within the main WNV genotypes that have circulated in the US. Within genotype SW/WN03, we have identified a cluster with strains derived from blood donors and birds from Idaho and North Dakota collected during 2006–2007, termed here MW/WN06. Using different codon-based and branch-site selection models, we detected a number of codons subjected to positive pressure in WNV genes. The mean nucleotide substitution rate for WNV isolates obtained from humans was calculated to be 5.06×10⁻⁴ substitutions/site/year (s/s/y). The Bayesian skyline plot shows that after a period of high genetic variability following the introduction of WNV into the US, the WNV population appears to have reached genetic stability. The establishment of WNV in the US represents a unique opportunity to understand how an arbovirus adapts and evolves in a naïve environment. We describe a novel, well-supported cluster of WNV formed by strains collected from humans and birds from Idaho and North Dakota. Adequate genetic surveillance is essential to public health since new mutants could potentially affect viral pathogenesis, decrease performance of diagnostic assays, and negatively impact the efficacy of vaccines and the development of specific therapies.

West Nile Virus (WNV) is a mosquito-borne virus of African origin that is widespread around the world. The WNV life-cycle involves mosquitoes and birds, but humans and other animals can be infected, although they are not considered to be important players in the transmission cycle. Clinically, most WNV infections are unapparent, but the virus can disseminate to the central nervous system causing a potentially fatal neurological disease, especially in susceptible populations including elderly and immunocompromised individuals. West Nile virus can also be transmitted by organ transplant and by transfusion of blood and blood components. Like other arboviruses, WNV has the extraordinary capacity of growing in the different microenvironments represented by the invertebrate vector and the vertebrate hosts. From an evolutionary standpoint, the arrival of WNV in the US in 1999 represents a unique opportunity to explore the processes involved in the adaptation and dissemination of an arbovirus in a naïve environment. From the study of WNV sequences, we can not only learn about the evolutionary mechanisms that govern arboviruses, but also update diagnostic tests that rely on the detection of the viral genome upon the occurrence of mutations and study the existence of genetic markers that may be responsible for increases in clinical cases and their severity.

Genetic characterization of Betacoronavirus lineage C viruses in bats reveals marked sequence divergence in the spike protein of pipistrellus bat coronavirus HKU5 in Japanese pipistrelle: implications for the origin of the novel Middle East respiratory syndrome coronavirus

While the novel Middle East respiratory syndrome coronavirus (MERS-CoV) is closely related to Tylonycteris bat CoV HKU4 (Ty-BatCoV HKU4) and Pipistrellus bat CoV HKU5 (Pi-BatCoV HKU5) in bats from Hong Kong, and other potential lineage C betacoronaviruses in bats from Africa, Europe, and America, its animal origin remains obscure. To better understand the role of bats in its origin, we examined the molecular epidemiology and evolution of lineage C betacoronaviruses among bats. Ty-BatCoV HKU4 and Pi-BatCoV HKU5 were detected in 29% and 25% of alimentary samples from lesser bamboo bat (Tylonycteris pachypus) and Japanese pipistrelle (Pipistrellus abramus), respectively. Sequencing of their RNA polymerase (RdRp), spike (S), and nucleocapsid (N) genes revealed that MERS-CoV is more closely related to Pi-BatCoV HKU5 in RdRp (92.1% to 92.3% amino acid [aa] identity) but is more closely related to Ty-BatCoV HKU4 in S (66.8% to 67.4% aa identity) and N (71.9% to 72.3% aa identity). Although both viruses were under purifying selection, the S of Pi-BatCoV HKU5 displayed marked sequence polymorphisms and more positively selected sites than that of Ty-BatCoV HKU4, suggesting that Pi-BatCoV HKU5 may generate variants to occupy new ecological niches along with its host in diverse habitats. Molecular clock analysis showed that they diverged from a common ancestor with MERS-CoV at least several centuries ago. Although MERS-CoV may have diverged from potential lineage C betacoronaviruses in European bats more recently, these bat viruses were unlikely to be the direct ancestor of MERS-CoV. Intensive surveillance for lineage C betaCoVs in Pipistrellus and related bats with diverse habitats and other animals in the Middle East may fill the evolutionary gap.

A 175 million year history of T cell regulatory molecules reveals widespread selection, with adaptive evolution of disease alleles

T cell activation plays a central role in immune response and in the maintenance of self-tolerance. We analyzed the evolutionary history of T cell regulatory molecules. Nine genes involved in triggering T cell activation or in regulating the ensuing response evolved adaptively in mammals. Several positively selected sites overlap with positions interacting with the binding partner or with cellular components. Population genetic analysis in humans revealed a complex scenario of local (FASLG, CD40LG, HAVCR2) and worldwide (FAS, ICOSLG) adaptation and H. sapiens-to-Neandertal gene flow (gene transfer between populations). Disease variants in these genes are preferential targets of pathogen-driven selection, and a Crohn's disease risk polymorphism targeted by bacterial-driven selection modulates the expression of ICOSLG in response to a bacterial superantigen. Therefore, we used evolutionary information to generate experimentally testable hypotheses concerning the function of specific genetic variants and indicate that adaptation to infection underlies the maintenance of autoimmune risk alleles.

The comparative genomics of human respiratory syncytial virus subgroups A and B: genetic variability and molecular evolutionary dynamics

Genomic variation and related evolutionary dynamics of human respiratory syncytial virus (RSV), a common causative agent of severe lower respiratory tract infections, may affect its transmission behavior. RSV evolutionary patterns are likely to be influenced by a precarious interplay between selection favoring variants with higher replicative fitness and variants that evade host immune responses. Studying RSV genetic variation can reveal both the genes and the individual codons within these genes that are most crucial for RSV survival. In this study, we conducted genetic diversity and evolutionary rate analyses on 36 RSV subgroup B (RSV-B) whole-genome sequences. The attachment protein, G, was the most variable protein; accordingly, the G gene had a higher substitution rate than other RSV-B genes. Overall, less genetic variability was found among the available RSV-B genome sequences than among RSV-A genome sequences in a comparable sample. The mean substitution rates of the two subgroups were, however, similar (for subgroup A, 6.47 × 10(-4) substitutions/site/year [95% credible interval {CI 95%}, 5.56 × 10(-4) to 7.38 × 10(-4)]; for subgroup B, 7.76 × 10(-4) substitutions/site/year [CI 95%, 6.89 × 10(-4) to 8.58 × 10(-4)]), with the time to their most recent common ancestors (TMRCAs) being much lower for RSV-B (19 years) than for RSV-A (46.8 years). The more recent RSV-B TMRCA is apparently the result of a genetic bottleneck that, over longer time scales, is still compatible with neutral population dynamics. Whereas the immunogenic G protein seems to require high substitution rates to ensure immune evasion, strong purifying selection in conserved proteins such as the fusion protein and nucleocapsid protein is likely essential to preserve RSV viability.

Citrus tristeza virus: Evolution of Complex and Varied Genotypic Groups

Amongst the Closteroviridae, Citrus tristeza virus (CTV) is almost unique in possessing a number of distinct and characterized strains, isolates of which produce a wide range of phenotype combinations among its different hosts. There is little understanding to connect genotypes to phenotypes, and to complicate matters more, these genotypes are found throughout the world as members of mixed populations within a single host plant. There is essentially no understanding of how combinations of genotypes affect symptom expression and disease severity. We know little about the evolution of the genotypes that have been characterized to date, little about the biological role of their diversity and particularly, about the effects of recombination. Additionally, genotype grouping has not been standardized. In this study we utilized an extensive array of CTV genomic information to classify the major genotypes, and to determine the major evolutionary processes that led to their formation and subsequent retention. Our analyses suggest that three major processes act on these genotypes: (1) ancestral diversification of the major CTV lineages, followed by (2) conservation and co-evolution of the major functional domains within, though not between CTV genotypes, and (3) extensive recombination between lineages that have given rise to new genotypes that have subsequently been retained within the global population. The effects of genotype diversity and host-interaction are discussed, as is a proposal for standardizing the classification of existing and novel CTV genotypes.

Evidence of recombination and genetic diversity in southern rice black-streaked dwarf virus

Southern rice black-streaked dwarf virus (SRBSDV) causes one of the most serious viral diseases of rice in China and Vietnam. Sequence identities of S10, encoding the major capsid protein, were 98.0 %-100 % and 98.3 %-100 % at the nucleotide and amino acid level, respectively. Our results suggest that the codon at position 550 of S10 is under positive selection, while most of the other codons are under neutral evolution. Putative recombination events were identified in genomic RNA segments S1, 2, 4, 5, 6 and 10, which are rare in plant-infecting dsRNA viruses. This study reveals the current state of SRBSDV evolution.

A phylomedicine approach to understanding the evolution of auditory sensory perception and disease in mammals

Hereditary deafness affects 0.1% of individuals globally and is considered as one of the most debilitating diseases of man. Despite recent advances, the molecular basis of normal auditory function is not fully understood and little is known about the contribution of single-nucleotide variations to the disease. Using cross-species comparisons of 11 ‘deafness’ genes (Myo15, Ush1 g, Strc, Tecta, Tectb, Otog, Col11a2, Gjb2, Cldn14, Kcnq4, Pou3f4) across 69 evolutionary and ecologically divergent mammals, we elucidated whether there was evidence for: (i) adaptive evolution acting on these genes across mammals with similar hearing capabilities; and, (ii) regions of long-term evolutionary conservation within which we predict disease-associated mutations should occur. We find evidence of adaptive evolution acting on the eutherian mammals in Myo15, Otog and Tecta. Examination of selection pressures in Tecta and Pou3f4 across a taxonomic sample that included a wide representation of auditory specialists, the bats, did not uncover any evidence for a role in echolocation. We generated ‘conservation indices’ based on selection estimates at nucleotide sites and found that known disease mutations fall within sites of high evolutionary conservation. We suggest that methods such as this, derived from estimates of evolutionary conservation using phylogenetically divergent taxa, will help to differentiate between deleterious and benign mutations.

tat Exon 1 exhibits functional diversity during HIV-1 subtype C primary infection

Human immunodeficiency virus type 1 (HIV-1) Tat is a mediator of viral transcription and is involved in the control of virus replication. However, associations between HIV-1 Tat diversity and functional effects during primary HIV-1 infection are still unclear. We estimated selection pressures in tat exon 1 using the mixed-effects model of evolution with 672 viral sequences generated from 20 patients infected with HIV-1 subtype C (HIV-1C) over 500 days postseroconversion. tat exon 1 residues 3, 4, 21, 24, 29, 39, and 68 were under positive selection, and we established that specific amino acid signature patterns were apparent in primary HIV-1C infection compared with chronic infection. We assessed the impact of these mutations on long terminal repeat (LTR) activity and found that Tat activity was negatively affected by the Ala(21) substitution identified in 13/20 (65%) of patients, which reduced LTR activity by 88% (± 1%) (P < 0.001). The greatest increase in Tat activity was seen with the Gln(35)/Lys(39) double mutant that resulted in an additional 49% (± 14%) production of LTR-driven luciferase (P = 0.012). There was a moderate positive correlation between Tat-mediated LTR activity and HIV-1 RNA in plasma (P = 0.026; r = 0.400) after 180 days postseroconversion that was reduced by 500 days postseroconversion (P = 0.043; r = 0.266). Although Tat activation of the LTR is not a strong predictor of these clinical variables, there are significant linear relationships between Tat transactivation and patients' plasma viral loads and CD4 counts, highlighting the complex interplay between Tat mutations in early HIV-1C infection.

HIV-1 conserved-element vaccines: relationship between sequence conservation and replicative capacity

To overcome the problem of HIV-1 variability, candidate vaccine antigens have been designed to be composed of conserved elements of the HIV-1 proteome. Such candidate vaccines could be improved with a better understanding of both HIV-1 evolutionary constraints and the fitness cost of specific mutations. We evaluated the in vitro fitness cost of 23 mutations engineered in the HIV-1 subtype B Gag-p24 Center-of-Tree (COT) protein through fitness competition assays. While some mutations at conserved sites exacted a high fitness cost, as expected under the assumption that the most conserved residue confers the highest fitness, there was no overall strong relationship between sequence conservation and replicative capacity. By comparing sites that have evolved since the beginning of the epidemic to those that have remain unchanged, we found that sites that have evolved over time were more likely to correspond to HLA-associated sites and that their mutation had limited fitness costs. Our data showed no transcendent link between high conservation and high fitness cost, indicating that merely focusing on conserved segments of HIV-1 would not be sufficient for a successful vaccine strategy. Nonetheless, a subset of sites exacted a high fitness cost upon mutation--these sites have been under selective pressure to change since the beginning of the epidemic but have proved virtually nonmutable and could constitute preferred targets for vaccine design.

Inference of functional divergence among proteins when the evolutionary process is non-stationary

Functional shifts during protein evolution are expected to yield shifts in substitution rate, and statistical methods can test for this at both codon and amino acid levels. Although methods based on models of sequence evolution serve as powerful tools for studying evolutionary processes, violating underlying assumptions can lead to false biological conclusions. It is not unusual for functional shifts to be accompanied by changes in other aspects of the evolutionary process, such as codon or amino acid frequencies. However, models used to test for functional divergence assume these frequencies remain constant over time. We employed simulation to investigate the impact of non-stationary evolution on functional divergence inference. We investigated three likelihood ratio tests based on codon models and found varying degrees of sensitivity. Joint effects of shifts in frequencies and selection pressures can be large, leading to false signals for positive selection. Amino acid-based tests (FunDi and Bivar) were also compromised when several aspects of the substitution process were not adequately modeled. We applied the same tests to a core genome "scan" for functional divergence between light-adapted ecotypes of the cyanobacteria Prochlorococcus, and carried out gene-specific simulations for ten genes. Results of those simulations illustrated how the inference of functional divergence at the genomic level can be seriously impacted by model misspecification. Although computationally costly, simulations motivated by data in hand are warranted when several aspects of the substitution process are either misspecified or not included in the models upon which the statistical tests were built.

Structure and evolution of the sea star egg receptor for sperm bindin

Selection on coevolving sperm- and egg-recognition molecules is a potent engine of population divergence leading to reproductive isolation and speciation. The study of receptor-ligand pairs can reveal co-evolution of male- and female-expressed genes or differences between their evolution in response to selective factors such as sperm competition and sexual conflict. Phylogeographical studies of these patterns have been limited by targeted gene methods that favour short protein-coding sequences amplifiable by PCR. Here, I use high-throughput transcriptomic methods to characterize the structure and divergence of full-length coding sequences for the gene encoding the protein component of a large complex egg surface glycopeptide receptor for the sperm acrosomal protein bindin from the sea star Patiria miniata. I used a simple but effective method for resolving nucleotide polymorphisms into haplotypes for phylogeny-based analyses of selection. The protein domain organization of sea star egg bindin receptor (EBR1) was similar to sea urchins and included a pair of protein-recognition domains plus a series of tandem repeat domains of two types. Two populations separated by a well-characterized phylogeographical break included lineages of EBR1 alleles under positive selection at several codons (similar to selection on sperm bindin in the same populations). However, these populations shared the same alleles that were under selection for amino acid differences at multiple codons (unlike the pattern of selection for population divergence in sperm bindin). The significance of positively selected EBR1 domains and alleles could be tested in functional analyses of fertilization rates associated with EBR1 (and bindin) polymorphisms.

FUBAR: a fast, unconstrained bayesian approximation for inferring selection

Model-based analyses of natural selection often categorize sites into a relatively small number of site classes. Forcing each site to belong to one of these classes places unrealistic constraints on the distribution of selection parameters, which can result in misleading inference due to model misspecification. We present an approximate hierarchical Bayesian method using a Markov chain Monte Carlo (MCMC) routine that ensures robustness against model misspecification by averaging over a large number of predefined site classes. This leaves the distribution of selection parameters essentially unconstrained, and also allows sites experiencing positive and purifying selection to be identified orders of magnitude faster than by existing methods. We demonstrate that popular random effects likelihood methods can produce misleading results when sites assigned to the same site class experience different levels of positive or purifying selection--an unavoidable scenario when using a small number of site classes. Our Fast Unconstrained Bayesian AppRoximation (FUBAR) is unaffected by this problem, while achieving higher power than existing unconstrained (fixed effects likelihood) methods. The speed advantage of FUBAR allows us to analyze larger data sets than other methods: We illustrate this on a large influenza hemagglutinin data set (3,142 sequences). FUBAR is available as a batch file within the latest HyPhy distribution (http://www.hyphy.org), as well as on the Datamonkey web server (http://www.datamonkey.org/).

Detection of positive selection eliminating effects of structural constraints in hemagglutinin of H3N2 human influenza A virus

In the evolutionary studies of proteins, the average effect of natural selection operating on amino acid mutations may be examined by comparing the numbers of synonymous (dS) and nonsynonymous (dN) substitutions that have accumulated during the same time period. In this method, destabilizing mutations occurring across protein molecules may interfere with detection of natural selection, particularly positive selection, operating on other mutations. Here an attempt to detect positive selection eliminating effects of structural constraints is demonstrated using hemagglutinin (HA) of H3N2 human influenza A virus as an example. Compatible and incompatible amino acids were inferred at each site from the computational analysis of three-dimensional structure using the thermodynamic stability as an indicator, and natural selection was examined by comparing dS and dN among compatible amino acids. In the analysis of 2701 nucleotide sequences for the entire coding region of HA, the new method identified twice as many positively selected amino acid sites as the ordinary method (16 and 4 sites in the former method without and with correction for multiple testing, respectively, and 8 and 2 sites in the latter method). Positively selected sites were involved in epitopes, receptor-binding pocket, epistasis, and stabilization, which appeared to be biologically reasonable. Nevertheless, there still appeared to be several problems, which may largely render this method conservative. It may be effective to analyze many densely sampled sequences in this method.

Origin and evolution of dishevelled

Dishevelled (Dsh or Dvl) is an important signaling protein, playing a key role in Wnt signaling and relaying cellular information for several developmental pathways. Dsh is highly conserved among metazoans and has expanded into a multigene family in most bilaterian lineages, including vertebrates, planarians, and nematodes. These orthologs, where explored, are known to have considerable overlap in function, but evidence for functional specialization continues to mount. We performed a comparative analysis of Dsh across animals to explore protein architecture and identify conserved and divergent features that could provide insight into functional specialization with an emphasis on invertebrates, especially nematodes. We find evidence of dynamic evolution of Dsh, particularly among nematodes, with taxa varying in ortholog number from one to three. We identify a new domain specific to some nematode lineages and find an unexpected nuclear localization signal conserved in many Dsh orthologs. Our findings raise questions of protein evolution in general and provide clues as to how animals have dealt with the complex intricacies of having a protein, such as Dsh, act as a central messenger hub connected to many different and vitally important pathways. We discuss our findings in the context of functional specialization and bring many testable hypotheses to light.

Convergence and divergence in the evolution of the APOBEC3G-Vif interaction reveal ancient origins of simian immunodeficiency viruses

Naturally circulating lentiviruses are abundant in African primate species today, yet their origins and history of transmitting between hosts remain obscure. As a means to better understand the age of primate lentiviruses, we analyzed primate genomes for signatures of lentivirus-driven evolution. Specifically, we studied the adaptive evolution of host restriction factor APOBEC3G (A3G) in Old World Monkey (OWM) species. We find recurrent mutation of A3G in multiple primate lineages at sites that determine susceptibility to antagonism by the lentiviral accessory protein Vif. Using a broad panel of SIV Vif isolates, we demonstrate that natural variation in OWM A3G confers resistance to Vif-mediated degradation, suggesting that adaptive variants of the host factor were selected upon exposure to pathogenic lentiviruses at least 5–6 million years ago (MYA). Furthermore, in members of the divergent Colobinae subfamily of OWM, a multi-residue insertion event in A3G that arose at least 12 MYA blocks the activity of Vif, suggesting an even more ancient origin of SIV. Moreover, analysis of the lentiviruses associated with Colobinae monkeys reveal that the interface of the A3G-Vif interaction has shifted and given rise to a second genetic conflict. Our analysis of virus-driven evolution describes an ancient yet ongoing genetic conflict between simian primates and lentiviruses on a million-year time scale.

The emergence of AIDS in the late 20^th century has provoked studies to better understand the evolutionary history of viruses and the factors that govern their spread. Pandemic human immunodeficiency virus-type 1 (HIV-1), which currently infects 34 million people worldwide, emerged following the transmission of a lentivirus between chimpanzees and humans. A growing list of apparently nonpathogenic, species-specific strains has now been characterized in dozens of African primates, suggesting that primate lentiviruses are older and more widespread than originally thought. To estimate the extent to which primates and lentiviruses have coexisted, we examined the interaction between host and virus on a molecular level and tracked its dynamics over evolutionary time. We report that the immunity factor APOBEC3G is evolving in tandem with the lentiviral accessory gene vif, allowing us to associate instances of host evolution with instances of lentivirus infection in deep and shallow timescales. Specifically, we show that the region of APOBEC3G targeted by Vif is adaptively diversifying in independent primate lineages in a manner that suggests that lentiviruses are millions of years old. Our study reveals that, while primate lentiviruses may have modern consequences for human health, they have ancient origins in our non-human primate relatives.

Residue mutations and their impact on protein structure and function: detecting beneficial and pathogenic changes

The present review focuses on the evolution of proteins and the impact of amino acid mutations on function from a structural perspective. Proteins evolve under the law of natural selection and undergo alternating periods of conservative evolution and of relatively rapid change. The likelihood of mutations being fixed in the genome depends on various factors, such as the fitness of the phenotype or the position of the residues in the three-dimensional structure. For example, co-evolution of residues located close together in three-dimensional space can occur to preserve global stability. Whereas point mutations can fine-tune the protein function, residue insertions and deletions (‘decorations’ at the structural level) can sometimes modify functional sites and protein interactions more dramatically. We discuss recent developments and tools to identify such episodic mutations, and examine their applications in medical research. Such tools have been tested on simulated data and applied to real data such as viruses or animal sequences. Traditionally, there has been little if any cross-talk between the fields of protein biophysics, protein structure–function and molecular evolution. However, the last several years have seen some exciting developments in combining these approaches to obtain an in-depth understanding of how proteins evolve. For example, a better understanding of how structural constraints affect protein evolution will greatly help us to optimize our models of sequence evolution. The present review explores this new synthesis of perspectives.

Genetic variability among complete human respiratory syncytial virus subgroup A genomes: bridging molecular evolutionary dynamics and epidemiology

Human respiratory syncytial virus (RSV) is an important cause of severe lower respiratory tract infections in infants and the elderly. In the vast majority of cases, however, RSV infections run mild and symptoms resemble those of a common cold. The immunological, clinical, and epidemiological profile of severe RSV infections suggests a disease caused by a virus with typical seasonal transmission behavior, lacking clear-cut virulence factors, but instead causing disease by modifying the host’s immune response in a way that stimulates pathogenesis. Yet, the interplay between RSV-evoked immune responses and epidemic behavior, and how this affects the genomic evolutionary dynamics of the virus, remains poorly understood. Here, we present a comprehensive collection of 33 novel RSV subgroup A genomes from strains sampled over the last decade, and provide the first measurement of RSV-A genomic diversity through time in a phylodynamic framework. In addition, we map amino acid substitutions per protein to determine mutational hotspots in specific domains. Using Bayesian genealogical inference, we estimated the genomic evolutionary rate to be 6.47×10⁻⁴ (credible interval: 5.56×10⁻⁴, 7.38×10⁻⁴) substitutions/site/year, considerably slower than previous estimates based on G gene sequences only. The G gene is however marked by elevated substitution rates compared to other RSV genes, which can be attributed to relaxed selective constraints. In line with this, site-specific selection analyses identify the G gene as the major target of diversifying selection. Importantly, statistical analysis demonstrates that the immune driven positive selection does not leave a measurable imprint on the genome phylogeny, implying that RSV lineage replacement mainly follows nonselective epidemiological processes. The roughly 50 years of RSV-A genomic evolution are characterized by a constant population size through time and general co-circulation of lineages over many epidemic seasons – a conclusion that might be taken into account when developing future therapeutic and preventive strategies.

Identification of positive selection in disease response genes within members of the Poaceae

Millions of years of coevolution between plants and pathogens can leave footprints on their genomes and genes involved on this interaction are expected to show patterns of positive selection in which novel, beneficial alleles are rapidly fixed within the population. Using information about upregulated genes in maize during Colletotrichum graminicola infection and resources available in the Phytozome database, we looked for evidence of positive selection in the Poaceae lineage, acting on protein coding sequences related with plant defense. We found six genes with evidence of positive selection and another eight with sites showing episodic selection. Some of them have already been described as evolving under positive selection, but others are reported here for the first time including genes encoding isocitrate lyase, dehydrogenases, a multidrug transporter, a protein containing a putative leucine-rich repeat and other proteins with unknown functions. Mapping positively selected residues onto the predicted 3-D structure of proteins showed that most of them are located on the surface, where proteins are in contact with other molecules. We present here a set of Poaceae genes that are likely to be involved in plant defense mechanisms and have evidence of positive selection. These genes are excellent candidates for future functional validation.

Evolutionary dynamics of the interferon-induced transmembrane gene family in vertebrates

Vertebrate interferon-induced transmembrane (IFITM) genes have been demonstrated to have extensive and diverse functions, playing important roles in the evolution of vertebrates. Despite observance of their functionality, the evolutionary dynamics of this gene family are complex and currently unknown. Here, we performed detailed evolutionary analyses to unravel the evolutionary history of the vertebrate IFITM family. A total of 174 IFITM orthologous genes and 112 pseudogenes were identified from 27 vertebrate genome sequences. The vertebrate IFITM family can be divided into immunity-related IFITM (IR-IFITM), IFITM5 and IFITM10 sub-families in phylogeny, implying origins from three different progenitors. In general, vertebrate IFITM genes are located in two loci, one containing the IFITM10 gene, and the other locus containing IFITM5 and various numbers of IR-IFITM genes. Conservation of evolutionary synteny was observed in these IFITM genes. Significant functional divergence was detected among the three IFITM sub-families. No gene duplication or positive selection was found in IFITM5 sub-family, implying the functional conservation of IFITM5 in vertebrate evolution, which is involved in bone formation. No IFITM5 locus was identified in the marmoset genome, suggesting a potential association with the tiny size of this monkey. The IFITM10 sub-family was divided into two groups: aquatic and terrestrial types. Functional divergence was detected between the two groups, and five IFITM10-like genes from frog were dispersed into the two groups. Both gene duplication and positive selection were observed in aquatic vertebrate IFITM10-like genes, indicating that IFITM10 might be associated with the adaptation to aquatic environments. A large number of lineage- and species-specific gene duplications were observed in IR-IFITM sub-family and positive selection was detected in IR-IFITM of primates and rodents. Because primates have experienced a long history of viral infection, such rapid expansion and positive selection suggests that the evolution of primate IR-IFITM genes is associated with broad-spectrum antiviral activity.

Variable evolutionary routes to host establishment across repeated rabies virus host shifts among bats

Determining the genetic pathways that viruses traverse to establish in new host species is crucial to predict the outcome of cross-species transmission but poorly understood for most host-virus systems. Using sequences encoding 78% of the rabies virus genome, we explored the extent, repeatability and dynamic outcome of evolution associated with multiple host shifts among New World bats. Episodic bursts of positive selection were detected in several viral proteins, including regions associated with host cell interaction and viral replication. Host shifts involved unique sets of substitutions, and few sites exhibited repeated evolution across adaptation to many bat species, suggesting diverse genetic determinants over host range. Combining these results with genetic reconstructions of the demographic histories of individual viral lineages revealed that although rabies viruses shared consistent three-stage processes of emergence in each new bat species, host shifts involving greater numbers of positively selected substitutions had longer delays between cross-species transmission and enzootic viral establishment. Our results point to multiple evolutionary routes to host establishment in a zoonotic RNA virus that may influence the speed of viral emergence.