Circadian-driven tissue specificity is constrained under caloric restricted feeding conditions.

Tissue specificity is a fundamental property of an organ that affects numerous biological processes, including aging and longevity, and is regulated by the circadian clock. However, the distinction between circadian-affected tissue specificity and other tissue specificities remains poorly understood. Here, using multi-omics data on circadian rhythms in mice, we discovered that approximately 35% of tissue-specific genes are directly affected by circadian regulation. These circadian-affected tissue-specific genes have higher expression levels and are associated with metabolism in hepatocytes. They also exhibit specific features in long-reads sequencing data. Notably, these genes are associated with aging and longevity at both the gene level and at the network module level. The expression of these genes oscillates in response to caloric restricted feeding regimens, which have been demonstrated to promote longevity. In addition, aging and longevity genes are disrupted in various circadian disorders. Our study indicates that the modulation of circadian-affected tissue specificity is essential for understanding the circadian mechanisms that regulate aging and longevity at the genomic level.

Chemical-induced phase transition and global conformational reorganization of chromatin.

Chemicals or drugs can accumulate within biomolecular condensates formed through phase separation in cells. Here, we use super-resolution imaging to search for chemicals that induce phase transition within chromatin at the microscale. This microscopic screening approach reveals that adriamycin (doxorubicin) - a widely used anticancer drug that is known to interact with chromatin - specifically induces visible local condensation and global conformational change of chromatin in cancer and primary cells. Hi-C and ATAC-seq experiments systematically and quantitatively demonstrate that adriamycin-induced chromatin condensation is accompanied by weakened chromatin interaction within topologically associated domains, compartment A/B switching, lower chromatin accessibility, and corresponding transcriptomic changes. Mechanistically, adriamycin complexes with histone H1 and induces phase transition of H1, forming fibrous aggregates in vitro. These results reveal a phase separation-driven mechanism for a chemotherapeutic drug.

Noncoding transcripts are linked to brain resting-state activity in non-human primates.

Brain-derived transcriptomes are known to correlate with resting-state brain activity in humans. Whether this association holds in nonhuman primates remains uncertain. Here, we search for such molecular correlates by integrating 757 transcriptomes derived from 100 macaque cortical regions with resting-state activity in separate conspecifics. We observe that 150 noncoding genes explain variations in resting-state activity at a comparable level with protein-coding genes. In-depth analysis of these noncoding genes reveals that they are connected to the function of nonneuronal cells such as oligodendrocytes. Co-expression network analysis finds that the modules of noncoding genes are linked to both autism and schizophrenia risk genes. Moreover, genes associated with resting-state noncoding genes are highly enriched in human resting-state functional genes and memory-effect genes, and their links with resting-state functional magnetic resonance imaging (fMRI) signals are altered in the brains of patients with autism. Our results highlight the potential for noncoding RNAs to explain resting-state activity in the nonhuman primate brain.

Brain-wide and cell-specific transcriptomic insights into MRI-derived cortical morphology in macaque monkeys.

Integrative analyses of transcriptomic and neuroimaging data have generated a wealth of information about biological pathways underlying regional variability in imaging-derived brain phenotypes in humans, but rarely in nonhuman primates due to the lack of a comprehensive anatomically-defined atlas of brain transcriptomics. Here we generate complementary bulk RNA-sequencing dataset of 819 samples from 110 brain regions and single-nucleus RNA-sequencing dataset, and neuroimaging data from 162 cynomolgus macaques, to examine the link between brain-wide gene expression and regional variation in morphometry. We not only observe global/regional expression profiles of macaque brain comparable to human but unravel a dorsolateral-ventromedial gradient of gene assemblies within the primate frontal lobe. Furthermore, we identify a set of 971 protein-coding and 34 non-coding genes consistently associated with cortical thickness, specially enriched for neurons and oligodendrocytes. These data provide a unique resource to investigate nonhuman primate models of human diseases and probe cross-species evolutionary mechanisms.

Significant Evolutionary Constraints on Neuron Cells Revealed by Single-Cell Transcriptomics.

Recent advances in single-cell RNA-sequencing technology have enabled us to characterize a variety of different cell types in each brain region. However, the evolutionary differences among these cell types remain unclear. Here, we analyzed single-cell RNA-seq data of >280,000 cells and developmental transcriptomes of bulk brain tissues. At the single-cell level, we found that the evolutionary constraints on the cell types of different organs significantly overlap with each other and the transcriptome of neuron cells is one of the most restricted evolutionarily. In addition, mature neurons are under more constraints than neuron stem cells as well as nascent neurons and the order of the constraints of various cell types of the brain is largely conserved in different subregions. We also found that although functionally similar brain regions have comparable evolutionary constraints, the early fetal brain is the least constrained and this pattern is conserved in the mouse, macaque, and humans. These results demonstrate the importance of maintaining the plasticity of early brain development during evolution. The delineation of evolutionary differences between brain cell types has great potential for an improved understanding of the pathogenesis of neurological diseases and drug development efforts aimed at the manipulation of molecular activities at the single-cell level.

Single-cell RNA-seq reveals distinct injury responses in different types of DRG sensory neurons.

Peripheral nerve injury leads to various injury-induced responses in sensory neurons including physiological pain, neuronal cell death, and nerve regeneration. In this study, we performed single-cell RNA-sequencing (scRNA-seq) analysis of mouse nonpeptidergic nociceptors (NP), peptidergic nociceptors (PEP), and large myelinated sensory neurons (LM) under both control and injury conditions at 3 days after sciatic nerve transection (SNT). After performing principle component and weighted gene co-expression network analysis, we categorized dorsal root ganglion (DRG) neurons into different subtypes and discovered co-regulated injury-response genes including novel regeneration associated genes (RAGs) in association with neuronal development, protein translation and cytoplasm transportation. In addition, we found significant up-regulation of the genes associated with cell death such as Pdcd2 in a subset of NP neurons after axotomy, implicating their actions in neuronal cell death upon nerve injury. Our study revealed the distinctive and sustained heterogeneity of transcriptomic responses to injury at single neuron level, implicating the involvement of different gene regulatory networks in nerve regeneration, neuronal cell death and neuropathy in different population of DRG neurons.

Simultaneous profiling of transcriptome and DNA methylome from a single cell.

BACKGROUND: Single-cell transcriptome and single-cell methylome technologies have become powerful tools to study RNA and DNA methylation profiles of single cells at a genome-wide scale. A major challenge has been to understand the direct correlation of DNA methylation and gene expression within single-cells. Due to large cell-to-cell variability and the lack of direct measurements of transcriptome and methylome of the same cell, the association is still unclear. RESULTS: Here, we describe a novel method (scMT-seq) that simultaneously profiles both DNA methylome and transcriptome from the same cell. In sensory neurons, we consistently identify transcriptome and methylome heterogeneity among single cells but the majority of the expression variance is not explained by proximal promoter methylation, with the exception of genes that do not contain CpG islands. By contrast, gene body methylation is positively associated with gene expression for only those genes that contain a CpG island promoter. Furthermore, using single nucleotide polymorphism patterns from our hybrid mouse model, we also find positive correlation of allelic gene body methylation with allelic expression. CONCLUSIONS: Our method can be used to detect transcriptome, methylome, and single nucleotide polymorphism information within single cells to dissect the mechanisms of epigenetic gene regulation.

Selective demethylation and altered gene expression are associated with ICF syndrome in human-induced pluripotent stem cells and mesenchymal stem cells.

Immunodeficiency, centromeric instability and facial anomalies type I (ICF1) syndrome is a rare genetic disease caused by mutations in DNA methyltransferase (DNMT) 3B, a de novo DNA methyltransferase. However, the molecular basis of how DNMT3B deficiency leads to ICF1 pathogenesis is unclear. Induced pluripotent stem cell (iPSC) technology facilitates the study of early human developmental diseases via facile in vitro paradigms. Here, we generate iPSCs from ICF Type 1 syndrome patient fibroblasts followed by directed differentiation of ICF1-iPSCs to mesenchymal stem cells (MSCs). By performing genome-scale bisulfite sequencing, we find that DNMT3B-deficient iPSCs exhibit global loss of non-CG methylation and select CG hypomethylation at gene promoters and enhancers. Further unbiased scanning of ICF1-iPSC methylomes also identifies large megabase regions of CG hypomethylation typically localized in centromeric and subtelomeric regions. RNA sequencing of ICF1 and control iPSCs reveals abnormal gene expression in ICF1-iPSCs relevant to ICF syndrome phenotypes, some directly associated with promoter or enhancer hypomethylation. Upon differentiation of ICF1 iPSCs to MSCs, we find virtually all CG hypomethylated regions remained hypomethylated when compared with either wild-type iPSC-derived MSCs or primary bone-marrow MSCs. Collectively, our results show specific methylome and transcriptome defects in both ICF1-iPSCs and differentiated somatic cell lineages, providing a valuable stem cell system for further in vitro study of the molecular pathogenesis of ICF1 syndrome. GEO accession number: GSE46030.

Integrated analysis of DNA methylation and RNA transcriptome during in vitro differentiation of human pluripotent stem cells into retinal pigment epithelial cells.

Using the paradigm of in vitro differentiation of hESCs/iPSCs into retinal pigment epithelial (RPE) cells, we have recently profiled mRNA and miRNA transcriptomes to define a set of RPE mRNA and miRNA signature genes implicated in directed RPE differentiation. In this study, in order to understand the role of DNA methylation in RPE differentiation, we profiled genome-scale DNA methylation patterns using the method of reduced representation bisulfite sequencing (RRBS). We found dynamic waves of de novo methylation and demethylation in four stages of RPE differentiation. Integrated analysis of DNA methylation and RPE transcriptomes revealed a reverse-correlation between levels of DNA methylation and expression of a subset of miRNA and mRNA genes that are important for RPE differentiation and function. Gene Ontology (GO) analysis suggested that genes undergoing dynamic methylation changes were related to RPE differentiation and maturation. We further compared methylation patterns among human ESC- and iPSC-derived RPE as well as primary fetal RPE (fRPE) cells, and discovered that specific DNA methylation pattern is useful to classify each of the three types of RPE cells. Our results demonstrate that DNA methylation may serve as biomarkers to characterize the cell differentiation process during the conversion of human pluripotent stem cells into functional RPE cells.

Identification of miRNA signatures during the differentiation of hESCs into retinal pigment epithelial cells.

Retinal pigment epithelium (RPE) cells can be obtained through in vitro differentiation of both embryonic stem cells (ESCs) and induced pluripotent stem cells (iPSCs). We have previously identified 87 signature genes relevant to RPE cell differentiation and function through transcriptome analysis of both human ESC- and iPSC-derived RPE as well as normal fetal RPE. Here, we profile miRNA expression through small RNA-seq in human ESCs and their RPE derivatives. Much like conclusions drawn from our previous transcriptome analysis, we find that the overall miRNA landscape in RPE is distinct from ESCs and other differentiated somatic tissues. We also profile miRNA expression during intermediate stages of RPE differentiation and identified unique subsets of miRNAs that are gradually up- or down-regulated, suggesting that dynamic regulation of these miRNAs is associated with the RPE differentiation process. Indeed, the down-regulation of a subset of miRNAs during RPE differentiation is associated with up-regulation of RPE-specific genes, such as RPE65, which is exclusively expressed in RPE. We conclude that miRNA signatures can be used to classify different degrees of in vitro differentiation of RPE from human pluripotent stem cells. We suggest that RPE-specific miRNAs likely contribute to the functional maturation of RPE in vitro, similar to the regulation of RPE-specific mRNA expression.