Translation initiation by non-AUG codons in Arabidopsis thaliana transgenic plants

Translation initiation at AUG and non-AUG triplets in plants

In plants and other eukaryotes, precise selection of translation initiation site (TIS) on mRNAs shapes the proteome in response to cellular events or environmental cues. The canonical translation of mRNAs initiates at a 5' proximal AUG codon in a favorable context. However, the coding and non-coding regions of plant genomes contain numerous unannotated alternative AUG and non-AUG TISs. Determining how and why these unexpected and prevalent TISs are activated in plants has emerged as an exciting research area. In this review, we focus on the selection of plant TISs and highlight studies that revealed previously unannotated TISs used in vivo via comparative genomics and genome-wide profiling of ribosome positioning and protein N-terminal ends. The biological signatures of non-AUG TIS-initiated open reading frames (ORFs) in plants are also discussed. We describe what is understood about cis-regulatory RNA elements and trans-acting eukaryotic initiation factors (eIFs) in the site selection for translation initiation by featuring the findings in plants along with supporting findings in non-plant species. The prevalent, unannotated TISs provide a hidden reservoir of ORFs that likely help reshape plant proteomes in response to developmental or environmental cues. These findings underscore the importance of understanding the mechanistic basis of TIS selection to functionally annotate plant genomes, especially for crops with large genomes.

Chloroplasts evolved an additional layer of translational regulation based on non-AUG start codons for proteins with different turnover rates

Chloroplasts have evolved from photosynthetic cyanobacteria-like progenitors through endosymbiosis. The chloroplasts of present-day land plants have their own transcription and translation systems that show several similarities with prokaryotic organisms. A remarkable feature of the chloroplast translation system is the use of non-AUG start codons in the protein synthesis of certain genes that are evolutionarily conserved from Algae to angiosperms. However, the biological significance of such use of non-AUG codons is not fully understood. The present study was undertaken to unravel the significance of non-AUG start codons in vivo using the chloroplast genetic engineering approach. For this purpose, stable transplastomic tobacco plants expressing a reporter gene i.e. uidA (GUS) under four different start codons (AUG/UUG/GUG/CUG) were generated and β-glucuronidase (GUS) expression was compared. To investigate further the role of promoter sequences proximal to the start codon, uidA was expressed under two different chloroplast gene promoters psbA and psbC that use AUG and a non-AUG (GUG) start codons, respectively, and also showed significant differences in the DNA sequence surrounding the start codon. Further, to delineate the role of RNA editing that creates AUG start codon by editing non-AUG codons, if any, which is another important feature of the chloroplast transcription and translation system, transcripts were sequenced. In addition, a proteomic approach was used to identify the translation initiation site(s) of GUS and the N-terminal amino acid encoded when expressed under different non-AUG start codons. The results showed that chloroplasts use non-AUG start codons in combination with the translation initiation site as an additional layer of gene regulation to over-express proteins that are required at high levels due to their high rates of turnover.

Evaluating Plant Gene Models Using Machine Learning

Gene models are regions of the genome that can be transcribed into RNA and translated to proteins, or belong to a class of non-coding RNA genes. The prediction of gene models is a complex process that can be unreliable, leading to false positive annotations. To help support the calling of confident conserved gene models and minimize false positives arising during gene model prediction we have developed Truegene, a machine learning approach to classify potential low confidence gene models using 14 gene and 41 protein-based characteristics. Amino acid and nucleotide sequence-based features were calculated for conserved (high confidence) and non-conserved (low confidence) annotated genes from the published Pisum sativum Cameor genome. These features were used to train eXtreme Gradient Boost (XGBoost) classifier models to predict whether a gene model is likely to be real. The optimized models demonstrated a prediction accuracy ranging from 87% to 90% and an F-1 score of 0.91–0.94. We used SHapley Additive exPlanations (SHAP) and feature importance plots to identify the features that contribute to the model predictions, and we show that protein and gene-based features can be used to build accurate models for gene prediction that have applications in supporting future gene annotation processes.

Codon usage bias

Codon usage bias is the preferential or non-random use of synonymous codons, a ubiquitous phenomenon observed in bacteria, plants and animals. Different species have consistent and characteristic codon biases. Codon bias varies not only with species, family or group within kingdom, but also between the genes within an organism. Codon usage bias has evolved through mutation, natural selection, and genetic drift in various organisms. Genome composition, GC content, expression level and length of genes, position and context of codons in the genes, recombination rates, mRNA folding, and tRNA abundance and interactions are some factors influencing codon bias. The factors shaping codon bias may also be involved in evolution of the universal genetic code. Codon-usage bias is critical factor determining gene expression and cellular function by influencing diverse processes such as RNA processing, protein translation and protein folding. Codon usage bias reflects the origin, mutation patterns and evolution of the species or genes. Investigations of codon bias patterns in genomes can reveal phylogenetic relationships between organisms, horizontal gene transfers, molecular evolution of genes and identify selective forces that drive their evolution. Most important application of codon bias analysis is in the design of transgenes, to increase gene expression levels through codon optimization, for development of transgenic crops. The review gives an overview of deviations of genetic code, factors influencing codon usage or bias, codon usage bias of nuclear and organellar genes, computational methods to determine codon usage and the significance as well as applications of codon usage analysis in biological research, with emphasis on plants.

RNA2-encoded VP37 protein of Broad bean wilt virus 1 is a determinant of pathogenicity, host susceptibility, and a suppressor of post-transcriptional gene silencing

Broad bean wilt virus 1 (BBWV-1, genus Fabavirus, family Secoviridae) is a bipartite, single-stranded positive-sense RNA virus infecting many horticultural and ornamental crops worldwide. RNA1 encodes proteins involved in viral replication whereas RNA2 encodes two coat proteins (the large and small coat proteins) and two putative movement proteins (MPs) of different sizes with overlapping C-terminal regions. In this work, we determined the role played by the small putative BBWV-1 MP (VP37) on virus pathogenicity, host specificity, and suppression of post-transcriptional gene silencing (PTGS). We engineered a BBWV-1 35S-driven full-length cDNA infectious clone corresponding to BBWV-1 RNA1 and RNA2 (pBBWV1-Wt) and generated a mutant knocking out VP37 (pBBWV1-G492C). Agroinfiltration assays showed that pBBWV1-Wt, as the original BBWV-1 isolate, infected broad bean, tomato, pepper, and Nicotiana benthamiana, whereas pBBWV1-G492C did not infect pepper and tomato systemically. Also, pBBWV1-G492C induced milder symptoms in broad bean and N. benthamiana than pBBWV1-Wt. Differential retrotranscription and amplification of the (+) and (-) strands showed that pBBWV1-G492C replicated in the agroinfiltrated leaves of pepper but not in tomato. All this suggests that VP37 is a determinant of pathogenicity and host specificity. Transient expression of VP37 through a potato virus X (PVX) vector enhanced PVX symptoms and induced systemic necrosis associated with programmed cell death in N. benthamiana plants. Finally, VP37 was identified as a viral suppressor of RNA silencing by transient expression in N. benthamiana 16c plants and movement complementation of a viral construct based on turnip crinkle virus (pTCV-GFP).

Novel pipeline identifies new upstream ORFs and non-AUG initiating main ORFs with conserved amino acid sequences in the 5' leader of mRNAs in Arabidopsis thaliana

Eukaryotic mRNAs contain a 5' leader sequence preceding the main open reading frame (mORF) and, depending on the species, 20%-50% of eukaryotic mRNAs harbor an upstream ORF (uORF) in the 5' leader. An unknown fraction of these uORFs encode sequence conserved peptides (conserved peptide uORFs, CPuORFs). Experimentally validated CPuORFs demonstrated to regulate the translation of downstream mORFs often do so in a metabolite concentration-dependent manner. Previous research has shown that most CPuORFs possess a start codon context suboptimal for translation initiation, which turns out to be favorable for translational regulation. The suboptimal initiation context may even include non-AUG start codons, which makes CPuORFs hard to predict. For this reason, we developed a novel pipeline to identify CPuORFs unbiased of start codon using well-annotated sequence data from 31 eudicot plant species and rice. Our new pipeline was able to identify 29 novel Arabidopsis thaliana (Arabidopsis) CPuORFs, conserved across a wide variety of eudicot species of which 15 do not initiate with an AUG start codon. In addition to CPuORFs, the pipeline was able to find 14 conserved coding regions directly upstream and in frame with the mORF, which likely initiate translation on a non-AUG start codon. Altogether, our pipeline identified highly conserved coding regions in the 5' leaders of Arabidopsis transcripts, including in genes with proven functional importance such as LHY, a key regulator of the circadian clock, and the RAPTOR1 subunit of the target of rapamycin (TOR) kinase.

Isolation and characterization of a wheat IF2 homolog required for innate immunity to stripe rust

The wheat eIF2 homolog, TaIF2, is induced by the stripe rust pathogen CYR 32 at an early stage of inoculation and is related to the innate immunity resistance level in wheat. The initiation of translation represents a critical control point in the regulation of gene expression in all organisms. We previously identified an upregulated EST S186 (EL773056) from an SSH-cDNA library of the Shaanmai 139 strain of wheat (Triticum aestivum) infected with Puccinia striiformis (Pst). In the present work, we isolated a cDNA clone and identified it as a wheat IF2 homolog. This cDNA consisted of 1,314 nucleotides and contained an open reading frame of 795 nucleotides encoding a polypeptide of 254 amino acids. The amino acids represent a conserved domain in EF-Tu, mtIF2-II, and mtIF2-Ivc. The alignment result showed that it maybe a partial cDNA of the initiation factor 2/eukaryotic initiation factor 5B (IF2/eIF5B) superfamily gene. Paradoxically, results of a Swiss-model analysis suggesting a low QMEAN Z-score implied that it was a membrane protein. Quantitative RT-PCR studies confirmed that the wheat eIF2 (TaIF2) homolog was differentially expressed in three near-isogenic lines. Critical time points for the induction of resistance by inoculation with Pst CYR32 in YrSM139-1B + YrSM139-2D immune resistance genotype occurred at 1 and 3 dpi (days post-infection). RNAi test showed that the inoculated BSMV-IF2 leaves of Shaanmai 139 showed obvious cell death after 15 days of inoculation with CYR 32. qRT-PCR analysis of the target gene in cDNA samples isolated from BSMV-IF2-Pst, BSMV-0-Pst and Pst infected leaves confirmed that the expression of TaIF2 is suppressed by BSMV-IF2 at 3 dpi. This suggested that TaIF2/eIF5B plays an important role in the mechanism of innate immunity to stripe rust pathogen.

Binding of more than one Tva800 molecule is required for ASLV-A entry

BACKGROUND

Understanding the mechanism by which viruses enter their target cell is an essential part of understanding their infectious cycle. Previous studies have focussed on the multiplicity of viral envelope proteins that need to bind to their cognate receptor to initiate entry. Avian sarcoma and leukosis virus Envelope protein (ASLV Env) mediates entry via a receptor, Tva, which can be attached to the cell surface either by a phospholipid anchor (Tva800) or a transmembrane domain (Tva950). In these studies, we have now investigated the number of target receptors necessary for entry of ASLV Env-pseudotyped virions.

RESULTS

Using titration and modelling experiments we provide evidence that binding of more than one receptor, probably two, is needed for entry of virions via Tva800. However, binding of just one Tva950 receptor is sufficient for successful entry.

CONCLUSIONS

The different modes of attachment of Tva800 and Tva950 to the cell membrane have important implications for the utilisation of these proteins as receptors for viral binding and/or uptake.

Introducing an RNA editing requirement into a plastid-localised transgene reduces but does not eliminate functional gene transfer to the nucleus

In higher plants, DNA transfer from the plastid (chloroplast) genome to the nucleus is a frequent, ongoing process. However, there has been uncertainty over whether this transfer occurs by a direct DNA mechanism or whether RNA intermediates are involved. Previous experiments utilising transplastomic Nicotiana tabacum (tp7 and tp17) enabled the detection of plastid-to-nucleus transfer in real time. To determine whether RNA intermediates are involved in this transfer, transplastomic lines (tpneoACG) were generated containing, in their plastid genomes, a nuclear promoter-driven kanamycin resistance gene (neo) with a start codon that required plastid RNA editing but otherwise identical to tp7 and tp17. Therefore it was expected that kanamycin resistance would only be acquired following RNA-mediated transfer of neo to the nucleus. Screening of tpneoACG progeny revealed several kanamycin-resistant plants, each of which contained the neo gene located in the nucleus. Surprisingly, neo was unedited in all these plants, indicating that neoACG was active in the absence of an edited start codon and suggesting that RNA intermediates were not involved in the transfers. However, analysis of tpneoACG revealed that only a low proportion of transcripts potentially able to mediate neo transfer were edited, thus precluding unequivocal conclusions regarding the role of RNA in plastid-to-nucleus transfer. The low proportion of edited transcripts was found to be due to predominant antisense neo transcripts, rather than to low editing efficiency of the sense transcripts. This study highlights a number of important considerations in the design of experiments utilising plastid RNA editing. The results also suggest that RNA editing sites reduce but do not eliminate functional plastid-to-nucleus gene transfer. This is relevant both in an evolutionary context and in placing RNA editing-dependent genes in the plastid genome as a means of transgene containment.

Evolution of plant phage-type RNA polymerases: the genome of the basal angiosperm Nuphar advena encodes two mitochondrial and one plastid phage-type RNA polymerases

Background

In mono- and eudicotyledonous plants, a small nuclear gene family (RpoT, RNA polymerase of the T3/T7 type) encodes mitochondrial as well as chloroplast RNA polymerases homologous to the T-odd bacteriophage enzymes. RpoT genes from angiosperms are well characterized, whereas data from deeper branching plant species are limited to the moss Physcomitrella and the spikemoss Selaginella. To further elucidate the molecular evolution of the RpoT polymerases in the plant kingdom and to get more insight into the potential importance of having more than one phage-type RNA polymerase (RNAP) available, we searched for the respective genes in the basal angiosperm Nuphar advena.

Results

By screening a set of BAC library filters, three RpoT genes were identified. Both genomic gene sequences and full-length cDNAs were determined. The NaRpoT mRNAs specify putative polypeptides of 996, 990 and 985 amino acids, respectively. All three genes comprise 19 exons and 18 introns, conserved in their positions with those known from RpoT genes of other land plants. The encoded proteins show a high degree of conservation at the amino acid sequence level, including all functional crucial regions and residues known from the phage T7 RNAP. The N-terminal transit peptides of two of the encoded polymerases, NaRpoTm1 and NaRpoTm2, conferred targeting of green fluorescent protein (GFP) exclusively to mitochondria, whereas the third polymerase, NaRpoTp, was targeted to chloroplasts. Remarkably, translation of NaRpoTp mRNA has to be initiated at a CUG codon to generate a functional plastid transit peptide. Thus, besides AGAMOUS in Arabidopsis and the Nicotiana RpoTp gene, N. advena RpoTp provides another example for a plant mRNA that is exclusively translated from a non-AUG codon. In contrast to the RpoT of the lycophyte Selaginella and those of the moss Physcomitrella, which are according to phylogenetic analyses in sister positions to all other phage-type polymerases of angiosperms, the Nuphar RpoTs clustered with the well separated clades of mitochondrial (NaRpoTm1 and NaRpoTm2) and plastid (NaRpoTp) polymerases.

Conclusions

Nuphar advena encodes two mitochondrial and one plastid phage-type RNAP. Identification of a plastid-localized phage-type RNAP in this basal angiosperm, orthologous to all other RpoTp enzymes of flowering plants, suggests that the duplication event giving rise to a nuclear gene-encoded plastid RNA polymerase, not present in lycopods, took place after the split of lycopods from all other tracheophytes. A dual-targeted mitochondrial and plastididal RNA polymerase (RpoTmp), as present in eudicots but not monocots, was not detected in Nuphar suggesting that its occurrence is an evolutionary novelty of eudicotyledonous plants like Arabidopsis.

Noncanonical translation initiation of the Arabidopsis flowering time and alternative polyadenylation regulator FCA

The RNA binding protein FCA regulates the floral transition and is required for silencing RNAs corresponding to specific noncoding sequences in the Arabidopsis thaliana genome. Through interaction with the canonical RNA 3' processing machinery, FCA affects alternative polyadenylation of many transcripts, including antisense RNAs at the locus encoding the floral repressor FLC. This potential for widespread alteration of gene regulation clearly needs to be tightly regulated, and we have previously shown that FCA expression is autoregulated through poly(A) site choice. Here, we show distinct layers of FCA regulation that involve sequences within the 5' region that regulate noncanonical translation initiation and alter the expression profile. FCA translation in vivo occurs exclusively at a noncanonical CUG codon upstream of the first in-frame AUG. We fully define the upstream flanking sequences essential for its selection, revealing features that distinguish this from other non-AUG start site mechanisms. Bioinformatic analysis identified 10 additional Arabidopsis genes that likely initiate translation at a CUG codon. Our findings reveal further unexpected complexity in the regulation of FCA expression with implications for its roles in regulating flowering time and gene expression and more generally show plant mRNA exceptions to AUG translation initiation.

Characterisation of a new reporter system allowing high throughput in planta screening for recombination events before and after controlled DNA double strand break induction

DNA double strand breaks (DSBs) are created either by DNA damaging reagents or in a programmed manner, for example during meiosis. Homologous recombination (HR) can be used to repair DSBs, a process vital both for cell survival and for genetic rearrangement during meiosis. In order to easily quantify this mechanism, a new HR reporter gene that is suitable for the detection of rare recombination events in high-throughput screens was developed in Arabidopsis thaliana. This reporter, pPNP, is composed of two mutated Pat genes and has also one restriction site for the meganuclease I-SceI. A functional Pat gene can be reconstituted by an HR event giving plants which are resistant to the herbicide glufosinate. The basal frequency of intra-chromosomal recombination is very low (10(-5)) and can be strongly increased by the expression of I-SceI which creates a DSB. Expression of I-SceI under the control of the 35S CaMV promoter dramatically increases HR frequency (10,000 fold); however the measured recombinant events are in majority somatic. In contrast only germinal recombination events were measured when the meganuclease was expressed from a floral-specific promoter. Finally, the reporter was used to test a dexamethasone inducible I-SceI which could produce up to 200x more HR events after induction. This novel inducible I-SceI should be useful in fundamental studies of the mechanism of repair of DSBs and for biotechnological applications.

Impact of the loss of AtMSH2 on double-strand break-induced recombination between highly diverged homeologous sequences in Arabidopsis thaliana germinal tissues

We experimented a novel reporter system to analyze intrachromosomal recombination between homeologous sequences in Arabidopsis germ cell lineages. The recombination substrates used are the BAR and PAT genes which diverge by about 13% at the nucleotide level and confer resistance to the herbicide glufosinate. DNA double-strand breaks (DSBs) were generated by the I-Sce1 endonuclease to induce recombination. Loss of AtMSH2 induces a 3-fold increase of the frequency of recombination events indicating that AtMSH2 is involved in the anti-recombination activity that prevents exchange between highly diverged sequences in Arabidopsis. Molecular analysis of recombined alleles indicates that in wild type plants the single strand annealing (SSA) pathway can process more efficiently homologous 3' ends than 3' ends generated by resection of non-homologous overhangs. The loss of AtMSH2 disturbs this process, leading to a modification of the distribution of the BAR/PAT junctions and therefore showing that the MSH2 function is also involved in determining the structure of the recombined alleles. In addition, conversion tracts were observed in some alleles. They are shorter in MSH2 deficient plants than in wild-type, suggesting that a short-patch mismatch repair, not controlled by MSH2, could exist in Arabidopsis.

A survey of mRNA sequences with a non-AUG start codon in RefSeq database

Alternative initiation in translation is one of the important mechanisms in which multiple proteins are synthesized from a single mRNA. In many cases, translation initiation occurring at a non-AUG codon has been reported by several experimental studies. We have analyzed all mRNA sequences in the RefSeq database and found that coding regions of about 0.1% of the total mRNA sequences begin with a non-AUG codon (nonAUG mRNAs). Major fraction of non-AUG mRNAs is predicted from genomic sequences. More than 100 non-AUG sequences are highly curated and 52 of them are explicitly annotated that they use alternate start codons for translation initiation. Analysis of these sequences reveals that majority of the protein products contain domains that are DNA/RNA-binding, kinases, growth factors, or involved in immune response or cell proliferation. Thus, the proteins translated from non-canonical codons seem to be implicated in regulatory role and/or signaling mechanism. The sequence context of the non-AUG start codons shows that purine at -3 position and/or G at +4 position are strongly conserved and the corresponding genes give rise to alternate transcripts and/or multiple isoforms. We have also developed a database "nonAUG" (http://bioinfo.iitk.ac.in) that contains a collection of all mRNA sequences whose coding regions start with a non-AUG codon. nonAUG database will be continuously updated and is freely available to the scientific community.