Evolutionary changes in the genetic code

Complete mitochondrial genome of Rhaphidophora duxiu (Orthoptera: Rhaphidophoridae: Rhaphidophorinae)

We present the complete mitochondrial genome of Rhaphidophora duxiu from China. The mitogenome of R. duxiu is circular, AT-rich (75.3%), and 15,898 bp in length. It comprises 13 protein-coding genes, 2 ribosomal RNA genes, and 22 transfer RNA genes. It is identical in gene content to Rhaphidophora quadrispina.

Virulence effectors encoded in the rice yellow dwarf phytoplasma genome participate in pathogenesis

Bacteria-like phytoplasmas alternate between plant and insect hosts, secreting proteins that disrupt host development. In this study, we sequenced the complete genome of "Candidatus Phytoplasma oryzae" strain HN2022, associated with rice yellow dwarf (RYD) disease, using PacBio HiFi technology. The strain was classified within the 16Sr XI-B subgroup. Through SignalP v5.0 for prediction and subsequent expression analysis of secreted proteins in Nicotiana benthamiana and rice (Oryza sativa L.), we identified the key virulence effector proteins RY348 and RY378. RY348, a homolog of Secreted Aster Yellows Phytoplasma Effector 54 (SAP54), targets and degrades the MADS-box transcription factors MADS1 and MADS15, causing pollen sterility. Meanwhile, RY378 impacts the strigolactone and auxin signaling pathways, substantially increasing tillering. These findings offer insights into the interactions between plants and phytoplasmas.

Complete mitochondrial genome of the introduced Indian walking stick Carausius morosus (Lonchodidae, Insecta) from California

We present the complete mitochondrial genome of Carausius morosus from Salinas, CA. The mitochondrial genome of C. morosus is circular, AT rich (78.1%), and 16,671 bp in length. It consists of 13 protein-coding, 22 transfer RNA, and 2 ribosomal RNA genes and is identical in gene content to Carausius sp.

Rare-event sampling analysis uncovers the fitness landscape of the genetic code

The genetic code refers to a rule that maps 64 codons to 20 amino acids. Nearly all organisms, with few exceptions, share the same genetic code, the standard genetic code (SGC). While it remains unclear why this universal code has arisen and been maintained during evolution, it may have been preserved under selection pressure. Theoretical studies comparing the SGC and numerically created hypothetical random genetic codes have suggested that the SGC has been subject to strong selection pressure for being robust against translation errors. However, these prior studies have searched for random genetic codes in only a small subspace of the possible code space due to limitations in computation time. Thus, how the genetic code has evolved, and the characteristics of the genetic code fitness landscape, remain unclear. By applying multicanonical Monte Carlo, an efficient rare-event sampling method, we efficiently sampled random codes from a much broader random ensemble of genetic codes than in previous studies, estimating that only one out of every 1020 random codes is more robust than the SGC. This estimate is significantly smaller than the previous estimate, one in a million. We also characterized the fitness landscape of the genetic code that has four major fitness peaks, one of which includes the SGC. Furthermore, genetic algorithm analysis revealed that evolution under such a multi-peaked fitness landscape could be strongly biased toward a narrow peak, in an evolutionary path-dependent manner.

Two-Step Epimerization of Deoxynivalenol by Quinone-Dependent Dehydrogenase and Candida parapsilosis ACCC 20221

Deoxynivalenol (DON), one of the main mycotoxins with enteric toxicity, genetic toxicity, and immunotoxicity, and is widely found in corn, barley, wheat, and rye. In order to achieve effective detoxification of DON, the least toxic 3-epi-DON (1/357th of the toxicity of DON) was chosen as the target for degradation. Quinone-dependent dehydrogenase (QDDH) reported from Devosia train D6-9 detoxifies DON by converting C3-OH to a ketone group with toxicity of less than 1/10 that of DON. In this study, the recombinant plasmid pPIC9K-QDDH was constructed and successfully expressed in Pichia pastoris GS115. Within 12 h, recombinant QDDH converted 78.46% of the 20 μg/mL DON to 3-keto-DON. Candida parapsilosis ACCC 20221 was screened for its activity in reducing 86.59% of 3-keto-DON within 48 h; its main products were identified as 3-epi-DON and DON. In addition, a two-step method was performed for epimerizing DON: 12 h catalysis by recombinant QDDH and 6 h transformation of the C. parapsilosis ACCC 20221 cell catalyst. The production rates of 3-keto-DON and 3-epi-DON were 51.59% and 32.57%, respectively, after manipulation. Through this study, effective detoxification of 84.16% of DON was achieved, with the products being mainly 3-keto-DON and 3-epi-DON.

A crescendo of competent coding (c3) contains the Standard Genetic Code

The Standard Genetic Code (SGC) can arise by fusion of partial codes evolved in different individuals, perhaps for differing prior tasks. Such code fragments can be unified into an SGC after later evolution of accurate third-position Crick wobble. Late wobble advent fills in the coding table, leaving only later development of translational initiation and termination to reach the SGC in separated domains of life. This code fusion mechanism is computationally implemented here. Late Crick wobble after C3 fusion (c3-lCw) is tested for its ability to evolve the SGC. Compared with previously studied isolated coding tables, or with increasing numbers of parallel, but nonfusing codes, c3-lCw reaches the SGC sooner, is successful in a smaller population, and presents accurate and complete codes more frequently. Notably, a long crescendo of SGC-like codes is exposed for selection of superior translation. c3-lCw also effectively suppresses varied disordered assignments, thus converging on a unified code. Such merged codes closely approach the SGC, making its selection plausible. For example: Under routine conditions, ≈1 of 22 c3-lCw environments evolves codes with ≥20 assignments and ≤3 differences from the SGC, notably including codes identical to the Standard Genetic Code.

Frameshift and wild-type proteins are often highly similar because the genetic code and genomes were optimized for frameshift tolerance

Frameshift mutations have been considered of significant importance for the molecular evolution of proteins and their coding genes, while frameshift protein sequences encoded in the alternative reading frames of coding genes have been considered to be meaningless. However, functional frameshifts have been found widely existing. It was puzzling how a frameshift protein kept its structure and functionality while substantial changes occurred in its primary amino-acid sequence. This study shows that the similarities among frameshifts and wild types are higher than random similarities and are determined at different levels. Frameshift substitutions are more conservative than random substitutions in the standard genetic code (SGC). The frameshift substitutions score of SGC ranks in the top 2.0-3.5% of alternative genetic codes, showing that SGC is nearly optimal for frameshift tolerance. In many genes and certain genomes, frameshift-resistant codons and codon pairs appear more frequently than expected, suggesting that frameshift tolerance is achieved through not only the optimality of the genetic code but, more importantly, the further optimization of a specific gene or genome through the usages of codons/codon pairs, which sheds light on the role of frameshift mutations in molecular and genomic evolution.

Enzymatic noncovalent synthesis of peptide assemblies generates multimolecular crowding in cells for biomedical applications

Enzymatic noncovalent synthesis enables the spatiotemporal control of multimolecular crowding in cells, thus offering a unique opportunity for modulating cellular functions. This article introduces some representative enzymes and molecular building blocks for generating peptide assemblies as multimolecular crowding in cells, highlights the relevant biomedical applications, such as anticancer therapy, molecular imaging, trafficking proteins, genetic engineering, artificial intracellular filaments, cell morphogenesis, and antibacterial, and briefly discusses the promises of ENS as a multistep molecular process in biology and medicine.

Complete Mitochondrial Genome of Pseudocaranx dentex (Carangidae, Perciformes) Provides Insight into Phylogenetic and Evolutionary Relationship among Carangidae Family

Pseudocaranx dentex (white trevally) which belongs to the Carangidae family, is an important commercial fishery and aquaculture resource in Asia. However, its evolution and population genetics have received little attention which was limited by the mitogenome information absence. Here, we sequenced and annotated the complete mitochondrial genome of P. dentex which was 16,569 bp in length, containing twenty-two tRNAs (transfer RNAs), thirteen PCGs (protein-coding genes), two rRNAs (ribosomal RNAs), and one non-coding region with conservative gene arrangement. The Ka/Ks ratio analysis among Carangidae fishes indicated the PCGs were suffering purify selection and the values were related to the taxonomic status and further influenced by their living habits. Phylogenetic analysis based on the PCGs sequences of mitogenomes among 36 species presented three major clades in Carangidae. According to the phylogenetic tree, we further analyzed the taxonomic confusion of Carangoides equula which was on the same branch with P. dentex but a different branch with Carangoides spp. We inferred Kaiwarinus equula should be the accepted name and belong to the independent Kaiwarinus genus which was the sister genus of Pseudocaranx. This work provides mitochondrial genetic information and verifies the taxonomic status of P. dentex, and further helps to recognize the phylogenetic relationship and evolutionary history of Carangidae.

Haplotype-resolved diverse human genomes and integrated analysis of structural variation

Long-read and strand-specific sequencing technologies together facilitate the de novo assembly of high-quality haplotype-resolved human genomes without parent-child trio data. We present 64 assembled haplotypes from 32 diverse human genomes. These highly contiguous haplotype assemblies (average minimum contig length needed to cover 50% of the genome: 26 million base pairs) integrate all forms of genetic variation, even across complex loci. We identified 107,590 structural variants (SVs), of which 68% were not discovered with short-read sequencing, and 278 SV hotspots (spanning megabases of gene-rich sequence). We characterized 130 of the most active mobile element source elements and found that 63% of all SVs arise through homology-mediated mechanisms. This resource enables reliable graph-based genotyping from short reads of up to 50,340 SVs, resulting in the identification of 1526 expression quantitative trait loci as well as SV candidates for adaptive selection within the human population.

Crick Wobble and Superwobble in Standard Genetic Code Evolution

Wobble coding is inevitable during evolution of the Standard Genetic Code (SGC). It ultimately splits half of NN U/C/A/G coding boxes with different assignments. Further, it contributes to pervasive SGC order by reinforcing close spacing for identical SGC assignments. But wobble cannot appear too soon, or it will inhibit encoding and more decisively, obstruct evolution of full coding tables. However, these prior results assumed Crick wobble, NN U/C and NN A/G, read by a single adaptor RNA. Superwobble translates NN U/C/A/G codons, using one adaptor RNA with an unmodified 5' anticodon U (appropriate to earliest coding) in modern mitochondria, plastids, and mycoplasma. Assuming the SGC was selected when evolving codes most resembled it, characteristics of the critical selection events can be calculated. For example, continuous superwobble infrequently evolves SGC-like coding tables. So, continuous superwobble is a very improbable origin hypothesis. In contrast, late-arising superwobble shares late Crick wobble's frequent resemblance to SGC order. Thus late superwobble is possible, but yields SGC-like assignments less frequently than late Crick wobble. Ancient coding ambiguity, most simply, arose from Crick wobble alone. This is consistent with SGC assignments to NAN codons.

Seaview Version 5: A Multiplatform Software for Multiple Sequence Alignment, Molecular Phylogenetic Analyses, and Tree Reconciliation

We present Seaview version 5, a multiplatform program to perform multiple alignment and phylogenetic tree building from molecular sequence data. Seaview provides network access to sequence databases, alignment with arbitrary algorithm, parsimony, distance and maximum likelihood tree building with PhyML, and display, printing, and copy-to-clipboard or to SVG files of rooted or unrooted, binary or multifurcating phylogenetic trees. While Seaview is primarily a program providing a graphical user interface to guide the user into performing desired analyses, Seaview possesses also a command-line mode adequate for user-provided scripts. Seaview version 5 introduces the ability to reconcile a gene tree with a reference species tree and use this reconciliation to root and rearrange the gene tree. Seaview is freely available at http://doua.prabi.fr/software/seaview .

Enzymatic Noncovalent Synthesis for Mitochondrial Genetic Engineering of Cancer Cells

Since mitochondria contribute to tumorigenesis and drug resistance in cancer, mitochondrial genetic engineering promises a new direction for cancer therapy. Here, we report the use of the perimitochondrial enzymatic noncovalent synthesis (ENS) of peptides for delivering genes selectively into the mitochondria of cancer cells for mitochondrial genetic engineering. Specifically, the micelles of peptides bind to the voltage-dependent anion channel (VDAC) on mitochondria for the proteolysis by enterokinase (ENTK), generating perimitochondrial nanofibers in cancer cells. This process, facilitating selective delivery of nucleic acid or gene vectors into mitochondria of cancer cells, enables the mitochondrial transgene expression of CRISPR/Cas9, FUNDC1, p53, and fluorescent proteins. Mechanistic investigation indicates that the interaction of the peptide assemblies with the VDAC and mitochondrial membrane potential are necessary for mitochondria targeting. This local enzymatic control of intermolecular noncovalent interactions enables selective mitochondrial genetic engineering, thus providing a strategy for targeting cancer cells.

Prediction and Characterization of Cationic Arginine-Rich Plant Antimicrobial Peptide SM-985 From Teosinte (Zea mays ssp. mexicana)

Antimicrobial peptides (AMPs) are effective against different plant pathogens and newly considered as part of plant defense systems. From prokaryotes to eukaryotes, AMPs can exist in all forms of life. SM-985 is a cationic AMP (CAMP) isolated from the cDNA library of Mexican teosinte (Zea mays ssp. mexicana). A computational prediction server running with different algorithms was used to screen the teosinte cDNA library for AMPs, and the SM-985 peptide was predicted as an AMP with high probability prediction values. SM-985 is an arginine-rich peptide and composed of 21 amino acids (MW: 2671.06 Da). The physicochemical properties of SM-985 are very promising as an AMP, including the net charge (+8), hydrophobicity ratio of 23%, Boman index of 5.19 kcal/mol, and isoelectric point of 12.95. The SM-985 peptide has amphipathic α-helix conformations. The antimicrobial activity of SM-985 was confirmed against six bacterial plant pathogens, and the MIC of SM-985 against Gram-positive indicators was 8 μM, while the MIC of SM-985 against Gram-negative indicators was 4 μM. The SM-985 interacting with the bacterial membrane and this interaction were examined by treatment of the bacterial indicators with FITC-SM-985 peptide, which showed a high binding affinity of SM-985 to the bacterial membrane (whether Gram-positive or Gram-negative). Scanning electron microscopy (SEM) and transmission electron microscopy (TEM) images of the treated bacteria with SM-985 demonstrated cell membrane damage and cell lysis. In vivo antimicrobial activity was examined, and SM-985 prevented leaf spot disease infection caused by Pst DC3000 on Solanum lycopersicum. Moreover, SM-985 showed sensitivity to calcium chloride salt, which is a common feature of CAMPs.

coil: an R package for cytochrome c oxidase I (COI) DNA barcode data cleaning, translation, and error evaluation

Biological conclusions based on DNA barcoding and metabarcoding analyses can be strongly influenced by the methods utilized for data generation and curation, leading to varying levels of success in the separation of biological variation from experimental error. The 5′ region of cytochrome c oxidase subunit I (COI-5P) is the most common barcode gene for animals, with conserved structure and function that allows for biologically informed error identification. Here, we present coil ( https://CRAN.R-project.org/package=coil ), an R package for the pre-processing and frameshift error assessment of COI-5P animal barcode and metabarcode sequence data. The package contains functions for placement of barcodes into a common reading frame, accurate translation of sequences to amino acids, and highlighting insertion and deletion errors. The analysis of 10 000 barcode sequences of varying quality demonstrated how coil can place barcode sequences in reading frame and distinguish sequences containing indel errors from error-free sequences with greater than 97.5% accuracy. Package limitations were tested through the analysis of COI-5P sequences from the plant and fungal kingdoms as well as the analysis of potential contaminants: nuclear mitochondrial pseudogenes and Wolbachia COI-5P sequences. Results demonstrated that coil is a strong technical error identification method but is not reliable for detecting all biological contaminants.

The transcriptome of the rumen ciliate Entodinium caudatum reveals some of its metabolic features

Background

Rumen ciliates play important roles in rumen function by digesting and fermenting feed and shaping the rumen microbiome. However, they remain poorly understood due to the lack of definitive direct evidence without influence by prokaryotes (including symbionts) in co-cultures or the rumen. In this study, we used RNA-Seq to characterize the transcriptome of Entodinium caudatum, the most predominant and representative rumen ciliate species.

Results

Of a large number of transcripts, > 12,000 were annotated to the curated genes in the NR, UniProt, and GO databases. Numerous CAZymes (including lysozyme and chitinase) and peptidases were represented in the transcriptome. This study revealed the ability of E. caudatum to depolymerize starch, hemicellulose, pectin, and the polysaccharides of the bacterial and fungal cell wall, and to degrade proteins. Many signaling pathways, including the ones that have been shown to function in E. caudatum, were represented by many transcripts. The transcriptome also revealed the expression of the genes involved in symbiosis, detoxification of reactive oxygen species, and the electron-transport chain. Overall, the transcriptomic evidence is consistent with some of the previous premises about E. caudatum. However, the identification of specific genes, such as those encoding lysozyme, peptidases, and other enzymes unique to rumen ciliates might be targeted to develop specific and effective inhibitors to improve nitrogen utilization efficiency by controlling the activity and growth of rumen ciliates. The transcriptomic data will also help the assembly and annotation in future genomic sequencing of E. caudatum.

Conclusion

As the first transcriptome of a single species of rumen ciliates ever sequenced, it provides direct evidence for the substrate spectrum, fermentation pathways, ability to respond to various biotic and abiotic stimuli, and other physiological and ecological features of E. caudatum. The presence and expression of the genes involved in the lysis and degradation of microbial cells highlight the dependence of E. caudatum on engulfment of other rumen microbes for its survival and growth. These genes may be explored in future research to develop targeted control of Entodinium species in the rumen. The transcriptome can also facilitate future genomic studies of E. caudatum and other related rumen ciliates.

Dynamic Characteristics and HIV Infection of Men who have Sex with Men from 2011 to 2017 in Yinchuan, Ningxia, China

OBJECTIVES

Ningxia Hui Autonomous Region, an important area for ethnic Hui settlement in Northwest China, is a low HIV prevalence region. However, HIV infection rates among men who have sex with men (MSM) in Ningxia have increased to an alarming level, despite scale-up of control measures in recent years. This study aimed to understand the demographical and sexual behavior dynamics of MSM and to explore the factors associated with HIV infection.

METHODS

Annual cross-sectional surveys were carried out among MSM during 2011~2017 in Yinchuan, the capital city of Ningxia. Information regarding social demographics, sexual behavior and HIV prevention knowledge was collected. Blood samples were taken for HIV, HCV serological and genetic analysis, and syphilis serological analysis. The dynamic trend was analyzed with trend χ2 test and factors associated with HIV infection were identified by multivariate logistic regression analysis.

RESULTS

The study found a decreasing trend for mean age of the MSM population over the study period. MSMs with a college education or higher increased significantly, while the proportions that were in a marriage significantly decreased over the study period. The rate of HIV positive among MSM increased during the study period (p<0.05), however, the rate of recently diagnosed infections decreased from 2012 (p<0.05). Overall, a very high proportion (98%) of MSM had basic knowledge of HIV prevention, however, only approximately 40% of them used condoms consistently during anal sex with male partners. Unprotected anal sex was identified as a risk factor associated with HIV infection, as was syphilis infection. Local residency status and MSM who received intervention and detection services were the factors that decreased HIV infection risk. Sequence analysis identified the HIV-1 CRF55_01B subtype from MSM for the first time in Yinchuan.

CONCLUSION

The reduction of recent HIV diagnoses is an encouraging sign of successful HIV control measures in MSM in Ningxia. The finding that a high proportion of MSM had knowledge of HIV prevention but still conducted unprotected sex highlights the need for further control measures to change unsafe sexual practices among MSM.

A degeneration-reducing criterion for optimal digital mapping of genetic codes

Bioinformatics may seem to be a scientific field processing primarily large string datasets, as nucleotides and amino acids are represented with dedicated characters. On the other hand, many computational tasks that bioinformatics challenges are mathematical problems understandable as operations with digits. In fact, many computational tasks are solved this way in the background. One of the most widely used digital representations is mapping of nucleotides and amino acids with integers 0–3 and 0–20, respectively. The limitation of this mapping occurs when the digital signal of nucleotides has to be translated into a digital signal of amino acids as the genetic code is degenerated. This causes non-monotonies in a mapping function. Although map for reducing this undesirable effect has already been proposed, it is defined theoretically and for standard genetic codes only. In this study, we derived a novel optimal criterion for reducing the influence of degeneration by utilizing a large dataset of real sequences with various genetic codes. As a result, we proposed a new robust global optimal map suitable for any genetic code as well as specialized optimal maps for particular genetic codes.

•

Optimization of 1D numerical representation for DNA to protein translation.

•

Reducing genetic code degeneracy in numerical representation of DNA sequences.

•

More robust numerical conversion used for genomic-proteomic analysis.

Genome Assembly

Genome assembly uses sequence similarity to go from sequencing reads to longer contiguous sequences (contigs). Scaffolds are contigs linked together by gaps where the order and orientation of the contigs is known but the exact sequence connecting two contigs is unknown, represented by Ns which estimate the gap length. Here we describe recommendations for genome assembly for different sequencing technologies, describe organelle assembly, and review how to perform assembly quality control.

Intraspecific variations in Cyt b and D-loop sequences of Testudine species, Lissemys punctata from south Karnataka

The freshwater Testudine species have gained importance in recent years, as most of their population is threatened due to exploitation for delicacy and pet trade. In this regard, Lissemys punctata, a freshwater terrapin, predominantly distributed in Asian countries has gained its significance for the study. A pilot study report on mitochondrial markers (Cyt b and D-loop) conducted on L. punctata species from southern Karnataka, India was presented in this investigation. A complete region spanning 1.14 kb and ∼1 kb was amplified by HotStart PCR and sequenced by Sanger sequencing. The Cyt b sequence revealed 85 substitution sites, no indels and 17 parsimony informative sites, whereas D-loop showed 189 variable sites, 51 parsimony informative sites with 5′ functional domains TAS, CSB-F, CSBs (1, 2, 3) preceding tandem repeat at 3′ end. Current data highlights the intraspecific variations in these target regions and variations validated using suitable evolutionary models points out that the overall point mutations observed in the region are transitions leading to no structural and functional alterations. The mitochondrial data generated uncover the genetic diversity within species and conservationist can utilize the data to estimate the effective population size or for forensic identification of animal or its seizures during unlawful trade activities.

Mitochondrial mutations in human cancer: Curation of translation

As a genetic disease, cancer is caused by the activation of oncogenes and the inhibition of tumor suppressor genes via genetic and epigenetic mechanisms. Given the important role of energy metabolism in tumors, we analyzed the cancer-derived mutations occurring in the DNA of the mitochondrion. Mutations in the mitochondrial DNA (mtDNA) compared to nuclear DNA are 62% decreased relative to the coding length per chromosome. We find that the majority of these mutations affects highly conserved nucleotides - significantly exceeding the conservation of the mtDNA - and are devoid of single nucleotide polymorphisms (SNPs). Surprisingly, the leading resources for tumor genetics information universally use the standard genetic code for translation of nucleotide into amino acid sequences in their online resources. However, the nuclear and mitochondrial genetic codes differ for four codons and the usage of incomplete STOP codons. Hence, we analyze and curate the consequences for all mutations in the mtDNA and comprehensively reclassify missense, nonsense and synonymous mutations accordingly. In total, 10% of the mutations are incorrectly translated leading to significant changes in the distribution of mutation types with tripling of nonsense and 69% loss of nonstop extension mutations. Lastly, we provide a curated dataset of coding and non-coding mitochondrial mutations in cancer merged, standardized, duplicate-free and aggregated from two databases as a resource including orthogonal data on their high conservation and SNPs. This study generally highlights the need to universally regard the important differences between the standard and mitochondrial genetic code in life science research.

A close-up view of codon selection in eukaryotic initiation

When given an option to choose among a set of alternatives and only one selection is right, one might stop and reflect over which one is best. However, the ribosome has no time to stop and make such reflections, proteins need to be produced and very fast. Eukaryotic translation initiation is an example of such a conundrum. Here, scanning for the correct codon match must be fast, efficient and accurate. We highlight our recent computational findings, which show how the initiation machinery manages to recognize one specific codon among many possible challengers, by fine-tuning the energetic landscape of base-pairing with the aid of the initiation factors eIF1 and eIF1A. Using a recent 3-dimensional structure of the eukaryotic initiation complex we have performed simulations of codon recognition in atomic detail. These calculations provide an in-depth energetic and structural view of how discrimination against near-cognate codons is achieved by the initiation complex.

Differential bicodon usage in lowly and highly abundant proteins

Degeneracy in the genetic code implies that different codons can encode the same amino acid. Usage preference of synonymous codons has been observed in all domains of life. There is much evidence suggesting that this bias has a major role on protein elongation rate, contributing to differential expression and to co-translational folding. In addition to codon usage bias, other preference variations have been observed such as codon pairs. In this paper, I report that codon pairs have significant different frequency usage for coding either lowly or highly abundant proteins. These usage preferences cannot be explained by the frequency usage of the single codons. The statistical analysis of coding sequences of nine organisms reveals that in many cases bicodon preferences are shared between related organisms. Furthermore, it is observed that misfolding in the drug-transport protein, encoded by MDR1 gene, is better explained by a big change in the pause propensity due to the synonymous bicodon variant, rather than by a relatively small change in codon usage. These findings suggest that codon pair usage can be a more powerful framework to understand translation elongation rate, protein folding efficiency, and to improve protocols to optimize heterologous gene expression.

How tRNAs dictate nuclear codon reassignments: Only a few can capture non-cognate codons

mRNA decoding by tRNAs and tRNA charging by aminoacyl-tRNA synthetases are biochemically separated processes that nevertheless in general involve the same nucleotides. The combination of charging and decoding determines the genetic code. Codon reassignment happens when a differently charged tRNA replaces a former cognate tRNA. The recent discovery of the polyphyly of the yeast CUG sense codon reassignment challenged previous mechanistic considerations and led to the proposal of the so-called tRNA loss driven codon reassignment hypothesis. Accordingly, codon capture is caused by loss of a tRNA or by mutations in the translation termination factor, subsequent reduction of the codon frequency through reduced translation fidelity and final appearance of a new cognate tRNA. Critical for codon capture are sequence and structure of the new tRNA, which must be compatible with recognition regions of aminoacyl-tRNA synthetases. The proposed hypothesis applies to all reported nuclear and organellar codon reassignments.

Two complete mitochondrial genomes from Praticolella mexicana Perez, 2011 (Polygyridae) and gene order evolution in Helicoidea (Mollusca, Gastropoda)

Helicoidea is a diverse group of land snails with a global distribution. While much is known regarding the relationships of helicoid taxa, comparatively little is known about the evolution of the mitochondrial genome in the superfamily. We sequenced two complete mitochondrial genomes from Praticolellamexicana Perez, 2011 representing the first such data from the helicoid family Polygyridae, and used them in an evolutionary analysis of mitogenomic gene order. We found the mitochondrial genome of Praticolellamexicana to be 14,008 bp in size, possessing the typical 37 metazoan genes. Multiple alternate stop codons are used, as are incomplete stop codons. Mitogenome size and nucleotide content is consistent with other helicoid species. Our analysis of gene order suggested that Helicoidea has undergone four mitochondrial rearrangements in the past. Two rearrangements were limited to tRNA genes only, and two involved protein coding genes.

Cross-species conservation of complementary amino acid-ribonucleobase interactions and their potential for ribosome-free encoding

The role of amino acid-RNA nucleobase interactions in the evolution of RNA translation and protein-mRNA autoregulation remains an open area of research. We describe the inference of pairwise amino acid-RNA nucleobase interaction preferences using structural data from known RNA-protein complexes. We observed significant matching between an amino acid’s nucleobase affinity and corresponding codon content in both the standard genetic code and mitochondrial variants. Furthermore, we showed that knowledge of nucleobase preferences allows statistically significant prediction of protein primary sequence from mRNA using purely physiochemical information. Interestingly, ribosomal primary sequences were more accurately predicted than non-ribosomal sequences, suggesting a potential role for direct amino acid-nucleobase interactions in the genesis of amino acid-based ribosomal components. Finally, we observed matching between amino acid-nucleobase affinities and corresponding mRNA sequences in 35 evolutionarily diverse proteomes. We believe these results have important implications for the study of the evolutionary origins of the genetic code and protein-mRNA cross-regulation.

Chemical Evolution of a Bacterial Proteome

We have changed the amino acid set of the genetic code of Escherichia coli by evolving cultures capable of growing on the synthetic noncanonical amino acid L-β-(thieno[3,2-b]pyrrolyl)alanine ([3,2]Tpa) as a sole surrogate for the canonical amino acid L-tryptophan (Trp). A long-term cultivation experiment in defined synthetic media resulted in the evolution of cells capable of surviving Trp→[3,2]Tpa substitutions in their proteomes in response to the 20,899 TGG codons of the E. coli W3110 genome. These evolved bacteria with new-to-nature amino acid composition showed robust growth in the complete absence of Trp. Our experimental results illustrate an approach for the evolution of synthetic cells with alternative biochemical building blocks.

Study of the Binding Energies between Unnatural Amino Acids and Engineered Orthogonal Tyrosyl-tRNA Synthetases

We utilized several computational approaches to evaluate the binding energies of tyrosine (Tyr) and several unnatural Tyr analogs, to several orthogonal aaRSes derived from Methanocaldococcus jannaschii and Escherichia coli tyrosyl-tRNA synthetases. The present study reveals the following: (1) AutoDock Vina and ROSETTA were able to distinguish binding energy differences for individual pairs of favorable and unfavorable aaRS-amino acid complexes, but were unable to cluster together all experimentally verified favorable complexes from unfavorable aaRS-Tyr complexes; (2) MD-MM/PBSA provided the best prediction accuracy in terms of clustering favorable and unfavorable enzyme-substrate complexes, but also required the highest computational cost; and (3) MM/PBSA based on single energy-minimized structures has a significantly lower computational cost compared to MD-MM/PBSA, but still produced sufficiently accurate predictions to cluster aaRS-amino acid interactions. Although amino acid-aaRS binding is just the first step in a complex series of processes to acylate a tRNA with its corresponding amino acid, the difference in binding energy, as shown by MD-MM/PBSA, is important for a mutant orthogonal aaRS to distinguish between a favorable unnatural amino acid (unAA) substrate from unfavorable natural amino acid substrates. Our computational study should assist further designing and engineering of orthogonal aaRSes for the genetic encoding of novel unAAs.

Ancient DNA sequence revealed by error-correcting codes

A previously described DNA sequence generator algorithm (DNA-SGA) using error-correcting codes has been employed as a computational tool to address the evolutionary pathway of the genetic code. The code-generated sequence alignment demonstrated that a residue mutation revealed by the code can be found in the same position in sequences of distantly related taxa. Furthermore, the code-generated sequences do not promote amino acid changes in the deviant genomes through codon reassignment. A Bayesian evolutionary analysis of both code-generated and homologous sequences of the Arabidopsis thaliana malate dehydrogenase gene indicates an approximately 1 MYA divergence time from the MDH code-generated sequence node to its paralogous sequences. The DNA-SGA helps to determine the plesiomorphic state of DNA sequences because a single nucleotide alteration often occurs in distantly related taxa and can be found in the alternative codon patterns of noncanonical genetic codes. As a consequence, the algorithm may reveal an earlier stage of the evolution of the standard code.

On models of the genetic code generated by binary dichotomic algorithms

In this paper we introduce the concept of a BDA-generated model of the genetic code which is based on binary dichotomic algorithms (BDAs). A BDA-generated model is based on binary dichotomic algorithms (BDAs). Such a BDA partitions the set of 64 codons into two disjoint classes of size 32 each and provides a generalization of known partitions like the Rumer dichotomy. We investigate what partitions can be generated when a set of different BDAs is applied sequentially to the set of codons. The search revealed that these models are able to generate code tables with very different numbers of classes ranging from 2 to 64. We have analyzed whether there are models that map the codons to their amino acids. A perfect matching is not possible. However, we present models that describe the standard genetic code with only few errors. There are also models that map all 64 codons uniquely to 64 classes showing that BDAs can be used to identify codons precisely. This could serve as a basis for further mathematical analysis using coding theory, for example. The hypothesis that BDAs might reflect a molecular mechanism taking place in the decoding center of the ribosome is discussed. The scan demonstrated that binary dichotomic partitions are able to model different aspects of the genetic code very well. The search was performed with our tool Beady-A. This software is freely available at http://mi.informatik.hs-mannheim.de/beady-a. It requires a JVM version 6 or higher.

Are "universal" DNA primers really universal?

"Universal" DNA primers LCO 1490 and HCO 2198 were originally designed from three coding and six anticoding strands by comparing highly conserved regions of mitochondrial cytochrome c oxidase subunit I (COI) genes across 15 taxa. These primers have been successful in amplifying a 710-bp fragment of highly conserved regions of the COI gene for more than 80 invertebrate species from 11 phyla. In the present study, 130,843 variations were reviewed in the primer region of mitochondrial molecular markers by comparing 725 COI sequences from the kingdom Animalia. It was found that, for 177 invertebrate species, the forward primer (LCO 1490) showed only four conserved regions, compared to 12 in the original study. For ascidians, fungi and vertebrates, it showed approximately 50 % conserved regions, dropping to one conserved region for echinoderms. However, the reverse primer (HCO 2198) was highly conserved across 725 COI primer sequences. A similar pattern was observed in amino acid distributions. There was a significant difference in the means of base pair differences from the level of family, genus and species for LCO 1490 [analysis of variance (ANOVA), F 6,188 = 8.193, P < 0.001] and at the level of genus and species for HCO 2198 (ANOVA, F 6,77 = 2.538, P < 0.027). We conclude that, at different taxonomic levels, it is possible to design forward primers from reference sequences belonging to the level of order (maximum 5 bp differences), family (maximum 6 bp differences) or genus (maximum 1 bp difference). Reverse primers can be designed from the level of family (maximum 5 bp differences) or genus (maximum 2 bp differences).

A symmetry model for genetic coding via a wallpaper group composed of the traditional four bases and an imaginary base E: towards category theory-like systematization of molecular/genetic biology

Background

Previously, we suggested prototypal models that describe some clinical states based on group postulates. Here, we demonstrate a group/category theory-like model for molecular/genetic biology as an alternative application of our previous model. Specifically, we focus on deoxyribonucleic acid (DNA) base sequences.

Results

We construct a wallpaper pattern based on a five-letter cruciform motif with letters C, A, T, G, and E. Whereas the first four letters represent the standard DNA bases, the fifth is introduced for ease in formulating group operations that reproduce insertions and deletions of DNA base sequences. A basic group Z₅ = {r, u, d, l, n} of operations is defined for the wallpaper pattern, with which a sequence of points can be generated corresponding to changes of a base in a DNA sequence by following the orbit of a point of the pattern under operations in group Z₅. Other manipulations of DNA sequence can be treated using a vector-like notation ‘D_j’ corresponding to a DNA sequence but based on the five-letter base set; also, ‘D_j’s are expressed graphically. Insertions and deletions of a series of letters ‘E’ are admitted to assist in describing DNA recombination. Likewise, a vector-like notation R_j can be constructed for sequences of ribonucleic acid (RNA). The wallpaper group B = {Z₅^×∞, ●} (an ∞-fold Cartesian product of Z₅) acts on D_j (or R_j) yielding changes to D_j (or R_j) denoted by ‘D_j◦B_(j→k) = D_k’ (or ‘R_j◦B_(j→k) = R_k’). Based on the operations of this group, two types of groups—a modulo 5 linear group and a rotational group over the Gaussian plane, acting on the five bases—are linked as parts of the wallpaper group for broader applications. As a result, changes, insertions/deletions and DNA (RNA) recombination (partial/total conversion) are described. As an exploratory study, a notation for the canonical “central dogma” via a category theory-like way is presented for future developments.

Conclusions

Despite the large incompleteness of our methodology, there is fertile ground to consider a symmetry model for genetic coding based on our specific wallpaper group. A more integrated formulation containing “central dogma” for future molecular/genetic biology remains to be explored.

Molecular codes in biological and chemical reaction networks

Shannon’s theory of communication has been very successfully applied for the analysis of biological information. However, the theory neglects semantic and pragmatic aspects and thus cannot directly be applied to distinguish between (bio-) chemical systems able to process “meaningful” information from those that do not. Here, we present a formal method to assess a system’s semantic capacity by analyzing a reaction network’s capability to implement molecular codes. We analyzed models of chemical systems (martian atmosphere chemistry and various combustion chemistries), biochemical systems (gene expression, gene translation, and phosphorylation signaling cascades), an artificial chemistry, and random reaction networks. Our study suggests that different chemical systems posses different semantic capacities. No semantic capacity was found in the model of the martian atmosphere chemistry, the studied combustion chemistries, and highly connected random networks, i.e. with these chemistries molecular codes cannot be implemented. High semantic capacity was found in the studied biochemical systems and in random reaction networks where the number of second order reactions is twice the number of species. We conclude that our approach can be applied to evaluate the information processing capabilities of a chemical system and may thus be a useful tool to understand the origin and evolution of meaningful information, e.g. in the context of the origin of life.

A universal trend among proteomes indicates an oily last common ancestor

Despite progresses in ancestral protein sequence reconstruction, much needs to be unraveled about the nature of the putative last common ancestral proteome that served as the prototype of all extant lifeforms. Here, we present data that indicate a steady decline (oil escape) in proteome hydrophobicity over species evolvedness (node number) evident in 272 diverse proteomes, which indicates a highly hydrophobic (oily) last common ancestor (LCA). This trend, obtained from simple considerations (free from sequence reconstruction methods), was corroborated by regression studies within homologous and orthologous protein clusters as well as phylogenetic estimates of the ancestral oil content. While indicating an inherent irreversibility in molecular evolution, oil escape also serves as a rare and universal reaction-coordinate for evolution (reinforcing Darwin's principle of Common Descent), and may prove important in matters such as (i) explaining the emergence of intrinsically disordered proteins, (ii) developing composition- and speciation-based "global" molecular clocks, and (iii) improving the statistical methods for ancestral sequence reconstruction.

Transfer RNA modifications: nature's combinatorial chemistry playground

Following synthesis, tRNAs are peppered by numerous chemical modifications which may differentially affect a tRNA's structure and function. Although modifications affecting the business ends of a tRNA are predictably important for cell viability, a majority of modifications play more subtle structural roles that can affect tRNA stability and folding. The current trend is that modifications act in concert and it is in the context of the specific sequence of a given tRNA that they impart their differing effects. Recent developments in the modification field have highlighted the diversity of modifications in tRNA. From these, the combinatorial nature of modifications in explaining previously described phenotypes derived from their absence has emerged as a growing theme.

Updating our view of organelle genome nucleotide landscape

Organelle genomes show remarkable variation in architecture and coding content, yet their nucleotide composition is relatively unvarying across the eukaryotic domain, with most having a high adenine and thymine (AT) content. Recent studies, however, have uncovered guanine and cytosine (GC)-rich mitochondrial and plastid genomes. These sequences come from a small but eclectic list of species, including certain green plants and animals. Here, I review GC-rich organelle DNAs and the insights they have provided into the evolution of nucleotide landscape. I emphasize that GC-biased mitochondrial and plastid DNAs are more widespread than once thought, sometimes occurring together in the same species, and suggest that the forces biasing their nucleotide content can differ both among and within lineages, and may be associated with specific genome architectural features and life history traits.

Biochemical analysis with the expanded genetic lexicon

The information used to build proteins is stored in the genetic material of every organism. In nature, ribosomes use 20 native amino acids to synthesize proteins in most circumstances. However, laboratory efforts to expand the genetic repertoire of living cells and organisms have successfully encoded more than 80 nonnative amino acids in E. coli, yeast, and other eukaryotic systems. The selectivity, fidelity, and site-specificity provided by the technology have enabled unprecedented flexibility in manipulating protein sequences and functions in cells. Various biophysical probes can be chemically conjugated or directly incorporated at specific residues in proteins, and corresponding analytical techniques can then be used to answer diverse biological questions. This review summarizes the methodology of genetic code expansion and its recent progress, and discusses the applications of commonly used analytical methods.

NovelSNPer: A Fast Tool for the Identification and Characterization of Novel SNPs and InDels

Typically, next-generation resequencing projects produce large lists of variants. NovelSNPer is a software tool that permits fast and efficient processing of such output lists. In a first step, NovelSNPer determines if a variant represents a known variant or a previously unknown variant. In a second step, each variant is classified into one of 15 SNP classes or 19 InDel classes. Beside the classes used by Ensembl, we introduce POTENTIAL_START_GAINED and START_LOST as new functional classes and present a classification scheme for InDels. NovelSNPer is based upon the gene structure information stored in Ensembl. It processes two million SNPs in six hours. The tool can be used online or downloaded.

Stops making sense: translational trade-offs and stop codon reassignment

Background

Efficient gene expression involves a trade-off between (i) premature termination of protein synthesis; and (ii) readthrough, where the ribosome fails to dissociate at the terminal stop. Sense codons that are similar in sequence to stop codons are more susceptible to nonsense mutation, and are also likely to be more susceptible to transcriptional or translational errors causing premature termination. We therefore expect this trade-off to be influenced by the number of stop codons in the genetic code. Although genetic codes are highly constrained, stop codon number appears to be their most volatile feature.

Results

In the human genome, codons readily mutable to stops are underrepresented in coding sequences. We construct a simple mathematical model based on the relative likelihoods of premature termination and readthrough. When readthrough occurs, the resultant protein has a tail of amino acid residues incorrectly added to the C-terminus. Our results depend strongly on the number of stop codons in the genetic code. When the code has more stop codons, premature termination is relatively more likely, particularly for longer genes. When the code has fewer stop codons, the length of the tail added by readthrough will, on average, be longer, and thus more deleterious. Comparative analysis of taxa with a range of stop codon numbers suggests that genomes whose code includes more stop codons have shorter coding sequences.

Conclusions

We suggest that the differing trade-offs presented by alternative genetic codes may result in differences in genome structure. More speculatively, multiple stop codons may mitigate readthrough, counteracting the disadvantage of a higher rate of nonsense mutation. This could help explain the puzzling overrepresentation of stop codons in the canonical genetic code and most variants.

Properties of biological mutation networks and their implications for ALife

We report a study of networks constructed from mutation patterns observed in biology. These networks form evolutionary trajectories, which allow for both frequent substitution of closely related structures, and a small evolutionary distance between any two structures. These two properties define the small-world phenomenon. The mutation behavior between tokens in an evolvable artificial chemistry determines its ability to explore evolutionary space. This concept is underrepresented in previous work on string-based chemistries. We argue that small-world mutation networks will confer better exploration of the evolutionary space than either random or fully regular mutation strategies. We calculate network statistics from two data sets: amino acid substitution matrices, and codon-level single point mutations. The first class are observed data from protein alignments; while the second class is defined by the standard genetic code that is used to translate RNA into amino acids. We report a methodology for creating small-world mutation networks for artificial chemistries with arbitrary node count and connectivity. We argue that ALife systems would benefit from this approach, as it delivers a more viable exploration of evolutionary space.

Complex phylogenetic distribution of a non-canonical genetic code in green algae

BACKGROUND

A non-canonical nuclear genetic code, in which TAG and TAA have been reassigned from stop codons to glutamine, has evolved independently in several eukaryotic lineages, including the ulvophycean green algal orders Dasycladales and Cladophorales. To study the phylogenetic distribution of the standard and non-canonical genetic codes, we generated sequence data of a representative set of ulvophycean green algae and used a robust green algal phylogeny to evaluate different evolutionary scenarios that may account for the origin of the non-canonical code.

RESULTS

This study demonstrates that the Dasycladales and Cladophorales share this alternative genetic code with the related order Trentepohliales and the genus Blastophysa, but not with the Bryopsidales, which is sister to the Dasycladales. This complex phylogenetic distribution whereby all but one representative of a single natural lineage possesses an identical deviant genetic code is unique.

CONCLUSIONS

We compare different evolutionary scenarios for the complex phylogenetic distribution of this non-canonical genetic code. A single transition to the non-canonical code followed by a reversal to the canonical code in the Bryopsidales is highly improbable due to the profound genetic changes that coincide with codon reassignment. Multiple independent gains of the non-canonical code, as hypothesized for ciliates, are also unlikely because the same deviant code has evolved in all lineages. Instead we favor a stepwise acquisition model, congruent with the ambiguous intermediate model, whereby the non-canonical code observed in these green algal orders has a single origin. We suggest that the final steps from an ambiguous intermediate situation to a non-canonical code have been completed in the Trentepohliales, Dasycladales, Cladophorales and Blastophysa but not in the Bryopsidales. We hypothesize that in the latter lineage an initial stage characterized by translational ambiguity was not followed by final reassignment of both stop codons to glutamine. Instead the standard code was retained by the disappearance of the ambiguously decoding tRNAs from the genome. We correlate the emergence of a non-canonical genetic code in the Ulvophyceae to their multinucleate nature.

Polyunsaturated fatty acids modulate prostaglandin synthesis by ovine amnion cells in vitro

Diets or supplements high in n-3 and n-6 polyunsaturated fatty acids (PUFAs) have been shown to influence the timing of parturition. PUFAs are substrates for prostaglandin (PG) synthesis, and PGs play central roles in parturition. Hence, the effects of altering PUFA composition may be mediated through alterations in the type and relative quantities of PGs synthesised. Therefore, we have investigated the effects of a range of n-3 and n-6 PUFAs in vitro on PG synthesis by amnion cells of late gestation ewes. The n-6 PUFA, arachidonic acid (20:4, n-6), increased synthesis of two-series PGs. Degree of stimulation induced by the n-6 PUFAs was dependent on the position of the PUFA in the PG synthetic pathway, i.e. PG production of the two-series (principally prostaglandin E(2):PGE(2)) increased progressively with longer chain PUFAs. Effects of n-3 PUFAs on output of PGE(2) were more modest and variable. The two shorter chain n-3 PUFAs, α-linolenic acid (18:3, n-3) and stearidonic acid (18:4, n-3), induced a small but significant increase in PGE(2) output, while the longest chain n-3 PUFA docosahexaenoic acid (22:6, n-3) inhibited PGE(2) synthesis. Dihomo-γ-linolenic acid (20:3, n-6), the PUFA substrate for synthesis of one-series PGs, induced an increase in PGE(1) generation and a decrease in PGE(2) and PGE(3) outputs. Hence, we have demonstrated that PUFA supplementation of ovine amnion cells in vitro affects the type and quantity of PGs synthesised.

Dual functions of codons in the genetic code

The discovery of the genetic code provided one of the basic foundations of modern molecular biology. Most organisms use the same genetic language, but there are also well-documented variations representing codon reassignments within specific groups of organisms (such as ciliates and yeast) or organelles (such as plastids and mitochondria). In addition, duality in codon function is known in the use of AUG in translation initiation and methionine insertion into internal protein positions as well as in the case of selenocysteine and pyrrolysine insertion (encoded by UGA and UAG, respectively) in competition with translation termination. Ambiguous meaning of CUG in coding for serine and leucine is also known. However, a recent study revealed that codons in any position within the open reading frame can serve a dual function and that a change in codon meaning can be achieved by availability of a specific type of RNA stem-loop structure in the 3'-untranslated region. Thus, duality of codon function is a more widely used feature of the genetic code than previously known, and this observation raises the possibility that additional recoding events and additional novel features have evolved in the genetic code.

A comparative genomics analysis of codon reassignments reveals a link with mitochondrial proteome size and a mechanism of genetic code change via suppressor tRNAs

Using a comparative genomics approach we demonstrate a negative correlation between the number of codon reassignments undergone by 222 mitochondrial genomes and the mitochondrial genome size, the number of mitochondrial ORFs, and the sizes of the large and small subunit mitochondrial rRNAs. In addition, we show that the TGA-to-tryptophan codon reassignment, which has occurred 11 times in mitochondrial genomes, is found in mitochondrial genomes smaller than those which have not undergone the reassignment. We therefore propose that mitochondrial codon reassignments occur in a wide range of phyla, particularly in Metazoa, due to a reduced "proteomic constraint" on the mitochondrial genetic code, compared to the nuclear genetic code. The reduced proteomic constraint reflects the small size of the mitochondrial-encoded proteome and allows codon reassignments to occur with less likelihood of lethality. In addition, we demonstrate a striking link between nonsense codon reassignments and the decoding properties of naturally occurring nonsense suppressor tRNAs. This suggests that natural preexisting nonsense suppression facilitated nonsense codon reassignments and constitutes a novel mechanism of genetic code change. These findings explain for the first time the identity of the stop codons and amino acids reassigned in mitochondrial and nuclear genomes. Nonsense suppressor tRNAs provided the raw material for nonsense codon reassignments, implying that the properties of the tRNA anticodon have dictated the identity of nonsense codon reassignments.

The yeast translation release factors Mrf1p and Sup45p (eRF1) are methylated, respectively, by the methyltransferases Mtq1p and Mtq2p

The translation release factors (RFs) RF1 and RF2 of Escherichia coli are methylated at the N5-glutamine of the GGQ motif by PrmC methyltransferase. This motif is conserved in organisms from bacteria to higher eukaryotes. The Saccharomyces cerevisiae RFs, mitochondrial Mrf1p and cytoplasmic Sup45p (eRF1), have sequence similarities to the bacterial RFs, including the potential site of glutamine methylation in the GGQ motif. A computational analysis revealed two yeast proteins, Mtq1p and Mtq2p, that have strong sequence similarity to PrmC. Mass spectrometric analysis demonstrated that Mtq1p and Mtq2p methylate Mrf1p and Sup45p, respectively, in vivo. A tryptic peptide of Mrf1p, GGQHVNTTDSAVR, containing the GGQ motif was found to be approximately 50% methylated at the glutamine residue in the normal strain but completely unmodified in the peptide from mtq1-Delta. Moreover, Mtq1p methyltransferase activity was observed in an in vitro assay. In similar experiments, it was determined that Mtq2p methylates Sup45p. The Sup45p methylation by Mtq2p was recently confirmed independently (Heurgue-Hamard, V., Champ, S., Mora, L., Merkulova-Rainon, T., Kisselev, L. L., and Buckingham, R. H. (2005) J. Biol. Chem. 280, 2439-2445). Analysis of the deletion mutants showed that although mtq1-Delta had only moderate growth defects on nonfermentable carbon sources, the mtq2-Delta had multiple phenotypes, including cold sensitivity and sensitivity to translation fidelity antibiotics paromomycin and geneticin, to high salt and calcium concentrations, to polymyxin B, and to caffeine. Also, the mitochondrial mit(-) mutation, cox2-V25, containing a premature stop mutation, was suppressed by mtq1-Delta. Most interestingly, the mtq2-Delta was significantly more resistant to the anti-microtubule drugs thiabendazole and benomyl, suggesting that Mtq2p may also methylate certain microtubule-related proteins.

Origins of the genetic code: the escaped triplet theory

There is very significant evidence that cognate codons and/or anticodons are unexpectedly frequent in RNA-binding sites for seven of eight biological amino acids that have been tested. This suggests that a substantial fraction of the genetic code has a stereochemical basis, the triplets having escaped from their original function in amino acid-binding sites to become modern codons and anticodons. We explicitly show that this stereochemical basis is consistent with subsequent optimization of the code to minimize the effect of coding mistakes on protein structure. These data also strengthen the argument for invention of the genetic code in an RNA world and for the RNA world itself.

Gene algebra from a genetic code algebraic structure

By considering two important factors involved in the codon-anticodon interactions, the hydrogen bond number and the chemical type of bases, a codon array of the genetic code table as an increasing code scale of interaction energies of amino acids in proteins was obtained. Next, in order to consecutively obtain all codons from the codon AAC, a sum operation has been introduced in the set of codons. The group obtained over the set of codons is isomorphic to the group (Z(64), +) of the integer module 64. On the Z(64)-algebra of the set of 64(N) codon sequences of length N, gene mutations are described by means of endomorphisms f:(Z(64))(N)-->(Z(64))(N). Endomorphisms and automorphisms helped us describe the gene mutation pathways. For instance, 77.7% mutations in 749 HIV protease gene sequences correspond to unique diagonal endomorphisms of the wild type strain HXB2. In particular, most of the reported mutations that confer drug resistance to the HIV protease gene correspond to diagonal automorphisms of the wild type. What is more, in the human beta-globin gene a similar situation appears where most of the single codon mutations correspond to automorphisms. Hence, in the analyses of molecular evolution process on the DNA sequence set of length N, the Z(64)-algebra will help us explain the quantitative relationships between genes.