Epigenetic memory of coronavirus infection in innate immune cells and their progenitors.

Inflammation can trigger lasting phenotypes in immune and non-immune cells. Whether and how human infections and associated inflammation can form innate immune memory in hematopoietic stem and progenitor cells (HSPC) has remained unclear. We found that circulating HSPC, enriched from peripheral blood, captured the diversity of bone marrow HSPC, enabling investigation of their epigenomic reprogramming following coronavirus disease 2019 (COVID-19). Alterations in innate immune phenotypes and epigenetic programs of HSPC persisted for months to 1 year following severe COVID-19 and were associated with distinct transcription factor (TF) activities, altered regulation of inflammatory programs, and durable increases in myelopoiesis. HSPC epigenomic alterations were conveyed, through differentiation, to progeny innate immune cells. Early activity of IL-6 contributed to these persistent phenotypes in human COVID-19 and a mouse coronavirus infection model. Epigenetic reprogramming of HSPC may underlie altered immune function following infection and be broadly relevant, especially for millions of COVID-19 survivors.

Precursors for Nonlymphoid-Tissue Treg Cells Reside in Secondary Lymphoid Organs and Are Programmed by the Transcription Factor BATF.

Specialized regulatory T (Treg) cells accumulate and perform homeostatic and regenerative functions in nonlymphoid tissues. Whether common precursors for nonlymphoid-tissue Treg cells exist and how they differentiate remain elusive. Using transcription factor nuclear factor, interleukin 3 regulated (Nfil3) reporter mice and single-cell RNA-sequencing (scRNA-seq), we identified two precursor stages of interleukin 33 (IL-33) receptor ST2-expressing nonlymphoid tissue Treg cells, which resided in the spleen and lymph nodes. Global chromatin profiling of nonlymphoid tissue Treg cells and the two precursor stages revealed a stepwise acquisition of chromatin accessibility and reprogramming toward the nonlymphoid-tissue Treg cell phenotype. Mechanistically, we identified and validated the transcription factor Batf as the driver of the molecular tissue program in the precursors. Understanding this tissue development program will help to harness regenerative properties of tissue Treg cells for therapy.

The spatial landscape of lung pathology during COVID-19 progression.

Recent studies have provided insights into the pathology of and immune response to COVID-191-8. However, a thorough investigation of the interplay between infected cells and the immune system at sites of infection has been lacking. Here we use high-parameter imaging mass cytometry9 that targets the expression of 36 proteins to investigate the cellular composition and spatial architecture of acute lung injury in humans (including injuries derived from SARS-CoV-2 infection) at single-cell resolution. These spatially resolved single-cell data unravel the disordered structure of the infected and injured lung, alongside the distribution of extensive immune infiltration. Neutrophil and macrophage infiltration are hallmarks of bacterial pneumonia and COVID-19, respectively. We provide evidence that SARS-CoV-2 infects predominantly alveolar epithelial cells and induces a localized hyperinflammatory cell state that is associated with lung damage. We leverage the temporal range of fatal outcomes of COVID-19 in relation to the onset of symptoms, which reveals increased macrophage extravasation and increased numbers of mesenchymal cells and fibroblasts concomitant with increased proximity between these cell types as the disease progresses-possibly as a result of attempts to repair the damaged lung tissue. Our data enable us to develop a biologically interpretable landscape of lung pathology from a structural, immunological and clinical standpoint. We use this landscape to characterize the pathophysiology of the human lung from its macroscopic presentation to the single-cell level, which provides an important basis for understanding COVID-19 and lung pathology in general.

A molecular single-cell lung atlas of lethal COVID-19.

Respiratory failure is the leading cause of death in patients with severe SARS-CoV-2 infection1,2, but the host response at the lung tissue level is poorly understood. Here we performed single-nucleus RNA sequencing of about 116,000 nuclei from the lungs of nineteen individuals who died of COVID-19 and underwent rapid autopsy and seven control individuals. Integrated analyses identified substantial alterations in cellular composition, transcriptional cell states, and cell-to-cell interactions, thereby providing insight into the biology of lethal COVID-19. The lungs from individuals with COVID-19 were highly inflamed, with dense infiltration of aberrantly activated monocyte-derived macrophages and alveolar macrophages, but had impaired T cell responses. Monocyte/macrophage-derived interleukin-1ß and epithelial cell-derived interleukin-6 were unique features of SARS-CoV-2 infection compared to other viral and bacterial causes of pneumonia. Alveolar type 2 cells adopted an inflammation-associated transient progenitor cell state and failed to undergo full transition into alveolar type 1 cells, resulting in impaired lung regeneration. Furthermore, we identified expansion of recently described CTHRC1+ pathological fibroblasts3 contributing to rapidly ensuing pulmonary fibrosis in COVID-19. Inference of protein activity and ligand-receptor interactions identified putative drug targets to disrupt deleterious circuits. This atlas enables the dissection of lethal COVID-19, may inform our understanding of long-term complications of COVID-19 survivors, and provides an important resource for therapeutic development.

Structural cells are key regulators of organ-specific immune responses.

The mammalian immune system implements a remarkably effective set of mechanisms for fighting pathogens1. Its main components are haematopoietic immune cells, including myeloid cells that control innate immunity, and lymphoid cells that constitute adaptive immunity2. However, immune functions are not unique to haematopoietic cells, and many other cell types display basic mechanisms of pathogen defence3-5. To advance our understanding of immunology outside the haematopoietic system, here we systematically investigate the regulation of immune genes in the three major types of structural cells: epithelium, endothelium and fibroblasts. We characterize these cell types across twelve organs in mice, using cellular phenotyping, transcriptome sequencing, chromatin accessibility profiling and epigenome mapping. This comprehensive dataset revealed complex immune gene activity and regulation in structural cells. The observed patterns were highly organ-specific and seem to modulate the extensive interactions between structural cells and haematopoietic immune cells. Moreover, we identified an epigenetically encoded immune potential in structural cells under tissue homeostasis, which was triggered in response to systemic viral infection. This study highlights the prevalence and organ-specific complexity of immune gene activity in non-haematopoietic structural cells, and it provides a high-resolution, multi-omics atlas of the epigenetic and transcriptional networks that regulate structural cells in the mouse.

Unsupervised discovery of tissue architecture in multiplexed imaging.

Multiplexed imaging and spatial transcriptomics enable highly resolved spatial characterization of cellular phenotypes, but still largely depend on laborious manual annotation to understand higher-order patterns of tissue organization. As a result, higher-order patterns of tissue organization are poorly understood and not systematically connected to disease pathology or clinical outcomes. To address this gap, we developed an approach called UTAG to identify and quantify microanatomical tissue structures in multiplexed images without human intervention. Our method combines information on cellular phenotypes with the physical proximity of cells to accurately identify organ-specific microanatomical domains in healthy and diseased tissue. We apply our method to various types of images across healthy and disease states to show that it can consistently detect higher-level architectures in human tissues, quantify structural differences between healthy and diseased tissue, and reveal tissue organization patterns at the organ scale.

Ultra-high-throughput single-cell RNA sequencing and perturbation screening with combinatorial fluidic indexing.

Cell atlas projects and high-throughput perturbation screens require single-cell sequencing at a scale that is challenging with current technology. To enable cost-effective single-cell sequencing for millions of individual cells, we developed 'single-cell combinatorial fluidic indexing' (scifi). The scifi-RNA-seq assay combines one-step combinatorial preindexing of entire transcriptomes inside permeabilized cells with subsequent single-cell RNA-seq using microfluidics. Preindexing allows us to load several cells per droplet and computationally demultiplex their individual expression profiles. Thereby, scifi-RNA-seq massively increases the throughput of droplet-based single-cell RNA-seq, and provides a straightforward way of multiplexing thousands of samples in a single experiment. Compared with multiround combinatorial indexing, scifi-RNA-seq provides an easy and efficient workflow. Compared to cell hashing methods, which flag and discard droplets containing more than one cell, scifi-RNA-seq resolves and retains individual transcriptomes from overloaded droplets. We benchmarked scifi-RNA-seq on various human and mouse cell lines, validated it for primary human T cells and applied it in a highly multiplexed CRISPR screen with single-cell transcriptome readout of T cell receptor activation.

A Unique Cellular Organization of Human Distal Airways and Its Disarray in Chronic Obstructive Pulmonary Disease.

Rationale: Remodeling and loss of distal conducting airways, including preterminal and terminal bronchioles (pre-TBs/TBs), underlie progressive airflow limitation in chronic obstructive pulmonary disease (COPD). The cellular basis of these structural changes remains unknown. Objectives: To identify biological changes in pre-TBs/TBs in COPD at single-cell resolution and determine their cellular origin. Methods: We established a novel method of distal airway dissection and performed single-cell transcriptomic profiling of 111,412 cells isolated from different airway regions of 12 healthy lung donors and pre-TBs of 5 patients with COPD. Imaging CyTOF and immunofluorescence analysis of pre-TBs/TBs from 24 healthy lung donors and 11 subjects with COPD were performed to characterize cellular phenotypes at a tissue level. Region-specific differentiation of basal cells isolated from proximal and distal airways was studied using an air-liquid interface model. Measurements and Main Results: The atlas of cellular heterogeneity along the proximal-distal axis of the human lung was assembled and identified region-specific cellular states, including SCGB3A2+ SFTPB+ terminal airway-enriched secretory cells (TASCs) unique to distal airways. TASCs were lost in COPD pre-TBs/TBs, paralleled by loss of region-specific endothelial capillary cells, increased frequency of CD8+ T cells normally enriched in proximal airways, and augmented IFN-γ signaling. Basal cells residing in pre-TBs/TBs were identified as a cellular origin of TASCs. Regeneration of TASCs by these progenitors was suppressed by IFN-γ. Conclusions: Altered maintenance of the unique cellular organization of pre-TBs/TBs, including loss of the region-specific epithelial differentiation in these bronchioles, represents the cellular manifestation and likely the cellular basis of distal airway remodeling in COPD.

Hyperglycemia Induces Trained Immunity in Macrophages and Their Precursors and Promotes Atherosclerosis.

BACKGROUND: Cardiovascular risk in diabetes remains elevated despite glucose-lowering therapies. We hypothesized that hyperglycemia induces trained immunity in macrophages, promoting persistent proatherogenic characteristics. METHODS: Bone marrow-derived macrophages from control mice and mice with diabetes were grown in physiological glucose (5 mmol/L) and subjected to RNA sequencing (n=6), assay for transposase accessible chromatin sequencing (n=6), and chromatin immunoprecipitation sequencing (n=6) for determination of hyperglycemia-induced trained immunity. Bone marrow transplantation from mice with (n=9) or without (n=6) diabetes into (normoglycemic) Ldlr-/- mice was used to assess its functional significance in vivo. Evidence of hyperglycemia-induced trained immunity was sought in human peripheral blood mononuclear cells from patients with diabetes (n=8) compared with control subjects (n=16) and in human atherosclerotic plaque macrophages excised by laser capture microdissection. RESULTS: In macrophages, high extracellular glucose promoted proinflammatory gene expression and proatherogenic functional characteristics through glycolysis-dependent mechanisms. Bone marrow-derived macrophages from diabetic mice retained these characteristics, even when cultured in physiological glucose, indicating hyperglycemia-induced trained immunity. Bone marrow transplantation from diabetic mice into (normoglycemic) Ldlr-/- mice increased aortic root atherosclerosis, confirming a disease-relevant and persistent form of trained innate immunity. Integrated assay for transposase accessible chromatin, chromatin immunoprecipitation, and RNA sequencing analyses of hematopoietic stem cells and bone marrow-derived macrophages revealed a proinflammatory priming effect in diabetes. The pattern of open chromatin implicated transcription factor Runt-related transcription factor 1 (Runx1). Similarly, transcriptomes of atherosclerotic plaque macrophages and peripheral leukocytes in patients with type 2 diabetes were enriched for Runx1 targets, consistent with a potential role in human disease. Pharmacological inhibition of Runx1 in vitro inhibited the trained phenotype. CONCLUSIONS: Hyperglycemia-induced trained immunity may explain why targeting elevated glucose is ineffective in reducing macrovascular risk in diabetes and suggests new targets for disease prevention and therapy.

Combined chemosensitivity and chromatin profiling prioritizes drug combinations in CLL.

The Bruton tyrosine kinase (BTK) inhibitor ibrutinib has substantially improved therapeutic options for chronic lymphocytic leukemia (CLL). Although ibrutinib is not curative, it has a profound effect on CLL cells and may create new pharmacologically exploitable vulnerabilities. To identify such vulnerabilities, we developed a systematic approach that combines epigenome profiling (charting the gene-regulatory basis of cell state) with single-cell chemosensitivity profiling (quantifying cell-type-specific drug response) and bioinformatic data integration. By applying our method to a cohort of matched patient samples collected before and during ibrutinib therapy, we identified characteristic ibrutinib-induced changes that provide a starting point for the rational design of ibrutinib combination therapies. Specifically, we observed and validated preferential sensitivity to proteasome, PLK1, and mTOR inhibitors during ibrutinib treatment. More generally, our study establishes a broadly applicable method for investigating treatment-specific vulnerabilities by integrating the complementary perspectives of epigenetic cell states and phenotypic drug responses in primary patient samples.

Mitochondria Are a Subset of Extracellular Vesicles Released by Activated Monocytes and Induce Type I IFN and TNF Responses in Endothelial Cells.

RATIONALE: Extracellular vesicles, including microvesicles, are increasingly recognized as important mediators in cardiovascular disease. The cargo and surface proteins they carry are considered to define their biological activity, including their inflammatory properties. Monocyte to endothelial cell signaling is a prerequisite for the propagation of inflammatory responses. However, the contribution of microvesicles in this process is poorly understood. OBJECTIVE: To elucidate the mechanisms by which microvesicles derived from activated monocytic cells exert inflammatory effects on endothelial cells. METHODS AND RESULTS: LPS (lipopolysaccharide)-stimulated monocytic cells release free mitochondria and microvesicles with mitochondrial content as demonstrated by flow cytometry, quantitative polymerase chain reaction, Western Blot, and transmission electron microscopy. Using RNAseq analysis and quantitative reverse transcription-polymerase chain reaction, we demonstrated that both mitochondria directly isolated from and microvesicles released by LPS-activated monocytic cells, as well as circulating microvesicles isolated from volunteers receiving low-dose LPS-injections, induce type I IFN (interferon), and TNF (tumor necrosis factor) responses in endothelial cells. Depletion of free mitochondria significantly reduced the ability of these microvesicles to induce type I IFN and TNF-dependent genes. We identified mitochondria-associated TNFα and RNA from stressed mitochondria as major inducers of these responses. Finally, we demonstrated that the proinflammatory potential of microvesicles and directly isolated mitochondria were drastically reduced when they were derived from monocytic cells with nonrespiring mitochondria or monocytic cells cultured in the presence of pyruvate or the mitochondrial reactive oxygen species scavenger MitoTEMPO. CONCLUSIONS: Mitochondria and mitochondria embedded in microvesicles constitute a major subset of extracellular vesicles released by activated monocytes, and their proinflammatory activity on endothelial cells is determined by the activation status of their parental cells. Thus, mitochondria may represent critical intercellular mediators in cardiovascular disease and other inflammatory settings associated with type I IFN and TNF signaling.

Pooled CRISPR screening with single-cell transcriptome readout.

CRISPR-based genetic screens are accelerating biological discovery, but current methods have inherent limitations. Widely used pooled screens are restricted to simple readouts including cell proliferation and sortable marker proteins. Arrayed screens allow for comprehensive molecular readouts such as transcriptome profiling, but at much lower throughput. Here we combine pooled CRISPR screening with single-cell RNA sequencing into a broadly applicable workflow, directly linking guide RNA expression to transcriptome responses in thousands of individual cells. Our method for CRISPR droplet sequencing (CROP-seq) enables pooled CRISPR screens with single-cell transcriptome resolution, which will facilitate high-throughput functional dissection of complex regulatory mechanisms and heterogeneous cell populations.

ChIPmentation: fast, robust, low-input ChIP-seq for histones and transcription factors.

Chromatin immunoprecipitation followed by sequencing (ChIP-seq) is widely used to map histone marks and transcription factor binding throughout the genome. Here we present ChIPmentation, a method that combines chromatin immunoprecipitation with sequencing library preparation by Tn5 transposase ('tagmentation'). ChIPmentation introduces sequencing-compatible adaptors in a single-step reaction directly on bead-bound chromatin, which reduces time, cost and input requirements, thus providing a convenient and broadly useful alternative to existing ChIP-seq protocols.

Malondialdehyde epitopes are sterile mediators of hepatic inflammation in hypercholesterolemic mice.

Diet-related health issues such as nonalcoholic fatty liver disease and cardiovascular disorders are known to have a major inflammatory component. However, the exact pathways linking diet-induced changes (e.g., hyperlipidemia) and the ensuing inflammation have remained elusive so far. We identified biological processes related to innate immunity and oxidative stress as prime response pathways in livers of low-density lipoprotein receptor-deficient mice on a Western-type diet using RNA sequencing and in silico functional analyses of transcriptome data. The observed changes were independent of the presence of microbiota and thus indicative of a role for sterile triggers. We further show that malondialdehyde (MDA) epitopes, products of lipid peroxidation and markers for enhanced oxidative stress, are detectable in hepatic inflammation predominantly on dying cells and stimulate cytokine secretion as well as leukocyte recruitment in vitro and in vivo. MDA-induced cytokine secretion in vitro was dependent on the presence of the scavenger receptors CD36 and MSR1. Moreover, in vivo neutralization of endogenously generated MDA epitopes by intravenous injection of a specific MDA antibody results in decreased hepatic inflammation in low-density lipoprotein receptor-deficient mice on a Western-type diet. CONCLUSION: Accumulation of MDA epitopes plays a major role during diet-induced hepatic inflammation and can be ameliorated by administration of an anti-MDA antibody. (Hepatology 2017;65:1181-1195).

Evolutionary conservation of the eumetazoan gene regulatory landscape.

Despite considerable differences in morphology and complexity of body plans among animals, a great part of the gene set is shared among Bilateria and their basally branching sister group, the Cnidaria. This suggests that the common ancestor of eumetazoans already had a highly complex gene repertoire. At present it is therefore unclear how morphological diversification is encoded in the genome. Here we address the possibility that differences in gene regulation could contribute to the large morphological divergence between cnidarians and bilaterians. To this end, we generated the first genome-wide map of gene regulatory elements in a nonbilaterian animal, the sea anemone Nematostella vectensis. Using chromatin immunoprecipitation followed by deep sequencing of five chromatin modifications and a transcriptional cofactor, we identified over 5000 enhancers in the Nematostella genome and could validate 75% of the tested enhancers in vivo. We found that in Nematostella, but not in yeast, enhancers are characterized by the same combination of histone modifications as in bilaterians, and these enhancers preferentially target developmental regulatory genes. Surprisingly, the distribution and abundance of gene regulatory elements relative to these genes are shared between Nematostella and bilaterian model organisms. Our results suggest that complex gene regulation originated at least 600 million yr ago, predating the common ancestor of eumetazoans.

Pooled multicolour tagging for visualizing subcellular protein dynamics.

Imaging-based methods are widely used for studying the subcellular localization of proteins in living cells. While routine for individual proteins, global monitoring of protein dynamics following perturbation typically relies on arrayed panels of fluorescently tagged cell lines, limiting throughput and scalability. Here, we describe a strategy that combines high-throughput microscopy, computer vision and machine learning to detect perturbation-induced changes in multicolour tagged visual proteomics cell (vpCell) pools. We use genome-wide and cancer-focused intron-targeting sgRNA libraries to generate vpCell pools and a large, arrayed collection of clones each expressing two different endogenously tagged fluorescent proteins. Individual clones can be identified in vpCell pools by image analysis using the localization patterns and expression level of the tagged proteins as visual barcodes, enabling simultaneous live-cell monitoring of large sets of proteins. To demonstrate broad applicability and scale, we test the effects of antiproliferative compounds on a pool with cancer-related proteins, on which we identify widespread protein localization changes and new inhibitors of the nuclear import/export machinery. The time-resolved characterization of changes in subcellular localization and abundance of proteins upon perturbation in a pooled format highlights the power of the vpCell approach for drug discovery and mechanism-of-action studies.

Transcriptomic signatures of progressive and regressive liver fibrosis and portal hypertension.

Persistent liver injury triggers a fibrogenic program that causes pathologic remodeling of the hepatic microenvironment (i.e., liver fibrosis) and portal hypertension. The dynamics of gene regulation during liver disease progression and early regression remain understudied. Here, we generated hepatic transcriptome profiles in two well-established liver disease models at peak fibrosis and during spontaneous regression after the removal of the inducing agents. We linked the dynamics of key disease readouts, such as portal pressure, collagen area, and transaminase levels, to differentially expressed genes, enabling the identification of transcriptomic signatures of progressive vs. regressive liver fibrosis and portal hypertension. These candidate biomarkers (e.g., Tcf4, Mmp7, Trem2, Spp1, Scube1, Islr) were validated in RNA sequencing datasets of patients with cirrhosis and portal hypertension, and those cured from hepatitis C infection. Finally, deconvolution identified major cell types and suggested an association of macrophage and portal hepatocyte signatures with portal hypertension and fibrosis area.