David G. Hendrickson

Differential analysis of gene regulation at transcript resolution with RNA-seq

Differential analysis of gene and transcript expression using high-throughput RNA sequencing (RNA-seq) is complicated by several sources of measurement variability and poses numerous statistical challenges. We present Cuffdiff 2, an algorithm that estimates expression at transcript-level resolution and controls for variability evident across replicate libraries. Cuffdiff 2 robustly identifies differentially expressed transcripts and genes and reveals differential splicing and promoter-preference changes. We demonstrate the accuracy of our approach through differential analysis of lung fibroblasts in response to loss of the developmental transcription factor HOXA1, which we show is required for lung fibroblast and HeLa cell cycle progression. Loss of HOXA1 results in significant expression level changes in thousands of individual transcripts, along with isoform switching events in key regulators of the cell cycle. Cuffdiff 2 performs robust differential analysis in RNA-seq experiments at transcript resolution, revealing a layer of regulation not readily observable with other high-throughput technologies.

Concordant regulation of translation and mRNA abundance for hundreds of targets of a human microRNA.

A specific microRNA reduces the synthesis of hundreds of proteins via concordant effects on the abundance and translation of the mRNAs that encode them.

Topological organization of multichromosomal regions by the long intergenic noncoding RNA Firre

RNA is known to be an abundant and important structural component of the nuclear matrix, including long noncoding RNAs (lncRNA). Yet the molecular identities, functional roles, and localization dynamics of lncRNAs that influence nuclear architecture remain poorly understood. Here, we describe one lncRNA, Firre, that interacts with the nuclear matrix factor hnRNPU, through a 156 bp repeating sequence and Firre localizes across a ~5 Mb domain on the X-chromosome. We further observed Firre localization across at least five distinct trans-chromosomal loci, which reside in spatial proximity to the Firre genomic locus on the X-chromosome. Both genetic deletion of the Firre locus or knockdown of hnRNPU resulted in loss of co-localization of these trans-chromosomal interacting loci. Thus, our data suggest a model in which lncRNAs such as Firre can interface with and modulate nuclear architecture across chromosomes.RNA, including long noncoding RNA (lncRNA), is known to be an abundant and important structural component of the nuclear matrix. However, the molecular identities, functional roles and localization dynamics of lncRNAs that influence nuclear architecture remain poorly understood. Here, we describe one lncRNA, Firre, that interacts with the nuclear-matrix factor hnRNPU through a 156-bp repeating sequence and localizes across an ~5-Mb domain on the X chromosome. We further observed Firre localization across five distinct trans-chromosomal loci, which reside in spatial proximity to the Firre genomic locus on the X chromosome. Both genetic deletion of the Firre locus and knockdown of hnRNPU resulted in loss of colocalization of these trans-chromosomal interacting loci. Thus, our data suggest a model in which lncRNAs such as Firre can interface with and modulate nuclear architecture across chromosomes.

Long noncoding RNAs regulate adipogenesis.

The prevalence of obesity has led to a surge of interest in understanding the detailed mechanisms underlying adipocyte development. Many protein-coding genes, mRNAs, and microRNAs have been implicated in adipocyte development, but the global expression patterns and functional contributions of long noncoding RNA (lncRNA) during adipogenesis have not been explored. Here we profiled the transcriptome of primary brown and white adipocytes, preadipocytes, and cultured adipocytes and identified 175 lncRNAs that are specifically regulated during adipogenesis. Many lncRNAs are adipose-enriched, strongly induced during adipogenesis, and bound at their promoters by key transcription factors such as peroxisome proliferator-activated receptor γ (PPARγ) and CCAAT/enhancer-binding protein α (CEBPα). RNAi-mediated loss of function screens identified functional lncRNAs with varying impact on adipogenesis. Collectively, we have identified numerous lncRNAs that are functionally required for proper adipogenesis.

Systematic Identification of mRNAs Recruited to Argonaute 2 by Specific microRNAs and Corresponding Changes in Transcript Abundance

microRNAs (miRNAs) are small non-coding RNAs that regulate mRNA stability and translation through the action of the RNAi-induced silencing complex (RISC). Our current understanding of miRNA function is inferred largely from studies of the effects of miRNAs on steady-state mRNA levels and from seed match conservation and context in putative targets. Here we have taken a more direct approach to these issues by comprehensively assessing the miRNAs and mRNAs that are physically associated with Argonaute 2 (Ago2), which is a core RISC component. We transfected HEK293T cells with epitope-tagged Ago2, immunopurified Ago2 together with any associated miRNAs and mRNAs, and quantitatively determined the levels of these RNAs by microarray analyses. We found that Ago2 immunopurified samples contained a representative repertoire of the cells miRNAs and a select subset of the cells total mRNAs. Transfection of the miRNAs miR-1 and miR-124 caused significant changes in the association of scores of mRNAs with Ago2. The mRNAs whose association with Ago2 increased upon miRNA expression were much more likely to contain specific miRNA seed matches and to have their overall mRNA levels decrease in response to the miRNA transfection than expected by chance. Hundreds of mRNAs were recruited to Ago2 by each miRNA via seed sequences in 3′-untranslated regions and coding sequences and a few mRNAs appear to be targeted via seed sequences in 5′-untranslated regions. Microarray analysis of Ago2 immunopurified samples provides a simple, direct method for experimentally identifying the targets of miRNAs and for elucidating roles of miRNAs in cellular regulation.

Programming human pluripotent stem cells into white and brown adipocytes

The utility of human pluripotent stem cells is dependent on efficient differentiation protocols that convert these cells into relevant adult cell types. Here we report the robust and efficient differentiation of human pluripotent stem cells into white or brown adipocytes. We found that inducible expression of PPARG2 alone or combined with CEBPB and/or PRDM16 in mesenchymal progenitor cells derived from pluripotent stem cells programmed their development towards a white or brown adipocyte cell fate with efficiencies of 85%–90%. These adipocytes retained their identity independent of transgene expression, could be maintained in culture for several weeks, expressed mature markers and had mature functional properties such as lipid catabolism and insulin-responsiveness. When transplanted into mice, the programmed cells gave rise to ectopic fat pads with the morphological and functional characteristics of white or brown adipose tissue. These results indicate that the cells could be used to faithfully model human disease.

miR-142 regulates the tumorigenicity of human breast cancer stem cells through the canonical WNT signaling pathway

MicroRNAs (miRNAs) are important regulators of stem and progenitor cell functions. We previously reported that miR-142 and miR-150 are upregulated in human breast cancer stem cells (BCSCs) as compared to the non-tumorigenic breast cancer cells. In this study, we report that miR-142 efficiently recruits the APC mRNA to an RNA-induced silencing complex, activates the canonical WNT signaling pathway in an APC-suppression dependent manner, and activates the expression of miR-150. Enforced expression of miR-142 or miR-150 in normal mouse mammary stem cells resulted in the regeneration of hyperproliferative mammary glands in vivo. Knockdown of endogenous miR-142 effectively suppressed organoid formation by BCSCs and slowed tumor growth initiated by human BCSCs in vivo. These results suggest that in some tumors, miR-142 regulates the properties of BCSCs at least in part by activating the WNT signaling pathway and miR-150 expression. DOI: http://dx.doi.org/10.7554/eLife.01977.001

Genome-wide RNA-Seq of Human Motor Neurons Implicates Selective ER Stress Activation in Spinal Muscular Atrophy.

Spinal muscular atrophy (SMA) is caused by mutations in the SMN1 gene. Because this gene is expressed ubiquitously, it remains poorly understood why motor neurons (MNs) are one of the most affected cell types. To address this question, we carried out RNA sequencing studies using fixed, antibody-labeled, and purified MNs produced from control and SMA patient-derived induced pluripotent stem cells (iPSCs). We found SMA-specific changes in MNs, including hyper-activation of the ER stress pathway. Functional studies demonstrated that inhibition of ER stress improves MN survival in vitro even in MNs expressing low SMN. In SMA mice, systemic delivery of an ER stress inhibitor that crosses the blood-brain barrier led to the preservation of spinal cord MNs. Therefore, our study implies that selective activation of ER stress underlies MN death in SMA. Moreover, the approach we have taken would be broadly applicable to the study of disease-prone human cells in heterogeneous cultures.

Transposable elements modulate human RNA abundance and splicing via specific RNA-protein interactions

BackgroundTransposable elements (TEs) have significantly influenced the evolution of transcriptional regulatory networks in the human genome. Post-transcriptional regulation of human genes by TE-derived sequences has been observed in specific contexts, but has yet to be systematically and comprehensively investigated. Here, we study a collection of 75 CLIP-Seq experiments mapping the RNA binding sites for a diverse set of 51 human proteins to explore the role of TEs in post-transcriptional regulation of human mRNAs and lncRNAs via RNA-protein interactions.ResultsWe detect widespread interactions between RNA binding proteins (RBPs) and many families of TE-derived sequence in the CLIP-Seq data. Further, alignment coverage peaks on specific positions of the TE consensus sequences, illuminating a diversity of TE-specific RBP binding motifs. Evidence of binding and conservation of these motifs in the nonrepetitive transcriptome suggests that TEs have generally appropriated existing sequence preferences of the RBPs. Depletion assays for numerous RBPs show that TE-derived binding sites affect transcript abundance and splicing similarly to nonrepetitive sites. However, in a few cases the effect of RBP binding depends on the specific TE family bound; for example, the ubiquitously expressed RBP HuR confers transcript stability unless bound to an Alu element.ConclusionsOur meta-analysis suggests a widespread role for TEs in shaping RNA-protein regulatory networks in the human genome.

LINC00520 is induced by Src, STAT3, and PI3K and plays a functional role in breast cancer

Long non-coding RNAs (lncRNAs) have been implicated in normal cellular homeostasis as well as pathophysiological conditions, including cancer. Here we performed global gene expression profiling of mammary epithelial cells transformed by oncogenic v-Src, and identified a large subset of uncharacterized lncRNAs potentially involved in breast cancer development. Specifically, our analysis revealed a novel lncRNA, LINC00520 that is upregulated upon ectopic expression of oncogenic v-Src, in a manner that is dependent on the transcription factor STAT3. Similarly, LINC00520 is also increased in mammary epithelial cells transformed by oncogenic PI3K and its expression is decreased upon knockdown of mutant PIK3CA. Additional expression profiling highlight that LINC00520 is elevated in a subset of human breast carcinomas, with preferential enrichment in the basal-like molecular subtype. ShRNA-mediated depletion of LINC00520 results in decreased cell migration and loss of invasive structures in 3D. RNA sequencing analysis uncovers several genes that are differentially expressed upon ectopic expression of LINC00520, a significant subset of which are also induced in v-Src-transformed MCF10A cells. Together, these findings characterize LINC00520 as a lncRNA that is regulated by oncogenic Src, PIK3CA and STAT3, and which may contribute to the molecular etiology of breast cancer.