Eukaryotic Transcription Research Articles

C-terminal determinants for RNA binding motif 7 protein stability and RNA recognition

The 7SK ribonucleoprotein (RNP) is a critical regulator of eukaryotic transcription. Recently, RNA binding motif 7 (RBM7) containing an RNA recognition motif (RRM) was reported to associate with 7SK RNA and core 7SK RNP protein components in response to DNA damage. However, little is known about the mode of RBM7-7SK RNA recognition. Here, we found that RRM constructs containing extended C-termini have increased solubility compared to a minimal RRM construct, although these constructs aggregate in a temperature and concentration-dependent manner. Using solution NMR dynamics experiments, we identified additional structural features observed previously in crystal but not in solution structures. To identify potential RBM7-7SK RNA binding sites, we analyzed deposited data from in cellulo crosslinking experiments and found that RBM7 primarily crosslinks to the distal region of 7SK stem-loop 3 (SL3). Electrophoretic mobility shift assays and NMR chemical shift perturbation experiments showed weak binding to 7SK SL3 constructs in vitro. Together, these results provide new insights into RBM7 RRM folding and recognition of 7SK RNA.

The new uORFdb: integrating literature, sequence, and variation data in a central hub for uORF research.

Upstream open reading frames (uORFs) are initiated by AUG or near-cognate start codons and have been identified in the transcript leader sequences of the majority of eukaryotic transcripts. Functionally, uORFs are implicated in downstream translational regulation of the main protein coding sequence and may serve as a source of non-canonical peptides. Genetic defects in uORF sequences have been linked to the development of various diseases, including cancer. To simplify uORF-related research, the initial release of uORFdb in 2014 provided a comprehensive and manually curated collection of uORF-related literature. Here, we present an updated sequence-based version of uORFdb, accessible at https://www.bioinformatics.uni-muenster.de/tools/uorfdb. The new uORFdb enables users to directly access sequence information, graphical displays, and genetic variation data for over 2.4 million human uORFs. It also includes sequence data of >4.2 million uORFs in 12 additional species. Multiple uORFs can be displayed in transcript- and reading-frame-specific models to visualize the translational context. A variety of filters, sequence-related information, and links to external resources (UCSC Genome Browser, dbSNP, ClinVar) facilitate immediate in-depth analysis of individual uORFs. The database also contains uORF-related somatic variation data obtained from whole-genome sequencing (WGS) analyses of 677 cancer samples collected by the TCGA consortium.

Structures of +1 nucleosome-bound PIC-Mediator complex.

RNA polymerase II-mediated eukaryotic transcription starts with the assembly of the preinitiation complex (PIC) on core promoters. The +1 nucleosome is well positioned about 40 base pairs downstream of the transcription start site (TSS) and is commonly known as a barrier of transcription. The +1 nucleosome-bound PIC-Mediator structures show that PIC-Mediator prefers binding to T40N nucleosome located 40 base pairs downstream of TSS and contacts T50N but not the T70N nucleosome. The nucleosome facilitates the organization of PIC-Mediator on the promoter by binding TFIIH subunit p52 and Mediator subunits MED19 and MED26 and may contribute to transcription initiation. PIC-Mediator exhibits multiple nucleosome-binding patterns, supporting a structural role of the +1 nucleosome in the coordination of PIC-Mediator assembly. Our study reveals the molecular mechanism of PIC-Mediator organization on chromatin and underscores the significance of the +1 nucleosome in regulating transcription initiation.

Plant transcripts with long or structured upstream open reading frames in the NDL2 5' UTR can escape nonsense-mediated mRNA decay in a reinitiation-independent manner.

Many eukaryotic transcripts contain upstream open reading frames (uORFs). Translated uORFs can inhibit the translation of main ORFs by imposing the need for reinitiation of translation. Translated uORFs can also lead to transcript degradation by the nonsense-mediated mRNA decay (NMD) pathway. In mammalian cells, translated uORFs were shown to target their transcripts to NMD if the uORFs were long (>23-32 amino acids), structured, or inhibit reinitiation. Reinitiation was shown to rescue uORF-containing mammalian transcripts from NMD. Much less is known about the significance of the length, structure, and reinitiation efficiency of translated uORFs for NMD targeting in plants. Although high-throughput studies suggested that uORFs do not globally reduce plant transcript abundance, it was not clear whether this was due to NMD-escape-permitting parameters of uORF recognition, length, structure, or reinitiation efficiency. We expressed in Arabidopsis reporter genes that included NDL2 5' untranslated region and various uORFs with modulation of the above parameters. We found that transcripts can escape NMD in plants even when they include efficiently translated uORFs up to 70 amino acids long, or structured uORFs, in the absence of reinitiation. These data highlight an apparent difference between the rules that govern the exposure of uORF-containing transcripts to NMD in mammalian and plant cells.

Genome-wide identification of Arabidopsis non-AUG-initiated upstream ORFs with evolutionarily conserved regulatory sequences that control protein expression levels.

This study identified four novel regulatory non-AUG-initiated upstream ORFs (uORFs) with evolutionarily conserved sequences in Arabidopsis and elucidated the mechanism by which a non-AUG-initiated uORF promotes main ORF translation. Upstream open reading frames (uORFs) are short ORFs found in the 5'-untranslated regions (5'-UTRs) of eukaryotic transcripts and can influence the translation of protein-coding main ORFs (mORFs). Recent genome-wide ribosome profiling studies have revealed that hundreds or thousands of uORFs initiate translation at non-AUG start codons. However, the physiological significance of these non-AUG uORFs has so far been demonstrated for only a few of them. In this study, to identify physiologically important regulatory non-AUG uORFs in Arabidopsis, we took an approach that combined bioinformatics and experimental analysis. Since physiologically important non-AUG uORFs are likely to be conserved across species, we first searched the Arabidopsis genome for non-AUG-initiated uORFs with evolutionarily conserved sequences. Then, we examined the effects of the conserved non-AUG uORFs on the expression of the downstream mORFs using transient expression assays. As a result, three inhibitory and one promotive non-AUG uORFs were identified. Among the inhibitory non-AUG uORFs, two exerted repressive effects on mORF expression in an amino acid sequence-dependent manner. These two non-AUG uORFs are likely to encode regulatory peptides that cause ribosome stalling, thereby enhancing their repressive effects. In contrast, one of the identified regulatory non-AUG uORFs promoted mORF expression by alleviating the inhibitory effect of a downstream AUG-initiated uORF. These findings provide insights into the mechanisms that enable non-AUG uORFs to play regulatory roles despite their low translation initiation efficiencies.

Molecular biology of the cell

Molecular biology of the cell

Suppressor mutations that make the essential transcription factor Spn1/Iws1 dispensable in Saccharomyces cerevisiae.

Spn1/Iws1 is an essential eukaryotic transcription elongation factor that is conserved from yeast to humans as an integral member of the RNA polymerase II elongation complex. Several studies have shown that Spn1 functions as a histone chaperone to control transcription, RNA splicing, genome stability, and histone modifications. However, the precise role of Spn1 is not understood, and there is little understanding of why it is essential for viability. To address these issues, we have isolated 8 suppressor mutations that bypass the essential requirement for Spn1 in Saccharomyces cerevisiae. Unexpectedly, the suppressors identify several functionally distinct complexes and activities, including the histone chaperone FACT, the histone methyltransferase Set2, the Rpd3S histone deacetylase complex, the histone acetyltransferase Rtt109, the nucleosome remodeler Chd1, and a member of the SAGA coactivator complex, Sgf73. The identification of these distinct groups suggests that there are multiple ways in which Spn1 bypass can occur, including changes in histone acetylation and alterations in other histone chaperones. Thus, Spn1 may function to overcome repressive chromatin by multiple mechanisms during transcription. Our results suggest that bypassing a subset of these functions allows viability in the absence of Spn1.

SHP-1 dephosphorylates histone H2B to facilitate its ubiquitination during transcription.

Dynamic regulation of phosphorylation and dephosphorylation of histones is essential for eukaryotic transcription, but the enzymes engaged in histone dephosphorylation are not fully explored. Here, we show that the tyrosine phosphatase SHP-1 dephosphorylates histone H2B and plays a critical role during transition from the initiation to the elongation stage of transcription. Nuclear-localized SHP-1 is associated with the Paf1 complex at chromatin and dephosphorylates H2B at tyrosine 121. Moreover, knockout of SHP-1, or expression of a mutant mimicking constitutive phosphorylation of H2B Y121, leads to a reduction in genome-wide H2B ubiquitination, which subsequently causes defects in RNA polymerase II-dependent transcription. Mechanistically, we demonstrate that Y121 phosphorylation precludes H2B's interaction with the E2 enzyme, indicating that SHP-1-mediated dephosphorylation of this residue may be a prerequisite for efficient H2B ubiquitination. Functionally, we find that SHP-1-mediated H2B dephosphorylation contributes to maintaining basal autophagic flux in cells through the efficient transcription of autophagy and lysosomal genes. Collectively, our study reveals an important modification of histone H2B regulated by SHP-1 that has a role during eukaryotic transcription.

Identification and functional annotation of long intergenic non-coding RNAs in Brassicaceae.

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.

Identification of aberrant transcription termination at specific gene loci with DNA hypomethylated transcription termination sites caused by DNA methyltransferase deficiency.

CpG methylation of genomic DNA is a well-known repressive epigenetic marker in eukaryotic transcription, and DNA methylation of promoter regions is correlated with gene silencing. In contrast to the promoter regions, the function of DNA methylation during transcription termination remains to be elucidated. A recent study revealed that mouse DNA methyltransferase 3a (Dnmt3a) mainly functions in de novo methylation in the promoter and gene body regions, including transcription termination sites (TTSs), during development. To investigate the relationship between DNA methylation overlapping the TTSs and transcription termination, we performed bioinformatics analysis using six pre-existing Dnmt-/- mouse cell datasets: four types of neurons (three Dnmt3a-/- and one Dnmt1-/- mutants) and two types of embryonic fibroblasts (MEFs) (Dnmt3a-/- and Dnmt3b-/- mutants). Combined analyses using methylome and transcriptome data revealed that read counts downstream of hypomethylated TTSs were increased in three types of neurons (two Dnmt3a-/- and one Dnmt1-/- mutants). Among these, an increase in chimeric transcripts downstream of the TTSs was observed in Dnmt3a-/- mature olfactory sensory neurons and Dnmt3a-/- agouti-related peptide (protein)-producing neurons, thereby indicating that read-through occurs in hypomethylated TTSs at specific gene loci in these two mutants. Conversely, in Dnmt3a-/- MEFs, we detected reductions in read counts downstream of hypomethylated TTSs. These results indicate that the hypomethylation of TTSs can both positively and negatively regulate transcription termination, dependent on Dnmt and cell types. This study is the first to identify the aberrant termination of transcription at specific gene loci with DNA hypomethylated TTSs attributable to Dnmt deficiency.

Glutamate-aspartate transporter 1 attenuates oxygen-glucose deprivation-induced injury by promoting glutamate metabolism in primary cortical neurons.

Ischemic stroke is a common cerebral disease. However, the treatment for the disease is limited. Daurian ground squirrel (GS;Spermophilus dauricus), a hibernating mammalian species, is highly tolerant to ischemia. In the present study, GS neurons in a non-hibernating state were found to be more resistant to oxygen-glucose deprivation (OGD), an ischemic model in vitro. We leveraged the differences in the endurance capacity of GS and rats to investigate the mechanisms of resistance to ischemia in GSneurons. We first identified glutamate-aspartate transporter 1 (GLAST) as a cytoprotective factor that contributed to tolerance against OGD injury of GS neurons. The expression of GLAST in GS neurons was much higher than that in rat neurons. Overexpression of GLAST rescued viability in rat neurons, and GS neurons exhibited decreased viability following GLAST knockdown under OGD conditions. Mechanistically, more glutamate was transported into neurons after GLAST overexpression and served as substrates for ATP production. Furthermore, eukaryotic transcription initiation factor 4E binding protein 1 was downregulated by GLAST to rescue neuronal viability. Our findings not only revealed an important molecular mechanism underlying the survival of hibernating mammals but also suggested that neuronal GLAST may be a potential target for ischemic stroke therapy.

Structural insights into nuclear transcription by eukaryotic DNA-dependent RNA polymerases.

The eukaryotic transcription apparatus synthesizes a staggering diversity of RNA molecules. The labour of nuclear gene transcription is, therefore, divided among multiple DNA-dependent RNA polymerases. RNA polymerase I (Pol I) transcribes ribosomal RNA, Pol II synthesizes messenger RNAs and various non-coding RNAs (including long non-coding RNAs, microRNAs and small nuclear RNAs) and Pol III produces transfer RNAs and other short RNA molecules. Pol I, Pol II and Pol III are large, multisubunit protein complexes that associate with a multitude of additional factors to synthesize transcripts that largely differ in size, structure and abundance. The three transcription machineries share common characteristics, but differ widely in various aspects, such as numbers of RNA polymerase subunits, regulatory elements and accessory factors, which allows them to specialize in transcribing their specific RNAs. Common to the three RNA polymerases is that the transcription process consists of three major steps: transcription initiation, transcript elongation and transcription termination. In this Review, we outline the common principles and differences between the Pol I, Pol II and Pol III transcription machineries and discuss key structural and functional insights obtained into the three stages of their transcription processes.

Abstract 333: Tnik: A Novel Regulator Of Endothelial Interferon Signaling And Activation

Background: Endothelial cells (ECs) are prominent cell types governing vascular tone and homeostasis. Activation of ECs is the central pathophysiological process in the initiation of atherosclerosis (AS). Interferon signaling (IFN Type I and II) is activated in ECs under inflammatory conditions such as systemic lupus erythematosus and is implicated in EC activation. But the underlying mechanisms are still unknown. Through RNA sequencing and pathway analysis, we have identified TNIK, a Ste20 family of MAP kinases (MAP4K) to regulate IFN signaling in ECs. Our study is the first to report the role and function of TNIK in IFN signaling and EC activation. Hypothesis: TNIK regulates IFNs through STAT mediated signaling and downstream interferon stimulated gene (ISG) sets Methods: Cultured Human aortic endothelial cells (HAECs) were treated with control or TNIK-targeted siRNA. Later cells were subjected to static, laminar or disturbed flow conditions. Total RNA was isolated and sequenced by mRNA-seq. The differential expression data were analyzed using CLC genomics and Ingenuity Pathway Analysis (IPA). Western blotting and qRT-PCR were performed to validate STAT activation and interferon stimulated genes (ISGs) Results: We found >450 genes differentially expressed in siTNIK vs siCTRL treated ECs (fold change>2, FDR<0.05). IPA revealed that under static (ST) and disturbed flow (DF) conditions TNIK depletion lead to inhibition IFN and upregulation of eukaryotic transcription elongation factor signaling (eIF2). In addition, under DF, TNIK depletion lead to downregulation of hyperctokinemia/hyperchemokinemia. Validation through western blotting and qRT-PCR confirmed that TNIK regulates STAT expression and IFN stimulated target genes. Conclusions: Our study elucidates the unique role of TNIK in regulating STAT mediated IFN signaling and subsequent ISGs in ECs.

Phosphorylation‐specific recruitment of human SCAF6 to RNA polymerase II during transcription regulation

The C‐terminal domain (CTD) of human RNA polymerase II (RNAPII) acts as a regulatory hub for RNA processing and progression of eukaryotic transcription. Transcriptional control is fine‐tuned through dynamic phosphorylation of the 52 repetitive heptads on the CTD to generate post‐translational modification patterns that precisely recruit cellular proteins according to transcriptional stages. However, it remains unclear what transcription regulators are differentially recruited by specific phospho‐marks. We assembled an enzymatic toolbox to recapitulate the dynamic phosphorylation patterns of CTD during transcription. We mapped the individual modification sites in high precision with advanced mass spectrometry strategies such as ultraviolet photodissociation. Utilizing a label‐free quantitative proteomics approach, we identified proteins that interact with different phosphoisoforms of Pol II CTD. Specifically, we identified RNA‐binding protein SCAF6 that are recruited for transcription by Ser2 and Thr4 phosphorylation. We characterized the interaction of SCAF6 with differentially phosphorylated CTD. The timely recruitment of SCAF6 is then investigated in cells using genome‐wide sequencing. In conclusion, we have expanded the repertoire of RNA‐binding proteins that directly interact with the CTD of RNA Pol II and illuminate the concerted interplay between splicing and transcription.

Structural insights into cooperative DNA recognition by the CCAAT-binding complex and its bZIP transcription factor HapX

The heterotrimeric CCAAT-binding complex (CBC) is a fundamental eukaryotic transcription factor recognizing the CCAAT box. In certain fungi, like Aspergilli, the CBC cooperates with the basic leucine zipper HapX to control iron metabolism. HapX functionally depends on the CBC, and the stable interaction of both requires DNA. To study this cooperative effect, X-ray structures of the CBC-HapX-DNA complex were determined. Downstream of the CCAAT box, occupied by the CBC, a HapX dimer binds to the major groove. The leash-like N terminus of the distal HapX subunit contacts the CBC, and via a flexible polyproline type II helix mediates minor groove interactions that stimulate sequence promiscuity. Invitro and invivo mutagenesis suggest that the structural and functional plasticity of HapX results from local asymmetry and its ability to target major and minor grooves simultaneously. The latter feature may also apply to related transcription factors such as yeast Hap4 and distinct Yap family members.

Nuclear and cytoplasmic poly(A) binding proteins (PABPs) favor distinct transcripts and isoforms.

The poly(A)-tail appended to the 3′-end of most eukaryotic transcripts plays a key role in their stability, nuclear transport, and translation. These roles are largely mediated by Poly(A) Binding Proteins (PABPs) that coat poly(A)-tails and interact with various proteins involved in the biogenesis and function of RNA. While it is well-established that the nuclear PABP (PABPN) binds newly synthesized poly(A)-tails and is replaced by the cytoplasmic PABP (PABPC) on transcripts exported to the cytoplasm, the distribution of transcripts for different genes or isoforms of the same gene on these PABPs has not been investigated on a genome-wide scale. Here, we analyzed the identity, splicing status, poly(A)-tail size, and translation status of RNAs co-immunoprecipitated with endogenous PABPN or PABPC in human cells. At steady state, many protein-coding and non-coding RNAs exhibit strong bias for association with PABPN or PABPC. While PABPN-enriched transcripts more often were incompletely spliced and harbored longer poly(A)-tails and PABPC-enriched RNAs had longer half-lives and higher translation efficiency, there are curious outliers. Overall, our study reveals the landscape of RNAs bound by PABPN and PABPC, providing new details that support and advance the current understanding of the roles these proteins play in poly(A)-tail synthesis, maintenance, and function.

Does the Pachytene Checkpoint, a Feature of Meiosis, Filter Out Mistakes in Double-Strand DNA Break Repair and as a side-Effect Strongly Promote Adaptive Speciation?

SynopsisThis essay aims to explain two biological puzzles: why eukaryotic transcription units are composed of short segments of coding DNA interspersed with long stretches of non-coding (intron) DNA, and the near ubiquity of sexual reproduction. As is well known, alternative splicing of its coding sequences enables one transcription unit to produce multiple variants of each encoded protein. Additionally, padding transcription units with non-coding DNA (often many thousands of base pairs long) provides a readily evolvable way to set how soon in a cell cycle the various mRNAs will begin being expressed and the total amount of mRNA that each transcription unit can make during a cell cycle. This regulation complements control via the transcriptional promoter and facilitates the creation of complex eukaryotic cell types, tissues, and organisms. However, it also makes eukaryotes exceedingly vulnerable to double-strand DNA breaks, which end-joining break repair pathways can repair incorrectly. Transcription units cover such a large fraction of the genome that any mis-repair producing a reorganized chromosome has a high probability of destroying a gene. During meiosis, the synaptonemal complex aligns homologous chromosome pairs and the pachytene checkpoint detects, selectively arrests, and in many organisms actively destroys gamete-producing cells with chromosomes that cannot adequately synapse; this creates a filter favoring transmission to the next generation of chromosomes that retain the parental organization, while selectively culling those with interrupted transcription units. This same meiotic checkpoint, reacting to accidental chromosomal reorganizations inflicted by error-prone break repair, can, as a side effect, provide a mechanism for the formation of new species in sympatry. It has been a long-standing puzzle how something as seemingly maladaptive as hybrid sterility between such new species can arise. I suggest that this paradox is resolved by understanding the adaptive importance of the pachytene checkpoint, as outlined above.

Effects of polyphenols on ncRNAs in cancer-An update.

In recent years, oncotherapy has received considerable attention concerning plant polyphenols. Increasing evidence suggests that because of the efficiency of polyphenols, they may have anti-tumour effects in various cancers. However, their regulatory structures remain elusive. Long non-coding RNAs (lncRNAs) have been identified in the regulation of various forms of tumorigenesis and tumour development. Long non-coding RNAs have recently emerged as regulatory eukaryotic transcripts and therapeutic targets with important and diverse functions in health and diseases. LncRNAs may be associated with the initiation, development, and progression of cancer. This review summarizes the research on the modulatory effects of IncRNAs and their roles in mediating cellular processes. The mechanisms of action of polyphenols underlying their therapeutic effects on cancers are also discussed. Based on our review, polyphenols might facilitate a significant epigenetic modification as part of their tissue- and/or cell-related biological effects. This finding may be attributed to their interaction with cellular signalling pathways involved in chronic diseases. Certain lncRNAs might be the target of specific polyphenols, and some critical signalling processes involved in the intervention of cancers might mediate the therapeutic roles of polyphenols.

LINC01535 Attenuates ccRCC Progression through Regulation of the miR-146b-5p/TRIM2 Axis and Inactivation of the PI3K/Akt Pathway.

lncRNAs, a group of eukaryotic cell genome-encoded transcripts, have been demonstrated to exert a notable impact on tumorigenesis. LINC01535, belonging to the lncRNA family, was reported to have an aberrant expression in certain types of cancers and thus affect cancer progression. Nevertheless, the expression pattern and potential roles of LINC01535 in clear cell renal cell carcinoma (ccRCC) remain to be elucidated. Here, LINC01535 expression was detected in ccRCC by RT-qPCR, cell proliferation by CCK-8 assays, and invasion by transwell assays. Besides, effects of LINC01535 on in vivo tumor growth were investigated by xenograft tumor models. The miR-146b-5p/LINC01535/TRIM2 interaction was evaluated via luciferase reporter assays. This study showed downregulation of LINC01535 in ccRCC. Moreover, LINC01535 upregulation attenuated in vitro ccRCC development and hindered in vivo tumor growth. Furthermore, LINC01535 sponged miR-146b-5p which had a negative correlation with LINC01535, and TRIM2 was a direct target of miR-146b-5p and mediated by LINC01535. Mechanically, LINC01535/miR-146b-5p/TRIM2 axis affected ccRCC progression by mediating the PI3K/Akt signaling. All in all, our observations suggest the LINC01535/miR-146b-5p/TRIM2 axis as a crucial role in ccRCC.

Identification and Characterization of circRNAs under Drought Stress in Moso Bamboo (Phyllostachys edulis)

Circular RNAs (circRNAs) are a class of endogenous noncoding RNAs formed by 3′-5′ ligation during splicing. They play an important role in the regulation of transcription and miRNA in eukaryotes. Drought is one of the detrimental abiotic stresses that limit plant growth and productivity. How circRNAs influence the response to drought stress in moso bamboo (Phyllostachys edulis) remains elusive. In this study, we investigate the expression pattern of circRNAs in moso bamboo at 6 h, 12 h, 24 h and 48 h after drought treatment by deep sequencing and bioinformatics analysis and identify 4931 circRNAs, 52 of which are differentially expressed (DEcircRNAs) in drought-treated and untreated moso bamboo. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses of the host genes that generate the DEcircRNAs indcate that these DEcircRNAs are predicted to be involved in biochemical processes in response to drought, such as ubiquitin-mediated proteolysis, calcium-dependent protein kinase phosphorylation, amino acid biosynthesis and plant hormone signal transduction including abscisic acid. In addition, some circRNAs are shown to act as sponges for 291 miRNAs. Taken together, our results characterize the transcriptome profiles of circRNAs in drought responses and provide new insights into resistance breeding of moso bamboo.