Codon usage bias controls mRNA and protein abundance in trypanosomatids

Strategies for improving the genome-editing efficiency of class 2 CRISPR/Cas system

Since its advent, gene-editing technology has been widely used in microorganisms, animals, plants, and other species. This technology shows remarkable application prospects, giving rise to a new biotechnological industry. In particular, third-generation gene editing technology, represented by the CRISPR/Cas9 system, has become the mainstream gene editing technology owing to its advantages of high efficiency, simple operation, and low cost. These systems can be widely used because they have been modified and optimized, leading to notable improvements in the efficiency of gene editing. This review introduces the characteristics of popular CRISPR/Cas systems and optimization methods aimed at improving the editing efficiency of class 2 CRISPR/Cas systems, providing a reference for the development of superior gene editing systems. Additionally, the review discusses the development and optimization of base editors, primer editors, gene activation and repression tools, as well as the advancement and refinement of compact systems such as IscB, TnpB, Fanzor, and Cas12f.

Post-transcriptional reprogramming by thousands of mRNA untranslated regions in trypanosomes

Although genome-wide polycistronic transcription places major emphasis on post-transcriptional controls in trypanosomatids, messenger RNA cis-regulatory untranslated regions (UTRs) have remained largely uncharacterised. Here, we describe a genome-scale massive parallel reporter assay coupled with 3'-UTR-seq profiling in the African trypanosome and identify thousands of regulatory UTRs. Increased translation efficiency was associated with dosage of adenine-rich poly-purine tracts (pPuTs). An independent assessment of native UTRs using machine learning based predictions confirmed the robust correspondence between pPuTs and positive control, as did an assessment of synthetic UTRs. Those 3'-UTRs associated with upregulated expression in bloodstream-stage cells were also enriched in uracil-rich poly-pyrimidine tracts, suggesting a mechanism for developmental activation through pPuT 'unmasking'. Thus, we describe a cis-regulatory UTR sequence 'code' that underpins gene expression control in the context of a constitutively transcribed genome. We conclude that thousands of UTRs post-transcriptionally reprogram gene expression profiles in trypanosomes.

A Modular Cloning Toolkit for the production of recombinant proteins in Leishmania tarentolae

Modular Cloning (MoClo) is based on libraries of standardized genetic parts that can be directionally assembled via Golden Gate cloning in one-pot reactions into transcription units and multigene constructs. Here, a team of bachelor students established a MoClo toolkit for the protist Leishmania tarentolae in the frame of the international Genetically Engineered Machine (iGEM) competition. Our modular toolkit is based on a domesticated version of a commercial LEXSY expression vector and comprises 34 genetic parts encoding various affinity tags, targeting signals as well as fluorescent and luminescent proteins. We demonstrated the utility of our kit by the successful production of 16 different tagged versions of the receptor binding domain (RBD) of the SARS-CoV-2 spike protein in L. tarentolae liquid cultures. While highest yields of secreted recombinant RBD were obtained for GST-tagged fusion proteins 48 h post induction, C-terminal peptide tags were often degraded and resulted in lower yields of secreted RBD. Fusing secreted RBD to a synthetic O-glycosylation SP20 module resulted in an apparent molecular mass shift around 10 kDa. No disadvantage regarding the production of RBD was detected when the three antibiotics of the LEXSY system were omitted during the 48-h induction phase. Furthermore, the successful purification of secreted RBD from the supernatant of L. tarentolae liquid cultures was demonstrated in pilot experiments. In summary, we established a MoClo toolkit and exemplified its application for the production of recombinant proteins in L. tarentolae.

Bioinformatic Analysis Reveals the Role of Translation Elongation Efficiency Optimisation in the Evolution of Ralstonia Genus

Translation efficiency modulates gene expression in prokaryotes. The comparative analysis of translation elongation efficiency characteristics of Ralstonia genus bacteria genomes revealed that these characteristics diverge in accordance with the phylogeny of Ralstonia. The first branch of this genus is a group of bacteria commonly found in moist environments such as soil and water that includes the species R. mannitolilytica, R. insidiosa, and R. pickettii, which are also described as nosocomial infection pathogens. In contrast, the second branch is plant pathogenic bacteria consisting of R. solanacearum, R. pseudosolanacearum, and R. syzygii. We found that the soil Ralstonia have a significantly lower number and energy of potential secondary structures in mRNA and an increased role of codon usage bias in the optimization of highly expressed genes' translation elongation efficiency, not only compared to phytopathogenic Ralstonia but also to Cupriavidus necator, which is closely related to the Ralstonia genus. The observed alterations in translation elongation efficiency of orthologous genes are also reflected in the difference of potentially highly expressed gene' sets' content among Ralstonia branches with different lifestyles. Analysis of translation elongation efficiency characteristics can be considered a promising approach for studying complex mechanisms that determine the evolution and adaptation of bacteria in various environments.

Specific recognition and ubiquitination of translating ribosomes by mammalian CCR4-NOT

Translation affects messenger RNA stability and, in yeast, this is mediated by the Ccr4-Not deadenylation complex. The details of this process in mammals remain unclear. Here, we use cryogenic electron microscopy (cryo-EM) and crosslinking mass spectrometry to show that mammalian CCR4-NOT specifically recognizes ribosomes that are stalled during translation elongation in an in vitro reconstituted system with rabbit and human components. Similar to yeast, mammalian CCR4-NOT inserts a helical bundle of its CNOT3 subunit into the empty E site of the ribosome. Our cryo-EM structure shows that CNOT3 also locks the L1 stalk in an open conformation to inhibit further translation. CCR4-NOT is required for stable association of the nonconstitutive subunit CNOT4, which ubiquitinates the ribosome, likely to signal stalled translation elongation. Overall, our work shows that human CCR4-NOT not only detects but also enforces ribosomal stalling to couple translation and mRNA decay.

Translation and mRNA Stability Control

Messenger RNA (mRNA) stability and translational efficiency are two crucial aspects of the post-transcriptional process that profoundly impact protein production in a cell. While it is widely known that ribosomes produce proteins, studies during the past decade have surprisingly revealed that ribosomes also control mRNA stability in a codon-dependent manner, a process referred to as codon optimality. Therefore, codons, the three-nucleotide words read by the ribosome, have a potent effect on mRNA stability and provide cis-regulatory information that extends beyond the amino acids they encode. While the codon optimality molecular mechanism is still unclear, the translation elongation rate appears to trigger mRNA decay. Thus, transfer RNAs emerge as potential master gene regulators affecting mRNA stability. Furthermore, while few factors related to codon optimality have been identified in yeast, the orthologous genes in vertebrates do not necessary share the same functions. Here, we discuss codon optimality findings and gene regulation layers related to codon composition in different eukaryotic species.

Transcriptomic, proteomic, and functional consequences of codon usage bias in human cells during heterologous gene expression

Differences in codon frequency between genomes, genes, or positions along a gene, modulate transcription and translation efficiency, leading to phenotypic and functional differences. Here, we present a multiscale analysis of the effects of synonymous codon recoding during heterologous gene expression in human cells, quantifying the phenotypic consequences of codon usage bias at different molecular and cellular levels, with an emphasis on translation elongation. Six synonymous versions of an antibiotic resistance gene were generated, fused to a fluorescent reporter, and independently expressed in HEK293 cells. Multiscale phenotype was analyzed by means of quantitative transcriptome and proteome assessment, as proxies for gene expression; cellular fluorescence, as a proxy for single-cell level expression; and real-time cell proliferation in absence or presence of antibiotic, as a proxy for the cell fitness. We show that differences in codon usage bias strongly impact the molecular and cellular phenotype: (i) they result in large differences in mRNA levels and protein levels, leading to differences of over 15 times in translation efficiency; (ii) they introduce unpredicted splicing events; (iii) they lead to reproducible phenotypic heterogeneity; and (iv) they lead to a trade-off between the benefit of antibiotic resistance and the burden of heterologous expression. In human cells in culture, codon usage bias modulates gene expression by modifying mRNA availability and suitability for translation, leading to differences in protein levels and eventually eliciting functional phenotypic changes.

Impact of Growth Rate on the Protein-mRNA Ratio in Pseudomonas aeruginosa

Our understanding of how bacterial pathogens colonize and persist during human infection has been hampered by the limited characterization of bacterial physiology during infection and a research bias toward in vitro, fast-growing bacteria. Recent research has begun to address these gaps in knowledge by directly quantifying bacterial mRNA levels during human infection, with the goal of assessing microbial community function at the infection site. However, mRNA levels are not always predictive of protein levels, which are the primary functional units of a cell. Here, we used carefully controlled chemostat experiments to examine the relationship between mRNA and protein levels across four growth rates in the bacterial pathogen Pseudomonas aeruginosa. We found a genome-wide positive correlation between mRNA and protein abundances across all growth rates, with genes required for P. aeruginosa viability having stronger correlations than nonessential genes. We developed a statistical method to identify genes whose mRNA abundances poorly predict protein abundances and calculated an RNA-to-protein (RTP) conversion factor to improve mRNA predictions of protein levels. The application of the RTP conversion factor to publicly available transcriptome data sets was highly robust, enabling the more accurate prediction of P. aeruginosa protein levels across strains and growth conditions. Finally, the RTP conversion factor was applied to P. aeruginosa human cystic fibrosis (CF) infection transcriptomes to provide greater insights into the functionality of this bacterium in the CF lung. This study addresses a critical problem in infection microbiology by providing a framework for enhancing the functional interpretation of bacterial human infection transcriptome data. IMPORTANCE Our understanding of bacterial physiology during human infection is limited by the difficulty in assessing bacterial function at the infection site. Recent studies have begun to address this question by quantifying bacterial mRNA levels in human-derived samples using transcriptomics. One challenge for these studies is the poor predictivity of mRNA for protein levels for some genes. Here, we addressed this challenge by measuring the transcriptomes and proteomes of P. aeruginosa grown at four growth rates. Our results revealed that the growth rate does not impact the genome-wide correlation of mRNA and protein levels. We used statistical methods to identify the genes for which mRNA and protein were poorly correlated and developed an RNA-to-protein (RTP) conversion factor that improved the predictivity of protein levels across strains and growth conditions. Our results provide new insights into mRNA-protein correlations and tools to enhance our understanding of bacterial physiology from transcriptome data.

Synonymous mutations in the phosphoglycerate kinase 1 gene induce an altered response to protein misfolding in Schizosaccharomyces pombe

Background

Proteostasis refers to the processes that regulate the biogenesis, folding, trafficking, and degradation of proteins. Any alteration in these processes can lead to cell malfunction. Protein synthesis, a key proteostatic process, is highly-regulated at multiple levels to ensure adequate adaptation to environmental and physiological challenges such as different stressors, proteotoxic conditions and aging, among other factors. Because alterations in protein translation can lead to protein misfolding, examining how protein translation is regulated may also help to elucidate in part how proteostasis is controlled. Codon usage bias has been implicated in the fine-tuning of translation rate, as more-frequent codons might be read faster than their less-frequent counterparts. Thus, alterations in codon usage due to synonymous mutations may alter translation kinetics and thereby affect the folding of the nascent polypeptide, without altering its primary structure. To date, it has been difficult to predict the effect of synonymous mutations on protein folding and cellular fitness due to a scarcity of relevant data. Thus, the purpose of this work was to assess the effect of synonymous mutations in discrete regions of the gene that encodes the highly-expressed enzyme 3-phosphoglycerate kinase 1 (pgk1) in the fission yeast Schizosaccharomyces pombe.

Results

By means of systematic replacement of synonymous codons along pgk1, we found slightly-altered protein folding and activity in a region-specific manner. However, alterations in protein aggregation, heat stress as well as changes in proteasome activity occurred independently of the mutated region. Concomitantly, reduced mRNA levels of the chaperones Hsp9 and Hsp16 were observed.

Conclusion

Taken together, these data suggest that codon usage bias of the gene encoding this highly-expressed protein is an important regulator of protein function and proteostasis.

Validation of Promoters and Codon Optimization on CRISPR/Cas9-Engineered Jurkat Cells Stably Expressing αRep4E3 for Interfering with HIV-1 Replication

Persistent and efficient therapeutic protein expression in the specific target cell is a significant concern in gene therapy. The controllable integration site, suitable promoter, and proper codon usage influence the effectiveness of the therapeutic outcome. Previously, we developed a non-immunoglobulin scaffold, alpha repeat protein (αRep4E3), as an HIV-1 RNA packaging interference system in SupT1 cells using the lentiviral gene transfer. Although the success of anti-HIV-1 activity was evidenced, the integration site is uncontrollable and may not be practical for clinical translation. In this study, we use the CRISPR/Cas9 gene editing technology to precisely knock-in αRep4E3 genes into the adeno-associated virus integration site 1 (AAVS1) safe harbor locus of the target cells. We compare the αRep4E3 expression under the regulation of three different promoters, including cytomegalovirus (CMV), human elongation factor-1 alpha (EF1α), and ubiquitin C (UbC) promoters with and without codon optimization in HEK293T cells. The results demonstrated that the EF1α promoter with codon-optimized αRep4E3mCherry showed higher protein expression than other promoters with non-optimized codons. We then performed a proof-of-concept study by knocking in the αRep4E3mCherry gene at the AAVS1 locus of the Jurkat cells. The results showed that the αRep4E3mCherry-expressing Jurkat cells exhibited anti-HIV-1 activities against HIV-1NL4-3 strain as evidenced by decreased capsid (p24) protein levels and viral genome copies as compared to the untransfected Jurkat control cells. Altogether, our study demonstrates that the αRep4E3 could interfere with the viral RNA packaging and suggests that the αRep4E3 scaffold protein could be a promising anti-viral molecule that offers a functional cure for people living with HIV-1.

Heterologous reporter expression in the planarian Schmidtea mediterranea through somatic mRNA transfection

Planarians have long been studied for their regenerative abilities. Moving forward, tools for ectopic expression of non-native proteins will be of substantial value. Using a luminescent reporter to overcome the strong autofluorescence of planarian tissues, we demonstrate heterologous protein expression in planarian cells and live animals. Our approach is based on the introduction of mRNA through several nanotechnological and chemical transfection methods. We improve reporter expression by altering untranslated region (UTR) sequences and codon bias, facilitating the measurement of expression kinetics in both isolated cells and whole planarians using luminescence imaging. We also examine protein expression as a function of variations in the UTRs of delivered mRNA, demonstrating a framework to investigate gene regulation at the post-transcriptional level. Together, these advances expand the toolbox for the mechanistic analysis of planarian biology and establish a foundation for the development and expansion of transgenic techniques in this unique model system.

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.

iCodon customizes gene expression based on the codon composition

Messenger RNA (mRNA) stability substantially impacts steady-state gene expression levels in a cell. mRNA stability is strongly affected by codon composition in a translation-dependent manner across species, through a mechanism termed codon optimality. We have developed iCodon (www.iCodon.org), an algorithm for customizing mRNA expression through the introduction of synonymous codon substitutions into the coding sequence. iCodon is optimized for four vertebrate transcriptomes: mouse, human, frog, and fish. Users can predict the mRNA stability of any coding sequence based on its codon composition and subsequently generate more stable (optimized) or unstable (deoptimized) variants encoding for the same protein. Further, we show that codon optimality predictions correlate with both mRNA stability using a massive reporter library and expression levels using fluorescent reporters and analysis of endogenous gene expression in zebrafish embryos and/or human cells. Therefore, iCodon will benefit basic biological research, as well as a wide range of applications for biotechnology and biomedicine.

Variation in upstream open reading frames contributes to allelic diversity in maize protein abundance

SignificanceProteins are the machinery which execute essential cellular functions. However, measuring their abundance within an organism can be difficult and resource-intensive. Cells use a variety of mechanisms to control protein synthesis from mRNA, including short open reading frames (uORFs) that lie upstream of the main coding sequence. Ribosomes can preferentially translate uORFs instead of the main coding sequence, leading to reduced translation of the main protein. In this study, we show that uORF sequence variation between individuals can lead to different rates of protein translation and thus variable protein abundances. We also demonstrate that natural variation in uORFs occurs frequently and can be linked to whole-plant phenotypes, indicating that uORF sequence variation likely contributes to plant adaptation.

The Codon Usage Bias Analysis of Free-Living Ciliates' Macronuclear Genomes and Clustered Regularly Interspaced Short Palindromic Repeats/Cas9 Vector Construction of Stylonychia lemnae

Ciliates represent higher unicellular animals, and several species are also important model organisms for molecular biology research. Analyses of codon usage bias (CUB) of the macronuclear (MAC) genome in ciliates can not only promote a better understanding of the genetic mode and evolution history of these organisms but also help optimize codons to improve the gene editing efficiency of model ciliates. In this study, macronuclear genome sequences of nine free-living ciliates were analyzed with CodonW software to calculate the following indices: the guanine-cytosine content (GC); the frequency of the nucleotides U, C, A, and G at the third position of codons (U3s, C3s, A3s, G3s); the effective number of codons (ENC); the correlation between GC at the first and second positions (GC12); the frequency of the nucleotides G + C at the third position of synonymous codons (GC3s); the relative synonymous codon usage (RSCU). Parity rule 2 plot analysis, neutrality plot analysis, and correlation analysis were performed to explore the factors that influence codon preference. The results showed that the GC contents in nine ciliates' MAC genomes were lower than 50% and appeared AT-rich. The base compositions of GC12 and GC3s are markedly distinct and the codon usage pattern and evolution of ciliates are affected by genetic mutation and natural selection. According to the synonymous codon analysis, the codons of most ciliates ended with A or U and eight codons were the general optimal codons of nine ciliates. A clustered regularly interspaced short palindromic repeats/Cas9 (CRISPR/Cas9) expression vector of Stylonychia lemnae was constructed by optimizing the macronuclear genome codon and was successfully used to knock out the Adss gene. This is the first such extensive investigation of the MAC genome CUB of ciliates and the initial successful application of the CRISPR/Cas9 technique in free-living ciliates.

Several different sequences are implicated in bloodstream-form-specific gene expression in Trypanosoma brucei

The parasite Trypanosoma brucei grows as bloodstream forms in mammalian hosts, and as procyclic forms in tsetse flies. In trypanosomes, gene expression regulation depends heavily on post-transcriptional mechanisms. Both the RNA-binding protein RBP10 and glycosomal phosphoglycerate kinase PGKC are expressed only in mammalian-infective forms. RBP10 targets procyclic-specific mRNAs for destruction, while PGKC is required for bloodstream-form glycolysis. Developmental regulation of both is essential: expression of either RBP10 or PGKC in procyclic forms inhibits their proliferation. We show that the 3’-untranslated region of the RBP10 mRNA is extraordinarily long—7.3kb—and were able to identify six different sequences, scattered across the untranslated region, which can independently cause bloodstream-form-specific expression. The 3’-untranslated region of the PGKC mRNA, although much shorter, still contains two different regions, of 125 and 153nt, that independently gave developmental regulation. No short consensus sequences were identified that were enriched either within these regulatory regions, or when compared with other mRNAs with similar regulation, suggesting that more than one regulatory RNA-binding protein is important for repression of mRNAs in procyclic forms. We also identified regions, including an AU repeat, that increased expression in bloodstream forms, or suppressed it in both forms. Trypanosome mRNAs that encode RNA-binding proteins often have extremely extended 3’-untranslated regions. We suggest that one function of this might be to act as a fail-safe mechanism to ensure correct regulation even if mRNA processing or expression of trans regulators is defective.

The parasite Trypanosoma brucei causes sleeping sickness in humans, and nagana in cattle, and is transmitted by Tsetse flies. It grows in the bloodstream and tissue fluids of mammalian hosts, as "bloodstream forms", and as "procyclic forms" in the midgut of tsetse flies. Several hundred proteins are expressed in a stage-specific fashion, and this is essential for parasite survival in the different environments. RBP10 is an RNA-binding protein that is expressed only in bloodstream forms. It binds to procyclic-specific mRNAs, and causes their destruction. PGKC is an enzyme that is also specifically expressed in bloodstream forms. Developmental regulation of both is essential: expression of either RBP10 or PGKC in procyclic forms prevents their growth. The mRNAs encoding both proteins are very unstable in procyclic forms, and the sequences responsible are in an "untranslated region" of the mRNA—sequences that follow the part that codes for protein. We here show that the mRNA encoding PGKC has two regions that independently cause developmental regulation, and that the very long untranslated region of the RBP10 mRNA has no fewer than six regulatory regions, but there were no obvious similarities between them. We suggest that the presence of several different regulatory sequences in trypanosome mRNAs might be a fail-safe mechanism to ensure correct regulation.

Characterization of ADAT2/3 molecules in Trypanosoma cruzi and regulation of mucin gene expression by tRNA editing

Adenosine-to-inosine conversion at position 34 (A34-to-I) of certain tRNAs is essential for expanding their decoding capacity. This reaction is catalyzed by the adenosine deaminase acting on tRNA (ADAT) complex, which in Eukarya is formed by two subunits: ADAT2 and ADAT3. We herein identified and thoroughly characterized the ADAT molecules from the protozoan pathogen Trypanosoma cruzi, the causative agent of Chagas Disease. TcADAT2 and TcADAT3 spontaneously form a catalytically active complex, as shown by expression in engineered bacteria and/or by the increased ex vivo tRNA A-to-I deamination activity of T. cruzi epimastigotes overexpressing TcADAT subunits. Importantly, enhanced TcADAT2/3 activity in transgenic parasites caused a shift in their in vivo tRNAThrAGU signature, which correlated with significant changes in the expression of the Thr-rich TcSMUG proteins. To our knowledge, this is the first evidence indicating that T. cruzi tRNA editing can be modulated in vivo, in turn post-transcriptionally changing the expression of specific genes. Our findings suggest tRNA editing/availability as a forcible step in controlling gene expression and driving codon adaptation in T. cruzi. Moreover, we unveil certain differences between parasite and mammalian host tRNA editing and processing, such as cytosine-to-uridine conversion at position 32 of tRNAThrAGU in T. cruzi, that may be exploited for the identification of novel druggable targets of intervention.

Identification of Oocyst-Driven Toxoplasma gondii Infections in Humans and Animals through Stage-Specific Serology-Current Status and Future Perspectives

The apicomplexan zoonotic parasite Toxoplasma gondii has three infective stages: sporozoites in sporulated oocysts, which are shed in unsporulated form into the environment by infected felids; tissue cysts containing bradyzoites, and fast replicating tachyzoites that are responsible for acute toxoplasmosis. The contribution of oocysts to infections in both humans and animals is understudied despite being highly relevant. Only a few diagnostic antigens have been described to be capable of discriminating which parasite stage has caused an infection. Here we provide an extensive overview of the antigens and serological assays used to detect oocyst-driven infections in humans and animals according to the literature. In addition, we critically discuss the possibility to exploit the increasing knowledge of the T. gondii genome and the various 'omics datasets available, by applying predictive algorithms, for the identification of new oocyst-specific proteins for diagnostic purposes. Finally, we propose a workflow for how such antigens and assays based on them should be evaluated to ensure reproducible and robust results.

In Vivo Structure-Function Analysis and Redox Interactomes of Leishmania tarentolae Erv

Import and oxidative folding of proteins in the mitochondrial intermembrane space differ among eukaryotic lineages. While opisthokonts such as yeast rely on the receptor and oxidoreductase Mia40 in combination with the Mia40:cytochrome c oxidoreductase Erv, kinetoplastid parasites and other Excavata/Discoba lack Mia40 but have a functional Erv homologue. Whether excavate Erv homologues rely on a Mia40 replacement or directly interact with imported protein substrates remains controversial. Here, we used the CRISPR-Cas9 system to generate a set of tagged and untagged homozygous mutants of LTERV from the kinetoplastid model parasite Leishmania tarentolae. Modifications of the shuttle cysteine motif of LtErv were lethal, whereas replacement of clamp residue Cys¹⁷ or removal of the kinetoplastida-specific second (KISS) domain had no impact on parasite viability under standard growth conditions. However, removal of the KISS domain rendered parasites sensitive to heat stress and led to the accumulation of homodimeric and mixed LtErv disulfides. We therefore determined and compared the redox interactomes of tagged wild-type LtErv and LtErv^ΔKISS using stable isotope labeling by amino acids in cell culture (SILAC) and quantitative mass spectrometry. While the Mia40-replacement candidate Mic20 and all but one typical substrate with twin Cx_3/9C-motifs were absent in both redox interactomes, we identified a small set of alternative potential interaction partners with putative redox-active cysteine residues. In summary, our study reveals parasite-specific intracellular structure-function relationships and redox interactomes of LtErv with implications for current hypotheses on mitochondrial protein import in nonopisthokonts. IMPORTANCE The discovery of the redox proteins Mia40/CHCHD4 and Erv1/ALR, as well as the elucidation of their relevance for oxidative protein folding in the mitochondrial intermembrane space of yeast and mammals, founded a new research topic in redox biology and mitochondrial protein import. The lack of Mia40/CHCHD4 in protist lineages raises fundamental and controversial questions regarding the conservation and evolution of this essential pathway. Do protist Erv homologues act alone, or do they use the candidate Mic20 or another protein as a Mia40 replacement? Furthermore, we previously showed that Erv homologues in L. tarentolae and the human pathogen L. infantum are not only essential but also differ structurally and mechanistically from yeast and human Erv1/ALR. Here, we analyzed the relevance of such structural differences in vivo and determined the first redox interactomes of a nonopisthokont Erv homologue. Our data challenge recent hypotheses on mitochondrial protein import in nonopisthokonts.

Differential regulation of mRNA fate by the human Ccr4-Not complex is driven by coding sequence composition and mRNA localization

BACKGROUND

Regulation of protein output at the level of translation allows for a rapid adaptation to dynamic changes to the cell's requirements. This precise control of gene expression is achieved by complex and interlinked biochemical processes that modulate both the protein synthesis rate and stability of each individual mRNA. A major factor coordinating this regulation is the Ccr4-Not complex. Despite playing a role in most stages of the mRNA life cycle, no attempt has been made to take a global integrated view of how the Ccr4-Not complex affects gene expression.

RESULTS

This study has taken a comprehensive approach to investigate post-transcriptional regulation mediated by the Ccr4-Not complex assessing steady-state mRNA levels, ribosome position, mRNA stability, and protein production transcriptome-wide. Depletion of the scaffold protein CNOT1 results in a global upregulation of mRNA stability and the preferential stabilization of mRNAs enriched for G/C-ending codons. We also uncover that mRNAs targeted to the ER for their translation have reduced translational efficiency when CNOT1 is depleted, specifically downstream of the signal sequence cleavage site. In contrast, translationally upregulated mRNAs are normally localized in p-bodies, contain disorder-promoting amino acids, and encode nuclear localized proteins. Finally, we identify ribosome pause sites that are resolved or induced by the depletion of CNOT1.

CONCLUSIONS

We define the key mRNA features that determine how the human Ccr4-Not complex differentially regulates mRNA fate and protein synthesis through a mechanism linked to codon composition, amino acid usage, and mRNA localization.

Quantitative elucidation of associations between nucleotide identity and physicochemical properties of amino acids and the functional insight

Studies on codon property would deepen our understanding of the origin of primitive life and enlighten biotechnical application. Here, we proposed a quantitative measurement of codon-amino acid association and found that seven out of 13 physicochemical properties have stronger associations with the nucleotide identity at the second codon position, indicating that protein structure and function may associate more closely with it than the other two sites. When extending the effect of codon-amino acid association to protein level, it was found that the correlation between the second codon position (measured by the relative frequencies of nucleobase T and A at this codon site) and hydrophobicity (by the form of GRAVY value) became stronger with 96% genomes having R > 0.90 and p < 1e-60. Furthermore, we revealed that informational genes encoding proteins have lower GRAVY values than operational proteins (p < 3e-37) in both prokaryotic and eukaryotic genomes. The above results reveal a complete link from codon identity (A₂ versus T₂) to amino acid property (hydrophilic versus hydrophobic) and then to protein functions (informational versus operational). Hence, our work may help to understand how the nucleotide sequence determines protein function.

Natural genetic variation in Drosophila melanogaster reveals genes associated with Coxiella burnetii infection

The gram-negative bacterium Coxiella burnetii is the causative agent of Query (Q) fever in humans and coxiellosis in livestock. Host genetics are associated with C. burnetii pathogenesis both in humans and animals; however, it remains unknown if specific genes are associated with severity of infection. We employed the Drosophila Genetics Reference Panel to perform a genome-wide association study to identify host genetic variants that affect host survival to C. burnetii infection. The genome-wide association study identified 64 unique variants (P < 10⁻⁵) associated with 25 candidate genes. We examined the role each candidate gene contributes to host survival during C. burnetii infection using flies carrying a null mutation or RNAi knockdown of each candidate. We validated 15 of the 25 candidate genes using at least one method. This is the first report establishing involvement of many of these genes or their homologs with C. burnetii susceptibility in any system. Among the validated genes, FER and tara play roles in the JAK/STAT, JNK, and decapentaplegic/TGF-β signaling pathways which are components of known innate immune responses to C. burnetii infection. CG42673 and DIP-ε play roles in bacterial infection and synaptic signaling but have no previous association with C. burnetii pathogenesis. Furthermore, since the mammalian ortholog of CG13404 (PLGRKT) is an important regulator of macrophage function, CG13404 could play a role in host susceptibility to C. burnetii through hemocyte regulation. These insights provide a foundation for further investigation regarding the genetics of C. burnetii susceptibility across a wide variety of hosts.

Genome-wide role of codon usage on transcription and identification of potential regulators

Codon usage bias is a fundamental feature of all genomes and plays an important role in determining gene expression levels. The codon usage was thought to influence gene expression mainly due to its impact on translation. Recently, however, codon usage was shown to affect transcription of fungal and mammalian genes, indicating the existence of a gene regulatory phenomenon with unknown mechanism. In Neurospora, codon usage biases strongly correlate with mRNA levels genome-wide, and here we show that the correlation between codon usage and RNA levels is maintained in the nucleus. In addition, codon optimality is tightly correlated with both total and nuclear RNA levels, suggesting that codon usage broadly influences mRNA levels through transcription in a translation-independent manner. A large-scale RNA sequencing-based genetic screen in Neurospora identified 18 candidate factors that when deleted decreased the genome-wide correlation between codon usage and RNA levels and reduced the codon usage effect on gene expression. Most of these factors, such as the H3K36 methyltransferase, are chromatin regulators or transcription factors. Together, our results suggest that the transcriptional effect of codon usage is mediated by multiple transcriptional regulatory mechanisms.

tRNA 3' shortening by LCCR4 as a response to stress in Trypanosoma brucei

Sensing of environmental cues is crucial for cell survival. To adapt to changes in their surroundings cells need to tightly control the repertoire of genes expressed at any time. Regulation of translation is key, especially in organisms in which transcription is hardly controlled, like Trypanosoma brucei. In this study, we describe the shortening of the bulk of the cellular tRNAs during stress at the expense of the conserved 3' CCA-tail. This tRNA shortening is specific for nutritional stress and renders tRNAs unsuitable substrates for translation. We uncovered the nuclease LCCR4 (Tb927.4.2430), a homologue of the conserved deadenylase Ccr4, as being responsible for tRNA trimming. Once optimal growth conditions are restored tRNAs are rapidly repaired by the trypanosome tRNA nucleotidyltransferase thus rendering the recycled tRNAs amenable for translation. This mechanism represents a fast and efficient way to repress translation during stress, allowing quick reactivation with a low energy input.

The effects of codon bias and optimality on mRNA and protein regulation

The central dogma of molecular biology entails that genetic information is transferred from nucleic acid to proteins. Notwithstanding retro-transcribing genetic elements, DNA is transcribed to RNA which in turn is translated into proteins. Recent advancements have shown that each stage is regulated to control protein abundances for a variety of essential physiological processes. In this regard, mRNA regulation is essential in fine-tuning or calibrating protein abundances. In this review, we would like to discuss one of several mRNA-intrinsic features of mRNA regulation that has been gaining traction of recent-codon bias and optimality. Specifically, we address the effects of codon bias with regard to codon optimality in several biological processes centred on translation, such as mRNA stability and protein folding among others. Finally, we examine how different organisms or cell types, through this system, are able to coordinate physiological pathways to respond to a variety of stress or growth conditions.

Effects of codon usage on gene expression are promoter context dependent

Codon usage bias is a universal feature of all genomes. Although codon usage has been shown to regulate mRNA and protein levels by influencing mRNA decay and transcription in eukaryotes, little or no genome-wide correlations between codon usage and mRNA levels are detected in mammalian cells, raising doubt on the significance of codon usage effect on gene expression. Here we show that gene-specific regulation reduces the genome-wide codon usage and mRNA correlations: Constitutively expressed genes exhibit much higher genome-wide correlations than differentially expressed genes from fungi to human cells. Using Drosophila S2 cells as a model system, we showed that the effect of codon usage on mRNA expression level is promoter-dependent. Regions downstream of the core promoters of differentially expressed genes can repress the codon usage effects on mRNA expression. An element in the Hsp70 promoter was identified to be necessary and sufficient for this inhibitory effect. The promoter-dependent codon usage effects on mRNA levels are regulated at the transcriptional level through modulation of histone modifications, nucleosome densities and premature termination. Together, our results demonstrate that promoters play a major role in determining whether codon usage influences gene expression and further establish the transcription-dependent codon usage effects on gene expression.

Modification of Transfer RNA Levels Affects Cyclin Aggregation and the Correct Duplication of Yeast Cells

Codon usage bias (the preferential use of certain synonymous codons (optimal) over others is found at the organism level (intergenomic) within specific genomes (intragenomic) and even in certain genes. Whether it is the result of genetic drift due to GC/AT content and/or natural selection is a topic of intense debate. Preferential codons are mostly found in genes encoding highly-expressed proteins, while lowly-expressed proteins usually contain a high proportion of rare (lowly-represented) codons. While optimal codons are decoded by highly expressed tRNAs, rare codons are usually decoded by lowly-represented tRNAs. Whether rare codons play a role in controlling the expression of lowly- or temporarily-expressed proteins is an open question. In this work we approached this question using two strategies, either by replacing rare glycine codons with optimal counterparts in the gene that encodes the cell cycle protein Cdc13, or by overexpression the tRNA ^Gly that decodes rare codons from the fission yeast, Schizosaccharomyces pombe. While the replacement of synonymous codons severely affected cell growth, increasing tRNA levels affected the aggregation status of Cdc13 and cell division. These lead us to think that rare codons in lowly-expressed cyclin proteins are crucial for cell division, and that the overexpression of tRNA that decodes rare codons affects the expression of proteins containing these rare codons. These codons may be the result of the natural selection of codons in genes that encode lowly-expressed proteins.

Synonymous but Not Silent: The Codon Usage Code for Gene Expression and Protein Folding

Codon usage bias, the preference for certain synonymous codons, is found in all genomes. Although synonymous mutations were previously thought to be silent, a large body of evidence has demonstrated that codon usage can play major roles in determining gene expression levels and protein structures. Codon usage influences translation elongation speed and regulates translation efficiency and accuracy. Adaptation of codon usage to tRNA expression determines the proteome landscape. In addition, codon usage biases result in nonuniform ribosome decoding rates on mRNAs, which in turn influence the cotranslational protein folding process that is critical for protein function in diverse biological processes. Conserved genome-wide correlations have also been found between codon usage and protein structures. Furthermore, codon usage is a major determinant of mRNA levels through translation-dependent effects on mRNA decay and translation-independent effects on transcriptional and posttranscriptional processes. Here, we discuss the multifaceted roles and mechanisms of codon usage in different gene regulatory processes.

Identification of Novel Interspersed DNA Repetitive Elements in the Trypanosoma cruzi Genome Associated with the 3'UTRs of Surface Multigenic Families

Trypanosoma cruzi is the etiological agent of Chagas disease, which affects millions of people in Latin America. No transcriptional control of gene expression has been demonstrated in this organism, and 50% of its genome consists of repetitive elements and members of multigenic families. In this study, we applied a novel bioinformatics approach to predict new repetitive elements in the genome sequence of T. cruzi. A new repetitive sequence measuring 241 nt was identified and found to be interspersed along the genome sequence from strains of different DTUs. This new repeat was mostly on intergenic regions, and upstream and downstream regions of the 241 nt repeat were enriched in surface protein genes. RNAseq analysis revealed that the repeat was part of processed mRNAs and was predominantly found in the 3′ untranslated regions (UTRs) of genes of multigenic families encoding surface proteins. Moreover, we detected a correlation between the presence of the repeat in the 3′UTR of multigenic family genes and the level of differential expression of these genes when comparing epimastigote and trypomastigote transcriptomes. These data suggest that this sequence plays a role in the posttranscriptional regulation of the expression of multigenic families.

Nonoptimal Codon Usage Is Critical for Protein Structure and Function of the Master General Amino Acid Control Regulator CPC-1

Under amino acid starvation conditions, eukaryotic organisms activate a general amino acid control response. In Neurospora crassa, Cross Pathway Control Protein 1 (CPC-1), the ortholog of the Saccharomyces cerevisiae bZIP transcription factor GCN4, functions as the master regulator of the general amino acid control response. Codon usage biases are a universal feature of eukaryotic genomes and are critical for regulation of gene expression. Although codon usage has also been implicated in the regulation of protein structure and function, genetic evidence supporting this conclusion is very limited. Here, we show that Neurospora cpc-1 has a nonoptimal NNU-rich codon usage profile that contrasts with the strong NNC codon preference in the genome. Although substitution of the cpc-1 NNU codons with synonymous NNC codons elevated CPC-1 expression in Neurospora, it altered the CPC-1 degradation rate and abolished its amino acid starvation-induced protein stabilization. The codon-manipulated CPC-1 protein also exhibited different sensitivity to limited protease digestion. Furthermore, CPC-1 functions in rescuing the cell growth of the cpc-1 deletion mutant and activation of the expression of its target genes were impaired by the synonymous codon changes. Together, these results reveal the critical role of codon usage in regulation of CPC-1 expression and function and establish a genetic example of the importance of codon usage in protein folding.IMPORTANCE The general amino acid control response is critical for adaptation of organisms to amino acid starvation conditions. The preference to use certain synonymous codons is a universal feature of all genomes. Synonymous codon changes were previously thought to be silent mutations. In this study, we showed that the Neurospora cpc-1 gene has an unusual codon usage profile compared to other genes in the genome. We found that codon optimization of the cpc-1 gene without changing its amino acid sequence resulted in elevated CPC-1 expression, an altered protein degradation rate, and impaired protein functions due to changes in protein structure. Together, these results reveal the critical role of synonymous codon usage in regulation of CPC-1 expression and function and establish a genetic example of the importance of codon usage in protein structure.

The Ongoing Quest to Crack the Genetic Code for Protein Production

Understanding the genetic design principles that determine protein production remains a major challenge. Although the key principles of gene expression were discovered 50 years ago, additional factors are still being uncovered. Both protein-coding and non-coding sequences harbor elements that collectively influence the efficiency of protein production by modulating transcription, mRNA decay, and translation. The influences of many contributing elements are intertwined, which complicates a full understanding of the individual factors. In natural genes, a functional balance between these factors has been obtained in the course of evolution, whereas for genetic-engineering projects, our incomplete understanding still limits optimal design of synthetic genes. However, notable advances have recently been made, supported by high-throughput analysis of synthetic gene libraries as well as by state-of-the-art biomolecular techniques. We discuss here how these advances further strengthen understanding of the gene expression process and how they can be harnessed to optimize protein production.

A code within the genetic code: codon usage regulates co-translational protein folding

The genetic code is degenerate, and most amino acids are encoded by two to six synonymous codons. Codon usage bias, the preference for certain synonymous codons, is a universal feature of all genomes examined. Synonymous codon mutations were previously thought to be silent; however, a growing body evidence now shows that codon usage regulates protein structure and gene expression through effects on co-translational protein folding, translation efficiency and accuracy, mRNA stability, and transcription. Codon usage regulates the speed of translation elongation, resulting in non-uniform ribosome decoding rates on mRNAs during translation that is adapted to co-translational protein folding process. Biochemical and genetic evidence demonstrate that codon usage plays an important role in regulating protein folding and function in both prokaryotic and eukaryotic organisms. Certain protein structural types are more sensitive than others to the effects of codon usage on protein folding, and predicted intrinsically disordered domains are more prone to misfolding caused by codon usage changes than other domain types. Bioinformatic analyses revealed that gene codon usage correlates with different protein structures in diverse organisms, indicating the existence of a codon usage code for co-translational protein folding. This review focuses on recent literature on the role and mechanism of codon usage in regulating translation kinetics and co-translational protein folding.

Interplay between Position-Dependent Codon Usage Bias and Hydrogen Bonding at the 5' End of ORFeomes

Codon usage bias exerts control over a wide variety of molecular processes. The positioning of synonymous codons within coding sequences (CDSs) dictates protein expression by mechanisms such as local translation efficiency, mRNA Gibbs free energy, and protein cotranslational folding. In this work, we explore how codon usage affects the position-dependent content of hydrogen bonding, which in turn influences energy requirements for unwinding double-stranded DNA (dsDNA). We categorized codons according to their hydrogen bond content and found differential effects on hydrogen bonding encoded by codon variants. The specific positional disposition of codon variants within CDSs creates a ramp of hydrogen bonding at the 5' end of the ORFeome in Escherichia coli CDSs occupying the first position of operons are subjected to selective pressure that reduces their hydrogen bonding compared to internal CDSs, and highly transcribed CDSs demand a lower maximum capacity of hydrogen bonds per codon, suggesting that the energetic requirement for unwinding the dsDNA in highly transcribed CDSs has evolved to be minimized in E. coli Subsequent analysis of over 14,000 ORFeomes showed a pervasive ramp of hydrogen bonding at the 5' end in Bacteria and Archaea that positively correlates with the probability of mRNA secondary structure formation. Both the ramp and the correlation were not found in Fungi The position-dependent hydrogen bonding might be part of the mechanism that contributes to the coordination between transcription and translation in Bacteria and Archaea A Web-based application to analyze the position-dependent hydrogen bonding of ORFeomes has been developed and is publicly available (https://juanvillada.shinyapps.io/hbonds/).IMPORTANCE Redundancy of the genetic code creates a vast space of alternatives to encode a protein. Synonymous codons exert control over a variety of molecular and physiological processes of cells mainly through influencing protein biosynthesis. Recent findings have shown that synonymous codon choice affects transcription by controlling mRNA abundance, mRNA stability, transcription termination, and transcript biosynthesis cost. In this work, by analyzing thousands of Bacteria, Archaea, and Fungi genomes, we extend recent findings by showing that synonymous codon choice, corresponding to the number of hydrogen bonds in a codon, can also have an effect on the energetic requirements for unwinding double-stranded DNA in a position-dependent fashion. This report offers new perspectives on the mechanism behind the transcription-translation coordination and complements previous hypotheses on the resource allocation strategies used by Bacteria and Archaea to manage energy efficiency in gene expression.

EIF2α phosphorylation is regulated in intracellular amastigotes for the generation of infective Trypanosoma cruzi trypomastigote forms

Trypanosomatids regulate gene expression mainly at the post-transcriptional level through processing, exporting and stabilising mRNA and control of translation. In most eukaryotes, protein synthesis is regulated by phosphorylation of eukaryotic initiation factor 2 (eIF2) at serine 51. Phosphorylation halts overall translation by decreasing availability of initiator tRNA^met to form translating ribosomes. In trypanosomatids, the N-terminus of eIF2α is extended with threonine 169 the homologous phosphorylated residue. Here, we evaluated whether eIF2α phosphorylation varies during the Trypanosoma cruzi life cycle, the etiological agent of Chagas' disease. Total levels of eIF2α are diminished in infective and non-replicative trypomastigotes compared with proliferative forms from the intestine of the insect vector or amastigotes from mammalian cells, consistent with decreased protein synthesis reported in infective forms. eIF2α phosphorylation increases in proliferative intracellular forms prior to differentiation into trypomastigotes. Parasites overexpressing eIF2α^T169A or with an endogenous CRISPR/Cas9-generated eIF2α^T169A mutation were created and analysis revealed alterations to the proteome, largely unrelated to the presence of μORF in epimastigotes. eIF2α^T169A mutant parasites produced fewer trypomastigotes with lower infectivity than wild type, with increased levels of sialylated mucins and oligomannose glycoproteins, and decreased galactofuranose epitopes and the surface protease GP63 on the cell surface. We conclude that eIF2α expression and phosphorylation levels affect proteins relevant for intracellular progression of T. cruzi.

The RNA-associated proteins MKT1 and MKT1L form alternative PBP1-containing complexes in Trypanosoma brucei

Control of gene expression in kinetoplastids such as trypanosomes depends heavily on RNA-binding proteins that influence mRNA decay and translation. We previously showed that the trypanosome protein MKT1 forms a multicomponent protein complex: MKT1 interacts with PBP1, which in turn recruits LSM12 and poly(A)-binding protein. MKT1 is recruited to mRNAs by sequence-specific RNA-binding proteins, resulting in stabilization of the bound mRNA. We here show that PBP1, LSM12, and a 117-residue protein, XAC1 (Tb927.7.2780), are present in complexes that contain either MKT1 or an MKT1-like protein, MKT1L (Tb927.10.1490). All five proteins are present predominantly in the complexes, and we found evidence for a minor subset of complexes containing both MKT1 and MKT1L. XAC1-containing complexes reproducibly contained RNA-binding proteins that were previously found associated with MKT1. Moreover, XAC1- or MKT1-containing complexes specifically recruited one of the two poly(A)-binding proteins, PABP2, and one of the six cap-binding translation initiation complexes, EIF4E6-EIF4G5. Yeast two-hybrid assay results indicated that MKT1 directly interacts with EIF4G5. MKT1-PBP1 complexes can therefore interact with mRNAs via their poly(A) tails and caps, as well as through sequence-specific RNA-binding proteins. Correspondingly, MKT1 is associated with many mRNAs, although not with those encoding ribosomal proteins. Meanwhile, MKT1L resembles MKT1 at the C terminus but additionally features an N-terminal extension with low-complexity regions. Although MKT1L depletion inhibited cell proliferation, we found no evidence that it specifically interacts with RNA-binding proteins or mRNA. We speculate that MKT1L may compete with MKT1 for PBP1 binding and thereby modulate the function of MKT1-containing complexes.

Adaptation of codon usage to tRNA I34 modification controls translation kinetics and proteome landscape

Codon usage bias is a universal feature of all genomes and plays an important role in regulating protein expression levels. Modification of adenosine to inosine at the tRNA anticodon wobble position (I34) by adenosine deaminases (ADATs) is observed in all eukaryotes and has been proposed to explain the correlation between codon usage and tRNA pool. However, how the tRNA pool is affected by I34 modification to influence codon usage-dependent gene expression is unclear. Using Neurospora crassa as a model system, by combining molecular, biochemical and bioinformatics analyses, we show that silencing of adat2 expression severely impaired the I34 modification levels for the ADAT-related tRNAs, resulting in major ADAT-related tRNA profile changes and reprogramming of translation elongation kinetics on ADAT-related codons. adat2 silencing also caused genome-wide codon usage-biased ribosome pausing on mRNAs and proteome landscape changes, leading to selective translational repression or induction of different mRNAs. The induced expression of CPC-1, the Neurospora ortholog of yeast GCN4p, mediates the transcriptional response after adat2 silencing and amino acid starvation. Together, our results demonstrate that the tRNA I34 modification by ADAT plays a major role in driving codon usage-biased translation to shape proteome landscape.

Modification of transfer RNA (tRNA) can have profound impacts on gene expression by shaping cellular tRNA pool. How codon usage bias and tRNA profiles synergistically regulate gene expression is unclear. By combining molecular, biochemical and bioinformatics analyses, we showed that the correlation between genome codon usage and tRNA I34 (inosine 34) modification modulates translation elongation kinetics and proteome landscape. Inhibition of tRNA I34 modification causes codon usage-dependent ribosome pausing on mRNAs during translation and changes cellular protein contents in a codon usage biased manner. Together, our results demonstrate that the tRNA I34 modification plays a major role in driving codon usage-dependent translation to determine proteome landscape in a eukaryotic organism.

GCN2-Like Kinase Modulates Stress Granule Formation During Nutritional Stress in Trypanosoma cruzi

The integrated stress response in eukaryotic cells is an orchestrated pathway that leads to eukaryotic Initiation Factor 2 alpha subunit (eIF2α) phosphorylation at ser51 and ultimately activates pathways to mitigate cellular damages. Three putative kinases (Tck1, Tck2, and Tck3) are found in the Trypanosoma cruzi genome, the flagellated parasite that causes Chagas disease. These kinases present similarities to other eukaryotic eIF2α kinases, exhibiting a typical insertion loop in the kinase domain of the protein. We found that this insertion loop is conserved among kinase 1 of several T. cruzi strains but differs among various Kinetoplastidae species, suggesting unique roles. Kinase 1 is orthologous of GCN2 of several eukaryotes, which have been implicated in the eIF2α ser51 phosphorylation in situations that mainly affects the nutrients levels. Therefore, we further investigated the responses to nutritional stress of T. cruzi devoid of TcK1 generated by CRISPR/Cas9 gene replacement. In nutrient-rich conditions, replicative T. cruzi epimastigotes depleted of TcK1 proliferate as wild type cells but showed increased levels of polysomes relative to monosomes. Upon nutritional deprivation, the polysomes decreased more than in TcK1 depleted line. However, eIF2α is still phosphorylated in TcK1 depleted line, as in wild type parasites. eIF2α phosphorylation increased at longer incubations times, but KO parasites showed less accumulation of ribonucleoprotein granules containing ATP-dependent RNA helicase involved in mRNA turnover (DHH1) and Poly-A binding protein (PABP1). Additionally, the formation of metacyclic-trypomastigotes is increased in the absence of Tck1 compared to controls. These metacyclics, as well as tissue culture trypomastigotes derived from the TcK1 knockout line, were less infective to mammalian host cells, although replicated faster inside mammalian cells. These results indicate that GCN2-like kinase in T. cruzi affects stress granule formation, independently of eIF2α phosphorylation upon nutrient deprivation. It also modulates the fate of the parasites during differentiation, invasion, and intracellular proliferation.

Genome-wide codon usage pattern analysis reveals the correlation between codon usage bias and gene expression in Cuscuta australis

The protein-coding genes and pseudogenes of Cuscuta australis had the diverse contribution to the formation and evolution of parasitism. The codon usage pattern analysis of these two type genes could be used to understand the gene transcription and translation. In this study, we systematically analyzed the codon usage patterns of protein-coding sequences and pseudogenes sequences in C. australis. The results showed that the high frequency codons of protein coding sequences and pseudogenes had the same A/U bias in the third position. However, these two sequences had converse bias at the third base in optimal codons: the protein coding sequences preferred G/C-ending codons while pseudogene sequences preferred A/U-ending codons. Neutrality plot and effective number of codons plot revealed that natural selection played a more important role than mutation pressure in two sequences codon usage bias. Furthermore, the gene expression level had a significant positive correlation with codon usage bias in C. australis. Highly-expressed protein coding genes exhibited a higher codon bias than lowly-expressed genes. Meanwhile, the high-expression genes tended to use G/C-ending synonymous codons. This result further verified the optimal codons usage bias and its correlation with the gene expression in C. australis.

The endoplasmic reticulum-associated mRNA-binding proteins ERBP1 and ERBP2 interact in bloodstream-form Trypanosoma brucei

Kinetoplastids rely heavily on post-transcriptional mechanisms for control of gene expression, and on RNA-binding proteins that regulate mRNA splicing, translation and decay. Trypanosoma brucei ERBP1 (Tb927.10.14150) and ERBP2 (Tb927.9.9550) were previously identified as mRNA binding proteins that lack canonical RNA-binding domains. We show here that ERBP1 is associated with the endoplasmic reticulum, like ERBP2, and that the two proteins interact in vivo. Loss of ERBP1 from bloodstream-form T. brucei initially resulted in a growth defect but proliferation was restored after more prolonged cultivation. Pull-down analysis of tagged ERBP1 suggests that it preferentially binds to ribosomal protein mRNAs. The ERBP1 sequence resembles that of Saccharomyces cerevisiae Bfr1, which also localises to the endoplasmic reticulum and binds to ribosomal protein mRNAs. However, unlike Bfr1, ERBP1 does not bind to mRNAs encoding secreted proteins, and it is also not recruited to stress granules after starvation.

eRF1 mediates codon usage effects on mRNA translation efficiency through premature termination at rare codons

Codon usage bias is a universal feature of eukaryotic and prokaryotic genomes and plays an important role in regulating gene expression levels. A major role of codon usage is thought to regulate protein expression levels by affecting mRNA translation efficiency, but the underlying mechanism is unclear. By analyzing ribosome profiling results, here we showed that codon usage regulates translation elongation rate and that rare codons are decoded more slowly than common codons in all codon families in Neurospora. Rare codons resulted in ribosome stalling in manners both dependent and independent of protein sequence context and caused premature translation termination. This mechanism was shown to be conserved in Drosophila cells. In both Neurospora and Drosophila cells, codon usage plays an important role in regulating mRNA translation efficiency. We found that the rare codon-dependent premature termination is mediated by the translation termination factor eRF1, which recognizes ribosomes stalled on rare sense codons. Silencing of eRF1 expression resulted in codon usage-dependent changes in protein expression. Together, these results establish a mechanism for how codon usage regulates mRNA translation efficiency.

Coding regions affect mRNA stability in human cells

A new paradigm has emerged that coding regions can regulate mRNA stability in model organisms. Here, due to differences in cognate tRNA abundance, synonymous codons are translated at different speeds, and slow codons then stimulate mRNA decay. To ask if this phenomenon also occurs in humans, we isolated RNA stability effects due to coding regions using the human ORFeome collection. We find that many open reading frame (ORF) characteristics, such as length and secondary structure, fail to provide explanations for how coding regions alter mRNA stability, and, instead, that the ORF relies on translation to impact mRNA stability. Consistent with what has been seen in other organisms, codon use is related to the effects of ORFs on transcript stability. Importantly, we found instability-associated codons have longer A-site dwell times, suggesting for the first time in humans a connection between elongation speed and mRNA decay. Thus, we propose that codon usage alters decoding speeds and so affects human mRNA stability.

Codon bias confers stability to human mRNAs

Codon bias has been implicated as one of the major factors contributing to mRNA stability in several model organisms. However, the molecular mechanisms of codon bias on mRNA stability remain unclear in humans. Here, we show that human cells possess a mechanism to modulate RNA stability through a unique codon bias. Bioinformatics analysis showed that codons could be clustered into two distinct groups-codons with G or C at the third base position (GC3) and codons with either A or T at the third base position (AT3): the former stabilizing while the latter destabilizing mRNA. Quantification of codon bias showed that increased GC3-content entails proportionately higher GC-content. Through bioinformatics, ribosome profiling, and in vitro analysis, we show that decoupling the effects of codon bias reveals two modes of mRNA regulation, one GC3- and one GC-content dependent. Employing an immunoprecipitation-based strategy, we identify ILF2 and ILF3 as RNA-binding proteins that differentially regulate global mRNA abundances based on codon bias. Our results demonstrate that codon bias is a two-pronged system that governs mRNA abundance.

Efficient and specific oligo-based depletion of rRNA

In most organisms, ribosomal RNA (rRNA) contributes to >85% of total RNA. Thus, to obtain useful information from RNA-sequencing (RNA-seq) analyses at reasonable sequencing depth, typically, mature polyadenylated transcripts are enriched or rRNA molecules are depleted. Targeted depletion of rRNA is particularly useful when studying transcripts lacking a poly(A) tail, such as some non-coding RNAs (ncRNAs), most bacterial RNAs and partially degraded or immature transcripts. While several commercially available kits allow effective rRNA depletion, their efficiency relies on a high degree of sequence homology between oligonucleotide probes and the target RNA. This restricts the use of such kits to a limited number of organisms with conserved rRNA sequences. In this study we describe the use of biotinylated oligos and streptavidin-coated paramagnetic beads for the efficient and specific depletion of trypanosomal rRNA. Our approach reduces the levels of the most abundant rRNA transcripts to less than 5% with minimal off-target effects. By adjusting the sequence of the oligonucleotide probes, our approach can be used to deplete rRNAs or other abundant transcripts independent of species. Thus, our protocol provides a useful alternative for rRNA removal where enrichment of polyadenylated transcripts is not an option and commercial kits for rRNA are not available.

A post-transcriptional respiratome regulon in trypanosomes

Post-transcriptional regulons coordinate the expression of groups of genes in eukaryotic cells, yet relatively few have been characterized. Parasitic trypanosomatids are particularly good models for studies on such mechanisms because they exhibit almost exclusive polycistronic, and unregulated, transcription. Here, we identify the Trypanosoma brucei ZC3H39/40 RNA-binding proteins as regulators of the respiratome; the mitochondrial electron transport chain (complexes I-IV) and the FoF1-ATP synthase (complex V). A high-throughput RNAi screen initially implicated both ZC3H proteins in variant surface glycoprotein (VSG) gene silencing. This link was confirmed and both proteins were shown to form a cytoplasmic ZC3H39/40 complex. Transcriptome and mRNA-interactome analyses indicated that the impact on VSG silencing was indirect, while the ZC3H39/40 complex specifically bound and stabilized transcripts encoding respiratome-complexes. Quantitative proteomic analyses revealed specific positive control of >20 components from complexes I, II and V. Our findings establish a link between the mitochondrial respiratome and VSG gene silencing in bloodstream form T. brucei. They also reveal a major respiratome regulon controlled by the conserved trypanosomatid ZC3H39/40 RNA-binding proteins.

Regulation of gene expression in trypanosomatids: living with polycistronic transcription

In trypanosomes, RNA polymerase II transcription is polycistronic and individual mRNAs are excised by trans-splicing and polyadenylation. The lack of individual gene transcription control is compensated by control of mRNA processing, translation and degradation. Although the basic mechanisms of mRNA decay and translation are evolutionarily conserved, there are also unique aspects, such as the existence of six cap-binding translation initiation factor homologues, a novel decapping enzyme and an mRNA stabilizing complex that is recruited by RNA-binding proteins. High-throughput analyses have identified nearly a hundred regulatory mRNA-binding proteins, making trypanosomes valuable as a model system to investigate post-transcriptional regulation.

African trypanosomes

African trypanosomes cause human African trypanosomiasis and animal African trypanosomiasis. They are transmitted by tsetse flies in sub-Saharan Africa. Although most famous for their mechanisms of immune evasion by antigenic variation, there have been recent important studies that illuminate important aspects of the biology of these parasites both in their mammalian host and during passage through their tsetse fly vector. This Primer overviews current research themes focused on these parasites and discusses how these biological insights and the development of new technologies to interrogate gene function are being used in the search for new approaches to control the parasite. The new insights into the biology of trypanosomes in their host and vector highlight that we are in a ‘golden age’ of discovery for these fascinating parasites.

Translation affects mRNA stability in a codon-dependent manner in human cells

mRNA translation decodes nucleotide into amino acid sequences. However, translation has also been shown to affect mRNA stability depending on codon composition in model organisms, although universality of this mechanism remains unclear. Here, using three independent approaches to measure exogenous and endogenous mRNA decay, we define which codons are associated with stable or unstable mRNAs in human cells. We demonstrate that the regulatory information affecting mRNA stability is encoded in codons and not in nucleotides. Stabilizing codons tend to be associated with higher tRNA levels and higher charged/total tRNA ratios. While mRNAs enriched in destabilizing codons tend to possess shorter poly(A)-tails, the poly(A)-tail is not required for the codon-mediated mRNA stability. This mechanism depends on translation; however, the number of ribosome loads into a mRNA modulates the codon-mediated effects on gene expression. This work provides definitive evidence that translation strongly affects mRNA stability in a codon-dependent manner in human cells.

Proteins are made by joining together building blocks called amino acids into strings. The proteins are ‘translated’ from genetic sequences called mRNA molecules. These sequences can be thought of as series of ‘letters’, which are read in groups of three known as codons. Molecules called tRNAs recognize the codons and add the matching amino acids to the end of the protein. Each tRNA can recognize one or several codons, and the levels of different tRNAs inside the cell vary.

There are 61 codons that code for amino acids, but only 20 amino acids. This means that some codons produce the same amino acid. Despite this, there is evidence to suggest that not all of the codons that produce the same amino acid are exactly equivalent. In bacteria, yeast and zebrafish, some codons seem to make the mRNA molecule more stable, and others make it less stable. This might help the cell to control how many proteins it makes. It was not clear whether the same is true for humans.

To find out, Wu et al. used three separate methods to examine mRNA stability in four types of human cell. Overall, the results revealed that some codons help to stabilize the mRNA, while others make the mRNA molecule break down faster. The effect seems to depend on the supply of tRNAs that have a charged amino acid; mRNA molecules were more likely to self-destruct in cells that contained codons with low levels of the tRNA molecules.

Wu et al. also found that conditions in the cell can alter how strongly the codons affect mRNA stability. For example, a cell that has been infected by a virus reduces translation. Under these conditions, the identity of the codons in the mRNA has less effect on the stability of the mRNA molecule.

Changes to protein production happen in many diseases. Understanding what controls these changes could help to reveal more about our fundamental biology, and what happens when it goes wrong.

Mechanistic Similarities between Antigenic Variation and Antibody Diversification during Trypanosoma brucei Infection

Trypanosoma brucei, which causes African trypanosomiasis, avoids immunity by periodically switching its surface composition. The parasite is coated by 10 million identical, monoallelically expressed variant surface glycoprotein (VSG) molecules. Multiple distinct parasites (with respect to their VSG coat) coexist simultaneously during each wave of parasitemia. This substantial antigenic load is countered by B cells whose antigen receptors (antibodies or immunoglobulins) are also monoallelically expressed, and that diversify dynamically to counter each variant antigen. Here we examine parallels between the processes that generate VSGs and antibodies. We also discuss current insights into VSG mRNA regulation that may inform the emerging field of Ig mRNA biology. We conclude by extending the parallels between VSG and Ig to the protein level.