Protective role for kidney TREM2high macrophages in obesity- and diabetes-induced kidney injury.

Diabetic kidney disease (DKD), the most common cause of kidney failure, is a frequent complication of diabetes and obesity, and yet to date, treatments to halt its progression are lacking. We analyze kidney single-cell transcriptomic profiles from DKD patients and two DKD mouse models at multiple time points along disease progression-high-fat diet (HFD)-fed mice aged to 90-100 weeks and BTBR ob/ob mice (a genetic model)-and report an expanding population of macrophages with high expression of triggering receptor expressed on myeloid cells 2 (TREM2) in HFD-fed mice. TREM2high macrophages are enriched in obese and diabetic patients, in contrast to hypertensive patients or healthy controls in an independent validation cohort. Trem2 knockout mice on an HFD have worsening kidney filter damage and increased tubular epithelial cell injury, all signs of worsening DKD. Together, our studies suggest that strategies to enhance kidney TREM2high macrophages may provide therapeutic benefits for DKD.

An atlas of healthy and injured cell states and niches in the human kidney.

Understanding kidney disease relies on defining the complexity of cell types and states, their associated molecular profiles and interactions within tissue neighbourhoods1. Here we applied multiple single-cell and single-nucleus assays (>400,000 nuclei or cells) and spatial imaging technologies to a broad spectrum of healthy reference kidneys (45 donors) and diseased kidneys (48 patients). This has provided a high-resolution cellular atlas of 51 main cell types, which include rare and previously undescribed cell populations. The multi-omic approach provides detailed transcriptomic profiles, regulatory factors and spatial localizations spanning the entire kidney. We also define 28 cellular states across nephron segments and interstitium that were altered in kidney injury, encompassing cycling, adaptive (successful or maladaptive repair), transitioning and degenerative states. Molecular signatures permitted the localization of these states within injury neighbourhoods using spatial transcriptomics, while large-scale 3D imaging analysis (around 1.2 million neighbourhoods) provided corresponding linkages to active immune responses. These analyses defined biological pathways that are relevant to injury time-course and niches, including signatures underlying epithelial repair that predicted maladaptive states associated with a decline in kidney function. This integrated multimodal spatial cell atlas of healthy and diseased human kidneys represents a comprehensive benchmark of cellular states, neighbourhoods, outcome-associated signatures and publicly available interactive visualizations.

FALCON systematically interrogates free fatty acid biology and identifies a novel mediator of lipotoxicity.

Cellular exposure to free fatty acids (FFAs) is implicated in the pathogenesis of obesity-associated diseases. However, there are no scalable approaches to comprehensively assess the diverse FFAs circulating in human plasma. Furthermore, assessing how FFA-mediated processes interact with genetic risk for disease remains elusive. Here, we report the design and implementation of fatty acid library for comprehensive ontologies (FALCON), an unbiased, scalable, and multimodal interrogation of 61 structurally diverse FFAs. We identified a subset of lipotoxic monounsaturated fatty acids associated with decreased membrane fluidity. Furthermore, we prioritized genes that reflect the combined effects of harmful FFA exposure and genetic risk for type 2 diabetes (T2D). We found that c-MAF-inducing protein (CMIP) protects cells from FFA exposure by modulating Akt signaling. In sum, FALCON empowers the study of fundamental FFA biology and offers an integrative approach to identify much needed targets for diverse diseases associated with disordered FFA metabolism.

FALCON systematically interrogates free fatty acid biology and identifies a novel mediator of lipotoxicity.

Cellular exposure to free fatty acids (FFA) is implicated in the pathogenesis of obesity-associated diseases. However, studies to date have assumed that a few select FFAs are representative of broad structural categories, and there are no scalable approaches to comprehensively assess the biological processes induced by exposure to diverse FFAs circulating in human plasma. Furthermore, assessing how these FFA- mediated processes interact with genetic risk for disease remains elusive. Here we report the design and implementation of FALCON (Fatty Acid Library for Comprehensive ONtologies) as an unbiased, scalable and multimodal interrogation of 61 structurally diverse FFAs. We identified a subset of lipotoxic monounsaturated fatty acids (MUFAs) with a distinct lipidomic profile associated with decreased membrane fluidity. Furthermore, we developed a new approach to prioritize genes that reflect the combined effects of exposure to harmful FFAs and genetic risk for type 2 diabetes (T2D). Importantly, we found that c-MAF inducing protein (CMIP) protects cells from exposure to FFAs by modulating Akt signaling and we validated the role of CMIP in human pancreatic beta cells. In sum, FALCON empowers the study of fundamental FFA biology and offers an integrative approach to identify much needed targets for diverse diseases associated with disordered FFA metabolism. Highlights: FALCON (Fatty Acid Library for Comprehensive ONtologies) enables multimodal profiling of 61 free fatty acids (FFAs) to reveal 5 FFA clusters with distinct biological effectsFALCON is applicable to many and diverse cell typesA subset of monounsaturated FAs (MUFAs) equally or more toxic than canonical lipotoxic saturated FAs (SFAs) leads to decreased membrane fluidityNew approach prioritizes genes that represent the combined effects of environmental (FFA) exposure and genetic risk for diseaseC-Maf inducing protein (CMIP) is identified as a suppressor of FFA-induced lipotoxicity via Akt-mediated signaling.

Single-nucleus cross-tissue molecular reference maps toward understanding disease gene function.

Understanding gene function and regulation in homeostasis and disease requires knowledge of the cellular and tissue contexts in which genes are expressed. Here, we applied four single-nucleus RNA sequencing methods to eight diverse, archived, frozen tissue types from 16 donors and 25 samples, generating a cross-tissue atlas of 209,126 nuclei profiles, which we integrated across tissues, donors, and laboratory methods with a conditional variational autoencoder. Using the resulting cross-tissue atlas, we highlight shared and tissue-specific features of tissue-resident cell populations; identify cell types that might contribute to neuromuscular, metabolic, and immune components of monogenic diseases and the biological processes involved in their pathology; and determine cell types and gene modules that might underlie disease mechanisms for complex traits analyzed by genome-wide association studies.

High-resolution Slide-seqV2 spatial transcriptomics enables discovery of disease-specific cell neighborhoods and pathways.

High-resolution spatial transcriptomics enables mapping of RNA expression directly from intact tissue sections; however, its utility for the elucidation of disease processes and therapeutically actionable pathways remains unexplored. We applied Slide-seqV2 to mouse and human kidneys, in healthy and distinct disease paradigms. First, we established the feasibility of Slide-seqV2 in tissue from nine distinct human kidneys, which revealed a cell neighborhood centered around a population of LYVE1+ macrophages. Second, in a mouse model of diabetic kidney disease, we detected changes in the cellular organization of the spatially restricted kidney filter and blood-flow-regulating apparatus. Third, in a mouse model of a toxic proteinopathy, we identified previously unknown, disease-specific cell neighborhoods centered around macrophages. In a spatially restricted subpopulation of epithelial cells, we discovered perturbations in 77 genes associated with the unfolded protein response. Our studies illustrate and experimentally validate the utility of Slide-seqV2 for the discovery of disease-specific cell neighborhoods.

WGS and RNA Studies Diagnose Noncoding DMD Variants in Males With High Creatine Kinase.

OBJECTIVE: To describe the diagnostic utility of whole-genome sequencing and RNA studies in boys with suspected dystrophinopathy, for whom multiplex ligation-dependent probe amplification and exomic parallel sequencing failed to yield a genetic diagnosis, and to use remnant normal DMD splicing in 3 families to define critical levels of wild-type dystrophin bridging clinical spectrums of Duchenne to myalgia. METHODS: Exome, genome, and/or muscle RNA sequencing was performed for 7 males with elevated creatine kinase. PCR of muscle-derived complementary DNA (cDNA) studied consequences for DMD premessenger RNA (pre-mRNA) splicing. Quantitative Western blot was used to determine levels of dystrophin, relative to control muscle. RESULTS: Splice-altering intronic single nucleotide variants or structural rearrangements in DMD were identified in all 7 families. Four individuals, with abnormal splicing causing a premature stop codon and nonsense-mediated decay, expressed remnant levels of normally spliced DMD mRNA. Quantitative Western blot enabled correlation of wild-type dystrophin and clinical severity, with 0%-5% dystrophin conferring a Duchenne phenotype, 10% ± 2% a Becker phenotype, and 15% ± 2% dystrophin associated with myalgia without manifesting weakness. CONCLUSIONS: Whole-genome sequencing relied heavily on RNA studies to identify DMD splice-altering variants. Short-read RNA sequencing was regularly confounded by the effectiveness of nonsense-mediated mRNA decay and low read depth of the giant DMD mRNA. PCR of muscle cDNA provided a simple, yet informative approach. Highly relevant to genetic therapies for dystrophinopathies, our data align strongly with previous studies of mutant dystrophin in Becker muscular dystrophy, with the collective conclusion that a fractional increase in levels of normal dystrophin between 5% and 20% is clinically significant.

Targeting a Braf/Mapk pathway rescues podocyte lipid peroxidation in CoQ-deficiency kidney disease.

Mutations affecting mitochondrial coenzyme Q (CoQ) biosynthesis lead to kidney failure due to selective loss of podocytes, essential cells of the kidney filter. Curiously, neighboring tubular epithelial cells are spared early in disease despite higher mitochondrial content. We sought to illuminate noncanonical, cell-specific roles for CoQ, independently of the electron transport chain (ETC). Here, we demonstrate that CoQ depletion caused by Pdss2 enzyme deficiency in podocytes results in perturbations in polyunsaturated fatty acid (PUFA) metabolism and the Braf/Mapk pathway rather than ETC dysfunction. Single-nucleus RNA-Seq from kidneys of Pdss2kd/kd mice with nephrotic syndrome and global CoQ deficiency identified a podocyte-specific perturbation of the Braf/Mapk pathway. Treatment with GDC-0879, a Braf/Mapk-targeting compound, ameliorated kidney disease in Pdss2kd/kd mice. Mechanistic studies in Pdss2-depleted podocytes revealed a previously unknown perturbation in PUFA metabolism that was confirmed in vivo. Gpx4, an enzyme that protects against PUFA-mediated lipid peroxidation, was elevated in disease and restored after GDC-0879 treatment. We demonstrate broader human disease relevance by uncovering patterns of GPX4 and Braf/Mapk pathway gene expression in tissue from patients with kidney diseases. Our studies reveal ETC-independent roles for CoQ in podocytes and point to Braf/Mapk as a candidate pathway for the treatment of kidney diseases.

Loss of sarcospan exacerbates pathology in mdx mice, but does not affect utrophin amelioration of disease.

The dystrophin-glycoprotein complex (DGC) is a membrane adhesion complex that provides structural stability at the sarcolemma by linking the myocyte's internal cytoskeleton and external extracellular matrix. In Duchenne muscular dystrophy (DMD), the absence of dystrophin leads to the loss of the DGC at the sarcolemma, resulting in sarcolemmal instability and progressive muscle damage. Utrophin (UTRN), an autosomal homolog of dystrophin, is upregulated in dystrophic muscle and partially compensates for the loss of dystrophin in muscle from patients with DMD. Here, we examine the interaction between Utr and sarcospan (SSPN), a small transmembrane protein that is a core component of both UTRN-glycoprotein complex (UGC) and DGC. We show that additional loss of SSPN causes an earlier onset of disease in dystrophin-deficient mdx mice by reducing the expression of the UGC at the sarcolemma. In order to further evaluate the role of SSPN in maintaining therapeutic levels of Utr at the sarcolemma, we tested the effect of Utr transgenic overexpression in mdx mice lacking SSPN (mdx:SSPN -/-:Utr-Tg). We found that overexpression of Utr restored SSPN to the sarcolemma in mdx muscle but that the ablation of SSPN in mdx muscle reduced Utr at the membrane. Nevertheless, Utr overexpression reduced central nucleation and improved grip strength in both lines. These findings demonstrate that high levels of Utr transgenic overexpression ameliorate the mdx phenotype independently of SSPN expression but that loss of SSPN may impair Utr-based mechanisms that rely on lower levels of Utr protein.

Principles of Spatial Transcriptomics Analysis: A Practical Walk-Through in Kidney Tissue.

Spatial transcriptomic technologies capture genome-wide readouts across biological tissue space. Moreover, recent advances in this technology, including Slide-seqV2, have achieved spatial transcriptomic data collection at a near-single cell resolution. To-date, a repertoire of computational tools has been developed to discern cell type classes given the transcriptomic profiles of tissue coordinates. Upon applying these tools, we can explore the spatial patterns of distinct cell types and characterize how genes are spatially expressed within different cell type contexts. The kidney is one organ whose function relies upon spatially defined structures consisting of distinct cellular makeup. Thus, the application of Slide-seqV2 to kidney tissue has enabled us to elucidate spatially characteristic cellular and genetic profiles at a scale that remains largely unexplored. Here, we review spatial transcriptomic technologies, as well as computational approaches for cell type mapping and spatial cell type and transcriptomic characterizations. We take kidney tissue as an example to demonstrate how the technologies are applied, while considering the nuances of this architecturally complex tissue.

Highly sensitive spatial transcriptomics at near-cellular resolution with Slide-seqV2.

Measurement of the location of molecules in tissues is essential for understanding tissue formation and function. Previously, we developed Slide-seq, a technology that enables transcriptome-wide detection of RNAs with a spatial resolution of 10 µm. Here we report Slide-seqV2, which combines improvements in library generation, bead synthesis and array indexing to reach an RNA capture efficiency ~50% that of single-cell RNA-seq data (~10-fold greater than Slide-seq), approaching the detection efficiency of droplet-based single-cell RNA-seq techniques. First, we leverage the detection efficiency of Slide-seqV2 to identify dendritically localized mRNAs in neurons of the mouse hippocampus. Second, we integrate the spatial information of Slide-seqV2 data with single-cell trajectory analysis tools to characterize the spatiotemporal development of the mouse neocortex, identifying underlying genetic programs that were poorly sampled with Slide-seq. The combination of near-cellular resolution and high transcript detection efficiency makes Slide-seqV2 useful across many experimental contexts.

HyPR-seq: Single-cell quantification of chosen RNAs via hybridization and sequencing of DNA probes.

Single-cell quantification of RNAs is important for understanding cellular heterogeneity and gene regulation, yet current approaches suffer from low sensitivity for individual transcripts, limiting their utility for many applications. Here we present Hybridization of Probes to RNA for sequencing (HyPR-seq), a method to sensitively quantify the expression of hundreds of chosen genes in single cells. HyPR-seq involves hybridizing DNA probes to RNA, distributing cells into nanoliter droplets, amplifying the probes with PCR, and sequencing the amplicons to quantify the expression of chosen genes. HyPR-seq achieves high sensitivity for individual transcripts, detects nonpolyadenylated and low-abundance transcripts, and can profile more than 100,000 single cells. We demonstrate how HyPR-seq can profile the effects of CRISPR perturbations in pooled screens, detect time-resolved changes in gene expression via measurements of gene introns, and detect rare transcripts and quantify cell-type frequencies in tissue using low-abundance marker genes. By directing sequencing power to genes of interest and sensitively quantifying individual transcripts, HyPR-seq reduces costs by up to 100-fold compared to whole-transcriptome single-cell RNA-sequencing, making HyPR-seq a powerful method for targeted RNA profiling in single cells.

The effect of LRRK2 loss-of-function variants in humans.

Human genetic variants predicted to cause loss-of-function of protein-coding genes (pLoF variants) provide natural in vivo models of human gene inactivation and can be valuable indicators of gene function and the potential toxicity of therapeutic inhibitors targeting these genes1,2. Gain-of-kinase-function variants in LRRK2 are known to significantly increase the risk of Parkinson's disease3,4, suggesting that inhibition of LRRK2 kinase activity is a promising therapeutic strategy. While preclinical studies in model organisms have raised some on-target toxicity concerns5-8, the biological consequences of LRRK2 inhibition have not been well characterized in humans. Here, we systematically analyze pLoF variants in LRRK2 observed across 141,456 individuals sequenced in the Genome Aggregation Database (gnomAD)9, 49,960 exome-sequenced individuals from the UK Biobank and over 4 million participants in the 23andMe genotyped dataset. After stringent variant curation, we identify 1,455 individuals with high-confidence pLoF variants in LRRK2. Experimental validation of three variants, combined with previous work10, confirmed reduced protein levels in 82.5% of our cohort. We show that heterozygous pLoF variants in LRRK2 reduce LRRK2 protein levels but that these are not strongly associated with any specific phenotype or disease state. Our results demonstrate the value of large-scale genomic databases and phenotyping of human loss-of-function carriers for target validation in drug discovery.

Disease-associated astrocytes in Alzheimer's disease and aging.

The role of non-neuronal cells in Alzheimer's disease progression has not been fully elucidated. Using single-nucleus RNA sequencing, we identified a population of disease-associated astrocytes in an Alzheimer's disease mouse model. These disease-associated astrocytes appeared at early disease stages and increased in abundance with disease progression. We discovered that similar astrocytes appeared in aged wild-type mice and in aging human brains, suggesting their linkage to genetic and age-related factors.

Efficient, continuous mutagenesis in human cells using a pseudo-random DNA editor.

Here we describe TRACE (T7 polymerase-driven continuous editing), a method that enables continuous, targeted mutagenesis in human cells using a cytidine deaminase fused to T7 RNA polymerase. TRACE induces high rates of mutagenesis over multiple cell generations in genes under the control of a T7 promoter integrated in the genome. We used TRACE in a MEK1 inhibitor-resistance screen, and identified functionally correlated mutations.

Single cell census of human kidney organoids shows reproducibility and diminished off-target cells after transplantation.

Human iPSC-derived kidney organoids have the potential to revolutionize discovery, but assessing their consistency and reproducibility across iPSC lines, and reducing the generation of off-target cells remain an open challenge. Here, we profile four human iPSC lines for a total of 450,118 single cells to show how organoid composition and development are comparable to human fetal and adult kidneys. Although cell classes are largely reproducible across time points, protocols, and replicates, we detect variability in cell proportions between different iPSC lines, largely due to off-target cells. To address this, we analyze organoids transplanted under the mouse kidney capsule and find diminished off-target cells. Our work shows how single cell RNA-seq (scRNA-seq) can score organoids for reproducibility, faithfulness and quality, that kidney organoids derived from different iPSC lines are comparable surrogates for human kidney, and that transplantation enhances their formation by diminishing off-target cells.

Identification of a Novel Deep Intronic Mutation in CAPN3 Presenting a Promising Target for Therapeutic Splice Modulation.

Calpainopathy, also known as limb girdle muscular dystrophy (LGMD) type 2A (LGMD2A) or LGMD R1 Calpain3-related, is one of the most common genetically characterized forms of limb-girdle muscular dystrophy with a wide range of phenotypic severity. We evaluated a consanguineous family with a clinical phenotype consistent with calpainopathy in whom conventional sequencing did not detect any mutations in the CAPN3 gene. Using whole exome sequencing paired with haplotype analysis, we identified a homozygous deep intronic single base pair deletion in CAPN3 (c.946-29delT). Familial segregation studies were consistent with recessive inheritance. Immunoblotting of muscle tissue from the patient showed complete absence of calpain 3. In silico analysis predicted the deletion to disrupt the branch point and subsequently alter splicing of exon 7. Studies of patient fibroblasts and muscle tissue confirmed altered splicing, resulting in an inclusion of a 389-bp intronic sequence upstream of exon 7, originating from a cryptic splice acceptor site in intron 6. This out-of-frame insertion results in a premature stop codon, leading to an apparent absence of protein likely due to degradation of the transcript via nonsense-mediated decay. We then designed phosphorodiamidate morpholino oligomers (PMOs) as splice modulators to block the new splice acceptor site. This approach successfully prevented the aberrant splicing - reverting the majority of the splice to the wildtype transcript. These results confirm the pathogenicity of this novel deep intronic mutation and provide a mutation-specific therapeutic strategy. Thus, deep intronic mutations in CAPN3 may be pathogenic and should be considered in the appropriate clinical setting. The identification of mutations which may be missed by traditional Sanger sequencing is essential as they may be excellent targets for individualized therapeutic strategies using RNA-directed splice modulation.

Stabilization of the cardiac sarcolemma by sarcospan rescues DMD-associated cardiomyopathy.

In the current preclinical study, we demonstrate the therapeutic potential of sarcospan (SSPN) overexpression to alleviate cardiomyopathy associated with Duchenne muscular dystrophy (DMD) utilizing dystrophin-deficient mdx mice with utrophin haploinsufficiency that more accurately represent the severe disease course of human DMD. SSPN interacts with dystrophin, the DMD disease gene product, and its autosomal paralog utrophin, which is upregulated in DMD as a partial compensatory mechanism. SSPN transgenic mice have enhanced abundance of fully glycosylated α-dystroglycan, which may further protect dystrophin-deficient cardiac membranes. Baseline echocardiography reveals SSPN improves systolic function and hypertrophic indices in mdx and mdx:utr-heterozygous mice. Assessment of SSPN transgenic mdx mice by hemodynamic pressure-volume methods highlights enhanced systolic performance compared to mdx controls. SSPN restores cardiac sarcolemma stability, the primary defect in DMD disease, reduces fibrotic response and improves contractile function. We demonstrate that SSPN ameliorates more advanced cardiac disease in the context of diminished sarcolemma expression of utrophin and ß1D integrin that mitigate disease severity and partially restores responsiveness to ß-adrenergic stimulation. Overall, our current and previous findings suggest SSPN overexpression in DMD mouse models positively impacts skeletal, pulmonary and cardiac performance by addressing the stability of proteins at the sarcolemma that protect the heart from injury, supporting SSPN and membrane stabilization as a therapeutic target for DMD.

A recurrent COL6A1 pseudoexon insertion causes muscular dystrophy and is effectively targeted by splice-correction therapies.

The clinical application of advanced next-generation sequencing technologies is increasingly uncovering novel classes of mutations that may serve as potential targets for precision medicine therapeutics. Here, we show that a deep intronic splice defect in the COL6A1 gene, originally discovered by applying muscle RNA sequencing in patients with clinical findings of collagen VI-related dystrophy (COL6-RD), inserts an in-frame pseudoexon into COL6A1 mRNA, encodes a mutant collagen α1(VI) protein that exerts a dominant-negative effect on collagen VI matrix assembly, and provides a unique opportunity for splice-correction approaches aimed at restoring normal gene expression. Using splice-modulating antisense oligomers, we efficiently skipped the pseudoexon in patient-derived fibroblast cultures and restored a wild-type matrix. Similarly, we used CRISPR/Cas9 to precisely delete an intronic sequence containing the pseudoexon and efficiently abolish its inclusion while preserving wild-type splicing. Considering that this splice defect is emerging as one of the single most frequent mutations in COL6-RD, the design of specific and effective splice-correction therapies offers a promising path for clinical translation.