Protein functional features are reflected in the patterns of mRNA translation speed

Optimizing interleukin-6 and 8 expression, clarification and purification in plant cell packs and plants for application in advanced therapy medicinal products and cellular agriculture

Healthcare and nutrition are facing a paradigm shift in light of advanced therapy medicinal products (ATMPs) and cellular agriculture options respectively. Both options heavily rely on some sort of animal cell culture, e.g. autologous stem cells. These cultures require various growth factors, such as interleukin-6 and 8 (IL-6/8), in a pure, safe and sustainable form that can be provided in a scalable manner. Plants seem well suited for this task because purification of small proteins can be readily achieved by membrane separation, human/animal pathogens do not replicate in plants and production can be scaled up using in-door farming or agricultural practices. Here, we illustrate this capacity by first optimizing the codon usage of IL-6/8 for translation in Nicotiana spp., as well as testing the effect of untranslated regions and product targeting to different sub-cellular compartments on expression in a high-throughput plant cell pack (PCP) assay. In the chloroplast, IL-6 accumulated up to 6.9±3.8 (SD, n=2) and 14.4±7.4 mg kg-1 (SD, n=5) were observed in case of IL-8. When transferring IL-8 expression into whole plants, accumulation was 12.3±1.5 mg kg-1 (SD, n=3). After extraction and clarification, IL-8 was purified using a two-stage process consisting of an ultrafiltration/diafiltration step with 100 kDa and 10 kDa cut off membranes followed by an IMAC polishing step. The purity, yield and recovery were 97.8%, 6.6 mg kg-1 and 38%, respectively. We evaluated the ability of the proposed purification process to remove endotoxins to ensure the compatibility of plant-made growth factors with cell culture.

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.

Multiplexing mechanical and translational cues on genes

The genetic code gives precise instructions on how to translate codons into amino acids. Due to the degeneracy of the genetic code-18 out of 20 amino acids are encoded for by more than one codon-more information can be stored in a basepair sequence. Indeed, various types of additional information have been discussed in the literature, e.g., the positioning of nucleosomes along eukaryotic genomes and the modulation of the translating efficiency in ribosomes to influence cotranslational protein folding. The purpose of this study is to show that it is indeed possible to carry more than one additional layer of information on top of a gene. In particular, we show how much translation efficiency and nucleosome positioning can be adjusted simultaneously without changing the encoded protein. We achieve this by mapping genes on weighted graphs that contain all synonymous genes, and then finding shortest paths through these graphs. This enables us, for example, to readjust the disrupted translational efficiency profile after a gene has been introduced from one organism (e.g., human) into another (e.g., yeast) without greatly changing the nucleosome landscape intrinsically encoded by the DNA molecule.

Identification of key non-coding RNAs and transcription factors regulators and their potential drugs for steroid-induced femoral head necrosis

Steroid-induced necrosis of femoral head (SINFH) is a femoral head necrotic disease caused by prolonged use of hormones. The detailed pathogenesis has not been fully demonstrated. In this study, we employed the bioinformatics approach to probe the roles of SINFH inhibitors. Core dysfunction modules related to SINFH was obtained. Meanwhile, GO and KEGG analysis of genes in dysfunction modules are carried out. Furthermore, the pivot prediction analysis of dysfunction modules related to ncRNA and transcription factor (TF) has been performed. The functions of the enriched modules were focused on multiple perspectives, including circulation, gland development, bone development and reconstruction, calcium production, and fatty acid metabolism regulation. The ncRNAs and TFs analysis showed that miR-322-5p, miR-124-3p, miR-125a-3p, and Ctnnb1 were important members of SINFH dysfunction. Drug targets suggested that Zinc and adenosine monophosphate may have an impact on SINFH dysfunction. SINFH was closely related to bone development and reconstruction.

The Correlation Between DsRed mRNA Levels and Transient DsRed Protein Expression in Plants Depends on Leaf Age and the 5' Untranslated Region

The yield of recombinant proteins in plants determines their economic competitiveness as a production platform compared to microbes and mammalian cells. The promoter, untranslated regions (UTRs) and codon usage can all contribute to the yield, but potential interactions among these components have not been examined in detail. Here the effect of two promoters (35SS and nos) and four 5'UTRs on the spatiotemporal expression of DsRed mRNA and the accumulation of DsRed protein during transient expression in tobacco (Nicotiana tabacum) mediated by Agrobacterium tumefaciens is investigated. The authors found that the mRNA levels peaked 2-3 days post-infiltration (dpi), and rapidly declined thereafter, whereas DsRed protein was first detected after ≈3 days and concentrations continued to increase until at least 5 dpi. This temporal decoupling of mRNA and protein expression was strongest in the older leaves, which also produced the lowest DsRed yields. The accumulation of DsRed linearly correlated with mRNA levels in all but the youngest leaves, where more DsRed was synthesized per mRNA molecule. This was the case for both promoters, although the nos promoter had a higher protein/mRNA ratio than the 35SS promoter. Furthermore, the type of 5'UTR affected DsRed protein accumulation by 50% starting from similar levels of mRNA. The authors concluded that DsRed mRNA levels are not the limiting factor for DsRed protein expression in plants, but that translation-associated processes such as initiation, elongation, and release are bottlenecks that should be addressed in future studies.

Genetic Code Optimization for Cotranslational Protein Folding: Codon Directional Asymmetry Correlates with Antiparallel Betasheets, tRNA Synthetase Classes

A new codon property, codon directional asymmetry in nucleotide content (CDA), reveals a biologically meaningful genetic code dimension: palindromic codons (first and last nucleotides identical, codon structure XZX) are symmetric (CDA = 0), codons with structures ZXX/XXZ are 5'/3' asymmetric (CDA = - 1/1; CDA = - 0.5/0.5 if Z and X are both purines or both pyrimidines, assigning negative/positive (-/+) signs is an arbitrary convention). Negative/positive CDAs associate with (a) Fujimoto's tetrahedral codon stereo-table; (b) tRNA synthetase class I/II (aminoacylate the 2'/3' hydroxyl group of the tRNA's last ribose, respectively); and (c) high/low antiparallel (not parallel) betasheet conformation parameters. Preliminary results suggest CDA-whole organism associations (body temperature, developmental stability, lifespan). Presumably, CDA impacts spatial kinetics of codon-anticodon interactions, affecting cotranslational protein folding. Some synonymous codons have opposite CDA sign (alanine, leucine, serine, and valine), putatively explaining how synonymous mutations sometimes affect protein function. Correlations between CDA and tRNA synthetase classes are weaker than between CDA and antiparallel betasheet conformation parameters. This effect is stronger for mitochondrial genetic codes, and potentially drives mitochondrial codon-amino acid reassignments. CDA reveals information ruling nucleotide-protein relations embedded in reversed (not reverse-complement) sequences (5'-ZXX-3'/5'-XXZ-3').

Comparative survey of the relative impact of mRNA features on local ribosome profiling read density

Ribosome profiling (Ribo-seq), a promising technology for exploring ribosome decoding rates, is characterized by the presence of infrequent high peaks in ribosome footprint density and by long alignment gaps. Here, to reduce the impact of data heterogeneity we introduce a simple normalization method, Ribo-seq Unit Step Transformation (RUST). RUST is robust and outperforms other normalization techniques in the presence of heterogeneous noise. We illustrate how RUST can be used for identifying mRNA sequence features that affect ribosome footprint densities globally. We show that a few parameters extracted with RUST are sufficient for predicting experimental densities with high accuracy. Importantly the application of RUST to 30 publicly available Ribo-seq data sets revealed a substantial variation in sequence determinants of ribosome footprint frequencies, questioning the reliability of Ribo-seq as an accurate representation of local ribosome densities without prior quality control. This emphasizes our incomplete understanding of how protocol parameters affect ribosome footprint densities.

Ribosome profiling data can suffer from uneven coverage which hampers estimation of elongation rates. Connor et al. present an enhanced data smoothing method for Ribo-seq data and highlight significant variability in sequence determinants of ribosome density in publicly available data sets.

Distinct interacting core taxa in co-occurrence networks enable discrimination of polymicrobial oral diseases with similar symptoms

Polymicrobial diseases, which can be life threatening, are caused by the presence and interactions of multiple microbes. Peri-implantitis and periodontitis are representative polymicrobial diseases that show similar clinical symptoms. To establish a means of differentiating between them, we compared microbial species and functional genes in situ by performing metatranscriptomic analyses of peri-implantitis and periodontitis samples obtained from the same subjects (n = 12 each). Although the two diseases differed in terms of 16S rRNA-based taxonomic profiles, they showed similarities with respect to functional genes and taxonomic and virulence factor mRNA profiles. The latter—defined as microbial virulence types—differed from those of healthy periodontal sites. We also showed that networks based on co-occurrence relationships of taxonomic mRNA abundance (co-occurrence networks) were dissimilar between the two diseases. Remarkably, these networks consisted mainly of taxa with a high relative mRNA-to-rRNA ratio, with some showing significant co-occurrence defined as interacting core taxa, highlighting differences between the two groups. Thus, peri-implantitis and periodontitis have shared as well as distinct microbiological characteristics. Our findings provide insight into microbial interactions in polymicrobial diseases with unknown etiologies.

Nucleoside modifications in the regulation of gene expression: focus on tRNA

Both, DNA and RNA nucleoside modifications contribute to the complex multi-level regulation of gene expression. Modified bases in tRNAs modulate protein translation rates in a highly dynamic manner. Synonymous codons, which differ by the third nucleoside in the triplet but code for the same amino acid, may be utilized at different rates according to codon-anticodon affinity. Nucleoside modifications in the tRNA anticodon loop can favor the interaction with selected codons by stabilizing specific base pairs. Similarly, weakening of base pairing can discriminate against binding to near-cognate codons. mRNAs enriched in favored codons are translated in higher rates constituting a fine-tuning mechanism for protein synthesis. This so-called codon bias establishes a basic protein level, but sometimes it is necessary to further adjust the production rate of a particular protein to actual requirements, brought by, e.g., stages in circadian rhythms, cell cycle progression or exposure to stress. Such an adjustment is realized by the dynamic change of tRNA modifications resulting in the preferential translation of mRNAs coding for example for stress proteins to facilitate cell survival. Furthermore, tRNAs contribute in an entirely different way to another, less specific stress response consisting in modification-dependent tRNA cleavage that contributes to the general down-regulation of protein synthesis. In this review, we summarize control functions of nucleoside modifications in gene regulation with a focus on recent findings on protein synthesis control by tRNA base modifications.