Opposite evolutionary effects between different alternative splicing patterns

Alternative splicing (AS) has been recognized as a mechanism of relaxing selection pressure on protein subsequences. Here, we show that AS may also yield contrary evolutionary effects. We compare the evolutionary rates of 2 types of alternatively spliced exons (ASEs)-simple and complex. The former does not change the boundaries of its flanking exons, whereas the latter does. By analyzing over 26,000 human-mouse orthologous exons, we demonstrate that complex ASEs have lower Ka and Ka/Ks ratio and higher Ks than constitutively spliced exons (CSEs), whereas simple ASEs have evolutionary rates to the opposite of CSEs. Our results indicate that complex ASEs are subject to stronger selection pressure than CSEs at the protein level, but the trend is reversed at the RNA level. Therefore, the previous view that ASEs accelerate evolution of protein subsequences needs to be modified.     

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.     

Evidence for purifying selection against synonymous mutations in mammalian exonic splicing enhancers

Silent sites in mammals have classically been assumed to be free from selective pressures. Consequently, the synonymous substitution rate (Ks) is often used as a proxy for the mutation rate. Although accumulating evidence demonstrates that the assumption is not valid, the mechanism by which selection acts remain unclear. Recent work has revealed that the presence of exonic splicing enhancers (ESEs) in coding sequence might influence synonymous evolution. ESEs are predominantly located near intron-exon junctions, which may explain the reduced single-nucleotide polymorphism (SNP) density in these regions. Here we show that synonymous sites in putative ESEs evolve more slowly than the remaining exonic sequence. Differential mutabilities of ESEs do not appear to explain this difference. We observe that substitution frequency at fourfold synonymous sites decreases as one approaches the ends of exons, consistent with the existing SNP data. This gradient is at least in part explained by ESEs being more abundant near junctions. Between-gene variation in Ks is hence partly explained by the proportion of the gene that acts as an ESE. Given the relative abundance of ESEs and the reduced rates of synonymous divergence within them, we estimate that constraints on synonymous evolution within ESEs causes the true mutation rate to be underestimated by not more than approximately 8%. We also find that Ks outside of ESEs is much lower in alternatively spliced exons than in constitutive exons, implying that other causes of selection on synonymous mutations exist. Additionally, selection on ESEs appears to affect nonsynonymous sites and may explain why amino acid usage near intron-exon junctions is nonrandom.     

Different alternative splicing patterns are subject to opposite selection pressure for protein reading frame preservation

Background  

Alternative splicing (AS) has been regarded capable of altering selection pressure on protein subsequences. Particularly, the frequency of reading frame preservation (FRFP), as a measure of selection pressure, has been reported to be higher in alternatively spliced exons (ASEs) than in constitutively spliced exons (CSEs). However, recently it has been reported that different ASE types – simple and complex ASEs – may be subject to opposite selection forces. Therefore, it is necessary to re-evaluate the evolutionary effects of such splicing patterns on frame preservation.     

Positive selection in alternatively spliced exons of human genes

Alternative splicing is a well-recognized mechanism of accelerated genome evolution. We have studied single-nucleotide polymorphisms and human-chimpanzee divergence in the exons of 6672 alternatively spliced human genes, with the aim of understanding the forces driving the evolution of alternatively spliced sequences. Here, we show that alternatively spliced exons and exon fragments (alternative exons) from minor isoforms experience lower selective pressure at the amino acid level, accompanied by selection against synonymous sequence variation. The results of the McDonald-Kreitman test suggest that alternatively spliced exons, unlike exons constitutively included in the mRNA, are also subject to positive selection, with up to 27% of amino acids fixed by positive selection.     

Nonneutral evolution and differential mutation rate of gender-associated mitochondrial DNA lineages in the marine mussel Mytilus.

Mussels have two types of mitochondrial DNA (mtDNA). The M type is transmitted paternally, and the F type is transmitted maternally. To test hypotheses of the molecular evolution of both mtDNA genomes, 50 nucleotide sequences were obtained for 396 bp of the COIII gene of European populations of Mytilus edulis and the Atlantic and Mediterranean forms of M. galloprovincialis. Analysis based on the proportion of synonymous and nonsynonymous substitutions indicate that mtDNA is evolving in a non-neutral and complex fashion. Previous studies on American mussels demonstrated that the F genome experiences a higher purifying selection and that the M genome evolves faster. Here we show that these patterns also hold in European populations. However, in contrast to American populations, where an excess of replacement substitution between F and M lineages has been reported, a significant excess of replacement polymorphism within mtDNA lineages is observed in European populations of M. galloprovincialis. European populations also show an excess of replacement polymorphism within the F but not within the M genome with respect to American M. trossulus, as well as a consistent pattern of excess of rare variants in both F and M genomes. These results are consistent with a nearly neutral model of molecular evolution and a recent relaxation of selective constraints on European mtDNA. Levels of diversity are significantly higher for the M than F genome, and the M genome also accumulates synonymous and nonsynonymous substitutions at a higher rate, in contrast with earlier reports where no difference for the synonymous rate was observed. It is suggested that a subtle balance between relaxed selection and a higher mutation rate explains the faster evolutionary rate of the M lineage.     

Molecular evolution in the gnd locus of Salmonella enterica

The gnd gene, the structural gene for 6-phosphogluconate dehydrogenase, was sequenced and analyzed in 34 isolates from different serovars of the seven subspecies of Salmonella enterica to provide comparative information on the evolution in this gene, which has been studied extensively in Escherichia coli. The gene tree obtained by the neighbor- joining method in general gave separate branches for each subspecies, with the few exceptions readily explained by recombination. There is evidence of recombination involving transfer of long (more than 400 bp) and short (30-150 bp) segments of DNA. Four of the six long-segment transfers detected are at the 5' end of the gene, and in all four cases a variant of the chi sequence is located close to the recombination junction and appears to have mediated the recombination events. We suggest that in these four cases and in a fifth case with intersubspecies transfer of the whole gnd gene, the adjacent rfb (O antigen) locus may have been transferred in the same event. The estimates of the number of synonymous substitutions per synonymous site, KS, and the number of nonsynonymous substitutions per nonsynonymous site, KA, within the E. coli and S. enterica gnd genes, and also between the two species show an interesting distribution, with KS being lower toward the ends of the gene and KA in particular being lower in the first than in the second domain. In S. enterica, synonymous sites also seem to be subjected to negative selection. The ratio of KA to KS was higher within S. enterica and E. coli than between them, which may indicate that intraspecies variation is essentially between clones and that mildly deleterious mutations can be fixed within clones, which would thus raise KA within species.      

Codon Usage Bias and Base Composition of Nuclear Genes in Drosophila

The nuclear genes of Drosophila evolve at various rates. This variation seems to correlate with codon-usage bias. In order to elucidate the determining factors of the various evolutionary rates and codon-usage bias in the Drosophila nuclear genome, we compared patterns of codon-usage bias with base compositions of exons and introns. Our results clearly show the existence of selective constraints at the translational level for synonymous (silent) sites and, on the other hand, the neutrality or near neutrality of long stretches of nucleotide sequence within noncoding regions. These features were found for comparisons among nuclear genes in a particular species (Drosophila melanogaster, Drosophila pseudoobscura and Drosophila virilis) as well as in a particular gene (alcohol dehydrogenase) among different species in the genus Drosophila. The patterns of evolution of synonymous sites in Drosophila are more similar to those in the prokaryotes than they are to those in mammals. If a difference in the level of expression of each gene is a main reason for the difference in the degree of selective constraint, the evolution of synonymous sites of Drosophila genes would be sensitive to the level of expression among genes and would change as the level of expression becomes altered in different species. Our analysis verifies these predictions and also identifies additional selective constraints at the translational level in Drosophila.     

Sex-linked mammalian sperm proteins evolve faster than autosomal ones

X-linked genes can evolve slower or faster depending on whether most recessive, or at least partially recessive alleles are deleterious or beneficial due to their hemizygous expression in males. Molecular studies of X chromosome divergence have provided conflicting evidence for both a higher and lower rate of nucleotide substitution at both synonymous and nonsynonymous sites, depending on the nucleotide sites sampled. Using human and mouse orthologous genes, we tested the hypothesis that genes encoding male-specific sperm proteins are evolving faster on the X chromosome compared with autosomes. X-linked sperm proteins have an average nonsynonymous mutation rate almost twice as high as sperm genes found on autosomes, unlike other tissue-specific genes, where no significant difference in the nonsynonymous mutation rate between the X chromosome and autosomes was found. However, no difference was found in the average synonymous mutation rate of X-linked versus autosomal sperm proteins, which along with corresponding higher values of Ka/Ks in X-linked sperm proteins suggest that differences in selective forces and not mutation rates are the underlying cause of higher X-linked mammalian sperm protein divergence.     

Rapid evolution of a primate sperm protein: relaxation of functional constraint or positive Darwinian selection?

Protamines are arginine-rich proteins that replace histones and bind sperm DNA during spermatogenesis in vertebrates. Previous studies have shown that protamine exons evolve faster than does the protamine intron. It has been suggested that this is a result of a relaxation of functional constraint. However, a more likely explanation is that the evolutionary rate of exons has been accelerated by positive Darwinian selection, because introns are generally believed to evolve in a neutral fashion. Therefore, we examined the possibility that positive selection has been acting on the protamine genes of three groups of placental mammals: primates (hominoids and Old World monkeys), rodents (mice, rats, and guinea pigs), and pecoran ruminants (deer and bovids). We found that the nucleotide substitution rate at nonsynonymous sites is significantly higher than the rate at synonymous and intron sites for protamine P1 of hominoids and Old World monkeys. This result suggests that positive selection has been operating on protamine P1 of these species. In contrast, no clear-cut evidence of positive selection was found for protamine P1 of ruminants and rodents or protamine P2 of primates. The agent of positive selection on primate protamine P1 remains unknown, though sperm competition is a possibility. Further investigations on the function and intraspecific polymorphism of this protein are needed in order to identify the selection agent.     

A test of translational selection at 'silent' sites in the human genome: base composition comparisons in alternatively spliced genes

Natural selection appears to discriminate among synonymous codons to enhance translational efficiency in a wide range of prokaryotes and eukaryotes. Codon bias is strongly related to gene expression levels in these species. In addition, between-gene variation in silent DNA divergence is inversely correlated with codon bias. However, in mammals, between-gene comparisons are complicated by distinctive nucleotide-content bias (isochores) throughout the genome. In this study, we attempted to identify translational selection by analyzing the DNA sequences of alternatively spliced genes in humans and in Drosophila melanogaster. Among codons in an alternatively spliced gene, those in constitutively expressed exons are translated more often than those in alternatively spliced exons. Thus, translational selection should act more strongly to bias codon usage and reduce silent divergence in constitutive than in alternative exons. By controlling for regional forces affecting base-composition evolution, this within-gene comparison makes it possible to detect codon selection at synonymous sites in mammals. We found that GC-ending codons are more abundant in constitutive than alternatively spliced exons in both Drosophila and humans. Contrary to our expectation, however, silent DNA divergence between mammalian species is higher in constitutive than in alternative exons.     

Selection in favor of nucleotides G and C diversifies evolution rates and levels of polymorphism at mammalian synonymous sites

The impact of synonymous nucleotide substitutions on fitness in mammals remains controversial. Despite some indications of selective constraint, synonymous sites are often assumed to be neutral, and the rate of their evolution is used as a proxy for mutation rate. We subdivide all sites into four classes in terms of the mutable CpG context, nonCpG, postC, preG, and postCpreG, and compare four-fold synonymous sites and intron sites residing outside transposable elements. The distribution of the rate of evolution across all synonymous sites is trimodal. Rate of evolution at nonCpG synonymous sites, not preceded by C and not followed by G, is approximately 10% below that at such intron sites. In contrast, rate of evolution at postCpreG synonymous sites is approximately 30% above that at such intron sites. Finally, synonymous and intron postC and preG sites evolve at similar rates. The relationship between the levels of polymorphism at the corresponding synonymous and intron sites is very similar to that between their rates of evolution. Within every class, synonymous sites are occupied by G or C much more often than intron sites, whose nucleotide composition is consistent with neutral mutation-drift equilibrium. These patterns suggest that synonymous sites are under weak selection in favor of G and C, with the average coefficient s approximately 0.25/Ne approximately 10(-5), where Ne is the effective population size. Such selection decelerates evolution and reduces variability at sites with symmetric mutation, but has the opposite effects at sites where the favored nucleotides are more mutable. The amino-acid composition of proteins dictates that many synonymous sites are CpGprone, which causes them, on average, to evolve faster and to be more polymorphic than intron sites. An average genotype carries approximately 10(7) suboptimal nucleotides at synonymous sites, implying synergistic epistasis in selection against them.     

Similarity of Synonymous Substitution Rates Across Mammalian Genomes

Given that a gene has a high (or low) synonymous substitution rate in one mammalian species, will it also have a high (or low) synonymous substitution rate in another mammalian species? Such similarities in the rate of synonymous substitution can reveal both selective pressures and neutral processes acting on mammalian gene sequences; however, the existence of such an effect has been a matter of disagreement. We resolve whether such synonymous substitution rate similarities exist using 7462 ortholog triplets aligned across rat, mouse, and human, a dataset two orders of magnitude larger than previous studies. We find that a gene’s synonymous substitution rate in the rat-mouse branch of the phylogeny is correlated with its rate in the branch connecting human and the rat-mouse ancestor. We confirm this for several different measures of synonymous substitution rate, including corrections for base composition and CpG dinucleotides, and we verify the results in the larger mouse-human-rat-dog phylogeny. This similarity of rates is most apparent for genes in which synonymous sites are well conserved across species, suggesting that a significant component of the effect is due to purifying selection. We observe rate correlations at a resolution as fine as a few hundred kilobases, and the genes with the most similar synonymous substitution rates are enriched for regulatory functions. Genes with above-average substitution rates also exhibit significant, though somewhat weaker, rate correlations, suggesting that some neutral processes may have persisted in the phylogeny as well.     

Nucleotide Substitution Rate of Mammalian Mitochondrial Genomes

We present here for the first time a comprehensive study based on the analysis of closely related organisms to provide an accurate determination of the nucleotide substitution rate in mammalian mitochondrial genomes. This study examines the evolutionary pattern of the different functional mtDNA regions as accurately as possible on the grounds of available data, revealing some important ``genomic laws.' The main conclusions can be summarized as follows. (1) High intragenomic variability in the evolutionary dynamic of mtDNA was found. The substitution rate is strongly dependent on the region considered, and slow- and fast-evolving regions can be identified. Nonsynonymous sites, the D-loop central domain, and tRNA and rRNA genes evolve much more slowly than synonymous sites and the two peripheral D-loop region domains. The synonymous rate is fairly uniform over the genome, whereas the rate of nonsynonymous sites depends on functional constraints and therefore differs considerably between genes. (2) The commonly accepted statement that mtDNA evolves more rapidly than nuclear DNA is valid only for some regions, thus it should be referred to specific mitochondrial components. In particular, nonsynonymous sites show comparable rates in mitochondrial and nuclear genes; synonymous sites and small rRNA evolve about 20 times more rapidly and tRNAs about 100 times more rapidly in mitochondria than in their nuclear counterpart. (3) A species-specific evolution is particularly evident in the D-loop region. As the divergence times of the organism pairs under consideration are known with sufficient accuracy, absolute nucleotide substitution rates are also provided. Received: 11 May 1998 / Accepted: 2 September 1998     

A comprehensive survey of human polymorphisms at conserved splice dinucleotides and its evolutionary relationship with alternative splicing

Background  

Alternative splicing (AS) is a key molecular process that endows biological functions with diversity and complexity. Generally, functional redundancy leads to the generation of new functions through relaxation of selective pressure in evolution, as exemplified by duplicated genes. It is also known that alternatively spliced exons (ASEs) are subject to relaxed selective pressure. Within consensus sequences at the splice junctions, the most conserved sites are dinucleotides at both ends of introns (splice dinucleotides). However, a small number of single nucleotide polymorphisms (SNPs) occur at splice dinucleotides. An intriguing question relating to the evolution of AS diversity is whether mutations at splice dinucleotides are maintained as polymorphisms and produce diversity in splice patterns within the human population. We therefore surveyed validated SNPs in the database dbSNP located at splice dinucleotides of all human genes that are defined by the H-Invitational Database.     

The problem of counting sites in the estimation of the synonymous and nonsynonymous substitution rates: implications for the correlation between the synonymous substitution rate and codon usage bias

Most methods for estimating the rate of synonymous and nonsynonymous substitution per site define a site as a mutational opportunity: the proportion of sites that are synonymous is equal to the proportion of mutations that would be synonymous under the model of evolution being considered. Here we demonstrate that this definition of a site can give misleading results and that a physical definition of site should be used in some circumstances. We illustrate our point by reexamining the relationship between codon usage bias and the synonymous substitution rate. It has recently been shown that the rate of synonymous substitution, calculated using the Goldman-Yang method, which encapsulates the mutational-opportunity definition of a site at a high level of sophistication, is either positively correlated or uncorrelated to synonymous codon bias in Drosophila. Using other methods, which account for synonymous codon bias but define a site physically, we show that there is a negative correlation between the synonymous substitution rate and codon bias and that the lack of a negative correlation using the Goldman-Yang method is due to the way in which the number of synonymous sites is counted. We also show that there is a positive correlation between the synonymous substitution rate and third position GC content in mammals, but that the relationship is considerably weaker than that obtained using the Goldman-Yang method. We argue that the Goldman-Yang method is misleading in this context and conclude that methods that rely on a mutational-opportunity definition of a site should be used with caution.     

Variation in synonymous codon use and DNA polymorphism within the Drosophila genome

A strong negative correlation between the rate of amino-acid substitution and codon usage bias in Drosophila has been attributed to interference between positive selection at nonsynonymous sites and weak selection on codon usage. To further explore this possibility we have investigated polymorphism and divergence at three kinds of sites: synonymous, nonsynonymous and intronic in relation to codon bias in D. melanogaster and D. simulans. We confirmed that protein evolution is one of the main explicative parameters for interlocus codon bias variation (r(2) approximately 40%). However, intron or synonymous diversities, which could have been expected to be good indicators of local interference [here defined as the additional increase of drift due to selection on tightly linked sites, also called 'genetic draft' by Gillespie (2000)] did not covary significantly with codon bias or with protein evolution. Concurrently, levels of polymorphism were reduced in regions of low recombination rates whereas codon bias was not. Finally, while nonsynonymous diversities were very well correlated between species, neither synonymous nor intron diversities observed in D. melanogaster were correlated with those observed in D. simulans. All together, our results suggest that the selective constraint on the protein is a stable component of gene evolution while local interference is not. The pattern of variation in genetic draft along the genome therefore seems to be instable through evolutionary times and should therefore be considered as a minor determinant of codon bias variance. We argue that selective constraints for optimal codon usage are likely to be correlated with selective constraints on the protein, both between codons within a gene, as previously suggested, and also between genes within a genome.