Enhlink infers distal and context-specific enhancer-promoter linkages.

Enhancers play a crucial role in regulating gene expression and their functional status can be queried with cell type precision using using single-cell (sc)ATAC-seq. To facilitate analysis of such data, we developed Enhlink, a novel computational approach that leverages single-cell signals to infer linkages between regulatory DNA sequences, such as enhancers and promoters. Enhlink uses an ensemble strategy that integrates cell-level technical covariates to control for batch effects and biological covariates to infer robust condition-specific links and their associated p-values. It can integrate simultaneous gene expression and chromatin accessibility measurements of individual cells profiled by multi-omic experiments for increased specificity. We evaluated Enhlink using simulated and real scATAC-seq data, including those paired with physical enhancer-promoter links enumerated by promoter capture Hi-C and with multi-omic scATAC-/RNA-seq data we generated from the mouse striatum. These examples demonstrated that our method outperforms popular alternative strategies. In conjunction with eQTL analysis, Enhlink revealed a putative super-enhancer regulating key cell type-specific markers of striatal neurons. Taken together, our analyses demonstrate that Enhlink is accurate, powerful, and provides features that can lead to novel biological insights.

Adult mouse fibroblasts retain organ-specific transcriptomic identity.

Organ fibroblasts are essential components of homeostatic and diseased tissues. They participate in sculpting the extracellular matrix, sensing the microenvironment, and communicating with other resident cells. Recent studies have revealed transcriptomic heterogeneity among fibroblasts within and between organs. To dissect the basis of interorgan heterogeneity, we compare the gene expression of murine fibroblasts from different tissues (tail, skin, lung, liver, heart, kidney, and gonads) and show that they display distinct positional and organ-specific transcriptome signatures that reflect their embryonic origins. We demonstrate that expression of genes typically attributed to the surrounding parenchyma by fibroblasts is established in embryonic development and largely maintained in culture, bioengineered tissues and ectopic transplants. Targeted knockdown of key organ-specific transcription factors affects fibroblast functions, in particular genes involved in the modulation of fibrosis and inflammation. In conclusion, our data reveal that adult fibroblasts maintain an embryonic gene expression signature inherited from their organ of origin, thereby increasing our understanding of adult fibroblast heterogeneity. The knowledge of this tissue-specific gene signature may assist in targeting fibrotic diseases in a more precise, organ-specific manner.

Protocol for Isolation of Cardiac Interstitial Cells from Adult Murine Hearts for Unbiased Single Cell Profiling.

Interstitial cells have a crucial role in cardiac fibrosis and repair of the mammalian heart. Single-cell profiling using droplet-based technology has revolutionized the investigation of cell states and identities. Here, we present a protocol for the efficient isolation of high-quality live nucleated non-cardiomyocytes from adult murine heart, for unbiased single-cell RNA sequencing using 10× Chromium technology. This protocol has been applied to homeostatic and injured hearts from different mouse strains. For complete details on the use and execution of this protocol, please refer to Forte et al. (2020).

Dynamic Interstitial Cell Response during Myocardial Infarction Predicts Resilience to Rupture in Genetically Diverse Mice.

Cardiac ischemia leads to the loss of myocardial tissue and the activation of a repair process that culminates in the formation of a scar whose structural characteristics dictate propensity to favorable healing or detrimental cardiac wall rupture. To elucidate the cellular processes underlying scar formation, here we perform unbiased single-cell mRNA sequencing of interstitial cells isolated from infarcted mouse hearts carrying a genetic tracer that labels epicardial-derived cells. Sixteen interstitial cell clusters are revealed, five of which were of epicardial origin. Focusing on stromal cells, we define 11 sub-clusters, including diverse cell states of epicardial- and endocardial-derived fibroblasts. Comparing transcript profiles from post-infarction hearts in C57BL/6J and 129S1/SvImJ inbred mice, which displays a marked divergence in the frequency of cardiac rupture, uncovers an early increase in activated myofibroblasts, enhanced collagen deposition, and persistent acute phase response in 129S1/SvImJ mouse hearts, defining a crucial time window of pathological remodeling that predicts disease outcome.

A haplolethal locus uncovered by deletions in the mouse T complex.

Proper levels of gene expression are important for normal mammalian development. Typically, altered gene dosage caused by karyotypic abnormalities results in embryonic lethality or birth defects. Segmental aneuploidy can be compatible with life but often results in contiguous gene syndromes. The ability to manipulate the mouse genome allows the systematic exploration of regions that are affected by alterations in gene dosage. To explore the effects of segmental haploidy in the mouse t complex on chromosome 17, radiation-induced deletion complexes centered at the Sod2 and D17Leh94 loci were generated in embryonic stem (ES) cells. A small interval was identified that, when hemizygous, caused specific embryonic lethal phenotypes (exencephaly and edema) in most fetuses. The penetrance of these phenotypes was background dependent. Additionally, evidence for parent-of-origin effects was observed. This genetic approach should be useful for identifying genes that are imprinted or whose dosage is critical for normal embryonic development.

A comparison of cDNA, oligonucleotide, and Affymetrix GeneChip gene expression microarray platforms.

We have conducted a study to compare the variability in measured gene expression levels associated with three types of microarray platforms. Total RNA samples were obtained from liver tissue of four male mice, two each from inbred strains A/J and C57BL/6J. The same four samples were assayed on Affymetrix Mouse Genome Expression Set 430 GeneChips (MOE430A and MOE430B), spotted cDNA microarrays, and spotted oligonucleotide microarrays using eight arrays of each type. Variances associated with measurement error were observed to be comparable across all microarray platforms. The MOE430A GeneChips and cDNA arrays had higher precision across technical replicates than the MOE430B GeneChips and oligonucleotide arrays. The Affymetrix platform showed the greatest range in the magnitude of expression levels followed by the oligonucleotide arrays. We observed good concordance in both estimated expression level and statistical significance of common genes between the Affymetrix MOE430A GeneChip and the oligonucleotide arrays. Despite their apparently high precision, cDNA arrays showed poor concordance with other platforms.