Accurate RNA quantification at the single-cell level is critical for understanding the dynamics of gene expression and regulation across space and time. Single molecule FISH (smFISH), such as RNAscope, provides spatial and quantitative measurements of individual transcripts, therefore, can be used to explore differential gene expression among a heterogeneous cell population if combined with cell identify information. However, such analysis is not straightforward, and existing image analysis pipelines cannot integrate both RNA transcripts and cellular staining information to automatically output cell type-specific gene expression. We developed an efficient and customizable analysis method, Single-Molecule Automatic RNA Transcription Quantification (SMART-Q), to enable the analysis of gene transcripts in a cell type-specific manner. SMART-Q efficiently infers cell identity information from multiplexed immuno-staining and quantifies cell type-specific transcripts using a 3D Gaussian fitting algorithm. Furthermore, we have optimized SMART-Q for user experiences, such as flexible parameters specification, batch data outputs, and visualization of analysis results. SMART-Q meets the demands for efficient quantification of single-molecule RNA and can be widely used for cell type-specific RNA transcript analysis.

Enhancers are important regulatory elements that can control gene activity across vast genetic distances. However, the underlying nature of this regulation remains obscured because it has been difficult to observe in living cells. Here, we visualize the spatial organization and transcriptional output of the key pluripotency regulator Sox2 and its essential enhancer Sox2 Control Region (SCR) in living embryonic stem cells (ESCs). We find that Sox2 and SCR show no evidence of enhanced spatial proximity and that spatial dynamics of this pair is limited over tens of minutes. Sox2 transcription occurs in short, intermittent bursts in ESCs and, intriguingly, we find this activity demonstrates no association with enhancer proximity, suggesting that direct enhancer-promoter contacts do not drive contemporaneous Sox2 transcription. Our study establishes a framework for interrogation of enhancer function in living cells and supports an unexpected mechanism for enhancer control of Sox2 expression that uncouples transcription from enhancer proximity.

With <2% of the human genome coding for proteins, a major challenge is to interpret the function of the noncoding DNA. Millions of regulatory sequences have been predicted in the human genome through analysis of DNA methylation, chromatin modification, hypersensitivity to nucleases, and transcription factor binding, but few have been shown to regulate transcription in their native contexts. We have developed a high-throughput CRISPR/Cas9-based genome-editing strategy and used it to interrogate 174 candidate regulatory sequences within the 1-Mbp POU5F1 locus in human embryonic stem cells (hESCs). We identified two classical regulatory elements, including a promoter and a proximal enhancer, that are essential for POU5F1 transcription in hESCs. Unexpectedly, we also discovered a new class of enhancers that contribute to POU5F1 transcription in an unusual way: Disruption of such sequences led to a temporary loss of POU5F1 transcription that is fully restored after a few rounds of cell division. These results demonstrate the utility of high-throughput screening for functional characterization of noncoding DNA and reveal a previously unrecognized layer of gene regulation in human cells.

A wide range of human diseases result from haploinsufficiency, where the function of one of the two gene copies is lost. Here, we targeted the remaining functional copy of a haploinsufficient gene using CRISPR-mediated activation (CRISPRa) in Sim1 and Mc4r heterozygous mouse models to rescue their obesity phenotype. Transgenic-based CRISPRa targeting of the Sim1 promoter or its distant hypothalamic enhancer up-regulated its expression from the endogenous functional allele in a tissue-specific manner, rescuing the obesity phenotype in Sim1 heterozygous mice. To evaluate the therapeutic potential of CRISPRa, we injected CRISPRa-recombinant adeno-associated virus into the hypothalamus, which led to reversal of the obesity phenotype in Sim1 and Mc4r haploinsufficient mice. Our results suggest that endogenous gene up-regulation could be a potential strategy to treat altered gene dosage diseases.

Acinar cells play an essential role in the secretory function of exocrine organs. Despite this requirement, how acinar cells are generated during organogenesis is unclear. Using the acini-ductal network of the developing human and murine salivary gland, we demonstrate an unexpected role for SOX2 and parasympathetic nerves in generating the acinar lineage that has broad implications for epithelial morphogenesis. Despite SOX2 being expressed by progenitors that give rise to both acinar and duct cells, genetic ablation of SOX2 results in a failure to establish acini but not ducts. Furthermore, we show that SOX2 targets acinar-specific genes and is essential for the survival of acinar but not ductal cells. Finally, we illustrate an unexpected and novel role for peripheral nerves in the creation of acini throughout development via regulation of SOX2. Thus, SOX2 is a master regulator of the acinar cell lineage essential to the establishment of a functional organ.

Droplet single-cell RNA-sequencing (dscRNA-seq) has enabled rapid, massively parallel profiling of transcriptomes. However, assessing differential expression across multiple individuals has been hampered by inefficient sample processing and technical batch effects. Here we describe a computational tool, demuxlet, that harnesses natural genetic variation to determine the sample identity of each droplet containing a single cell (singlet) and detect droplets containing two cells (doublets). These capabilities enable multiplexed dscRNA-seq experiments in which cells from unrelated individuals are pooled and captured at higher throughput than in standard workflows. Using simulated data, we show that 50 single-nucleotide polymorphisms (SNPs) per cell are sufficient to assign 97% of singlets and identify 92% of doublets in pools of up to 64 individuals. Given genotyping data for each of eight pooled samples, demuxlet correctly recovers the sample identity of >99% of singlets and identifies doublets at rates consistent with previous estimates. We apply demuxlet to assess cell-type-specific changes in gene expression in 8 pooled lupus patient samples treated with interferon (IFN)-β and perform eQTL analysis on 23 pooled samples.

Systemic lupus erythematosus (SLE) is an autoimmune disease in which outcomes vary among different racial groups. We leverage cell-sorted RNA-seq data (CD14+ monocytes, B cells, CD4+ T cells, and NK cells) from 120 SLE patients (63 Asian and 57 White individuals) and apply a four-tier approach including unsupervised clustering, differential expression analyses, gene co-expression analyses, and machine learning to identify SLE subgroups within this multiethnic cohort. K-means clustering on each cell-type resulted in three clusters for CD4 and CD14, and two for B and NK cells. To understand the identified clusters, correlation analysis revealed significant positive associations between the clusters and clinical parameters including disease activity as well as ethnicity. We then explored differentially expressed genes between Asian and White groups for each cell-type. The shared differentially expressed genes across cells were involved in SLE or other autoimmune-related pathways. Co-expression analysis identified similarly regulated genes across samples and grouped these genes into modules. Finally, random forest classification of disease activity in the White and Asian cohorts showed the best classification in CD4+ T cells in White individuals. The results from these analyses will help stratify patients based on their gene expression signatures to enable SLE precision medicine.

Objective

There is an established yet unexplained link between interferon (IFN) and systemic lupus erythematosus (SLE). The expression of sequences derived from transposable elements (TEs) may contribute to SLE phenotypes, specifically production of type I IFNs and generation of autoantibodies.

Methods

We profiled cell-sorted RNA-sequencing data (CD4⁺ T cells, CD14⁺ monocytes, CD19⁺ B cells, and natural killer cells) from peripheral blood mononuclear cells of 120 patients with SLE and quantified TE expression identifying 27,135 TEs. We tested for differential TE expression across 10 SLE phenotypes, including autoantibody production and disease activity.

Results

We found 731 differentially expressed (DE) TEs across all SLE phenotypes that were mostly cell specific and phenotype specific. DE TEs were enriched for specific families and open reading frames of viral genes encoded in TE sequences. Increased expression of DE TEs was associated with genes involved in antiviral activity, such as LY6E, ISG15, and TRIM22, and pathways such as IFN signaling.

Conclusion

These findings suggest that expression of TEs contributes to activation of SLE-related mechanisms in a cell-specific manner, which can impact disease diagnostics and therapeutics.

Current pooled CRISPR screens for cis-regulatory elements (CREs), based on transcriptional output changes, are typically limited to characterizing CREs of only one gene. Here, we describe CRISPRpath, a scalable screening strategy for parallelly characterizing CREs of genes linked to the same biological pathway and converging phenotypes. We demonstrate the ability of CRISPRpath for simultaneously identifying functional enhancers of six genes in the 6-thioguanine–induced DNA mismatch repair pathway using both CRISPR interference (CRISPRi) and CRISPR nuclease (CRISPRn) approaches. Sixty percent of the identified enhancers are known promoters with distinct epigenomic features compared to other active promoters, including increased chromatin accessibility and interactivity. Furthermore, by imposing different levels of selection pressure, CRISPRpath can distinguish enhancers exerting strong impact on gene expression from those exerting weak impact. Our results offer a nuanced view of cis-regulation and demonstrate that CRISPRpath can be leveraged for understanding the complex gene regulatory program beyond transcriptional output at scale.

Lineage-specific epigenomic changes during human corticogenesis have been difficult to study owing to challenges with sample availability and tissue heterogeneity. For example, previous studies using single-cell RNA sequencing identified at least 9 major cell types and up to 26 distinct subtypes in the dorsal cortex alone^1,2. Here we characterize cell-type-specific cis-regulatory chromatin interactions, open chromatin peaks, and transcriptomes for radial glia, intermediate progenitor cells, excitatory neurons, and interneurons isolated from mid-gestational samples of the human cortex. We show that chromatin interactions underlie several aspects of gene regulation, with transposable elements and disease-associated variants enriched at distal interacting regions in a cell-type-specific manner. In addition, promoters with increased levels of chromatin interactivity-termed super-interactive promoters-are enriched for lineage-specific genes, suggesting that interactions at these loci contribute to the fine-tuning of transcription. Finally, we develop CRISPRview, a technique that integrates immunostaining, CRISPR interference, RNAscope, and image analysis to validate cell-type-specific cis-regulatory elements in heterogeneous populations of primary cells. Our findings provide insights into cell-type-specific gene expression patterns in the developing human cortex and advance our understanding of gene regulation and lineage specification during this crucial developmental window.