ABSTRACT
The 3' end of a gene, often called a terminator, modulates mRNA stability, localization, translation, and polyadenylation. Here, we adapted Plant STARR-seq, a massively parallel reporter assay, to measure the activity of over 50,000 terminators from the plants Arabidopsis thaliana and Zea mays. We characterize thousands of plant terminators, including many that outperform bacterial terminators commonly used in plants. Terminator activity is species-specific, differing in tobacco leaf and maize protoplast assays. While recapitulating known biology, our results reveal the relative contributions of polyadenylation motifs to terminator strength. We built a computational model to predict terminator strength and used it to conduct in silico evolution that generated optimized synthetic terminators. Additionally, we discover alternative polyadenylation sites across tens of thousands of terminators; however, the strongest terminators tend to have a dominant cleavage site. Our results establish features of plant terminator function and identify strong naturally occurring and synthetic terminators.
Subject(s)
Arabidopsis , Polyadenylation , Zea mays , Zea mays/genetics , Zea mays/metabolism , Arabidopsis/genetics , Arabidopsis/metabolism , Gene Expression Regulation, Plant , Terminator Regions, Genetic/genetics , Nicotiana/genetics , Nicotiana/metabolism , RNA, Messenger/genetics , RNA, Messenger/metabolismABSTRACT
The transcription termination process is an important part of the gene expression process in the cell. It has been studied extensively, but many aspects of the mechanism are not well understood. The widespread availability of experimental RNA-seq data from high-throughput experiments provides a unique opportunity to infer the end of the transcription units genome wide. This data is available for both Rho-dependent and Rho-independent termination pathways that drive transcription termination in bacteria. Our book chapter gives an overview of the current knowledge of Rho-independent transcription termination mechanisms and the prediction approaches currently deployed to infer the termination sites. Thereafter, we describe our method that uses cluster hairpins to detect Rho-independent transcription termination sites. These clusters are a group of hairpins that lies at <15 bp from each other and are together capable of enforcing the termination process. The idea of a group of hairpins being extensively used for transcription termination is new, and results show that at least 52% of the total cases are of this type, while in the remaining cases, a single strong hairpin is capable of driving transcription termination. The reads derived from the RNA-seq data for corresponding bacteria have been used to validate the predicted sites. The predictions that match these RNA-seq derived sites have higher confidence, and we find almost 98% of the predicted sites, including alternate termination sites, to match the RNA-seq data. We discuss the features of predicted hairpins in detail for a better understanding of the Rho-independent transcription termination mechanism in bacteria. We also explain how users can use the tools developed by us to do transcription terminator predictions and design their experiments through genome-level visualization of the transcription termination sites from the precomputed INTERPIN database.
Subject(s)
RNA-Seq , Transcription Termination, Genetic , RNA-Seq/methods , Software , Computational Biology/methods , RNA, Bacterial/genetics , Bacteria/genetics , Sequence Analysis, RNA/methods , Terminator Regions, Genetic/genetics , Gene Expression Regulation, BacterialABSTRACT
A riboswitch generally regulates the expression of its downstream genes through conformational change in its expression platform (EP) upon ligand binding. The cyclic diguanosine monophosphate (c-di-GMP) class I riboswitch Bc1 is widespread and conserved among Bacillus cereus group species. In this study, we revealed that Bc1 has a long EP with two typical ρ-independent terminator sequences 28 bp apart. The upstream terminator T1 is dominant in vitro, while downstream terminator T2 is more efficient in vivo. Through mutation analysis, we elucidated that Bc1 exerts a rare and incoherent "transcription-translation" dual regulation with T2 playing a crucial role. However, we found that Bc1 did not respond to c-di-GMP under in vitro transcription conditions, and the expressions of downstream genes did not change with fluctuation in intracellular c-di-GMP concentration. To explore this puzzle, we conducted SHAPE-MaP and confirmed the interaction of Bc1 with c-di-GMP. This shows that as c-di-GMP concentration increases, T1 unfolds but T2 remains almost intact and functional. The presence of T2 masks the effect of T1 unwinding, resulting in no response of Bc1 to c-di-GMP. The high Shannon entropy values of EP region imply the potential alternative structures of Bc1. We also found that zinc uptake regulator can specifically bind to the dual terminator coding sequence and slightly trigger the response of Bc1 to c-di-GMP. This work will shed light on the dual-regulation riboswitch and enrich our understanding of the RNA world.IMPORTANCEIn nature, riboswitches are involved in a variety of metabolic regulation, most of which preferentially regulate transcription termination or translation initiation of downstream genes in specific ways. Alternatively, the same or different riboswitches can exist in tandem to enhance regulatory effects or respond to multiple ligands. However, many putative conserved riboswitches have not yet been experimentally validated. Here, we found that the c-di-GMP riboswitch Bc1 with a long EP could form a dual terminator and exhibit non-canonical and incoherent "transcription-translation" dual regulation. Besides, zinc uptake regulator specifically bound to the coding sequence of the Bc1 EP and slightly mediated the action of Bc1. The application of SHAPE-MaP to the dual regulation mechanism of Bc1 may establish the foundation for future studies of such complex untranslated regions in other bacterial genomes.
Subject(s)
Bacillus thuringiensis , Cyclic GMP , Gene Expression Regulation, Bacterial , Riboswitch , Riboswitch/genetics , Cyclic GMP/analogs & derivatives , Cyclic GMP/metabolism , Cyclic GMP/genetics , Bacillus thuringiensis/genetics , Bacillus thuringiensis/metabolism , Nucleic Acid Conformation , Transcription, Genetic , Terminator Regions, Genetic/genetics , Bacterial Proteins/genetics , Bacterial Proteins/metabolism , RNA, Bacterial/genetics , RNA, Bacterial/metabolismABSTRACT
The objectives of the present study were to identify the cis-elements of the promoter absolutely required for the efficient rat NHE3 gene transcription and to locate positive and negative regulatory elements in the 5’-flanking sequence (5’FS), which might modulate the gene expression in proximal tubules, and to compare this result to those reported for intestinal cell lines. We analyzed the promoter activity of different 5’FS segments of the rat NHE3 gene, in the OKP renal proximal tubule cell line by measuring the activity of the reporter gene luciferase. Because the segment spanning the first 157 bp of 5’FS was the most active it was studied in more detail by sequential deletions, point mutations, and gel shift assays. The essential elements for gene transcription are in the region -85 to -33, where we can identify consensual binding sites for Sp1 and EGR-1, which are relevant to NHE3 gene basal transcription. Although a low level of transcription is still possible when the first 25 bp of the 5’FS are used as promoter, efficient transcription only occurs with 44 bp of 5’FS. There are negative regulatory elements in the segments spanning -1196 to -889 and -467 to -152, and positive enhancers between -889 and -479 bp of 5’FS. Transcription factors in the OKP cell nuclear extract efficiently bound to DNA elements of rat NHE3 promoter as demonstrated by gel shift assays, suggesting a high level of similarity between transcription factors of both species, including Sp1 and EGR-1.