Spatial proteomics of human diabetic kidney disease, from health to class III.

AIMS/HYPOTHESIS: Diabetic kidney disease (DKD) is the leading cause of chronic and end-stage kidney disease in the USA and worldwide. Animal models have taught us much about DKD mechanisms, but translation of this knowledge into treatments for human disease has been slowed by the lag in our molecular understanding of human DKD. METHODS: Using our Spatial TissuE Proteomics (STEP) pipeline (comprising curated human kidney tissues, multiplexed immunofluorescence and powerful analysis tools), we imaged and analysed the expression of 21 proteins in 23 tissue sections from individuals with diabetes and healthy kidneys (n=5), compared to those with DKDIIA, IIA-B and IIB (n=2 each) and DKDIII (n=1). RESULTS: These analyses revealed the existence of 11 cellular clusters (kidney compartments/cell types): podocytes, glomerular endothelial cells, proximal tubules, distal nephron, peritubular capillaries, blood vessels (endothelial cells and vascular smooth muscle cells), macrophages, myeloid cells, other CD45+ inflammatory cells, basement membrane and the interstitium. DKD progression was associated with co-localised increases in inflammatory cells and collagen IV deposition, with concomitant loss of native proteins of each nephron segment. Cell-type frequency and neighbourhood analyses highlighted a significant increase in inflammatory cells and their adjacency to tubular and αSMA+ (α-smooth muscle actin-positive) cells in DKD. Finally, DKD progression showed marked regional variability within single tissue sections, as well as inter-individual variability within each DKD class. CONCLUSIONS/INTERPRETATION: Using the STEP pipeline, we found alterations in protein expression, cellular phenotypic composition and microenvironment structure with DKD progression, demonstrating the power of this pipeline to reveal the pathophysiology of human DKD.

PEPSI: Polarity measurements from spatial proteomics imaging suggest immune cell engagement.

Subcellular protein localization is important for understanding functional states of cells, but measuring and quantifying this information can be difficult and typically requires high-resolution microscopy. In this work, we develop a metric to define surface protein polarity from immunofluorescence (IF) imaging data and use it to identify distinct immune cell states within tumor microenvironments. We apply this metric to characterize over two million cells across 600 patient samples and find that cells identified as having polar expression exhibit characteristics relating to tumor-immune cell engagement. Additionally, we show that incorporating these polarity-defined cell subtypes improves the performance of deep learning models trained to predict patient survival outcomes. This method provides a first look at using subcellular protein expression patterns to phenotype immune cell functional states with applications to precision medicine.

7-UP: Generating in silico CODEX from a small set of immunofluorescence markers.

Multiplex immunofluorescence (mIF) assays multiple protein biomarkers on a single tissue section. Recently, high-plex CODEX (co-detection by indexing) systems enable simultaneous imaging of 40+ protein biomarkers, unlocking more detailed molecular phenotyping, leading to richer insights into cellular interactions and disease. However, high-plex data can be slower and more costly to collect, limiting its applications, especially in clinical settings. We propose a machine learning framework, 7-UP, that can computationally generate in silico 40-plex CODEX at single-cell resolution from a standard 7-plex mIF panel by leveraging cellular morphology. We demonstrate the usefulness of the imputed biomarkers in accurately classifying cell types and predicting patient survival outcomes. Furthermore, 7-UP's imputations generalize well across samples from different clinical sites and cancer types. 7-UP opens the possibility of in silico CODEX, making insights from high-plex mIF more widely available.

CasKAS: direct profiling of genome-wide dCas9 and Cas9 specificity using ssDNA mapping.

Detecting and mitigating off-target activity is critical to the practical application of CRISPR-mediated genome and epigenome editing. While numerous methods have been developed to map Cas9 binding specificity genome-wide, they are generally time-consuming and/or expensive, and not applicable to catalytically dead CRISPR enzymes. We have developed CasKAS, a rapid, inexpensive, and facile assay for identifying off-target CRISPR enzyme binding and cleavage by chemically mapping the unwound single-stranded DNA structures formed upon binding of a sgRNA-loaded Cas9 protein. We demonstrate this method in both in vitro and in vivo contexts.

Nociceptor neurons affect cancer immunosurveillance.

Solid tumours are innervated by nerve fibres that arise from the autonomic and sensory peripheral nervous systems1-5. Whether the neo-innervation of tumours by pain-initiating sensory neurons affects cancer immunosurveillance remains unclear. Here we show that melanoma cells interact with nociceptor neurons, leading to increases in their neurite outgrowth, responsiveness to noxious ligands and neuropeptide release. Calcitonin gene-related peptide (CGRP)-one such nociceptor-produced neuropeptide-directly increases the exhaustion of cytotoxic CD8+ T cells, which limits their capacity to eliminate melanoma. Genetic ablation of the TRPV1 lineage, local pharmacological silencing of nociceptors and antagonism of the CGRP receptor RAMP1 all reduced the exhaustion of tumour-infiltrating leukocytes and decreased the growth of tumours, nearly tripling the survival rate of mice that were inoculated with B16F10 melanoma cells. Conversely, CD8+ T cell exhaustion was rescued in sensory-neuron-depleted mice that were treated with local recombinant CGRP. As compared with wild-type CD8+ T cells, Ramp1-/- CD8+ T cells were protected against exhaustion when co-transplanted into tumour-bearing Rag1-deficient mice. Single-cell RNA sequencing of biopsies from patients with melanoma revealed that intratumoral RAMP1-expressing CD8+ T cells were more exhausted than their RAMP1-negative counterparts, whereas overexpression of RAMP1 correlated with a poorer clinical prognosis. Overall, our results suggest that reducing the release of CGRP from tumour-innervating nociceptors could be a strategy to improve anti-tumour immunity by eliminating the immunomodulatory effects of CGRP on cytotoxic CD8+ T cells.

Graph deep learning for the characterization of tumour microenvironments from spatial protein profiles in tissue specimens.

Multiplexed immunofluorescence imaging allows the multidimensional molecular profiling of cellular environments at subcellular resolution. However, identifying and characterizing disease-relevant microenvironments from these rich datasets is challenging. Here we show that a graph neural network that leverages spatial protein profiles in tissue specimens to model tumour microenvironments as local subgraphs captures distinctive cellular interactions associated with differential clinical outcomes. We applied this spatial cellular-graph strategy to specimens of human head-and-neck and colorectal cancers assayed with 40-plex immunofluorescence imaging to identify spatial motifs associated with cancer recurrence and with patient survival after treatment. The graph deep learning model was substantially more accurate in predicting patient outcomes than deep learning approaches that model spatial data on the basis of the local composition of cell types, and it generated insights into the effect of the spatial compartmentalization of tumour cells and granulocytes on patient prognosis. Local graphs may also aid in the analysis of disease-relevant motifs in histology samples characterized via spatial transcriptomics and other -omics techniques.

Single-Cell Transcriptomic Census of Endothelial Changes Induced by Matrix Stiffness and the Association with Atherosclerosis.

Vascular endothelial cell (EC) plasticity plays a critical role in the progression of atherosclerosis by giving rise to mesenchymal phenotypes in the plaque lesion. Despite the evidence for arterial stiffening as a major contributor to atherosclerosis, the complex interplay among atherogenic stimuli in vivo has hindered attempts to determine the effects of extracellular matrix (ECM) stiffness on endothelial-mesenchymal transition (EndMT). To study the regulatory effects of ECM stiffness on EndMT, an in vitro model is developed in which human coronary artery ECs are cultured on physiological or pathological stiffness substrates. Leveraging single-cell RNA sequencing, cell clusters with mesenchymal transcriptional features are identified to be more prevalent on pathological substrates than physiological substrates. Trajectory inference analyses reveal a novel mesenchymal-to-endothelial reverse transition, which is blocked by pathological stiffness substrates, in addition to the expected EndMT trajectory. ECs pushed to a mesenchymal character by pathological stiffness substrates are enriched in transcriptional signatures of atherosclerotic ECs from human and murine plaques. This study characterizes at single-cell resolution the transcriptional programs that underpin EC plasticity in both physiological or pathological milieus, and thus serves as a valuable resource for more precisely defining EndMT and the transcriptional programs contributing to atherosclerosis.

Chromatin and gene-regulatory dynamics of the developing human cerebral cortex at single-cell resolution.

Genetic perturbations of cortical development can lead to neurodevelopmental disease, including autism spectrum disorder (ASD). To identify genomic regions crucial to corticogenesis, we mapped the activity of gene-regulatory elements generating a single-cell atlas of gene expression and chromatin accessibility both independently and jointly. This revealed waves of gene regulation by key transcription factors (TFs) across a nearly continuous differentiation trajectory, distinguished the expression programs of glial lineages, and identified lineage-determining TFs that exhibited strong correlation between linked gene-regulatory elements and expression levels. These highly connected genes adopted an active chromatin state in early differentiating cells, consistent with lineage commitment. Base-pair-resolution neural network models identified strong cell-type-specific enrichment of noncoding mutations predicted to be disruptive in a cohort of ASD individuals and identified frequently disrupted TF binding sites. This approach illustrates how cell-type-specific mapping can provide insights into the programs governing human development and disease.

Transcription-dependent domain-scale three-dimensional genome organization in the dinoflagellate Breviolum minutum.

Dinoflagellate chromosomes represent a unique evolutionary experiment, as they exist in a permanently condensed, liquid crystalline state; are not packaged by histones; and contain genes organized into tandem gene arrays, with minimal transcriptional regulation. We analyze the three-dimensional genome of Breviolum minutum, and find large topological domains (dinoflagellate topologically associating domains, which we term 'dinoTADs') without chromatin loops, which are demarcated by convergent gene array boundaries. Transcriptional inhibition disrupts dinoTADs, implicating transcription-induced supercoiling as the primary topological force in dinoflagellates.

Diverse lncRNA mechanisms in brain development and disease.

Long noncoding RNAs (lncRNAs) are a diverse and pervasive class of genes. Recent studies in the mammalian brain have uncovered several novel mechanisms. LncRNA loci are often located in proximity to developmental transcriptional factors. The lncRNA product may act like a transcription factor to control distantly located genes, or in other instances, the lncRNA loci contain DNA regulatory elements that act locally on neighboring genes. Circular RNAs are covalently closed single-stranded RNAs that can control neuronal function by acting as microRNA sponges and additional mechanisms. LncRNAs can also engage in target-directed microRNA degradation to shape the pool of microRNAs and translation. Thus, diverse mechanisms allow lncRNAs to act in the nucleus and cytoplasm to control neuronal fate and function.

Chromatin accessibility dynamics in a model of human forebrain development.

Forebrain development is characterized by highly synchronized cellular processes, which, if perturbed, can cause disease. To chart the regulatory activity underlying these events, we generated a map of accessible chromatin in human three-dimensional forebrain organoids. To capture corticogenesis, we sampled glial and neuronal lineages from dorsal or ventral forebrain organoids over 20 months in vitro. Active chromatin regions identified in human primary brain tissue were observed in organoids at different developmental stages. We used this resource to map genetic risk for disease and to explore evolutionary conservation. Moreover, we integrated chromatin accessibility with transcriptomics to identify putative enhancer-gene linkages and transcription factors that regulate human corticogenesis. Overall, this platform brings insights into gene-regulatory dynamics at previously inaccessible stages of human forebrain development, including signatures of neuropsychiatric disorders.

Mitigation of off-target toxicity in CRISPR-Cas9 screens for essential non-coding elements.

Pooled CRISPR-Cas9 screens are a powerful method for functionally characterizing regulatory elements in the non-coding genome, but off-target effects in these experiments have not been systematically evaluated. Here, we investigate Cas9, dCas9, and CRISPRi/a off-target activity in screens for essential regulatory elements. The sgRNAs with the largest effects in genome-scale screens for essential CTCF loop anchors in K562 cells were not single guide RNAs (sgRNAs) that disrupted gene expression near the on-target CTCF anchor. Rather, these sgRNAs had high off-target activity that, while only weakly correlated with absolute off-target site number, could be predicted by the recently developed GuideScan specificity score. Screens conducted in parallel with CRISPRi/a, which do not induce double-stranded DNA breaks, revealed that a distinct set of off-targets also cause strong confounding fitness effects with these epigenome-editing tools. Promisingly, filtering of CRISPRi libraries using GuideScan specificity scores removed these confounded sgRNAs and enabled identification of essential regulatory elements.

An improved ATAC-seq protocol reduces background and enables interrogation of frozen tissues.

We present Omni-ATAC, an improved ATAC-seq protocol for chromatin accessibility profiling that works across multiple applications with substantial improvement of signal-to-background ratio and information content. The Omni-ATAC protocol generates chromatin accessibility profiles from archival frozen tissue samples and 50-µm sections, revealing the activities of disease-associated DNA elements in distinct human brain structures. The Omni-ATAC protocol enables the interrogation of personal regulomes in tissue context and translational studies.

Genome-scale transcriptional activation by an engineered CRISPR-Cas9 complex.

Systematic interrogation of gene function requires the ability to perturb gene expression in a robust and generalizable manner. Here we describe structure-guided engineering of a CRISPR-Cas9 complex to mediate efficient transcriptional activation at endogenous genomic loci. We used these engineered Cas9 activation complexes to investigate single-guide RNA (sgRNA) targeting rules for effective transcriptional activation, to demonstrate multiplexed activation of ten genes simultaneously, and to upregulate long intergenic non-coding RNA (lincRNA) transcripts. We also synthesized a library consisting of 70,290 guides targeting all human RefSeq coding isoforms to screen for genes that, upon activation, confer resistance to a BRAF inhibitor. The top hits included genes previously shown to be able to confer resistance, and novel candidates were validated using individual sgRNA and complementary DNA overexpression. A gene expression signature based on the top screening hits correlated with markers of BRAF inhibitor resistance in cell lines and patient-derived samples. These results collectively demonstrate the potential of Cas9-based activators as a powerful genetic perturbation technology.

Genome editing using Cas9 nickases.

The RNA-guided, sequence-specific endonuclease Cas9 has been widely adopted as genome engineering tool due to its efficiency and ease of use. Derived from the microbial CRISPR (clustered regularly interspaced short palindromic repeat) type II adaptive immune system, Cas9 has now been successfully engineered for genome editing applications in a variety of animal and plant species. To reduce potential off-target mutagenesis by wild-type Cas9, homology- and structure-guided mutagenesis of Streptococcus pyogenes Cas9 catalytic domains has produced "nicking" enzymes (Cas9n) capable of inducing single-strand nicks rather than double-strand breaks. Since nicks are generally repaired with high fidelity in eukaryotic cells, Cas9n can be leveraged to mediate highly specific genome editing, either via nonhomologous end-joining or homology-directed repair. Here we describe the preparation, testing, and application of Cas9n reagents for precision mammalian genome engineering.

Genome-wide binding of the CRISPR endonuclease Cas9 in mammalian cells.

Bacterial type II CRISPR-Cas9 systems have been widely adapted for RNA-guided genome editing and transcription regulation in eukaryotic cells, yet their in vivo target specificity is poorly understood. Here we mapped genome-wide binding sites of a catalytically inactive Cas9 (dCas9) from Streptococcus pyogenes loaded with single guide RNAs (sgRNAs) in mouse embryonic stem cells (mESCs). Each of the four sgRNAs we tested targets dCas9 to between tens and thousands of genomic sites, frequently characterized by a 5-nucleotide seed region in the sgRNA and an NGG protospacer adjacent motif (PAM). Chromatin inaccessibility decreases dCas9 binding to other sites with matching seed sequences; thus 70% of off-target sites are associated with genes. Targeted sequencing of 295 dCas9 binding sites in mESCs transfected with catalytically active Cas9 identified only one site mutated above background levels. We propose a two-state model for Cas9 binding and cleavage, in which a seed match triggers binding but extensive pairing with target DNA is required for cleavage.

Double nicking by RNA-guided CRISPR Cas9 for enhanced genome editing specificity.

Targeted genome editing technologies have enabled a broad range of research and medical applications. The Cas9 nuclease from the microbial CRISPR-Cas system is targeted to specific genomic loci by a 20 nt guide sequence, which can tolerate certain mismatches to the DNA target and thereby promote undesired off-target mutagenesis. Here, we describe an approach that combines a Cas9 nickase mutant with paired guide RNAs to introduce targeted double-strand breaks. Because individual nicks in the genome are repaired with high fidelity, simultaneous nicking via appropriately offset guide RNAs is required for double-stranded breaks and extends the number of specifically recognized bases for target cleavage. We demonstrate that using paired nicking can reduce off-target activity by 50- to 1,500-fold in cell lines and to facilitate gene knockout in mouse zygotes without sacrificing on-target cleavage efficiency. This versatile strategy enables a wide variety of genome editing applications that require high specificity.