Evolution of MHC class IIB in the genome of wild and ornamental guppies, Poecilia reticulata

This is the first study to quantify genomic sequence variation of the major histocompatibility complex (MHC) in wild and ornamental guppies, Poecilia reticulata. We sequenced 196-219 bp of exon 2 MHC class IIB (DAB) in 56 wild Trinidadian guppies and 14 ornamental strain guppies. Each of two natural populations possessed high allelic richness (15-16 alleles), whereas only three or fewer DAB alleles were amplified from ornamental guppies. The disparity in allelic richness between wild and ornamental fish cannot be fully explained by fixation of alleles by inbreeding, nor by the presence of non-amplified sequences (ie null alleles). Rather, we suggest that the same allele is fixed at duplicated MHC DAB loci owing to gene conversion. Alternatively, the number of loci in the ornamental strains has contracted during >100 generations in captivity, a hypothesis consistent with the accordion model of MHC evolution. We furthermore analysed the substitution patterns by making pairwise comparisons of sequence variation at the putative peptide binding region (PBR). The rate of non-synonymous substitutions (dN) only marginally exceeded synonymous substitutions (dS) in PBR codons. Highly diverged sequences showed no evidence for diversifying selection, possibly because synonymous substitutions have accumulated since their divergence. Also, the substitution pattern of similar alleles did not show evidence for diversifying selection, plausibly because advantageous non-synonymous substitutions have not yet accumulated. Intermediately diverged sequences showed the highest relative rate of non-synonymous substitutions, with dN/dS>14 in some pairwise comparisons. Consequently, a curvilinear relationship was observed between the dN/dS ratio and the level of sequence divergence.     

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.     

Apparent Selection Intensity for the Cytochrome Oxidase Subunit I Gene Varies with Mode of Reproduction in Echinoderms

When most amino acid substitutions in protein-coding genes are slightly deleterious rather than selectively neutral, life history differences can potentially modify the effective population size or the selective regime, resulting in altered ratios of non-synonymous to synonymous substitutions among taxa. We studied substitution patterns for the mitochondrial cytochrome oxidase subunit I (COI) gene in a sea star genus (Leptasterias spp.) with an obligate brood-protecting mode of reproduction and small-scale population genetic subdivision, and compared the results to available COI sequences in nine other genera of echinoderms with pelagic larvae: three sea stars, five sea urchins and one brittle star. We predicted that this life history difference would be associated with differences in the ratio of non-synonymous (dN) to synonymous (dS) substitution rates. Leptasterias had a significantly greater dN/dS ratio (both between species and within species), a significantly smaller transition/transversion rate ratio, and a significantly lower average nucleotide diversity within species, than did the non-brooding genera. Other explanations for the results, such as altered mutation rates or selective sweeps, were not supported by the data analysis. These findings highlight the potential influence of reproductive traits and other life history factors on patterns of nucleotide substitution within and between species.     

Gene expression,chromosome heterogeneity and the fast-X effect in mammals

The higher rate of non-synonymous over synonymous substitutions (dN/dS) of the X chromosome compared with autosomes is often interpreted as a consequence of X hemizygosity. However, other factors, such as gene expression, are also known to vary between X and autosomes. Analysing 4800 orthologues in six mammals, we found that gene expression levels, associated with GC content, fully account for the variation in dN/dS between X and autosomes with no detectable effect of hemizygosity. We also report an extensive variance in dN/dS and gene expression between autosomes.     

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.     

Estimating diversifying selection and functional constraint in the presence of recombination

Models of molecular evolution that incorporate the ratio of nonsynonymous to synonymous polymorphism (dN/dS ratio) as a parameter can be used to identify sites that are under diversifying selection or functional constraint in a sample of gene sequences. However, when there has been recombination in the evolutionary history of the sequences, reconstructing a single phylogenetic tree is not appropriate, and inference based on a single tree can give misleading results. In the presence of high levels of recombination, the identification of sites experiencing diversifying selection can suffer from a false-positive rate as high as 90%. We present a model that uses a population genetics approximation to the coalescent with recombination and use reversible-jump MCMC to perform Bayesian inference on both the dN/dS ratio and the recombination rate, allowing each to vary along the sequence. We demonstrate that the method has the power to detect variation in the dN/dS ratio and the recombination rate and does not suffer from a high false-positive rate. We use the method to analyze the porB gene of Neisseria meningitidis and verify the inferences using prior sensitivity analysis and model criticism techniques.     

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.      

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.     

Phylogenetic analysis suggests that sociality is associated with reduced effectiveness of selection

The evolution of sociality in spiders is associated with female bias, reproductive skew and an inbreeding mating system, factors that cause a reduction in effective population size and increase effects of genetic drift. These factors act to decrease the effectiveness of selection, thereby increasing the fixation probability of deleterious mutations. Comparative studies of closely related species with contrasting social traits and mating systems provide the opportunity to test consequences of low effective population size on the effectiveness of selection empirically. We used phylogenetic analyses of three inbred social spider species and seven outcrossing subsocial species of the genus Stegodyphus, and compared dN/dS ratios and codon usage bias between social Inbreeding and subsocial outcrossing mating systems to assess the effectiveness of selection. The overall results do not differ significantly between the social inbreeding and outcrossing species, but suggest a tendency for lower codon usage bias and higher dN/dS ratios in the social inbreeding species compared with their outcrossing congeners. The differences in dN/dS ratio and codon usage bias between social and subsocial species are modest but consistent with theoretical expectations of reduced effectiveness of selection in species with relatively low effective population size. The modest differences are consistent with relatively recent evolution of social mating systems. Additionally, the short terminal branches and lack of speciation of the social lineages, together with low genetic diversity lend support for the transient state of permanent sociality in spiders.     

The impact of recombination on dN/dS within recently emerged bacterial clones

The development of next-generation sequencing platforms is set to reveal an unprecedented level of detail on short-term molecular evolutionary processes in bacteria. Here we re-analyse genome-wide single nucleotide polymorphism (SNP) datasets for recently emerged clones of methicillin resistant Staphylococcus aureus (MRSA) and Clostridium difficile. We note a highly significant enrichment of synonymous SNPs in those genes which have been affected by recombination, i.e. those genes on mobile elements designated "non-core" (in the case of S. aureus), or those core genes which have been affected by homologous replacements (S. aureus and C. difficile). This observation suggests that the previously documented decrease in dN/dS over time in bacteria applies not only to genomes of differing levels of divergence overall, but also to horizontally acquired genes of differing levels of divergence within a single genome. We also consider the role of increased drift acting on recently emerged, highly specialised clones, and the impact of recombination on selection at linked sites. This work has implications for a wide range of genomic analyses.     

Preponderance of slightly deleterious polymorphism in mitochondrial DNA: nonsynonymous/synonymous rate ratio is much higher within species than between species

We estimated synonymous (dN) and nonsynonymous (dS) substitution rates for protein-coding genes of the mitochondrial genome from two individuals each of the species human, chimpanzee, and gorilla. The genes were analyzed both separately and in a combined data set. Pairwise sequence comparisons suggest that the dN/dS rate ratios are about 5-10 times higher in within-species comparisons than in between-species comparisons. This result is confirmed by a more rigorous likelihood ratio test, which rejected the null hypothesis that the dN/dS rate ratios are identical within and between species. The likelihood models account for the genetic code structure, transition/transversion rate ratio, and codon usage bias and are expected to produce more reliable results than the commonly used contingency test. Separate analyses of different genes show that the dN/dS rate ratios are higher within species than between species for all 13 mitochondrial genes, with the difference being statistically significant for all except three small or slowly evolving genes. Furthermore, in conserved genes, nonsynonymous rates within species tend to be higher than the between-species rates by a greater proportion than in fast-changing genes. Our findings confirm and extend earlier results obtained from smaller data sets and suggest the operation of slightly deleterious mutations throughout the mitochondrial genome in the hominoids. Implications of the results for evolutionary studies and, in particular, for studies of the origin of modern humans, are discussed.     

Selection on protein-coding genes of natural cyanobacterial populations

We examined the distribution of synonymous and non-synonymous changes in 12 protein-coding genes of natural populations of cyanobacteria to infer changes in gene functionality. By comparing mutation distributions within and across species using the McDonald-Kreitman test, we found data sets to contain evidence for purifying selection (hetR of Trichodesmium, nifH of Cylindrospermopsis raceborskii and rpoC1 of Anabaena lemmermannii) and positive selection (kaiC of Microcoleus chthonoplastes and rbcX of Anabaena and Aphanizomenon sp.). Other genes from the same set of clonal isolates (petB and rbcL in M. chthonoplastes and Anabaena/Aphanizomenon, respectively) did not harbour evidence for either form of selection. The results of branch models of codon evolution agreed fully with the results of the McDonald-Kreitman test in terms of significance and absolute value of the dN/dS estimates. The high frequency of gene-specific mutation patterns and their association with branches that separate closely related cyanobacterial genera suggest that evolutionary tests are suited to uncover gene-specific selective differentiation in cyanobacterial genomes. At the same time, given the lack of information about the history of cyanobacteria, analysis of larger numbers of protein-coding genes of clonal cyanobacterial isolates will produce more detailed pictures of the effects of natural selection.     

Physicochemical evolution and molecular adaptation of the cetacean and artiodactyl cytochrome b proteins

Cetaceans have most likely experienced metabolic shifts since evolutionarily diverging from their terrestrial ancestors, shifts that may be reflected in the proteins such as cytochrome b that are responsible for metabolic efficiency. However, accepted statistical methods for detecting molecular adaptation are largely biased against even moderately conservative proteins because the primary criterion involves a comparison of nonsynonymous and synonymous substitution rates (dN/dS); they do not allow for the possibility that adaptation may come in the form of very few amino acid changes. We apply the MM01 model to the possible molecular adaptation of cytochrome b among cetaceans because it does not rely on a dN/dS ratio, instead evaluating positive selection in terms of the amino acid properties that comprise protein phenotypes that selection at the molecular level may act upon. We also apply the codon-degeneracy model (CDM), which focuses on evaluating overall patterns of nucleotide substitution in terms of base exchange, codon position, and synonymy to estimate the overall effect of selection. Using these relatively new models, we characterize the molecular adaptation that has occurred in the cetacean cytochrome b protein by comparing revealed amino acid replacement patterns to those found among artiodactyls, the modern terrestrial mammals found to be most closely related to cetaceans. Our findings suggest that several regions of the cetacean cytochrome b protein have experienced molecular adaptation. Also, these adaptations are spatially associated with domain structure, protein function, and the structure and function of the cytochrome bc(1) complex and its constituents. We also have found a general correlation between the results of the analytical software programs TreeSAAP (which implements the MM01 model) and CDM (which implements the codon-degeneracy model).     

Cytochrome b evolution in birds and mammals: an evaluation of the avian constraint hypothesis.

Patterns of molecular evolution in birds have long been considered anomalous. Compared with other vertebrates, birds have reduced levels of genetic divergence between groups of similar taxonomic ranks for a variety of nuclear and mitochondrial markers. This observation led to the avian constraint hypothesis, which identifies increased functional constraint on avian proteins as the cause for the reduction in genetic divergence. Subsequent investigations provided additional support for the avian constraint hypothesis when rates of molecular evolution were found to be slower in birds than in mammals in a variety of independent calibrations. It is possible to test the avian constraint hypothesis as an explanation for this avian slowdown by comparing DNA sequence data from protein-coding regions in birds and homologous regions in mammals. The increased selective constraints should lead to a reduction in the proportion of amino acid replacement substitutions. To test for such a decrease, we calculated the numbers of amino acid replacement substitutions per replacement site (dN) and silent substitutions per silent site (dS) for the complete mitochondrial cytochrome b gene using 38 avian and 43 mammalian comparisons that were phylogenetically independent. We find that dN/dS is significantly smaller in birds than in mammals. This difference cannot be explained by differences in codon bias affecting dS values. We suggest that the avian slowdown can be explained, at least in part, by a decreased tolerance for amino acid substitutions in avian species relative to mammalian species.     

Independent effects of alternative splicing and structural constraint on the evolution of mammalian coding exons

Alternative splicing (AS) is known to significantly affect exon-level protein evolutionary rates in mammals. Particularly, alternatively spliced exons (ASEs) have a higher nonsynonymous-to-synonymous substitution rate (dN/dS) ratio than constitutively spliced exons (CSEs), possibly because the former are required only occasionally for normal biological functions. Meanwhile, intrinsically disordered regions (IDRs), the protein regions lacking fixed 3D structures, are also reported to have an increased evolutionary rate due to lack of structural constraint. Interestingly, IDRs tend to be located in alternative protein regions. Yet which of these two factors is the major determinant of the increased dN/dS in mammalian ASEs remains unclear. By comparing human-macaque and human-mouse one-to-one orthologous genes, we demonstrate that AS and protein structural disorder have independent effects on mammalian exon evolution. We performed analyses of covariance to demonstrate that the slopes of the (dN/dS-percentage of IDR) regression lines differ significantly between CSEs and ASEs. In other words, the dN/dS ratios of both ASEs and CSEs increase with the proportion of IDR (PIDR), whereas ASEs have higher dN/dS ratios than CSEs when they have similar PIDRs. Since ASEs and IDRs may less frequently overlap with protein domains (which also affect dN/dS), we also examined the correlations between dN/dS ratio and exon type/PIDR by controlling for the density of protein domain. We found that the effects of exon type and PIDR on dN/dS are both independent of domain density. Our results imply that nature can select for different biological features with regard to ASEs and IDRs, even though the two biological features tend to be localized in the same protein regions.     

Molecular evolution of hepatitis B virus over 25 years

Determining the longitudinal molecular evolution of hepatitis B virus (HBV) is difficult due to HBV's genomic complexity and the need to study paired samples collected over long periods of time. In this study, serial samples were collected from eight hepatitis B virus e antigen-negative asymptomatic carriers of HBV genotype B in 1979 and 2004, thus providing a 25-year period to document the long-term molecular evolution of HBV. The rate, nature, and distribution of mutations that emerged over 25 years were determined by phylogenetic and linear regression analysis of full-length HBV genome sequences. Nucleotide hypervariability was observed within the polymerase and pre-S/S overlap region and within the core gene. The calculated mean number of nucleotide substitutions/site/year (7.9 x 10(-5)) was slightly higher than published estimates (1.5 x 10(-5) to 5 x 10(-5)). Nucleotide changes in the quasispecies population did not significantly alter the molecular evolutionary rate, based on linear regression analysis of evolutionary distances among serial clone pre-S region sequences. Therefore, the directly amplified or dominant sequence was sufficient to estimate the putative molecular evolutionary rate for these long-term serial samples. On average, the ratio of synonymous (dS) to nonsynonymous (dN) substitutions was highest for the polymerase-coding region and lowest for the core-coding region. The low dS/dN ratios observed within the core suggest that selection occurs within this gene region, possibly as an immune evasion strategy. The results of this study suggest that HBV sequence divergence may occur more rapidly than previously estimated, in a host immune phase-dependent manner.     

Detection of site-specific positive Darwinian selection on pandemic influenza A/H1N1 virus genome: integrative approaches

In the twenty-first century, the first pandemic novel human influenza A/H1N1virus (NIV) outbreak was reported at Mexico and USA on March and early April, 2009 respectively. The outbreak occurred among human populations due to the presence of meager or no immune response against newly emerged viruses. The success of vaccines and drugs depends on their low susceptibility to the formation of escape mutants in virus. Identification of excess, non-synonymous substitutions over synonymous ones is a main indicator of positive Darwinian selection in protein-coding genes of NIVs. The positive Darwinian selection operating on each site of proteins were inferred by computing ω, the ratio of the non-synonymous/synonymous substitutions [dN/dS (or) K_a/K_s], which was calculated by three different methods in terms of codon-based maximum likelihood, branch-site and empirical Bayesian methods under various models. Totally, nine sites from PB2, PB1, HA, M2 and NS1 are inferred as positively selected. The function for amino acid sites of NIVs proteins under positive selection are inferred by comparing the sites with experimentally determined functionally known amino acid sites. Completely 4 positively selected sites of PB1, HA and M2 are found to be involved in B-cell epitopes (BCEs). Interestingly, most of these sites are also involving in T-cell epitopes (TCEs). However, more sites under positive selection forces are involved in TCEs than those of BCEs. Amino acid sites engaged in both BCEs and TCEs should be measured as highly suitable targets, because these sites could induce the strong humoral and cellular immune responses against targets.     

Evolution of mitochondrial-encoded cytochrome oxidase subunits in endothermic fish: the importance of taxon-sampling in codon-based models

While endothermy is ubiquitous in birds and mammals, it is not exclusive to these most recently arisen vertebrate classes. The ability to warm specific organs and/or tissues above ambient temperature (regional endothermy) has evolved at least three times in phylogentically discrete fish lineages: lamnid sharks (Lamnidae), tunas (Scombridae) and billfishes (Istiophoridae and Xiphidae). Given the links between endothermy and metabolic rate, we looked for evidence of convergent molecular evolution in mtDNA-encoded cytochrome c oxidase (COX) subunits in each of these discrete lineages. We found no evidence that the endothermic phenotype in fishes is driven or accompanied by molecular convergence. Though we found little evidence for positively-selected sites in any of the lineages in any subunit, the conclusions were sensitive to the choice of maximum-likelihood model. Several sites identified by Na?ve Empirical Bayes (NEB) were not found when Bayes Empirical Bayes (BEB) was employed. As well, conclusions were profoundly influenced by taxon-sampling. Several of the putative sites of positive selection in COX II were no longer apparent as we augmented taxon sampling. The lack of convergent molecular evolution in these remarkable taxa, combined with the profound influence of model choice and taxon sampling provide a cautionary note on the use of rates of non-synonymous to synonymous mutations (dN/dS) to explore questions of the evolution of physiological function.     

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.