Codon Length Research Articles

Comparative analysis of codon usage patterns of Plasmodium helical interspersed subtelomeric (PHIST) proteins.

Plasmodium falciparum is a protozoan parasite that causes the most severe form of malaria in humans worldwide, which is predominantly found in sub-Saharan Africa, where it is responsible for the majority of malaria-related deaths. Plasmodium helical interspersed subtelomeric (PHIST) proteins are a family of proteins, with a conserved PHIST domain, which are typically located at the subtelomeric regions of the Plasmodium falciparum chromosomes and play crucial roles in the interaction between the parasite and its human host, such as cytoadherence, immune evasion, and host cell remodeling. However, the specific utilization of synonymous codons by PHIST proteins in Plasmodium falciparum is still unknown. Codon usage bias (CUB) refers to the unequal usage of synonymous codons during translation, resulting in over- or underrepresentation of certain nucleotide patterns. This imbalance in CUB can impact various cellular processes, including protein expression levels and genetic variation. To investigate this, the CUB of 88 PHIST protein coding sequences (CDSs) from 5 subgroups were analyzed in this study. The results showed that both codon base composition and relative synonymous codon usage (RSCU) analysis identified a higher occurrence of AT-ended codons (AGA and UUA) in PHIST proteins of Plasmodium falciparum. The average effective number of codons (ENC) for these PHIST proteins was 36.69, indicating a weak codon preference among them, as it was greater than 35. Additionally, the correlation analysis among codon base composition (GC1, GC2, GC3, GCs), codon adaptation index (CAI), codon bias index (CBI), frequency of optimal codons (FOP), ENC, general average hydropathicity (GRAVY), aromaticity (AROMO), length of synonymous codons (L_sym), and length of amino acids (L_aa) revealed the influence of base composition and codon usage indices on codon usage bias, with GC1 having a significant impact in this study. Furthermore, the neutrality plot analysis, PR2-bias plot analysis, and ENC-GC3 plot analysis provided additional evidence that natural selection plays a crucial role in determining codon bias in PHIST proteins. In conclusion, this study has enhanced our understanding of the characteristics of codon usage and genetic evolution in PHIST proteins, thereby providing data foundation for further research on antimalarial drugs or vaccines.

Analysis of codon usage bias of thioredoxin in apicomplexan protozoa

BackgroundApicomplexan protozoa are a diverse group of obligate intracellular parasites causing many diseases that affect humans and animals, such as malaria, toxoplasmosis, and cryptosporidiosis. Apicomplexan protozoa possess unique thioredoxins (Trxs) that have been shown to regulate various cellular processes including metabolic redox regulation, parasite survival, and host immune evasion. However, it is still unknown how synonymous codons are used by apicomplexan protozoa Trxs.MethodsCodon usage bias (CUB) is the unequal usage of synonymous codons during translation which leads to the over- or underrepresentation of certain nucleotide patterns. This imbalance in CUB can impact a variety of cellular processes including protein expression levels and genetic variation. This study analyzed the CUB of 32 Trx coding sequences (CDS) from 11 apicomplexan protozoa.ResultsThe results showed that both codon base composition and relative synonymous codon usage (RSCU) analysis revealed that AT-ended codons were more frequently used in Cryptosporidium spp. and Plasmodium spp., while the Eimeria spp., Babesia spp., Hammondia hammondi, Neospora caninum, and Toxoplasma gondii tended to end in G/C. The average effective number of codon (ENC) value of these apicomplexan protozoa is 46.59, which is > 35, indicating a weak codon preference among apicomplexan protozoa Trxs. Furthermore, the correlation analysis among codon base composition (GC1, GC2, GC3, GCs), codon adaptation index (CAI), codon bias index (CBI), frequency of optimal codons (FOP), ENC, general average hydropathicity (GRAVY), aromaticity (AROMO), length of synonymous codons (L_sym), and length of amino acids (L_aa) indicated the influence of base composition and codon usage indices on CUB. Additionally, the neutrality plot analysis, PR2-bias plot analysis, and ENC-GC3 plot analysis further demonstrated that natural selection plays an important role in apicomplexan protozoa Trxs codon bias.ConclusionsIn conclusion, this study increased the understanding of codon usage characteristics and genetic evolution of apicomplexan protozoa Trxs, which expanded new ideas for vaccine and drug research.Graphical

An Adaptive Mapping Method Using Spectral Envelope Approach for DNA Spectral Analysis.

The digital signal processing approaches were investigated as a preliminary indicator for discriminating between the protein coding and non-coding regions of DNA. This is because a three-base periodicity (TBP) has already been proven to exist in protein-coding regions arising from the length of codons (three nucleic acids). This demonstrates that there is a prominent peak in the energy spectrum of a DNA coding sequence at frequency rad/sample. However, because DNA sequences are symbolic sequences, these should be mapped into one or more signals such that the hidden information is highlighted. We propose, therefore, two new algorithms for computing adaptive mappings and, by using them, finding periodicities. Both such algorithms are based on the spectral envelope approach. This adaptive approach is essentially important since a single mapping for any DNA sequence may ignore its intrinsic properties. Finally, the improved performance of the new methods is verified by using them with synthetic and real DNA sequences as compared to the classical methods, especially the minimum entropy mapping (MEM) spectrum, which is also an adaptive method. We demonstrated that our method is both more accurate and more responsive than all its counterparts. This is especially important in this application since it reduces the risks of a coding sequence being missed.

Trinity of G-tetrads and origin of translation

BackgroundThe RNA world hypothesis cannot address most of the questions of the origin of life without violating the continuity principle (small Darwinian steps without foresight and miracles). Moreover, the RNA world is an isolated system incapable of accommodating the genetic code and evolving into extant biochemistry. All these problems are rooted in the central assumption of the hypothesis: de novo appearance of the ribozymes, production of which represents a multistep reaction requiring the complementarity principle. Thus, even the basis of the RNA world is at odds with the continuity principle—it uses foresight (multistep reaction) and a miracle (complementarity principle). Can a three-dimensional (3D) architecture, capable of molecular recognition and catalysis, be formed in a single-step reaction without the complementarity or any other preexisting rules?HypothesisAt first glance, the above question sounds rhetoric since the complementarity principle is the essential feature of the RNA world; it turns an RNA polymer into a genetic material. Without it, the RNA world becomes as shapeless and unconvincing as other hypotheses based on the non-hereditary molecules (i.e., protein world). However, it was suggested recently that the quadruplexes could initiate life and take necessary evolutionary steps before the arrival of the complementarity rules. The hypothesis relies on the unique properties of guanines (Gs) to self-assemble into G-tetrads and efficiently polymerize without any external help or preexisting rules. Interestingly, polyG folds into an unusually stable and well-structured monomolecular architecture that uses the quadruplex domain (QD) assembly. The QD has a strictly defined zigzag-like building pattern to accommodate only three G-tetrads. Since both QD architecture and codon length are based on triplets, the inevitable question arises: are they related? Or could QD play the role of the early adapter and determine the codon length? The current paper is an attempt to answer this question.ConclusionWhile without translation apparatus most of the steps of the extant translation are physically impossible, the QD-mediated translation is sterically feasible and can be explained by physicochemical properties of the QD and the amino acids without violating the continuity principle. Astonishingly, the quadruplex world hypothesis can address all the shortcomings of the RNA world, including its most significant challenge—step-by-step evolution from the polymerization of the first polynucleotide to the extant biochemistry.

Identifying A- and P-site locations on ribosome-protected mRNA fragments using Integer Programming

Identifying the A- and P-site locations on ribosome-protected mRNA fragments from Ribo-Seq experiments is a fundamental step in the quantitative analysis of transcriptome-wide translation properties at the codon level. Many analyses of Ribo-Seq data have utilized heuristic approaches applied to a narrow range of fragment sizes to identify the A-site. In this study, we use Integer Programming to identify the A-site by maximizing an objective function that reflects the fact that the ribosome’s A-site on ribosome-protected fragments must reside between the second and stop codons of an mRNA. This identifies the A-site location as a function of the fragment’s size and its 5′ end reading frame in Ribo-Seq data generated from S. cerevisiae and mouse embryonic stem cells. The correctness of the identified A-site locations is demonstrated by showing that this method, as compared to others, yields the largest ribosome density at established stalling sites. By providing greater accuracy and utilization of a wider range of fragment sizes, our approach increases the signal-to-noise ratio of underlying biological signals associated with translation elongation at the codon length scale.

Phylogenetic tree construction using trinucleotide usage profile (TUP).

BackgroundIt has been a challenging task to build a genome-wide phylogenetic tree for a large group of species containing a large number of genes with long nucleotides sequences. The most popular method, called feature frequency profile (FFP-k), finds the frequency distribution for all words of certain length k over the whole genome sequence using (overlapping) windows of the same length. For a satisfactory result, the recommended word length (k) ranges from 6 to 15 and it may not be a multiple of 3 (codon length). The total number of possible words needed for FFP-k can range from 46=4096 to 415.ResultsWe propose a simple improvement over the popular FFP method using only a typical word length of 3. A new method, called Trinucleotide Usage Profile (TUP), is proposed based only on the (relative) frequency distribution using non-overlapping windows of length 3. The total number of possible words needed for TUP is 43=64, which is much less than the total count for the recommended optimal “resolution” for FFP. To build a phylogenetic tree, we propose first representing each of the species by a TUP vector and then using an appropriate distance measure between pairs of the TUP vectors for the tree construction. In particular, we propose summarizing a DNA sequence by a matrix of three rows corresponding to three reading frames, recording the frequency distribution of the non-overlapping words of length 3 in each of the reading frame. We also provide a numerical measure for comparing trees constructed with various methods.ConclusionsCompared to the FFP method, our empirical study showed that the proposed TUP method is more capable of building phylogenetic trees with a stronger biological support. We further provide some justifications on this from the information theory viewpoint. Unlike the FFP method, the TUP method takes the advantage that the starting of the first reading frame is (usually) known. Without this information, the FFP method could only rely on the frequency distribution of overlapping words, which is the average (or mixture) of the frequency distributions of three possible reading frames. Consequently, we show (from the entropy viewpoint) that the FFP procedure could dilute important gene information and therefore provides less accurate classification.

Stochastic alternative splicing is prevalent in mungbean (Vigna radiata)

SummaryAlternative splicing (AS) can produce multiple mature mRNAs from the same primary transcript, thereby generating diverse proteins and phenotypes from the same gene. To assess the prevalence of AS in mungbean (Vigna radiata), we analysed whole‐genome RNA sequencing data from root, leaf, flower and pod tissues and found that at least 37.9% of mungbean genes are subjected to AS. The number of AS transcripts exhibited a strong correlation with exon number and thus resembled a uniform probabilistic event rather than a specific regulatory function. The proportion of frameshift splicing was close to the expected frequency of random splicing. However, alternative donor and acceptor AS events tended to occur at multiples of three nucleotides (i.e. the codon length) from the main splice site. Genes with high exon number and expression level, which should have the most AS if splicing is purely stochastic, exhibited less AS, implying the existence of negative selection against excessive random AS. Functional AS is probably rare: a large proportion of AS isoforms exist at very low copy per cell on average or are expressed at much lower levels than default transcripts. Conserved AS was only detected in 629 genes (2.8% of all genes in the genome) when compared to Vigna angularis, and in 16 genes in more distant species like soya bean. These observations highlight the challenges of finding and cataloguing candidates for experimentally proven AS isoforms in a crop genome.

The complete mitochondrial genome of Fulica atra (Avian, Gruiformes, Rallidae)

To analyze the gene structure and the evolutionary roadmap of Fulica atra, the complete mitogenome is sequenced. It is composed of 37 genes and 1 control region, and the structure and arrangement of all genes are identical to other Rallidae. The comparative mitogenome revealed that the start codon and stop codon varied in Rallidae, and the gene lengths are different in ND2, COX1, ND3, ND5 and CYTB due to incompleteness of stop codon, frameshift mutation and various numbers of amino acids. We analyzed the correlation between phylogeny and gene characteristic in Rallidae with respect to the usage of start/stop codon and gene length, but no correlation was found. It indicates these discrepancies might happen independently. This work can afford an in-depth insight of phyletic evolution in Rallidae.

EF-G and EF4: translocation and back-translocation on the bacterial ribosome

Ribosomes translate the codon sequence of an mRNA into the amino acid sequence of the corresponding protein. One of the most crucial events is the translocation reaction, which involves movement of both the mRNA and the attached tRNAs by one codon length and is catalysed by the GTPase elongation factor G (EF-G). Interestingly, recent studies have identified a structurally related GTPase, EF4, that catalyses movement of the tRNA2-mRNA complex in the opposite direction when the ribosome stalls, which is known as back-translocation. In this Review, we describe recent insights into the mechanistic basis of both translocation and back-translocation.

DNA sequence watermarking based on random circular angle

This paper discusses DNA watermarking for copyright protection and authentication of a DNA sequence. We then propose a DNA watermarking method that confers mutation resistance, amino acid residue conservation, and watermark security. Our method allocates codons to random circular angles using a random mapping table and selects a number of codons for embedding targets using the Lipschitz regularity that is measured from the evolution across scales of local modulus maxima of codon circular angles. We then embed the watermark into random circular angles of codons without changing the amino acid residue. The length and location of target codons depend on the random mapping table and the singularity of detection of Lipschitz regularity. This table is used as the watermark key and can be applied to any codon sequence regardless of sequence length. Without knowledge of this table, it is very difficult to detect the length and location of sequences for extracting the watermark. From experimental results on the suitability of similar watermark capacities, we verified that our method has a lower bit-rate error for point mutations compared with previous methods. Further, we established that the entropies of the random mapping table and the location of target codons are high, indicating that the watermark is secure.

DNA as a medium for hiding data

Background Steganography is the science of hiding information by transmitting secret messages through unsuspicious cover carriers in a way that makes the presence of any embedded messages undetectable. The term originates in the Greek language and means, “covered writing”. While the goal of cryptography is to make a message unreadable, steganography aims at avoiding suspicion to the existence of a hidden message. A coding scheme is a set of rules that determines which symbol of the source alphabet is represented by which symbol in the target alphabet. Several different coding schemes have been developed within the past decade by different groups of researchers. The coding schemes differ in codon length, detectability, number of characters that can be encoded, and the number of steps involved in encoding a message. Recent advances in genetic engineering have allowed the insertion of artificial DNA strands into the living cells of organisms. Instead of expressing a message as a series of ones and zeros, it is represented in DNA code as a series of As, Cs, Gs, and Ts. A number of algorithms have been developed to encode a message in DNA characters and either disguise these messages as novel DNA sequences or encapsulate them within existing ones. In this work we describe methods to insert information into a DNA sequence for the purpose of data storage, watermarking, or communication of secret messages.

Gene Prediction Based on DNA Spectral Analysis: A Literature Review

The identification of regions of DNA sequences that code for proteins is one of the most fundamental applications in bioinformatics. These protein-coding regions are in contrast to other DNA regions that encode functional RNA molecules, provide structural stability of chromosomes, serve as genetic raw materials, represent molecular fossils, or have no known purpose (sometimes called "junk DNA"). A number of approaches have been suggested for differentiating between the protein-coding and non-protein-coding regions of DNA. A selection of these approaches is based on digital signal processing (DSP) techniques. These DSP techniques rely on the phenomenon that protein-coding regions have a prominent power spectrum peak at frequency f=⅓ arising from the length of codons (three nucleic acids). This article partitions the identification of protein-coding regions into four discrete steps. Based on this partitioning, DSP techniques can be easily described and compared based on their unique implementations of the processing steps. We compare the approaches, and discuss strengths and weaknesses of each in the context of different applications. Our work provides an accessible introduction and comparative review of DSP methods for the identification of protein-coding regions. Additionally, by breaking down the approaches into four steps, we suggest new combinations that may be worthy of future study.

Codon Size Reduction as the Origin of the Triplet Genetic Code

The genetic code appears to be optimized in its robustness to missense errors and frameshift errors. In addition, the genetic code is near-optimal in terms of its ability to carry information in addition to the sequences of encoded proteins. As evolution has no foresight, optimality of the modern genetic code suggests that it evolved from less optimal code variants. The length of codons in the genetic code is also optimal, as three is the minimal nucleotide combination that can encode the twenty standard amino acids. The apparent impossibility of transitions between codon sizes in a discontinuous manner during evolution has resulted in an unbending view that the genetic code was always triplet. Yet, recent experimental evidence on quadruplet decoding, as well as the discovery of organisms with ambiguous and dual decoding, suggest that the possibility of the evolution of triplet decoding from living systems with non-triplet decoding merits reconsideration and further exploration. To explore this possibility we designed a mathematical model of the evolution of primitive digital coding systems which can decode nucleotide sequences into protein sequences. These coding systems can evolve their nucleotide sequences via genetic events of Darwinian evolution, such as point-mutations. The replication rates of such coding systems depend on the accuracy of the generated protein sequences. Computer simulations based on our model show that decoding systems with codons of length greater than three spontaneously evolve into predominantly triplet decoding systems. Our findings suggest a plausible scenario for the evolution of the triplet genetic code in a continuous manner. This scenario suggests an explanation of how protein synthesis could be accomplished by means of long RNA-RNA interactions prior to the emergence of the complex decoding machinery, such as the ribosome, that is required for stabilization and discrimination of otherwise weak triplet codon-anticodon interactions.

Predicting mutation frequencies in stem-loop structures of derepressed genes: implications for evolution.

This work provides evidence that, during transcription, the mutability (propensity to mutate) of a base in a DNA secondary structure depends both on the stability of the structure and on the extent to which the base is unpaired. Zuker's DNA folding computer program reveals the most probable stem-loop structures (SLSs) and negative energies of folding (-DeltaG) for any given nucleotide sequence. We developed an interfacing program that calculates (i) the percentage of folds in which each base is unpaired during transcription; and (ii) the mutability index (MI) for each base, expressed as an absolute value and defined as -follows: MI = (% total folds in which the base is unpaired) x (highest -DeltaG of all folds in which it is unpaired). Thus, MIs predict the relative mutation or reversion frequencies of unpaired bases in SLSs. MIs for 16 mutable bases in auxotrophs, selected during starvation in derepressed genes, are compared with 70 background mutations in lacI and ebgR that were not derepressed during mutant selection. All the results are consistent with the location of known mutable bases in SLSs. Specific conclusions are: (i) Of 16 mutable bases in transcribing genes, 87% have higher MIs than the average base of the sequence analysed, compared with 50% for the 70 background mutations. (ii) In 15 of the mutable bases of transcribing genes, the correlation between MIs and relative mutation frequencies determined experimentally is good. There is no correlation for 35 mutable bases in the lacI gene. (iii) In derepressed auxotrophs, 100% of the codons containing the mutable bases are within one codon's length of a stem, compared with 53% for the background mutable bases in lacI. (iv) The data suggest that environmental stressors may cause as well as select mutations in derepressed genes. The implications of these results for evolution are discussed.

Distances between Completely Matched Tuples in the Modern Gene Sequences

Modern genome sequences are the products as a consequence of gene evolution. When we consider the gene evolution, it is necessary to investigate a universal feature of base sequences of the modern genes. We have investigated the characteristics of base sequences of modern genes using statistical analysis. In our previous works, we revealed that the base sequences of the modern genes universally have statistically significant repetitive short tuples [1]. These significant repetitive tuples (SRTs) exist not only in coding region but also in noncoding region and the lengths of SRTs are not multiple of the codon length. Each gene has a specific set of the base sequences of SRTs in the whole gene consisting of the coding and noncoding regions [2]. A set of genes encoding proteins with different functions contains few common SRTs. On the other hand, a gene family has a number of common SRTs. It is possible to say that the base sequences of SRTs are influenced by the codons used frequently in the coding region. If the base sequences of the repetitive short tuples are related to the sequences of codons, our results reflect the preference of codon usage. We have not yet studied the actual location of SRTs in the reading frame. Instead, in this work, we analysed distances of completely matched tuples to investigate whether or not the location of tuples with the same base sequence in coding region is related to the reading frame.

Improving the Efficiency of the Genetic Code by Varying the Codon Length—The Perfect Genetic Code

The function of DNA is to specify protein sequences. The four-base “alphabet” used in nucleic acids is translated to the 20 base alphabet of proteins (plus a stop signal) via the genetic code. The code is neither overlapping nor punctuated, but has mRNA sequences read in successive triplet codons until reaching a stop codon. The true genetic code uses three bases for every amino acid. The efficiency of the genetic code can be significantly increased if the requirement for a fixed codon length is dropped so that the more common amino acids have shorter codon lengths and rare amino acids have longer codon lengths. More efficient codes can be derived using the Shannon–Fano and Huffman coding algorithms. The compression achieved using a Huffman code cannot be improved upon. I have used these algorithms to derive efficient codes for representing protein sequences using both two and four bases. The length of DNA required to specify the complete set of protein sequences could be significantly shorter if transcription used a variable codon length. The restriction to a fixed codon length of three bases means that it takes 42% more DNA than the minimum necessary, and the genetic code is 70% efficient. One can think of many reasons why this maximally efficient code has not evolved: there is very little redundancy so almost any mutation causes an amino acid change. Many mutations will be potentially lethal frame-shift mutations, if the mutation leads to a change in codon length. It would be more difficult for the machinery of transcription to cope with a variable codon length. Nevertheless, in the strict and narrow sense of coding for protein sequences using the minimum length of DNA possible, the Huffman code derived here is perfect.

A putative beta-glucanase pseudogene behind the potato GBSS gene.

We identified an open reading frame (ORF) which is located closely behind the gene encoding granule-bound starch synthase (GBSS) of potato (Solanum tuberosum L.). The ORF ends with a perfect 43 bp direct repeat, which carries the stop triplet precisely at the beginning of the second repeat. The deduced protein shows homology with all known isoforms of plant beta-1,3-glucanases and beta-1,3-1,4-glucanases. Although the DNA sequence is unique in potato and tomato (Lycopersicon esculentum L.), no expression of the gene was found in these species. Taken together with the unusual codon usage and length of the predicted protein, this sequence could be a pseudogene.

Codon recognition in polypeptide chain termination:site directed crosslinking of termination codon toEscherichia colirelease factor 2

An RNA synthesized in vitro was positioned on the Escherichia coli ribosome at the P site with tRNAala, and with a termination codon, UAA, as the next codon in the A site. Such a complex bound stoichiometric amounts of release factor 2 (RF-2); a corresponding RNA with UAC in place of UAA was not a template for the factor. An RNA containing 4-thio-UAA in place of the UAA supported binding of RF-2, and this has allowed site-directed crosslinking from the first position of the termination codon to answer two long standing questions about the termination of protein biosynthesis, the position of the termination codon and its proximity to the release factor during codon recognition. An RF-2.mRNA crosslinked product was detected, indicating the release factor and the termination codon are in close physical contact during the codon recognition event of termination. The 4-thio-U crosslinked also to the ribosome but only to the 30S subunit, and the proteins and the rRNA site concerned were identified. RF-2 decreased significantly the crosslinking to the ribosomal components, but no new crosslink sites were found. If the stop codon was deliberately displaced from the decoding site by one codon's length then a different pattern of crosslinking in particular to the rRNA resulted. These observations are consistent with a model of codon recognition by RF-2 at the decoding site, without a major shift in position of the codon.

Origin of noncoding DNA sequences: molecular fossils of genome evolution.

The total amount of noncoding sequences on chromosomes of contemporary organisms varies significantly from species to species. We propose a hypothesis for the origin of these noncoding sequences that assumes that (i) an approximately equal to 0.55-kilobase (kb)-long reading frame composed the primordial gene and (ii) a 20-kb-long single-stranded polynucleotide is the longest molecule (as a genome) that was polymerized at random and without a specific template in the primordial soup/cell. The statistical distribution of stop codons allows examination of the probability of generating reading frames of approximately equal to 0.55 kb in this primordial polynucleotide. This analysis reveals that with three stop codons, a run of at least 0.55-kb equivalent length of nonstop codons would occur in 4.6% of 20-kb-long polynucleotide molecules. We attempt to estimate the total amount of noncoding sequences that would be present on the chromosomes of contemporary species assuming that present-day chromosomes retain the prototype primordial genome structure. Theoretical estimates thus obtained for most eukaryotes do not differ significantly from those reported for these specific organisms, with only a few exceptions. Furthermore, analysis of possible stop-codon distributions suggests that life on earth would not exist, at least in its present form, had two or four stop codons been selected early in evolution.

Improved direct expression of prochymosin cDNA through changing the SD-ATG codon length.

Improved direct expression of prochymosin cDNA through changing the SD-ATG codon length.