The impact of metabolic engineering on nontarget pathways and outcomes of metabolic engineering from different genomes are poorly understood questions. Therefore, squalene biosynthesis genes FARNESYL DIPHOSPHATE SYNTHASE (FPS) and SQUALENE SYNTHASE (SQS) were engineered via the Nicotiana tabacum chloroplast (C), nuclear (N) or both (CN) genomes to promote squalene biosynthesis. SQS levels were ~4300-fold higher in C and CN lines than in N, but all accumulated ~150-fold higher squalene due to substrate or storage limitations. Abnormal leaf and flower phenotypes, including lower pollen production and reduced fertility, were observed regardless of the compartment or level of transgene expression. Substantial changes in metabolomes of all lines were observed: levels of 65-120 unrelated metabolites, including the toxic alkaloid nicotine, changed by as much as 32-fold. Profound effects of transgenesis on nontarget gene expression included changes in the abundance of 19 076 transcripts by up to 2000-fold in CN; 7784 transcripts by up to 1400-fold in N; and 5224 transcripts by as much as 2200-fold in C. Transporter-related transcripts were induced, and cell cycle-associated transcripts were disproportionally repressed in all three lines. Transcriptome changes were validated by qRT-PCR. The mechanism underlying these large changes likely involves metabolite-mediated anterograde and/or retrograde signalling irrespective of the level of transgene expression or end product, due to imbalance of metabolic pools, offering new insight into both anticipated and unanticipated consequences of metabolic engineering.

The functions of the major mammalian cytoplasmic poly(A) binding protein, PABPC1, have been characterized predominantly in the context of its binding to the 3' poly(A) tails of mRNAs. These interactions play important roles in post-transcriptional gene regulation by enhancing translation and mRNA stability. Here, we performed transcriptome-wide CLIP-seq analysis to identify additional PABPC1 binding sites within genomically encoded mRNA sequences that may impact on gene regulation. From this analysis, we found that PABPC1 binds directly to the canonical polyadenylation signal in thousands of mRNAs in the mouse transcriptome. PABPC1 binding also maps to translation initiation and termination sites bracketing open reading frames, exemplified most dramatically in replication-dependent histone mRNAs. Additionally, a more restricted subset of PABPC1 interaction sites comprised A-rich sequences within the 5' UTRs of mRNAs, including Pabpc1 mRNA itself. Functional analyses revealed that these PABPC1 interactions in the 5' UTR mediate both auto- and trans-regulatory translational control. In total, these findings reveal a repertoire of PABPC1 binding that is substantially broader than previously recognized with a corresponding potential to impact and coordinate post-transcriptional controls critical to a broad array of cellular functions.

Posttranscriptional regulation in eukaryotes requires cis- and trans-acting features and factors including RNA secondary structure and RNA-binding proteins (RBPs). However, a comprehensive view of the structural and RBP interaction landscape of nuclear RNAs has yet to be compiled for any organism. Here, we use our ribonuclease-mediated structure and RBP-binding site mapping approaches to globally profile these features in Arabidopsis seedling nuclei in vivo. We reveal anticorrelated patterns of secondary structure and RBP binding throughout nuclear mRNAs that demarcate sites of alternative splicing and polyadenylation. We also uncover a collection of protein-bound sequence motifs, and identify their structural contexts, co-occurrences in transcripts encoding functionally related proteins, and interactions with putative RBPs. Finally, using these motifs, we find that the chloroplast RBP CP29A also interacts with nuclear mRNAs. In total, we provide a simultaneous view of the RNA secondary structure and RBP interaction landscapes in a eukaryotic nucleus.

The MSI2 RNA-binding protein is a potent oncogene playing key roles in haematopoietic stem cell homeostasis and malignant haematopoiesis. Here we demonstrate that MSI2 is expressed in the intestinal stem cell compartment, that its expression is elevated in colorectal adenocarcinomas, and that MSI2 loss-of-function abrogates colorectal cancer cell growth. MSI2 gain-of-function in the intestinal epithelium in a drug-inducible mouse model is sufficient to phenocopy many of the morphological and molecular consequences of acute loss of the APC tumour suppressor in the intestinal epithelium in a Wnt-independent manner. Transcriptome-wide RNA-binding analysis indicates that MSI2 acts as a pleiotropic inhibitor of known intestinal tumour suppressors including Lrig1, Bmpr1a, Cdkn1a and Pten. Finally, we demonstrate that inhibition of the PDK-AKT-mTORC1 axis rescues oncogenic consequences of MSI2 induction. Taken together, our findings identify MSI2 as a central component in an unappreciated oncogenic pathway promoting intestinal transformation.

As the relevant literature and the number of experiments increase at a super linear rate, databases that curate and collect experimentally verified microRNA (miRNA) targets have gradually emerged. These databases attempt to provide efficient access to this wealth of experimental data, which is scattered in thousands of manuscripts. Aim of TarBase 6.0 (http://www.microrna.gr/tarbase) is to face this challenge by providing a significant increase of available miRNA targets derived from all contemporary experimental techniques (gene specific and high-throughput), while incorporating a powerful set of tools in a user-friendly interface. TarBase 6.0 hosts detailed information for each miRNA-gene interaction, ranging from miRNA- and gene-related facts to information specific to their interaction, the experimental validation methodologies and their outcomes. All database entries are enriched with function-related data, as well as general information derived from external databases such as UniProt, Ensembl and RefSeq. DIANA microT miRNA target prediction scores and the relevant prediction details are available for each interaction. TarBase 6.0 hosts the largest collection of manually curated experimentally validated miRNA-gene interactions (more than 65,000 targets), presenting a 16.5-175-fold increase over other available manually curated databases.

microRNAs (miRNAs) are small endogenous RNA molecules that are implicated in many biological processes through post-transcriptional regulation of gene expression. The DIANA-microT Web server provides a user-friendly interface for comprehensive computational analysis of miRNA targets in human and mouse. The server has now been extended to support predictions for two widely studied species: Drosophila melanogaster and Caenorhabditis elegans. In the updated version, the Web server enables the association of miRNAs to diseases through bibliographic analysis and provides insights for the potential involvement of miRNAs in biological processes. The nomenclature used to describe mature miRNAs along different miRBase versions has been extensively analyzed, and the naming history of each miRNA has been extracted. This enables the identification of miRNA publications regardless of possible nomenclature changes. User interaction has been further refined allowing users to save results that they wish to analyze further. A connection to the UCSC genome browser is now provided, enabling users to easily preview predicted binding sites in comparison to a wide array of genomic tracks, such as single nucleotide polymorphisms. The Web server is publicly accessible in www.microrna.gr/microT-v4.

RNA silencing pathways play critical roles in gene regulation, virus infection, and transposon control. RNA interference (RNAi) is mediated by small interfering RNAs (siRNAs), which are liberated from double-stranded (ds)RNA precursors by Dicer and guide the RNA-induced silencing complex (RISC) to targets. Although principles governing small RNA sorting into RISC have been uncovered, the spectrum of RNA species that can be targeted by Dicer proteins, particularly the viral RNAs present during an infection, are poorly understood. Dicer-2 potently restricts viral infection in insects by generating virus-derived siRNAs from viral RNA. To better characterize the substrates of Dicer-2, we examined the virus-derived siRNAs produced during the Drosophila antiviral RNAi response to four different viruses using high-throughput sequencing. We found that each virus was uniquely targeted by the RNAi pathway; dicing substrates included dsRNA replication intermediates and intramolecular RNA stem loops. For instance, a putative intergenic RNA hairpin encoded by Rift Valley Fever virus generates abundant small RNAs in both Drosophila and mosquito cells, while repetitive sequences within the genomic termini of Vaccinia virus, which give rise to abundant small RNAs in Drosophila, were found to be transcribed in both insect and mammalian cells. Moreover, we provide evidence that the RNA species targeted by Dicer-2 can be modulated by the presence of a viral suppressor of RNAi. This study uncovered several novel, heavily targeted features within viral genomes, offering insight into viral replication, viral immune evasion strategies, and the mechanism of antiviral RNAi.

The functional structure of all biologically active molecules is dependent on intra- and inter-molecular interactions. This is especially evident for RNA molecules whose functionality, maturation, and regulation require formation of correct secondary structure through encoded base-pairing interactions. Unfortunately, intra- and inter-molecular base-pairing information is lacking for most RNAs. Here, we marry classical nuclease-based structure mapping techniques with high-throughput sequencing technology to interrogate all base-paired RNA in Arabidopsis thaliana and identify ∼200 new small (sm)RNA-producing substrates of RNA-DEPENDENT RNA POLYMERASE6. Our comprehensive analysis of paired RNAs reveals conserved functionality within introns and both 5' and 3' untranslated regions (UTRs) of mRNAs, as well as a novel population of functional RNAs, many of which are the precursors of smRNAs. Finally, we identify intra-molecular base-pairing interactions to produce a genome-wide collection of RNA secondary structure models. Although our methodology reveals the pairing status of RNA molecules in the absence of cellular proteins, previous studies have demonstrated that structural information obtained for RNAs in solution accurately reflects their structure in ribonucleoprotein complexes. Furthermore, our identification of RNA-DEPENDENT RNA POLYMERASE6 substrates and conserved functional RNA domains within introns and both 5' and 3' untranslated regions (UTRs) of mRNAs using this approach strongly suggests that RNA molecules are correctly folded into their secondary structure in solution. Overall, our findings highlight the importance of base-paired RNAs in eukaryotes and present an approach that should be widely applicable for the analysis of this key structural feature of RNA.

PIWI family proteins bind to small RNAs known as PIWI-interacting RNAs (piRNAs) and play essential roles in the germline by silencing transposons and by promoting germ cell specification and function. Here we report that the widely used Kc167 cell line, derived from Drosophila melanogaster embryos, expresses piRNAs that are loaded to Aub and Piwi. Kc167 piRNAs are produced by a canonical, primary piRNA biogenesis pathway, from phased processing of precursor transcripts by the Zuc endonuclease, Armi helicase, and dGasz mitochondrial scaffold protein. Kc167 piRNAs derive from cytoplasmic transcripts, notably tRNAs and mRNAs, and their abundance correlates with that of parent transcripts. The expression of Aub is robust in Kc167, that of Piwi is modest, while Ago3 is undetectable, explaining the lack of transposon-related piRNA amplification by the Aub-Ago3, ping-pong mechanism. We propose that the default state of the primary piRNA biogenesis machinery is random transcript sampling to allow generation of piRNAs from any transcript, including newly acquired retrotransposons. This state is unmasked in Kc167, likely because they do not express piRNA cluster transcripts in sufficient amounts and do not amplify transposon piRNAs. We use Kc167 to characterize an inactive isoform of Aub protein. Since most Kc167 piRNAs are genic, they can be mapped uniquely to the genome, facilitating computational analyses. Furthermore, because Kc167 is a widely used and well-characterized cell line that is easily amenable to experimental manipulations, we expect that it will serve as an excellent system to study piRNA biogenesis and piRNA-related factors.

After transcription, a messenger RNA (mRNA) is further post-transcriptionally regulated by several features including RNA secondary structure and covalent RNA modifications (specifically N6-methyladenosine, m6A). Both RNA secondary structure and m6A have been demonstrated to regulate mRNA stability and translation and have been independently linked to plant responses to soil salinity levels. However, the effect of m6A on regulating RNA secondary structure and the combinatorial interplay between these two RNA features during salt stress response has yet to be studied. Here, we globally identify RNA-protein interactions and RNA secondary structure during systemic salt stress. This analysis reveals that RNA secondary structure changes significantly during salt stress, and that it is independent of global changes in RNA-protein interactions. Conversely, we find that m6A is anti-correlated with RNA secondary structure in a condition-dependent manner, with salt-specific m6A correlated with a decrease in mRNA secondary structure during salt stress. Taken together, we suggest that salt-specific m6A deposition and the associated loss of RNA secondary structure results in increases in mRNA stability for transcripts encoding abiotic stress response proteins and ultimately increases in protein levels from these stabilized transcripts. In total, our comprehensive analyses reveal important post-transcriptional regulatory mechanisms involved in plant long-term salt stress response and adaptation.

Recent high-throughput structure-sensitive genome-wide sequencing-based assays have enabled large-scale studies of RNA structure, and robust transcriptome-wide computational prediction of individual RNA structures across RNA classes from these assays has potential to further improve the prediction accuracy. Here, we describe HiPR, a novel method for RNA structure prediction at single-nucleotide resolution that combines high-throughput structure probing data (DMS-seq, DMS-MaPseq) with a novel probabilistic folding algorithm. On validation data spanning a variety of RNA classes, HiPR often increases accuracy for predicting RNA structures, giving researchers new tools to study RNA structure.

Cerebral malaria (CM), coma caused by Plasmodium falciparum-infected red blood cells (iRBCs), is the deadliest complication of malaria. The mechanisms that lead to CM development are incompletely understood. Here we report on the identification of activation and inhibition pathways leading to mouse CM with supporting evidence from the analysis of human specimens. We find that CM suppression can be induced by vascular injury when sporozoites exit the circulation to infect the liver and that CM suppression is mediated by the release of soluble factors into the circulation. Among these factors is insulin like growth factor 1 (IGF1), administration of which inhibits CM development in mice. IMPORTANCE Liver infection by Plasmodium sporozoites is a required step for infection of the organism. We found that alternate pathways of sporozoite liver infection differentially influence cerebral malaria (CM) development. CM is one of the primary causes of death following malaria infection. To date, CM research has focused on how CM phenotypes develop but no successful therapeutic treatment or prognostic biomarkers are available. Here we show for the first time that sporozoite liver invasion can trigger CM-inhibitory immune responses. Importantly, we identified a number of early-stage prognostic CM inhibitory biomarkers, many of which had never been associated with CM development. Serological markers identified using a mouse model are directly relevant to human CM.

Polyadenylation plays a key role in producing mature mRNAs in eukaryotes. It is widely believed that the poly(A)-binding proteins (PABs) uniformly bind to poly(A)-tailed mRNAs, regulating their stability and translational efficiency.

Long intergenic noncoding RNAs (lincRNAs) are a large yet enigmatic class of eukaryotic transcripts that can have critical biological functions. The wealth of RNA-sequencing (RNA-seq) data available for plants provides the opportunity to implement a harmonized identification and annotation effort for lincRNAs that enables cross-species functional and genomic comparisons as well as prioritization of functional candidates. In this study, we processed >24 Tera base pairs of RNA-seq data from >16,000 experiments to identify ∼130,000 lincRNAs in four Brassicaceae: Arabidopsis thaliana, Camelina sativa, Brassica rapa, and Eutrema salsugineum. We used nanopore RNA-seq, transcriptome-wide structural information, peptide data, and epigenomic data to characterize these lincRNAs and identify conserved motifs. We then used comparative genomic and transcriptomic approaches to highlight lincRNAs in our data set with sequence or transcriptional conservation. Finally, we used guilt-by-association analyses to assign putative functions to lincRNAs within our data set. We tested this approach on a subset of lincRNAs associated with germination and seed development, observing germination defects for Arabidopsis lines harboring T-DNA insertions at these loci. LincRNAs with Brassicaceae-conserved putative miRNA binding motifs, small open reading frames, or abiotic-stress modulated expression are a few of the annotations that will guide functional analyses into this cryptic portion of the transcriptome.

The evolutionarily conserved DNA repair complex Ku serves as the primary sensor of free DNA ends in eukaryotic cells. Its rapid association with DNA ends is crucial for several cellular processes, including non-homologous end joining (NHEJ) DNA repair and telomere protection. In this study, we conducted a transient kinetic analysis to investigate the impact of the SAP domain on individual phases of the Ku-DNA interaction. Specifically, we examined the initial binding, the subsequent docking of Ku onto DNA, and sliding of Ku along DNA. Our findings revealed that the C-terminal SAP domain of Ku70 facilitates the initial phases of the Ku-DNA interaction but does not affect the sliding process. This suggests that the SAP domain may either establish the first interactions with DNA, or stabilize these initial interactions during loading. To assess the biological role of the SAP domain, we generated Arabidopsis plants expressing Ku lacking the SAP domain. Intriguingly, despite the decreased efficiency of the ΔSAP Ku complex in loading onto DNA, the mutant plants exhibited full proficiency in classical NHEJ and telomere maintenance. This indicates that the speed with which Ku loads onto telomeres or DNA double-strand breaks is not the decisive factor in stabilizing these DNA structures.

Expansins facilitate cell expansion by mediating pH-dependent cell wall (CW) loosening. However, the role of expansins in controlling CW biomechanical properties in specific tissues and organs remains elusive. We monitored hormonal responsiveness and spatial specificity of expression and localization of expansins predicted to be the direct targets of cytokinin signaling in Arabidopsis (Arabidopsis thaliana). We found EXPANSIN1 (EXPA1) homogenously distributed throughout the CW of columella/lateral root cap, while EXPA10 and EXPA14 localized predominantly at 3-cell boundaries in the epidermis/cortex in various root zones. EXPA15 revealed cell-type-specific combination of homogenous vs. 3-cell boundaries localization. By comparing Brillouin frequency shift and AFM-measured Young's modulus, we demonstrated Brillouin light scattering (BLS) as a tool suitable for non-invasive in vivo quantitative assessment of CW viscoelasticity. Using both BLS and AFM, we showed that EXPA1 overexpression upregulated CW stiffness in the root transition zone (TZ). The dexamethasone-controlled EXPA1 overexpression induced fast changes in the transcription of numerous CW-associated genes, including several EXPAs and XYLOGLUCAN:XYLOGLUCOSYL TRANSFERASEs (XTHs), and associated with rapid pectin methylesterification determined by in situ Fourier-transform infrared spectroscopy in the root TZ. The EXPA1-induced CW remodeling is associated with the shortening of the root apical meristem, leading to root growth arrest. Based on our results, we propose that expansins control root growth by a delicate orchestration of CW biomechanical properties, possibly regulating both CW loosening and CW remodeling.

Recently, post-transcriptional regulation of mRNA mediated by N6-methyladenosine (m6A) has been found to have profound effects on transcriptome regulation during plant responses to various abiotic stresses. However, whether this RNA modification can affect an oxidative stress response in plants has not been studied. To assess the role of m6A modifications during copper-induced oxidative stress responses, m6A-IP-seq was performed in Arabidopsis seedlings exposed to high levels of copper sulfate. This analysis revealed large-scale shifts in this modification on the transcripts most relevant for oxidative stress. This altered epitranscriptomic mark is known to influence transcript abundance and translation; therefore we scrutinized these possibilities. We found an increased abundance of copper-enriched m6A-containing transcripts. Similarly, we also found increased ribosome occupancy of copper-enriched m6A-containing transcripts, specifically those encoding proteins involved with stress responses relevant to oxidative stressors. Furthermore, the significance of the m6A epitranscriptome on plant oxidative stress tolerance was uncovered by assessing germination and seedling development of the mta (N6-methyladenosine RNA methyltransferase A mutant complemented with ABI3:MTA) mutant exposed to high copper treatment. These analyses suggested hypersensitivity of the mta mutant compared to the wild-type plants in response to copper-induced oxidative stress. Overall, our findings suggest an important role for m6A in the oxidative stress response of Arabidopsis.

Germ cells implement elaborate mechanisms to protect their genetic material and to regulate gene expression during differentiation. Piwi proteins bind Piwi-interacting RNAs (piRNAs), small germline RNAs whose biogenesis and functions are still largely elusive. We used high-throughput sequencing after cross-linking and immunoprecipitation (HITS-CLIP) coupled with RNA-sequencing (RNA-seq) to characterize the genome-wide target RNA repertoire of Mili (Piwil2) and Miwi (Piwil1), two Piwi proteins expressed in mouse postnatal testis. We report the in vivo pathway of primary piRNA biogenesis and implicate distinct nucleolytic activities that process Piwi-bound precursor transcripts. Our studies indicate that pachytene piRNAs are the end products of RNA processing. HITS-CLIP demonstrated that Miwi binds spermiogenic mRNAs directly, without using piRNAs as guides, and independent biochemical analyses of testis mRNA ribonucleoproteins (mRNPs) established that Miwi functions in the formation of mRNP complexes that stabilize mRNAs essential for spermiogenesis.

Epigenetic mechanisms, including histone post-translational modifications, control longevity in diverse organisms. Relatedly, loss of proper transcriptional regulation on a global scale is an emerging phenomenon of shortened life span, but the specific mechanisms linking these observations remain to be uncovered. Here, we describe a life span screen in Saccharomyces cerevisiae that is designed to identify amino acid residues of histones that regulate yeast replicative aging. Our results reveal that lack of sustained histone H3K36 methylation is commensurate with increased cryptic transcription in a subset of genes in old cells and with shorter life span. In contrast, deletion of the K36me2/3 demethylase Rph1 increases H3K36me3 within these genes, suppresses cryptic transcript initiation, and extends life span. We show that this aging phenomenon is conserved, as cryptic transcription also increases in old worms. We propose that epigenetic misregulation in aging cells leads to loss of transcriptional precision that is detrimental to life span, and, importantly, this acceleration in aging can be reversed by restoring transcriptional fidelity.

MicroRNAs (miRNAs) and small interfering RNAs (siRNAs) are abundant endogenous small RNAs (smRNAs) that control transcript expression through posttranscriptional gene silencing. Here, we show that concomitant loss of XRN4/EIN5, a 5'-3' exoribonuclease, and ABH1/CBP80, a subunit of the mRNA cap binding complex, results in Arabidopsis plants manifesting myriad developmental defects. We find that ABH1/CBP80 is necessary to obtain proper mature miRNA levels, which suggests this protein affects the miRNA-mediated RNA silencing pathway. Additionally, we show that XRN4/EIN5 affects the levels of a smRNA class that is processed from both sense and antisense strands of approximately 130 endogenous transcripts that apparently are converted to double-stranded RNA (dsRNA) and subsequently processed. We find that the parent transcripts of these smRNAs accumulate in an uncapped form upon loss of XRN4/EIN5, which suggests that uncapped endogenous transcripts can become smRNA biogenesis substrates. Overall, our results reveal unexpected connections between RNA metabolism and silencing pathways.