A step-by-step guide to analyzing CAGE data using R/Bioconductor

CAGE sequencing reveals CFTR-dependent dysregulation of type I IFN signaling in activated cystic fibrosis macrophages

An intense, nonresolving airway inflammatory response leads to destructive lung disease in cystic fibrosis (CF). Dysregulation of macrophage immune function may be a key facet governing the progression of CF lung disease, but the underlying mechanisms are not fully understood. We used 5' end centered transcriptome sequencing to profile P. aeruginosa LPS-activated human CF macrophages, showing that CF and non-CF macrophages deploy substantially distinct transcriptional programs at baseline and following activation. This includes a significantly blunted type I IFN signaling response in activated patient cells relative to healthy controls that was reversible upon in vitro treatment with CFTR modulators in patient cells and by CRISPR-Cas9 gene editing to correct the F508del mutation in patient-derived iPSC macrophages. These findings illustrate a previously unidentified immune defect in human CF macrophages that is CFTR dependent and reversible with CFTR modulators, thus providing new avenues in the search for effective anti-inflammatory interventions in CF.

Polysome-CAGE of TCL1-driven chronic lymphocytic leukemia revealed multiple N-terminally altered epigenetic regulators and a translation stress signature

The transformation of normal to malignant cells is accompanied by substantial changes in gene expression programs through diverse mechanisms. Here, we examined the changes in the landscape of transcription start sites and alternative promoter (AP) usage and their impact on the translatome in TCL1-driven chronic lymphocytic leukemia (CLL). Our findings revealed a marked elevation of APs in CLL B cells from Eµ-Tcl1 transgenic mice, which are particularly enriched with intra-genic promoters that generate N-terminally truncated or modified proteins. Intra-genic promoter activation is mediated by (1) loss of function of 'closed chromatin' epigenetic regulators due to the generation of inactive N-terminally modified isoforms or reduced expression; (2) upregulation of transcription factors, including c-Myc, targeting the intra-genic promoters and their associated enhancers. Exogenous expression of Tcl1 in MEFs is sufficient to induce intra-genic promoters of epigenetic regulators and promote c-Myc expression. We further found a dramatic translation downregulation of transcripts bearing CNY cap-proximal trinucleotides, reminiscent of cells undergoing metabolic stress. These findings uncovered the role of Tcl1 oncogenic function in altering promoter usage and mRNA translation in leukemogenesis.

Global Analysis of Transcription Start Sites in the New Ovine Reference Genome (Oar rambouillet v1.0)

The overall aim of the Ovine FAANG project is to provide a comprehensive annotation of the new highly contiguous sheep reference genome sequence (Oar rambouillet v1.0). Mapping of transcription start sites (TSS) is a key first step in understanding transcript regulation and diversity. Using 56 tissue samples collected from the reference ewe Benz2616, we have performed a global analysis of TSS and TSS-Enhancer clusters using Cap Analysis Gene Expression (CAGE) sequencing. CAGE measures RNA expression by 5' cap-trapping and has been specifically designed to allow the characterization of TSS within promoters to single-nucleotide resolution. We have adapted an analysis pipeline that uses TagDust2 for clean-up and trimming, Bowtie2 for mapping, CAGEfightR for clustering, and the Integrative Genomics Viewer (IGV) for visualization. Mapping of CAGE tags indicated that the expression levels of CAGE tag clusters varied across tissues. Expression profiles across tissues were validated using corresponding polyA+ mRNA-Seq data from the same samples. After removal of CAGE tags with <10 read counts, 39.3% of TSS overlapped with 5' ends of 31,113 transcripts that had been previously annotated by NCBI (out of a total of 56,308 from the NCBI annotation). For 25,195 of the transcripts, previously annotated by NCBI, no TSS meeting stringent criteria were identified. A further 14.7% of TSS mapped to within 50 bp of annotated promoter regions. Intersecting these predicted TSS regions with annotated promoter regions (±50 bp) revealed 46% of the predicted TSS were "novel" and previously un-annotated. Using whole-genome bisulfite sequencing data from the same tissues, we were able to determine that a proportion of these "novel" TSS were hypo-methylated (32.2%) indicating that they are likely to be reproducible rather than "noise". This global analysis of TSS in sheep will significantly enhance the annotation of gene models in the new ovine reference assembly. Our analyses provide one of the highest resolution annotations of transcript regulation and diversity in a livestock species to date.

Characterization of Arabidopsis thaliana Promoter Bidirectionality and Antisense RNAs by Inactivation of Nuclear RNA Decay Pathways

In animals, RNA polymerase II initiates transcription bidirectionally from gene promoters to produce pre-mRNAs on the forward strand and promoter upstream transcripts (PROMPTs) on the reverse strand. PROMPTs are degraded by the nuclear exosome. Previous studies based on nascent RNA approaches concluded that Arabidopsis (Arabidopsis thaliana) does not produce PROMPTs. Here, we used steady-state RNA sequencing in mutants defective in nuclear RNA decay including the exosome to reassess the existence of Arabidopsis PROMPTs. While they are rare, we identified ∼100 cases of exosome-sensitive PROMPTs in Arabidopsis. Such PROMPTs are sources of small interfering RNAs in exosome-deficient mutants, perhaps explaining why plants have evolved mechanisms to suppress PROMPTs. In addition, we found ∼200 long, unspliced and exosome-sensitive antisense RNAs that arise from transcription start sites within parts of the genome encoding 3'-untranslated regions on the sense strand. The previously characterized noncoding RNA that regulates expression of the key seed dormancy regulator, DELAY OF GERMINATION1, is a typical representative of this class of RNAs. Transcription factor genes are overrepresented among loci with exosome-sensitive antisense RNAs, suggesting a potential for widespread control of gene expression via this class of noncoding RNAs. Lastly, we assess the use of alternative promoters in Arabidopsis and compare the accuracy of existing TSS annotations.

CAGEfightR: analysis of 5'-end data using R/Bioconductor

BACKGROUND

5'-end sequencing assays, and Cap Analysis of Gene Expression (CAGE) in particular, have been instrumental in studying transcriptional regulation. 5'-end methods provide genome-wide maps of transcription start sites (TSSs) with base pair resolution. Because active enhancers often feature bidirectional TSSs, such data can also be used to predict enhancer candidates. The current availability of mature and comprehensive computational tools for the analysis of 5'-end data is limited, preventing efficient analysis of new and existing 5'-end data.

RESULTS

We present CAGEfightR, a framework for analysis of CAGE and other 5'-end data implemented as an R/Bioconductor-package. CAGEfightR can import data from BigWig files and allows for fast and memory efficient prediction and analysis of TSSs and enhancers. Downstream analyses include quantification, normalization, annotation with transcript and gene models, TSS shape statistics, linking TSSs to enhancers via co-expression, identification of enhancer clusters, and genome-browser style visualization. While built to analyze CAGE data, we demonstrate the utility of CAGEfightR in analyzing nascent RNA 5'-data (PRO-Cap). CAGEfightR is implemented using standard Bioconductor classes, making it easy to learn, use and combine with other Bioconductor packages, for example popular differential expression tools such as limma, DESeq2 and edgeR.

CONCLUSIONS

CAGEfightR provides a single, scalable and easy-to-use framework for comprehensive downstream analysis of 5'-end data. CAGEfightR is designed to be interoperable with other Bioconductor packages, thereby unlocking hundreds of mature transcriptomic analysis tools for 5'-end data. CAGEfightR is freely available via Bioconductor: bioconductor.org/packages/CAGEfightR .