Pan-cancer tRNA-derived fragment CAT1 coordinates RBPMS to stabilize NOTCH2 mRNA to promote tumorigenesis.

Transfer RNA-derived fragments (tRFs) are a class of small non-coding regulatory RNAs that are involved in the pathophysiology of many diseases. However, the role of tRFs in cancer progression remains largely elusive. Here, we demonstrate that a pan-cancer 3'-tRF, CAT1 (cancer associated tRF 1), is ubiquitously upregulated in tumors and associated with poor prognosis of a variety of cancers, including lung cancer. The upregulated CAT1 in cancer cells binds to RNA-binding protein with multiple splicing (RBPMS) and displaces NOTCH2 association from RBPMS, thereby inhibiting the subsequent CCR4-NOT deadenylation-complex-mediated NOTCH2 mRNA decay. The CAT1-enhanced NOTCH2 expression promotes lung cancer cell proliferation and metastasis in vitro and in vivo. In addition, plasma CAT1 levels are substantially increased in patients with lung cancer compared to non-cancer control subjects. Our findings reveal an intrinsic connection between cancer-specific upregulation of CAT1 and cancer progression, show the regulation of NOTCH signaling in cancer by a 3'-tRF, and highlight its great clinical potential.

OncotRF: an online resource for exploration of tRNA-derived fragments in human cancers.

Transfer RNA-derived fragments (tRFs) are a new class of small non-coding RNAs whose biological roles in cancers are not well understood. Emerging evidence suggests that tRFs are involved in gene regulation at multiple levels. In this study, we constructed an integrative database, OncotRF (http://bioinformatics.zju.edu.cn/OncotRF), for in silico exploration of tRF functions, and identification of diagnostic and prognostic biomarkers in cancers. The database contains an analysis pipeline for tRF identification and characterization, analysis results of 11,211 small RNA sequencing samples and 8,776 RNA sequencing samples, and clinicopathologic annotation data from The Cancer Genome Atlas (TCGA). The results include: tRF identification and quantification across 33 cancers, abnormally expressed tRFs and genes, tRF-gene correlations, tRF-gene networks, survival analyses, and tRF-related functional enrichment analyses. Users are also able to identify differentially expressed tRFs, predict their functions, and assess the relevance of the tRF expression levels to the clinical outcome according to user-defined groups. Additionally, an online Kaplan-Meier plotter is available in OncotRF for plotting survival curves according to user-defined groups. OncotRF will be a valuable online database and functional annotation tool for researchers studying the roles, functions, and mechanisms of tRFs in human cancers.

Gossypium hirsutum Salt Tolerance Is Enhanced by Overexpression of G. arboreum JAZ1.

Gossypium arboreum possesses many favorable traits including robust defense against biotic and abiotic stress although it has been withdrawn from the market because of lower yield and fiber quality compared to G. hirsutum (upland cotton). It is therefore important to explore and utilize the beneficial genes of G. arboretum for G. hirsutum cultivar breeding. Here, the function of G. arboreum JAZ1 in tolerance to salt stress was determined through loss-of-function analysis. GaJAZ1can interact with GaMYC2 to repress expression of downstream genes whose promoters contain a G-box cis element, affecting plant tolerance to salinity stress. The experimental data from NaCl treatments and a 2 year continuous field trial with natural saline-alkaline soil showed that the ectopically overexpressed GaJAZ1 significantly increased salt tolerance in upland cotton compared to the wild type, showing higher growth vigor with taller plants, increased fresh weight, and more bolls, which is due to reprogrammed expression of tolerance-related genes and promotion of root development. High-throughput RNA sequencing of GaJAZ1 transgenic and wild-type plants showed many differentially expressed genes involved in JA signaling and biosynthesis, salt stress-related genes, and hormone-related genes, suggesting that overexpressing GaJAZ1 can reprogram the expression of defense-related genes in G. hirsutum plants to increase tolerance to salt stress. The research provides a foundation to explore and utilize favorable genes from Gossypium species for upland cotton cultivar breeding.

Global identification and characterization of tRNA-derived RNA fragment landscapes across human cancers.

Transfer RNA-derived RNA fragments (tRFs) are a class of small non-coding RNAs that are abundant in many organisms, but their role in cancer has not been fully explored. Here, we report a functional genomic landscape of tRFs in 8118 specimens across 15 cancer types from The Cancer Genome Atlas. These tRFs exhibited characteristics of widespread expression, high sequence conservation, cytoplasmic localization, specific patterns of tRNA cleavage and conserved cleavage in tissues. A cross-tumor analysis revealed significant commonality among tRF expression subtypes from distinct tissues of origins, characterized by upregulation of a group of tRFs with similar size and activation of cancer-associated signaling. One of the largest superclusters was composed of 22 nt 3'-tRFs upregulated in 13 cancer types, all of which share the activation of Ras/MAPK, RTK and TSC/mTOR signaling. tRF-based subgrouping provided clinically relevant stratifications and significantly improved outcome prediction by incorporating clinical variables. Additionally, we discovered 11 cancer driver tRFs using an effective approach for accurately exploring cross-tumor and platform trends. As a proof of concept, we performed comprehensive functional assays on a non-microRNA driver tRF, 5'-IleAAT-8-1-L20, and validated its oncogenic roles in lung cancer in vitro and in vivo. Our study also provides a valuable tRF resource for identifying diagnostic and prognostic biomarkers, developing cancer therapy and studying cancer pathogenesis.

Circ2Disease: a manually curated database of experimentally validated circRNAs in human disease.

Circular RNAs (circRNAs), a new class of regulatory noncoding RNAs, play important roles in human diseases. While a growing number of circRNAs have been characterized with biological functions, it is necessary to integrate all the information to facilitate studies on circRNA functions and regulatory networks in human diseases. Circ2Disease database contains 273 manually curated associations between 237 circRNAs and 54 human diseases with strong experimental evidence from 120 studies. Each association includes circRNA name, disease name, expression pattern, experimental method, a brief functional description of the circRNA-disease relationship, and other detailed information. The experimentally validated miRNAs that may be 'sponged up' by these circRNAs and their validated targets were also integrated to form a comprehensive regulatory network. Circ2Disease provides a user-friendly interface to browse, search, analyze regulatory network and download data. With the rapidly increasing interest in circRNAs, Circ2Disease will significantly improve our understanding of circRNA deregulation in diseases and is a useful resource for studying posttranscriptional regulatory roles of circRNAs in human diseases.

ccNET: Database of co-expression networks with functional modules for diploid and polyploid Gossypium.

Plant genera with both diploid and polyploid species are a common evolutionary occurrence. Polyploids, especially allopolyploids such as cotton and wheat, are a great model system for heterosis research. Here, we have integrated genome sequences and transcriptome data of Gossypium species to construct co-expression networks and identified functional modules from different cotton species, including 1155 and 1884 modules in G. arboreum and G. hirsutum, respectively. We overlayed the gene expression results onto the co-expression network. We further provided network comparison analysis for orthologous genes across the diploid and allotetraploid Gossypium We also constructed miRNA-target networks and predicted PPI networks for both cotton species. Furthermore, we integrated in-house ChIP-seq data of histone modification (H3K4me3) together with cis-element analysis and gene sets enrichment analysis tools for studying possible gene regulatory mechanism in Gossypium species. Finally, we have constructed an online ccNET database (http://structuralbiology.cau.edu.cn/gossypium) for comparative gene functional analyses at a multi-dimensional network and epigenomic level across diploid and polyploid Gossypium species. The ccNET database will be beneficial for community to yield novel insights into gene/module functions during cotton development and stress response, and might be useful for studying conservation and diversity in other polyploid plants, such as T. aestivum and Brassica napus.

Co-expression network analyses identify functional modules associated with development and stress response in Gossypium arboreum.

Cotton is an economically important crop, essential for the agriculture and textile industries. Through integrating transcriptomic data, we discovered that multi-dimensional co-expression network analysis was powerful for predicting cotton gene functions and functional modules. Here, the recently available transcriptomic data on Gossypium arboreum, including data on multiple growth stages of tissues and stress treatment samples were applied to construct a co-expression network exploring multi-dimensional expression (development and stress) through multi-layered approaches. Based on differential gene expression and network analysis, a fibre development regulatory module of the gene GaKNL1 was found to regulate the second cell wall through repressing the activity of REVOLUTA, and a tissue-selective module of GaJAZ1a was examined in response to water stress. Moreover, comparative genomics analysis of the JAZ1-related regulatory module revealed high conservation across plant species. In addition, 1155 functional modules were identified through integrating the co-expression network, module classification and function enrichment tools, which cover functions such as metabolism, stress responses, and transcriptional regulation. In the end, an online platform was built for network analysis (http://structuralbiology.cau.edu.cn/arboreum), which could help to refine the annotation of cotton gene function and establish a data mining system to identify functional genes or modules with important agronomic traits.

mRNA-seq analysis of the Gossypium arboreum transcriptome reveals tissue selective signaling in response to water stress during seedling stage.

The cotton diploid species, Gossypium arboreum, shows important properties of stress tolerance and good genetic stability. In this study, through mRNA-seq, we de novo assembled the unigenes of multiple samples with 3h H(2)O, NaCl, or PEG treatments in leaf, stem and root tissues and successfully obtained 123,579 transcripts of G. arboreum, 89,128 of which were with hits through BLAST against known cotton ESTs and draft genome of G. raimondii. About 36,961 transcripts (including 1,958 possible transcription factor members) were identified with differential expression under water stresses. Principal component analysis of differential expression levels in multiple samples suggested tissue selective signalling responding to water stresses. Venn diagram analysis showed the specificity and intersection of transcripts' response to NaCl and PEG treatments in different tissues. Self-organized mapping and hierarchical cluster analysis of the data also revealed strong tissue selectivity of transcripts under salt and osmotic stresses. In addition, the enriched gene ontology (GO) terms for the selected tissue groups were differed, including some unique enriched GO terms such as photosynthesis and tetrapyrrole binding only in leaf tissues, while the stem-specific genes showed unique GO terms related to plant-type cell wall biogenesis, and root-specific genes showed unique GO terms such as monooxygenase activity. Furthermore, there were multiple hormone cross-talks in response to osmotic and salt stress. In summary, our multidimensional mRNA sequencing revealed tissue selective signalling and hormone crosstalk in response to salt and osmotic stresses in G. arboreum. To our knowledge, this is the first such report of spatial resolution of transcriptome analysis in G. arboreum. Our study will potentially advance understanding of possible transcriptional networks associated with water stress in cotton and other crop species.

Comparative genomic analysis of NAC transcriptional factors to dissect the regulatory mechanisms for cell wall biosynthesis.

BACKGROUND: NAC domain transcription factors are important transcriptional regulators involved in plant growth, development and stress responses. Recent studies have revealed several classes of NAC transcriptional factors crucial for controlling secondary cell wall biosynthesis. These transcriptional factors mainly include three classes, SND, NST and VND. Despite progress, most current analysis is carried out in the model plant Arabidopsis. Moreover, many downstream genes regulated by these transcriptional factors are still not clear. METHODS: In order to identify the key homologue genes across species and discover the network controlling cell wall biosynthesis, we carried out comparative genome analysis of NST, VND and SND genes across 19 higher plant species along with computational modelling of genes regulated or co-regulated with these transcriptional factors. RESULTS: The comparative genome analysis revealed that evolutionarily the secondary-wall-associated NAC domain transcription factors first appeared in Selaginella moellendorffii. In fact, among the three groups, only VND genes appeared in S. moellendorffii, which is evolutionarily earlier than the other two groups. The Arabidopsis and rice gene expression analysis showed specific patterns of the secondary cell wall-associated NAC genes (SND, NST and VND). Most of them were preferentially expressed in the stem, especially the second internodes. Furthermore, comprehensive co-regulatory network analysis revealed that the SND and MYB genes were co-regulated, which indicated the coordinative function of these transcriptional factors in modulating cell wall biosynthesis. In addition, the co-regulatory network analysis revealed many novel genes and pathways that could be involved in cell wall biosynthesis and its regulation. The gene ontology analysis also indicated that processes like carbohydrate synthesis, transport and stress response, are coordinately regulated toward cell wall biosynthesis. CONCLUSIONS: Overall, we provided a new insight into the evolution and the gene regulatory network of a subgroup of the NAC gene family controlling cell wall composition through bioinformatics data mining and bench validation. Our work might benefit to elucidate the possible molecular mechanism underlying the regulation network of secondary cell wall biosynthesis.

Transcriptome analysis reveals salt-stress-regulated biological processes and key pathways in roots of cotton (Gossypium hirsutum L.).

High salinity is one of the main factors limiting cotton growth and productivity. The genes that regulate salt stress in TM-1 upland cotton were monitored using microarray and real-time PCR (RT-PCR) with samples taken from roots. Microarray analysis showed that 1503 probe sets were up-regulated and 1490 probe sets were down-regulated in plants exposed for 3h to 100mM NaCl, and RT-PCR analysis validated 42 relevant/related genes. The distribution of enriched gene ontology terms showed such important processes as the response to water stress and pathways of hormone metabolism and signal transduction were induced by the NaCl treatment. Some key regulatory gene families involved in abiotic and biotic sources of stress such as WRKY, ERF, and JAZ were differentially expressed. Our transcriptome analysis might provide some useful insights into salt-mediated signal transduction pathways in cotton and offer a number of candidate genes as potential markers of tolerance to salt stress.

Gene expression profiles deciphering rice phenotypic variation between Nipponbare (Japonica) and 93-11 (Indica) during oxidative stress.

Rice is a very important food staple that feeds more than half the world's population. Two major Asian cultivated rice (Oryza sativa L.) subspecies, japonica and indica, show significant phenotypic variation in their stress responses. However, the molecular mechanisms underlying this phenotypic variation are still largely unknown. A common link among different stresses is that they produce an oxidative burst and result in an increase of reactive oxygen species (ROS). In this study, methyl viologen (MV) as a ROS agent was applied to investigate the rice oxidative stress response. We observed that 93-11 (indica) seedlings exhibited leaf senescence with severe lesions under MV treatment compared to Nipponbare (japonica). Whole-genome microarray experiments were conducted, and 1,062 probe sets were identified with gene expression level polymorphisms between the two rice cultivars in addition to differential expression under MV treatment, which were assigned as Core Intersectional Probesets (CIPs). These CIPs were analyzed by gene ontology (GO) and highlighted with enrichment GO terms related to toxin and oxidative stress responses as well as other responses. These GO term-enriched genes of the CIPs include glutathine S-transferases (GSTs), P450, plant defense genes, and secondary metabolism related genes such as chalcone synthase (CHS). Further insertion/deletion (InDel) and regulatory element analyses for these identified CIPs suggested that there may be some eQTL hotspots related to oxidative stress in the rice genome, such as GST genes encoded on chromosome 10. In addition, we identified a group of marker genes individuating the japonica and indica subspecies. In summary, we developed a new strategy combining biological experiments and data mining to study the possible molecular mechanism of phenotypic variation during oxidative stress between Nipponbare and 93-11. This study will aid in the analysis of the molecular basis of quantitative traits.