Single-mitochondrion sequencing uncovers distinct mutational patterns and heteroplasmy landscape in mouse astrocytes and neurons.

BACKGROUND: Mitochondrial (mt) heteroplasmy can cause adverse biological consequences when deleterious mtDNA mutations accumulate disrupting "normal" mt-driven processes and cellular functions. To investigate the heteroplasmy of such mtDNA changes, we developed a moderate throughput mt isolation procedure to quantify the mt single-nucleotide variant (SNV) landscape in individual mouse neurons and astrocytes. In this study, we amplified mt-genomes from 1645 single mitochondria isolated from mouse single astrocytes and neurons to (1) determine the distribution and proportion of mt-SNVs as well as mutation pattern in specific target regions across the mt-genome, (2) assess differences in mtDNA SNVs between neurons and astrocytes, and (3) study co-segregation of variants in the mouse mtDNA. RESULTS: (1) The data show that specific sites of the mt-genome are permissive to SNV presentation while others appear to be under stringent purifying selection. Nested hierarchical analysis at the levels of mitochondrion, cell, and mouse reveals distinct patterns of inter- and intra-cellular variation for mt-SNVs at different sites. (2) Further, differences in the SNV incidence were observed between mouse neurons and astrocytes for two mt-SNV 9027:G > A and 9419:C > T showing variation in the mutational propensity between these cell types. Purifying selection was observed in neurons as shown by the Ka/Ks statistic, suggesting that neurons are under stronger evolutionary constraint as compared to astrocytes. (3) Intriguingly, these data show strong linkage between the SNV sites at nucleotide positions 9027 and 9461. CONCLUSIONS: This study suggests that segregation as well as clonal expansion of mt-SNVs is specific to individual genomic loci, which is important foundational data in understanding of heteroplasmy and disease thresholds for mutation of pathogenic variants.

Stepwise developmental mimicry generates proximal-biased kidney organoids.

The kidney maintains body fluid homeostasis by reabsorbing essential compounds and excreting waste. Proximal tubule cells, crucial for renal reabsorption of a range of sugars, ions, and amino acids, are highly susceptible to damage, leading to pathologies necessitating dialysis and kidney transplants. While human pluripotent stem cell-derived kidney organoids are used for modeling renal development, disease, and injury, the formation of proximal nephron cells in these 3D structures is incomplete. Here, we describe how to drive the development of proximal tubule precursors in kidney organoids by following a blueprint of in vivo human nephrogenesis. Transient manipulation of the PI3K signaling pathway activates Notch signaling in the early nephron and drives nephrons toward a proximal precursor state. These "proximal-biased" (PB) organoid nephrons proceed to generate proximal nephron precursor cells. Single-cell transcriptional analyses across the organoid nephron differentiation, comparing control and PB types, confirm the requirement of transient Notch signaling for proximal development. Indicative of functional maturity, PB organoids demonstrate dextran and albumin uptake, akin to in vivo proximal tubules. Moreover, PB organoids are highly sensitive to nephrotoxic agents, display an injury response, and drive expression of HAVCR1 / KIM1 , an early proximal-specific marker of kidney injury. Injured PB organoids show evidence of collapsed tubules, DNA damage, and upregulate the injury-response marker SOX9 . The PB organoid model therefore has functional relevance and potential for modeling mechanisms underpinning nephron injury. These advances improve the use of iPSC-derived kidney organoids as tools to understand developmental nephrology, model disease, test novel therapeutics, and for understanding human renal physiology.

Single-Mitochondrion Sequencing Uncovers Distinct Mutational Patterns and Heteroplasmy Landscape in Mouse Astrocytes and Neurons.

Background: Mitochondrial (mt) heteroplasmy can cause adverse biological consequences when deleterious mtDNA mutations accumulate disrupting 'normal' mt-driven processes and cellular functions. To investigate the heteroplasmy of such mtDNA changes we developed a moderate throughput mt isolation procedure to quantify the mt single-nucleotide variant (SNV) landscape in individual mouse neurons and astrocytes In this study we amplified mt-genomes from 1,645 single mitochondria (mts) isolated from mouse single astrocytes and neurons to 1. determine the distribution and proportion of mt-SNVs as well as mutation pattern in specific target regions across the mt-genome, 2. assess differences in mtDNA SNVs between neurons and astrocytes, and 3. Study cosegregation of variants in the mouse mtDNA. Results: 1. The data show that specific sites of the mt-genome are permissive to SNV presentation while others appear to be under stringent purifying selection. Nested hierarchical analysis at the levels of mitochondrion, cell, and mouse reveals distinct patterns of inter- and intra-cellular variation for mt-SNVs at different sites. 2. Further, differences in the SNV incidence were observed between mouse neurons and astrocytes for two mt-SNV 9027:G>A and 9419:C>T showing variation in the mutational propensity between these cell types. Purifying selection was observed in neurons as shown by the Ka/Ks statistic, suggesting that neurons are under stronger evolutionary constraint as compared to astrocytes. 3. Intriguingly, these data show strong linkage between the SNV sites at nucleotide positions 9027 and 9461. Conclusion: This study suggests that segregation as well as clonal expansion of mt-SNVs is specific to individual genomic loci, which is important foundational data in understanding of heteroplasmy and disease thresholds for mutation of pathogenic variants.

SCO-pH: Microfluidic dynamic phenotyping platform for high-throughput screening of single cell acidification.

Studies on the dynamics of single cell phenotyping have been hampered by the lack of quantitative high-throughput metabolism assays. Extracellular acidification, a prominent phenotype, yields significant insights into cellular metabolism, including tumorigenicity. Here, we develop a versatile microfluidic system for single cell optical pH analysis (SCO-pH), which compartmentalizes single cells in 140-pL droplets and immobilizes approximately 40,000 droplets in a two-dimensional array for temporal extracellular pH analysis. SCO-pH distinguishes cells undergoing hyperglycolysis induced by oligomycin A from untreated cells by monitoring their extracellular acidification. To facilitate pH sensing in each droplet, we encapsulate a cell-impermeable pH probe whose fluorescence intensities are quantified. Using this approach, we can differentiate hyperglycolytic cells and concurrently observe single cell heterogeneity in extracellular acidification dynamics. This high-throughput system will be useful in applications that require dynamic phenotyping of single cells with significant heterogeneity.

Uniform quantification of single-nucleus ATAC-seq data with Paired-Insertion Counting (PIC) and a model-based insertion rate estimator.

Existing approaches to scoring single-nucleus assay for transposase-accessible chromatin with sequencing (snATAC-seq) feature matrices from sequencing reads are inconsistent, affecting downstream analyses and displaying artifacts. We show that, even with sparse single-cell data, quantitative counts are informative for estimating the regulatory state of a cell, which calls for a consistent treatment. We propose Paired-Insertion Counting as a uniform method for snATAC-seq feature characterization and provide a probability model for inferring latent insertion dynamics from snATAC-seq count matrices.

PACS allows comprehensive dissection of multiple factors governing chromatin accessibility from snATAC-seq data.

Single nucleus ATAC-seq (snATAC-seq) experimental designs have become increasingly complex with multiple factors that might affect chromatin accessibility, including genotype, cell type, tissue of origin, sample location, batch, etc., whose compound effects are difficult to test by existing methods. In addition, current snATAC-seq data present statistical difficulties due to their sparsity and variations in individual sequence capture. To address these problems, we present a zero-adjusted statistical model, Probability model of Accessible Chromatin of Single cells (PACS), that can allow complex hypothesis testing of factors that affect accessibility while accounting for sparse and incomplete data. For differential accessibility analysis, PACS controls the false positive rate and achieves on average a 17% to 122% higher power than existing tools. We demonstrate the effectiveness of PACS through several analysis tasks including supervised cell type annotation, compound hypothesis testing, batch effect correction, and spatiotemporal modeling. We apply PACS to several datasets from a variety of tissues and show its ability to reveal previously undiscovered insights in snATAC-seq data.

CHEX-seq detects single-cell genomic single-stranded DNA with catalytical potential.

Genomic DNA (gDNA) undergoes structural interconversion between single- and double-stranded states during transcription, DNA repair and replication, which is critical for cellular homeostasis. We describe "CHEX-seq" which identifies the single-stranded DNA (ssDNA) in situ in individual cells. CHEX-seq uses 3'-terminal blocked, light-activatable probes to prime the copying of ssDNA into complementary DNA that is sequenced, thereby reporting the genome-wide single-stranded chromatin landscape. CHEX-seq is benchmarked in human K562 cells, and its utilities are demonstrated in cultures of mouse and human brain cells as well as immunostained spatially localized neurons in brain sections. The amount of ssDNA is dynamically regulated in response to perturbation. CHEX-seq also identifies single-stranded regions of mitochondrial DNA in single cells. Surprisingly, CHEX-seq identifies single-stranded loci in mouse and human gDNA that catalyze porphyrin metalation in vitro, suggesting a catalytic activity for genomic ssDNA. We posit that endogenous DNA enzymatic activity is a function of genomic ssDNA.

HRS mediates tumor immune evasion by regulating proteostasis-associated interferon pathway activation.

By sorting receptor tyrosine kinases into endolysosomes, the endosomal sorting complexes required for transport (ESCRTs) are thought to attenuate oncogenic signaling in tumor cells. Paradoxically, ESCRT members are upregulated in tumors. Here, we show that disruption of hepatocyte growth factor-regulated tyrosine kinase substrate (HRS), a pivotal ESCRT component, inhibited tumor growth by promoting CD8+ T cell infiltration in melanoma and colon cancer mouse models. HRS ablation led to misfolded protein accumulation and triggered endoplasmic reticulum (ER) stress, resulting in the activation of the type I interferon pathway in an inositol-requiring enzyme-1α (IRE1α)/X-box binding protein 1 (XBP1)-dependent manner. HRS was upregulated in tumor cells with high tumor mutational burden (TMB). HRS expression associates with the response to PD-L1/PD-1 blockade therapy in melanoma patients with high TMB tumors. HRS ablation sensitized anti-PD-1 treatment in mouse melanoma models. Our study shows a mechanism by which tumor cells with high TMB evade immune surveillance and suggests HRS as a promising target to improve immunotherapy.

Upregulation of exosome secretion from tumor-associated macrophages plays a key role in the suppression of anti-tumor immunity.

Macrophages play a pivotal role in tumor immunity. We report that reprogramming of macrophages to tumor-associated macrophages (TAMs) promotes the secretion of exosomes. Mechanistically, increased exosome secretion is driven by MADD, which is phosphorylated by Akt upon TAM induction and activates Rab27a. TAM exosomes carry high levels of programmed death-ligand 1 (PD-L1) and potently suppress the proliferation and function of CD8+ T cells. Analysis of patient melanoma tissues indicates that TAM exosomes contribute significantly to CD8+ T cell suppression. Single-cell RNA sequencing analysis showed that exosome-related genes are highly expressed in macrophages in melanoma; TAM-specific RAB27A expression inversely correlates with CD8+ T cell infiltration. In a murine melanoma model, lipid nanoparticle delivery of small interfering RNAs (siRNAs) targeting macrophage RAB27A led to better T cell activation and sensitized tumors to anti-programmed cell death protein 1 (PD-1) treatment. Our study demonstrates tumors use TAM exosomes to combat CD8 T cells and suggests targeting TAM exosomes as a potential strategy to improve immunotherapies.

Transcriptional profiling of dental sensory and proprioceptive trigeminal neurons using single-cell RNA sequencing.

Dental primary afferent (DPA) neurons and proprioceptive mesencephalic trigeminal nucleus (MTN) neurons, located in the trigeminal ganglion and the brainstem, respectively, are essential for controlling masticatory functions. Despite extensive transcriptomic studies on various somatosensory neurons, there is still a lack of knowledge about the molecular identities of these populations due to technical challenges in their circuit-validated isolation. Here, we employed high-depth single-cell RNA sequencing (scRNA-seq) in combination with retrograde tracing in mice to identify intrinsic transcriptional features of DPA and MTN neurons. Our transcriptome analysis revealed five major types of DPA neurons with cell type-specific gene enrichment, some of which exhibit unique mechano-nociceptive properties capable of transmitting nociception in response to innocuous mechanical stimuli in the teeth. Furthermore, we discovered cellular heterogeneity within MTN neurons that potentially contribute to their responsiveness to mechanical stretch in the masseter muscle spindles. Additionally, DPA and MTN neurons represented sensory compartments with distinct molecular profiles characterized by various ion channels, receptors, neuropeptides, and mechanoreceptors. Together, our study provides new biological insights regarding the highly specialized mechanosensory functions of DPA and MTN neurons in pain and proprioception.

Direct androgen receptor control of sexually dimorphic gene expression in the mammalian kidney.

Mammalian organs exhibit distinct physiology, disease susceptibility, and injury responses between the sexes. In the mouse kidney, sexually dimorphic gene activity maps predominantly to proximal tubule (PT) segments. Bulk RNA sequencing (RNA-seq) data demonstrated that sex differences were established from 4 and 8 weeks after birth under gonadal control. Hormone injection studies and genetic removal of androgen and estrogen receptors demonstrated androgen receptor (AR)-mediated regulation of gene activity in PT cells as the regulatory mechanism. Interestingly, caloric restriction feminizes the male kidney. Single-nuclear multiomic analysis identified putative cis-regulatory regions and cooperating factors mediating PT responses to AR activity in the mouse kidney. In the human kidney, a limited set of genes showed conserved sex-linked regulation, whereas analysis of the mouse liver underscored organ-specific differences in the regulation of sexually dimorphic gene expression. These findings raise interesting questions on the evolution, physiological significance, disease, and metabolic linkage of sexually dimorphic gene activity.

Systematic Sampling of the Female Reproductive System for Molecular Characterization.

As part of the National Institutes of Health Human BioMolecular Atlas Program to develop a global platform to map the 37 trillion cells in the adult human body, we are generating a comprehensive molecular characterization of the female reproductive system. Data gathered from multiple single-cell/single-nucleus and spatial molecular assays will be used to build a 3D molecular atlas. Herein, we describe our multistep protocol, beginning with an optimized organ procurement workflow that maintains functional characteristics of the uterus, ovaries, and fallopian tubes by perfusing these organs with preservation solution. We have also developed a structured tissue sampling procedure that retains information on individual-level anatomic, physiologic, and individual diversity of the female reproductive system, toward full exploration of the function and structure of female reproductive cells. © 2023 Wiley Periodicals LLC. Basic Protocol 1: Preparation and preservation of the female reproductive system (ovaries, fallopian tubes, and uterus) prior to procurement Basic Protocol 2: Removal of the female reproductive system en bloc Basic Protocol 3: Postsurgical dissection of ovaries Basic Protocol 4: Postsurgical dissection of fallopian tubes Basic Protocol 5: Postsurgical dissection of cervix Basic Protocol 6: Postsurgical dissection of uterine body Support Protocol 1: OCT-embedded tissue protocol Support Protocol 2: Tissue fixation protocol Support Protocol 3: Snap-frozen tissue protocol Basic Protocol 7: Tissue slice preparation for Visium analysis Support Protocol 4: Hematoxylin and eosin staining for 10X Visium imaging Basic Protocol 8: Manual tissue dissociation for Multiome analysis Basic Protocol 9: Tissue dissociation for Multiome analysis using S2 Singulator.

Direct androgen receptor regulation of sexually dimorphic gene expression in the mammalian kidney.

Mammalian organs exhibit distinct physiology, disease susceptibility and injury responses between the sexes. In the mouse kidney, sexually dimorphic gene activity maps predominantly to proximal tubule (PT) segments. Bulk RNA-seq data demonstrated sex differences were established from 4 and 8 weeks after birth under gonadal control. Hormone injection studies and genetic removal of androgen and estrogen receptors demonstrated androgen receptor (AR) mediated regulation of gene activity in PT cells as the regulatory mechanism. Interestingly, caloric restriction feminizes the male kidney. Single-nuclear multiomic analysis identified putative cis-regulatory regions and cooperating factors mediating PT responses to AR activity in the mouse kidney. In the human kidney, a limited set of genes showed conserved sex-linked regulation while analysis of the mouse liver underscored organ-specific differences in the regulation of sexually dimorphic gene expression. These findings raise interesting questions on the evolution, physiological significance, and disease and metabolic linkage, of sexually dimorphic gene activity.

Astrocyte ethanol exposure reveals persistent and defined calcium response subtypes and associated gene signatures.

Astrocytes play a critical role in brain function, but their contribution during ethanol (EtOH) consumption remains largely understudied. In light of recent findings on the heterogeneity of astrocyte physiology and gene expression, an approach with the ability to identify subtypes and capture this heterogeneity is necessary. Here, we combined measurements of calcium signaling and gene expression to define EtOH-induced astrocyte subtypes. In the absence of a demonstrated EtOH receptor, EtOH is believed to have effects on the function of many receptors and downstream biological cascades that underlie calcium responsiveness. This mechanism of EtOH-induced calcium signaling is unknown and this study provides the first step in understanding the characteristics of cells displaying these observed responses. To characterize underlying astrocyte subtypes, we assessed the correlation between calcium signaling and astrocyte gene expression signature in response to EtOH. We found that various EtOH doses increased intracellular calcium levels in a subset of astrocytes, distinguishing three cellular response types and one nonresponsive subtype as categorized by response waveform properties. Furthermore, single-cell RNA-seq analysis of astrocytes from the different response types identified type-enriched discriminatory gene expression signatures. Combining single-cell calcium responses and gene expression analysis identified specific astrocyte subgroups among astrocyte populations defined by their response to EtOH. This result provides a basis for identifying the relationship between astrocyte susceptibility to EtOH and corresponding measurable markers of calcium signaling and gene expression, which will be useful to investigate potential subgroup-specific influences of astrocytes on the physiology and pathology of EtOH exposure in the brain.

Learning Proteome Domain Folding Using LSTMs in an Empirical Kernel Space.

The recognition of protein structural folds is the starting point for protein function inference and for many structural prediction tools. We previously introduced the idea of using empirical comparisons to create a data-augmented feature space called PESS (Protein Empirical Structure Space)1 as a novel approach for protein structure prediction. Here, we extend the previous approach by generating the PESS feature space over fixed-length subsequences of query peptides, and applying a sequential neural network model, with one long short-term memory cell layer followed by a fully connected layer. Using this approach, we show that only a small group of domains as a training set is needed to achieve near state-of-the-art accuracy on fold recognition. Our method improves on the previous approach by reducing the training set required and improving the model's ability to generalize across species, which will help fold prediction for newly discovered proteins.

Single-cell analysis identifies the interaction of altered renal tubules with basophils orchestrating kidney fibrosis.

Inflammation is an important component of fibrosis but immune processes that orchestrate kidney fibrosis are not well understood. Here we apply single-cell sequencing to a mouse model of kidney fibrosis. We identify a subset of kidney tubule cells with a profibrotic-inflammatory phenotype characterized by the expression of cytokines and chemokines associated with immune cell recruitment. Receptor-ligand interaction analysis and experimental validation indicate that CXCL1 secreted by profibrotic tubules recruits CXCR2+ basophils. In mice, these basophils are an important source of interleukin-6 and recruitment of the TH17 subset of helper T cells. Genetic deletion or antibody-based depletion of basophils results in reduced renal fibrosis. Human kidney single-cell, bulk gene expression and immunostaining validate a function for basophils in patients with kidney fibrosis. Collectively, these studies identify basophils as contributors to the development of renal fibrosis and suggest that targeting these cells might be a useful clinical strategy to manage chronic kidney disease.

Multi-omics integration in the age of million single-cell data.

An explosion in single-cell technologies has revealed a previously underappreciated heterogeneity of cell types and novel cell-state associations with sex, disease, development and other processes. Starting with transcriptome analyses, single-cell techniques have extended to multi-omics approaches and now enable the simultaneous measurement of data modalities and spatial cellular context. Data are now available for millions of cells, for whole-genome measurements and for multiple modalities. Although analyses of such multimodal datasets have the potential to provide new insights into biological processes that cannot be inferred with a single mode of assay, the integration of very large, complex, multimodal data into biological models and mechanisms represents a considerable challenge. An understanding of the principles of data integration and visualization methods is required to determine what methods are best applied to a particular single-cell dataset. Each class of method has advantages and pitfalls in terms of its ability to achieve various biological goals, including cell-type classification, regulatory network modelling and biological process inference. In choosing a data integration strategy, consideration must be given to whether the multi-omics data are matched (that is, measured on the same cell) or unmatched (that is, measured on different cells) and, more importantly, the overall modelling and visualization goals of the integrated analysis.

Parallel pathways for recruiting effector proteins determine centromere drive and suppression.

Selfish centromere DNA sequences bias their transmission to the egg in female meiosis. Evolutionary theory suggests that centromere proteins evolve to suppress costs of this "centromere drive." In hybrid mouse models with genetically different maternal and paternal centromeres, selfish centromere DNA exploits a kinetochore pathway to recruit microtubule-destabilizing proteins that act as drive effectors. We show that such functional differences are suppressed by a parallel pathway for effector recruitment by heterochromatin, which is similar between centromeres in this system. Disrupting the kinetochore pathway with a divergent allele of CENP-C reduces functional differences between centromeres, whereas disrupting heterochromatin by CENP-B deletion amplifies the differences. Molecular evolution analyses using Murinae genomes identify adaptive evolution in proteins in both pathways. We propose that centromere proteins have recurrently evolved to minimize the kinetochore pathway, which is exploited by selfish DNA, relative to the heterochromatin pathway that equalizes centromeres, while maintaining essential functions.

Single cell regulatory landscape of the mouse kidney highlights cellular differentiation programs and disease targets.

Determining the epigenetic program that generates unique cell types in the kidney is critical for understanding cell-type heterogeneity during tissue homeostasis and injury response. Here, we profile open chromatin and gene expression in developing and adult mouse kidneys at single cell resolution. We show critical reliance of gene expression on distal regulatory elements (enhancers). We reveal key cell type-specific transcription factors and major gene-regulatory circuits for kidney cells. Dynamic chromatin and expression changes during nephron progenitor differentiation demonstrates that podocyte commitment occurs early and is associated with sustained Foxl1 expression. Renal tubule cells follow a more complex differentiation, where Hfn4a is associated with proximal and Tfap2b with distal fate. Mapping single nucleotide variants associated with human kidney disease implicates critical cell types, developmental stages, genes, and regulatory mechanisms. The single cell multi-omics atlas reveals key chromatin remodeling events and gene expression dynamics associated with kidney development.