Genomic analysis and phylogenetic characterization of Himalayan snow trout, Schizothorax esocinus based on mitochondrial protein-coding genes.

BACKGROUND: Mitochondrial DNA (mtDNA) has become a significant tool for exploring genetic diversity and delineating evolutionary links across diverse taxa. Within the group of cold-water fish species that are native to the Indian Himalayan region, Schizothorax esocinus holds particular importance due to its ecological significance and is potentially vulnerable to environmental changes. This research aims to clarify the phylogenetic relationships within the Schizothorax genus by utilizing mitochondrial protein-coding genes. METHODS: Standard protocols were followed for the isolation of DNA from S. esocinus. For the amplification of mtDNA, overlapping primers were used, and then subsequent sequencing was performed. The genetic features were investigated by the application of bioinformatic approaches. These approaches covered the evaluation of nucleotide composition, codon usage, selective pressure using nonsynonymous substitution /synonymous substitution (Ka/Ks) ratios, and phylogenetic analysis. RESULTS: The study specifically examined the 13 protein-coding genes of Schizothorax species which belongs to the Schizothoracinae subfamily. Nucleotide composition analysis showed a bias towards A + T content, consistent with other cyprinid fish species, suggesting evolutionary conservation. Relative Synonymous Codon Usage highlighted leucine as the most frequent (5.18%) and cysteine as the least frequent (0.78%) codon. The positive AT-skew and the predominantly negative GC-skew indicated the abundance of A and C. Comparative analysis revealed significant conservation of amino acids in multiple genes. The majority of amino acids were hydrophobic rather than polar. The purifying selection was revealed by the genetic distance and Ka/Ks ratios. Phylogenetic study revealed a significant genetic divergence between S. esocinus and other Schizothorax species with interspecific K2P distances ranging from 0.00 to 8.87%, with an average of 5.76%. CONCLUSION: The present study provides significant contributions to the understanding of mitochondrial genome diversity and genetic evolution mechanisms in Schizothoracinae, hence offering vital insights for the development of conservation initiatives aimed at protecting freshwater fish species.

Analysis of chloroplast genome characteristics and codon usage bias in 14 species of Annonaceae.

For the study of species evolution, chloroplast gene expression, and transformation, the chloroplast genome is an invaluable resource. Codon usage bias (CUB) analysis is a tool that is utilized to improve gene expression and investigate evolutionary connections in genetic transformation. In this study, we analysed chloroplast genome differences, codon usage patterns and the sources of variation on CUB in 14 Annonaceae species using bioinformatics tools. The study showed that there was a significant variation in both gene sizes and numbers between the 14 species, but conservation was still maintained. It's worth noting that there were noticeable differences in the IR/SC sector boundary and the types of SSRs among the 14 species. The mono-nucleotide repeat type was the most common, with A/T repeats being more prevalent than G/C repeats. Among the different types of repeats, forward and palindromic repeats were the most abundant, followed by reverse repeats, and complement repeats were relatively rare. Codon composition analysis revealed that all 14 species had a frequency of GC lower than 50%. Additionally, it was observed that the proteins in-coding sequences of chloroplast genes tend to end with A/T at the third codon position. Among these species, 21 codons exhibited bias (RSCU > 1), and there were 8 high-frequency (HF) codons and 5 optimal codons that were identical across the species. According to the ENC-plot and Neutrality plot analysis, natural selection had less impact on the CUB of A. muricate and A. reticulata. Based on the PR2-plot, it was evident that base G had a higher frequency than C, and T had a higher frequency A. The correspondence analysis (COA) revealed that codon usage patterns different in Annonaceae.

Decoding Evolution of Rubioideae: Plastomes Reveal Sweet Secrets of Codon Usage, Diagnostides, and Superbarcoding.

Galium genus belongs to the Rubiaceae family, which consists of approximately 14,000 species. In comparison to its well-known relatives, the plastomes of the Galium genus have not been explored so far. The plastomes of this genus have a typical, quadripartite structure, but differ in gene content, since the infA gene is missing in Galium palustre and Galium trfidum. An evaluation of the effectiveness of using entire chloroplast genome sequences as superbarcodes for accurate plant species identification revealed the high potential of this method for molecular delimitation within the genus and tribe. The trnE-UUC-psbD region showed the biggest number of diagnostides (diagnostic nucleotides) which might be new potential barcodes, not only in Galium, but also in other closely related genera. Relative synonymous codon usage (RSCU) appeared to be connected with the phylogeny of the Rubiaceae family, showing that during evolution, plants started preferring specific codons over others.

An important consideration when expressing mAbs in Escherichiacoli.

Monoclonal antibodies (mAbs) are a driving force in the biopharmaceutical industry. Therapeutic mAbs are usually produced in mammalian cells, but there has been a push towards the use of alternative production hosts, such as Escherichia coli. When the genes encoding for a mAb heavy and light chains are codon-optimized for E. coli expression, a truncated form of the heavy chain can form along with the full-length product. In this work, the role of codon optimization in the formation of a truncated product was investigated. This study used the amino acid sequences of several therapeutic mAbs and multiple optimization algorithms. It was found that several algorithms incorporate sequences that lead to a truncated product. Approaches to avoid this truncated form are discussed.

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

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

Robust genetic codes enhance protein evolvability.

The standard genetic code defines the rules of translation for nearly every life form on Earth. It also determines the amino acid changes accessible via single-nucleotide mutations, thus influencing protein evolvability-the ability of mutation to bring forth adaptive variation in protein function. One of the most striking features of the standard genetic code is its robustness to mutation, yet it remains an open question whether such robustness facilitates or frustrates protein evolvability. To answer this question, we use data from massively parallel sequence-to-function assays to construct and analyze 6 empirical adaptive landscapes under hundreds of thousands of rewired genetic codes, including those of codon compression schemes relevant to protein engineering and synthetic biology. We find that robust genetic codes tend to enhance protein evolvability by rendering smooth adaptive landscapes with few peaks, which are readily accessible from throughout sequence space. However, the standard genetic code is rarely exceptional in this regard, because many alternative codes render smoother landscapes than the standard code. By constructing low-dimensional visualizations of these landscapes, which each comprise more than 16 million mRNA sequences, we show that such alternative codes radically alter the topological features of the network of high-fitness genotypes. Whereas the genetic codes that optimize evolvability depend to some extent on the detailed relationship between amino acid sequence and protein function, we also uncover general design principles for engineering nonstandard genetic codes for enhanced and diminished evolvability, which may facilitate directed protein evolution experiments and the bio-containment of synthetic organisms, respectively.

Massively parallel CRISPR-assisted homologous recombination enables saturation editing of full-length endogenous genes in yeast.

Performing saturation editing of chromosomal genes will enable the study of genetic variants in situ and facilitate protein and cell engineering. However, current in vivo editing of endogenous genes either lacks flexibility or is limited to discrete codons and short gene fragments, preventing a comprehensive exploration of genotype-phenotype relationships. To enable facile saturation editing of full-length genes, we used a protospacer adjacent motif-relaxed Cas9 variant and homology-directed repair to achieve above 60% user-defined codon replacement efficiencies in Saccharomyces cerevisiae genome. Coupled with massively parallel DNA design and synthesis, we developed a saturation gene editing method termed CRISPR-Cas9- and homology-directed repair-assisted saturation editing (CHASE) and achieved highly saturated codon swapping of long genomic regions. By applying CHASE to massively edit a well-studied global transcription factor gene, we found known and unreported genetic variants affecting an industrially relevant microbial trait. The user-defined codon editing capability and wide targeting windows of CHASE substantially expand the scope of saturation gene editing.

In-silico Codon Context and Synonymous Usage Analysis of Genes for Molecular Mechanisms Inducing Autophagy and Apoptosis with Reference to Neurodegenerative Disorders.

Background: Autophagy and apoptosis are cellular processes that maintain cellular homeostasis and remove damaged or aged organelles or aggregated and misfolded proteins. Stress factors initiate the signaling pathways common to autophagy and apoptosis. An imbalance in the autophagy and apoptosis, led by cascade of molecular mechanism prior to both processes culminate into neurodegeneration. Objective: In present study, we urge to investigate the codon usage pattern of genes which are common before initiating autophagy and apoptosis. Methods: In the present study, we took up eleven genes (DAPK1, BECN1, PIK3C3 (VPS34), BCL2, MAPK8, BNIP3 L (NIX), PMAIP1, BAD, BID, BBC3, MCL1) that are part of molecular signaling mechanism prior to autophagy and apoptosis. We analyzed dinucleotide odds ratio, codon bias, usage, context, and rare codon analysis. Results: CpC and GpG dinucleotides were abundant, with the dominance of G/C ending codons as preferred codons. Clustering analysis revealed that MAPK8 had a distinct codon usage pattern compared to other envisaged genes. Both positive and negative contexts were observed, and GAG-GAG followed by CTG-GCC was the most abundant codon pair. Of the six synonymous arginine codons, two codons CGT and CGA were the rarest. Conclusions: The information presented in the study may be used to manipulate the process of autophagy and apoptosis and to check the pathophysiology associated with their dysregulation.

Comprehensive analysis of two hotspot codons in the TUBB4B gene and associated phenotypes.

Our purpose was to elucidate the genotype and ophthalmological and audiological phenotype in TUBB4B-associated inherited retinal dystrophy (IRD) and sensorineural hearing loss (SNHL), and to model the effects of all possible amino acid substitutions at the hotspot codons Arg390 and Arg391. Six patients from five families with heterozygous missense variants in TUBB4B were included in this observational study. Ophthalmological testing included best-corrected visual acuity, fundus examination, optical coherence tomography, fundus autofluorescence imaging, and full-field electroretinography (ERG). Audiological examination included pure-tone and speech audiometry in adult patients and auditory brainstem response testing in a child. Genetic testing was performed by disease gene panel analysis based on genome sequencing. The molecular consequences of the substitutions of residues 390 and 391 on TUBB4B and its interaction with α-tubulin were predicted in silico on its three-dimensional structure obtained by homology modelling. Two independent patients had amino acid exchanges at position 391 (p.(Arg391His) or p.(Arg391Cys)) of the TUBB4B protein. Both had a distinct IRD phenotype with peripheral round yellowish lesions with pigmented spots and mild or moderate SNHL, respectively. Yet the phenotype was milder with a sectorial pattern of bone spicules in one patient, likely due to a genetically confirmed mosaicism for p.(Arg391His). Three patients were heterozygous for an amino acid exchange at position 390 (p.(Arg390Gln) or p.(Arg390Trp)) and presented with another distinct retinal phenotype with well demarcated pericentral retinitis pigmentosa. All showed SNHL ranging from mild to severe. One additional patient showed a variant distinct from codon 390 or 391 (p.(Tyr310His)), and presented with congenital profound hearing loss and reduced responses in ERG. Variants at codon positions 390 and 391 were predicted to decrease the structural stability of TUBB4B and its complex with α-tubulin, as well as the complex affinity. In conclusion, the twofold larger reduction in heterodimer affinity exhibited by Arg391 substitutions suggested an association with the more severe retinal phenotype, compared to the substitution at Arg390.

Intraspecific and interspecific variations in the synonymous codon usage in mitochondrial genomes of 8 pleurotus strains.

In this study, we investigated the codon bias of twelve mitochondrial core protein coding genes (PCGs) in eight Pleurotus strains, two of which are from the same species. The results revealed that the codons of all Pleurotus strains had a preference for ending in A/T. Furthermore, the correlation between codon base compositions and codon adaptation index (CAI), codon bias index (CBI) and frequency of optimal codons (FOP) indices was also detected, implying the influence of base composition on codon bias. The two P. ostreatus species were found to have differences in various base bias indicators. The average effective number of codons (ENC) of mitochondrial core PCGs of Pleurotus was found to be less than 35, indicating strong codon preference of mitochondrial core PCGs of Pleurotus. The neutrality plot analysis and PR2-Bias plot analysis further suggested that natural selection plays an important role in Pleurotus codon bias. Additionally, six to ten optimal codons (ΔRSCU > 0.08 and RSCU > 1) were identified in eight Pleurotus strains, with UGU and ACU being the most widely used optimal codons in Pleurotus. Finally, based on the combined mitochondrial sequence and RSCU value, the genetic relationship between different Pleurotus strains was deduced, showing large variations between them. This research has improved our understanding of synonymous codon usage characteristics and evolution of this important fungal group.

Generation of genetic codes with 2-adic codon algebra and adaptive dynamics.

This is a brief review on modeling genetic codes with the aid of 2-adic dynamical systems. In this model amino acids are encoded by the attractors of such dynamical systems. Each genetic code is coupled to the special class of 2-adic dynamics. We consider the discrete dynamical systems, These are the iterations of a function F:Z2→Z2, where Z2 is the ring of 2-adic numbers (2-adic tree). A genetic code is characterized by the set of attractors of a function belonging to the code generating functional class. The main mathematical problem is to reduce degeneration of dynamic representation and select the optimal generating function. Here optimality can be treated in many ways. One possibility is to consider the Lipschitz functions playing the crucial role in general theory of iterations. Then we minimize the Lip-constant. The main issue is to find the proper biological interpretation of code-functions. One can speculate that the evolution of the genetic codes can be described in information space of the nucleotide-strings endowed with ultrametric (treelike) geometry. A code-function is a fitness function; the solutions of the genetic code optimization problem are attractors of the code-function. We illustrate this approach by generation of the standard nuclear and (vertebrate) mitochondrial genetics codes.

Comparative analysis of codon usage bias in chloroplast genomes of ten medicinal species of Rutaceae.

Rutaceae family comprises economically important plants due to their extensive applications in spices, food, oil, medicine, etc. The Rutaceae plants is able to better utilization through biotechnology. Modern biotechnological approaches primarily rely on the heterologous expression of functional proteins in different vectors. However, several proteins are difficult to express outside their native environment. The expression potential of functional genes in heterologous systems can be maximized by replacing the rare synonymous codons in the vector with preferred optimal codons of functional genes. Codon usage bias plays a critical role in biogenetic engineering-based research and development. In the current study, 727 coding sequences (CDSs) obtained from the chloroplast genomes of ten Rutaceae plant family members were analyzed for codon usage bias. The nucleotide composition analysis of codons showed that these codons were rich in A/T(U) bases and preferred A/T(U) endings. Analyses of neutrality plots, effective number of codons (ENC) plots, and correlations between ENC and codon adaptation index (CAI) were conducted, which revealed that natural selection is a major driving force for the Rutaceae plant family's codon usage bias, followed by base mutation. In the ENC vs. CAI plot, codon usage bias in the Rutaceae family had a negligible relationship with gene expression level. For each sample, we screened 12 codons as preferred and high-frequency codons simultaneously, of which GCU encoding Ala, UUA encoding Leu, and AGA encoding Arg were the most preferred codons. Taken together, our study unraveled the synonymous codon usage pattern in the Rutaceae family, providing valuable information for the genetic engineering of Rutaceae plant species in the future.

A new look at the fluorescent protein-based approach for identifying optimal coding sequence for recombinant protein expression in E. coli.

Due to the degeneracy of the genetic code, most amino acids are encoded by several codons. The choice among synonymous codons at the N-terminus of genes has a profound effect on protein expression in Escherichia coli. This is often explained by the different contributions of synonymous codons to mRNA secondary structure formation. Strong secondary structures at the 5'-end of mRNA interfere with ribosome binding and affect the process of translation initiation. In silico optimization of the gene 5'-end can significantly increase the level of protein expression; however, this method is not always effective due to the uncertainty of the exact mechanism by which synonymous substitutions affect expression; thus, it may produce nonoptimal variants as well as miss some of the best producers. In this paper, an alternative approach is proposed based on screening a partially randomized library of expression constructs comprising hundreds of selected synonymous variants. The effect of such substitutions was evaluated using the gene of interest fused to the reporter gene of the fluorescent protein with subsequent screening for the most promising candidates according to the reporter's signal intensity. The power of the approach is demonstrated by a significant increase in the prokaryotic expression of three proteins: canine cystatin C, human BCL2-associated athanogene 3 and human cardiac troponin I. This simple approach was suggested which may provide an efficient, easy, and inexpensive optimization method for poorly expressed proteins in bacteria.

Translation velocity determines the efficacy of engineered suppressor tRNAs on pathogenic nonsense mutations.

Nonsense mutations - the underlying cause of approximately 11% of all genetic diseases - prematurely terminate protein synthesis by mutating a sense codon to a premature stop or termination codon (PTC). An emerging therapeutic strategy to suppress nonsense defects is to engineer sense-codon decoding tRNAs to readthrough and restore translation at PTCs. However, the readthrough efficiency of the engineered suppressor tRNAs (sup-tRNAs) largely varies in a tissue- and sequence context-dependent manner and has not yet yielded optimal clinical efficacy for many nonsense mutations. Here, we systematically analyze the suppression efficacy at various pathogenic nonsense mutations. We discover that the translation velocity of the sequence upstream of PTCs modulates the sup-tRNA readthrough efficacy. The PTCs most refractory to suppression are embedded in a sequence context translated with an abrupt reversal of the translation speed leading to ribosomal collisions. Moreover, modeling translation velocity using Ribo-seq data can accurately predict the suppression efficacy at PTCs. These results reveal previously unknown molecular signatures contributing to genotype-phenotype relationships and treatment-response heterogeneity, and provide the framework for the development of personalized tRNA-based gene therapies.

Theories of the origin of the genetic code: Strong corroboration for the coevolution theory.

I analyzed all the theories and models of the origin of the genetic code, and over the years, I have considered the main suggestions that could explain this origin. The conclusion of this analysis is that the coevolution theory of the origin of the genetic code is the theory that best captures the majority of observations concerning the organization of the genetic code. In other words, the biosynthetic relationships between amino acids would have heavily influenced the origin of the organization of the genetic code, as supported by the coevolution theory. Instead, the presence in the genetic code of physicochemical properties of amino acids, which have also been linked to the physicochemical properties of anticodons or codons or bases by stereochemical and physicochemical theories, would simply be the result of natural selection. More explicitly, I maintain that these correlations between codons, anticodons or bases and amino acids are in fact the result not of a real correlation between amino acids and codons, for example, but are only the effect of the intervention of natural selection. Specifically, in the genetic code table we expect, for example, that the most similar codons - that is, those that differ by only one base - will have more similar physicochemical properties. Therefore, the 64 codons of the genetic code table ordered in a certain way would also represent an ordering of some of their physicochemical properties. Now, a study aimed at clarifying which physicochemical property of amino acids has influenced the allocation of amino acids in the genetic code has established that the partition energy of amino acids has played a role decisive in this. Indeed, under some conditions, the genetic code was found to be approximately 98% optimized on its columns. In this same work, it was shown that this was most likely the result of the action of natural selection. If natural selection had truly allocated the amino acids in the genetic code in such a way that similar amino acids also have similar codons - this, not through a mechanism of physicochemical interaction between, for example, codons and amino acids - then it might turn out that even different physicochemical properties of codons (or anticodons or bases) show some correlation with the physicochemical properties of amino acids, simply because the partition energy of amino acids is correlated with other physicochemical properties of amino acids. It is very likely that this would inevitably lead to a correlation between codons (or anticodons or bases) and amino acids. In other words, since the codons (anticodons or bases) are ordered in the genetic code, that is to say, some of their physicochemical properties should also be ordered by a similar order, and given that the amino acids would also appear to have been ordered in the genetic code by selection natural, then it should inevitably turn out that there is a correlation between, for example, the hydrophobicity of anticodons and that of amino acids. Instead, the intervention of natural selection in organizing the genetic code would appear to be highly compatible with the main mechanism of structuring the genetic code as supported by the coevolution theory. This would make the coevolution theory the only plausible explanation for the origin of the genetic code.

Elucidating the multichromosomal structure within the Brasenia schreberi mitochondrial genome through assembly and analysis.

Brasenia schreberi, a plant species traditionally utilized in Chinese medicine and cuisine, represents an early evolutionary stage among flowering plants (angiosperms). While the plastid genome of this species has been published, its mitochondrial genome (mitogenome) has not been extensively explored, with a notable absence of thorough comparative analyses of its organellar genomes. In our study, we had assembled the entire mitogenome of B. schreberi utilizing the sequencing data derived from both Illumina platform and Oxford Nanopore. The B. schreberi mitogenome mostly exists as six circular DNA molecules, with the largest being 628,257 base pairs (bp) and the smallest 110,220 bp, amounting to 1.49 megabases (Mb). Then we annotated the mitogenome of B. schreberi. The mitogenome encompasses a total of 71 genes: 40 of these are coding proteins genes (PCGs), 28 are genes for transfer RNA (tRNA), and the remaining 3 are genes for ribosomal RNA (rRNA). In the analysis of codon usage, we noted a unique codon preference specific to each amino acid. The most commonly used codons exhibited an average RSCU of 1.36, indicating a noticeable bias in codon selection. In the repeat sequence analysis, a total of 553 simple sequence repeats (SSRs) were identified, 1,822 dispersed repeats (comprising 1,015 forward and 807 palindromic repeats), and 608 long terminal repeats (LTRs). Additionally, in the analysis of homologous sequences between organelle genomes, we detected 38 homologous sequences derived from the plastid genome, each exceeding 500 bp, within the B. schreberi mitochondrial genome. Notably, ten tRNA genes (trnC-GCA, trnM-CAU, trnI-CAU, trnQ-UUG, trnN-GUU, trnT-GGU, trnW-CCA, trnA-UGC, trnI-GAU, and trnV-GAC) appear to have been completely transferred from the chloroplast to the mitogenome. Utilizing the Deepred-mt to predict the RNA editing sites in the mitogenome, we have identified 675 high-quality RNA editing sites in the 40 mitochondrial PCGs. In the final stage of our study, we performed an analysis of colinearity and inferred the phylogenetic relationship of B. schreberi with other angiosperms, utilizing the mitochondrial PCGs as a basis. The results showed that the non-coding regions of the B. schreberi mitogenome are characterized by an abundance of repetitive sequences and exogenous sequences, and B. schreberi is more closely related with Euryale ferox.

The calcium-binding photoprotein clytin II: Expression of the preferred human codon-optimized clytin II gene in Chinese hamster ovary-K1 cells and its use in the G-protein-coupled receptor assays.

Clytin II (CLII) is a Ca2+-binding photoprotein and has been identified as an isotype of clytin I (CLI). CLII consists of apoCLII (an apoprotein) and 2-peroxide of coelenterazine (an adduct of molecular oxygen to coelenterazine), which is identical to the widely used Ca2+-binding photoprotein, aequorin (AQ). However, CLII triggered by Ca2+ exhibits a 4.5-fold higher maximum luminescence intensity (Imax) compared to both AQ and CLI, and it is approximately 5 times less sensitive to Ca2+ than AQ. To confirm the suitability of the preferred human codon-optimized CLII (pCLII) gene for cell-based G-protein-coupled receptor (GPCR) assays, a transformant stably expressing apoprotein of pCLII using the pCLII gene in the mitochondria of CHO-K1 cells was established and in situ regenerated pCLII in the cells were applied to the high-throughput screening system. An ATP-stimulated GPCR assay for endogenous P2Y purinergic receptors was confirmed using the established stable transformant.

Codon usage bias analysis in the mitochondrial genomes of five Rhingia Scopoli (Diptera, Syrphidae, Eristalinae) species.

This study presents the sequencing and annotation of mitochondrial genomes from five Rhingia species of the family Syrphidae, focusing on codon bias. Each species possessed 22 tRNAs genes, 13 protein-coding genes, 2 rRNAs genes, and a control region, without any observed gene rearrangements. Nucleotide composition analysis revealed a higher AT content compared with GC content, indicating AT enrichment. Neutrality plot, Parity rule 2 bias, and effective number of codons plot analyses collectively indicated that natural selection primarily influences the codon usage bias in the five Rhingia species. Relative synonymous codon usage analysis identified the optimal codons for Rhingia binotata, R. fromosana, R. campestris, R. louguanensis, and R. xanthopoda as 10, 14, 10, 11, and 12, respectively, all ending with A/U and exhibiting AT preference. Phylogenetic analysis, based on maximum likelihood and Bayesian inference methods applied to three datasets, confirmed the monophyly of Rhingia. In conclusion, this research establishes a foundation for understanding the phylogenetic evolution and codon usage patterns in Rhingia, offering valuable for future studies.

Initial Amino Acid:Codon Assignments and Strength of Codon:Anticodon Binding.

The ribosome brings 3'-aminoacyl-tRNA and 3'-peptidyl-tRNAs together to enable peptidyl transfer by binding them in two major ways. First, their anticodon loops are bound to mRNA, itself anchored at the ribosomal subunit interface, by contiguous anticodon:codon pairing augmented by interactions with the decoding center of the small ribosomal subunit. Second, their acceptor stems are bound by the peptidyl transferase center, which aligns the 3'-aminoacyl- and 3'-peptidyl-termini for optimal interaction of the nucleophilic amino group and electrophilic ester carbonyl group. Reasoning that intrinsic codon:anticodon binding might have been a major contributor to bringing tRNA 3'-termini into proximity at an early stage of ribosomal peptide synthesis, we wondered if primordial amino acids might have been assigned to those codons that bind the corresponding anticodon loops most tightly. By measuring the binding of anticodon stem loops to short oligonucleotides, we determined that family-box codon:anticodon pairings are typically tighter than split-box codon:anticodon pairings. Furthermore, we find that two family-box anticodon stem loops can tightly bind a pair of contiguous codons simultaneously, whereas two split-box anticodon stem loops cannot. The amino acids assigned to family boxes correspond to those accessible by what has been termed cyanosulfidic chemistry, supporting the contention that these limited amino acids might have been the first used in primordial coded peptide synthesis.

An Integer Linear Programming Model to Optimize Coding DNA Sequences By Joint Control of Transcript Indicators.

A Coding DNA Sequence (CDS) is a fraction of DNA whose nucleotides are grouped into consecutive triplets called codons, each one encoding an amino acid. Because most amino acids can be encoded by more than one codon, the same amino acid chain can be obtained by a very large number of different CDSs. These synonymous CDSs show different features that, also depending on the organism the transcript is expressed in, could affect translational efficiency and yield. The identification of optimal CDSs with respect to given transcript indicators is in general a challenging task, but it has been observed in recent literature that integer linear programming (ILP) can be a very flexible and efficient way to achieve it. In this article, we add evidence to this observation by proposing a new ILP model that simultaneously optimizes different well-grounded indicators. With this model, we efficiently find solutions that dominate those returned by six existing codon optimization heuristics.