Increasing evidence suggests that interactions between regulatory genomic elements play an important role in regulating gene expression. We generated a genome-wide interaction map of regulatory elements in human cells (ENCODE tier 1 cells, K562, GM12878) using Chromatin Interaction Analysis by Paired-End Tag sequencing (ChIA-PET) experiments targeting six broadly distributed factors. Bound regions covered 80% of DNase I hypersensitive sites including 99.7% of TSS and 98% of enhancers. Correlating this map with ChIP-seq and RNA-seq data sets revealed cohesin, CTCF, and ZNF143 as key components of three-dimensional chromatin structure and revealed how the distal chromatin state affects gene transcription. Comparison of interactions between cell types revealed that enhancer-promoter interactions were highly cell-type-specific. Construction and comparison of distal and proximal regulatory networks revealed stark differences in structure and biological function. Proximal binding events are enriched at genes with housekeeping functions, while distal binding events interact with genes involved in dynamic biological processes including response to stimulus. This study reveals new mechanistic and functional insights into regulatory region organization in the nucleus.

Alternative splicing shapes mammalian transcriptomes, with many RNA molecules undergoing multiple distant alternative splicing events. Comprehensive transcriptome analysis, including analysis of exon co-association in the same molecule, requires deep, long-read sequencing. Here we introduce an RNA sequencing method, synthetic long-read RNA sequencing (SLR-RNA-seq), in which small pools (≤1,000 molecules/pool, ≤1 molecule/gene for most genes) of full-length cDNAs are amplified, fragmented and short-read-sequenced. We demonstrate that these RNA sequences reconstructed from the short reads from each of the pools are mostly close to full length and contain few insertion and deletion errors. We report many previously undescribed isoforms (human brain: ∼13,800 affected genes, 14.5% of molecules; mouse brain ∼8,600 genes, 18% of molecules) and up to 165 human distant molecularly associated exon pairs (dMAPs) and distant molecularly and mutually exclusive pairs (dMEPs). Of 16 associated pairs detected in the mouse brain, 9 are conserved in human. Our results indicate conserved mechanisms that can produce distant but phased features on transcript and proteome isoforms.

Author contributions L.X. and A.M contributed equally to this study. S.M., E.M. and A.B. planned the study, with input from M.P.S. and J.W. S.M. and E.M. performed the reprogramming experiments, and analysed and interpreted data. S.M., E.M. and L.X. wrote the manuscript with the help of A.B. S.M. generated, processed and analysed bulk and single-cell RNA-seq datasets, analysed the metabolomics data, and performed most THY1-related and conditioned medium experiments. E.M. generated and propagated transgene-free iPS cell lines. All other studies were done by both E.M and S.M., unless otherwise noted. A.M. and S.M. performed wound healing experiments under the supervision of M.T.L. L.X. helped with reprogramming, FACS, and immunofluorescence experiments. F.J. generated the in vitro single-cell RNA-seq data under the supervision of M.P.S. R.S. generated the ChIP-seq libraries under the supervision of J.W. K.H. helped with statistics and PAGODA analysis. X.L. performed metabolomics experiments and helped with metabolomics data analysis and validation under the supervision of M.P.S. K.D. helped with reprogramming and western blotting experiments. L.P. helped with reprogramming and RT-qPCR experiments. C.E.A. and Y.S. performed the induced neuron reprogramming experiment under the supervision of M.W. B.A.B. helped with analysis of the epigenomic data. A.L.S.C. identified and collected the human samples. All authors discussed the results and commented on the manuscript. Data availabilityAll raw sequencing reads for population RNA-seq, ChIP-seq and single-cell RNA-seq data can be found under BioProject PRJNA316110. The command and configuration files, in addition to a list of all versioned dependencies present in the running environment, are available on the Github repository for this paper (https://github.com/brunetlab/Mahmoudi_et_al_2018) (except for the code for the processing of metabolomics data, which is available upon request).

Left ventricular non-compaction (LVNC) is the third most prevalent cardiomyopathy in children and its pathogenesis has been associated with the developmental defect of the embryonic myocardium. We show that patient-specific induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs) generated from LVNC patients carrying a mutation in the cardiac transcription factor TBX20 recapitulate a key aspect of the pathological phenotype at the single-cell level and was associated with perturbed transforming growth factor beta (TGFβ) signaling. LVNC iPSC-CMs have decreased proliferative capacity due to abnormal activation of TGFβ signaling. TBX20 regulates the expression of TGFβ signaling modifiers including a known genetic cause of LVNC, PRDM16, and genome editing of PRDM16 caused proliferation defects in iPSC-CMs. Inhibition of TGFβ signaling and genome correction of the TBX20 mutation were sufficient to reverse the disease phenotype. Our study demonstrates that iPSC-CMs are a useful tool for the exploration of pathological mechanisms underlying poorly understood cardiomyopathies including LVNC.

The ubiquitous presence of long noncoding RNAs (lncRNAs) in eukaryotes points to the importance of understanding how their sequences impact function. As many lncRNAs regulate nuclear events and thus must localize to nuclei, we analyzed the sequence requirements for nuclear localization in an intergenic lncRNA named BORG (BMP2-OP1-responsive gene), which is both spliced and polyadenylated but is strictly localized in nuclei. Subcellular localization of BORG was not dependent on the context or level of its expression or decay but rather depended on the sequence of the mature, spliced transcript. Mutational analyses indicated that nuclear localization of BORG was mediated through a novel RNA motif consisting of the pentamer sequence AGCCC with sequence restrictions at positions ؊8 (T or A) and ؊3 (G or C) relative to the first nucleotide of the pentamer. Mutation of the motif to a scrambled sequence resulted in complete loss of nuclear localization, while addition of even a single copy of the motif to a cytoplasmically localized RNA was sufficient to impart nuclear localization. Further, the presence of this motif in other cellular RNAs showed a direct correlation with nuclear localization, suggesting that the motif may act as a general nuclear localization signal for cellular RNAs.