Complementation testing identifies genes mediating effects at quantitative trait loci underlying fear-related behavior

Knowing the genes involved in quantitative traits provides an entry point to understanding the biological bases of behavior, but there are very few examples where the pathway from genetic locus to behavioral change is known. To explore the role of specific genes in fear behavior, we mapped three fear-related traits, tested fourteen genes at six quantitative trait loci (QTLs) by quantitative complementation, and identified six genes. Four genes, Lamp, Ptprd, Nptx2, and Sh3gl, have known roles in synapse function; the fifth, Psip1, was not previously implicated in behavior; and the sixth is a long non-coding RNA, 4933413L06Rik, of unknown function. Variation in transcriptome and epigenetic modalities occurred preferentially in excitatory neurons, suggesting that genetic variation is more permissible in excitatory than inhibitory neuronal circuits. Our results relieve a bottleneck in using genetic mapping of QTLs to uncover biology underlying behavior and prompt a reconsideration of expected relationships between genetic and functional variation.

Polyploidisation pleiotropically buffers ageing in hepatocytes

BACKGROUND & AIMS

Polyploidy in hepatocytes has been proposed as a genetic mechanism to buffer against transcriptional dysregulation. Here, we aim to demonstrate the role of polyploidy in modulating gene regulatory networks in hepatocytes during ageing.

METHODS

We performed single-nucleus RNA-sequencing in hepatocyte nuclei of different ploidy levels isolated from young and old wild-type mice. Changes in the gene expression and regulatory network were compared to three independent haploinsufficient strains for HNF4A, CEBPA or CTCF, representing non-deleterious perturbations. Phenotypic characteristics of the liver section were additionally evaluated histologically, whereas the genomic allele composition of hepatocytes was analysed by BaseScope.

RESULTS

We observed that ageing in wild-type mice results in nuclei polyploidy and marked increase in steatosis. Haploinsufficiency of liver-specific master regulators (HFN4A or CEBPA) results in the enrichment of hepatocytes with tetraploid nuclei at a young age, affecting the genomic regulatory network, and dramatically suppressing ageing-related steatosis tissue-wide. Notably, these phenotypes are not the result of subtle disruption to liver-specific transcriptional networks, since haploinsufficiency in CTCF insulator protein resulted in the same phenotype. Further quantification of genotypes of tetraploid hepatocytes in young and old HFN4A haploinsufficient mice revealed that during ageing, tetraploid hepatocytes lead to the selection of wild-type alleles, restoring non-deleterious genetic perturbation. ConclusionsOur results suggest a model whereby polyploidisation leads to fundamentally different cell states. Polyploid conversion enables pleiotropic buffering against age-related decline via non-random allelic segregation to restore a wild-type genome.

Whole organism snRNA-seq reveals systemic peripheral changes in Alzheimer's Disease fly models

Peripheral tissues become disrupted in Alzheimer's Disease (AD). However, a comprehensive understanding of how the expression of AD-associated toxic proteins, Aβ42 and Tau, in neurons impacts the periphery is lacking. Using Drosophila, a prime model organism for studying aging and neurodegeneration, we generated the Alzheimer's Disease Fly Cell Atlas (AD-FCA): whole-organism single-nucleus transcriptomes of 219 cell types from adult flies neuronally expressing human Aβ42 or Tau. In-depth analyses and functional data reveal impacts on peripheral sensory neurons by Aβ42 and on various non-neuronal peripheral tissues by Tau, including the gut, fat body, and reproductive system. This novel AD atlas provides valuable insights into potential biomarkers and the intricate interplay between the nervous system and peripheral tissues in response to AD-associated proteins.

Transplanted human neural stem cells rescue phenotypes in zQ175 Huntington's disease mice and innervate the striatum

Huntington's disease (HD), a genetic neurodegenerative disorder, primarily affects the striatum and cortex with progressive loss of medium-sized spiny neurons (MSNs) and pyramidal neurons, disrupting cortico-striatal circuitry. A promising regenerative therapeutic strategy of transplanting human neural stem cells (hNSCs) is challenged by the need for long-term functional integration. We previously described that, with short-term hNSC transplantation into the striatum of HD R6/2 mice, human cells differentiated into electrophysiologically active immature neurons, improving behavior and biochemical deficits. Here, we show that long-term (8 months) implantation of hNSCs into the striatum of HD zQ175 mice ameliorates behavioral deficits, increases brain-derived neurotrophic factor (BDNF) levels, and reduces mutant huntingtin (mHTT) accumulation. Patch clamp recordings, immunohistochemistry, single-nucleus RNA sequencing (RNA-seq), and electron microscopy demonstrate that hNSCs differentiate into diverse neuronal populations, including MSN- and interneuron-like cells, and form connections. Single-nucleus RNA-seq analysis also shows restoration of several mHTT-mediated transcriptional changes of endogenous striatal HD mouse cells. Remarkably, engrafted cells receive synaptic inputs, innervate host neurons, and improve membrane and synaptic properties. Overall, the findings support hNSC transplantation for further evaluation and clinical development for HD.

Determination of key events in mouse hepatocyte maturation at the single-cell level

Hepatocytes, the liver's predominant cells, perform numerous essential biological functions. However, crucial events and regulators during hepatocyte maturation require in-depth investigation. In this study, we performed single-cell RNA sequencing (scRNA-seq) and single-nucleus RNA sequencing (snRNA-seq) to explore the precise hepatocyte development process in mice. We defined three maturation stages of postnatal hepatocytes, each of which establishes specific metabolic functions and exhibits distinct proliferation rates. Hepatic zonation is gradually formed during hepatocyte maturation. Hepatocytes or their nuclei with distinct ploidies exhibit zonation preferences in distribution and asynchrony in maturation. Moreover, by combining gene regulatory network analysis with in vivo genetic manipulation, we identified critical maturation- and zonation-related transcription factors. This study not only delineates the comprehensive transcriptomic profiles of hepatocyte maturation but also presents a paradigm to identify genes that function in the development of hepatocyte maturation and zonation by combining genetic manipulation and measurement of coordinates in a single-cell developmental trajectory.

Opportunities and challenges in the application of single-cell and spatial transcriptomics in plants

Single-cell and spatial transcriptomics have diverted researchers' attention from the multicellular level to the single-cell level and spatial information. Single-cell transcriptomes provide insights into the transcriptome at the single-cell level, whereas spatial transcriptomes help preserve spatial information. Although these two omics technologies are helpful and mature, further research is needed to ensure their widespread applicability in plant studies. Reviewing recent research on plant single-cell or spatial transcriptomics, we compared the different experimental methods used in various plants. The limitations and challenges are clear for both single-cell and spatial transcriptomic analyses, such as the lack of applicability, spatial information, or high resolution. Subsequently, we put forth further applications, such as cross-species analysis of roots at the single-cell level and the idea that single-cell transcriptome analysis needs to be combined with other omics analyses to achieve superiority over individual omics analyses. Overall, the results of this review suggest that combining single-cell transcriptomics, spatial transcriptomics, and spatial element distribution can provide a promising research direction, particularly for plant research.

Integrative single-cell analysis of cardiac and pulmonary sarcoidosis using publicly available cardiac and bronchoalveolar lavage fluid sequencing datasets

Introduction

Cardiac presentation of autoimmune sarcoidosis, known as cardiac sarcoidosis (CS), is a poorly understood disease with high mortality and low diagnosis rate. While CS is an immunological syndrome, little is known about how cardiac parenchymal and stromal cells mediate its pathogenesis. Moreover, while most current sarcoidosis research is based on research in pulmonary sarcoidosis (PS), it remains unclear how much both presentations of sarcoidosis overlap. To tackle these concerns, we leveraged publicly available sarcoidosis transcriptomic datasets.

Methods

Two publicly available bronchoalveolar lavage single-cell RNA sequencing datasets were integrated to analyze PS relative to control. Additionally, two publicly available cardiac single-nucleus RNA sequencing datasets were integrated to analyze CS relative to control. Following integration, we ran cell-cell communication, transcription factor, and differential expression analyses on parenchymal, stromal, and immune subsets identified in our analysis.

Results

Our analysis revealed that there was an expansion of stromal and immune cells in PS and CS. We also observed upregulation of Th17.1 and attenuated activation transcriptional profiles in the immune cells of CS and PS relative to control. Additionally, we found upregulation of pro-inflammatory and pro-fibrotic transcriptional profiles in the cardiac stromal cells of CS relative to control. We also found that cardiomyocytes exhibited upregulated cardiac stress and proliferation transcriptional profiles in CS relative to control.

Conclusions

Our integrative transcriptomic analysis shows that despite tissue-specific differences, there are shared transcriptional trends between CS and PS. It also shows that stromal and parenchymal populations exhibit transcriptional trends that could explain their pathogenic role in CS.

Impact of non-neuronal cells in Alzheimer's disease from a single-nucleus profiling perspective

The role of non-neuronal cells has been relatively overlooked in Alzheimer's disease (AD) neuropathogenesis compared to neuronal cells since the first characterization of the disease. Genome wide-association studies (GWAS) performed in the last few decades have greatly contributed to highlighting the critical impact of non-neuronal cells in AD by uncovering major genetic risk factors that are found largely in these cell types. The recent development of single cell or single nucleus technologies has revolutionized the way we interrogate the transcriptomic and epigenetic profiles of neurons, microglia, astrocytes, oligodendrocytes, pericytes, and endothelial cells simultaneously in the same sample and in an individual manner. Here, we review the latest advances in single-cell/nucleus RNA sequencing and Assay for Transposase-Accessible Chromatin (ATAC) sequencing to more accurately understand the function of non-neuronal cells in AD. We conclude by giving an overview of what still needs to be achieved to better appreciate the interconnected roles of each cell type in the context of AD.

Transcriptomic atlas and interaction networks of brain cells in mouse CNS demyelination and remyelination

Demyelination is a hallmark of multiple sclerosis, leukoencephalopathies, cerebral vasculopathies, and several neurodegenerative diseases. The cuprizone mouse model is widely used to simulate demyelination and remyelination occurring in these diseases. Here, we present a high-resolution single-nucleus RNA sequencing (snRNA-seq) analysis of gene expression changes across all brain cells in this model. We define demyelination-associated oligodendrocytes (DOLs) and remyelination-associated MAFBhi microglia, as well as astrocytes and vascular cells with signatures of altered metabolism, oxidative stress, and interferon response. Furthermore, snRNA-seq provides insights into how brain cell types connect and interact, defining complex circuitries that impact demyelination and remyelination. As an explicative example, perturbation of microglia caused by TREM2 deficiency indirectly impairs the induction of DOLs. Altogether, this study provides a rich resource for future studies investigating mechanisms underlying demyelinating diseases.

Transcriptome analysis reveals tumor microenvironment changes in glioblastoma

A better understanding of transcriptional evolution of IDH-wild-type glioblastoma may be crucial for treatment optimization. Here, we perform RNA sequencing (RNA-seq) (n = 322 test, n = 245 validation) on paired primary-recurrent glioblastoma resections of patients treated with the current standard of care. Transcriptional subtypes form an interconnected continuum in a two-dimensional space. Recurrent tumors show preferential mesenchymal progression. Over time, hallmark glioblastoma genes are not significantly altered. Instead, tumor purity decreases over time and is accompanied by co-increases in neuron and oligodendrocyte marker genes and, independently, tumor-associated macrophages. A decrease is observed in endothelial marker genes. These composition changes are confirmed by single-cell RNA-seq and immunohistochemistry. An extracellular matrix-associated gene set increases at recurrence and bulk, single-cell RNA, and immunohistochemistry indicate it is expressed mainly by pericytes. This signature is associated with significantly worse survival at recurrence. Our data demonstrate that glioblastomas evolve mainly by microenvironment (re-)organization rather than molecular evolution of tumor cells.

Cell Type-Specific Nuclei Markers: The Need for Human Brain Research to Go Nuclear

Identifying and interrogating cell type-specific populations within the heterogeneous milieu of the human brain is paramount to resolving the processes of normal brain homeostasis and the pathogenesis of neurological disorders. While brain cell type-specific markers are well established, most are localized on cellular membranes or within the cytoplasm, with limited literature describing those found in the nucleus. Due to the complex cytoarchitecture of the human brain, immunohistochemical studies require well-defined cell-specific nuclear markers for more precise and efficient quantification of the cellular populations. Furthermore, efficient nuclear markers are required for cell type-specific purification and transcriptomic interrogation of archived human brain tissue through nuclei isolation-based RNA sequencing. To sate the growing demand for robust cell type-specific nuclear markers, we thought it prudent to comprehensively review the current literature to identify and consolidate a novel series of robust cell type-specific nuclear markers that can assist researchers across a range of neuroscientific disciplines. The following review article collates and discusses several key and prospective cell type-specific nuclei markers for each of the major human brain cell types; it then concludes by discussing the potential applications of cell type-specific nuclear workflows and the power of nuclear-based neuroscientific research.

Understanding and evaluating ambiguity in single-cell and single-nucleus RNA-sequencing

Recently, a new modification has been proposed by Hjörleifsson and Sullivan et al. to the model used to classify the splicing status of reads (as spliced (mature), unspliced (nascent), or ambiguous) in single-cell and single-nucleus RNA-seq data. Here, we evaluate both the theoretical basis and practical implementation of the proposed method. The proposed method is highly-conservative, and therefore, unlikely to mischaracterize reads as spliced (mature) or unspliced (nascent) when they are not. However, we find that it leaves a large fraction of reads classified as ambiguous, and, in practice, allocates these ambiguous reads in an all-or-nothing manner, and differently between single-cell and single-nucleus RNA-seq data. Further, as implemented in practice, the ambiguous classification is implicit and based on the index against which the reads are mapped, which leads to several drawbacks compared to methods that consider both spliced (mature) and unspliced (nascent) mapping targets simultaneously - for example, the ability to use confidently assigned reads to rescue ambiguous reads based on shared UMIs and gene targets. Nonetheless, we show that these conservative assignment rules can be obtained directly in existing approaches simply by altering the set of targets that are indexed. To this end, we introduce the spliceu reference and show that its use with alevin-fry recapitulates the more conservative proposed classification. We also observe that, on experimental data, and under the proposed allocation rules for ambiguous UMIs, the difference between the proposed classification scheme and existing conventions appears much smaller than previously reported. We demonstrate the use of the new piscem index for mapping simultaneously against spliced (mature) and unspliced (nascent) targets, allowing classification against the full nascent and mature transcriptome in human or mouse in <3GB of memory. Finally, we discuss the potential of incorporating probabilistic evidence into the inference of splicing status, and suggest that it may provide benefits beyond what can be obtained from discrete classification of UMIs as splicing-ambiguous.

Single-cell transcriptomics: A new tool for studying diabetic kidney disease

The kidney is a complex organ comprising various functional partitions and special cell types that play important roles in maintaining homeostasis in the body. Diabetic kidney disease (DKD) is the leading cause of end-stage renal disease and is an independent risk factor for cardiovascular diseases. Owing to the complexity and heterogeneity of kidney structure and function, the mechanism of DKD development has not been fully elucidated. Single-cell sequencing, including transcriptomics, epigenetics, metabolomics, and proteomics etc., is a powerful technology that enables the analysis of specific cell types and states, specifically expressed genes or pathways, cell differentiation trajectories, intercellular communication, and regulation or co-expression of genes in various diseases. Compared with other omics, RNA sequencing is a more developed technique with higher utilization of tissues or samples. This article reviewed the application of single-cell transcriptomics in the field of DKD and highlighted the key signaling pathways in specific tissues or cell types involved in the occurrence and development of DKD. The comprehensive understanding of single-cell transcriptomics through single-cell RNA-seq and single-nucleus RNA-seq will provide us new insights into the pathogenesis and treatment strategy of various diseases including DKD.

Molecular signatures underlying neurofibrillary tangle susceptibility in Alzheimer's disease

Tau aggregation in neurofibrillary tangles (NFTs) is closely associated with neurodegeneration and cognitive decline in Alzheimer's disease (AD). However, the molecular signatures that distinguish between aggregation-prone and aggregation-resistant cell states are unknown. We developed methods for the high-throughput isolation and transcriptome profiling of single somas with NFTs from the human AD brain, quantified the susceptibility of 20 neocortical subtypes for NFT formation and death, and identified both shared and cell-type-specific signatures. NFT-bearing neurons shared a marked upregulation of synaptic transmission-related genes, including a core set of 63 genes enriched for synaptic vesicle cycling. Oxidative phosphorylation and mitochondrial dysfunction were highly cell-type dependent. Apoptosis was only modestly enriched, and the susceptibilities of NFT-bearing and NFT-free neurons for death were highly similar. Our analysis suggests that NFTs represent cell-type-specific responses to stress and synaptic dysfunction. We provide a resource for biomarker discovery and the investigation of tau-dependent and tau-independent mechanisms of neurodegeneration.

Heterogeneous Skeletal Muscle Cell and Nucleus Populations Identified by Single-Cell and Single-Nucleus Resolution Transcriptome Assays

Single-cell RNA-seq (scRNA-seq) has revolutionized modern genomics, but the large size of myotubes and myofibers has restricted use of scRNA-seq in skeletal muscle. For the study of muscle, single-nucleus RNA-seq (snRNA-seq) has emerged not only as an alternative to scRNA-seq, but as a novel method providing valuable insights into multinucleated cells such as myofibers. Nuclei within myofibers specialize at junctions with other cell types such as motor neurons. Nuclear heterogeneity plays important roles in certain diseases such as muscular dystrophies. We survey current methods of high-throughput single cell and subcellular resolution transcriptomics, including single-cell and single-nucleus RNA-seq and spatial transcriptomics, applied to satellite cells, myoblasts, myotubes and myofibers. We summarize the major myonuclei subtypes identified in homeostatic and regenerating tissue including those specific to fiber type or at junctions with other cell types. Disease-specific nucleus populations were found in two muscular dystrophies, FSHD and Duchenne muscular dystrophy, demonstrating the importance of performing transcriptome studies at the single nucleus level in muscle.

Vision-dependent specification of cell types and function in the developing cortex

The role of postnatal experience in sculpting cortical circuitry, while long appreciated, is poorly understood at the level of cell types. We explore this in the mouse primary visual cortex (V1) using single-nucleus RNA sequencing, visual deprivation, genetics, and functional imaging. We find that vision selectively drives the specification of glutamatergic cell types in upper layers (L) (L2/3/4), while deeper-layer glutamatergic, GABAergic, and non-neuronal cell types are established prior to eye opening. L2/3 cell types form an experience-dependent spatial continuum defined by the graded expression of ∼200 genes, including regulators of cell adhesion and synapse formation. One of these genes, Igsf9b, a vision-dependent gene encoding an inhibitory synaptic cell adhesion molecule, is required for the normal development of binocular responses in L2/3. In summary, vision preferentially regulates the development of upper-layer glutamatergic cell types through the regulation of cell-type-specific gene expression programs.

Single-Nucleus RNA-Seq Is Not Suitable for Detection of Microglial Activation Genes in Humans

Single-nucleus RNA sequencing (snRNA-seq) is used as an alternative to single-cell RNA-seq, as it allows transcriptomic profiling of frozen tissue. However, it is unclear whether snRNA-seq is able to detect cellular state in human tissue. Indeed, snRNA-seq analyses of human brain samples have failed to detect a consistent microglial activation signature in Alzheimer’s disease. Our comparison of microglia from single cells and single nuclei of four human subjects reveals that, although most genes show similar relative abundances in cells and nuclei, a small population of genes (∼1%) is depleted in nuclei compared to whole cells. This population is enriched for genes previously implicated in microglial activation, including APOE, CST3, SPP1, and CD74, comprising 18% of previously identified microglial-disease-associated genes. Given the low sensitivity of snRNA-seq to detect many activation genes, we conclude that snRNA-seq is not suited for detecting cellular activation in microglia in human disease.

•

A small set of genes is depleted in microglial nuclei relative to single cells

•

This set is enriched for microglial activation genes, including APOE and SPP1

•

This depletion is confirmed in publicly available datasets

•

Single-nucleus sequencing is not suited for the detection of human microglial activation

Single-Nucleus RNA-Seq Reveals Dysregulation of Striatal Cell Identity Due to Huntington's Disease Mutations

Huntington's disease (HD) is a dominantly inherited neurodegenerative disorder caused by a trinucleotide expansion in exon 1 of the huntingtin (HTT) gene. Cell death in HD occurs primarily in striatal medium spiny neurons (MSNs), but the involvement of specific MSN subtypes and of other striatal cell types remains poorly understood. To gain insight into cell type-specific disease processes, we studied the nuclear transcriptomes of 4524 cells from the striatum of a genetically precise knock-in mouse model of the HD mutation, HttQ175/+, and from wild-type controls. We used 14- to 15-month-old male mice, a time point at which multiple behavioral, neuroanatomical, and neurophysiological changes are present but at which there is no known cell death. Thousands of differentially expressed genes (DEGs) were distributed across most striatal cell types, including transcriptional changes in glial populations that are not apparent from RNA-seq of bulk tissue. Reconstruction of cell type-specific transcriptional networks revealed a striking pattern of bidirectional dysregulation for many cell type-specific genes. Typically, these genes were repressed in their primary cell type, yet de-repressed in other striatal cell types. Integration with existing epigenomic and transcriptomic data suggest that partial loss-of-function of the polycomb repressive complex 2 (PRC2) may underlie many of these transcriptional changes, leading to deficits in the maintenance of cell identity across virtually all cell types in the adult striatum.SIGNIFICANCE STATEMENT Huntington's disease (HD) is a dominantly inherited neurodegenerative disorder characterized by specific loss of medium spiny neurons (MSNs) in the striatum, accompanied by more subtle changes in many other cell types. It is thought that changes in transcriptional regulation are an important underlying mechanism, but cell type-specific gene expression changes are not well understood, particularly at time points relevant to the onset of disease-related symptoms. Single-nucleus (sn)RNA-seq in a genetically precise mouse model enabled us to identify novel patterns of transcriptional dysregulation because of HD mutations, including bidirectional dysregulation of many cell type identity genes that may be driven by partial loss-of-function of the polycomb repressive complex (PRC). Identifying these regulators of transcriptional dysregulation in HD can be leveraged to design novel disease-modifying therapeutics.

Single-cell transcriptomes of developing and adult olfactory receptor neurons in Drosophila

Recognition of environmental cues is essential for the survival of all organisms. Transcriptional changes occur to enable the generation and function of the neural circuits underlying sensory perception. To gain insight into these changes, we generated single-cell transcriptomes of Drosophila olfactory- (ORNs), thermo-, and hygro-sensory neurons at an early developmental and adult stage using single-cell and single-nucleus RNA sequencing. We discovered that ORNs maintain expression of the same olfactory receptors across development. Using receptor expression and computational approaches, we matched transcriptomic clusters corresponding to anatomically and physiologically defined neuron types across multiple developmental stages. We found that cell-type-specific transcriptomes partly reflected axon trajectory choices in development and sensory modality in adults. We uncovered stage-specific genes that could regulate the wiring and sensory responses of distinct ORN types. Collectively, our data reveal transcriptomic features of sensory neuron biology and provide a resource for future studies of their development and physiology.

Increased APOE ε4 expression is associated with the difference in Alzheimer's disease risk from diverse ancestral backgrounds

INTRODUCTION

Apolipoprotein E (APOE) ε4 confers less risk for Alzheimer's disease (AD) in carriers with African local genomic ancestry (ALA) than APOE ε4 carriers with European local ancestry (ELA). Cell type specific transcriptional variation between the two local ancestries (LAs) could contribute to this disease risk differences.

METHODS

Single-nucleus RNA sequencing was performed on frozen frontal cortex of homozygous APOE ε4/ε4 AD patients: seven with ELA, four with ALA.

RESULTS

A total of 60,908 nuclei were sequenced. Within the LA region (chr19:44-46Mb), APOE was the gene most differentially expressed, with ELA carriers having significantly more expression (overall P < 1.8E^-317 ) in 24 of 32 cell clusters. The transcriptome of one astrocyte cluster, with high APOE ε4 expression and specific to ELA, is suggestive of A1 reactive astrocytes.

DISCUSSION

AD patients with ELA expressed significantly greater levels of APOE than ALA APOE ε4 carriers. These differences in APOE expression could contribute to the reduced risk for AD seen in African APOE ε4 carriers.

Single-nucleus transcriptomic survey of cell diversity and functional maturation in postnatal mammalian hearts

A fundamental challenge in understanding cardiac biology and disease is that the remarkable heterogeneity in cell type composition and functional states have not been well characterized at single-cell resolution in maturing and diseased mammalian hearts. Massively parallel single-nucleus RNA sequencing (snRNA-seq) has emerged as a powerful tool to address these questions by interrogating the transcriptome of tens of thousands of nuclei isolated from fresh or frozen tissues. snRNA-seq overcomes the technical challenge of isolating intact single cells from complex tissues, including the maturing mammalian hearts; reduces biased recovery of easily dissociated cell types; and minimizes aberrant gene expression during the whole-cell dissociation. Here we applied sNucDrop-seq, a droplet microfluidics-based massively parallel snRNA-seq method, to investigate the transcriptional landscape of postnatal maturing mouse hearts in both healthy and disease states. By profiling the transcriptome of nearly 20,000 nuclei, we identified major and rare cardiac cell types and revealed significant heterogeneity of cardiomyocytes, fibroblasts, and endothelial cells in postnatal developing hearts. When applied to a mouse model of pediatric mitochondrial cardiomyopathy, we uncovered profound cell type-specific modifications of the cardiac transcriptional landscape at single-nucleus resolution, including changes of subtype composition, maturation states, and functional remodeling of each cell type. Furthermore, we employed sNucDrop-seq to decipher the cardiac cell type-specific gene regulatory network (GRN) of GDF15, a heart-derived hormone and clinically important diagnostic biomarker of heart disease. Together, our results present a rich resource for studying cardiac biology and provide new insights into heart disease using an approach broadly applicable to many fields of biomedicine.