Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “synonymous codon”

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 757 records · Page 42Linked to original sources

Protein evolution and codon usage bias on the neo-sex chromosomes of Drosophila miranda.

The neo-sex chromosomes of Drosophila miranda constitute an ideal system to study the effects of recombination on patterns of genome evolution. Due to a fusion of an autosome with the Y chromosome, one homolog is transmitted clonally. Here, I compare patterns of molecular evolution of 18 protein-coding genes located on the recombining neo-X and their homologs on the nonrecombining neo-Y chromosome. The rate of protein evolution has significantly increased on the neo-Y lineage since its formation. Amino acid substitutions are accumulating uniformly among neo-Y-linked genes, as expected if all loci on the neo-Y chromosome suffer from a reduced effectiveness of natural selection. In contrast, there is significant heterogeneity in the rate of protein evolution among neo-X-linked genes, with most loci being under strong purifying selection and two genes showing evidence for adaptive evolution. This observation agrees with theory predicting that linkage limits adaptive protein evolution. Both the neo-X and the neo-Y chromosome show an excess of unpreferred codon substitutions over preferred ones and no difference in this pattern was observed between the chromosomes. This suggests that there has been little or no selection maintaining codon bias in the D. miranda lineage. A change in mutational bias toward AT substitutions also contributes to the decline in codon bias. The contrast in patterns of molecular evolution between amino acid mutations and synonymous mutations on the neo-sex-linked genes can be understood in terms of chromosome-specific differences in effective population size and the distribution of selective effects of mutations.

Animals↗

Characterization and phylogenetic relationship of prosimian MHC class I genes.

MHC class I cDNA sequences from the most divergent primate group of extant primates compared to human, the suborder Strepsirrhini (prosimians), are described. The sequences are derived from the gray mouse lemur (Microcebus murinus) and the ring-tailed lemur (Lemur catta), which are members of the malagasy Lemuriformes, as well as from the pygmy slow loris (Nycticebus pygmaeus), a prosimian from East Asia. The M. murinus sequences have been analyzed in detail. Analysis of the expression level, G/C content, and synonymous vs. nonsynonymous substitution rates in the peptide-binding region codons suggests that these cDNA clones represent classical class I (class Ia) genes. According to Southern blot analysis, the genome of the gray mouse lemur might contain about 10 class I genes. In gene tree analysis, the strepsirrhine class Ia genes described here cluster significantly separately from the known class I genes of Catarrhini (humans, apes, Old World monkeys) and Platyrrhini (New World monkeys) species, suggesting that the class I loci of Simiiformes arose by gene duplications which occurred after the divergence of prosimians.

Amino Acid Sequence↗

Correlations of nucleotide substitution rates and base composition of mammalian coding sequences with protein structure.

We investigated the relationships between the nucleotide substitution rates and the predicted secondary structures in the three states representation (alpha-helix, beta-sheet, and coil). The analysis was carried out on 34 alignments, each of which comprised sequences belonging to at least four different mammalian orders. The rates of synonymous substitution were found to be significantly different in regions predicted to be alpha-helix, beta-sheet, or coil. Likewise, the nonsynonymous rates also differ, although expectedly at a lower extent, in the three types of secondary structure, suggesting that different selective constraints associated with the different structures are affecting in a similar way the synonymous and nonsynonymous rates. Moreover, the base composition of the third codon positions is different in coding sequence regions corresponding to different secondary structures of proteins.

Animals↗

Elevated evolutionary rates in the laboratory strain of Saccharomyces cerevisiae.

By using the maximum likelihood method, we made a genome-wide comparison of the evolutionary rates in the lineages leading to the laboratory strain (S288c) and a wild strain (YJM789) of Saccharomyces cerevisiae and found that genes in the laboratory strain tend to evolve faster than in the wild strain. The pattern of elevated evolution suggests that relaxation of selection intensity is the dominant underlying reason, which is consistent with recurrent bottlenecks in the S. cerevisiae laboratory strain population. Supporting this conclusion are the following observations: (i) the increases in nonsynonymous evolutionary rate occur for genes in all functional categories; (ii) most of the synonymous evolutionary rate increases in S288c occur in genes with strong codon usage bias; (iii) genes under stronger negative selection have a larger increase in nonsynonymous evolutionary rate; and (iv) more genes with adaptive evolution were detected in the laboratory strain, but they do not account for the majority of the increased evolution. The present discoveries suggest that experimental and possible industrial manipulations of the laboratory strain of yeast could have had a strong effect on the genetic makeup of this model organism. Furthermore, they imply an evolution of laboratory model organisms away from their wild counterparts, questioning the relevancy of the models especially when extensive laboratory cultivation has occurred. In addition, these results shed light on the evolution of livestock and crop species that have been under human domestication for years.

Biological Evolution↗

Transfer RNAs with complementary anticodons: could they reflect early evolution of discriminative genetic code adaptors?

In accordance with the hypercycle theory of M. Eigen and P. Schuster [(1979) Hypercycle: A Principle of Natural Self-Organization (Springer, New York)], the ancestors of modern tRNAs appear to have emerged via the shortest possible way, both complementary strands of a short symmetrical double helix serving as pre-tRNAs with complementary anticodons. This conclusion is based upon results of comparative sequence analysis of the 17-base-long anticodon loop and stem of tRNAs totaling 896 and especially of 22 pairs of consensus tRNAs with complementary or quasi-complementary anticodons. With regard to the anticodon loop and stem of pairs of consensus tRNAs, complementary distances were considerably less than direct distances--i.e., antiparallel pairing invariably yielded fewer mismatches than direct pairing. Furthermore, the smallest complementary distance was detected when two antiparallel sequences formed irregular G-U bonds in their anticodon triplets. The above implies that pre-tRNAs in peribiotic times were long hairpin structures having 73 bases or more, the middle base of an anticodon being the center of symmetry. Accordingly, each pair of pre-tRNAs with complementary anticodons should have been almost identical with each other except for their three central bases. The above situation appears to have dictated the early establishment of direct links between anticodons and the type of amino acids with which tRNAs are to be charged. This direct link is still maintained between modern aminoacyl-tRNA synthetases and anticodons. Replication of the double helices concertedly generated new codons for the same pair of amino acids. Thus, occurrence of synonymous as well as certain "palindromic" features of the genetic code table might have been determined by this mechanism.

Anticodon↗

The population genetics of the origin and divergence of the Drosophila simulans complex species.

The origins and divergence of Drosophila simulans and close relatives D. mauritiana and D. sechellia were examined using the patterns of DNA sequence variation found within and between species at 14 different genes. D. sechellia consistently revealed low levels of polymorphism, and genes from D. sechellia have accumulated mutations at a rate that is approximately 50% higher than the same genes from D. simulans. At synonymous sites, D. sechellia has experienced a significant excess of unpreferred codon substitutions. Together these observations suggest that D. sechellia has had a reduced effective population size for some time, and that it is accumulating slightly deleterious mutations as a result. D. simulans and D. mauritiana are both highly polymorphic and the two species share many polymorphisms, probably since the time of common ancestry. A simple isolation speciation model, with zero gene flow following incipient species separation, was fitted to both the simulans/mauritiana divergence and the simulans/sechellia divergence. In both cases the model fit the data quite well, and the analyses revealed little evidence of gene flow between the species. The exception is one gene copy at one locus in D. sechellia, which closely resembled other D. simulans sequences. The overall picture is of two allopatric speciation events that occurred quite near one another in time.

Animals↗

Identification of a novel HLA-DPB1 allele, DPB1*9701, by sequence-based typing.

This report describes the identification of a novel DPB1 allele, DPB *9701, found in an Italian Caucasian individual. The new allele was detected by human leukocyte antigen sequence-based typing carried out to investigate the role of genetic factors in determining the outcome of hepatitis C virus infection. DPB1*9701 was identical to DPB1*0501 except for a single-nucleotide substitution at codon 43 (GGG --> TGG). This nucleotide change is a non-synonymous mutation and results in the amino acid substitution glycine (G) --> tryptophan (W). The nucleotide sequence has been deposited in GenBank under the accession number AY033075, and denominated DPB1*9701 by the official World Health Organization Nomenclature Committee.

Base Sequence↗

Compartmentalization of hepatitis C virus quasispecies in blood mononuclear cells of patients with mixed cryoglobulinemic syndrome.

The aim of this study was to investigate the quasispecies heterogeneity of hepatitis C virus (HCV) in the plasma, cryoprecipitate, and peripheral lymphocytes of chronically infected HCV patients with mixed cryoglobulinemia (MC). We studied 360 clones from 10 HCV-positive patients with MC and 8 age-, gender- and HCV genotype-matched subjects with chronic HCV infection but without MC. A partial nucleotide sequence encompassing the E1/E2 region, including hypervariable region 1 (HVR1), was amplified and cloned from plasma, cryoprecipitates, and peripheral blood mononuclear cells (PBMC), and the genetic diversity and complexity and synonymous and nonsynonymous substitution rates were determined. Heterogeneous selection pressure at codon sites was evaluated. Compartmentalization was estimated by phylogenetic and phenetic (Mantel's test) approaches. The patients with MC had 3.3 times lower nonsynonymous substitution rates (1.7 versus 5.7 substitutions/100 sites). Among the subjects with HCV genotype 1, the MC patients had significantly less complexity than the controls, whereas the diversity and complexity were similar in the genotype 2 patients and controls. Site-specific selection analysis confirmed the low frequency of MC patients showing positive selection. There was a significant correlation between positive selection and the infecting HCV genotype. The quasispecies were less heterogeneous in PBMC than in plasma. Significant compartmentalization of HCV quasispecies was observed in the PBMC of four of nine subjects (three with MC) and seven of nine cryoprecipitates. In one subject with MC, we detected a 5-amino-acid insertion at codons 385 to 389 of HVR1. Our results suggest reduced quasispecies heterogeneity in MC patients that is related to a low selection pressure which is probably due to an impaired immune response, the HCV genotype, and/or the duration of the infection. The frequent HCV quasispecies compartmentalization in patients' PBMC suggests a possible pathogenetic significance.

Aged↗

The 'evolutionary signal' of homoplasy in protein-coding gene sequences and its consequences for a priori weighting in phylogeny.

To analyse independently homoplasy for the six possible types of substitution (i.e., A-G, C-T, A-C, A-T, C-G and G-T) at each of the three codon-positions of the cytochrome b gene, two approaches were used: the first is based on the consistency index which measures the amount of homoplasy, and the second is based on the saturation analysis which describes graphically the distribution of homoplasy within the taxonomic sampling. The results obtained from a data set of 32 sequences of Artiodactyla indicate that evolution of the cytochrome b is governed by differential constraints: 1) between the six substitutions-types, 2) between the three codon-positions, and 3) between the two mtDNA strands. Moreover, we find that non-synonymous sites can be more homoplastic than synonymous sites when the possibilities of substitutions are severely restricted because of the functional requirements of hydrophobicity. Most weighting schemes applied to protein-coding genes are elaborated from unjustified assumptions. We propose to weight each substitution-type at each codon-position according to its homoplasy content evaluated either with the consistency index or with an index representing the level of mutational saturation.

Amino Acid Sequence↗

Hepadnavirus evolution and molecular strategy of adaptation in a new host.

In order to elucidate the mechanisms of hepadnavirus evolution in vivo and to trace the fate of known quasispecies in a single animal during the acute phase of infection, a woodchuck (Marmota monax) was infected with the hepadnavirus woodchuck hepatitis B virus (WHV). Woodchuck 197 (W197) was injected intravenously with pooled sera collected from a chronic carrier that had been infected originally with a molecular clone of known genome sequence (WHV7). Viral genome variants from both the inoculum and the follow-up sera from W197 were characterized for the presence of quasispecies related to the WHV7 sequence. Interestingly, WHV7-related genomes were predominant 6 weeks post-infection (p.i.), whereas a highly heterogeneous virus population was present in the first viraemic serum (4 weeks p.i.). Using WHV7 as the prototype, the variability of the Pol and PreS/S regions in the first 11 weeks p.i. has been calculated. The sequence population in serum collected 6 weeks p.i. was highly homogeneous, with a mean variability of 0.36% in the region analysed. Mean variability values ranging from 0.82% to 1.61% were found in quasispecies from the other sera. The presence of possible selective pressure was analysed by means of the non-synonymous versus synonymous variation ratio (dn/d5). We found that the dn/d5 values were stable for the S ORF (ranging from 2.6 to 3.0), whereas a wider range was observed for the Pol ORF (from 1.4 to 3.0). Furthermore, from the analysis of the variability of the codon positions for the two overlapping ORFs it was found that, in most cases, non-synonymous mutations at position 1 of the Pol ORF (position 3 of the S ORF) corresponded to synonymous variation in the S (Pol) ORF, indicating independent evolution of the encoded proteins.

Acute Disease↗

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↗

On the varied pattern of evolution of 2 fungal genomes: a critique of Hughes and Friedman.

A number of statistical tests have been proposed to detect positive Darwinian selection affecting a few amino acid sites in a protein, exemplified by an excess of nonsynonymous nucleotide substitutions. These tests are often more powerful than pairwise sequence comparison, which averages synonymous (d(S)) and nonsynonymous (d(N)) rates over the whole gene. In a recent study, however, Hughes AL and Friedman R (2005. Variation in the pattern of synonymous and nonsynonymous difference between two fungal genomes. Mol Bio Evol. 22: 1320-1324) argue that d(S) and d(N) are expected to fluctuate along the sequence by chance and that an excess of nonsynonymous differences in individual codons is no evidence for positive selection. The authors compared codons in protein-coding genes from the genomes of 2 yeast species, Saccharomyces cerevisiae and Saccharomyces paradoxus. They calculated the proportions of synonymous and nonsynonymous differences per site (p(S) and p(N)) in every codon and discovered that p(N) is often greater than p(S) and that among some codons p(S) and p(N) are negatively correlated. The authors argued that these results invalidate previous tests of codons under positive selection. Here I discuss several errors of statistics in the analysis of Hughes and Friedman, including confusion of statistics with parameters, arbitrary data filtering, and derivation of hypotheses from data. I also apply likelihood ratio tests of positive selection to the yeast data and illustrate empirically that Hughes and Friedman's criticisms on such tests are not valid.

Bias↗

Codon-substitution models for heterogeneous selection pressure at amino acid sites.

Comparison of relative fixation rates of synonymous (silent) and nonsynonymous (amino acid-altering) mutations provides a means for understanding the mechanisms of molecular sequence evolution. The nonsynonymous/synonymous rate ratio (omega = d(N)d(S)) is an important indicator of selective pressure at the protein level, with omega = 1 meaning neutral mutations, omega < 1 purifying selection, and omega > 1 diversifying positive selection. Amino acid sites in a protein are expected to be under different selective pressures and have different underlying omega ratios. We develop models that account for heterogeneous omega ratios among amino acid sites and apply them to phylogenetic analyses of protein-coding DNA sequences. These models are useful for testing for adaptive molecular evolution and identifying amino acid sites under diversifying selection. Ten data sets of genes from nuclear, mitochondrial, and viral genomes are analyzed to estimate the distributions of omega among sites. In all data sets analyzed, the selective pressure indicated by the omega ratio is found to be highly heterogeneous among sites. Previously unsuspected Darwinian selection is detected in several genes in which the average omega ratio across sites is <1, but in which some sites are clearly under diversifying selection with omega > 1. Genes undergoing positive selection include the beta-globin gene from vertebrates, mitochondrial protein-coding genes from hominoids, the hemagglutinin (HA) gene from human influenza virus A, and HIV-1 env, vif, and pol genes. Tests for the presence of positively selected sites and their subsequent identification appear quite robust to the specific distributional form assumed for omega and can be achieved using any of several models we implement. However, we encountered difficulties in estimating the precise distribution of omega among sites from real data sets.

Amino Acid Substitution↗

Rates of synonymous substitution and base composition of nuclear genes in Drosophila.

We compared the rates of synonymous (silent) substitution among various genes in a number of species of Drosophila. First, we found that even for a particular gene, the rate of synonymous substitution varied considerably with Drosophila lineages. Second, we showed a large variation in synonymous substitution rates among nuclear genes in Drosophila. These rates of synonymous substitution were correlated negatively with C content and positively with A content at the third codon positions. Nucleotide sequences were also compared between pseudogenes and their functional homologs. The C content of the pseudogenes was lower than that of the functional genes and the A content of the former was higher than that of the latter. Because the synonymous substitution for functional genes and the nucleotide substitution for pseudogenes are exempted from any selective constraint at the protein level, these observations could be explained by a biased pattern of mutation in the Drosophila nuclear genome. Such a bias in the mutation pattern may affect the molecular clock (local clock) of each nuclear gene of each species. Finally, we obtained the average rates of synonymous substitution for three gene groups in Drosophila; 11.0 x 10(-9), 17.5 x 10(-9) and 27.1 x 10(-9)/site/year.

Animals↗

Molecular evolution of dinoflagellate luciferases, enzymes with three catalytic domains in a single polypeptide.

Enzymes with multiple catalytic sites are rare, and their evolutionary significance remains to be established. This study of luciferases from seven dinoflagellate species examines the previously undescribed evolution of such proteins. All these enzymes have the same unique structure: three homologous domains, each with catalytic activity, preceded by an N-terminal region of unknown function. Both pairwise comparison and phylogenetic inference indicate that the similarity of the corresponding individual domains between species is greater than that between the three different domains of each polypeptide. Trees constructed from each of the three individual domains are congruent with the tree of the full-length coding sequence. Luciferase and ribosomal DNA trees both indicate that the Lingulodinium polyedrum luciferase diverged early from the other six. In all species, the amino acid sequence in the central regions of the three domains is strongly conserved, suggesting it as the catalytic site. Synonymous substitution rates also are greatly reduced in the central regions of two species but not in the other five. This lineage-specific difference in synonymous substitution rates in the central region of the domains correlates inversely with the content of GC3, which can be accounted for by the biased usage toward C-ending codons at the degenerate sites. RNA modeling of the central region of the L. polyedrum luciferase domain suggests a function of the constrained synonymous substitutions in the circadian-controlled protein synthesis.

Amino Acid Sequence↗

Dynamic structures and functions of transfer ribonucleic acids from extreme thermophiles.

tRNA species from an extreme thermophile T. thermophilus that grows up to 85 degrees C have been found to be more thermostable than those from moderate thermophiles and mesophiles. Such thermostability of T. thermophilus tRNA species is partly due to the high contents of G.C base pairs in the stem regions. In addition, a novel modified nucleoside s2T has been found that substitutes T in position 54. The extent of 2-thiolation of T(54) has been found to depend on environmental temperatures from 50 to 80 degrees C. Two tRNA(Ile) species have been isolated from T. thermophilus HB8, tRNA(1aIle) with s2T(54) and tRNA(1bIle) with T(54), which have the identical nucleotide sequence except for position 54. However, the melting temperature of tRNA(1aIle) is higher by 3 degrees C than that of tRNA(1bIle). This clearly indicates that the 2-thiolation of T(54) contributes directly to the thermostability of T. thermophilus tRNA species. Proton NMR analyses have shown that the nucleoside s2T is "rigid" and predominantly takes the C3'-endo-gg-anti form of A-RNA, because of the steric effect of the bulky 2-thiocarbonyl groups and the 2'-hydroxyl group. Thus, the inherent rigidity of s2T in position 54 significantly enhances the stability of the tertiary structure of tRNA. In protein synthesis of T. thermophilus, s2T(54)-bearing tRNA and T(54)-bearing tRNA species are selectively utilized depending on environmental temperature. In the anticodons of major tRNA species from T. thermophilus, G or C exclusively appears in the first position, and GGN and CCN are favored over synonymous GCN or CGN. These characteristic anticodon sequences correspond to the characteristic codon usage in thermophile genes.

Anticodon↗

Adaptive evolution in the SRZ chemoreceptor families of Caenorhabditis elegans and Caenorhabditis briggsae.

We investigated the possibility of positive selection acting on members of the putative seven-pass chemoreceptor superfamily in Caenorhabditis elegans, which comprises approximately 1,300 genes encoding seven-pass G protein-coupled receptors (GPCRs). Using a maximum-likelihood approach, we conducted statistical tests for evidence of codon sites where the ratio of nonsynonymous substitutions per nonsynonymous site to synonymous substitutions per synonymous site (d(N)/d(S)) was >1. Evidence for positive selection was found only for the srz family, about which virtually nothing specific is known. We extended the annotation of the srz gene family, establishing gene models for 60 srz genes in C. elegans and 28 srz genes in Caenorhabditis briggsae. d(N)/d(S) ratios varied dramatically in different regions of the SRZ proteins, peaking in predicted extracellular regions. These regions included 23 sites where evidence of positive selection was highly significant, corresponding remarkably well with regions implicated in ligand binding in other GPCR family members. We interpret these results as indicating that the srz family is under positive selection, probably driven by ligand binding.

Adaptation, Physiological↗

The atypical codon usage of the plant psbA gene may be the remnant of an ancestral bias.

The psbA gene of the chloroplast genome has a codon usage that is unusual for plant chloroplast genes. In the present study the evolutionary status of this codon usage is tested by reconstructing putative ancestral psbA sequences to determine the pattern of change in codon bias during angiosperm divergence. It is shown that the codon biases of the ancestral genes are much stronger than all extant flowering plant psbA genes. This is related to previous work that demonstrated a significant increase in synonymous substitution in psbA relative to other chloroplast genes. It is suggested, based on the two lines of evidence, that the codon bias of this gene currently is not being maintained by selection. Rather, the atypical codon bias simply may be a remnant of an ancestral codon bias that now is being degraded by the mutation bias of the chloroplast genome, in other words, that the psbA gene is not at equilibrium. A model for the evolution of selective pressure on the codon usage of plant chloroplast genes is discussed.

Base Sequence↗