Codon optimality, bias and usage in translation and mRNA decay

YMLA: A comparative platform to carry out functional enrichment analysis for multiple gene lists in yeast

Comparative analysis among multiple gene lists on their functional features is now a routine task due to the advancement of high-throughput experiments. Several enrichment analysis tools were developed in the past. However, these tools mainly focus on one gene list and contain only gene ontology or interaction features. What makes it worse, comparative investigation and customized feature set reanalysis are still unavailable. Therefore, we constructed the YMLA (Yeast Multiple List Analyzer) platform in this research. YMLA includes 39 yeast features and facilitates comparative analysis among multiple gene lists via tabular views, heatmaps, and network plots. Moreover, the customized feature set reanalysis function was implemented in YMLA to help form mechanism hypotheses based on a selected enriched feature subset. We demonstrated the biological applicability of YMLA via example lists consisting of genes with top/bottom translation efficiency values. The analysis results provided by YMLA reveal novel facts consistent with previous experiments. YMLA is available at https://cosbi7.ee.ncku.edu.tw/YMLA/.

Comprehensive Analysis of Codon Usage in Quercus Chloroplast Genome and Focus on psbA Gene

Quercus (oak) is an important economic and ecological tree species in the world, and it is the necessary feed for oak silkworm feeding. Chloroplasts play an important role in green plants but the codon usage of oak chloroplast genomes is not fully studied. We examined the codon usage of the oak chloroplast genomes in detail to facilitate the understanding of their biology and evolution. We downloaded all the protein coding genes of 26 non-redundant chloroplast reference genomes, removed short ones and those containing internal stop codons, and finally retained 50 genes shared by all genomes for comparative analyses. The base composition, codon bias, and codon preference are not significantly different between genomes but are significantly different among genes within these genomes. Oak chloroplast genomes prefer T/A-ending codons and avoid C/G-ending codons, and the psbA gene has the same preference except for the codons encoding amino acid Phe. Complex factors such as context-dependent mutations are the major factors affecting codon usage in these genomes, while selection plays an important role on the psbA gene. Our study provided an important understanding of codon usage in the oak chloroplast genomes and found that the psbA gene has nearly the same codon usage preference as other genes in the oak chloroplasts.

A nascent peptide code for translational control of mRNA stability in human cells

Stability of eukaryotic mRNAs is associated with their codon, amino acid, and GC content. Yet, coding sequence motifs that predictably alter mRNA stability in human cells remain poorly defined. Here, we develop a massively parallel assay to measure mRNA effects of thousands of synthetic and endogenous coding sequence motifs in human cells. We identify several families of simple dipeptide repeats whose translation triggers mRNA destabilization. Rather than individual amino acids, specific combinations of bulky and positively charged amino acids are critical for the destabilizing effects of dipeptide repeats. Remarkably, dipeptide sequences that form extended β strands in silico and in vitro slowdown ribosomes and reduce mRNA levels in vivo. The resulting nascent peptide code underlies the mRNA effects of hundreds of endogenous peptide sequences in the human proteome. Our work suggests an intrinsic role for the ribosome as a selectivity filter against the synthesis of bulky and aggregation-prone peptides.

Functional Characterization of Laccase Isozyme (PoLcc1) from the Edible Mushroom Pleurotus ostreatus Involved in Lignin Degradation in Cotton Straw

Fungal laccases play important roles in the degradation of lignocellulose. In this study, the laccase producing cotton straw medium for Pleurotus ostreatus was optimized by single-factor and orthogonal experiments, and to investigate the role of Lacc1 gene, one of the laccase-encoding genes, in the degradation of cotton straw lignin, an overexpression strain of Lacc1 gene was constructed, which was analyzed for the characteristics of lignin degradation. The results demonstrated that the culture conditions with the highest lignin degradation efficiency of the P. ostreatus were the cotton straw particle size of 0.75 mm, a solid-liquid ratio of 1:3 and containing 0.25 g/L of Tween in the medium, as well as an incubation temperature of 26 °C. Two overexpression strains (OE L1-1 and OE L1-4) of Lacc1 gene were obtained, and the gene expression increased 12.08- and 33.04-fold, respectively. The results of 1H-NMR and FTIR analyses of significant changes in lignin structure revealed that Lacc1 gene accelerated the degradation of lignin G-units and involved in the cleavage of β-O-4 linkages and the demethylation of lignin units. These findings will help to improve the efficiency of biodelignification and expand our understanding of its mechanism.

Unlocking the promise of mRNA therapeutics

The extraordinary success of mRNA vaccines against coronavirus disease 2019 (COVID-19) has renewed interest in mRNA as a means of delivering therapeutic proteins. Early clinical trials of mRNA therapeutics include studies of paracrine vascular endothelial growth factor (VEGF) mRNA for heart failure and of CRISPR-Cas9 mRNA for a congenital liver-specific storage disease. However, a series of challenges remains to be addressed before mRNA can be established as a general therapeutic modality with broad relevance to both rare and common diseases. An array of new technologies is being developed to surmount these challenges, including approaches to optimize mRNA cargos, lipid carriers with inherent tissue tropism and in vivo percutaneous delivery systems. The judicious integration of these advances may unlock the promise of biologically targeted mRNA therapeutics, beyond vaccines and other immunostimulatory agents, for the treatment of diverse clinical indications.

Evolution and codon usage bias of mitochondrial and nuclear genomes in Aspergillus section Flavi

The fungal genus Aspergillus contains a diversity of species divided into taxonomic sections of closely related species. Section Flavi contains 33 species, many of industrial, agricultural, or medical relevance. Here, we analyze the mitochondrial genomes (mitogenomes) of 20 Flavi species-including 18 newly assembled mitogenomes-and compare their evolutionary history and codon usage bias patterns to their nuclear counterparts. Codon usage bias refers to variable frequencies of synonymous codons in coding DNA and is shaped by a balance of neutral processes and natural selection. All mitogenomes were circular DNA molecules with highly conserved gene content and order. As expected, genomic content, including GC content, and genome size differed greatly between mitochondrial and nuclear genomes. Phylogenetic analysis based on 14 concatenated mitochondrial genes predicted evolutionary relationships largely consistent with those predicted by a phylogeny constructed from 2,422 nuclear genes. Comparing similarities in interspecies patterns of codon usage bias between mitochondrial and nuclear genomes showed that species grouped differently by patterns of codon usage bias depending on whether analyses were performed using mitochondrial or nuclear relative synonymous usage values. We found that patterns of codon usage bias at gene level are more similar between mitogenomes of different species than the mitogenome and nuclear genome of the same species. Finally, we inferred that, although most genes-both nuclear and mitochondrial-deviated from the neutral expectation for codon usage, mitogenomes were not under translational selection while nuclear genomes were under moderate translational selection. These results contribute to the study of mitochondrial genome evolution in filamentous fungi.

Stability and expression of SARS-CoV-2 spike-protein mutations

Protein fold stability likely plays a role in SARS-CoV-2 S-protein evolution, together with ACE2 binding and antibody evasion. While few thermodynamic stability data are available for S-protein mutants, many systematic experimental data exist for their expression. In this paper, we explore whether such expression levels relate to the thermodynamic stability of the mutants. We studied mutation-induced SARS-CoV-2 S-protein fold stability, as computed by three very distinct methods and eight different protein structures to account for method- and structure-dependencies. For all methods and structures used (24 comparisons), computed stability changes correlate significantly (99% confidence level) with experimental yeast expression from the literature, such that higher expression is associated with relatively higher fold stability. Also significant, albeit weaker, correlations were seen between stability and ACE2 binding effects. The effect of thermodynamic fold stability may be direct or a correlate of amino acid or site properties, notably the solvent exposure of the site. Correlation between computed stability and experimental expression and ACE2 binding suggests that functional properties of the SARS-CoV-2 S-protein mutant space are largely determined by a few simple features, due to underlying correlations. Our study lends promise to the development of computational tools that may ideally aid in understanding and predicting SARS-CoV-2 S-protein evolution.

The Shift in Synonymous Codon Usage Reveals Similar Genomic Variation during Domestication of Asian and African Rice

The domestication of wild rice occurred together with genomic variation, including the synonymous nucleotide substitutions that result in synonymous codon usage bias (SCUB). SCUB mirrors the evolutionary specialization of plants, but its characteristics during domestication were not yet addressed. Here, we found cytosine- and guanidine-ending (NNC and NNG) synonymous codons (SCs) were more pronounced than adenosine- and thymine-ending SCs (NNA and NNT) in both wild and cultivated species of Asian and African rice. The ratios of NNC/G to NNA/T codons gradually decreased following the rise in the number of introns, and the preference for NNA/T codons became more obvious in genes with more introns in cultivated rice when compared with those in wild rice. SCUB frequencies were heterogeneous across the exons, with a higher preference for NNA/T in internal exons than in terminal exons. The preference for NNA/T in internal but not terminal exons was more predominant in cultivated rice than in wild rice, with the difference between wild and cultivated rice becoming more remarkable with the rise in exon numbers. The difference in the ratios of codon combinations representing DNA methylation-mediated conversion from cytosine to thymine between wild and cultivated rice coincided with their difference in SCUB frequencies, suggesting that SCUB reveals the possible association between genetic and epigenetic variation during the domestication of rice. Similar patterns of SCUB shift in Asian and African rice indicate that genomic variation occurs in the same non-random manner. SCUB representing non-neutral synonymous mutations can provide insight into the mechanism of genomic variation in domestication and can be used for the genetic dissection of agricultural traits in rice and other crops.

Potassium viroporins as model systems for understanding eukaryotic ion channel behaviour

Ion channels are membrane proteins essential for a plethora of cellular functions including maintaining cell shape, ion homeostasis, cardiac rhythm and action potential in neurons. The complexity and often extensive structure of eukaryotic membrane proteins makes it difficult to understand their basic biological regulation. Therefore, this article suggests, viroporins - the miniature versions of eukaryotic protein homologs from viruses - might serve as model systems to provide insights into behaviour of eukaryotic ion channels in general. The structural requirements for correct assembly of the channel along with the basic functional properties of a K+ channel exist in the minimal design of the viral K+ channels from two viruses, Chlorella virus (Kcv) and Ectocarpus siliculosus virus (Kesv). These small viral proteins readily assemble into tetramers and they sort in cells to distinct target membranes. When these viruses-encoded channels are expressed into the mammalian cells, they utilise their protein machinery and hence can serve as excellent tools to study the cells protein sorting machinery. This combination of small size and robust function makes viral K+ channels a valuable model system for detection of basic structure-function correlations. It is believed that molecular and physiochemical analyses of these viroporins may serve as basis for the development of inhibitors or modulators to ion channel activity for targeting ion channel diseases - so called channelopathies. Therefore, it may provide a potential different scope for molecular pharmacology studies aiming at novel and innovative therapeutics associated with channel related diseases. This article reviews the structural and functional properties of Kcv and Kesv upon expression in mammalian cells and Xenopus oocytes. The mechanisms behind differential protein sorting in Kcv and Kesv are also thoroughly discussed.

Bioinformatic Assessment of Factors Affecting the Correlation between Protein Abundance and Elongation Efficiency in Prokaryotes

Protein abundance is crucial for the majority of genetically regulated cell functions to act properly in prokaryotic organisms. Therefore, developing bioinformatic methods for assessing the efficiency of different stages of gene expression is of great importance for predicting the actual protein abundance. One of these steps is the evaluation of translation elongation efficiency based on mRNA sequence features, such as codon usage bias and mRNA secondary structure properties. In this study, we have evaluated correlation coefficients between experimentally measured protein abundance and predicted elongation efficiency characteristics for 26 prokaryotes, including non-model organisms, belonging to diverse taxonomic groups The algorithm for assessing elongation efficiency takes into account not only codon bias, but also number and energy of secondary structures in mRNA if those demonstrate an impact on predicted elongation efficiency of the ribosomal protein genes. The results show that, for a number of organisms, secondary structures are a better predictor of protein abundance than codon usage bias. The bioinformatic analysis has revealed several factors associated with the value of the correlation coefficient. The first factor is the elongation efficiency optimization type-the organisms whose genomes are optimized for codon usage only have significantly higher correlation coefficients. The second factor is taxonomical identity-bacteria that belong to the class Bacilli tend to have higher correlation coefficients among the analyzed set. The third is growth rate, which is shown to be higher for the organisms with higher correlation coefficients between protein abundance and predicted translation elongation efficiency. The obtained results can be useful for further improvement of methods for protein abundance prediction.

TREAT: Therapeutic RNAs exploration inspired by artificial intelligence technology

Recent advances in RNA engineering have enabled the development of RNA-based therapeutics for a broad spectrum of applications. Developing RNA therapeutics start with targeted RNA screening and move to the drug design and optimization. However, existing target screening tools ignore noncoding RNAs and their disease-relevant regulatory relationships. And designing therapeutic RNAs encounters high computational complexity of multi-objective optimization to overcome the immunogenicity, instability and inefficient translational production. To unlock the therapeutic potential of noncoding RNAs and enable one-stop screening and design of therapeutic RNAs, we have built the platform TREAT. It incorporates 43,087,953 regulatory relationships between coding and noncoding genes from 81 biological networks under different physiological conditions. TREAT introduces graph representation learning with Random Walk Diffusions to perform disease-relevant target screening, in addition to the commonly utilized Topological Degree and PageRank algorithms. Design and optimization of large RNAs or interfering RNAs are both available. To reduce the computational complexity of multi-objective optimization for large RNA, we stratified the features into local and global features. The local features are evaluated on the fixed-length or dynamic-length local bins, whereas the latter are inspired by AI language models of protein sequence. Then the global assessment is performed on refined candidates, thus reducing the enormous search space. Overall, TREAT is a one-stop platform for the screening and designing of therapeutic RNAs, with particular attention to noncoding RNAs and cutting-edge AI technology embedded, leading the progress of innovative therapeutics for challenging diseases. TREAT is freely accessible at https://rna.org.cn/treat.

The High Mutational Sensitivity of ccdA Antitoxin Is Linked to Codon Optimality

Deep mutational scanning studies suggest that synonymous mutations are typically silent and that most exposed, nonactive-site residues are tolerant to mutations. Here, we show that the ccdA antitoxin component of the Escherichia coli ccdAB toxin-antitoxin system is unusually sensitive to mutations when studied in the operonic context. A large fraction (∼80%) of single-codon mutations, including many synonymous mutations in the ccdA gene shows inactive phenotype, but they retain native-like binding affinity towards cognate toxin, CcdB. Therefore, the observed phenotypic effects are largely not due to alterations in protein structure/stability, consistent with a large region of CcdA being intrinsically disordered. E. coli codon preference and strength of ribosome-binding associated with translation of downstream ccdB gene are found to be major contributors of the observed ccdA mutant phenotypes. In select cases, proteomics studies reveal altered ratios of CcdA:CcdB protein levels in vivo, suggesting that the ccdA mutations likely alter relative translation efficiencies of the two genes in the operon. We extend these results by studying single-site synonymous mutations that lead to loss of function phenotypes in the relBE operon upon introduction of rarer codons. Thus, in their operonic context, genes are likely to be more sensitive to both synonymous and nonsynonymous point mutations than inferred previously.

Relative specificity as an important consideration in the big data era

Technological breakthroughs such as high-throughput methods, genomics, single-cell studies, and machine learning have fundamentally transformed research and ushered in the big data era of biology. Nevertheless, current data collections, analyses, and modeling frequently overlook relative specificity, a crucial property of molecular interactions in biochemical systems. Relative specificity describes how, for example, an enzyme reacts with its many substrates at different rates, and how this discriminatory action alone is sufficient to modulate the substrates and downstream events. As a corollary, it is not only important to comprehensively identify an enzyme’s substrates, but also critical to quantitatively determine how the enzyme interacts with the substrates and to evaluate how it shapes subsequent biological outcomes. Genomics and high-throughput techniques have greatly facilitated the studies of relative specificity in the 21st century, and its functional significance has been demonstrated in complex biochemical systems including transcription, translation, protein kinases, RNA-binding proteins, and animal microRNAs (miRNAs), although it remains ignored in most work. Here we analyze recent findings in big data and relative specificity studies and explain how the incorporation of relative specificity concept might enhance our mechanistic understanding of gene functions, biological phenomena, and human diseases.

CPEB1 regulates the inflammatory immune response, phagocytosis, and alternative polyadenylation in microglia

Microglia are myeloid cells of the central nervous system that perform tasks essential for brain development, neural circuit homeostasis, and neural disease. Microglia react to inflammatory stimuli by upregulating inflammatory signaling through several different immune cell receptors such as the Toll-like receptor 4 (TLR4), which signals to several downstream effectors including transforming growth factor beta-activated kinase 1 (TAK1). Here, we show that TAK1 levels are regulated by CPEB1, a sequence-specific RNA binding protein that controls translation as well as RNA splicing and alternative poly(A) site selection in microglia. Lipopolysaccharide (LPS) binds the TLR4 receptor, which in CPEB1-deficient mice leads to elevated expression of ionized calcium binding adaptor molecule 1 (Iba1), a microglial protein that increases with inflammation, and increased levels of the cytokine IL6. This LPS-induced IL6 response is blocked by inhibitors of JNK, p38, ERK, NFκB, and TAK1. In contrast, phagocytosis, which is elevated in CPEB1-deficient microglia, is unaffected by LPS treatment or ERK inhibition, but is blocked by TAK1 inhibition. These data indicate that CPEB1 regulates microglial inflammatory responses and phagocytosis. RNA-seq indicates that these changes in inflammation and phagocytosis are accompanied by changes in RNA levels, splicing, and alternative poly(A) site selection. Thus, CPEB1 regulation of RNA expression plays a role in microglial function.

A critical overview of current progress for COVID-19: development of vaccines, antiviral drugs, and therapeutic antibodies

The novel coronavirus disease (COVID-19) pandemic remains a global public health crisis, presenting a broad range of challenges. To help address some of the main problems, the scientific community has designed vaccines, diagnostic tools and therapeutics for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection. The rapid pace of technology development, especially with regard to vaccines, represents a stunning and historic scientific achievement. Nevertheless, many challenges remain to be overcome, such as improving vaccine and drug treatment efficacies for emergent mutant strains of SARS-CoV-2. Outbreaks of more infectious variants continue to diminish the utility of available vaccines and drugs. Thus, the effectiveness of vaccines and drugs against the most current variants is a primary consideration in the continual analyses of clinical data that supports updated regulatory decisions. The first two vaccines granted Emergency Use Authorizations (EUAs), BNT162b2 and mRNA-1273, still show more than 60% protection efficacy against the most widespread current SARS-CoV-2 variant, Omicron. This variant carries more than 30 mutations in the spike protein, which has largely abrogated the neutralizing effects of therapeutic antibodies. Fortunately, some neutralizing antibodies and antiviral COVID-19 drugs treatments have shown continued clinical benefits. In this review, we provide a framework for understanding the ongoing development efforts for different types of vaccines and therapeutics, including small molecule and antibody drugs. The ripple effects of newly emergent variants, including updates to vaccines and drug repurposing efforts, are summarized. In addition, we summarize the clinical trials supporting the development and distribution of vaccines, small molecule drugs, and therapeutic antibodies with broad-spectrum activity against SARS-CoV-2 strains.

Control of immediate early gene expression by CPEB4-repressor complex-mediated mRNA degradation

BACKGROUND

Cytoplasmic polyadenylation element-binding protein 4 (CPEB4) is known to associate with cytoplasmic polyadenylation elements (CPEs) located in the 3' untranslated region (UTR) of specific mRNAs and assemble an activator complex promoting the translation of target mRNAs through cytoplasmic polyadenylation.

RESULTS

Here, we find that CPEB4 is part of an alternative repressor complex that mediates mRNA degradation by associating with the evolutionarily conserved CCR4-NOT deadenylase complex. We identify human CPEB4 as an RNA-binding protein (RBP) with enhanced association to poly(A) RNA upon inhibition of class I histone deacetylases (HDACs), a condition known to cause widespread degradation of poly(A)-containing mRNA. Photoactivatable ribonucleoside-enhanced crosslinking and immunoprecipitation (PAR-CLIP) analysis using endogenously tagged CPEB4 in HeLa cells reveals that CPEB4 preferentially binds to the 3'UTR of immediate early gene mRNAs, at G-containing variants of the canonical U- and A-rich CPE located in close proximity to poly(A) sites. By transcriptome-wide mRNA decay measurements, we find that the strength of CPEB4 binding correlates with short mRNA half-lives and that loss of CPEB4 expression leads to the stabilization of immediate early gene mRNAs. Akin to CPEB4, we demonstrate that CPEB1 and CPEB2 also confer mRNA instability by recruitment of the CCR4-NOT complex.

CONCLUSIONS

While CPEB4 was previously known for its ability to stimulate cytoplasmic polyadenylation, our findings establish an additional function for CPEB4 as the RNA adaptor of a repressor complex that enhances the degradation of short-lived immediate early gene mRNAs.

Changes in the genomes and methylomes of three Salmonella enterica serovars after long-term storage in ground black pepper

Low moisture foods (LMFs) have traditionally been recognized as safe for consumption, as most bacteria require higher water content to grow. However, outbreaks due to LMF foods are increasing, and the microbial pathogen Salmonella enterica is frequently implicated. S. enterica can survive in LMFs for years, but few serovars have been studied, and the mechanisms which underlie this longevity are not well understood. Here, we determine that S. enterica serovars S. Tennessee, S. Anatum, and S. Reading but not S. Oranienburg can survive in the ground black pepper for 6 years. S. Reading was not previously associated with any LMF. Using both Illumina and Pacific Biosciences sequencing technologies, we also document changes in the genomes and methylomes of the surviving serovars over this 6-year period. The three serovars acquired a small number of single nucleotide polymorphisms (SNPs) including seven substitutions (four synonymous, two non-synonymous, and one substitution in a non-coding region), and two insertion-deletions. Nine distinct N6-methyladenine (m6A) methylated motifs across the three serovars were identified including five which were previously known, G^m6ATC, CAG^m6AG, BATGC^m6AT, CRT^m6AYN⁶CTC, and CC^m6AN⁷TGAG, and four novel serovar-specific motifs, GRT^m6AN⁸TTYG, GA^m6ACN⁷GTA, GAA m^6ACY, and CAA^m6ANCC. Interestingly, the BATGCAT motif was incompletely methylated (35–64% sites across the genome methylated), suggesting a possible role in gene regulation. Furthermore, the number of methylated BATGC^m6AT motifs increased after storage in ground black pepper for 6 years from 475 to 657 (S. Tennessee), 366 to 608 (S. Anatum), and 525 to 570 (S. Reading), thus warranting further study as an adaptive mechanism. This is the first long-term assessment of genomic changes in S. enterica in a low moisture environment, and the first study to examine the methylome of any bacteria over a period of years, to our knowledge. These data contribute to our understanding of S. enterica survival in LMFs, and coupled with further studies, will provide the information necessary to design effective interventions which reduce S. enterica in LMFs and maintain a healthy, safe food supply.

Plastomes of Bletilla (Orchidaceae) and Phylogenetic Implications

The genus Bletilla is a small genus of only five species distributed across Asia, including B. chartacea, B. foliosa, B. formosana, B. ochracea and B. striata, which is of great medicinal importance. Furthermore, this genus is a member of the key tribe Arethuseae (Orchidaceae), harboring an extremely complicated taxonomic history. Recently, the monophyletic status of Bletilla has been challenged, and the phylogenetic relationships within this genus are still unclear. The plastome, which is rich in both sequence and structural variation, has emerged as a powerful tool for understanding plant evolution. Along with four new plastomes, this work is committed to exploring plastomic markers to elucidate the phylogeny of Bletilla. Our results reveal considerable plastomic differences between B. sinensis and the other three taxa in many aspects. Most importantly, the specific features of the IR junction patterns, novel pttRNA structures and codon aversion motifs can serve as useful molecular markers for Bletilla phylogeny. Moreover, based on maximum likelihood and Bayesian inference methods, our phylogenetic analyses based on two datasets of Arethuseae strongly imply that Bletilla is non-monophyletic. Accordingly, our findings from this study provide novel potential markers for species identification, and shed light on the evolution of Bletilla and Arethuseae.

Plasmid-Based Donor Templates for Nonviral CRISPR/Cas9-Mediated Gene Knock-In in Human T Cells

Effective and precise gene editing of T lymphocytes is critical for advancing the understanding of T cell biology and the development of next-generation cellular therapies. Although methods for effective CRISPR/Cas9-mediated gene knock-out in primary human T cells have been developed, complementary techniques for nonviral gene knock-in can be cumbersome and inefficient. Here, we report a simple and efficient method for nonviral CRISPR/Cas9-based gene knock-in utilizing plasmid-based donor DNA templates. © 2022 Wiley Periodicals LLC. Basic Protocol 1: Purification of human CD4⁺ or CD8⁺ T cells from blood Basic Protocol 2: Activation of purified CD4⁺ or CD8⁺ T cells using TransAct CD3/CD28 agonist-conjugated nanomatrix Basic Protocol 3: Preparation of Cas9/sgRNA RNPs Basic Protocol 4: Transfection of CAS9-RNP and knock-in template into human T cells Support Protocol 1: Purity check following magnetic T cell isolation Support Protocol 2: Dextramer staining of TCR-edited T cells Support Protocol 3: Functional characterization of TCR knock-in T cells Support Protocol 4: Detection of knock-in reporter activity in CRISPR/CAS9-edited T cells.

Genetic features of bovine viral diarrhea virus subgenotype 1c in newborn calves at nucleotide and synonymous codon usages

Bovine viral diarrhea virus (BVDV), serving as an important pathogen for newborn calves, poses threat to reproductive and economic losses in the cattle industry. To survey the infection rate and genetic diversity of BVDV in newborn calves in northern China, a total of 676 sera samples of newborn calves were collected from four provinces between 2021 and 2022. All sera samples were individually detected for BVDV infection by RT-PCR and ELISA. Our results showed that the overall serological rate was 9.76% (66/676) and the average positive rate of BVDV RNA was 8.14% (55/676) in the newborn calves. Eight BVDV strains were successfully isolated from RT-PCR positive sera samples, and four isolates displayed the cytopathic effect (CPE). Based on phylogenetic tree at the genome level, the eight strains were classified into subgenotype 1c. Moreover, the BVDV isolates had a close genetic relationship with the GSTZ strain at either nucleotide or codon usage level. Interestingly, in comparison of synonymous codon usage patterns between the BVDV isolates with CPE and ones without CPE, there were four synonymous codons (UCG, CCC, GCA, and AAC) which displayed the significant differences (p < 0.05) at codon usage pattern, suggesting that synonymous codon usage bias might play a role in BVDV-1c biotypes. In addition, the usage of synonymous codons containing CpG dinucleotides was suppressed by the BVDV-1c isolates, reflecting one of strategies of immune evasion of BVDV to its host. Taken together, our study provided data for monitoring and vaccination strategies of BVDV for newborn calves in northern China.

Oligodendrocyte differentiation alters tRNA modifications and codon optimality-mediated mRNA decay

Oligodendrocytes are specialized cells that confer neuronal myelination in the central nervous system. Leukodystrophies associated with oligodendrocyte deficits and hypomyelination are known to result when a number of tRNA metabolism genes are mutated. Thus, for unknown reasons, oligodendrocytes may be hypersensitive to perturbations in tRNA biology. In this study, we survey the tRNA transcriptome in the murine oligodendrocyte cell lineage and find that specific tRNAs are hypomodified in oligodendrocytes within or near the anticodon compared to oligodendrocyte progenitor cells (OPCs). This hypomodified state may be the result of differential expression of key modification enzymes during oligodendrocyte differentiation. Moreover, we observe a concomitant relationship between tRNA hypomodification and tRNA decoding potential; observing oligodendrocyte specific alterations in codon optimality-mediated mRNA decay and ribosome transit. Our results reveal that oligodendrocytes naturally maintain a delicate, hypersensitized tRNA/mRNA axis. We suggest this axis is a potential mediator of pathology in leukodystrophies and white matter disease when further insult to tRNA metabolism is introduced.

N4-acetylcytidine regulates the replication and pathogenicity of enterovirus 71

Chemical modifications are important for RNA function and metabolism. N4-acetylcytidine (ac4C) is critical for the translation and stability of mRNA. Although ac4C is found in RNA viruses, the detailed mechanisms through which ac4C affects viral replication are unclear. Here, we reported that the 5' untranslated region of the enterovirus 71 (EV71) genome was ac4C modified by the host acetyltransferase NAT10. Inhibition of NAT10 and mutation of the ac4C sites within the internal ribosomal entry site (IRES) suppressed EV71 replication. ac4C enhanced viral RNA translation via selective recruitment of PCBP2 to the IRES and boosted RNA stability. Additionally, ac4C increased the binding of RNA-dependent RNA polymerase (3D) to viral RNA. Notably, ac4C-deficient mutant EV71 showed reduced pathogenicity in vivo. Our findings highlighted the essential role of ac4C in EV71 infection and provided insights into potential antiviral treatments.

Evolutionary dynamics of codon usages for peste des petits ruminants virus

Peste des petits ruminants virus (PPRV) is an important agent of contagious, acute and febrile viral diseases in small ruminants, while its evolutionary dynamics related to codon usage are still lacking. Herein, we adopted information entropy, the relative synonymous codon usage values and similarity indexes and codon adaptation index to analyze the viral genetic features for 45 available whole genomes of PPRV. Some universal, lineage-specific, and gene-specific genetic features presented by synonymous codon usages of the six genes of PPRV that encode N, P, M, F, H and L proteins reflected evolutionary plasticity and independence. The high adaptation of PPRV to hosts at codon usages reflected high viral gene expression, but some synonymous codons that are rare in the hosts were selected in high frequencies in the viral genes. Another obvious genetic feature was that the synonymous codons containing CpG dinucleotides had weak tendencies to be selected in viral genes. The synonymous codon usage patterns of PPRV isolated during 2007–2008 and 2013–2014 in China displayed independent evolutionary pathway, although the overall codon usage patterns of these PPRV strains matched the universal codon usage patterns of lineage IV. According to the interplay between nucleotide and synonymous codon usages of the six genes of PPRV, the evolutionary dynamics including mutation pressure and natural selection determined the viral survival and fitness to its host.

Cotranslational Mechanisms of Protein Biogenesis and Complex Assembly in Eukaryotes

The formation of protein complexes is crucial to most biological functions. The cellular mechanisms governing protein complex biogenesis are not yet well understood, but some principles of cotranslational and posttranslational assembly are beginning to emerge. In bacteria, this process is favored by operons encoding subunits of protein complexes. Eukaryotic cells do not have polycistronic mRNAs, raising the question of how they orchestrate the encounter of unassembled subunits. Here we review the constraints and mechanisms governing eukaryotic co- and posttranslational protein folding and assembly, including the influence of elongation rate on nascent chain targeting, folding, and chaperone interactions. Recent evidence shows that mRNAs encoding subunits of oligomeric assemblies can undergo localized translation and form cytoplasmic condensates that might facilitate the assembly of protein complexes. Understanding the interplay between localized mRNA translation and cotranslational proteostasis will be critical to defining protein complex assembly in vivo.

Mitotic checkpoint gene expression is tuned by codon usage bias

The mitotic checkpoint (also called spindle assembly checkpoint, SAC) is a signaling pathway that safeguards proper chromosome segregation. Correct functioning of the SAC depends on adequate protein concentrations and appropriate stoichiometries between SAC proteins. Yet very little is known about the regulation of SAC gene expression. Here, we show in the fission yeast Schizosaccharomyces pombe that a combination of short mRNA half-lives and long protein half-lives supports stable SAC protein levels. For the SAC genes mad2+ and mad3+ , their short mRNA half-lives are caused, in part, by a high frequency of nonoptimal codons. In contrast, mad1+ mRNA has a short half-life despite a higher frequency of optimal codons, and despite the lack of known RNA-destabilizing motifs. Hence, different SAC genes employ different strategies of expression. We further show that Mad1 homodimers form co-translationally, which may necessitate a certain codon usage pattern. Taken together, we propose that the codon usage of SAC genes is fine-tuned to ensure proper SAC function. Our work shines light on gene expression features that promote spindle assembly checkpoint function and suggests that synonymous mutations may weaken the checkpoint.

A large-scale analysis of codon usage bias in 4868 bacterial genomes shows association of codon adaptation index with GC content, protein functional domains and bacterial phenotypes

Multiple synonymous codons code for the same amino acid, resulting in the degeneracy of the genetic code and in the preferred used of some codons called codon bias usage (CBU). We performed a large-scale analysis of codon usage bias analysing the distribution of the codon adaptation index (CAI) and the codon relative adaptiveness index (RA) in 4868 bacterial genomes. We found that CAI values differ significantly between protein functional domains and part of the protein outside domains and show how CAI, GC content and preferred usage of polymerase III alpha subunits are related. Additionally, we give evidence of the association between CAI and bacterial phenotypes.

Pan-cancer analyses of synonymous mutations based on tissue-specific codon optimality

•

Developed tissue-specific codon optimality in 29 human tissues.

•

Applied these to analyze synonymous mutations in ∼10,000 tumor and normal samples.

•

Synonymous mutations frequently increase optimal codons in most cancer types.

•

Synonymous mutations frequently increase optimal codons cell cycle-related genes.

•

Frequency of optimal codon gain relates to proliferation, DDR deficiency, and survival.

Codon optimality has been demonstrated to be an important determinant of mRNA stability and expression levels in multiple model organisms and human cell lines. However, tissue-specific codon optimality has not been developed to investigate how codon optimality is usually perturbed by somatic synonymous mutations in human cancers. Here, we determined tissue-specific codon optimality in 29 human tissues based on mRNA expression data from the Genotype-Tissue Expression project. We found that optimal codons were associated with differentiation, whereas non-optimal codons were correlated with proliferation. Furthermore, codons biased toward differentiation displayed greater tissue specificity in codon optimality, and the tissue specificity of codon optimality was primarily present in amino acids with high degeneracy of the genetic code. By applying tissue-specific codon optimality to somatic synonymous mutations in 8532 tumor samples across 24 cancer types and to those in 416 normal cells across six human tissues, we found that synonymous mutations frequently increased optimal codons in tumor cells and cancer-related genes (e.g., genes involved in cell cycle). Furthermore, an elevated frequency of optimal codon gain was found to promote tumor cell proliferation in three cancer types characterized by DNA damage repair deficiency and could act as a prognostic biomarker for patients with triple-negative breast cancer. In summary, this study profiled tissue-specific codon optimality in human tissues, revealed alterations in codon optimality caused by synonymous mutations in human cancers, and highlighted the non-negligible role of optimal codon gain in tumorigenesis and therapeutics.

DNA methylation differences during development distinguish sympatric morphs of Arctic charr (Salvelinus alpinus)

Changes in DNA methylation in specific coding or non-coding regions can influence development and potentially divergence in traits within species and groups. While the impact of epigenetic variation on developmental pathways associated with evolutionary divergence is the focus of intense investigation, few studies have looked at recently diverged systems. Phenotypic diversity between closely related populations of Arctic charr (Salvelinus alpinus), which diverged within the last 10 000 years, offers an interesting ecological model to address such effects. Using bisulfite sequencing, we studied general DNA methylation patterns during development in the four sympatric morphs of Arctic charr from Lake Thingvallavatn. The data revealed strong differences between developmental timepoints and between morphs (mainly along the benthic - limnetic axis), both at single CpG sites and in 1,000bp-regions. Genes located close to differentially methylated CpG sites were involved in nucleosome assembly, regulation of osteoclast differentiation, and cell-matrix adhesion. Differentially methylated regions were enriched in tRNA and rRNA sequences, and half of them were located close to transcription start sites. The expression of 14 genes showing methylation differences over time or between morphs was further investigated by qPCR and nine of these were found to be differentially expressed between morphs. Four genes (ARHGEF37-like, H3-like, MPP3 and MEGF9) showed a correlation between methylation and expression. Lastly, histone gene clusters displayed interesting methylation differences between timepoints and morphs, as well as intragenic methylation variation. The results presented here provide a motivation for further studies on the contribution of epigenetic traits, such as DNA methylation, to phenotypic diversity and developmental mechanisms.

Tracing genetic signatures of bat-to-human coronaviruses and early transmission of North American SARS-CoV-2

Highly pathogenic coronaviruses, including SARS-CoV-2, SARS-CoV and MERS-CoV, are thought to be transmitted from bats to humans, but the viral genetic signatures that contribute to bat-to-human transmission remain largely obscure. In this study, we identified an identical ribosomal frameshift motif among the three bat-human pairs of viruses and strong purifying selection after jumping from bats to humans. This represents genetic signatures of coronaviruses that are related to bat-to-human transmission. To further trace the early human-to-human transmission of SARS-CoV-2 in North America, a geographically stratified genome-wide association study (North American isolates and the remaining isolates) and a retrospective study were conducted. We determined that the single nucleotide polymorphisms (SNPs) 1,059.C > T and 25,563.G > T were significantly associated with approximately half of the North American SARS-CoV-2 isolates that accumulated largely during March 2020. Retrospectively tracing isolates with these two SNPs was used to reconstruct the early, reliable transmission history of North American SARS-CoV-2, and European isolates (February 26, 2020) showed transmission 3 days earlier than North American isolates and 17 days earlier than Asian isolates. Collectively, we identified the genetic signatures of the three pairs of coronaviruses and reconstructed an early transmission history of North American SARS-CoV-2. We envision that these genetic signatures are possibly diagnosable and predic markers for public health surveillance.

Gene Therapy in Orthopaedics: Progress and Challenges in Pre-Clinical Development and Translation

In orthopaedics, gene-based treatment approaches are being investigated for an array of common -yet medically challenging- pathologic conditions of the skeletal connective tissues and structures (bone, cartilage, ligament, tendon, joints, intervertebral discs etc.). As the skeletal system protects the vital organs and provides weight-bearing structural support, the various tissues are principally composed of dense extracellular matrix (ECM), often with minimal cellularity and vasculature. Due to their functional roles, composition, and distribution throughout the body the skeletal tissues are prone to traumatic injury, and/or structural failure from chronic inflammation and matrix degradation. Due to a mixture of environment and endogenous factors repair processes are often slow and fail to restore the native quality of the ECM and its function. In other cases, large-scale lesions from severe trauma or tumor surgery, exceed the body’s healing and regenerative capacity. Although a wide range of exogenous gene products (proteins and RNAs) have the potential to enhance tissue repair/regeneration and inhibit degenerative disease their clinical use is hindered by the absence of practical methods for safe, effective delivery. Cumulatively, a large body of evidence demonstrates the capacity to transfer coding sequences for biologic agents to cells in the skeletal tissues to achieve prolonged delivery at functional levels to augment local repair or inhibit pathologic processes. With an eye toward clinical translation, we discuss the research progress in the primary injury and disease targets in orthopaedic gene therapy. Technical considerations important to the exploration and pre-clinical development are presented, with an emphasis on vector technologies and delivery strategies whose capacity to generate and sustain functional transgene expression in vivo is well-established.

Design of typical genes for heterologous gene expression

Heterologous protein expression is an important method for analysing cellular functions of proteins, in genetic circuit engineering and in overexpressing proteins for biopharmaceutical applications and structural biology research. The degeneracy of the genetic code, which enables a single protein to be encoded by a multitude of synonymous gene sequences, plays an important role in regulating protein expression, but substantial uncertainty exists concerning the details of this phenomenon. Here we analyse the influence of a profiled codon usage adaptation approach on protein expression levels in the eukaryotic model organism Saccharomyces cerevisiae. We selected green fluorescent protein (GFP) and human α-synuclein (αSyn) as representatives for stable and intrinsically disordered proteins and representing a benchmark and a challenging test case. A new approach was implemented to design typical genes resembling the codon usage of any subset of endogenous genes. Using this approach, synthetic genes for GFP and αSyn were generated, heterologously expressed and evaluated in yeast. We demonstrate that GFP is expressed at high levels, and that the toxic αSyn can be adapted to endogenous, low-level expression. The new software is publicly available as a web-application for performing host-specific protein adaptations to a set of the most commonly used model organisms ( https://odysseus.motorprotein.de ).

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.

Intragenomic variation in non-adaptive nucleotide biases causes underestimation of selection on synonymous codon usage

Patterns of non-uniform usage of synonymous codons vary across genes in an organism and between species across all domains of life. This codon usage bias (CUB) is due to a combination of non-adaptive (e.g. mutation biases) and adaptive (e.g. natural selection for translation efficiency/accuracy) evolutionary forces. Most models quantify the effects of mutation bias and selection on CUB assuming uniform mutational and other non-adaptive forces across the genome. However, non-adaptive nucleotide biases can vary within a genome due to processes such as biased gene conversion (BGC), potentially obfuscating signals of selection on codon usage. Moreover, genome-wide estimates of non-adaptive nucleotide biases are lacking for non-model organisms. We combine an unsupervised learning method with a population genetics model of synonymous coding sequence evolution to assess the impact of intragenomic variation in non-adaptive nucleotide bias on quantification of natural selection on synonymous codon usage across 49 Saccharomycotina yeasts. We find that in the absence of a priori information, unsupervised learning can be used to identify genes evolving under different non-adaptive nucleotide biases. We find that the impact of intragenomic variation in non-adaptive nucleotide bias varies widely, even among closely-related species. We show that the overall strength and direction of translational selection can be underestimated by failing to account for intragenomic variation in non-adaptive nucleotide biases. Interestingly, genes falling into clusters identified by machine learning are also physically clustered across chromosomes. Our results indicate the need for more nuanced models of sequence evolution that systematically incorporate the effects of variable non-adaptive nucleotide biases on codon frequencies.

High-throughput functional evaluation of human cancer-associated mutations using base editors

Comprehensive phenotypic characterization of the many mutations found in cancer tissues is one of the biggest challenges in cancer genomics. In this study, we evaluated the functional effects of 29,060 cancer-related transition mutations that result in protein variants on the survival and proliferation of non-tumorigenic lung cells using cytosine and adenine base editors and single guide RNA (sgRNA) libraries. By monitoring base editing efficiencies and outcomes using surrogate target sequences paired with sgRNA-encoding sequences on the lentiviral delivery construct, we identified sgRNAs that induced a single primary protein variant per sgRNA, enabling linking those mutations to the cellular phenotypes caused by base editing. The functions of the vast majority of the protein variants (28,458 variants, 98%) were classified as neutral or likely neutral; only 18 (0.06%) and 157 (0.5%) variants caused outgrowing and likely outgrowing phenotypes, respectively. We expect that our approach can be extended to more variants of unknown significance and other tumor types.

Global SLAM-seq for accurate mRNA decay determination and identification of NMD targets

Gene expression analysis requires accurate measurements of global RNA degradation rates, earlier problematic with methods disruptive to cell physiology. Recently, metabolic RNA labeling emerged as an efficient and minimally invasive technique applied in mammalian cells. Here, we have adapted SH-linked alkylation for the metabolic sequencing of RNA (SLAM-seq) for a global mRNA stability study in yeast using 4-thiouracil pulse-chase labeling. We assign high-confidence half-life estimates for 67.5% of expressed ORFs, and measure a median half-life of 9.4 min. For mRNAs where half-life estimates exist in the literature, their ranking order was in good agreement with previous data, indicating that SLAM-seq efficiently classifies stable and unstable transcripts. We then leveraged our yeast protocol to identify targets of the nonsense-mediated decay (NMD) pathway by measuring the change in RNA half-lives, instead of steady-state RNA level changes. With SLAM-seq, we assign 580 transcripts as putative NMD targets, based on their measured half-lives in wild-type and upf3Δ mutants. We find 225 novel targets, and observe a strong agreement with previous reports of NMD targets, 61.2% of our candidates being identified in previous studies. This indicates that SLAM-seq is a simpler and more economic method for global quantification of mRNA half-lives. Our adaptation for yeast yielded global quantitative measures of the NMD effect on transcript half-lives, high correlation with RNA half-lives measured previously with more technically challenging protocols, and identification of novel NMD regulated transcripts that escaped prior detection.

Consistent Clustering Pattern of Prokaryotic Genes Based on Base Frequency at the Second Codon Position and its Association with Functional Category Preference

In 2002, our research group observed a gene clustering pattern based on the base frequency of A versus T at the second codon position in the genome of Vibrio cholera and found that the functional category distribution of genes in the two clusters was different. With the availability of a large number of sequenced genomes, we performed a systematic investigation of A2–T2 distribution and found that 2694 out of 2764 prokaryotic genomes have an optimal clustering number of two, indicating a consistent pattern. Analysis of the functional categories of the coding genes in each cluster in 1483 prokaryotic genomes indicated, that 99.33% of the genomes exhibited a significant difference (p < 0.01) in function distribution between the two clusters. Specifically, functional category P was overrepresented in the small cluster of 98.65% of genomes, whereas categories J, K, and L were overrepresented in the larger cluster of over 98.52% of genomes. Lineage analysis uncovered that these preferences appear consistently across all phyla. Overall, our work revealed an almost universal clustering pattern based on the relative frequency of A2 versus T2 and its role in functional category preference. These findings will promote the understanding of the rationality of theoretical prediction of functional classes of genes from their nucleotide sequences and how protein function is determined by DNA sequence.

Graphical abstract

T-G-A Deficiency Pattern in Protein-Coding Genes and Its Potential Reason

If a stop codon appears within one gene, then its translation will be terminated earlier than expected. False folding of premature protein will be adverse to the host; hence, all functional genes would tend to avoid the intragenic stop codons. Therefore, we hypothesize that there will be less frequency of nucleotides corresponding to stop codons at each codon position of genes. Here, we validate this inference by investigating the nucleotide frequency at a large scale and results from 19,911 prokaryote genomes revealed that nucleotides coinciding with stop codons indeed have the lowest frequency in most genomes. Interestingly, genes with three types of stop codons all tend to follow a T-G-A deficiency pattern, suggesting that the property of avoiding intragenic termination pressure is the same and the major stop codon TGA plays a dominant role in this effect. Finally, a positive correlation between the TGA deficiency extent and the base length was observed in start-experimentally verified genes of Escherichia coli (E. coli). This strengthens the proof of our hypothesis. The T-G-A deficiency pattern observed would help to understand the evolution of codon usage tactics in extant organisms.

The Chloroplast Genome of Wild Saposhnikovia divaricata: Genomic Features, Comparative Analysis, and Phylogenetic Relationships

Saposhnikovia divaricata, a well-known Chinese medicinal herb, is the sole species under the genus Saposhnikovia of the Apiaceae subfamily Apioideae Drude. However, information regarding its genetic diversity and evolution is still limited. In this study, the first complete chloroplast genome (cpDNA) of wild S. divaricata was generated using de novo sequencing technology. Similar to the characteristics of Ledebouriella seseloides, the 147,834 bp-long S. divaricata cpDNA contained a large single copy, a small single copy, and two inverted repeat regions. A total of 85 protein-coding, 8 ribosomal RNA, and 36 transfer RNA genes were identified. Compared with five other species, the non-coding regions in the S. divaricata cpDNA exhibited greater variation than the coding regions. Several repeat sequences were also discovered, namely, 33 forward, 14 reverse, 3 complement, and 49 microsatellite repeats. Furthermore, phylogenetic analysis using 47 cpDNA sequences of Apioideae members revealed that L. seseloides and S. divaricata clustered together with a 100% bootstrap value, thereby supporting the validity of renaming L. seseloides to S. divaricata at the genomic level. Notably, S. divaricata was most closely related to Libanotis buchtormensis, which contradicts previous reports. Therefore, these findings provide a valuable foundation for future studies on the genetic diversity and evolution of S. divaricata.

Senecavirus A Enhances Its Adaptive Evolution via Synonymous Codon Bias Evolution

Synonymous codon bias in the viral genome affects protein translation and gene expression, suggesting that the synonymous codon mutant plays an essential role in influencing virulence and evolution. However, how the recessive mutant form contributes to virus evolvability remains elusive. In this paper, we characterize how the Senecavirus A (SVA), a picornavirus, utilizes synonymous codon mutations to influence its evolution, resulting in the adaptive evolution of the virus to adverse environments. The phylogenetic tree and Median-joining (MJ)-Network of these SVA lineages worldwide were constructed to reveal SVA three-stage genetic development clusters. Furthermore, we analyzed the codon bias of the SVA genome of selected strains and found that SVA could increase the GC content of the third base of some amino acid synonymous codons to enhance the viral RNA adaptive evolution. Our results highlight the impact of recessive mutation of virus codon bias on the evolution of the SVA and uncover a previously underappreciated evolutionary strategy for SVA. They also underline the importance of understanding the genetic evolution of SVA and how SVA adapts to the adverse effects of external stress.

Translation Reprogramming as a Novel Therapeutic Target in MAFLD

Approved pharmacotherapies for metabolic-dysfunction-associated fatty liver disease (MAFLD) are lacking. Novel approaches and therapeutic targets that are likely to translate to clinical benefit are required. Targeting components of the translation machinery hold promise as a novel therapeutic approach that can overcome the well-known disease heterogeneity, as dysregulation of mRNA translation is a common feature independent of the MAFLD drivers. In this perspective, recent advances in understanding the role of mRNA translation in MAFLD are discussed, with a particular focus on the potential implications and challenges to "translate" these findings to the clinic, and an overview of similar recent efforts in other diseases is provided.

RBM33 directs the nuclear export of transcripts containing GC-rich elements

Although splicing is a major driver of RNA nuclear export, many intronless RNAs are efficiently exported to the cytoplasm through poorly characterized mechanisms. For example, GC-rich sequences promote nuclear export in a splicing-independent manner, but how GC content is recognized and coupled to nuclear export is unknown. Here, we developed a genome-wide screening strategy to investigate the mechanism of export of NORAD, an intronless cytoplasmic long noncoding RNA (lncRNA). This screen revealed an RNA binding protein, RBM33, that directs the nuclear export of NORAD and numerous other transcripts. RBM33 directly binds substrate transcripts and recruits components of the TREX-NXF1/NXT1 RNA export pathway. Interestingly, high GC content emerged as the feature that specifies RBM33-dependent nuclear export. Accordingly, RBM33 directly binds GC-rich elements in target transcripts. These results provide a broadly applicable strategy for the genetic dissection of nuclear export mechanisms and reveal a long-sought nuclear export pathway for transcripts with GC-rich sequences.

Codon optimality-mediated mRNA degradation: Linking translational elongation to mRNA stability

Messenger RNA (mRNA) translation by the ribosome represents the final step of a complicated molecular dance from DNA to protein. Although classically considered a decipherer that translates a 64-word genetic code into a proteome of astonishing complexity, the ribosome can also shape the transcriptome by controlling mRNA stability. Recent work has discovered that the ribosome is an arbiter of the general mRNA degradation pathway, wherein the ribosome transit rate serves as a major determinant of transcript half-lives. Specifically, members of the degradation complex sense ribosome translocation rates as a function of ribosome elongation rates. Central to this notion is the concept of codon optimality: although all codons impact translation rates, some are deciphered quickly, whereas others cause ribosome hesitation as a consequence of relative cognate tRNA concentration. These transient pauses induce a unique ribosome conformational state that is probed by the deadenylase complex, thereby inducing an orchestrated set of events that enhance both poly(A) shortening and cap removal. Together, these data imply that the coding region of an mRNA not only encodes for protein content but also impacts protein levels through determining the transcript's fate.

Characterization of tRNA expression profiles in large offspring syndrome

Background

Assisted Reproductive Technologies (ART) use can increase the risk of congenital overgrowth syndromes, such as large offspring syndrome (LOS) in ruminants. Epigenetic variations are known to influence gene expression and differentially methylated regions (DMRs) were previously determined to be associated with LOS in cattle. We observed DMRs overlapping tRNA clusters which could affect tRNA abundance and be associated with tissue specificity or overgrowth. Variations in tRNA expression have been identified in several disease pathways suggesting an important role in the regulation of biological processes. Understanding the role of tRNA expression in cattle offers an opportunity to reveal mechanisms of regulation at the translational level. We analyzed tRNA expression in the skeletal muscle and liver tissues of day 105 artificial insemination-conceived, ART-conceived with a normal body weight, and ART-conceived bovine fetuses with a body weight above the 97^th percentile compared to Control-AI.

Results

Despite the centrality of tRNAs to translation, in silico predictions have revealed dramatic differences in the number of tRNA genes between humans and cattle (597 vs 1,659). Consistent with reports in human, only a fraction of predicted tRNA genes are expressed. We detected the expression of 474 and 487 bovine tRNA genes in the muscle and liver with the remainder being unexpressed. 193 and 198 unique tRNA sequences were expressed in all treatment groups within muscle and liver respectively. In addition, an average of 193 tRNA sequences were expressed within the same treatment group in different tissues. Some tRNA isodecoders were differentially expressed between treatment groups. In the skeletal muscle and liver, we categorized 11 tRNA isoacceptors with undetected expression as well as an isodecoder that was unexpressed in the liver (Ser^GGA). Our results identified variation in the proportion of tRNA gene copies expressed between tissues and differences in the highest contributing tRNA anticodon within an amino acid family due to treatment and tissue type. Out of all amino acid families, roughly half of the most highly expressed tRNA isoacceptors correlated to their most frequent codon in the bovine genome.

Conclusion

Although the number of bovine tRNA genes is nearly triple of that of the tRNA genes in human, there is a shared occurrence of transcriptionally inactive tRNA genes in both species. We detected differential expression of tRNA genes as well as tissue- and treatment- specific tRNA transcripts with unique sequence variations that could modulate translation during protein homeostasis or cellular stress, and give rise to regulatory products targeting genes related to overgrowth in the skeletal muscle and/or tumor development in the liver of LOS individuals. While the absence of certain isodecoders may be relieved by wobble base pairing, missing tRNA species could increase the likelihood of mistranslation or mRNA degradation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08496-7.

mRNA vaccines: the most recent clinical applications of synthetic mRNA

Synthetic mRNA has been considered as an emerging biotherapeutic agent for the past decades. Recently, the SARS-CoV-2 pandemic has led to the first clinical use of synthetic mRNA. mRNA vaccines showed far surpassing influences on the public as compared to other vaccine platforms such as viral vector vaccines and recombinant protein vaccines. It allowed rapid development and production of vaccines that have never been achieved in history. Synthetic mRNA, called in vitro transcribed (IVT) mRNA, is the key component of mRNA vaccines. It has several advantages over conventional gene-expressing systems such as plasmid DNA and viral vectors. It can translate proteins in the cytoplasm by structurally resembling natural mRNA and exhibit various protein expression patterns depending on how it is engineered. Another advantage is that synthetic mRNA enables fast, scalable, and cost-effective production. Therefore, starting with the mRNA vaccine, synthetic mRNA is now in the spotlight as a promising new drug development agent. In this review, we will summarize the latest IVT mRNA technology such as new mRNA structures or large-scale production. In addition, the nature of the innate immunogenicity of IVT mRNA will be discussed along with its roles in the development of vaccines. Finally, the principles of the mRNA vaccine and the future direction of synthetic mRNA will be provided.

Synonymous Codon Pattern of Cowpea Mild Mottle Virus Sheds Light on Its Host Adaptation and Genome Evolution

Cowpea mild mottle virus (CpMMV) is an economically significant virus that causes severe disease on several legume crops. Aside from recombination, other factors driving its rapid evolution are elusive. In this study, the synonymous codon pattern of CpMMV and factors shaping it were analyzed. Phylogeny and nucleotide composition analyses showed that isolates of different geography or hosts had very similar nucleotide compositions. Relative synonymous codon usage (RSCU) and neutrality analyses suggest that CpMMV prefers A/U-ending codons and natural selection is the dominative factor that affects its codon bias. Dinucleotide composition and codon adaptation analyses indicate that the codon pattern of CpMMV is mainly shaped by the requirement of escaping of host dinucleotide-associated antiviral responses and translational efficiency.

Eukaryotic mRNA Decapping Activation

The 5'-terminal cap is a fundamental determinant of eukaryotic gene expression which facilitates cap-dependent translation and protects mRNAs from exonucleolytic degradation. Enzyme-directed hydrolysis of the cap (decapping) decisively affects mRNA expression and turnover, and is a heavily regulated event. Following the identification of the decapping holoenzyme (Dcp1/2) over two decades ago, numerous studies revealed the complexity of decapping regulation across species and cell types. A conserved set of Dcp1/2-associated proteins, implicated in decapping activation and molecular scaffolding, were identified through genetic and molecular interaction studies, and yet their exact mechanisms of action are only emerging. In this review, we discuss the prevailing models on the roles and assembly of decapping co-factors, with considerations of conservation across species and comparison across physiological contexts. We next discuss the functional convergences of decapping machineries with other RNA-protein complexes in cytoplasmic P bodies and compare current views on their impact on mRNA stability and translation. Lastly, we review the current models of decapping activation and highlight important gaps in our current understanding.

Advances in COVID-19 mRNA vaccine development

To date, the coronavirus disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has determined 399,600,607 cases and 5,757,562 deaths worldwide. COVID-19 is a serious threat to human health globally. The World Health Organization (WHO) has declared COVID-19 pandemic a major public health emergency. Vaccination is the most effective and economical intervention for controlling the spread of epidemics, and consequently saving lives and protecting the health of the population. Various techniques have been employed in the development of COVID-19 vaccines. Among these, the COVID-19 messenger RNA (mRNA) vaccine has been drawing increasing attention owing to its great application prospects and advantages, which include short development cycle, easy industrialization, simple production process, flexibility to respond to new variants, and the capacity to induce better immune response. This review summarizes current knowledge on the structural characteristics, antigen design strategies, delivery systems, industrialization potential, quality control, latest clinical trials and real-world data of COVID-19 mRNA vaccines as well as mRNA technology. Current challenges and future directions in the development of preventive mRNA vaccines for major infectious diseases are also discussed.

Combinatorial optimization of mRNA structure, stability, and translation for RNA-based therapeutics

Therapeutic mRNAs and vaccines are being developed for a broad range of human diseases, including COVID-19. However, their optimization is hindered by mRNA instability and inefficient protein expression. Here, we describe design principles that overcome these barriers. We develop an RNA sequencing-based platform called PERSIST-seq to systematically delineate in-cell mRNA stability, ribosome load, as well as in-solution stability of a library of diverse mRNAs. We find that, surprisingly, in-cell stability is a greater driver of protein output than high ribosome load. We further introduce a method called In-line-seq, applied to thousands of diverse RNAs, that reveals sequence and structure-based rules for mitigating hydrolytic degradation. Our findings show that highly structured "superfolder" mRNAs can be designed to improve both stability and expression with further enhancement through pseudouridine nucleoside modification. Together, our study demonstrates simultaneous improvement of mRNA stability and protein expression and provides a computational-experimental platform for the enhancement of mRNA medicines.