Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “synonymous substitution”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 271 records · Page 15Linked to original sources

Pattern of nucleotide substitution at major histocompatibility complex class I loci reveals overdominant selection.

The major histocompatibility complex (MHC) loci are known to be highly polymorphic in humans, mice and certain other mammals, with heterozygosity as high as 80-90% (ref. 1). Four different hypotheses have been proposed to explain this high degree of polymorphism: (1) a high mutation rate, (2) gene conversion or interlocus genetic exchange, (3) over dominant selection and (4) frequency-dependent selection. In an attempt to establish which of these hypotheses is correct, we examined the pattern of nucleotide substitution between polymorphic alleles in the region of the antigen recognition site (ARS) and other regions of human and mouse class I MHC genes. The results indicate that in ARS the rate of nonsynonymous (amino acid altering) substitution is significantly higher than that of synonymous substitution in both humans and mice, whereas in other regions the reverse is true. This observation, together with a theoretical study and other considerations, supports the hypothesis of overdominant selection (heterozygote advantage).

Animals↗

Structure, function, and evolution of mouse TL genes, nonclassical class I genes of the major histocompatibility complex.

In contrast to well-studied "classical" class I genes of the major histocompatibility complex (MHC), the biology of nonclassical class I genes remains largely unexamined. The mouse TL genes constitute one of the best defined systems among nonclassical class I genes in the T region of the MHC. To elucidate the function and the evolution of TL genes and their relationship to classical class I genes, seven TL DNA sequences, including one from a Japanese wild mouse, were examined and compared with those of several mouse and human classical class I genes. The TL genes differ from either classical class I genes or pseudogenes in the extent and pattern of nucleotide substitutions. Natural selection appears to have operated so as to preserve the function of TL, which might have been acquired in an early stage of its evolution. In a putative peptide-binding region encoded by TL genes, the rate of nonsynonymous (amino acid replacing) substitution is considerably lower than that of synonymous substitution. This conservation is completely opposite that in classical class I genes, in which the peptide-binding region has evolved to diversify amino acid sequences so as to recognize a variety of antigens. Thus, it is suggested that the function of TL antigens is distinct from that of classical class I antigens and is related to the recognition of a relatively restricted repertoire of antigens and their presentation to T-cell receptors.

Amino Acid Sequence↗

Positive Darwinian selection after gene duplication in primate ribonuclease genes.

Evolutionary mechanisms of origins of new gene function have been a subject of long-standing debate. Here we report a convincing case in which positive Darwinian selection operated at the molecular level during the evolution of novel function by gene duplication. The genes for eosinophil cationic protein (ECP) and eosinophil-derived neurotoxin (EDN) in primates belong to the ribonuclease gene family, and the ECP gene, whose product has an anti-pathogen function not displayed by EDN, was generated by duplication of the EDN gene about 31 million years ago. Using inferred nucleotide sequences of ancestral organisms, we showed that the rate of nonsynonymous nucleotide substitution was significantly higher than that of synonymous substitution for the ECP gene. This strongly suggests that positive Darwinian selection operated in the early stage of evolution of the ECP gene. It was also found that the number of arginine residues increased substantially in a short period of evolutionary time after gene duplication, and these amino acid changes probably produced the novel anti-pathogen function of ECP.

Amino Acid Sequence↗

Circumsporozoite protein genes of malaria parasites (Plasmodium spp.): evidence for positive selection on immunogenic regions.

The circumsporozoite (CS) protein is a cell surface protein of the sporozoite, the stage of the life cycle of malaria parasites (Plasmodium spp.) that infects the vertebrate host. Analysis of DNA sequences supports the hypothesis that in Plasmodium falciparum, positive Darwinian selection favors diversity in the T-cell epitopes (peptides presented to T cells by host MHC molecules) of the CS protein. In gene regions encoding T cell epitopes of P. falciparum, the rate of nonsynonymous nucleotide substitution is significantly higher than that of synonymous substitution, whereas this is not true of other gene regions. Furthermore nonsynonymous nucleotide substitutions in these regions cause a change of amino acid residue charge significantly more frequently than expected by chance. By contrast, in Plasmodium cynomolgi, the same regions show no evidence of positive selection, and residue charge is conserved. The CS protein has a central repeat region, which is the target of host antibodies. In P. falciparum, the amino acid sequence of the repeat region is conserved within and between alleles. In P. cynomolgi, on the other hand, there is evidence that positive selection has favored evolution of two different repeat types within a given allele.

Amino Acid Sequence↗

Higher frequency of concerted evolutionary events in rodents than in man at the polyubiquitin gene VNTR locus.

The polyubiquitin gene is an evolutionarily conserved eukaryotic gene, encoding tandemly repeated multiple ubiquitins, and is considered to be subject to concerted evolution. Here, we present the nucleotide sequences of new alleles of the polyubiquitin gene UbC in humans and CHUB2 in Chinese hamster, which encode a different number of ubiquitin units from those of previously reported genes. And we analyze the concerted evolution of these genes on the basis of their orthologous relationship. That the mean of the synonymous sequence difference Ks which is defined as the number of synonymous substitution relative to the total number of synonymous sites, within the UbC and CHUB2 genes (0.192 +/- 0.096) is significantly less than Ks between these genes (0.602 +/- 0.057) provides direct evidence for concerted evolution. Moreover, it also appears that concerted evolutionary events have been much more frequent in CHUB2 than in UbC, because Ks within CHUB2 (0.022 +/- 0.018) is much less than that within UbC (0.362 +/- 0.192). By a numerical simulation, postulating that the major mechanism of concerted evolution in polyubiquitin genes is unequal crossing over, we estimated the frequency of concerted evolutionary events of CHUB2 at 3.3 x 10(-5) per year and that of UbC at no more than 5.0 x 10(-7) per year.

Amino Acid Sequence↗

Adaptive diversification of bitter taste receptor genes in Mammalian evolution.

The diversity and evolution of bitter taste perception in mammals is not well understood. Recent discoveries of bitter taste receptor (T2R) genes provide an opportunity for a genetic approach to this question. We here report the identification of 10 and 30 putative T2R genes from the draft human and mouse genome sequences, respectively, in addition to the 23 and 6 previously known T2R genes from the two species. A phylogenetic analysis of the T2R genes suggests that they can be classified into three main groups, which are designated A, B, and C. Interestingly, while the one-to-one gene orthology between the human and mouse is common to group B and C genes, group A genes show a pattern of species- or lineage-specific duplication. It is possible that group B and C genes are necessary for detecting bitter tastants common to both humans and mice, whereas group A genes are used for species-specific bitter tastants. The analysis also reveals that phylogenetically closely related T2R genes are close in their chromosomal locations, demonstrating tandem gene duplication as the primary source of new T2Rs. For closely related paralogous genes, a rate of nonsynonymous nucleotide substitution significantly higher than the rate of synonymous substitution was observed in the extracellular regions of T2Rs, which are presumably involved in tastant-binding. This suggests the role of positive selection in the diversification of newly duplicated T2R genes. Because many natural poisonous substances are bitter, we conjecture that the mammalian T2R genes are under diversifying selection for the ability to recognize a diverse array of poisons that the organisms may encounter in exploring new habitats and diets.

Amino Acid Sequence↗

Adaptive evolution of the water stress-induced gene Asr2 in Lycopersicon species dwelling in arid habitats.

The Asr2 gene encodes a putative transcription factor that is up-regulated in leaves and roots of tomato plants exposed to water-deficit stress. This gene was first cloned and characterized in a cultivar of commercial tomato (Lycopersicon esculentum cv. Ailsa Craig). In this work, we report the complete coding sequences of the orthologous Asr2 genes in six wild tomato lineages: L. hirsutum, L. cheesmanii, L. esculentum v. cerasiforme, L. chilense, L. peruvianum v. humifusum and L. peruvianum f. glandulosum. Estimates of the Ka/Ks ratio (omega) in pairwise comparisons within the genus Lycopersicon were equal or greater than 1 (a signature of adaptive evolution) when involving L. chilense and L. peruvianum v. humifusum. Interestingly, these two species are distinct from the others in their adaptation to dry habitats. We also mapped the detected substitutions onto a phylogenetic tree of the genus Lycopersicon. Remarkably, there are two and three amino acid substitutions, which contrast with the absence of synonymous substitutions along the terminal branches leading to L. chilense and L. peruvianum v. humifusum, respectively. Likelihood ratio tests confirmed that omega values in the branches leading to these species are significantly different from the remaining branches of the tree. Moreover, inferred changes in the branches leading to these species that inhabit dry areas are nonconservative and may be associated with dramatic alterations in ASR2 protein conformation. In this work, we demonstrate accelerated rates of amino acid substitutions in the Asr2 gene of tomato lineages living in dry habitats, thus giving support to the hypothesis of adaptive Darwinian evolution.

Amino Acid Sequence↗

A Bayesian model comparison approach to inferring positive selection.

A popular approach to detecting positive selection is to estimate the parameters of a probabilistic model of codon evolution and perform inference based on its maximum likelihood parameter values. This approach has been evaluated intensively in a number of simulation studies and found to be robust when the available data set is large. However, uncertainties in the estimated parameter values can lead to errors in the inference, especially when the data set is small or there is insufficient divergence between the sequences. We introduce a Bayesian model comparison approach to infer whether the sequence as a whole contains sites at which the rate of nonsynonymous substitution is greater than the rate of synonymous substitution. We incorporated this probabilistic model comparison into a Bayesian approach to site-specific inference of positive selection. Using simulated sequences, we compared this approach to the commonly used empirical Bayes approach and investigated the effect of tree length on the performance of both methods. We found that the Bayesian approach outperforms the empirical Bayes method when the amount of sequence divergence is small and is less prone to false-positive inference when the sequences are saturated, while the results are indistinguishable for intermediate levels of sequence divergence.

Bayes Theorem↗

Variation in positive selection in termite GNBPs and Relish.

Social insects are model organisms for investigating molecular evolution in the innate immune system. Their diversity affords comparative analysis among closely related species, and group living is likely to contribute to the pathogen stress imposed on the immune system. We used different models of nucleotide substitution at nonsynonymous (amino acid altering) and synonymous (silent) sites to compare the different levels and type of selection among three immunity genes in 13 Australian termite species (Nasutitermes). The immunity genes include two encoding pathogen recognition proteins (gram-negative bacterial-binding proteins) that duplicated and diverged before or soon after the evolution of the termites and a transcription factor (Relish), which induces the production of antimicrobial peptides. A comparison of evolutionary models that assign four unrestricted classes of dN/dS (the ratio of the nonsynonymous to synonymous substitution rate) to different Nasutitermes lineages revealed that the occurrence of positive selection (dN/dS > 1) varies among lineages and the three genes. Positive selection appears to have driven the evolution of all three genes in an ancestral lineage of three subterranean termites. It had previously been suggested that there was a transition along this ancestral lineage to termite morphology and ecology associated with a diet of decayed wood, a diet that may expose termites to elevated levels of fungal and bacterial pathogens. Relish appears to have experienced the highest levels of selective pressure for change among all three genes. Positively selected sites in the molecule are located in regions that are important for its activation, which suggests that amino acid substitutions at these sites are a counter response to pathogen mechanisms that disrupt the activation of Relish.

Amino Acid Substitution↗

Molecular evolution of mammalian lactate dehydrogenase-A genes and pseudogenes: association of a mouse processed pseudogene with a B1 repetitive sequence.

A mouse genomic clone containing a lactate dehydrogenase-A (LDH-A) processed pseudogene and a B1 repetitive element was isolated, and a nucleotide sequence of approximately 3 kb was determined. The pseudogene and B1 element are flanked by perfect 13-bp repeats, and the B1 sequence starts at 14 nucleotides 3' to the presumptive polyadenylation signal of the pseudogene. The nucleotide sequences of the LDH-A genes and processed pseudogenes from mouse, rat, and human were compared, and a phylogenetic tree was constructed. The rate and pattern of nucleotide substitutions in the LDH-A pseudogenes are similar to previously reported results (Li et al. 1984). The average rate of nucleotide substitutions in the LDH-A pseudogenes is 4.3 X 10(-9)/site/year. The substitutions of C----T and G----A are most frequent, and A----G substitutions are relatively high. The rate of synonymous substitutions in the LDH-A genes is 5.3 X 10(-9), which is not significantly higher than the average rate of 4.7 X 10(-9) for 35 mammalian genes. The rate of nonsynonymous substitutions in the LDH-A genes is 0.20 X 10(-9), which is considerably lower than the average rate of 0.88 X 10(-9) for 35 mammalian genes. Thus, the mammalian LDH-A gene appears to be highly conserved in evolution.

Animals↗

Origin and evolution of new exons in rodents.

Gene number difference among organisms demonstrates that new gene origination is a fundamental biological process in evolution. Exon shuffling has been universally observed in the formation of new genes. Yet to be learned are the ways new exons originate and evolve, and how often new exons appear. To address these questions, we identified 2695 newly evolved exons in the mouse and rat by comparing the expressed sequences of 12,419 orthologous genes between human and mouse, using 743,856 pig ESTs as the outgroup. The new exon origination rate is about 2.71 x 10(-3) per gene per million years. These new exons have markedly accelerated rates both of nonsynonymous substitutions and of insertions/deletions (indels). A much higher proportion of new exons have K(a)/K(s) ratios >1 (where K(a) is the nonsynonymous substitution rate and K(s) is the synonymous substitution rate) than do the old exons shared by human and mouse, implying a role of positive selection in the rapid evolution. The majority of these new exons have sequences unique in the genome, suggesting that most new exons might originate through "exonization" of intronic sequences. Most of the new exons appear to be alternative exons that are expressed at low levels.

Alternative Splicing↗

Analyses of P-like transposable element sequences from the genome of Anopheles gambiae.

We have identified 50 P element-homologous sequences in the genome of Anopheles gambiae by performing homology searches against the public genome database of A. gambiae using the canonical P element from Drosophila melanogaster as a query sequence. While most of these sequences belong to P subfamilies previously described from anopheline mosquitoes, at least four new subfamilies were identified. One of these A. gambiae P elements, which we termed AgPLS, was analysed in detail. AgPLS consists of three exons and does not have inverted terminal repeats. This element retains several of the structural features of other P-encoded peptides, such as motifs involved in DNA-protein and protein-protein interaction, and a motif involved in GTP utilization. Strong sequence and structural similarity to functional P elements, a number of nonsynonymous substitutions that is smaller than that of synonymous substitutions and the presence of putative nuclear localization signals suggest that the A. gambiae elements may retain the capacity for transposition or its repression. These sequences seem to be most closely related to P elements described from Musca domestica and Lucilia cuprina, the only P element hosts known outside the family Drosophilidae.

Amino Acid Sequence↗

A new algorithm for analysis of within-host HIV-1 evolution.

A new algorithm for inferring the evolution of within-host viral sequences is presented. A sequential-linking approach is developed so that a longitudinal phylogenetic tree can be reconstructed from sequential molecular data that are obtained at different time points from the same host. The algorithm employs a codon-based model, which uses a Markov process to describe substitutions between codons, to calculate nonsynonymous and synonymous substitution rates and to distinguish positive selection and neutral evolution. The algorithm is applied to a data set of the V3 region of the HIV-1 envelope genes sequenced at different years after the infection of a single patient. The results suggest that this algorithm may provide a more realistic description of viral evolution than traditional evolutionary models, because it accounts for both neutral and adaptive evolution, and reconstructs a longitudinal phylogenetic tree that describes the dynamic process of viral evolution.

Algorithms↗

Overdispersed molecular clock at the major histocompatibility complex loci.

The extent of amino acid differences of major histocompatibility complex molecules within species is unusually high, consistent with the finding that some pairs of alleles have persisted for more than ten million years and the view that the polymorphism has been maintained by natural selection. The disparity between synonymous and non-synonymous substitutions in the antigen recognition site, however, suggests that some non-synonymous sites have undergone a number of substitutions whereas others have little or none. To describe statistically such an overdispersed underlying process, commonly used Poisson processes are inadequate. An alternative process leads to the surprising conclusion that each non-synonymous site has accumulated as many as 2.6 substitutions, on the average, in the two lineages leading to humans and mice. The standard deviation is also very large (6.6) and the dispersion index (the ratio of the variance to the mean) is at least 17. The substitution process thus inferred qualitatively agrees with the disposition (a boomerang pattern) of substitutions between HLA-A2 and Aw68 alleles, and quantitatively agrees well with that expected where the evolution of major histocompatibility complex molecules has long been driven mostly by balancing selection.

Alleles↗

Evolutionary conservation and selection of human disease gene orthologs in the rat and mouse genomes.

BACKGROUND: Model organisms have contributed substantially to our understanding of the etiology of human disease as well as having assisted with the development of new treatment modalities. The availability of the human, mouse and, most recently, the rat genome sequences now permit the comprehensive investigation of the rodent orthologs of genes associated with human disease. Here, we investigate whether human disease genes differ significantly from their rodent orthologs with respect to their overall levels of conservation and their rates of evolutionary change. RESULTS: Human disease genes are unevenly distributed among human chromosomes and are highly represented (99.5%) among human-rodent ortholog sets. Differences are revealed in evolutionary conservation and selection between different categories of human disease genes. Although selection appears not to have greatly discriminated between disease and non-disease genes, synonymous substitution rates are significantly higher for disease genes. In neurological and malformation syndrome disease systems, associated genes have evolved slowly whereas genes of the immune, hematological and pulmonary disease systems have changed more rapidly. Amino-acid substitutions associated with human inherited disease occur at sites that are more highly conserved than the average; nevertheless, 15 substituting amino acids associated with human disease were identified as wild-type amino acids in the rat. Rodent orthologs of human trinucleotide repeat-expansion disease genes were found to contain substantially fewer of such repeats. Six human genes that share the same characteristics as triplet repeat-expansion disease-associated genes were identified; although four of these genes are expressed in the brain, none is currently known to be associated with disease. CONCLUSIONS: Most human disease genes have been retained in rodent genomes. Synonymous nucleotide substitutions occur at a higher rate in disease genes, a finding that may reflect increased mutation rates in the chromosomal regions in which disease genes are found. Rodent orthologs associated with neurological function exhibit the greatest evolutionary conservation; this suggests that rodent models of human neurological disease are likely to most faithfully represent human disease processes. However, with regard to neurological triplet repeat expansion-associated human disease genes, the contraction, relative to human, of rodent trinucleotide repeats suggests that rodent loci may not achieve a 'critical repeat threshold' necessary to undergo spontaneous pathological repeat expansions. The identification of six genes in this study that have multiple characteristics associated with repeat expansion-disease genes raises the possibility that not all human loci capable of facilitating neurological disease by repeat expansion have as yet been identified.

Animals↗

The effect of tandem substitutions on the correlation between synonymous and nonsynonymous rates in rodents.

Nonsynonymous substitutions in DNA cause amino acid substitutions while synonymous substitutions in DNA leave amino acids unchanged. The cause of the correlation between the substitution rates at nonsynonymous (K(A)) and synonymous (K(S)) sites in mammals is a contentious issue, and one that impacts on many aspects of molecular evolution. Here we use a large set of orthologous mammalian genes to investigate the causes of the K(A)-K(S) correlation in rodents. The strength of the K(A)-K(S) correlation exceeds the neutral theory expectation when substitution rates are estimated using algorithmic methods, but not when substitution rates are estimated by maximum likelihood. Irrespective of this methodological uncertainty the strength of the K(A)-K(S) correlation appears mostly due to tandem substitutions, an excess of which is generated by substitutional nonindependence. Doublet mutations cannot explain the excess of tandem synonymous-nonsynonymous substitutions, and substitution patterns indicate that selection on silent sites is the likely cause. We find no evidence for selection on codon usage. The nature of the relationship between synonymous divergence and base composition is unclear because we find a significant correlation if we use maximum-likelihood methods but not if we use algorithmic methods. Finally, we find that K(S) is reduced at the start of genes, which suggests that selection for RNA structure may affect silent sites in mammalian protein-coding genes.

Algorithms↗

A dissection of volatility in yeast.

It has been suggested that volatility, the proportion of mutations which change an amino acid, can be used to infer the level of natural selection acting upon a gene. This conjecture is supported by a correlation between volatility and the rate of nonsynonymous substitution (dN), or the ratio of nonsynonymous and synonymous substitution rates, in a variety of organisms. These organisms include yeast, in which the correlations are quite strong. Here we show that these correlations are a by-product of a correlation between synonymous codon bias toward translationally optimal codons and dN. Although this analysis suggests that volatility is not a good measure of the selection, we suggest that it might be possible to infer something about the level of natural selection, from a single genome sequence, using translational codon bias.

Base Sequence↗

Differentiated evolutionary rates in alternative exons and the implications for splicing regulation.

BACKGROUND: Alternatively spliced exons play an important role in the diversification of gene function in most metazoans and are highly regulated by conserved motifs in exons and introns. Two contradicting properties have been associated to evolutionary conserved alternative exons: higher sequence conservation and higher rate of non-synonymous substitutions, relative to constitutive exons. In order to clarify this issue, we have performed an analysis of the evolution of alternative and constitutive exons, using a large set of protein coding exons conserved between human and mouse and taking into account the conservation of the transcript exonic structure. Further, we have also defined a measure of the variation of the arrangement of exonic splicing enhancers (ESE-conservation score) to study the evolution of splicing regulatory sequences. We have used this measure to correlate the changes in the arrangement of ESEs with the divergence of exon and intron sequences. RESULTS: We find evidence for a relation between the lack of conservation of the exonic structure and the weakening of the sequence evolutionary constraints in alternative and constitutive exons. Exons in transcripts with non-conserved exonic structures have higher synonymous (dS) and non-synonymous (dN) substitution rates than exons in conserved structures. Moreover, alternative exons in transcripts with non-conserved exonic structure are the least constrained in sequence evolution, and at high EST-inclusion levels they are found to be very similar to constitutive exons, whereas alternative exons in transcripts with conserved exonic structure have a dS significantly lower than average at all EST-inclusion levels. We also find higher conservation in the arrangement of ESEs in constitutive exons compared to alternative ones. Additionally, the sequence conservation at flanking introns remains constant for constitutive exons at all ESE-conservation values, but increases for alternative exons at high ESE-conservation values. CONCLUSION: We conclude that most of the differences in dN observed between alternative and constitutive exons can be explained by the conservation of the transcript exonic structure. Low dS values are more characteristic of alternative exons with conserved exonic structure, but not of those with non-conserved exonic structure. Additionally, constitutive exons are characterized by a higher conservation in the arrangement of ESEs, and alternative exons with an ESE-conservation similar to that of constitutive exons are characterized by a conservation of the flanking intron sequences higher than average, indicating the presence of more intronic regulatory signals.

Alternative Splicing↗