Symbolic recording of signalling and cis-regulatory element activity to DNA.

Measurements of gene expression or signal transduction activity are conventionally performed using methods that require either the destruction or live imaging of a biological sample within the timeframe of interest. Here we demonstrate an alternative paradigm in which such biological activities are stably recorded to the genome. Enhancer-driven genomic recording of transcriptional activity in multiplex (ENGRAM) is based on the signal-dependent production of prime editing guide RNAs that mediate the insertion of signal-specific barcodes (symbols) into a genomically encoded recording unit. We show how this strategy can be used for multiplex recording of the cell-type-specific activities of dozens to hundreds of cis-regulatory elements with high fidelity, sensitivity and reproducibility. Leveraging signal transduction pathway-responsive cis-regulatory elements, we also demonstrate time- and concentration-dependent genomic recording of WNT, NF-κB and Tet-On activities. By coupling ENGRAM to sequential genome editing via DNA Typewriter1, we stably record information about the temporal dynamics of two orthogonal signalling pathways to genomic DNA. Finally we apply ENGRAM to integratively record the transient activity of nearly 100 transcription factor consensus motifs across daily windows spanning the differentiation of mouse embryonic stem cells into gastruloids, an in vitro model of early mammalian development. Although these are proof-of-concept experiments and much work remains to fully realize the possibilities, the symbolic recording of biological signals or states within cells, to the genome and over time, has broad potential to complement contemporary paradigms for how we make measurements in biological systems.

Meningioma transcriptomic landscape demonstrates novel subtypes with regional associated biology and patient outcome.

Meningiomas, although mostly benign, can be recurrent and fatal. World Health Organization (WHO) grading of the tumor does not always identify high-risk meningioma, and better characterizations of their aggressive biology are needed. To approach this problem, we combined 13 bulk RNA sequencing (RNA-seq) datasets to create a dimension-reduced reference landscape of 1,298 meningiomas. The clinical and genomic metadata effectively correlated with landscape regions, which led to the identification of meningioma subtypes with specific biological signatures. The time to recurrence also correlated with the map location. Further, we developed an algorithm that maps new patients onto this landscape, where the nearest neighbors predict outcome. This study highlights the utility of combining bulk transcriptomic datasets to visualize the complexity of tumor populations. Further, we provide an interactive tool for understanding the disease and predicting patient outcomes. This resource is accessible via the online tool Oncoscape, where the scientific community can explore the meningioma landscape.

Multi-condition and multi-modal temporal profile inference during mouse embryonic development.

The emergence of single-cell time-series datasets enables modeling of changes in various types of cellular profiles over time. However, due to the disruptive nature of single-cell measurements, it is impossible to capture the full temporal trajectory of a particular cell. Furthermore, single-cell profiles can be collected at mismatched time points across different conditions (e.g., sex, batch, disease) and data modalities (e.g., scRNA-seq, scATAC-seq), which makes modeling challenging. Here we propose a joint modeling framework, Sunbear, for integrating multi-condition and multi-modal single-cell profiles across time. Sunbear can be used to impute single-cell temporal profile changes, align multi-dataset and multi-modal profiles across time, and extrapolate single-cell profiles in a missing modality. We applied Sunbear to reveal sex-biased transcription during mouse embryonic development and predict dynamic relationships between epigenetic priming and transcription for cells in which multi-modal profiles are unavailable. Sunbear thus enables the projection of single-cell time-series snapshots to multi-modal and multi-condition views of cellular trajectories.

A single-cell time-lapse of mouse prenatal development from gastrula to birth.

The house mouse (Mus musculus) is an exceptional model system, combining genetic tractability with close evolutionary affinity to humans1,2. Mouse gestation lasts only 3 weeks, during which the genome orchestrates the astonishing transformation of a single-cell zygote into a free-living pup composed of more than 500 million cells. Here, to establish a global framework for exploring mammalian development, we applied optimized single-cell combinatorial indexing3 to profile the transcriptional states of 12.4 million nuclei from 83 embryos, precisely staged at 2- to 6-hour intervals spanning late gastrulation (embryonic day 8) to birth (postnatal day 0). From these data, we annotate hundreds of cell types and explore the ontogenesis of the posterior embryo during somitogenesis and of kidney, mesenchyme, retina and early neurons. We leverage the temporal resolution and sampling depth of these whole-embryo snapshots, together with published data4-8 from earlier timepoints, to construct a rooted tree of cell-type relationships that spans the entirety of prenatal development, from zygote to birth. Throughout this tree, we systematically nominate genes encoding transcription factors and other proteins as candidate drivers of the in vivo differentiation of hundreds of cell types. Remarkably, the most marked temporal shifts in cell states are observed within one hour of birth and presumably underlie the massive physiological adaptations that must accompany the successful transition of a mammalian fetus to life outside the womb.

Induction and in silico staging of human gastruloids with neural tube, segmented somites & advanced cell types.

Embryonic organoids are emerging as powerful models for studying early mammalian development. For example, stem cell-derived 'gastruloids' form elongating structures containing all three germ layers1-4. However, although elongated, human gastruloids do not morphologically resemble post-implantation embryos. Here we show that a specific, discontinuous regimen of retinoic acid (RA) robustly induces human gastruloids with embryo-like morphological structures, including a neural tube and segmented somites. Single cell RNA-seq (sc-RNA-seq) further reveals that these human 'RA-gastruloids' contain more advanced cell types than conventional gastruloids, including neural crest cells, renal progenitor cells, skeletal muscle cells, and, rarely, neural progenitor cells. We apply a new approach to computationally stage human RA-gastruloids relative to somite-resolved mouse embryos, early human embryos and other gastruloid models, and find that the developmental stage of human RA-gastruloids is comparable to that of E9.5 mouse embryos, although some cell types show greater or lesser progression. We chemically perturb WNT and BMP signaling in human RA-gastruloids and find that these signaling pathways regulate somite patterning and neural tube length, respectively, while genetic perturbation of the transcription factors PAX3 and TBX6 markedly compromises the formation of neural crest and somites/renal cells, respectively. Human RA-gastruloids complement other embryonic organoids in serving as a simple, robust and screenable model for decoding early human embryogenesis.

Single-cell, whole-embryo phenotyping of mammalian developmental disorders.

Mouse models are a critical tool for studying human diseases, particularly developmental disorders1. However, conventional approaches for phenotyping may fail to detect subtle defects throughout the developing mouse2. Here we set out to establish single-cell RNA sequencing of the whole embryo as a scalable platform for the systematic phenotyping of mouse genetic models. We applied combinatorial indexing-based single-cell RNA sequencing3 to profile 101 embryos of 22 mutant and 4 wild-type genotypes at embryonic day 13.5, altogether profiling more than 1.6 million nuclei. The 22 mutants represent a range of anticipated phenotypic severities, from established multisystem disorders to deletions of individual regulatory regions4,5. We developed and applied several analytical frameworks for detecting differences in composition and/or gene expression across 52 cell types or trajectories. Some mutants exhibit changes in dozens of trajectories whereas others exhibit changes in only a few cell types. We also identify differences between widely used wild-type strains, compare phenotyping of gain- versus loss-of-function mutants and characterize deletions of topological associating domain boundaries. Notably, some changes are shared among mutants, suggesting that developmental pleiotropy might be 'decomposable' through further scaling of this approach. Overall, our findings show how single-cell profiling of whole embryos can enable the systematic molecular and cellular phenotypic characterization of mouse mutants with unprecedented breadth and resolution.

A single-cell transcriptional timelapse of mouse embryonic development, from gastrula to pup.

The house mouse, Mus musculus, is an exceptional model system, combining genetic tractability with close homology to human biology. Gestation in mouse development lasts just under three weeks, a period during which its genome orchestrates the astonishing transformation of a single cell zygote into a free-living pup composed of >500 million cells. Towards a global framework for exploring mammalian development, we applied single cell combinatorial indexing (sci-*) to profile the transcriptional states of 12.4 million nuclei from 83 precisely staged embryos spanning late gastrulation (embryonic day 8 or E8) to birth (postnatal day 0 or P0), with 2-hr temporal resolution during somitogenesis, 6-hr resolution through to birth, and 20-min resolution during the immediate postpartum period. From these data (E8 to P0), we annotate dozens of trajectories and hundreds of cell types and perform deeper analyses of the unfolding of the posterior embryo during somitogenesis as well as the ontogenesis of the kidney, mesenchyme, retina, and early neurons. Finally, we leverage the depth and temporal resolution of these whole embryo snapshots, together with other published data, to construct and curate a rooted tree of cell type relationships that spans mouse development from zygote to pup. Throughout this tree, we systematically nominate sets of transcription factors (TFs) and other genes as candidate drivers of the in vivo differentiation of hundreds of mammalian cell types. Remarkably, the most dramatic shifts in transcriptional state are observed in a restricted set of cell types in the hours immediately following birth, and presumably underlie the massive changes in physiology that must accompany the successful transition of a placental mammal to extrauterine life.

Optimized single-nucleus transcriptional profiling by combinatorial indexing.

Single-cell combinatorial indexing RNA sequencing (sci-RNA-seq) is a powerful method for recovering gene expression data from an exponentially scalable number of individual cells or nuclei. However, sci-RNA-seq is a complex protocol that has historically exhibited variable performance on different tissues, as well as lower sensitivity than alternative methods. Here, we report a simplified, optimized version of the sci-RNA-seq protocol with three rounds of split-pool indexing that is faster, more robust and more sensitive and has a higher yield than the original protocol, with reagent costs on the order of 1 cent per cell or less. The total hands-on time from nuclei isolation to final library preparation takes 2-3 d, depending on the number of samples sharing the experiment. The improvements also allow RNA profiling from tissues rich in RNases like older mouse embryos or adult tissues that were problematic for the original method. We showcase the optimized protocol via whole-organism analysis of an E16.5 mouse embryo, profiling ~380,000 nuclei in a single experiment. Finally, we introduce a 'Tiny-Sci' protocol for experiments in which input material is very limited.

Embryo model completes gastrulation to neurulation and organogenesis.

Embryonic stem (ES) cells can undergo many aspects of mammalian embryogenesis in vitro1-5, but their developmental potential is substantially extended by interactions with extraembryonic stem cells, including trophoblast stem (TS) cells, extraembryonic endoderm stem (XEN) cells and inducible XEN (iXEN) cells6-11. Here we assembled stem cell-derived embryos in vitro from mouse ES cells, TS cells and iXEN cells and showed that they recapitulate the development of whole natural mouse embryo in utero up to day 8.5 post-fertilization. Our embryo model displays headfolds with defined forebrain and midbrain regions and develops a beating heart-like structure, a trunk comprising a neural tube and somites, a tail bud containing neuromesodermal progenitors, a gut tube, and primordial germ cells. This complete embryo model develops within an extraembryonic yolk sac that initiates blood island development. Notably, we demonstrate that the neurulating embryo model assembled from Pax6-knockout ES cells aggregated with wild-type TS cells and iXEN cells recapitulates the ventral domain expansion of the neural tube that occurs in natural, ubiquitous Pax6-knockout embryos. Thus, these complete embryoids are a powerful in vitro model for dissecting the roles of diverse cell lineages and genes in development. Our results demonstrate the self-organization ability of ES cells and two types of extraembryonic stem cells to reconstitute mammalian development through and beyond gastrulation to neurulation and early organogenesis.

Systematic reconstruction of cellular trajectories across mouse embryogenesis.

Mammalian embryogenesis is characterized by rapid cellular proliferation and diversification. Within a few weeks, a single-cell zygote gives rise to millions of cells expressing a panoply of molecular programs. Although intensively studied, a comprehensive delineation of the major cellular trajectories that comprise mammalian development in vivo remains elusive. Here, we set out to integrate several single-cell RNA-sequencing (scRNA-seq) datasets that collectively span mouse gastrulation and organogenesis, supplemented with new profiling of ~150,000 nuclei from approximately embryonic day 8.5 (E8.5) embryos staged in one-somite increments. Overall, we define cell states at each of 19 successive stages spanning E3.5 to E13.5 and heuristically connect them to their pseudoancestors and pseudodescendants. Although constructed through automated procedures, the resulting directed acyclic graph (TOME (trajectories of mammalian embryogenesis)) is largely consistent with our contemporary understanding of mammalian development. We leverage TOME to systematically nominate transcription factors (TFs) as candidate regulators of each cell type's specification, as well as 'cell-type homologs' across vertebrate evolution.

Genome-wide association studies identify the role of caspase-9 in kidney disease.

Genome-wide association studies (GWAS) have identified hundreds of genetic risk regions for kidney dysfunction [estimated glomerular filtration rate (eGFR)]; however, the causal genes, cell types, and pathways are poorly understood. Integration of GWAS and human kidney expression of quantitative trait analysis using Bayesian colocations, transcriptome-wide association studies, and summary-based Mendelian randomization studies prioritized caspase-9 (CASP9) as a kidney disease risk gene. Human kidney single-cell epigenetic and immunostaining studies indicated kidney tubule cells as a disease-causing cell type. Mice with genetic deletion or pharmacological inhibition of CASP9 showed lower apoptosis while having improved mitophagy, resulting in dampened activation of cytosolic nucleotide sensing pathways (cGAS-STING), reduction of inflammation, and protection from acute kidney disease or renal fibrosis. In summary, here, we prioritized CASP9 as an eGFR GWAS target gene and demonstrated the causal role of CASP9 in kidney disease development via improving mitophagy and lowering inflammation and apoptosis.

Mapping the genetic architecture of human traits to cell types in the kidney identifies mechanisms of disease and potential treatments.

The functional interpretation of genome-wide association studies (GWAS) is challenging due to the cell-type-dependent influences of genetic variants. Here, we generated comprehensive maps of expression quantitative trait loci (eQTLs) for 659 microdissected human kidney samples and identified cell-type-eQTLs by mapping interactions between cell type abundances and genotypes. By partitioning heritability using stratified linkage disequilibrium score regression to integrate GWAS with single-cell RNA sequencing and single-nucleus assay for transposase-accessible chromatin with high-throughput sequencing data, we prioritized proximal tubules for kidney function and endothelial cells and distal tubule segments for blood pressure pathogenesis. Bayesian colocalization analysis nominated more than 200 genes for kidney function and hypertension. Our study clarifies the mechanism of commonly used antihypertensive and renal-protective drugs and identifies drug repurposing opportunities for kidney disease.

Transcriptome-wide association analysis identifies DACH1 as a kidney disease risk gene that contributes to fibrosis.

Genome-wide association studies (GWAS) for kidney function identified hundreds of risk regions; however, the causal variants, target genes, cell types, and disease mechanisms remain poorly understood. Here, we performed transcriptome-wide association studies (TWAS), summary Mendelian randomization, and MetaXcan to identify genes whose expression mediates the genotype effect on the phenotype. Our analyses identified Dachshund homolog 1 (DACH1), a cell-fate determination factor. GWAS risk variant was associated with lower DACH1 expression in human kidney tubules. Human and mouse kidney single-cell open chromatin data (snATAC-Seq) prioritized estimated glomerular filtration rate (eGFR) GWAS variants located on an intronic regulatory region in distal convoluted tubule cells. CRISPR-Cas9-mediated gene editing confirmed the role of risk variants in regulating DACH1 expression. Mice with tubule-specific Dach1 deletion developed more severe renal fibrosis both in folic acid and diabetic kidney injury models. Mice with tubule-specific Dach1 overexpression were protected from folic acid nephropathy. Single-cell RNA sequencing, chromatin immunoprecipitation, and functional analysis indicated that DACH1 controls the expression of cell cycle and myeloid chemotactic factors, contributing to macrophage infiltration and fibrosis development. In summary, integration of GWAS, TWAS, single-cell epigenome, expression analyses, gene editing, and functional validation in different mouse kidney disease models identified DACH1 as a kidney disease risk gene.

Kidney disease genetic risk variants alter lysosomal beta-mannosidase (MANBA) expression and disease severity.

More than 800 million people in the world suffer from chronic kidney disease (CKD). Genome-wide association studies (GWAS) have identified hundreds of loci where genetic variants are associated with kidney function; however, causal genes and pathways for CKD remain unknown. Here, we performed integration of kidney function GWAS and human kidney-specific expression quantitative trait analysis and identified that the expression of beta-mannosidase (MANBA) was lower in kidneys of subjects with CKD risk genotype. We also show an increased incidence of renal failure in subjects with rare heterozygous loss-of-function coding variants in MANBA using phenome-wide association analysis of 40,963 subjects with exome sequencing data. MANBA is a lysosomal gene highly expressed in kidney tubule cells. Deep phenotyping revealed structural and functional lysosomal alterations in human kidneys from subjects with CKD risk alleles and mice with genetic deletion of Manba Manba heterozygous and knockout mice developed more severe kidney fibrosis when subjected to toxic injury induced by cisplatin or folic acid. Manba loss altered multiple pathways, including endocytosis and autophagy. In the absence of Manba, toxic acute tubule injury induced inflammasome activation and fibrosis. Together, these results illustrate the convergence of common noncoding and rare coding variants in MANBA in kidney disease development and demonstrate the role of the endolysosomal system in kidney disease development.

Systematic integrated analysis of genetic and epigenetic variation in diabetic kidney disease.

Poor metabolic control and host genetic predisposition are critical for diabetic kidney disease (DKD) development. The epigenome integrates information from sequence variations and metabolic alterations. Here, we performed a genome-wide methylome association analysis in 500 subjects with DKD from the Chronic Renal Insufficiency Cohort for DKD phenotypes, including glycemic control, albuminuria, kidney function, and kidney function decline. We show distinct methylation patterns associated with each phenotype. We define methylation variations that are associated with underlying nucleotide variations (methylation quantitative trait loci) and show that underlying genetic variations are important drivers of methylation changes. We implemented Bayesian multitrait colocalization analysis (moloc) and summary data-based Mendelian randomization to systematically annotate genomic regions that show association with kidney function, methylation, and gene expression. We prioritized 40 loci, where methylation and gene-expression changes likely mediate the genotype effect on kidney disease development. Functional annotation suggested the role of inflammation, specifically, apoptotic cell clearance and complement activation in kidney disease development. Our study defines methylation changes associated with DKD phenotypes, the key role of underlying genetic variations driving methylation variations, and prioritizes methylome and gene-expression changes that likely mediate the genotype effect on kidney disease pathogenesis.

Loss of IL-27Rα Results in Enhanced Tubulointerstitial Fibrosis Associated with Elevated Th17 Responses.

Clinical and experimental studies have established that immune cells such as alternatively activated (M2) macrophages and Th17 cells play a role in the progression of chronic kidney disease, but the endogenous pathways that limit these processes are not well understood. The cytokine IL-27 has been shown to limit immune-mediated pathology in other systems by effects on these cell types, but this has not been thoroughly investigated in the kidney. Unilateral ureteral obstruction was performed on wild-type and IL-27Rα-/- mice. After 2 wk, kidneys were extracted, and the degree of injury was measured by hydroxyproline assay and quantification of neutrophil gelatinase-associated lipocalin mRNA. Immune cell infiltrate was evaluated by immunohistochemistry and flow cytometry. An anti-IL-17A mAb was subsequently administered to IL-27Rα-/- mice every 2 d from day of surgery with evaluation as described after 2 wk. After unilateral ureteral obstruction, IL-27 deficiency resulted in increased tissue injury and collagen deposition associated with higher levels of chemokine mRNA and increased numbers of M2 macrophages. Loss of the IL-27Rα led to increased infiltration of activated CD4+ T cells that coproduced IL-17A and TNF-α, and blockade of IL-17A partially ameliorated kidney injury. Patients with chronic kidney disease had elevated serum levels of IL-27 and IL-17A, whereas expression of transcripts for the IL-27RA and the IL-17RA in the tubular epithelial cells of patients with renal fibrosis correlated with disease severity. These data suggest that endogenous IL-27 acts at several points in the inflammatory cascade to limit the magnitude of immune-mediated damage to the kidney.

Phenome-wide association analysis suggests the APOL1 linked disease spectrum primarily drives kidney-specific pathways.

The relationship between commonly occurring genetic variants (G1 and G2) in the APOL1 gene in African Americans and different disease traits, such as kidney disease, cardiovascular disease, and pre-eclampsia, remains the subject of controversy. Here we took a genotype-first approach, a phenome-wide association study, to define the spectrum of phenotypes associated with APOL1 high-risk variants in 1,837 African American participants of Penn Medicine Biobank and 4,742 African American participants of Vanderbilt BioVU. In the Penn Medicine Biobank, outpatient creatinine measurement-based estimated glomerular filtration rate and multivariable regression models were used to evaluate the association between high-risk APOL1 status and renal outcomes. In meta-analysis of both cohorts, the strongest phenome-wide association study associations were for the high-risk APOL1 variants and diagnoses codes were highly significant for "kidney dialysis" (odds ratio 3.75) and "end stage kidney disease" (odds ratio 3.42). A number of phenotypes were associated with APOL1 high-risk genotypes in an analysis adjusted only for demographic variables. However, no associations were detected with non-renal phenotypes after controlling for chronic/end stage kidney disease status. Using calculated estimated glomerular filtration rate -based phenotype analysis in the Penn Medicine Biobank, APOL1 high-risk status was associated with prevalent chronic/end stage kidney disease /kidney transplant (odds ratio 2.27, 95% confidence interval 1.67-3.08). In high-risk participants, the estimated glomerular filtration rate was 15.4 mL/min/1.73m2; significantly lower than in low-risk participants. Thus, although APOL1 high-risk variants are associated with a range of phenotypes, the risks for other associated phenotypes appear much lower and in our dataset are driven by a primary effect on renal disease.

Target genes, variants, tissues and transcriptional pathways influencing human serum urate levels.

Elevated serum urate levels cause gout and correlate with cardiometabolic diseases via poorly understood mechanisms. We performed a trans-ancestry genome-wide association study of serum urate in 457,690 individuals, identifying 183 loci (147 previously unknown) that improve the prediction of gout in an independent cohort of 334,880 individuals. Serum urate showed significant genetic correlations with many cardiometabolic traits, with genetic causality analyses supporting a substantial role for pleiotropy. Enrichment analysis, fine-mapping of urate-associated loci and colocalization with gene expression in 47 tissues implicated the kidney and liver as the main target organs and prioritized potentially causal genes and variants, including the transcriptional master regulators in the liver and kidney, HNF1A and HNF4A. Experimental validation showed that HNF4A transactivated the promoter of ABCG2, encoding a major urate transporter, in kidney cells, and that HNF4A p.Thr139Ile is a functional variant. Transcriptional coregulation within and across organs may be a general mechanism underlying the observed pleiotropy between urate and cardiometabolic traits.

Genome-Wide Association Study of Diabetic Kidney Disease Highlights Biology Involved in Glomerular Basement Membrane Collagen.

BACKGROUND: Although diabetic kidney disease demonstrates both familial clustering and single nucleotide polymorphism heritability, the specific genetic factors influencing risk remain largely unknown. METHODS: To identify genetic variants predisposing to diabetic kidney disease, we performed genome-wide association study (GWAS) analyses. Through collaboration with the Diabetes Nephropathy Collaborative Research Initiative, we assembled a large collection of type 1 diabetes cohorts with harmonized diabetic kidney disease phenotypes. We used a spectrum of ten diabetic kidney disease definitions based on albuminuria and renal function. RESULTS: Our GWAS meta-analysis included association results for up to 19,406 individuals of European descent with type 1 diabetes. We identified 16 genome-wide significant risk loci. The variant with the strongest association (rs55703767) is a common missense mutation in the collagen type IV alpha 3 chain (COL4A3) gene, which encodes a major structural component of the glomerular basement membrane (GBM). Mutations in COL4A3 are implicated in heritable nephropathies, including the progressive inherited nephropathy Alport syndrome. The rs55703767 minor allele (Asp326Tyr) is protective against several definitions of diabetic kidney disease, including albuminuria and ESKD, and demonstrated a significant association with GBM width; protective allele carriers had thinner GBM before any signs of kidney disease, and its effect was dependent on glycemia. Three other loci are in or near genes with known or suggestive involvement in this condition (BMP7) or renal biology (COLEC11 and DDR1). CONCLUSIONS: The 16 diabetic kidney disease-associated loci may provide novel insights into the pathogenesis of this condition and help identify potential biologic targets for prevention and treatment.