Genes enriched in A/T-ending codons are co-regulated and conserved across mammals

Uncovering codon usage patterns during murine embryogenesis and tissue-specific developmental diseases

Introduction

Mouse models share significant genetic similarities with humans and have expanded our understanding of how embryonic tissue-specific genes influence disease states. By improved analyses of temporal, transcriptional data from these models, we can capture unique tissue codon usage patterns and determine how deviations from these patterns can influence developmental disorders.

Methods

We analyzed transcriptomic-weighted data from four mouse strains across three different germ layer tissues (liver, heart, and eye) and through embryonic stages. Applying a multifaceted approach, we calculated relative synonymous codon usage, reduced the dimensionality, and employed machine learning clustering techniques.

Results and discussion

These techniques identified relative synonymous codon usage differences/similarities among strains and deviations in codon usage patterns between healthy and disease-linked genes. Original transcriptomic mouse data and RefSeq gene sequences can be found at the associated Mouse Embryo CoCoPUTs (codon and codon pair usage tables) website. Future studies can leverage this resource to uncover further insights into the dynamics of embryonic development and the corresponding codon usage biases that are paramount to understanding disease processes of embryologic origin.

Genes and Pathways Comprising the Human and Mouse ORFeomes Display Distinct Codon Bias Signatures that Can Regulate Protein Levels

Arginine, glutamic acid and selenocysteine based codon bias has been shown to regulate the translation of specific mRNAs for proteins that participate in stress responses, cell cycle and transcriptional regulation. Defining codon-bias in gene networks has the potential to identify other pathways under translational control. Here we have used computational methods to analyze the ORFeome of all unique human (19,711) and mouse (22,138) open-reading frames (ORFs) to characterize codon-usage and codon-bias in genes and biological processes. We show that ORFeome-wide clustering of gene-specific codon frequency data can be used to identify ontology-enriched biological processes and gene networks, with developmental and immunological programs well represented for both humans and mice. We developed codon over-use ontology mapping and hierarchical clustering to identify multi-codon bias signatures in human and mouse genes linked to signaling, development, mitochondria and metabolism, among others. The most distinct multi-codon bias signatures were identified in human genes linked to skin development and RNA metabolism, and in mouse genes linked to olfactory transduction and ribosome, highlighting species-specific pathways potentially regulated by translation. Extreme codon bias was identified in genes that included transcription factors and histone variants. We show that re-engineering extreme usage of C- or U-ending codons for aspartic acid, asparagine, histidine and tyrosine in the transcription factors CEBPB and MIER1, respectively, significantly regulates protein levels. Our study highlights that multi-codon bias signatures can be linked to specific biological pathways and that extreme codon bias with regulatory potential exists in transcription factors for immune response and development.

Queuosine tRNA Modification: Connecting the Microbiome to the Translatome

Transfer RNA (tRNA) modifications play an important role in regulating mRNA translation at the codon level. tRNA modifications can influence codon selection and optimality, thus shifting translation toward specific sets of mRNAs in a dynamic manner. Queuosine (Q) is a tRNA modification occurring at the wobble position. In eukaryotes, queuosine is synthesized by the tRNA-guanine trans-glycosylase (TGT) complex, which incorporates the nucleobase queuine (or Qbase) into guanine of the GUN anticodons. Queuine is sourced from gut bacteria and dietary intake. Q was recently shown to be critical for cellular responses to oxidative and mitochondrial stresses, as well as its potential role in neurodegenerative diseases and brain health. These unique features of Q provide an interesting insight into the regulation of mRNA translation by gut bacteria, and the potential health implications. In this review, Q biology is examined in the light of recent literature and nearly 4 decades of research. Q's role in neuropsychiatric diseases and cancer is highlighted and discussed. Given the recent interest in Q, and the new findings, more research is needed to fully comprehend its biological function and disease relevance, especially in neurobiology.

Decoding Codon Bias: The Role of tRNA Modifications in Tissue-Specific Translation

The tRNA epitranscriptome has been recognized as an important player in mRNA translation regulation. Our knowledge of the role of the tRNA epitranscriptome in fine-tuning translation via codon decoding at tissue or cell levels remains incomplete. We analyzed tRNA expression and modifications as well as codon optimality across seven mouse tissues. Our analysis revealed distinct enrichment patterns of tRNA modifications in different tissues. Queuosine (Q) tRNA modification was most enriched in the brain compared to other tissues, while mitochondrial tRNA modifications and tRNA expression were highest in the heart. Using this observation, we synthesized, and delivered in vivo, codon-mutated EGFP for Q-codons, where the C-ending Q-codons were replaced with U-ending codons. The protein levels of mutant EGFP were downregulated in liver, which is poor in Q, while in brain EGFP, levels did not change. These data show that understanding tRNA modification enrichments across tissues is not only essential for understanding codon decoding and bias but can also be utilized for optimizing gene and mRNA therapeutics to be more tissue-, cell-, or condition-specific.

Decoding the general role of tRNA queuosine modification in eukaryotes

Transfer RNA (tRNA) contains modified nucleosides essential for modulating protein translation. One of these modifications is queuosine (Q), which affects NAU codons translation rate. For decades, multiple studies have reported a wide variety of species-specific Q-related phenotypes in different eukaryotes, hindering the identification of a general underlying mechanism behind that phenotypic diversity. Here, through bioinformatics analysis of representative eukaryotic genomes we have predicted: i) the genes enriched in NAU codons, whose translation would be affected by tRNA Q-modification (Q-genes); and ii) the specific biological processes of each organism enriched in Q-genes, which generally in eukaryotes would be related to ubiquitination, phosphatidylinositol metabolism, splicing, DNA repair or cell cycle. These bioinformatics results provide evidence to support for the first time in eukaryotes that the wide diversity of phenotypes associated with tRNA Q-modification previously described in various species would directly depend on the control of Q-genes translation, and would allow prediction of unknown Q-dependent processes, such as Akt activation and p53 expression, which we have tested in human cancer cells. Considering the relevance of the Q-related processes, our findings may support further exploration of the role of Q in cancer and other pathologies. Moreover, since eukaryotes must salvage Q from bacteria, we suggest that changes in Q supply by the microbiome would affect the expression of host Q-genes, altering its physiology.

Dynamic mRNA Stability Buffer Transcriptional Activation During Neuronal Differentiation and Is Regulated by SAMD4A

Neurons are exceptionally sensitive to oxidative stress, which is the basis for many neurodegenerative disease pathophysiologies. The posttranscriptional basis for neuronal differentiation and behavior is not well characterized. The steady-state levels of mRNA are outcomes of an interplay between RNA transcription and decay. However, the correlation between mRNA transcription, translation, and stability remains elusive. We utilized a SH-SY5Y-based neural differentiation model that is widely used to study neurodegenerative diseases. After neuronal differentiation, we observed enhanced sensitivity of mature neurons to mitochondrial stresses and ferroptosis induction. We employed a newly developed simplified mRNA stability profiling technique to explore the role of mRNA stability in SH-SY5Y neuronal differentiation model. Transcriptome-wide mRNA stability analysis revealed neural-specific RNA stability kinetics. Our analysis revealed that mRNA stability could either exert the buffering effect on gene products or change in the same direction as transcription. Importantly, we observed that changes in mRNA stability corrected over or under transcription of mRNAs to maintain mRNA translation dynamics. Furthermore, we conducted integrative analysis of our mRNA stability data set, and a published CRISPR-i screen focused on neuronal oxidative stress responses. Our analysis unveiled novel neuronal stress response genes that were not evident at the transcriptional or translational levels. SEPHS2 emerged as an important neuronal stress regulator based on this integrative analysis. Motif analysis unveiled SAMD4A as a major regulator of the dynamic changes in mRNA stability observed during differentiation. Knockdown of SAMD4A impaired neuronal differentiation and influenced the response to oxidative stress. Mechanistically, SAMD4A was found to alter the stability of several mRNAs. The novel insights into the interplay between mRNA stability and cellular behaviors provide a foundation for understanding neurodevelopmental processes and neurodegenerative disorders and highlight dynamic mRNA stability as an important layer of gene expression.

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

Background

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.

Methods

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.

Results

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.

Conclusion

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

BACKGROUND

Apicomplexan 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.

METHODS

Codon 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.

RESULTS

The 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.

CONCLUSIONS

In 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.

Subfunctionalisation of paralogous genes and evolution of differential codon usage preferences: The showcase of polypyrimidine tract binding proteins

Gene paralogs are copies of an ancestral gene that appear after gene or full genome duplication. When two sister gene copies are maintained in the genome, redundancy may release certain evolutionary pressures, allowing one of them to access novel functions. Here, we focused our study on gene paralogs on the evolutionary history of the three polypyrimidine tract binding protein genes (PTBP) and their concurrent evolution of differential codon usage preferences (CUPrefs) in vertebrate species. PTBP1-3 show high identity at the amino acid level (up to 80%) but display strongly different nucleotide composition, divergent CUPrefs and, in humans and in many other vertebrates, distinct tissue-specific expression levels. Our phylogenetic inference results show that the duplication events leading to the three extant PTBP1-3 lineages predate the basal diversification within vertebrates, and genomic context analysis illustrates that local synteny has been well preserved over time for the three paralogs. We identify a distinct evolutionary pattern towards GC3-enriching substitutions in PTBP1, concurrent with enrichment in frequently used codons and with a tissue-wide expression. In contrast, PTBP2s are enriched in AT-ending, rare codons, and display tissue-restricted expression. As a result of this substitution trend, CUPrefs sharply differ between mammalian PTBP1s and the rest of PTBPs. Genomic context analysis suggests that GC3-rich nucleotide composition in PTBP1s is driven by local substitution processes, while the evidence in this direction is thinner for PTBP2-3. An actual lack of co-variation between the observed GC composition of PTBP2-3 and that of the surrounding non-coding genomic environment would raise an interrogation on the origin of CUPrefs, warranting further research on a putative tissue-specific translational selection. Finally, we communicate an intriguing trend for the use of the UUG-Leu codon, which matches the trends of AT-ending codons. Our results are compatible with a scenario in which a combination of directional mutation-selection processes would have differentially shaped CUPrefs of PTBPs in vertebrates: the observed GC-enrichment of PTBP1 in placental mammals may be linked to genomic location and to the strong and broad tissue-expression, while AT-enrichment of PTBP2 and PTBP3 would be associated with rare CUPrefs and thus, possibly to specialized spatio-temporal expression. Our interpretation is coherent with a gene subfunctionalisation process by differential expression regulation associated with the evolution of specific CUPrefs.