Single-cell multiplex chromatin and RNA interactions in ageing human brain.

Dynamically organized chromatin complexes often involve multiplex chromatin interactions and sometimes chromatin-associated RNA1-3. Chromatin complex compositions change during cellular differentiation and ageing, and are expected to be highly heterogeneous among terminally differentiated single cells4-7. Here we introduce the multinucleic acid interaction mapping in single cells (MUSIC) technique for concurrent profiling of multiplex chromatin interactions, gene expression and RNA-chromatin associations within individual nuclei. When applied to 14 human frontal cortex samples from older donors, MUSIC delineated diverse cortical cell types and states. We observed that nuclei exhibiting fewer short-range chromatin interactions were correlated with both an 'older' transcriptomic signature and Alzheimer's disease pathology. Furthermore, the cell type exhibiting chromatin contacts between cis expression quantitative trait loci and a promoter tends to be that in which these cis expression quantitative trait loci specifically affect the expression of their target gene. In addition, female cortical cells exhibit highly heterogeneous interactions between XIST non-coding RNA and chromosome X, along with diverse spatial organizations of the X chromosomes. MUSIC presents a potent tool for exploration of chromatin architecture and transcription at cellular resolution in complex tissues.

Regulation of nuclear transcription by mitochondrial RNA in endothelial cells.

Chromatin-associated RNAs (caRNAs) form a relatively poorly recognized layer of the epigenome. The caRNAs reported to date are transcribed from the nuclear genome. Here, leveraging a recently developed assay for detection of caRNAs and their genomic association, we report that mitochondrial RNAs (mtRNAs) are attached to the nuclear genome and constitute a subset of caRNA, thus termed mt-caRNA. In four human cell types analyzed, mt-caRNAs preferentially attach to promoter regions. In human endothelial cells (ECs), the level of mt-caRNA-promoter attachment changes in response to environmental stress that mimics diabetes. Suppression of a non-coding mt-caRNA in ECs attenuates stress-induced nascent RNA transcription from the nuclear genome, including that of critical genes regulating cell adhesion, and abolishes stress-induced monocyte adhesion, a hallmark of dysfunctional ECs. Finally, we report increased nuclear localization of multiple mtRNAs in the ECs of human diabetic donors, suggesting many mtRNA translocate to the nucleus in a cell stress and disease-dependent manner. These data nominate mt-caRNAs as messenger molecules responsible for mitochondrial-nuclear communication and connect the immediate product of mitochondrial transcription with the transcriptional regulation of the nuclear genome.

Single-cell multiplex chromatin and RNA interactions in aging human brain.

The dynamically organized chromatin complexes often involve multiplex chromatin interactions and sometimes chromatin-associated RNA (caRNA) 1-3. Chromatin complex compositions change during cellular differentiation and aging, and are expected to be highly heterogeneous among terminally differentiated single cells 4-7. Here we introduce the Multi-Nucleic Acid Interaction Mapping in Single Cell (MUSIC) technique for concurrent profiling of multiplex chromatin interactions, gene expression, and RNA-chromatin associations within individual nuclei. Applied to 14 human frontal cortex samples from elderly donors, MUSIC delineates diverse cortical cell types and states. We observed the nuclei exhibiting fewer short-range chromatin interactions are correlated with an "older" transcriptomic signature and with Alzheimer's pathology. Furthermore, the cell type exhibiting chromatin contacts between cis expression quantitative trait loci (cis eQTLs) and a promoter tends to be the cell type where these cis eQTLs specifically affect their target gene's expression. Additionally, the female cortical cells exhibit highly heterogeneous interactions between the XIST non-coding RNA and Chromosome X, along with diverse spatial organizations of the X chromosomes. MUSIC presents a potent tool for exploring chromatin architecture and transcription at cellular resolution in complex tissues.

Genome-wide analysis of the interplay between chromatin-associated RNA and 3D genome organization in human cells.

The interphase genome is dynamically organized in the nucleus and decorated with chromatin-associated RNA (caRNA). It remains unclear whether the genome architecture modulates the spatial distribution of caRNA and vice versa. Here, we generate a resource of genome-wide RNA-DNA and DNA-DNA contact maps in human cells. These maps reveal the chromosomal domains demarcated by locally transcribed RNA, hereafter termed RNA-defined chromosomal domains. Further, the spreading of caRNA is constrained by the boundaries of topologically associating domains (TADs), demonstrating the role of the 3D genome structure in modulating the spatial distribution of RNA. Conversely, stopping transcription or acute depletion of RNA induces thousands of chromatin loops genome-wide. Activation or suppression of the transcription of specific genes suppresses or creates chromatin loops straddling these genes. Deletion of a specific caRNA-producing genomic sequence promotes chromatin loops that straddle the interchromosomal target sequences of this caRNA. These data suggest a feedback loop where the 3D genome modulates the spatial distribution of RNA, which in turn affects the dynamic 3D genome organization.

High-Resolution Characterization of Human Brain Cortex with High-Fidelity Spatial Transcriptomic Slides (HiFi-Slides).

Spatial transcriptomic tools and platforms help researchers to inspect tissues and cells with fine details of how they differentiate in expressions and how they orient themselves. With the higher resolution we get and higher throughput of expression targets, spatial analysis can truly become the core player for cell clustering, migration study, and, eventually, the novel model for pathological study. We present the demonstration of HiFi-slide, a whole transcriptomic sequencing technique that recycles used sequenced-by-synthesis flow cell surfaces to a high-resolution spatial mapping tool that can be directly applied to tissue cell gradient analysis, gene expression analysis, cell proximity analysis, and other cellular-level spatial studies.

Spatial and temporal organization of the genome: Current state and future aims of the 4D nucleome project.

The four-dimensional nucleome (4DN) consortium studies the architecture of the genome and the nucleus in space and time. We summarize progress by the consortium and highlight the development of technologies for (1) mapping genome folding and identifying roles of nuclear components and bodies, proteins, and RNA, (2) characterizing nuclear organization with time or single-cell resolution, and (3) imaging of nuclear organization. With these tools, the consortium has provided over 2,000 public datasets. Integrative computational models based on these data are starting to reveal connections between genome structure and function. We then present a forward-looking perspective and outline current aims to (1) delineate dynamics of nuclear architecture at different timescales, from minutes to weeks as cells differentiate, in populations and in single cells, (2) characterize cis-determinants and trans-modulators of genome organization, (3) test functional consequences of changes in cis- and trans-regulators, and (4) develop predictive models of genome structure and function.

Genetic Deletion of the LINC00520 Homolog in Mouse Aggravates Angiotensin II-Induced Hypertension.

(1) Background: Hypertension is a complex, multifactorial disease that is caused by genetic and environmental factors. Apart from genetic predisposition, the mechanisms involved in this disease have yet to be fully understood. We previously reported that LEENE (lncRNA enhancing endothelial nitric oxide expression, transcribed from LINC00520 in the human genome) regulates endothelial cell (EC) function by promoting the expression of endothelial nitric oxide synthase (eNOS) and vascular growth factor receptor 2 (VEGFR2). Mice with genetic deletion of the LEENE/LINC00520 homologous region exhibited impaired angiogenesis and tissue regeneration in a diabetic hindlimb ischemia model. However, the role of LEENE in blood pressure regulation is unknown. (2) Methods: We subjected mice with genetic ablation of leene and wild-type littermates to Angiotensin II (AngII) and monitored their blood pressure and examined their hearts and kidneys. We used RNA-sequencing to identify potential leene-regulated molecular pathways in ECs that contributed to the observed phenotype. We further performed in vitro experiments with murine and human ECs and ex vivo experiments with murine aortic rings to validate the select mechanism. (3) Results: We identified an exacerbated hypertensive phenotype of leene-KO mice in the AngII model, evidenced by higher systolic and diastolic blood pressure. At the organ level, we observed aggravated hypertrophy and fibrosis in the heart and kidney. Moreover, the overexpression of human LEENE RNA, in part, restored the signaling pathways impaired by leene deletion in murine ECs. Additionally, Axitinib, a tyrosine kinase inhibitor that selectively inhibits VEGFR suppresses LEENE in human ECs. (4) Conclusions: Our study suggests LEENE as a potential regulator in blood pressure control, possibly through its function in ECs.

Long noncoding RNA LEENE promotes angiogenesis and ischemic recovery in diabetes models.

Impaired angiogenesis in diabetes is a key process contributing to ischemic diseases such as peripheral arterial disease. Epigenetic mechanisms, including those mediated by long noncoding RNAs (lncRNAs), are crucial links connecting diabetes and the related chronic tissue ischemia. Here we identify the lncRNA that enhances endothelial nitric oxide synthase (eNOS) expression (LEENE) as a regulator of angiogenesis and ischemic response. LEENE expression was decreased in diabetic conditions in cultured endothelial cells (ECs), mouse hind limb muscles, and human arteries. Inhibition of LEENE in human microvascular ECs reduced their angiogenic capacity with a dysregulated angiogenic gene program. Diabetic mice deficient in Leene demonstrated impaired angiogenesis and perfusion following hind limb ischemia. Importantly, overexpression of human LEENE rescued the impaired ischemic response in Leene-knockout mice at tissue functional and single-cell transcriptomic levels. Mechanistically, LEENE RNA promoted transcription of proangiogenic genes in ECs, such as KDR (encoding VEGFR2) and NOS3 (encoding eNOS), potentially by interacting with LEO1, a key component of the RNA polymerase II-associated factor complex and MYC, a crucial transcription factor for angiogenesis. Taken together, our findings demonstrate an essential role for LEENE in the regulation of angiogenesis and tissue perfusion. Functional enhancement of LEENE to restore angiogenesis for tissue repair and regeneration may represent a potential strategy to tackle ischemic vascular diseases.

PHGDH expression increases with progression of Alzheimer's disease pathology and symptoms.

Chen et al. reveal an increase of phosphoglycerate dehydrogenase (PHGDH) mRNA and protein levels in two mouse models and four human cohorts in Alzheimer's disease brains compared to age- and sex-matched control brains. The increase of PHGDH expression in human brain correlates with symptomatic development and disease pathology.

Isolation and Profiling of Human Primary Mesenteric Arterial Endothelial Cells at the Transcriptome Level.

Endothelial cells (ECs) are crucial for vascular and whole-body function through their dynamic response to environmental cues. Elucidating the transcriptome and epigenome of ECs is paramount to understanding their roles in development, health, and disease, but is limited in the availability of isolated primary cells. Recent technologies have enabled the high-throughput profiling of EC transcriptome and epigenome, leading to the identification of previously unknown EC cell subpopulations and developmental trajectories. While EC cultures are a useful tool in the exploration of EC function and dysfunction, the culture conditions and multiple passages can introduce external variables that alter the properties of native EC, including morphology, epigenetic state, and gene expression program. To overcome this limitation, the present paper demonstrates a method of isolating human primary ECs from donor mesenteric arteries aiming to capture their native state. ECs in the intimal layer are dissociated mechanically and biochemically with the use of particular enzymes. The resultant cells can be directly used for bulk RNA or single-cell RNA-sequencing or plated for culture. In addition, a workflow is described for the preparation of human arterial tissue for spatial transcriptomics, specifically for a commercially available platform, although this method is also suitable for other spatial transcriptome profiling techniques. This methodology can be applied to different vessels collected from a variety of donors in health or disease states to gain insights into EC transcriptional and epigenetic regulation, a pivotal aspect of endothelial cell biology.

Stress-induced RNA-chromatin interactions promote endothelial dysfunction.

Chromatin-associated RNA (caRNA) has been proposed as a type of epigenomic modifier. Here, we test whether environmental stress can induce cellular dysfunction through modulating RNA-chromatin interactions. We induce endothelial cell (EC) dysfunction with high glucose and TNFα (H + T), that mimic the common stress in diabetes mellitus. We characterize the H + T-induced changes in gene expression by single cell (sc)RNA-seq, DNA interactions by Hi-C, and RNA-chromatin interactions by iMARGI. H + T induce inter-chromosomal RNA-chromatin interactions, particularly among the super enhancers. To test the causal relationship between H + T-induced RNA-chromatin interactions and the expression of EC dysfunction-related genes, we suppress the LINC00607 RNA. This suppression attenuates the expression of SERPINE1, a critical pro-inflammatory and pro-fibrotic gene. Furthermore, the changes of the co-expression gene network between diabetic and healthy donor-derived ECs corroborate the H + T-induced RNA-chromatin interactions. Taken together, caRNA-mediated dysregulation of gene expression modulates EC dysfunction, a crucial mechanism underlying numerous diseases.

Extracellular RNA in a single droplet of human serum reflects physiologic and disease states.

Extracellular RNAs (exRNAs) are present in human serum. It remains unclear to what extent these circulating exRNAs may reflect human physiologic and disease states. Here, we developed SILVER-seq (Small Input Liquid Volume Extracellular RNA Sequencing) to efficiently sequence both integral and fragmented exRNAs from a small droplet (5 µL to 7 µL) of liquid biopsy. We calibrated SILVER-seq in reference to other RNA sequencing methods based on milliliters of input serum and quantified droplet-to-droplet and donor-to-donor variations. We carried out SILVER-seq on more than 150 serum droplets from male and female donors ranging from 18 y to 48 y of age. SILVER-seq detected exRNAs from more than a quarter of the human genes, including small RNAs and fragments of mRNAs and long noncoding RNAs (lncRNAs). The detected exRNAs included those derived from genes with tissue (e.g., brain)-specific expression. The exRNA expression levels separated the male and female samples and were correlated with chronological age. Noncancer and breast cancer donors exhibited pronounced differences, whereas donors with or without cancer recurrence exhibited moderate differences in exRNA expression patterns. Even without using differentially expressed exRNAs as features, nearly all cancer and noncancer samples and a large portion of the recurrence and nonrecurrence samples could be correctly classified by exRNA expression values. These data suggest the potential of using exRNAs in a single droplet of serum for liquid biopsy-based diagnostics.

GITAR: An Open Source Tool for Analysis and Visualization of Hi-C Data.

Interactions between chromatin segments play a large role in functional genomic assays and developments in genomic interaction detection methods have shown interacting topological domains within the genome. Among these methods, Hi-C plays a key role. Here, we present the Genome Interaction Tools and Resources (GITAR), a software to perform a comprehensive Hi-C data analysis, including data preprocessing, normalization, and visualization, as well as analysis of topologically-associated domains (TADs). GITAR is composed of two main modules: (1) HiCtool, a Python library to process and visualize Hi-C data, including TAD analysis; and (2) processed data library, a large collection of human and mouse datasets processed using HiCtool. HiCtool leads the user step-by-step through a pipeline, which goes from the raw Hi-C data to the computation, visualization, and optimized storage of intra-chromosomal contact matrices and TAD coordinates. A large collection of standardized processed data allows the users to compare different datasets in a consistent way, while saving time to obtain data for visualization or additional analyses. More importantly, GITAR enables users without any programming or bioinformatic expertise to work with Hi-C data. GITAR is publicly available at http://genomegitar.org as an open-source software.

Rainbow-Seq: Combining Cell Lineage Tracing with Single-Cell RNA Sequencing in Preimplantation Embryos.

We developed the Rainbow-seq technology to trace cell division history and reveal single-cell transcriptomes. With distinct fluorescent protein genes as lineage markers, Rainbow-seq enables each single-cell RNA sequencing (RNA-seq) experiment to simultaneously decode the lineage marker genes and read single-cell transcriptomes. We triggered lineage tracking in each blastomere at the 2-cell stage, observed microscopically inequivalent contributions of the progeny to the two embryonic poles at the blastocyst stage, and analyzed every single cell at either 4- or 8-cell stage with deep paired-end sequencing of full-length transcripts. Although lineage difference was not marked unequivocally at a single-gene level, it became clear when the transcriptome was analyzed as a whole. Moreover, several groups of novel transcript isoforms with embedded repeat sequences exhibited lineage difference, suggesting a possible link between DNA demethylation and cell fate decision. Rainbow-seq bridged a critical gap between division history and single-cell RNA-seq assays.