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 343 records · Page 19Linked to original sources
PubMed2001

Bushbaby growth hormone is much more similar to nonprimate growth hormones than to rhesus monkey and human growth hormones.

Unlike other mammals, Old World primates have five growth hormone-like genes that are highly divergent at the amino acid level from the single growth hormone genes found in nonprimates. Additionally, there is a change in the interaction of growth hormone with its receptor in humans such that human growth hormone functions in nonprimates, whereas nonprimate growth hormone is ineffective in humans. A Southern blotting analysis of the genome of a prosimian, Galago senegalensis, revealed a single growth hormone locus. This single gene was PCR-amplified from genomic DNA and sequenced. It has a rate of nonsynonymous nucleotide substitution less than one fourth that of the human growth hormone gene, while the rates of synonymous substitution in the two species are less different. Human and rhesus monkey growth hormones exhibit variation at a number of amino acid residues that can affect receptor binding. The galago growth hormone is conservative at each of these sites, indicating that this growth hormone is functionally like nonprimate growth hormones. These observations indicate that the amplification and rapid divergence of primate growth hormones occurred after the separation of the higher primate lineage from the galago lineage.

Amino Acid Sequence↗
PubMed2000

Estimating synonymous and nonsynonymous substitution rates under realistic evolutionary models.

Approximate methods for estimating the numbers of synonymous and nonsynonymous substitutions between two DNA sequences involve three steps: counting of synonymous and nonsynonymous sites in the two sequences, counting of synonymous and nonsynonymous differences between the two sequences, and correcting for multiple substitutions at the same site. We examine complexities involved in those steps and propose a new approximate method that takes into account two major features of DNA sequence evolution: transition/transversion rate bias and base/codon frequency bias. We compare the new method with maximum likelihood, as well as several other approximate methods, by examining infinitely long sequences, performing computer simulations, and analyzing a real data set. The results suggest that when there are transition/transversion rate biases and base/codon frequency biases, previously described approximate methods for estimating the nonsynonymous/synonymous rate ratio may involve serious biases, and the bias can be both positive and negative. The new method is, in general, superior to earlier approximate methods and may be useful for analyzing large data sets, although maximum likelihood appears to always be the method of choice.

Animals↗
PubMed2004

Host-symbiont stability and fast evolutionary rates in an ant-bacterium association: cospeciation of camponotus species and their endosymbionts, candidatus blochmannia.

Bacterial endosymbionts are widespread across several insect orders and are involved in interactions ranging from obligate mutualism to reproductive parasitism. Candidatus Blochmannia gen. nov. (Blochmannia) is an obligate bacterial associate of Camponotus and related ant genera (Hymenoptera: Formicidae). The occurrence of Blochmannia in all Camponotus species sampled from field populations and its maternal transmission to host offspring suggest that this bacterium is engaged in a long-term, stable association with its ant hosts. However, evidence for cospeciation in this system is equivocal because previous phylogenetic studies were based on limited gene sampling, lacked statistical analysis of congruence, and have even suggested host switching. We compared phylogenies of host genes (the nuclear EF-1alphaF2 and mitochondrial COI/II) and Blochmannia genes (16S ribosomal DNA [rDNA], groEL, gidA, and rpsB), totaling more than 7 kilobases for each of 16 Camponotus species. Each data set was analyzed using maximum likelihood and Bayesian phylogenetic reconstruction methods. We found minimal conflict among host and symbiont phylogenies, and the few areas of discordance occurred at deep nodes that were poorly supported by individual data sets. Concatenated protein-coding genes produced a very well-resolved tree that, based on the Shimodaira-Hasegawa test, did not conflict with any host or symbiont data set. Correlated rates of synonymous substitution (d(S)) along corresponding branches of host and symbiont phylogenies further supported the hypothesis of cospeciation. These findings indicate that Blochmannia-Camponotus symbiosis has been evolutionarily stable throughout tens of millions of years. Based on inferred divergence times among the ant hosts, we estimated rates of sequence evolution of Blochmannia to be approximately 0.0024 substitutions per site per million years (s/s/MY) for the 16S rDNA gene and approximately 0.1094 s/s/MY at synonymous positions of the genes sampled. These rates are several-fold higher than those for related bacteria Buchnera aphidicola and Escherichia coli. Phylogenetic congruence among Blochmannia genes indicates genome stability that typifies primary endosymbionts of insects.

Animals↗
PubMed1993

Pattern of nucleotide substitutions in growth hormone-prolactin gene family: a paradigm for evolution by gene duplication.

The growth hormone-prolactin gene family in mammals is an interesting example of evolution by gene duplication. Divergence among members of duplicated gene families and among species was examined by using reported gene sequences of growth hormone, prolactin and their receptors. Sequence divergence among species was found to show a general tendency in which a generation-time effect is pronounced for synonymous substitutions but not so for nonsynonymous substitutions. Divergence among duplicated genes is characterized by the relatively high rate of nonsynonymous substitutions, i.e., the rate is close to that of synonymous ones. In view of the stage- and tissue-specific expression of duplicated genes, some of the amino acid substitutions among duplicated genes is likely to be caused by positive Darwinian selection.

Amino Acid Sequence↗
PubMed1998

Evolution and phylogenetic utility of the period gene in Lepidoptera.

Evolution and phylogenetic utility of the period gene are explored through sequence analysis of a relatively conserved 909-bp fragment in 26 lepidopteran species. Taxa range from tribes to superfamilies, primarily within the putative clade Macrolepidotera plus near outgroups, and include both strongly established and problematic groupings. Their divergence dates probably range from the late Cretaceous through much of the Tertiary. Comparisons within the same set of closely related species show that amino acid substitutions in period occur 4.9 and 44 times as frequently as they do in two other nuclear genes--dopa decarboxylase and elongation factor-1 alpha, respectively. In contrast, rates of observed synonymous substitution are within 60% of each other for these three genes. Synonymous changes in period approach saturation by the family level, whereas nonsynonymous and amino acid divergences across the Macrolepidoptera are less than half the maximal values reported for this gene. Phylogenetic analyses of period strongly supported groupings at the family level and below. In contrast to previous analyses at this level with other nuclear genes, much of the information lies in nonsynonymous change. Relationships up to the superfamily level were recovered with decreasing effectiveness, and little, if any, signal was apparent regarding relationships among superfamilies. This could reflect rapid radiation of the superfamilies, however, rather than saturation in the period locus; thus, period, in combination with other genes, remains a plausible candidate for approaching the difficult problems of lepidopteran family and superfamily relationships.

Amino Acid Sequence↗
PubMed2001

The acetylcholinesterase gene and organophosphorus resistance in the Australian sheep blowfly, Lucilia cuprina.

Acetylcholinesterase (AChE), encoded by the Ace gene, is the primary target of organophosphorous (OP) and carbamate insecticides. Ace mutations have been identified in OP resistants strains of Drosophila melanogaster. However, in the Australian sheep blowfly, Lucilia cuprina, resistance in field and laboratory generated strains is determined by point mutations in the Rop-1 gene, which encodes a carboxylesterase, E3. To investigate the apparent bias for the Rop-1/E3 mechanism in the evolution of OP resistance in L. cuprina, we have cloned the Ace gene from this species and characterized its product. Southern hybridization indicates the existence of a single Ace gene in L. cuprina. The amino acid sequence of L. cuprina AChE shares 85.3% identity with D. melanogaster and 92.4% with Musca domestica AChE. Five point mutations in Ace associated with reduced sensitivity to OP insecticides have been previously detected in resistant strains of D. melanogaster. These residues are identical in susceptible strains of D. melanogaster and L. cuprina, although different codons are used. Each of the amino acid substitutions that confer OP resistance in D. melanogaster could also occur in L. cuprina by a single non-synonymous substitution. These data suggest that the resistance mechanism used in L. cuprina is determined by factors other than codon bias. The same point mutations, singly and in combination, were introduced into the Ace gene of L. cuprina by site-directed mutagenesis and the resulting AChE enzymes expressed using a baculovirus system to characterise their kinetic properties and interactions with OP insecticides. The K(m) of wild type AChE for acetylthiocholine (ASCh) is 23.13 microM and the point mutations change the affinity to the substrate. The turnover number of Lucilia AChE for ASCh was estimated to be 1.27x10(3) min(-1), similar to Drosophila or housefly AChE. The single amino acid replacements reduce the affinities of the AChE for OPs and give up to 8.7-fold OP insensitivity, while combined mutations give up to 35-fold insensitivity. However, other published studies indicate these same mutations yield higher levels of OP insensitivity in D. melanogaster and A. aegypti. The inhibition data indicate that the wild type form of AChE of L. cuprina is 12.4-fold less sensitive to OP inhibition than the susceptible form of E3, suggesting that the carboxylesterases may have a role in the protection of AChE via a sequestration mechanism. This provides a possible explanation for the bias towards the evolution of resistance via the Rop-1/E3 mechanism in L. cuprina.

Acetylcholinesterase↗
PubMed1988

Molecular evolution of the Escherichia coli chromosome. I. Analysis of structure and natural variation in a previously uncharacterized region between trp and tonB.

We present the sequence of a 3500-bp region of the Escherichia coli strain K12 chromosome lying between the tryptophan operon and the tonB gene. Analysis of the sequence yields six open reading frames that have properties characteristic of genes for proteins. The reading frames are closely spaced, and putative transcription units and control sites compose over 95% of the DNA. The sequences of several wild strains of E. coli have been determined for a large segment of the region described. Comparison of these sequences reveals the effects of base substitutions, DNA rearrangements, and recombination. In the regions presumably expressed as polypeptides, most of the natural variation results from synonymous substitutions. However, the DNA rearrangements identified have end points within the open reading frames and disrupt them in a variety of ways. The effects of genetic recombination between strains, recently found to be significant on a large scale in E. coli, are also apparent in the region between trp and tonB.

Base Sequence↗
PubMed2001

Noncoding RNA gene detection using comparative sequence analysis.

BACKGROUND: Noncoding RNA genes produce transcripts that exert their function without ever producing proteins. Noncoding RNA gene sequences do not have strong statistical signals, unlike protein coding genes. A reliable general purpose computational genefinder for noncoding RNA genes has been elusive. RESULTS: We describe a comparative sequence analysis algorithm for detecting novel structural RNA genes. The key idea is to test the pattern of substitutions observed in a pairwise alignment of two homologous sequences. A conserved coding region tends to show a pattern of synonymous substitutions, whereas a conserved structural RNA tends to show a pattern of compensatory mutations consistent with some base-paired secondary structure. We formalize this intuition using three probabilistic "pair-grammars": a pair stochastic context free grammar modeling alignments constrained by structural RNA evolution, a pair hidden Markov model modeling alignments constrained by coding sequence evolution, and a pair hidden Markov model modeling a null hypothesis of position-independent evolution. Given an input pairwise sequence alignment (e.g. from a BLASTN comparison of two related genomes) we classify the alignment into the coding, RNA, or null class according to the posterior probability of each class. CONCLUSIONS: We have implemented this approach as a program, QRNA, which we consider to be a prototype structural noncoding RNA genefinder. Tests suggest that this approach detects noncoding RNA genes with a fair degree of reliability.

Algorithms↗
PubMed2004

Moderate mutation rate in the SARS coronavirus genome and its implications.

BACKGROUND: The outbreak of severe acute respiratory syndrome (SARS) caused a severe global epidemic in 2003 which led to hundreds of deaths and many thousands of hospitalizations. The virus causing SARS was identified as a novel coronavirus (SARS-CoV) and multiple genomic sequences have been revealed since mid-April, 2003. After a quiet summer and fall in 2003, the newly emerged SARS cases in Asia, particularly the latest cases in China, are reinforcing a wide-spread belief that the SARS epidemic would strike back. With the understanding that SARS-CoV might be with humans for years to come, knowledge of the evolutionary mechanism of the SARS-CoV, including its mutation rate and emergence time, is fundamental to battle this deadly pathogen. To date, the speed at which the deadly virus evolved in nature and the elapsed time before it was transmitted to humans remains poorly understood. RESULTS: Sixteen complete genomic sequences with available clinical histories during the SARS outbreak were analyzed. After careful examination of multiple-sequence alignment, 114 single nucleotide variations were identified. To minimize the effects of sequencing errors and additional mutations during the cell culture, three strategies were applied to estimate the mutation rate by 1) using the closely related sequences as background controls; 2) adjusting the divergence time for cell culture; or 3) using the common variants only. The mutation rate in the SARS-CoV genome was estimated to be 0.80 - 2.38 x 10-3 nucleotide substitution per site per year which is in the same order of magnitude as other RNA viruses. The non-synonymous and synonymous substitution rates were estimated to be 1.16 - 3.30 x 10-3 and 1.67 - 4.67 x 10-3 per site per year, respectively. The most recent common ancestor of the 16 sequences was inferred to be present as early as the spring of 2002. CONCLUSIONS: The estimated mutation rates in the SARS-CoV using multiple strategies were not unusual among coronaviruses and moderate compared to those in other RNA viruses. All estimates of mutation rates led to the inference that the SARS-CoV could have been with humans in the spring of 2002 without causing a severe epidemic.

Coronavirus↗
PubMed2006

Characterization of two paralogous muscleblind-like genes from the tiger pufferfish (Takifugu rubripes).

Muscleblind-like (Mbnl) proteins are required for terminal muscle differentiation in mammals. In this study we have identified two mbnl paralogues from the tiger pufferfish, tmbnl2a and tmbnl3, which are the first examples of non-mammalian mbnl genes. Tmbnl2a and tmbnl3 were found in regions of conserved synteny and had a high degree of global conservation with their mammalian homologues. Phylogenetic analysis showed that the T. rubripes genome contains one mbnl3 gene and two copies of mbnl1 and mbnl2. Moreover, the mbnl1 and mbnl3 paralogues are derived from duplication of a common ancestral gene. The average rates of synonymous substitutions between T. rubripes, mouse and human mbnl2 and mbnl3 genes were much higher than the corresponding rates of non-synonymous mutations, suggesting that Mbnl2 and Mbnl3 are subjected to strong purifying selection. Quantitation of tmbnl2a and tmbnl3 transcripts by real-time PCR revealed that these two paralogues are differentially expressed in fast and slow myotomal muscle, heart, liver, skin, brain and testes. Tmbnl2a was expressed at similar levels in all tissues examined, as was the mouse orthologue. Tmbnl3 was expressed at higher levels than tmbnl2a, with a ubiquitous tissue distribution. Expression of tmbnl3 remained high in adult pufferfish muscle whereas the mouse orthologue was down-regulated in adults, perhaps reflecting the indeterminate and determinate growth patterns of these taxa, respectively.

Amino Acid Sequence↗
PubMed2006

SWAKK: a web server for detecting positive selection in proteins using a sliding window substitution rate analysis.

We present a bioinformatic web server (SWAKK) for detecting amino acid sites or regions of a protein under positive selection. It estimates the ratio of non-synonymous to synonymous substitution rates (K(A)/K(S)) between a pair of protein-coding DNA sequences, by sliding a 3D window, or sphere, across one reference structure. The program displays the results on the 3D protein structure. In addition, for comparison or when a reference structure is unavailable, the server can also perform a sliding window analysis on the primary sequence. The SWAKK web server is available at http://oxytricha.princeton.edu/SWAKK/.

Amino Acid Substitution↗
PubMed2002

Gene location and bacterial sequence divergence.

Previous comparison of a relatively small set of homologous genes from Escherichia coli and Salmonella typhimurium revealed that genes nearer to the origin of replication had substitution rates lower than genes closer to the replication terminus. The recently completed sequences of numerous bacterial genomes have allowed us to test whether this effect of distance from the replication origin on substitution rates, as observed for the E. coli-S. typhimurium comparison, is a general feature of bacterial genomes. Extending the analysis to all 3,000 E. coli-S. typhimurium homologs confirmed the significant association between chromosomal position and synonymous site divergence. However, the effect, though still significant, is not as dramatic as originally thought. A similar association between relative chromosomal location and synonymous substitution rate was detected in the majority of other bacterial species comparisons within alpha- and gamma- Proteobacteria, and Firmicutes but was absent in Chlamydiales. The opposite trend, i.e., a decrease in synonymous divergence with distance from the replication origin, was detected in Mycobacteria. Analysis of the patterns of nucleotide substitutions revealed that the distance effect is not affected by gene orientation and is mainly caused by an increase in rates of transversions, suggesting that this effect may not be caused by recombinational repair or biased gene conversion, as originally suggested.

Bacteria↗
PubMed1994

Interspecific sequence comparison of the muscle-myosin heavy-chain genes from Drosophila hydei and Drosophila melanogaster.

The muscle-myosin heavy-chain (mMHC) gene of Drosophila hydei has been sequenced completely (size 23.3 kb). The sequence comparison with the D. melanogaster mMHC gene revealed that the exon-intron pattern is identical. The protein coding regions show a high degree of conservation (97%). The alternatively spliced exons (3a-b, 7a-d, 9a-c, 11a-e, and 15a-b) display more variations in the number of nonsynonymous and synonymous substitutions than the common exons (2, 4, 5, 6, 8, 10, 12, 13, 14, 16, 17, and 19). The base composition at synonymous sites of fourfold degenerate codons (third position) is not biased in the alternative exons. In the common exons there exists a bias for C and against A. These findings imply that the alternative exons of the Drosophila mMHC gene evolve at a different, in several cases higher, rate than the common ones. The 5' splice junctions and 5' and 3' untranslated regions show a high level of similarity, indicating a functional constraint on these sequences. The intron regions vary considerably in length within one species, but the corresponding introns are very similar in length between the two species and all contain stretches of sequence similarity. A particular example is the first intron, which contains multiple regions of similarity. In the conserved regions of intron 12 (head-tail border) sequences were found which have the potential to direct another smaller mMHC transcript.

Amino Acid Sequence↗
PubMed1993

Structure, function, and evolution of the family of superoxide dismutase proteins from halophilic archaebacteria.

The protein sequences of seven members of the superoxide dismutase (SOD) family from halophilic archaebacteria have been aligned and compared with each other and with the homologous Mn and Fe SOD sequences from eubacteria and the methanogenic archaebacterium Methanobacterium thermoautotrophicum. Of 199 common residues in the SOD proteins from halophilic archaebacteria, 125 are conserved in all seven sequences, and 64 of these are encoded by single unique triplets. The 74 remaining positions exhibit a high degree of variability, and for almost half of these, the encoding triplets are connected by at least two nonsynonymous nucleotide substitutions. The majority of nucleotide substitutions within the seven genes are nonsynonymous and result in amino acid replacement in the respective protein; silent third-codon-position (synonymous) substitutions are unexpectedly rare. Halophilic SODs contain 30 specific residues that are not found at the corresponding positions of the methanogenic or eubacterial SOD proteins. Seven of these are replacements of highly conserved amino acids in eubacterial SODs that are believed to play an important role in the three-dimensional structure of the protein. Residues implicated in formation of the active site, catalysis, and metal ion binding are conserved in all Mn and Fe SODs. Molecular phylogenies based on parsimony and neighbor-joining methods coherently group the halophile sequences but surprisingly fail to distinguish between the Mn SOD of Escherichia coli and the Fe SOD of M. thermoautotrophicum as the outgroup. These comparisons indicate that as a group, the SODs of halophilic archaebacteria have many unique and characteristic features. At the same time, the patterns of nucleotide substitution and amino acid replacement indicate that these genes and the proteins that they encode continue to be subject to strong and changing selection. This selection may be related to the presence of oxygen radicals and the inter- and intracellular composition and concentration of metal cations.

Amino Acid Sequence↗
PubMed1997

Role of random genetic drift in the evolution of interactive systems.

Through comparative studies of DNA sequences it has become possible to test the neutral and the selection theories of molecular evolution. The separate estimation of the numbers of synonymous and non-synonymous substitutions is one of the most powerful tools for detecting selection. The patterns on the average and variance of these two types of substitutions of mammalian genes turned out to be in accord with the slightly deleterious or the nearly neutral mutation theory for non-synonymous changes. Interactive systems at the amino acid level were suggested to be responsible for such nearly neutral or very weak selection. An attractive model is the NK model of Kauffman, which assumes that each amino acid makes a fitness contribution that depends upon the amino acid and upon K other amino acids among the N that make the protein. It is known that the fitness landscape is very rugged for K > or = 2. Population genetic analysis of this model suggest that protein evolution obeys the nearly neutral theory and that random genetic drift is important. In other words, evolution becomes rapid in small populations because proportion of near-neutrality increases among new mutations, and proteins as interactive systems evolve by shifting through random genetic drift on the multipeaked fitness landscape.

Animals↗
PubMed1995

The MHC E locus in macaques is polymorphic and is conserved between macaques and humans.

Although the functions of the molecules encoded by the classical MHC class I loci are well defined, no function has been ascribed to the molecules encoded by the non-classical MHC class I loci. To investigate the evolution and conservation of the non-classical loci, we cloned and sequenced HLA-E homologues in macaques. We isolated four E locus alleles from five rhesus monkeys and two E locus alleles from one cynomolgus monkey, which indicated that the E locus in macaques is polymorphic. We also compared the rate of nucleotide substitution in the second intron of the macaque and human E locus alleles with that of exons two and three. The rate of nucleotide substitution was significantly higher in the introns, which suggested that the E locus has evolved under selective pressure. Additionally, comparison of the rates of synonymous and non-synonymous substitutions in the peptide binding region versus the remainder of the molecule suggested that the codons encoding the amino acids in the peptide binding region had been conserved in macaques and humans over the 36 million years since macaques and humans last shared a common ancestor.

Amino Acid Sequence↗
PubMed2003

Dramatically elevated rate of mitochondrial substitution in lice (Insecta: Phthiraptera).

Few estimates of relative substitution rates, and the underlying mutation rates, exist between mitochondrial and nuclear genes in insects. Previous estimates for insects indicate a 2-9 times faster substitution rate in mitochondrial genes relative to nuclear genes. Here we use novel methods for estimating relative rates of substitution, which incorporate multiple substitutions, and apply these methods to a group of insects (lice, Order: Phthiraptera). First, we use a modification of copath analysis (branch length regression) to construct independent comparisons of rates, consisting of each branch in a phylogenetic tree. The branch length comparisons use maximum likelihood models to correct for multiple substitution. In addition, we estimate codon-specific rates under maximum likelihood for the different genes and compare these values. Estimates of the relative synonymous substitution rates between a mitochondrial (COI) and nuclear (EF-1alpha) gene in lice indicate a relative rate of several 100 to 1. This rapid relative mitochondrial rate (>100 times) is at least an order of magnitude faster than previous estimates for any group of organisms. Comparisons using the same methods for another group of insects (aphids) reveals that this extreme relative rate estimate is not simply attributable to the methods we used, because estimates from aphids are substantially lower. Taxon sampling affects the relative rate estimate, with comparisons involving more closely related taxa resulting in a higher estimate. Relative rate estimates also increase with model complexity, indicating that methods accounting for more multiple substitution estimate higher relative rates.

Animals↗
PubMed1987

cDNA sequence of a human skeletal muscle ADP/ATP translocator: lack of a leader peptide, divergence from a fibroblast translocator cDNA, and coevolution with mitochondrial DNA genes.

We have characterized a 1400-nucleotide cDNA for the human skeletal muscle ADP/ATP translocator. The deduced amino acid sequence is 94% homologous to the beef heart ADP/ATP translocator protein and contains only a single additional amino-terminal methionine. This implies that the human translocator lacks an amino-terminal targeting peptide, a conclusion substantiated by measuring the molecular weight of the protein synthesized in vitro. A 1400-nucleotide transcript encoding the skeletal muscle translocator was detected on blots of total RNA from human heart, kidney, skeletal muscle, and HeLa cells by hybridization with oligonucleotide probes homologous to the coding region and 3' noncoding region of the cDNA. However, the level of this mRNA varied substantially among tissues. Comparison of our skeletal muscle translocator sequence with that of a recently published human fibroblast translocator cognate revealed that the two proteins are 88% identical and diverged about 275 million years ago. Hence, tissues vary both in the level of expression of individual translocator genes and in differential expression of cognate translocator genes. Comparison of the base substitution rates of the ADP/ATP translocator and the oxidative phosphorylation genes encoded by mitochondrial DNA revealed that the mitochondrial DNA genes fix 10 times more synonymous substitutions and 12 times more replacement substitutions; yet, these nuclear and cytoplasmic respiration genes experience comparable evolutionary constraints. This suggests that the mitochondrial DNA genes are highly prone to deleterious mutations.

Base Sequence↗