Functional Divergence and Evolutionary Turnover in Mammalian Phosphoproteomes

Protein phosphorylation is a key mechanism to regulate protein functions. However, the contribution of this protein modification to species divergence is still largely unknown. Here, we studied the evolution of mammalian phosphoregulation by comparing the human and mouse phosphoproteomes. We found that 84% of the positions that are phosphorylated in one species or the other are conserved at the residue level. Twenty percent of these conserved sites are phosphorylated in both species. This proportion is 2.5 times more than expected by chance alone, suggesting that purifying selection is preserving phosphoregulation. However, we show that the majority of the sites that are conserved at the residue level are differentially phosphorylated between species. These sites likely result from false-negative identifications due to incomplete experimental coverage, false-positive identifications and non-functional sites. In addition, our results suggest that at least 5% of them are likely to be true differentially phosphorylated sites and may thus contribute to the divergence in phosphorylation networks between mouse and humans and this, despite residue conservation between orthologous proteins. We also showed that evolutionary turnover of phosphosites at adjacent positions (in a distance range of up to 40 amino acids) in human or mouse leads to an over estimation of the divergence in phosphoregulation between these two species. These sites tend to be phosphorylated by the same kinases, supporting the hypothesis that they are functionally redundant. Our results support the hypothesis that the evolutionary turnover of phosphorylation sites contributes to the divergence in phosphorylation profiles while preserving phosphoregulation. Overall, our study provides advanced analyses of mammalian phosphoproteomes and a framework for the study of their contribution to phenotypic evolution.     

Comparison of diverse protein sequences of the nuclear-encoded subunits of cytochrome C oxidase suggests conservation of structure underlies evolving functional sites

Interspecific comparisons of protein sequences can reveal regions of evolutionary conservation that are under purifying selection because of functional constraints. Interpreting these constraints requires combining evolutionary information with structural, biochemical, and physiological data to understand the biological function of conserved regions. We take this integrative approach to investigate the evolution and function of the nuclear-encoded subunits of cytochrome c oxidase (COX). We find that the nuclear-encoded subunits evolved subsequent to the origin of mitochondria and the subunit composition of the holoenzyme varies across diverse taxa that include animals, yeasts, and plants. By mapping conserved amino acids onto the crystal structure of bovine COX, we show that conserved residues are structurally organized into functional domains. These domains correspond to some known functional sites as well as to other uncharacterized regions. We find that amino acids that are important for structural stability are conserved at frequencies higher than expected within each taxon, and groups of conserved residues cluster together at distances of less than 5 A more frequently than do randomly selected residues. We, therefore, suggest that selection is acting to maintain the structural foundation of COX across taxa, whereas active sites vary or coevolve within lineages.     

Evidence that purifying selection acts on promoter sequences

We tested whether functionally important sites in bacterial, yeast, and animal promoters are more conserved than their neighbors. We found that substitutions are predominantly seen in less important sites and that those that occurred tended to have less impact on gene expression than possible alternatives. These results suggest that purifying selection operates on promoter sequences.     

Evidence for heterogeneous selective pressures in the evolution of the env gene in different human immunodeficiency virus type 1 subtypes

Recent studies have demonstrated the emergence of human immunodeficiency virus type 1 (HIV-1) subtypes with various levels of fitness. Using heterogeneous maximum-likelihood models of adaptive evolution implemented in the PAML software package, with env sequences representing each HIV-1 group M subtype, we examined the various intersubtype selective pressures operating across the env gene. We found heterogeneity of evolutionary mechanisms between the different subtypes with a category of amino acid sites observed that had undergone positive selection for subtypes C, F1, and G, while these sites had undergone purifying selection in all other subtypes. Also, amino acid sites within subtypes A and K that had undergone purifying selection were observed, while these sites had undergone positive selection in all other subtypes. The presence of such sites indicates heterogeneity of selective pressures within HIV-1 group M subtype evolution that may account for the various levels of fitness of the subtypes.     

Episodic evolution of growth hormone in primates and emergence of the species specificity of human growth hormone receptor

Growth hormone (GH) evolution is very conservative among mammals, except for primates and ruminant artiodactyls. In fact, most known mammalian GH sequences differ from the inferred ancestral mammalian sequence by only a few amino acids. In contrast, the human GH sequence differs from the inferred ancestral sequence by 59 amino acids. However, it is not known when this rapid evolution of GH occurred during primate evolution or whether it was due to positive selection. Also, human growth hormone receptor (GHR) displays species specificity; i.e., it can interact only with human (or rhesus monkey) GH, not with nonprimate GHS: The species specificity of human GHR is largely due to the Leu-->Arg change at position 43, and it has been hypothesized that this change must have been preceded by the His-->Asp change at position 171 of GH. Is this hypothesis true? And when did these changes occur? To address the above issues, we sequenced GH and GHR genes in prosimians and simians. Our data supported the above hypothesis and revealed that the species specificity of human GHR actually emerged in the common ancestor of Old World primates, but the transitional phase still persists in New World monkeys. Our data showed that the rapid evolution of primate GH occurred during a relatively short period (in the common ancestor of higher primates) and that the rate of change was especially high at functionally important sites, suggesting positive selection. However, the nonsynonymous rate/synonymous rate ratio at these sites was <1, so relaxation of purifying selection might have played a role in the rapid evolution of the GH gene in simians, possibly as a result of multiple gene duplications. Similar to GH, GHR displayed an accelerated rate of evolution in primates. Our data revealed proportionally more amino acid replacements at the functionally important sites in both GH and GHR in simians but, surprisingly, showed few coincidental replacements of amino acids forming the same intermolecular contacts between the two proteins.     

Phosphoglucose isomerases of hagfish, zebrafish, gray mullet, toad, and snake, with reference to the evolution of the genes in vertebrates

Phosphoglucose isomerase (PGI) is a protein with multiple functions. To infer its structure changes and evolution in vertebrates, we cloned cDNAs encoding PGI genes from hagfish (Paramyxine yangi), gray mullet (Mugil cephalus), zebrafish (Danio rerio), toad (Bufo melanosticus), and snake (Boiga kraepelini). Only one PGI gene was cloned in each of hagfish, toad, and snake, but two PGI genes were found in zebrafish and gray mullet, respectively. The PGI of hagfish encodes 554 amino acids, in contrast to the PGIs of bonyfishes, toad, and snake which encode 553 amino acids and the PGIs of mammals which encode 558 amino acids. Among 558 aligned amino acid sites, there are 314 sites (56.27%) totally conserved. To see if diversifying selection acts on PGI amino acids of vertebrates, we calculated the pairwise ratio of nonsynonymous versus synonymous substitution per site (Ka/Ks) and the ratio of radical amino acid changes versus conservative amino acid changes per sites (dR/dC) between PGI sequences. The average pairwise ratio between nonsynonymous substitutions per nucleotide (Ka) and synonymous substitutions per nucleotide (Ks) among vertebrate PGI sequences equals 0.047 +/- 0.019. The average pairwise ratio between radical amino acid changes and conservative amino acid changes (dR/dC) among the vertebrate PGIs equal 0.938 +/- 0.158 for charge changes, 0.558 +/- 0.085 for polarity changes, and 0.465 +/- 0.0714 when both polarity and volume are considered. There is no amino acid within the vertebrate PGIs under diversifying selection as analyzed by the method of Yang et al. (2000b). The results suggest that the present vertebrate PGIs are at evolutionary stasis and are being subjected to intense purifying selection. The purifying selection is to maintain polarity and volume of the protein but not the charge groups of amino acids. Phylogenetic analysis reveals that vertebrate PGIs can be classified into three major groups: the mammalian, amphibian-reptilian, and teleostean PGIs. The gene tree suggests that the gene duplication event of PGI in bonyfishes occurred before diversification of Acanthopterygii but after the split of bonyfishes and tetrapods. The evolution of multiple functions of PGI is discussed.     

Phylogenetic Gaussian Process Model for the Inference of Functionally Important Regions in Protein Tertiary Structures

A critical question in biology is the identification of functionally important amino acid sites in proteins. Because functionally important sites are under stronger purifying selection, site-specific substitution rates tend to be lower than usual at these sites. A large number of phylogenetic models have been developed to estimate site-specific substitution rates in proteins and the extraordinarily low substitution rates have been used as evidence of function. Most of the existing tools, e.g. Rate4Site, assume that site-specific substitution rates are independent across sites. However, site-specific substitution rates may be strongly correlated in the protein tertiary structure, since functionally important sites tend to be clustered together to form functional patches. We have developed a new model, GP4Rate, which incorporates the Gaussian process model with the standard phylogenetic model to identify slowly evolved regions in protein tertiary structures. GP4Rate uses the Gaussian process to define a nonparametric prior distribution of site-specific substitution rates, which naturally captures the spatial correlation of substitution rates. Simulations suggest that GP4Rate can potentially estimate site-specific substitution rates with a much higher accuracy than Rate4Site and tends to report slowly evolved regions rather than individual sites. In addition, GP4Rate can estimate the strength of the spatial correlation of substitution rates from the data. By applying GP4Rate to a set of mammalian B7-1 genes, we found a highly conserved region which coincides with experimental evidence. GP4Rate may be a useful tool for the in silico prediction of functionally important regions in the proteins with known structures.     

The molecular evolution of acrosin in placental mammals

Acrosin is thought to fulfill several different roles in fertilization including that of a serine protease and in secondary zona pellucida (ZP) binding. However, acrosin's importance as a fertilization protein has been questioned. Especially since it was discovered that acrosin knockout mice are fertile. In this study, we explored the sites involved in serine protease activity and secondary binding. We also assessed conservation in functional sites across species and examined whether amino acid changes present in the human population have the potential to affect fertility. In addition, since many mammalian reproduction proteins have been found to evolve rapidly, we tested for positive selection. Sequences from 43 mammals from all 19 placental orders, which included a total of 828 nucleotides from acrosin exons 2, 3, 4, and a portion of exon 5, were obtained. We found that all sites of the serine catalytic triad as well as three other sites linked to catalytic activity were completely conserved. Five of six sites proposed to play a role in secondary binding were 100% conserved as basic residues. These results support an evolutionary conserved role for acrosin as a serine protease and secondary binding protein across placental mammals. We found statistically significant support for positive selection within acrosin, but no single amino acid site reached the significance level of P > 0.95 for inclusion within the category omega > 1. Based upon two amino acid mutation scoring systems, three out of seven human residue changing single nucleotide polymorphisms (SNPs) were found to be potentially protein-altering mutations.     

竹亚科lea3基因的进化分析

选取竹亚科中两个超族、六个族和三个亚族的10个竹种为材料,分别是泰竹、凤尾竹、青皮竹、大叶慈、慈竹、野龙竹、毛竹、香竹、苦竹、菲白竹,分离克隆了它们的lea3基因,并将它们与外类群物种水稻进行序列比对和进化分析。结果发现在分支模型与分支位点模型的检测中,不同竹种所含lea3基因承受了不同的正选择压力,清除选择作用在lea3基因编码区中占主导地位（ω<1）。在位点模型的检测中,共检测出了18个显著性正选择位点,占总氨基酸数目的111%。对这18个显著性正选择位点进行定位后,发现其中的15个位于11个氨基酸串联重复序列附近。这说明lea3基因中的11个氨基酸串联重复序列区比基因其它区域更容易受自然选择作用影响。同时,在位点模型检测结果的基础上,通过对强烈清除选择位点的定位,发现在11个氨基酸串联重复序列区内存在一长段无强烈清除位点的序列区。     

Molecular evolution of <Emphasis Type="Italic">PKD2</Emphasis> gene family in mammals

PKD2 gene encodes a critical cation channel protein that plays important roles in various developmental processes and is usually evolutionarily conserved. In the present study, we analyzed the evolutionary patterns of PKD2 and its homologous genes (PKD2L1, PKD2L2) from nine mammalian species. In this study, we demonstrated the orthologs of PKD2 gene family evolved under a dominant purifying selection force. Our results in combination with the reported evidences from functional researches suggested the entire PKD2 gene family are conserved and perform essential biological roles during mammalian evolution. In rodents, PKD2 gene family members appeared to have evolved more rapidly than other mammalian lineages, probably resulting from relaxation of purifying selection. However, positive selection imposed on synonymous sites also potentially contributed to this case. For the paralogs, our results implied that PKD2L2 genes evolved under a weaker purifying selection constraint than PKD2 and PKD2L1 genes. Interestingly, some loop regions of transmembrane domain of PKD2L2 exhibited higher P _N/P _S ratios than expected, suggesting these regions are more functional divergent in organisms and worthy of special attention. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users. Chun Ye, Huan Sun have contributed equally to this work.     

ACP5 (Uteroferrin): phylogeny of an ancient and conserved gene expressed in the endometrium of mammals

Type 5 acid phosphatase (ACP5; also known as tartrate-resistant acid phosphatase or uteroferrin) is a metalloprotein secreted by the endometrial glandular epithelium of pigs, mares, sheep, and water buffalo. In this paper, we describe the phylogenetic distribution of endometrial expression of ACP5 and demonstrate that endometrial expression arose early in evolution (i.e., before divergence of prototherian and therian mammals ~166 million years ago). To determine expression of ACP5 in the pregnant endometrium, RNA was isolated from rhesus, mouse, rat, dog, sheep, cow, horse, armadillo, opossum, and duck-billed platypus. Results from RT-PCR and RNA-Seq experiments confirmed that ACP5 is expressed in all species examined. ACP5 was also demonstrated immunochemically in endometrium of rhesus, marmoset, sheep, cow, goat, and opossum. Alignment of inferred amino acid sequences shows a high conservation of ACP5 throughout speciation, with species-specific differences most extensive in the N-terminal and C-terminal regions of the protein. Analysis by Selecton indicated that most of the sites in ACP5 are undergoing purifying selection, and no sites undergoing positive selection were found. In conclusion, endometrial expression of ACP5 is a common feature in all orders of mammals and has been subjected to purifying selection. Expression of ACP5 in the uterus predates the divergence of therians and prototherians. ACP5 is an evolutionary conserved gene that likely exerts a common function important for pregnancy in mammals using a wide range of reproductive strategies.     

Purifying selection maintains highly conserved noncoding sequences in Drosophila

The majority of metazoan genomes consist of nonprotein-coding regions, although the functional significance of most noncoding DNA sequences remains unknown. Highly conserved noncoding sequences (CNSs) have proven to be reliable indicators of functionally constrained sequences such as cis-regulatory elements and noncoding RNA genes. However, CNSs may arise from nonselective evolutionary processes such as genomic regions with extremely low mutation rates known as mutation "cold spots." Here we combine comparative genomic data from recently completed insect genome projects with population genetic data in Drosophila melanogaster to test predictions of the mutational cold spot model of CNS evolution in the genus Drosophila. We find that point mutations in intronic and intergenic CNSs exhibit a significant reduction in levels of divergence relative to levels of polymorphism, as well as a significant excess of rare derived alleles, compared with either the nonconserved spacer regions between CNSs or with 4-fold silent sites in coding regions. Controlling for the effects of purifying selection, we find no evidence of positive selection acting on Drosophila CNSs, although we do find evidence for the action of recurrent positive selection in the spacer regions between CNSs. We estimate that approximately 85% of sites in Drosophila CNSs are under constraint with selection coefficients (N(e)s) on the order of 10-100, and thus, the estimated strength and number of sites under purifying selection is greater for Drosophila CNSs relative to those in the human genome. These patterns of nonneutral molecular evolution are incompatible with the mutational cold spot hypothesis to explain the existence of CNSs in Drosophila and, coupled with similar findings in mammals, argue against the general likelihood that CNSs are generated by mutational cold spots in any metazoan genome.     

Frequent and widespread parallel evolution of protein sequences

Understanding the patterns and causes of protein sequence evolution is a major challenge in evolutionary biology. One of the critical unresolved issues is the relative contribution of selection and genetic drift to the fixation of amino acid sequence differences between species. Molecular homoplasy, the independent evolution of the same amino acids at orthologous sites in different taxa, is one potential signature of selection; however, relatively little is known about its prevalence in eukaryotic proteomes. To quantify the extent and type of homoplasy among evolving proteins, we used phylogenetic methodology to analyze 8 genome-scale data matrices from clades of different evolutionary depths that span the eukaryotic tree of life. We found that the frequency of homoplastic amino acid substitutions in eukaryotic proteins was more than 2-fold higher than expected under neutral models of protein evolution. The overwhelming majority of homoplastic substitutions were parallelisms that involved the most frequently exchanged amino acids with similar physicochemical properties and that could be reached by a single-mutational step. We conclude that the role of homoplasy in shaping the protein record is much larger than generally assumed, and we suggest that its high frequency can be explained by both weak positive selection for certain substitutions and purifying selection that constrains substitutions to a small number of functionally equivalent amino acids.     

Widespread positive selection in synonymous sites of mammalian genes

Evolution of protein sequences is largely governed by purifying selection, with a small fraction of proteins evolving under positive selection. The evolution at synonymous positions in protein-coding genes is not nearly as well understood, with the extent and types of selection remaining, largely, unclear. A statistical test to identify purifying and positive selection at synonymous sites in protein-coding genes was developed. The method compares the rate of evolution at synonymous sites (Ks) to that in intron sequences of the same gene after sampling the aligned intron sequences to mimic the statistical properties of coding sequences. We detected purifying selection at synonymous sites in approximately 28% of the 1,562 analyzed orthologous genes from mouse and rat, and positive selection in approximately 12% of the genes. Thus, the fraction of genes with readily detectable positive selection at synonymous sites is much greater than the fraction of genes with comparable positive selection at nonsynonymous sites, i.e., at the level of the protein sequence. Unlike other genes, the genes with positive selection at synonymous sites showed no correlation between Ks and the rate of evolution in nonsynonymous sites (Ka), indicating that evolution of synonymous sites under positive selection is decoupled from protein evolution. The genes with purifying selection at synonymous sites showed significant anticorrelation between Ks and expression level and breadth, indicating that highly expressed genes evolve slowly. The genes with positive selection at synonymous sites showed the opposite trend, i.e., highly expressed genes had, on average, higher Ks. For the genes with positive selection at synonymous sites, a significantly lower mRNA stability is predicted compared to the genes with negative selection. Thus, mRNA destabilization could be an important factor driving positive selection in nonsynonymous sites, probably, through regulation of expression at the level of mRNA degradation and, possibly, also translation rate. So, unexpectedly, we found that positive selection at synonymous sites of mammalian genes is substantially more common than positive selection at the level of protein sequences. Positive selection at synonymous sites might act through mRNA destabilization affecting mRNA levels and translation.     

Unusual mammalian usage of TGA stop codons reveals that sequence conservation need not imply purifying selection

The assumption that conservation of sequence implies the action of purifying selection is central to diverse methodologies to infer functional importance. GC-biased gene conversion (gBGC), a meiotic mismatch repair bias strongly favouring GC over AT, can in principle mimic the action of selection, this being thought to be especially important in mammals. As mutation is GC→AT biased, to demonstrate that gBGC does indeed cause false signals requires evidence that an AT-rich residue is selectively optimal compared to its more GC-rich allele, while showing also that the GC-rich alternative is conserved. We propose that mammalian stop codon evolution provides a robust test case. Although in most taxa TAA is the optimal stop codon, TGA is both abundant and conserved in mammalian genomes. We show that this mammalian exceptionalism is well explained by gBGC mimicking purifying selection and that TAA is the selectively optimal codon. Supportive of gBGC, we observe (i) TGA usage trends are consistent at the focal stop codon and elsewhere (in UTR sequences); (ii) that higher TGA usage and higher TAA→TGA substitution rates are predicted by a high recombination rate; and (iii) across species the difference in TAA <-> TGA substitution rates between GC-rich and GC-poor genes is largest in genomes that possess higher between-gene GC variation. TAA optimality is supported both by enrichment in highly expressed genes and trends associated with effective population size. High TGA usage and high TAA→TGA rates in mammals are thus consistent with gBGC’s predicted ability to “drive” deleterious mutations and supports the hypothesis that sequence conservation need not be indicative of purifying selection. A general trend for GC-rich trinucleotides to reside at frequencies far above their mutational equilibrium in high recombining domains supports the generality of these results.

Is sequence conservation a sign of purifying selection and hence functional importance? This analysis of why mammals use and conserve the most error-prone stop codon suggests not, consistent with GC-biased gene conversion’s predicted ability to “drive” deleterious mutations and supporting the hypothesis that sequence conservation need not be indicative of purifying selection.     

ApoE: The role of conserved residues in defining function

The amino acid sequences of apolipoprotein E (apoE) from 63 different mammalian species have been downloaded from the protein database. The sequences were compared to human apoE4 to determine conserved and non‐conserved sequences of amino acids. ApoE4 is the major risk factor for the development of late onset Alzheimer's disease while apoE3, which differs from apoE4 by a single amino acid change at position 112, poses little or no risk for the development of this disease. Thus, the two proteins appear to be structurally and functionally different. Seven highly conserved regions, representing approximately 47 amino acids (of 299) have been found. These regions are distributed throughout the protein and reflect ligand binding sites as well as regions proposed to be involved in the propagation of the cysteine–arginine change at position 112 to distant regions of the protein in the N‐ and C‐terminal domains. Highly non‐conserved regions are at the N‐ and C‐terminal ends of the apoE protein.     

Molecular evolution of the mammalian prion protein

Prion protein (PrP) sequences are until now available for only six of the 18 orders of placental mammals. A broader comparison of mammalian prions might help to understand the enigmatic functional and pathogenic properties of this protein. We therefore determined PrP coding sequences in 26 mammalian species to include all placental orders and major subordinal groups. Glycosylation sites, cysteines forming a disulfide bridge, and a hydrophobic transmembrane region are perfectly conserved. Also, the sequences responsible for secondary structure elements, for N- and C-terminal processing of the precursor protein, and for attachment of the glycosyl-phosphatidylinositol membrane anchor are well conserved. The N-terminal region of PrP generally contains five or six repeats of the sequence P(Q/H)GGG(G/-)WGQ, but alleles with two, four, and seven repeats were observed in some species. This suggests, together with the pattern of amino acid replacements in these repeats, the regular occurrence of repeat expansion and contraction. Histidines implicated in copper ion binding and a proline involved in 4-hydroxylation are lacking in some species, which questions their importance for normal functioning of cellular PrP. The finding in certain species of two or seven repeats, and of amino acid substitutions that have been related to human prion diseases, challenges the relevance of such mutations for prion pathology. The gene tree deduced from the PrP sequences largely agrees with the species tree, indicating that no major deviations occurred in the evolution of the prion gene in different placental lineages. In one species, the anteater, a prion pseudogene was present in addition to the active gene.     

Contrasting modes of natural selection acting on pigmentation genes in the Drosophila dunni subgroup

Genes that encode for divergent adaptive traits may have genealogies that contrast with those from loci that are not functionally involved in differentiation. Here, we examine DNA sequence variation among the species of the eastern Caribbean Drosophila dunni subgroup at two loci, yellow and dopa decaboxylase (Ddc), which both play integral roles in pigmentation patterning of adult Drosophila. Phylogenetic analyses of these loci produce gene genealogies with topologies that mirror those described for other nuclear genes: the six morphologically distinct species within the subgroup are divided into only three lineages, with one lineage containing four species that share extensive ancestral polymorphism. At the Ddc locus these major lineages are delineated only by silent site variation. We observe a significantly higher rate of synonymous site divergence than non-synonymous divergence, consistent with strong purifying selection acting on the locus. In contrast, the yellow locus exhibits patterns of amino acid divergence and nucleotide diversity that are consistent with recent diversifying selection acting in two different lineages. This selection appears to be targeting amino acid variants in the signal sequence of the Yellow protein, a region which is tightly constrained among members of the larger D. cardini radiation. This result highlights not only the potential importance of yellow in the evolution of divergent pigmentation patterns among members of the D. dunni subgroup, but also hints that variation in signal peptide sequences may play a role in phenotypic diversification.     

Assessing the fraction of short-distance tandem splice sites under purifying selection

Many alternative splice events result in subtle mRNA changes, and most of them occur at short-distance tandem donor and acceptor sites. The splicing mechanism of such tandem sites likely involves the stochastic selection of either splice site. While tandem splice events are frequent, it is unknown how many are functionally important. Here, we use phylogenetic conservation to address this question, focusing on tandems with a distance of 3-9 nucleotides. We show that previous contradicting results on whether alternative or constitutive tandem motifs are more conserved between species can be explained by a statistical paradox (Simpson's paradox). Applying methods that take biases into account, we found higher conservation of alternative tandems in mouse, dog, and even chicken, zebrafish, and Fugu genomes. We estimated a lower bound for the number of alternative sites that are under purifying (negative) selection. While the absolute number of conserved tandem motifs decreases with the evolutionary distance, the fraction under selection increases. Interestingly, a number of frameshifting tandems are under selection, suggesting a role in regulating mRNA and protein levels via nonsense-mediated decay (NMD). An analysis of the intronic flanks shows that purifying selection also acts on the intronic sequence. We propose that stochastic splice site selection can be an advantageous mechanism that allows constant splice variant ratios in situations where a deviation in this ratio is deleterious.     

The relation between the neutrality index for mitochondrial genes and the distribution of mutational effects on fitness

We explore factors affecting patterns of polymorphism and divergence (as captured by the neutrality index) at mammalian mitochondrial loci. To do this, we develop a population genetic model that incorporates a fraction of neutral amino acid sites, mutational bias, and a probability distribution of selection coefficients against new nonsynonymous mutations. We confirm, by reanalyzing publicly available datasets, that the mitochondrial cyt-b gene shows a broad range of neutrality indices across mammalian taxa, and explore the biological factors that can explain this observation. We find that observed patterns of differences in the neutrality index, polymorphism, and divergence are not caused by differences in mutational bias. They can, however, be explained by a combination of a small fraction of neutral amino acid sites, weak selection acting on most amino acid mutations, and differences in effective population size among taxa.