CGN Codons Research Articles

Genomic factors shaping codon usage across the Saccharomycotina subphylum.

Codon usage bias, or the unequal use of synonymous codons, is observed across genes, genomes, and between species. It has been implicated in many cellular functions, such as translation dynamics and transcript stability, but can also be shaped by neutral forces. We characterized codon usage across 1,154 strains from 1,051 species from the fungal subphylum Saccharomycotina to gain insight into the biases, molecular mechanisms, evolution, and genomic features contributing to codon usage patterns. We found a general preference for A/T-ending codons and correlations between codon usage bias, GC content, and tRNA-ome size. Codon usage bias is distinct between the 12 orders to such a degree that yeasts can be classified with an accuracy greater than 90% using a machine-learning algorithm. We also characterized the degree to which codon usage bias is impacted by translational selection. We found it was influenced by a combination of features, including the number of coding sequences, BUSCO count, and genome length. Our analysis also revealed an extreme bias in codon usage in the Saccharomycodales associated with a lack of predicted arginine tRNAs that decode CGN codons, leaving only the AGN codons to encode arginine. Analysis of Saccharomycodales gene expression, tRNA sequences, and codon evolution suggests that avoidance of the CGN codons is associated with a decline in arginine tRNA function. Consistent with previous findings, codon usage bias within the Saccharomycotina is shaped by genomic features and GC bias. However, we find cases of extreme codon usage preference and avoidance along yeast lineages, suggesting additional forces may be shaping the evolution of specific codons.

Arginines of the CGN codon family are Achilles' heels of cancer genes.

Recent studies have revealed that arginine is the most favorable target of amino acid alteration in most cancer types and it has been suggested that the high preference for arginine mutations reflects the critical roles of this amino acid in the function of proteins. High rates of mutations of arginine residues in cancer, however, might also be due to increased mutability of arginine codons of the CGN family as the CpG dinucleotides of these codons may be methylated. In the present work we have analyzed spectra of single base substitutions of cancer genes (oncogenes, tumor suppressor genes) and passenger genes in cancer tissues to assess the contributions of CpG hypermutability and selection to arginine mutations. Our studies have shown that arginines encoded by the CGN codon family display higher rates of mutation in both cancer genes and passenger genes than arginine codons AGA and AGG that are devoid of CpG dinucleotide, suggesting that the predominance of arginine mutations in cancer is primarily due to CpG hypermutability, rather than selection for arginine replacement. Nevertheless, our results also suggest that CGN codons for arginines may serve as Achilles' heels of cancer genes. CpG hypermutability of key arginines of proto-oncogenes, leading to high rates of recurrence of driver mutations, contributes significantly to carcinogenesis. Similarly, our results indicate that hypermutability of the CpG dinucleotide of CGA codons (converting them to TGA stop codons) contributes significantly to recurrent truncation and inactivation of tumor suppressor genes.

YaaJ, the tRNA-Specific Adenosine Deaminase, Is Dispensable in Bacillus subtilis.

Post-transcriptional modifications of tRNA are crucial for their core function. The inosine (I; 6-deaminated adenosine) at the first position in the anticodon of tRNAArg(ICG) modulates the decoding capability and is generally considered essential for reading CGU, CGC, and CGA codons in eubacteria. We report here that the Bacillus subtilis yaaJ gene encodes tRNA-specific adenosine deaminase and is non-essential for viability. A β-galactosidase reporter assay revealed that the translational activity of CGN codons was not impaired in the yaaJ-deletion mutant. Furthermore, tRNAArg(CCG) responsible for decoding the CGG codon was dispensable, even in the presence or absence of yaaJ. These results strongly suggest that tRNAArg with either the anticodon ICG or ACG has an intrinsic ability to recognize all four CGN codons, providing a fundamental concept of non-canonical wobbling mediated by adenosine and inosine nucleotides in the anticodon. This is the first example of the four-way wobbling by inosine nucleotide in bacterial cells. On the other hand, the absence of inosine modification induced +1 frameshifting, especially at the CGA codon. Additionally, the yaaJ deletion affected growth and competency. Therefore, the inosine modification is beneficial for translational fidelity and proper growth-phase control, and that is why yaaJ has been actually conserved in B. subtilis.

Biases in arginine codon usage correlate with genetic disease risk

PurposeThe persistence of hypermutable CGN (CGG, CGA, CGC, CGU) arginine codons at high frequency suggests the possibility of negative selective pressure at these sites and that arginine codon usage could be a predictive indicator of human disease genes. MethodsWe analyzed arginine codons (CGN, AGG, AGA) from all canonical Ensembl protein coding gene transcripts before comparing the frequency of CGN codons between genes with and without human disease associations and with gnomAD constraint metrics. ResultsThe frequency of CGN codons among a gene’s total arginine codon count was higher in genes linked to syndromic autism spectrum disorder (ASD) compared with genes not associated with ASD. A comparison of genes annotated as dominant or recessive with control genes not matching either classification revealed a progressive increase in CGN codon frequency. Moreover, CGN frequency was positively correlated with a gene’s probability of loss-of-function intolerance (pLI) score and negatively correlated with observed-over-expected ratios for both loss-of-function and missense variants. ConclusionOur findings indicate that genes utilizing CGN arginine codons rather than AGG or AGA are more likely to underlie single-gene disorders, particularly for dominant phenotypes, and thus constitute candidate genes for the study of human genetic disease.

Complete mitochondrial genome and evolutionary analysis of Turritopsis dohrnii, the “immortal” jellyfish with a reversible life-cycle

Turritopsis dohrnii (Cnidaria, Hydrozoa, Hydroidolina, Anthoathecata) is the only known metazoan that is capable of reversing its life cycle via morph rejuvenation from the adult medusa stage to the juvenile polyp stage. Here, we present a complete mitochondrial (mt) genome sequence of T. dohrnii, which harbors genes for 13 proteins, two transfer RNAs, and two ribosomal RNAs. The T. dohrnii mt genome is characterized by typical features of species in the Hydroidolina subclass, such as a high A+T content (71.5%), reversed transcriptional orientation for the large rRNA subunit gene, and paucity of CGN codons. An incomplete complementary duplicate of the cox1 gene was found at the 5′ end of the T. dohrnii mt chromosome, as were variable repeat regions flanking the chromosome. We identified species-specific variations (nad5, nad6, cob, and cox1 genes) and putative selective constraints (atp8, nad1, nad2, and nad5 genes) in the mt genes of T. dohrnii, and predicted alterations in tertiary structures of respiratory chain proteins (NADH4, NADH5, and COX1 proteins) of T. dohrnii. Based on comparative analyses of available hydrozoan mt genomes, we also determined the taxonomic relationships of T. dohrnii, recovering Filifera IV as a paraphyletic taxon, and assessed intraspecific diversity of various Hydrozoa species.

Where have all the inosines gone? Conflicting evidence for A-to-I editing of the anticodon of higher eukaryotic tRNAACGArg questions the dogma of a universal wobble-mediated decoding of CGN codons.

Codon-anticodon recognition between triplets of an mRNA and a specific tRNA is the key element in the translation of the genetic code. In general, the precision of this process is dominated by a strict Watson-Crick base-pairing scheme. However, the degeneracy of the genetic code led Crick to propose the Wobble Hypothesis, permitting a less restraining interaction with the third base of the codon and involving the participation of inosine for decoding C-ending codons. The concept that the anticodon base A34 of tRNAACGArg in all eukaryotes, eubacteria, and plant chloroplasts is converted to I34 is firmly anchored in the literature despite conflicting evidence for its existence in higher eukaryote cytoplasmic tRNAACGArg. Here, we provide additional data and summarize the arguments favoring and contradicting post-transcriptional deamination of this position. A hypothesis that resolves the apparent conflict is proposed. © 2016 IUBMB Life, 68(6):419-422, 2016.

Life without tRNAArg–adenosine deaminase TadA: evolutionary consequences of decoding the four CGN codons as arginine in Mycoplasmas and other Mollicutes

In most bacteria, two tRNAs decode the four arginine CGN codons. One tRNA harboring a wobble inosine (tRNAArgICG) reads the CGU, CGC and CGA codons, whereas a second tRNA harboring a wobble cytidine (tRNAArgCCG) reads the remaining CGG codon. The reduced genomes of Mycoplasmas and other Mollicutes lack the gene encoding tRNAArgCCG. This raises the question of how these organisms decode CGG codons. Examination of 36 Mollicute genomes for genes encoding tRNAArg and the TadA enzyme, responsible for wobble inosine formation, suggested an evolutionary scenario where tadA gene mutations first occurred. This allowed the temporary accumulation of non-deaminated tRNAArgACG, capable of reading all CGN codons. This hypothesis was verified in Mycoplasma capricolum, which contains a small fraction of tRNAArgACG with a non-deaminated wobble adenosine. Subsets of Mollicutes continued to evolve by losing both the mutated tRNAArgCCG and tadA, and then acquired a new tRNAArgUCG. This permitted further tRNAArgACG mutations with tRNAArgGCG or its disappearance, leaving a single tRNAArgUCG to decode the four CGN codons. The key point of our model is that the A-to-I deamination activity had to be controlled before the loss of the tadA gene, allowing the stepwise evolution of Mollicutes toward an alternative decoding strategy.

Mitochondrial DNA of Clathrina clathrus (Calcarea, Calcinea): Six Linear Chromosomes, Fragmented rRNAs, tRNA Editing, and a Novel Genetic Code

Sponges (phylum Porifera) are a large and ancient group of morphologically simple but ecologically important aquatic animals. Although their body plan and lifestyle are relatively uniform, sponges show extensive molecular and genetic diversity. In particular, mitochondrial genomes from three of the four previously studied classes of Porifera (Demospongiae, Hexactinellida, and Homoscleromorpha) have distinct gene contents, genome organizations, and evolutionary rates. Here, we report the mitochondrial genome of Clathrina clathrus (Calcinea, Clathrinidae), a representative of the fourth poriferan class, the Calcarea, which proves to be the most unusual. Clathrina clathrus mitochondrial DNA (mtDNA) consists of six linear chromosomes 7.6-9.4 kb in size and encodes at least 37 genes: 13 protein codings, 2 ribosomal RNAs (rRNAs), and 24 transfer RNAs (tRNAs). Protein genes include atp9, which has now been found in all major sponge lineages, but no atp8. Our analyses further reveal the presence of a novel genetic code that involves unique reassignments of the UAG codons from termination to tyrosine and of the CGN codons from arginine to glycine. Clathrina clathrus mitochondrial rRNAs are encoded in three (srRNA) and ≥6 (lrRNA) fragments distributed out of order and on several chromosomes. The encoded tRNAs contain multiple mismatches in the aminoacyl acceptor stems that are repaired posttranscriptionally by 3'-end RNA editing. Although our analysis does not resolve the phylogenetic position of calcareous sponges, likely due to their high rates of mitochondrial sequence evolution, it confirms mtDNA as a promising marker for population studies in this group. The combination of unusual mitochondrial features in C. clathrus redefines the extremes of mtDNA evolution in animals and further argues against the idea of a "typical animal mtDNA."

Efficient expression of gene variants that harbour AGA codons next to the initiation codon

In an effort to improve the knowledge about the rules which direct the effect of the early ORF sequences on translation efficiency, we have analyzed the effect of pairs of the six arginine codons at the second and third positions on the expression of lacZ variants. Whereas the pairs of identical AGA or AGG codons were favorable for the gene expression, identical pairs of each of the four CGN codons were very inefficient. This result was unexpected because tandems of AGA or AGG codons located in more internal gene positions provoke deficient expression whilst internally located CGU and CGC are the most abundant and efficiently translated arginine codons. The mixed combinations of AGA and each of the CGN codons usually resulted in efficient rates of lacZ expression independently of the peptidyl-tRNA propensity to dissociate from the ribosome. Thus, the variant harboring the pair of AGA codons was expressed as efficiently as the variant carrying a pair of AAA codons in the same positions, a configuration reported as one of the most common and efficient for gene expression. We explain these results assuming that the presence of adenines in these early positions enhance gene expression. As expected, specific mRNA levels correlated with the intensity of lacZ expression for each variant. However, the induction of lacZ AGA AGA gene in pth cells accumulated peptidyl-tRNAArg4 as well as a short 5′-proximal lacZ mRNA fragment suggesting ribosome stalling due to depletion of aminoacylated-tRNAArg4.

CpG islands are the second main factor shaping codon usage in human genes

A correspondence analysis of codon usage in human genes revealed, as expected, that the first axis is strongly correlated with the base composition at synonymous third codon positions. At one extreme of the second axis were localized genes with a high frequency of NCG and CGN codons. The great majority of these sequences were embedded in CpG islands, while the opposite is true for the genes placed at the other extreme. The two main conclusions of this paper are: (1) the influence of CpG islands on codon usage, and (2) since the second axis is orthogonal (and therefore independent) of the first, GC3-rich genes are not necessarily associated with CpG islands.

The RNA polymerase III-dependent family of genes in hemiascomycetes: comparative RNomics, decoding strategies, transcription and evolutionary implications.

We present the first comprehensive analysis of RNA polymerase III (Pol III) transcribed genes in ten yeast genomes. This set includes all tRNA genes (tDNA) and genes coding for SNR6 (U6), SNR52, SCR1 and RPR1 RNA in the nine hemiascomycetes Saccharomyces cerevisiae, Saccharomyces castellii, Candida glabrata, Kluyveromyces waltii, Kluyveromyces lactis, Eremothecium gossypii, Debaryomyces hansenii, Candida albicans, Yarrowia lipolytica and the archiascomycete Schizosaccharomyces pombe. We systematically analysed sequence specificities of tRNA genes, polymorphism, variability of introns, gene redundancy and gene clustering. Analysis of decoding strategies showed that yeasts close to S.cerevisiae use bacterial decoding rules to read the Leu CUN and Arg CGN codons, in contrast to all other known Eukaryotes. In D.hansenii and C.albicans, we identified a novel tDNA-Leu (AAG), reading the Leu CUU/CUC/CUA codons with an unusual G at position 32. A systematic ‘p-distance tree’ using the 60 variable positions of the tRNA molecule revealed that most tDNAs cluster into amino acid-specific sub-trees, suggesting that, within hemiascomycetes, orthologous tDNAs are more closely related than paralogs. We finally determined the bipartite A- and B-box sequences recognized by TFIIIC. These minimal sequences are nearly conserved throughout hemiascomycetes and were satisfactorily retrieved at appropriate locations in other Pol III genes.

Complete nucleotide sequence of the mitochondrial DNA from

The total nucleotide sequence of the mitochondrial genome of the yeast Kluyveromyces lactis was determined. The DNA is a circular molecule of 40,291 base pairs, with 26.1% GC. It contains a set of protein- and RNA-coding genes equivalent to those of the Saccharomyces cerevisiae mitochondrial genome. The genome size is about one half of that of S. cerevisiae mitochondrial DNA. The difference in size is due essentially to a reduced proportion of intergenic and intronic sequences. The coding sequences occupy about one third of the genome, the rest being composed of AT-rich sequences and numerous short GC-rich clusters that are dispersed mostly in the non-coding regions and a few within coding sequences. The presence of these GC clusters is a characteristic feature common to K. lactis and S. cerevisiae mitochondrial DNA, although their sequence patterns are different. The absence of the NADH dehydrogenase subunit genes distinguishes this yeast and S. cerevisiae from the typically aerobic species. The genetic code appears to be that of the standard fungal mitochondrial genomes, with UGA as a tryptophan codon. There are only 22 transfer RNA genes, those corresponding to CUN and CGN codons being missing. CUN codons are absent in the protein-coding sequences. There are five CGN codons within the open reading frames, but they are located exclusively in the introns, rendering them untranslatable. Introns are found only the genes in KlCOX1 and LrRNA. The transcription promoter motif known in S. cerevisiae and several other yeast species is also present. All genes are transcribed from the same strand, except those on a single 7-kilobase pairs segment (EMBL Accession No. AY654900).

The conserved Cys-X1-X2-Cys motif present in the TtcA protein is required for the thiolation of cytidine in position 32 of tRNA from Salmonella enterica serovar Typhimurium.

The modified nucleoside 2-thiocytidine (s(2)C) has so far been found in tRNA from organisms belonging to the phylogenetic domains Archaea and Bacteria. In the bacteria Escherichia coli and Salmonella enterica serovar Typhimurium, s(2)C is present in position 32 of only four tRNA species-, and. An in-frame deletion of an S. enterica gene (designated ttcA, for "two-thio-cytidine") was constructed, and such a mutant has no detectable s(2)C in its tRNA. The TtcA protein family is characterized by the existence of both a PP-loop and a Cys-X(1)-X(2)-Cys motif in the central region of the protein but can be divided into two distinct groups based on the presence and location of additional Cys-X(1)-X(2)-Cys motifs in terminal regions of the sequence. Mutant analysis showed that both cysteines in this central conserved Cys-X(1)-X(2)-Cys motif are required for the formation of s(2)C. The DeltattcA1 mutant grows at the same rate as the congenic wild-type strain, and no growth disadvantage caused by the lack of s(2)C was observed in a mixed-population experiment. Lack of s(2)C32 did not reduce the selection rate at the ribosomal aminoacyl-tRNA site (A-site) for at any of its cognate CGN codons, whereas A-site selection at AGG by was dependent on the presence of s(2)C32. The presence of s(2)C32 in peptidyl- or in peptidyl- interfered with decoding in the A-site. The presence of s(2)C32 in decreased the rate of translation of the CGA codon but not that of the CGU codon.

The complete mitochondrial genome of Yarrowia lipolytica

We here report the complete nucleotide sequence of the 47.9 kb mitochondrial (mt) genome from the obligate aerobic yeast Yarrowia lipolytica. It encodes, all on the same strand, seven subunits of NADH: ubiquinone oxidoreductase (ND1-6, ND4L), apocytochrome b (COB), three subunits of cytochrome oxidase (COX1, 2, 3), three subunits of ATP synthetase (ATP6, 8 and 9), small and large ribosomal RNAs and an incomplete set of tRNAs. The Y. lipolytica mt genome is very similar to the Hansenula wingei mt genome, as judged from blocks of conserved gene order and from sequence homology. The extra DNA in the Y. lipolytica mt genome consists of 17 group 1 introns and stretches of A+Trich sequence, interspersed with potentially transposable GC clusters. The usual mould mt genetic code is used. Interestingly, there is no tRNA able to read CGN (arginine) codons. CGN codons could not be found in exonic open reading frames, whereas they do occur in intronic open reading frames. However, several of the intronic open reading frames have accumulated mutations and must be regarded as pseudogenes. We propose that this may have been triggered by the presence of untranslatable CGN codons. This sequence is available under EMBL Accession No. AJ307410.

Primary sequence of mitochondrial tRNA Arg of a nematode Ascaris suum: occurrence of unmodified adenosine at the first position of the anticodon

Mitochondrial tRNA Arg from a nematode, Ascaris suum, was purified and sequenced at the RNA level. An unmodified adenosine was found to exist at the anticodon first position, suggesting that, contrary to the conventional wobble rule, the anticodon ACG of the tRNA can translate all the CGN codons.

The genetic code in mitochondria and chloroplasts

The universal genetic code is used without changes in chloroplasts and in mitochondria of green plants. Non-plant mitochondria use codes that include changes from the universal code. Chloroplasts use 31 anticodons in translating the code; a number smaller than that used by bacteria, because chloroplasts have eliminated 10 CNN anticodons that are found in bacteria. Green plant mitochondria (mt) obtain some tRNAs from the cytosol, and genes for some other tRNAs have been acquired from chloroplast DNA. The code in non-plant mt differs from the universal code in the following usages found in various organisms: UGA for Trp, AUA for Met, AGR for Ser and stop, AAA for Asn, CUN for Thr, and possibly UAA for Tyr. CGN codons are not used by Torulopsis yeast mt. Non-plant mt, e.g. in vertebrates, may use a minimum of 22 anticodons for complete translation of mRNA sequences. The following possible causes are regarded as contributing to changes in the non-plant mt: directional mutation pressure, genomic economization, changes in charging specificity of tRNAs, loss of release factor RF2, changes in RF1, changes in anticodons, loss of lysidine-forming enzyme system, and disappearance of codons from coding sequences.

Codon recognition rules in yeast mitochondria.

The mitochondrial genome of Saccharomyces cerevisiae codes for 24 tRNAs. The nucleotide sequences of the tRNA genes suggest a unique set of rules that govern the decoding of the mitochondrial genetic code. The four codons of unmixed fmilies are recognized by single tRNAs that always have a U in the wobble position of the anticodon. The codons of the mixed families are read by two different tRNAs. Codons terminating in a C or U are recognized by tRNAs with a G and codons terminating in a G or A are recognized by tRNAs with a U in the corresponding positions of the anticodons. There are two exceptions to these rules. In the AUN family for isoleucine and methionine, the isoleucine tRNA has a G and the methionine tRNA has a C in the wobble position. The tRNA for the arginine CGN family also has an A in the wobble position of the anticodon. It is of interest that the CGN codons have not been found in the mitochondrial genes sequenced to date. The simplified decoding system of yeast mitochondria allows all the codons to be recognized by only 24 tRNAs.