Molecular hallmarks of heterochronic parabiosis at single-cell resolution.

The ability to slow or reverse biological ageing would have major implications for mitigating disease risk and maintaining vitality1. Although an increasing number of interventions show promise for rejuvenation2, their effectiveness on disparate cell types across the body and the molecular pathways susceptible to rejuvenation remain largely unexplored. Here we performed single-cell RNA sequencing on 20 organs to reveal cell-type-specific responses to young and aged blood in heterochronic parabiosis. Adipose mesenchymal stromal cells, haematopoietic stem cells and hepatocytes are among those cell types that are especially responsive. On the pathway level, young blood invokes new gene sets in addition to reversing established ageing patterns, with the global rescue of genes encoding electron transport chain subunits pinpointing a prominent role of mitochondrial function in parabiosis-mediated rejuvenation. We observed an almost universal loss of gene expression with age that is largely mimicked by parabiosis: aged blood reduces global gene expression, and young blood restores it in select cell types. Together, these data lay the groundwork for a systemic understanding of the interplay between blood-borne factors and cellular integrity.

Ageing hallmarks exhibit organ-specific temporal signatures.

Ageing is the single greatest cause of disease and death worldwide, and understanding the associated processes could vastly improve quality of life. Although major categories of ageing damage have been identified-such as altered intercellular communication, loss of proteostasis and eroded mitochondrial function1-these deleterious processes interact with extraordinary complexity within and between organs, and a comprehensive, whole-organism analysis of ageing dynamics has been lacking. Here we performed bulk RNA sequencing of 17 organs and plasma proteomics at 10 ages across the lifespan of Mus musculus, and integrated these findings with data from the accompanying Tabula Muris Senis2-or 'Mouse Ageing Cell Atlas'-which follows on from the original Tabula Muris3. We reveal linear and nonlinear shifts in gene expression during ageing, with the associated genes clustered in consistent trajectory groups with coherent biological functions-including extracellular matrix regulation, unfolded protein binding, mitochondrial function, and inflammatory and immune response. Notably, these gene sets show similar expression across tissues, differing only in the amplitude and the age of onset of expression. Widespread activation of immune cells is especially pronounced, and is first detectable in white adipose depots during middle age. Single-cell RNA sequencing confirms the accumulation of T cells and B cells in adipose tissue-including plasma cells that express immunoglobulin J-which also accrue concurrently across diverse organs. Finally, we show how gene expression shifts in distinct tissues are highly correlated with corresponding protein levels in plasma, thus potentially contributing to the ageing of the systemic circulation. Together, these data demonstrate a similar yet asynchronous inter- and intra-organ progression of ageing, providing a foundation from which to track systemic sources of declining health at old age.

An evolutionarily conserved ribosome-rescue pathway maintains epidermal homeostasis.

Ribosome-associated mRNA quality control mechanisms ensure the fidelity of protein translation1,2. Although these mechanisms have been extensively studied in yeast, little is known about their role in mammalian tissues, despite emerging evidence that stem cell fate is controlled by translational mechanisms3,4. One evolutionarily conserved component of the quality control machinery, Dom34 (in higher eukaryotes known as Pelota (Pelo)), rescues stalled ribosomes 5 . Here we show that Pelo is required for mammalian epidermal homeostasis. Conditional deletion of Pelo in mouse epidermal stem cells that express Lrig1 results in hyperproliferation and abnormal differentiation of these cells. By contrast, deletion of Pelo in Lgr5-expressing stem cells has no effect and deletion in Lgr6-expressing stem cells induces only a mild phenotype. Loss of Pelo results in accumulation of short ribosome footprints and global upregulation of translation, rather than affecting the expression of specific genes. Translational inhibition by rapamycin-mediated downregulation of mTOR (mechanistic target of rapamycin kinase) rescues the epidermal phenotype. Our study reveals that the ribosome-rescue machinery is important for mammalian tissue homeostasis and that it has specific effects on different stem cell populations.

MARS: discovering novel cell types across heterogeneous single-cell experiments.

Although tremendous effort has been put into cell-type annotation, identification of previously uncharacterized cell types in heterogeneous single-cell RNA-seq data remains a challenge. Here we present MARS, a meta-learning approach for identifying and annotating known as well as new cell types. MARS overcomes the heterogeneity of cell types by transferring latent cell representations across multiple datasets. MARS uses deep learning to learn a cell embedding function as well as a set of landmarks in the cell embedding space. The method has a unique ability to discover cell types that have never been seen before and annotate experiments that are as yet unannotated. We apply MARS to a large mouse cell atlas and show its ability to accurately identify cell types, even when it has never seen them before. Further, MARS automatically generates interpretable names for new cell types by probabilistically defining a cell type in the embedding space.

Inhibition of ß-catenin signalling in dermal fibroblasts enhances hair follicle regeneration during wound healing.

New hair follicles (HFs) do not form in adult mammalian skin unless epidermal Wnt signalling is activated genetically or within large wounds. To understand the postnatal loss of hair forming ability we monitored HF formation at small circular (2 mm) wound sites. At P2, new HFs formed in back skin, but HF formation was markedly decreased by P21. Neonatal tail also formed wound-associated HFs, albeit in smaller numbers. Postnatal loss of HF neogenesis did not correlate with wound closure rate but with a reduction in Lrig1-positive papillary fibroblasts in wounds. Comparative gene expression profiling of back and tail dermis at P1 and dorsal fibroblasts at P2 and P50 showed a correlation between loss of HF formation and decreased expression of genes associated with proliferation and Wnt/ß-catenin activity. Between P2 and P50, fibroblast density declined throughout the dermis and clones of fibroblasts became more dispersed. This correlated with a decline in fibroblasts expressing a TOPGFP reporter of Wnt activation. Surprisingly, between P2 and P50 there was no difference in fibroblast proliferation at the wound site but Wnt signalling was highly upregulated in healing dermis of P21 compared with P2 mice. Postnatal ß-catenin ablation in fibroblasts promoted HF regeneration in neonatal and adult mouse wounds, whereas ß-catenin activation reduced HF regeneration in neonatal wounds. Our data support a model whereby postnatal loss of hair forming ability in wounds reflects elevated dermal Wnt/ß-catenin activation in the wound bed, increasing the abundance of fibroblasts that are unable to induce HF formation.

Fibroblast state switching orchestrates dermal maturation and wound healing.

Murine dermis contains functionally and spatially distinct fibroblast lineages that cease to proliferate in early postnatal life. Here, we propose a model in which a negative feedback loop between extracellular matrix (ECM) deposition and fibroblast proliferation determines dermal architecture. Virtual-tissue simulations of our model faithfully recapitulate dermal maturation, predicting a loss of spatial segregation of fibroblast lineages and dictating that fibroblast migration is only required for wound healing. To test this, we performed in vivo live imaging of dermal fibroblasts, which revealed that homeostatic tissue architecture is achieved without active cell migration. In contrast, both fibroblast proliferation and migration are key determinants of tissue repair following wounding. The results show that tissue-scale coordination is driven by the interdependence of cell proliferation and ECM deposition, paving the way for identifying new therapeutic strategies to enhance skin regeneration.

Integration of Kinase and Calcium Signaling at the Level of Chromatin Underlies Inducible Gene Activation in T Cells.

TCR signaling pathways cooperate to activate the inducible transcription factors NF-κB, NFAT, and AP-1. In this study, using the calcium ionophore ionomycin and/or PMA on Jurkat T cells, we show that the gene expression program associated with activation of TCR signaling is closely related to specific chromatin landscapes. We find that calcium and kinase signaling cooperate to induce chromatin remodeling at ∼2100 chromatin regions, which demonstrate enriched binding motifs for inducible factors and correlate with target gene expression. We found that these regions typically function as inducible enhancers. Many of these elements contain composite NFAT/AP-1 sites, which typically support cooperative binding, thus further reinforcing the need for cooperation between calcium and kinase signaling in the activation of genes in T cells. In contrast, treatment with PMA or ionomycin alone induces chromatin remodeling at far fewer regions (∼600 and ∼350, respectively), which mostly represent a subset of those induced by costimulation. This suggests that the integration of TCR signaling largely occurs at the level of chromatin, which we propose plays a crucial role in regulating T cell activation.

How Mycobacterium tuberculosis builds a home: Single-cell analysis reveals M. tuberculosis ESX-1-mediated accumulation of anti-inflammatory macrophages in infected mouse lungs.

Mycobacterium tuberculosis (MTB) infects and replicates in lung mononuclear phagocytes (MNPs) with astounding ability to evade elimination. ESX-1, a type VII secretion system, acts as a virulence determinant that contributes to MTB's ability to survive within MNPs, but its effect on MNP recruitment and/or differentiation remains unknown. Here, using single-cell RNA sequencing, we studied the role of ESX-1 in MNP heterogeneity and response in mice and murine bone marrow-derived macrophages (BMDM). We found that ESX-1 is required for MTB to recruit diverse MNP subsets with high MTB burden. Further, MTB induces an anti-inflammatory signature in MNPs and BMDM in an ESX-1 dependent manner. Similarly, spatial transcriptomics revealed an upregulation of anti-inflammatory signals in MTB lesions, where monocyte-derived macrophages concentrate near MTB-infected cells. Together, our findings suggest that MTB ESX-1 mediates the recruitment and differentiation of anti-inflammatory MNPs, which MTB can infect and manipulate for survival.

Dysregulation of CD4+ and CD8+ resident memory T, myeloid, and stromal cells in steroid-experienced, checkpoint inhibitor colitis.

BACKGROUND: Colitis caused by checkpoint inhibitors (CPI) is frequent and is treated with empiric steroids, but CPI colitis mechanisms in steroid-experienced or refractory disease are unclear. METHODS: Using colon biopsies and blood from predominantly steroid-experienced CPI colitis patients, we performed multiplexed single-cell transcriptomics and proteomics to nominate contributing populations. RESULTS: CPI colitis biopsies showed enrichment of CD4+resident memory (RM) T cells in addition to CD8+ RM and cytotoxic CD8+ T cells. Matching T cell receptor (TCR) clonotypes suggested that both RMs are progenitors that yield cytotoxic effectors. Activated, CD38+ HLA-DR+ CD4+ RM and cytotoxic CD8+ T cells were enriched in steroid-experienced and a validation data set of steroid-naïve CPI colitis, underscoring their pathogenic potential across steroid exposure. Distinct from ulcerative colitis, CPI colitis exhibited perturbed stromal metabolism (NAD+, tryptophan) impacting epithelial survival and inflammation. Endothelial cells in CPI colitis after anti-TNF and anti-cytotoxic T-lymphocyte-associated antigen 4 (anti-CTLA-4) upregulated the integrin α4ß7 ligand molecular vascular addressin cell adhesion molecule 1 (MAdCAM-1), which may preferentially respond to vedolizumab (anti-α4ß7). CONCLUSIONS: These findings nominate CD4+ RM and MAdCAM-1+ endothelial cells for targeting in specific subsets of CPI colitis patients.

Defining and benchmarking open problems in single-cell analysis.

With the growing number of single-cell analysis tools, benchmarks are increasingly important to guide analysis and method development. However, a lack of standardisation and extensibility in current benchmarks limits their usability, longevity, and relevance to the community. We present Open Problems, a living, extensible, community-guided benchmarking platform including 10 current single-cell tasks that we envision will raise standards for the selection, evaluation, and development of methods in single-cell analysis.

Single-cell and spatial multi-omics highlight effects of anti-integrin therapy across cellular compartments in ulcerative colitis.

Ulcerative colitis (UC) is driven by immune and stromal subsets, culminating in epithelial injury. Vedolizumab (VDZ) is an anti-integrin antibody that is effective for treating UC. VDZ is known to inhibit lymphocyte trafficking to the intestine, but its broader effects on other cell subsets are less defined. To identify the inflammatory cells that contribute to colitis and are affected by VDZ, we performed single-cell transcriptomic and proteomic analyses of peripheral blood and colonic biopsies in healthy controls and patients with UC on VDZ or other therapies. Here we show that VDZ treatment is associated with alterations in circulating and tissue mononuclear phagocyte (MNP) subsets, along with modest shifts in lymphocytes. Spatial multi-omics of formalin-fixed biopsies demonstrates trends towards increased abundance and proximity of MNP and fibroblast subsets in active colitis. Spatial transcriptomics of archived specimens pre-treatment identifies epithelial-, MNP-, and fibroblast-enriched genes related to VDZ responsiveness, highlighting important roles for these subsets in UC.

Single-cell and spatial multi-omics highlight effects of anti-integrin therapy across cellular compartments in ulcerative colitis.

Ulcerative colitis (UC) is driven by immune and stromal subsets, culminating in epithelial injury. Vedolizumab (VDZ) is an anti-integrin antibody that is effective for treating UC. VDZ is known to inhibit lymphocyte trafficking to the intestine, but its broader effects on other cell subsets are less defined. To identify the inflammatory cells that contribute to colitis and are affected by VDZ, we perform single-cell transcriptomic and proteomic analyses of peripheral blood and colonic biopsies in healthy controls and patients with UC on VDZ or other therapies. Here we show that VDZ treatment is associated with alterations in circulating and tissue mononuclear phagocyte (MNP) subsets, along with modest shifts in lymphocytes. Spatial multi-omics of formalin-fixed biopsies demonstrates trends towards increased abundance and proximity of MNP and fibroblast subsets in active colitis. Spatial transcriptomics of archived specimens pre-treatment identifies epithelial-, MNP-, and fibroblast-enriched genes related to VDZ responsiveness, highlighting important roles for these subsets in UC.

scDiffCom: a tool for differential analysis of cell-cell interactions provides a mouse atlas of aging changes in intercellular communication.

Dysregulation of intercellular communication is a hallmark of aging. To better quantify and explore changes in intercellular communication, we present scDiffCom and scAgeCom. scDiffCom is an R package, relying on approximately 5,000 curated ligand-receptor interactions, that performs differential intercellular communication analysis between two conditions from single-cell transcriptomics data. Built upon scDiffCom, scAgeCom is an atlas of age-related cell-cell communication changes covering 23 mouse tissues from 58 single-cell RNA sequencing datasets from Tabula Muris Senis and the Calico murine aging cell atlas. It offers a comprehensive resource of tissue-specific and sex-specific aging dysregulations and highlights age-related intercellular communication changes widespread across the whole body, such as the upregulation of immune system processes and inflammation, the downregulation of developmental processes, angiogenesis and extracellular matrix organization and the deregulation of lipid metabolism. Our analysis emphasizes the relevance of the specific ligands, receptors and cell types regulating these processes. The atlas is available online ( https://scagecom.org ).

Single-nuclei characterization of pervasive transcriptional signatures across organs in response to COVID-19.

Background: Infection by coronavirus SARS-CoV2 is a severe and often deadly disease that has implications for the respiratory system and multiple organs across the human body. While the effects in the lung have been extensively studied, less is known about the impact COVID-19 has across other organs. Methods: Here, we contribute a single-nuclei RNA-sequencing atlas comprising six human organs across 20 autopsies where we analyzed the transcriptional changes due to COVID-19 in multiple cell types. The integration of data from multiple organs enabled the identification of systemic transcriptional changes. Results: Computational cross-organ analysis for endothelial cells and macrophages identified systemic transcriptional changes in these cell types in COVID-19 samples. In addition, analysis of gene modules showed enrichment of specific signaling pathways across multiple organs in COVID-19 autopsies. Conclusions: Altogether, the COVID Tissue Atlas enables the investigation of both cell type-specific and cross-organ transcriptional responses to COVID-19, providing insights into the molecular networks affected by the disease and highlighting novel potential targets for therapies and drug development. Funding: The Chan-Zuckerberg Initiative, The Chan-Zuckerberg Biohub.

A 2-Gene Host Signature for Improved Accuracy of COVID-19 Diagnosis Agnostic to Viral Variants.

The continued emergence of SARS-CoV-2 variants is one of several factors that may cause false-negative viral PCR test results. Such tests are also susceptible to false-positive results due to trace contamination from high viral titer samples. Host immune response markers provide an orthogonal indication of infection that can mitigate these concerns when combined with direct viral detection. Here, we leverage nasopharyngeal swab RNA-seq data from patients with COVID-19, other viral acute respiratory illnesses, and nonviral conditions (n = 318) to develop support vector machine classifiers that rely on a parsimonious 2-gene host signature to diagnose COVID-19. We find that optimal classifiers include an interferon-stimulated gene that is strongly induced in COVID-19 compared with nonviral conditions, such as IFI6, and a second immune-response gene that is more strongly induced in other viral infections, such as GBP5. The IFI6+GBP5 classifier achieves an area under the receiver operating characteristic curve (AUC) greater than 0.9 when evaluated on an independent RNA-seq cohort (n = 553). We further provide proof-of-concept demonstration that the classifier can be implemented in a clinically relevant RT-qPCR assay. Finally, we show that its performance is robust across common SARS-CoV-2 variants and is unaffected by cross-contamination, demonstrating its utility for improved accuracy of COVID-19 diagnostics. IMPORTANCE In this work, we study upper respiratory tract gene expression to develop and validate a 2-gene host-based COVID-19 diagnostic classifier and then demonstrate its implementation in a clinically practical qPCR assay. We find that the host classifier has utility for mitigating false-negative results, for example due to SARS-CoV-2 variants harboring mutations at primer target sites, and for mitigating false-positive viral PCR results due to laboratory cross-contamination. Both types of error carry serious consequences of either unrecognized viral transmission or unnecessary isolation and contact tracing. This work is directly relevant to the ongoing COVID-19 pandemic given the continued emergence of viral variants and the continued challenges of false-positive PCR assays. It also suggests the feasibility of pan-respiratory virus host-based diagnostics that would have value in congregate settings, such as hospitals and nursing homes, where unrecognized respiratory viral transmission is of particular concern.

The mitochondrial unfolded protein response regulates hippocampal neural stem cell aging.

Aging results in a decline in neural stem cells (NSCs), neurogenesis, and cognitive function, and evidence is emerging to demonstrate disrupted adult neurogenesis in the hippocampus of patients with several neurodegenerative disorders. Here, single-cell RNA sequencing of the dentate gyrus of young and old mice shows that the mitochondrial protein folding stress is prominent in activated NSCs/neural progenitors (NPCs) among the neurogenic niche, and it increases with aging accompanying dysregulated cell cycle and mitochondrial activity in activated NSCs/NPCs in the dentate gyrus. Increasing mitochondrial protein folding stress results in compromised NSC maintenance and reduced neurogenesis in the dentate gyrus, neural hyperactivity, and impaired cognitive function. Reducing mitochondrial protein folding stress in the dentate gyrus of old mice improves neurogenesis and cognitive function. These results establish the mitochondrial protein folding stress as a driver of NSC aging and suggest approaches to improve aging-associated cognitive decline.

Developmental origin and local signals cooperate to determine septal astrocyte identity.

Astrocyte specification during development is influenced by both intrinsic and extrinsic factors, but the precise contribution of each remains poorly understood. Here we show that septal astrocytes from Nkx2.1 and Zic4 expressing progenitor zones are allocated into non-overlapping domains of the medial (MS) and lateral septal nuclei (LS) respectively. Astrocytes in these areas exhibit distinctive molecular and morphological features tailored to the unique cellular and synaptic circuit environment of each nucleus. Using single-nucleus (sn) RNA sequencing, we trace the developmental trajectories of cells in the septum and find that neurons and astrocytes undergo region and developmental stage-specific local cell-cell interactions. We show that expression of the classic morphogens Sonic hedgehog (Shh) and Fibroblast growth factors (Fgfs) by MS and LS neurons respectively, functions to promote the molecular specification of local astrocytes in each region. Finally, using heterotopic cell transplantation, we show that both morphological and molecular specifications of septal astrocytes are highly dependent on the local microenvironment, regardless of developmental origins. Our data highlights the complex interplay between intrinsic and extrinsic factors shaping astrocyte identities and illustrates the importance of the local environment in determining astrocyte functional specialization.

Concordance of MERFISH spatial transcriptomics with bulk and single-cell RNA sequencing.

Spatial transcriptomics extends single-cell RNA sequencing (scRNA-seq) by providing spatial context for cell type identification and analysis. Imaging-based spatial technologies such as multiplexed error-robust fluorescence in situ hybridization (MERFISH) can achieve single-cell resolution, directly mapping single-cell identities to spatial positions. MERFISH produces a different data type than scRNA-seq, and a technical comparison between the two modalities is necessary to ascertain how to best integrate them. We performed MERFISH on the mouse liver and kidney and compared the resulting bulk and single-cell RNA statistics with those from the Tabula Muris Senis cell atlas and from two Visium datasets. MERFISH quantitatively reproduced the bulk RNA-seq and scRNA-seq results with improvements in overall dropout rates and sensitivity. Finally, we found that MERFISH independently resolved distinct cell types and spatial structure in both the liver and kidney. Computational integration with the Tabula Muris Senis atlas did not enhance these results. We conclude that MERFISH provides a quantitatively comparable method for single-cell gene expression and can identify cell types without the need for computational integration with scRNA-seq atlases.

Aging Fly Cell Atlas identifies exhaustive aging features at cellular resolution.

Aging is characterized by a decline in tissue function, but the underlying changes at cellular resolution across the organism remain unclear. Here, we present the Aging Fly Cell Atlas, a single-nucleus transcriptomic map of the whole aging Drosophila. We characterized 163 distinct cell types and performed an in-depth analysis of changes in tissue cell composition, gene expression, and cell identities. We further developed aging clock models to predict fly age and show that ribosomal gene expression is a conserved predictive factor for age. Combining all aging features, we find distinctive cell type-specific aging patterns. This atlas provides a valuable resource for studying fundamental principles of aging in complex organisms.

Indian Hedgehog release from TNF-activated renal epithelia drives local and remote organ fibrosis.

Progressive fibrosis is a feature of aging and chronic tissue injury in multiple organs, including the kidney and heart. Glioma-associated oncogene 1 expressing (Gli1+) cells are a major source of activated fibroblasts in multiple organs, but the links between injury, inflammation, and Gli1+ cell expansion and tissue fibrosis remain incompletely understood. We demonstrated that leukocyte-derived tumor necrosis factor (TNF) promoted Gli1+ cell proliferation and cardiorenal fibrosis through induction and release of Indian Hedgehog (IHH) from renal epithelial cells. Using single-cell-resolution transcriptomic analysis, we identified an "inflammatory" proximal tubular epithelial (iPT) population contributing to TNF- and nuclear factor κB (NF-κB)-induced IHH production in vivo. TNF-induced Ubiquitin D (Ubd) expression was observed in human proximal tubular cells in vitro and during murine and human renal disease and aging. Studies using pharmacological and conditional genetic ablation of TNF-induced IHH signaling revealed that IHH activated canonical Hedgehog signaling in Gli1+ cells, which led to their activation, proliferation, and fibrosis within the injured and aging kidney and heart. These changes were inhibited in mice by Ihh deletion in Pax8-expressing cells or by pharmacological blockade of TNF, NF-κB, or Gli1 signaling. Increased amounts of circulating IHH were associated with loss of renal function and higher rates of cardiovascular disease in patients with chronic kidney disease. Thus, IHH connects leukocyte activation to Gli1+ cell expansion and represents a potential target for therapies to inhibit inflammation-induced fibrosis.