Molecular Evolution of the Genomic RNA of Apple Stem Grooving Capillovirus

The complete genome of the German isolate AC of Apple stem grooving virus (ASGV) was sequenced. It encodes two overlapping open reading frames (ORFs), similarly to previously described ASGV isolates. Two regions of high variability were detected between the ASGV isolates, variable region 1 (V1, from amino acids (aa) 532 to 570), and variable region 2 (V2, from aa 1,583 to 1,868). The phylogenetic analysis of the V1 and V2 regions suggested that the ASGV diversity was structured by host plant species rather than geographical origin. The dN/dS ratio between nonsynonymous and synonymous nucleotide substitution rates varied greatly along the ASGV genome. Most of ORF1 showed predominant negative selection except for the two regions V1 and V2. V1 showed an elevated dN and an average dS when compared to the ORF1 background but no significant positive selection was detected. The V2 region of ORF1 showed an elevated dN and a low dS when compared to the ORF1 background with an average dN/dS????3.0 indicative of positive selection. However, the V2 area includes overlapping ORFs, making the dN/dS estimate biased. Joint estimates of the selection intensity in the different ORFs by a recent method indicated that this region of ORF1 was in fact evolving close to neutrality. This was convergent with previous results showing that introduction of stop codons in this region of ORF1 did not impair plant infection. These data suggest that the elimination of a stop codon caused the overprinting of a novel coding region over the ancestral ORF.     

Evolution of duplicate genes in a tetraploid animal, Xenopus laevis

To understand the evolution of duplicate genes, we compared rates of nucleotide substitution between 17 pairs of nonallelic duplicated genes in the tetraploid frog Xenopus laevis with rates between the orthologous loci of human and rodent. For all duplicated X. laevis genes, the number of synonymous substitutions per site (dS) was greater than the number of nonsynonymous substitutions per site (dN), indicating that these genes are subject to purifying selection. There was also a significant positive correlation (r = 0.915) between dN for the X. laevis genes and dN for the mammalian genes, suggesting that, at the amino acid level, the X. laevis genes and the mammalian genes are under similar constraints. Results of relative-rate tests showed nearly equal rates of nonsynonymous substitution in each copy of the X. laevis genes; apparently there are similar constraints on both copies. No correlation was found between dS for the X. laevis genes and dS for the mammalian genes. There was a significant positive correlation both between members of pairs of duplicated X. laevis genes (r = 0.951) and between human and rodent orthologues (r = 0.854) with respect to third- position G+C content but no such relationship between the X. laevis genes and either of their mammalian orthologues. The results indicate that both copies of a duplicate gene can be subject to purifying selection and thus support the hypothesis of selection against all genotypes containing a null allele at either of two duplicate loci.      

Positive selection during the diversification of class I vomeronasal receptor-like (V1RL) genes,putative pheromone receptor genes,in human and primate evolution

Vomeronasal receptors are the major receptors for pheromones in vertebrates, and five putative type 1 vomeronasal receptors (V1RL) have been identified in humans. The evolution of the V1RL1 gene in non-human primates, and patterns of selection on V1RL genes, were investigated. The presumed ortholog of V1RL1 was sequenced from 13 species of nonhuman primate, and in eight of these species V1RL1 was a pseudogene. Phylogenetic reconstructions reveal that V1RL1 pseudogene formation occurred independently in multiple primate lineages. Using maximum likelihood estimates of dN/dS ratios in PAML, we show that V1RL genes have evolved under neutral evolution in lineages in which they became a pseudogene. In contrast, among lineages in which V1RL genes contain an open reading frame, the majority of sites are under purifying selection and a minority are under significant positive selection. These results provide an interesting case where all three categories of selection can be teased apart in the same data set using maximum likelihood methods. The finding of positive selection on V1RL genes during primate evolution provides indirect support for the hypothesis that V1RL genes have a function in species-specific pheromone detection in primates.     

A Sliding Window-Based Method to Detect Selective Constraints in Protein-Coding Genes and Its Application to RNA Viruses

Here we present a new sliding window-based method specially designed to detect selective constraints in specific regions of a multiple protein-coding sequence alignment. In contrast to previous window-based procedures, our method is based on a nonarbitrary statistical approach to find the appropriate codon-window size to test deviations of synonymous (dS) and nonsynonymous (dN) nucleotide substitutions from the expectation. The probabilities of dN and dS are obtained from simulated data and used to detect significant deviations of dN and dS in a specific window region of the real sequence alignment. The nonsynonymous-to-synonymous rate ratio (w = dN/dS) was used to highlight selective constraints in any window wherein dS or dN was significantly different from the expectation. In these significant windows, w and its variance [V(w)] were calculated and used to test the neutral hypothesis. Computer simulations showed that the method is accurate even for highly divergent sequences. The main advantages of the new method are that it (i) uses a statistically appropriate window size to detect different selective patterns, (ii) is computationally less intensive than maximum likelihood methods, and (iii) detects saturation of synonymous sites, which can give deviations from neutrality. Hence, it allows the analysis of highly divergent sequences and the test of different alternative hypothesis as well. The application of the method to different human immunodeficiency virus type 1 and to foot-and-mouth disease virus genes confirms the action of positive selection on previously described regions as well as on new regions.     

Two mitochondrial genes under episodic positive selection in subterranean octodontoid rodents

Tuco-tucos (Ctenomys) and related coruros (Spalacopus) are South American subterranean rodents. An energetically demanding lifestyle within the hypoxic, underground atmosphere may change the selective regime on oxidative phosphorylation. We examined whether weak and/or episodic positive directional selection affected the evolution of two mitochondrial genes (COX2, CytB), in a background of purifying selection in these lineages. We estimated rates of synonymous (dS) and non-synonymous (dN) substitutions and found: 1) significantly higher dN/dS ratio in subterranean groups relative to non-subterranean related species, and 2) two codons in each gene under episodic selection: 94 and 277 of COX2 and 269 and 307 of CytB.     

Positive darwinian selection observed at the variable-region genes of immunoglobulins

To examine Gojobori and Nei's hypothesis that the immunoglobulin heavy- chain variable-region (VH) genes in mammals are subject to diversity- enhancing selection, we studied the rates of synonymous and nonsynonymous nucleotide substitution in the complementarity- determining regions (CDRs) and in the framework regions (FRs) of mouse and human VH genes. The results obtained indicate that the non- synonymous rate is higher than the synonymous rate in CDRs, whereas the reverse is true in FRs. This observation supports Gojobori and Nei's hypothesis and suggests that diversity-enhancing selection (similar to overdominant selection) operates mainly in CDRs and that this is one of the evolutionary factors that increase antibody diversity.      

The Population Genetics of dN/dS

Evolutionary pressures on proteins are often quantified by the ratio of substitution rates at non-synonymous and synonymous sites. The dN/dS ratio was originally developed for application to distantly diverged sequences, the differences among which represent substitutions that have fixed along independent lineages. Nevertheless, the dN/dS measure is often applied to sequences sampled from a single population, the differences among which represent segregating polymorphisms. Here, we study the expected dN/dS ratio for samples drawn from a single population under selection, and we find that in this context, dN/dS is relatively insensitive to the selection coefficient. Moreover, the hallmark signature of positive selection over divergent lineages, dN/dS>1, is violated within a population. For population samples, the relationship between selection and dN/dS does not follow a monotonic function, and so it may be impossible to infer selection pressures from dN/dS. These results have significant implications for the interpretation of dN/dS measurements among population-genetic samples.     

Evolution of immunoglobulin VH pseudogenes in chickens

In chickens, there is a single functional gene (VH1) coding for the heavy chain variable region of immunoglobulins, and immunoglobulin diversity is generated by gene conversion of the VH1 gene by many variable region pseudogenes (psi VH's) that exist on the 5' side of the VH1 gene. To understand the evolution of this unique genetic system, we conducted statistical analyses of VH1 and psi VH genes together with functional VH genes from other higher vertebrate species. The results indicate, first, that chicken VH genes are all closely related to one another and were derived relatively recently from an ancestral gene belonging to one of the three major groups of VH genes in higher vertebrates. Second, the rate of nonsynonymous substitution is slightly higher than that of synonymous substitution in the complementarity- determining regions (CDRs), which suggests that diversity-enhancing selection has operated in the CDRs even for pseudogenes. However, both the rates of synonymous and nonsynonymous substitution are higher in the CDRs than in the framework regions (FRs), apparently because of an interaction between positive selection and meiotic gene conversion in the CDRs. Third, a dot matrix analysis of the psi VH genes and genomic diversity (D) genes has indicated that the 3' end of psi VH genes is attached by D-gene-like sequences, and this region of psi VH genes has high similarity with D gene sequences. This suggests that V and D genes were fused at some point of evolutionary time and this fused element multiplied by gene duplication. Finally, two alternative hypotheses of explaining the evolution of the chicken VH gene system are presented.      

Identification of fast-evolving genes in the scleractinian coral Acropora using comparative EST analysis

To identify fast-evolving genes in reef-building corals, we performed direct comparative sequence analysis with expressed sequence tag (EST) datasets from two acroporid species: Acropora palmata from the Caribbean Sea and A. millepora from the Great Barrier Reef in Australia. Comparison of 589 independent sequences from 1,421 A. palmata contigs, with 10,247 A. millepora contigs resulted in the identification of 196 putative homologues. Most of the homologous pairs demonstrated high amino acid similarities (over 90%). Comparisons of putative homologues showing low amino acid similarities (under 90%) among the Acropora species to the near complete datasets from two other cnidarians (Hydra magnipapillata and Nematostella vectensis) implied that some were non-orthologous. Within 86 homologous pairs, 39 exhibited dN/dS ratios significantly less than 1, suggesting that these genes are under purifying selection associated with functional constraints. Eight independent genes showed dN/dS ratios exceeding 1, while three deviated significantly from 1, suggesting that these genes may play important roles in the adaptive evolution of Acropora. Our results also indicated that CEL-III lectin was under positive selection, consistent with a possible role in immunity or symbiont recognition. Further studies are needed to clarify the possible functions of the genes under positive selection to provide insight into the evolutionary process of corals.     

Selection in coastal Synechococcus (cyanobacteria) populations evaluated from environmental metagenomes

Environmental metagenomics provides snippets of genomic sequences from all organisms in an environmental sample and are an unprecedented resource of information for investigating microbial population genetics. Current analytical methods, however, are poorly equipped to handle metagenomic data, particularly of short, unlinked sequences. A custom analytical pipeline was developed to calculate dN/dS ratios, a common metric to evaluate the role of selection in the evolution of a gene, from environmental metagenomes sequenced using 454 technology of flow-sorted populations of marine Synechococcus, the dominant cyanobacteria in coastal environments. The large majority of genes (98%) have evolved under purifying selection (dN/dS<1). The metagenome sequence coverage of the reference genomes was not uniform and genes that were highly represented in the environment (i.e. high read coverage) tended to be more evolutionarily conserved. Of the genes that may have evolved under positive selection (dN/dS>1), 77 out of 83 (93%) were hypothetical. Notable among annotated genes, ribosomal protein L35 appears to be under positive selection in one Synechococcus population. Other annotated genes, in particular a possible porin, a large-conductance mechanosensitive channel, an ATP binding component of an ABC transporter, and a homologue of a pilus retraction protein had regions of the gene with elevated dN/dS. With the increasing use of next-generation sequencing in metagenomic investigations of microbial diversity and ecology, analytical methods need to accommodate the peculiarities of these data streams. By developing a means to analyze population diversity data from these environmental metagenomes, we have provided the first insight into the role of selection in the evolution of Synechococcus, a globally significant primary producer.     

The evolution of mitochondrial genomes in subterranean caviomorph rodents: adaptation against a background of purifying selection

South American tuco-tucos (Ctenomys) and the related coruro (Spalacopus) are two rodent lineages that have independently colonised the subterranean niche. The energetically demanding lifestyles of these species, coupled with the hypoxic atmospheres characteristic of subterranean environments, may have altered the selective regimes on genes encoding proteins related to cellular respiration. Here, we examined the molecular evolution of 13 protein-coding genes in the mitochondrial genome of seven caviomorph rodents, including these two subterranean genera and their above-ground relatives. Using maximum-likelihood and Bayesian approaches, we estimated rates of synonymous (dS) and nonsynonymous (dN) substitutions. We found a significantly higher ω ratio (dN/dS) in subterranean groups as compared to their non-subterranean counterparts in 11 of 13 genes, although no ω ratio was larger than 1. Additionally, we applied a method based on quantitative physicochemical properties to test for positive selection. Amino acid changes implicated in radical structural or functional shifts in the protein property were found to be ubiquitous across the phylogeny, but concentrated in the subterranean lineages. Convergent changes were also found between the subterranean genera used in this study and other mammals adapted to hypoxia. The results of this study suggest a link between niche shifts and weak directional (or episodic) selection at the molecular level against a background of purifying selection.     

Lineage‐specific sequence evolution and exon edge conservation partially explain the relationship between evolutionary rate and expression level in A. thaliana

Rapidly evolving proteins can aid the identification of genes underlying phenotypic adaptation across taxa, but functional and structural elements of genes can also affect evolutionary rates. In plants, the ‘edges’ of exons, flanking intron junctions, are known to contain splice enhancers and to have a higher degree of conservation compared to the remainder of the coding region. However, the extent to which these regions may be masking indicators of positive selection or account for the relationship between dN/dS and other genomic parameters is unclear. We investigate the effects of exon edge conservation on the relationship of dN/dS to various sequence characteristics and gene expression parameters in the model plant Arabidopsis thaliana. We also obtain lineage‐specific dN/dS estimates, making use of the recently sequenced genome of Thellungiella parvula, the second closest sequenced relative after the sister species Arabidopsis lyrata. Overall, we find that the effect of exon edge conservation, as well as the use of lineage‐specific substitution estimates, upon dN/dS ratios partly explains the relationship between the rates of protein evolution and expression level. Furthermore, the removal of exon edges shifts dN/dS estimates upwards, increasing the proportion of genes potentially under adaptive selection. We conclude that lineage‐specific substitutions and exon edge conservation have an important effect on dN/dS ratios and should be considered when assessing their relationship with other genomic parameters.     

Vertebrate extracellular calcium-sensing receptor evolution: Selection in relation to life history and habitat

Ionic calcium (Ca^2 +) supports essential functions within physiological systems, and consequently its concentration is homeostatically regulated within narrow bounds in the body fluids of animals through endocrine effects at ion-transporting osmoregulatory tissues. In vertebrates, extracellular Ca^2 + is detected at the cell surface by the extracellular calcium-sensing receptor (CaSR), a member of the G protein-coupled receptor (GPCR) superfamily. Interestingly, the taxonomic distribution of CaSRs is restricted to vertebrates, with some CaSR-like receptors apparently present in non-vertebrate chordates. Since bone is a known Ca^2 + storage site and is characteristically restricted to the vertebrate lineage, we hypothesized a functional association of CaSR with vertebrate skeleton that may have an ancient origin. Protein sequence alignment and phylogenetic analysis of vertebrate CaSRs and related GPCRs of the glutamate receptor-like family expose similarities and indel differences among these receptors, and reveal the evolutionary history of CaSRs. Evolutionary selection was tested statistically by evaluating the relationship between non-synonymous (replacement, dN) versus synonymous (silent, dS) amino acid substitution rates (as dN/dS) of protein-coding DNA sequences among branches of the estimated protein phylogeny. On a background of strong purifying selection (dN/dS < 1) in the CaSR phylogeny, statistical evidence for adaptive evolution (dN/dS > 1) was detected on some branches to major clades in the CaSR phylogeny, especially to the tetrapod vertebrate CaSRs and chordate CaSR-like branches. Testing also revealed overall purifying selection at the codon level. At some sites relaxation from strong purifying selection was seen, but evidence for adaptive evolution was not detected for individual sites. The results suggest purifying selection of CaSRs, and of adaptive evolution among some major vertebrate clades, reflecting clade specific differences in natural history and organismal biology, including skeletal involvement in calcium homeostasis.     

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.     

Likelihood ratio tests for detecting positive selection and application to primate lysozyme evolution

An excess of nonsynonymous substitutions over synonymous ones is an important indicator of positive selection at the molecular level. A lineage that underwent Darwinian selection may have a nonsynonymous/synonymous rate ratio (dN/dS) that is different from those of other lineages or greater than one. In this paper, several codon-based likelihood models that allow for variable dN/dS ratios among lineages were developed. They were then used to construct likelihood ratio tests to examine whether the dN/dS ratio is variable among evolutionary lineages, whether the ratio for a few lineages of interest is different from the background ratio for other lineages in the phylogeny, and whether the dN/dS ratio for the lineages of interest is greater than one. The tests were applied to the lysozyme genes of 24 primate species. The dN/dS ratios were found to differ significantly among lineages, indicating that the evolution of primate lysozymes is episodic, which is incompatible with the neutral theory. Maximum- likelihood estimates of parameters suggested that about nine nonsynonymous and zero synonymous nucleotide substitutions occurred in the lineage leading to hominoids, and the dN/dS ratio for that lineage is significantly greater than one. The corresponding estimates for the lineage ancestral to colobine monkeys were nine and one, and the dN/dS ratio for the lineage is not significantly greater than one, although it is significantly higher than the background ratio. The likelihood analysis thus confirmed most, but not all, conclusions Messier and Stewart reached using reconstructed ancestral sequences to estimate synonymous and nonsynonymous rates for different lineages.      

Pattern of the evolution of HIV-1 enν gene in C?te d?Ivoire

Cête d׳Ivoire continues to have the highest HIV-1 prevalence rate in West Africa, although the infection number is in constant decline. The external envelope protein of the viruses is a likely site of selection, and responsible for receptor binding and entry into host cells, and therefore constitutes an ideal region with which to investigate the evolutionary processes acting on HIV-1. In this study, we analyse 189 envelope glycoprotein V3 loop region sequences of viruse isolates from 1995 to 2009, from HIV-1 untreated patients living in Cête d׳Ivoire, to decipher the temporal relationship between disease diversity, divergence and selection. Our analyses show that the nonsynonymous and synonymous ratio (dN/dS) was lower than 1 for viral populations analysed within 15 years, which showed the sequences did not undergo adequate immune pressure. The phylogenetic tree of the sequences analysed demonstrated distinctly long internal branches and short external branches, suggesting that only a small number of viruses infected the new host cell at each transmission. In addition to identifying sites under purifying selection, we also identified neutral sites that can cause false positive inference of selection. These sites presented form a resource for future studies of selection pressures acting on HIV-1 enν gene in Cête d׳Ivoire and other West African countries.