Best practices for single-cell analysis across modalities.

Recent advances in single-cell technologies have enabled high-throughput molecular profiling of cells across modalities and locations. Single-cell transcriptomics data can now be complemented by chromatin accessibility, surface protein expression, adaptive immune receptor repertoire profiling and spatial information. The increasing availability of single-cell data across modalities has motivated the development of novel computational methods to help analysts derive biological insights. As the field grows, it becomes increasingly difficult to navigate the vast landscape of tools and analysis steps. Here, we summarize independent benchmarking studies of unimodal and multimodal single-cell analysis across modalities to suggest comprehensive best-practice workflows for the most common analysis steps. Where independent benchmarks are not available, we review and contrast popular methods. Our article serves as an entry point for novices in the field of single-cell (multi-)omic analysis and guides advanced users to the most recent best practices.

CD201+ fascia progenitors choreograph injury repair.

Optimal tissue recovery and organismal survival are achieved by spatiotemporal tuning of tissue inflammation, contraction and scar formation1. Here we identify a multipotent fibroblast progenitor marked by CD201 expression in the fascia, the deepest connective tissue layer of the skin. Using skin injury models in mice, single-cell transcriptomics and genetic lineage tracing, ablation and gene deletion models, we demonstrate that CD201+ progenitors control the pace of wound healing by generating multiple specialized cell types, from proinflammatory fibroblasts to myofibroblasts, in a spatiotemporally tuned sequence. We identified retinoic acid and hypoxia signalling as the entry checkpoints into proinflammatory and myofibroblast states. Modulating CD201+ progenitor differentiation impaired the spatiotemporal appearances of fibroblasts and chronically delayed wound healing. The discovery of proinflammatory and myofibroblast progenitors and their differentiation pathways provide a new roadmap to understand and clinically treat impaired wound healing.

Spatial single-cell mass spectrometry defines zonation of the hepatocyte proteome.

Single-cell proteomics by mass spectrometry is emerging as a powerful and unbiased method for the characterization of biological heterogeneity. So far, it has been limited to cultured cells, whereas an expansion of the method to complex tissues would greatly enhance biological insights. Here we describe single-cell Deep Visual Proteomics (scDVP), a technology that integrates high-content imaging, laser microdissection and multiplexed mass spectrometry. scDVP resolves the context-dependent, spatial proteome of murine hepatocytes at a current depth of 1,700 proteins from a cell slice. Half of the proteome was differentially regulated in a spatial manner, with protein levels changing dramatically in proximity to the central vein. We applied machine learning to proteome classes and images, which subsequently inferred the spatial proteome from imaging data alone. scDVP is applicable to healthy and diseased tissues and complements other spatial proteomics and spatial omics technologies.

CD151 Maintains Endolysosomal Protein Quality to Inhibit Vascular Inflammation.

BACKGROUND: Tetraspanin CD151 is highly expressed in endothelia and reinforces cell adhesion, but its role in vascular inflammation remains largely unknown. METHODS: In vitro molecular and cellular biological analyses on genetically modified endothelial cells, in vivo vascular biological analyses on genetically engineered mouse models, and in silico systems biology and bioinformatics analyses on CD151-related events. RESULTS: Endothelial ablation of Cd151 leads to pulmonary and cardiac inflammation, severe sepsis, and perilous COVID-19, and endothelial CD151 becomes downregulated in inflammation. Mechanistically, CD151 restrains endothelial release of proinflammatory molecules for less leukocyte infiltration. At the subcellular level, CD151 determines the integrity of multivesicular bodies/lysosomes and confines the production of exosomes that carry cytokines such as ANGPT2 (angiopoietin-2) and proteases such as cathepsin-D. At the molecular level, CD151 docks VCP (valosin-containing protein)/p97, which controls protein quality via mediating deubiquitination for proteolytic degradation, onto endolysosomes to facilitate VCP/p97 function. At the endolysosome membrane, CD151 links VCP/p97 to (1) IFITM3 (interferon-induced transmembrane protein 3), which regulates multivesicular body functions, to restrain IFITM3-mediated exosomal sorting, and (2) V-ATPase, which dictates endolysosome pH, to support functional assembly of V-ATPase. CONCLUSIONS: Distinct from its canonical function in strengthening cell adhesion at cell surface, CD151 maintains endolysosome function by sustaining VCP/p97-mediated protein unfolding and turnover. By supporting protein quality control and protein degradation, CD151 prevents proteins from (1) buildup in endolysosomes and (2) discharge through exosomes, to limit vascular inflammation. Also, our study conceptualizes that balance between degradation and discharge of proteins in endothelial cells determines vascular information. Thus, the IFITM3/V-ATPase-tetraspanin-VCP/p97 complexes on endolysosome, as a protein quality control and inflammation-inhibitory machinery, could be beneficial for therapeutic intervention against vascular inflammation.

Inhibition of LTßR signalling activates WNT-induced regeneration in lung.

Lymphotoxin ß-receptor (LTßR) signalling promotes lymphoid neogenesis and the development of tertiary lymphoid structures1,2, which are associated with severe chronic inflammatory diseases that span several organ systems3-6. How LTßR signalling drives chronic tissue damage particularly in the lung, the mechanism(s) that regulate this process, and whether LTßR blockade might be of therapeutic value have remained unclear. Here we demonstrate increased expression of LTßR ligands in adaptive and innate immune cells, enhanced non-canonical NF-κB signalling, and enriched LTßR target gene expression in lung epithelial cells from patients with smoking-associated chronic obstructive pulmonary disease (COPD) and from mice chronically exposed to cigarette smoke. Therapeutic inhibition of LTßR signalling in young and aged mice disrupted smoking-related inducible bronchus-associated lymphoid tissue, induced regeneration of lung tissue, and reverted airway fibrosis and systemic muscle wasting. Mechanistically, blockade of LTßR signalling dampened epithelial non-canonical activation of NF-κB, reduced TGFß signalling in airways, and induced regeneration by preventing epithelial cell death and activating WNT/ß-catenin signalling in alveolar epithelial progenitor cells. These findings suggest that inhibition of LTßR signalling represents a viable therapeutic option that combines prevention of tertiary lymphoid structures1 and inhibition of apoptosis with tissue-regenerative strategies.

CellRank for directed single-cell fate mapping.

Computational trajectory inference enables the reconstruction of cell state dynamics from single-cell RNA sequencing experiments. However, trajectory inference requires that the direction of a biological process is known, largely limiting its application to differentiating systems in normal development. Here, we present CellRank ( https://cellrank.org ) for single-cell fate mapping in diverse scenarios, including regeneration, reprogramming and disease, for which direction is unknown. Our approach combines the robustness of trajectory inference with directional information from RNA velocity, taking into account the gradual and stochastic nature of cellular fate decisions, as well as uncertainty in velocity vectors. On pancreas development data, CellRank automatically detects initial, intermediate and terminal populations, predicts fate potentials and visualizes continuous gene expression trends along individual lineages. Applied to lineage-traced cellular reprogramming data, predicted fate probabilities correctly recover reprogramming outcomes. CellRank also predicts a new dedifferentiation trajectory during postinjury lung regeneration, including previously unknown intermediate cell states, which we confirm experimentally.

COPD basal cells are primed towards secretory to multiciliated cell imbalance driving increased resilience to environmental stressors.

INTRODUCTION: Environmental pollutants injure the mucociliary elevator, thereby provoking disease progression in chronic obstructive pulmonary disease (COPD). Epithelial resilience mechanisms to environmental nanoparticles in health and disease are poorly characterised. METHODS: We delineated the impact of prevalent pollutants such as carbon and zinc oxide nanoparticles, on cellular function and progeny in primary human bronchial epithelial cells (pHBECs) from end-stage COPD (COPD-IV, n=4), early disease (COPD-II, n=3) and pulmonary healthy individuals (n=4). After nanoparticle exposure of pHBECs at air-liquid interface, cell cultures were characterised by functional assays, transcriptome and protein analysis, complemented by single-cell analysis in serial samples of pHBEC cultures focusing on basal cell differentiation. RESULTS: COPD-IV was characterised by a prosecretory phenotype (twofold increase in MUC5AC+) at the expense of the multiciliated epithelium (threefold reduction in Ac-Tub+), resulting in an increased resilience towards particle-induced cell damage (fivefold reduction in transepithelial electrical resistance), as exemplified by environmentally abundant doses of zinc oxide nanoparticles. Exposure of COPD-II cultures to cigarette smoke extract provoked the COPD-IV characteristic, prosecretory phenotype. Time-resolved single-cell transcriptomics revealed an underlying COPD-IV unique basal cell state characterised by a twofold increase in KRT5+ (P=0.018) and LAMB3+ (P=0.050) expression, as well as a significant activation of Wnt-specific (P=0.014) and Notch-specific (P=0.021) genes, especially in precursors of suprabasal and secretory cells. CONCLUSION: We identified COPD stage-specific gene alterations in basal cells that affect the cellular composition of the bronchial elevator and may control disease-specific epithelial resilience mechanisms in response to environmental nanoparticles. The identified phenomena likely inform treatment and prevention strategies.

Sfrp1 inhibits lung fibroblast invasion during transition to injury-induced myofibroblasts.

BACKGROUND: Fibroblast-to-myofibroblast conversion is a major driver of tissue remodelling in organ fibrosis. Distinct lineages of fibroblasts support homeostatic tissue niche functions, yet their specific activation states and phenotypic trajectories during injury and repair have remained unclear. METHODS: We combined spatial transcriptomics, multiplexed immunostainings, longitudinal single-cell RNA-sequencing and genetic lineage tracing to study fibroblast fates during mouse lung regeneration. Our findings were validated in idiopathic pulmonary fibrosis patient tissues in situ as well as in cell differentiation and invasion assays using patient lung fibroblasts. Cell differentiation and invasion assays established a function of SFRP1 in regulating human lung fibroblast invasion in response to transforming growth factor (TGF)ß1. MEASUREMENTS AND MAIN RESULTS: We discovered a transitional fibroblast state characterised by high Sfrp1 expression, derived from both Tcf21-Cre lineage positive and negative cells. Sfrp1 + cells appeared early after injury in peribronchiolar, adventitial and alveolar locations and preceded the emergence of myofibroblasts. We identified lineage-specific paracrine signals and inferred converging transcriptional trajectories towards Sfrp1 + transitional fibroblasts and Cthrc1 + myofibroblasts. TGFß1 downregulated SFRP1 in noninvasive transitional cells and induced their switch to an invasive CTHRC1+ myofibroblast identity. Finally, using loss-of-function studies we showed that SFRP1 modulates TGFß1-induced fibroblast invasion and RHOA pathway activity. CONCLUSIONS: Our study reveals the convergence of spatially and transcriptionally distinct fibroblast lineages into transcriptionally uniform myofibroblasts and identifies SFRP1 as a modulator of TGFß1-driven fibroblast phenotypes in fibrogenesis. These findings are relevant in the context of therapeutic interventions that aim at limiting or reversing fibroblast foci formation.

RIPK1 kinase-dependent inflammation and cell death contribute to the pathogenesis of COPD.

BACKGROUND: Receptor-interacting protein kinase 1 (RIPK1) is a key mediator of regulated cell death (including apoptosis and necroptosis) and inflammation, both drivers of COPD pathogenesis. We aimed to define the contribution of RIPK1 kinase-dependent cell death and inflammation in the pathogenesis of COPD. METHODS: We assessed RIPK1 expression in single-cell RNA sequencing (RNA-seq) data from human and mouse lungs, and validated RIPK1 levels in lung tissue of COPD patients via immunohistochemistry. Next, we assessed the consequences of genetic and pharmacological inhibition of RIPK1 kinase activity in experimental COPD, using Ripk1 S25D/S25D kinase-deficient mice and the RIPK1 kinase inhibitor GSK'547. RESULTS: RIPK1 expression increased in alveolar type 1 (AT1), AT2, ciliated and neuroendocrine cells in human COPD. RIPK1 protein levels were significantly increased in airway epithelium of COPD patients compared with never-smokers and smokers without airflow limitation. In mice, exposure to cigarette smoke (CS) increased Ripk1 expression similarly in AT2 cells, and further in alveolar macrophages and T-cells. Genetic and/or pharmacological inhibition of RIPK1 kinase activity significantly attenuated airway inflammation upon acute and subacute CS exposure, as well as airway remodelling, emphysema, and apoptotic and necroptotic cell death upon chronic CS exposure. Similarly, pharmacological RIPK1 kinase inhibition significantly attenuated elastase-induced emphysema and lung function decline. Finally, RNA-seq on lung tissue of CS-exposed mice revealed downregulation of cell death and inflammatory pathways upon pharmacological RIPK1 kinase inhibition. CONCLUSIONS: RIPK1 kinase inhibition is protective in experimental models of COPD and may represent a novel promising therapeutic approach.

FGFR2b signalling restricts lineage-flexible alveolar progenitors during mouse lung development and converges in mature alveolar type 2 cells.

The specification, characterization, and fate of alveolar type 1 and type 2 (AT1 and AT2) progenitors during embryonic lung development are poorly defined. Current models of distal epithelial lineage formation fail to capture the heterogeneity and dynamic contribution of progenitor pools present during early development. Furthermore, few studies explore the pathways involved in alveolar progenitor specification and fate. In this paper, we build upon our previously published work on the regulation of airway epithelial progenitors by fibroblast growth factor receptor 2b (FGFR2b) signalling during early (E12.5) and mid (E14.5) pseudoglandular stage lung development. Our results suggest that a significant proportion of AT2 and AT1 progenitors are lineage-flexible during late pseudoglandular stage development, and that lineage commitment is regulated in part by FGFR2b signalling. We have characterized a set of direct FGFR2b targets at E16.5 which are likely involved in alveolar lineage formation. These signature genes converge on a subpopulation of AT2 cells later in development and are downregulated in AT2 cells transitioning to the AT1 lineage during repair after injury in adults. Our findings highlight the extensive heterogeneity of pneumocytes by elucidating the role of FGFR2b signalling in these cells during early airway epithelial lineage formation, as well as during repair after injury.

The discovAIR project: a roadmap towards the Human Lung Cell Atlas.

The Human Cell Atlas (HCA) consortium aims to establish an atlas of all organs in the healthy human body at single-cell resolution to increase our understanding of basic biological processes that govern development, physiology and anatomy, and to accelerate diagnosis and treatment of disease. The Lung Biological Network of the HCA aims to generate the Human Lung Cell Atlas as a reference for the cellular repertoire, molecular cell states and phenotypes, and cell-cell interactions that characterise normal lung homeostasis in healthy lung tissue. Such a reference atlas of the healthy human lung will facilitate mapping the changes in the cellular landscape in disease. The discovAIR project is one of six pilot actions for the HCA funded by the European Commission in the context of the H2020 framework programme. discovAIR aims to establish the first draft of an integrated Human Lung Cell Atlas, combining single-cell transcriptional and epigenetic profiling with spatially resolving techniques on matched tissue samples, as well as including a number of chronic and infectious diseases of the lung. The integrated Human Lung Cell Atlas will be available as a resource for the wider respiratory community, including basic and translational scientists, clinical medicine, and the private sector, as well as for patients with lung disease and the interested lay public. We anticipate that the Human Lung Cell Atlas will be the founding stone for a more detailed understanding of the pathogenesis of lung diseases, guiding the design of novel diagnostics and preventive or curative interventions.

IL-37 regulates allergic inflammation by counterbalancing pro-inflammatory IL-1 and IL-33.

BACKGROUND: Children with asthma have impaired production of interleukin (IL) 37; in mice, IL-37 reduces hallmarks of experimental allergic asthma (EAA). However, it remains unclear how IL-37 exerts its inhibitory properties in asthma. This study aimed to identify the mechanism(s) by which IL-37 controls allergic inflammation. METHODS: IL-37 target cells were identified by single-cell RNA-seq of IL-1R5 and IL-1R8. Airway tissues were isolated by laser-capture microdissection and examined by microarray-based gene expression analysis. Mononuclear cells (MNC) and airway epithelial cells (AECs) were isolated and stimulated with allergen, IL-1ß, or IL-33 together with recombinant human (rh) IL-37. Wild-type, IL-1R1- and IL-33-deficient mice with EAA were treated with rhIL-37. IL-1ß, IL-33, and IL-37 levels were determined in sputum and nasal secretions from adult asthma patients without glucocorticoid therapy. RESULTS: IL-37 target cells included AECs, T cells, and dendritic cells. In mice with EAA, rhIL-37 led to differential expression of >90 genes induced by IL-1ß and IL-33. rhIL-37 reduced production of Th2 cytokines in allergen-activated MNCs from wild-type but not from IL-1R1-deficient mice and inhibited IL-33-induced Th2 cytokine release. Furthermore, rhIL-37 attenuated IL-1ß- and IL-33-induced pro-inflammatory mediator expression in murine AEC cultures. In contrast to wild-type mice, hIL-37 had no effect on EAA in IL-1R1- or IL-33-deficient mice. We also observed that expression/production ratios of both IL-1ß and IL-33 to IL-37 were dramatically increased in asthma patients compared to healthy controls. CONCLUSION: IL-37 downregulates allergic airway inflammation by counterbalancing the disease-amplifying effects of IL-1ß and IL-33.

Rivaroxaban Reduces Arterial Thrombosis by Inhibition of FXa-Driven Platelet Activation via Protease Activated Receptor-1.

RATIONALE: A reduced rate of myocardial infarction has been reported in patients with atrial fibrillation treated with FXa (factor Xa) inhibitors including rivaroxaban compared with vitamin K antagonists. At the same time, low-dose rivaroxaban has been shown to reduce mortality and atherothrombotic events in patients with coronary artery disease. Yet, the mechanisms underlying this reduction remain unknown. OBJECTIVE: In this study, we hypothesized that rivaroxaban's antithrombotic potential is linked to a hitherto unknown rivaroxaban effect that impacts on platelet reactivity and arterial thrombosis. METHODS AND RESULTS: In this study, we identified FXa as potent, direct agonist of the PAR-1 (protease-activated receptor 1), leading to platelet activation and thrombus formation, which can be inhibited by rivaroxaban. We found that rivaroxaban reduced arterial thrombus stability in a mouse model of arterial thrombosis using intravital microscopy. For in vitro studies, atrial fibrillation patients on permanent rivaroxaban treatment for stroke prevention, respective controls, and patients with new-onset atrial fibrillation before and after first intake of rivaroxaban (time series analysis) were recruited. Platelet aggregation responses, as well as thrombus formation under arterial flow conditions on collagen and atherosclerotic plaque material, were attenuated by rivaroxaban. We show that rivaroxaban's antiplatelet effect is plasma dependent but independent of thrombin and rivaroxaban's anticoagulatory capacity. CONCLUSIONS: Here, we identified FXa as potent platelet agonist that acts through PAR-1. Therefore, rivaroxaban exerts an antiplatelet effect that together with its well-known potent anticoagulatory capacity might lead to reduced frequency of atherothrombotic events and improved outcome in patients.

Structural model of a CRISPR RNA-silencing complex reveals the RNA-target cleavage activity in Cmr4.

The Cmr complex is an RNA-guided endonuclease that cleaves foreign RNA targets as part of the CRISPR prokaryotic defense system. We investigated the molecular architecture of the P. furiosus Cmr complex using an integrative structural biology approach. We determined crystal structures of P. furiosus Cmr1, Cmr2, Cmr4, and Cmr6 and combined them with known structural information to interpret the cryo-EM map of the complex. To support structure determination, we obtained residue-specific interaction data using protein crosslinking and mass spectrometry. The resulting pseudoatomic model reveals how the superhelical backbone of the complex is defined by the polymerizing principles of Cmr4 and Cmr5 and how it is capped at the extremities by proteins of similar folds. The inner surface of the superhelix exposes conserved residues of Cmr4 that we show are required for target-cleavage activity. The structural and biochemical data thus identify Cmr4 as the conserved endoribonuclease of the Cmr complex.

Acute Phase Response as a Biological Mechanism-of-Action of (Nano)particle-Induced Cardiovascular Disease.

Inhaled nanoparticles constitute a potential health hazard due to their size-dependent lung deposition and large surface to mass ratio. Exposure to high levels contributes to the risk of developing respiratory and cardiovascular diseases, as well as of lung cancer. Particle-induced acute phase response may be an important mechanism of action of particle-induced cardiovascular disease. Here, the authors review new important scientific evidence showing causal relationships between inhalation of particle and nanomaterials, induction of acute phase response, and risk of cardiovascular disease. Particle-induced acute phase response provides a means for risk assessment of particle-induced cardiovascular disease and underscores cardiovascular disease as an occupational disease.

Variability of forced vital capacity in progressive interstitial lung disease: a prospective observational study.

BACKGROUND: Fibrotic interstitial lung disease (ILD) is often associated with poor outcomes, but has few predictors of progression. Daily home spirometry has been proposed to provide important information about the clinical course of idiopathic pulmonary disease (IPF). However, experience is limited, and home spirometry is not a routine component of patient care in ILD. Using home spirometry, we aimed to investigate the predictive potential of daily measurements of forced vital capacity (FVC) in fibrotic ILD. METHODS: In this prospective observational study, patients with fibrotic ILD and clinical progression were provided with home spirometers for daily measurements over 6 months. Hospital based spirometry was performed after three and 6 months. Disease progression, defined as death, lung transplantation, acute exacerbation or FVC decline > 10% relative was assessed in the cohort. RESULTS: From May 2017 until August 2018, we included 47 patients (IPF n = 20; non-IPF n = 27). Sufficient daily measurements were performed by 85.1% of the study cohort. Among these 40 patients (IPF n = 17; non-IPF n = 23), who had a mean ± SD age of 60.7 ± 11.3 years and FVC 64.7 ± 21.7% predicted (2.4 ± 0.8 L), 12 patients experienced disease progression (death: n = 2; lung transplantation: n = 3; acute exacerbation: n = 1; FVC decline > 10%: n = 6). Within the first 28 days, a group of patients had high daily variability in FVC, with 60.0% having a variation ≥5%. Patients with disease progression had significantly higher FVC variability than those in the stable group (median variability 8.6% vs. 4.8%; p = 0.002). Cox regression identified FVC variability as independently associated with disease progression when controlling for multiple confounding variables (hazard ratio: 1.203; 95% CI:1.050-1.378; p = 0.0076). CONCLUSIONS: Daily home spirometry is feasible in IPF and non-IPF ILD and facilitates the identification of FVC variability, which was associated with disease progression.

Compartment-resolved Proteomic Analysis of Mouse Aorta during Atherosclerotic Plaque Formation Reveals Osteoclast-specific Protein Expression.

Atherosclerosis leads to vascular lesions that involve major rearrangements of the vascular proteome, especially of the extracellular matrix (ECM). Using single aortas from ApoE knock out mice, we quantified formation of plaques by single-run, high-resolution mass spectrometry (MS)-based proteomics. To probe localization on a proteome-wide scale we employed quantitative detergent solubility profiling. This compartment- and time-resolved resource of atherogenesis comprised 5117 proteins, 182 of which changed their expression status in response to vessel maturation and atherosclerotic plaque development. In the insoluble ECM proteome, 65 proteins significantly changed, including relevant collagens, matrix metalloproteinases and macrophage derived proteins. Among novel factors in atherosclerosis, we identified matrilin-2, the collagen IV crosslinking enzyme peroxidasin as well as the poorly characterized MAM-domain containing 2 (Mamdc2) protein as being up-regulated in the ECM during atherogenesis. Intriguingly, three subunits of the osteoclast specific V-ATPase complex were strongly increased in mature plaques with an enrichment in macrophages thus implying an active de-mineralization function.