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.

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.

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.

Integrative Transcriptomic and Single-Cell Protein Characterization of Colorectal Carcinoma Delineates Distinct Tumor Immune Microenvironments Associated with Overall Survival.

Colorectal carcinoma (CRC) is a heterogeneous group of tumors with varying therapeutic response and prognosis, and evidence suggests the tumor immune microenvironment (TIME) plays a pivotal role. Using advanced molecular and spatial biology technologies, we aimed to evaluate the TIME in patients with CRC to determine whether specific alterations in the immune composition correlated with prognosis. We identified primary and metastatic tumor samples from 31 consented patients, which were profiled with whole-exome sequencing and bulk RNA-seq. Immune cell deconvolution followed by gene set enrichment analysis and unsupervised clustering was performed. A subset of tumors underwent in situ analysis of the TIME spatial composition at single-cell resolution through Imaging Mass Mass Cytometry. Gene set enrichment analysis revealed two distinct groups of advanced CRC, one with an immune activated phenotype and the other with a suppressed immune microenvironment. The activated TIME phenotype contained increased Th1 cells, activated dendritic cells, tertiary lymphoid structures, and higher counts of CD8+ T cells whereas the inactive or suppressed TIME contained increased macrophages and a higher M2/M1 ratio. Our findings were further supported by RNA-seq data analysis from the TCGA CRC database, in which unsupervised clustering also identified two separate groups. The immunosuppressed CRC TIME had a lower overall survival probability (HR 1.66, p=0.007). This study supports the pertinent role of the CRC immune microenvironment in tumor progression and patient prognosis. We characterized the immune cell composition to better understand the complexity and vital role that immune activity states of the TIME play in determining patient outcome.

The evolution of metastatic upper tract urothelial carcinoma through genomic-transcriptomic and single-cell protein markers analysis.

The molecular characteristics of metastatic upper tract urothelial carcinoma (UTUC) are not well understood, and there is a lack of knowledge regarding the genomic and transcriptomic differences between primary and metastatic UTUC. To address these gaps, we integrate whole-exome sequencing, RNA sequencing, and Imaging Mass Cytometry using lanthanide metal-conjugated antibodies of 44 tumor samples from 28 patients with high-grade primary and metastatic UTUC. We perform a spatially-resolved single-cell analysis of cancer, immune, and stromal cells to understand the evolution of primary to metastatic UTUC. We discover that actionable genomic alterations are frequently discordant between primary and metastatic UTUC tumors in the same patient. In contrast, molecular subtype membership and immune depletion signature are stable across primary and matched metastatic UTUC. Molecular and immune subtypes are consistent between bulk RNA-sequencing and mass cytometry of protein markers from 340,798 single cells. Molecular subtypes at the single-cell level are highly conserved between primary and metastatic UTUC tumors within the same patient.

Multi-parametric atlas of the pre-metastatic liver for prediction of metastatic outcome in early-stage pancreatic cancer.

Metastasis occurs frequently after resection of pancreatic cancer (PaC). In this study, we hypothesized that multi-parametric analysis of pre-metastatic liver biopsies would classify patients according to their metastatic risk, timing and organ site. Liver biopsies obtained during pancreatectomy from 49 patients with localized PaC and 19 control patients with non-cancerous pancreatic lesions were analyzed, combining metabolomic, tissue and single-cell transcriptomics and multiplex imaging approaches. Patients were followed prospectively (median 3 years) and classified into four recurrence groups; early (<6 months after resection) or late (>6 months after resection) liver metastasis (LiM); extrahepatic metastasis (EHM); and disease-free survivors (no evidence of disease (NED)). Overall, PaC livers exhibited signs of augmented inflammation compared to controls. Enrichment of neutrophil extracellular traps (NETs), Ki-67 upregulation and decreased liver creatine significantly distinguished those with future metastasis from NED. Patients with future LiM were characterized by scant T cell lobular infiltration, less steatosis and higher levels of citrullinated H3 compared to patients who developed EHM, who had overexpression of interferon target genes (MX1 and NR1D1) and an increase of CD11B+ natural killer (NK) cells. Upregulation of sortilin-1 and prominent NETs, together with the lack of T cells and a reduction in CD11B+ NK cells, differentiated patients with early-onset LiM from those with late-onset LiM. Liver profiles of NED closely resembled those of controls. Using the above parameters, a machine-learning-based model was developed that successfully predicted the metastatic outcome at the time of surgery with 78% accuracy. Therefore, multi-parametric profiling of liver biopsies at the time of PaC diagnosis may determine metastatic risk and organotropism and guide clinical stratification for optimal treatment selection.

Identification of a stem cell mediating osteoblast versus adipocyte lineage selection.

Most skeletal fragility disorders are characterized by bone loss with a concurrent gain in marrow adipocytes 1-8. This suggests that a cell that forms adipocytes at the expense of osteoblasts is central to the pathogenesis of skeletal disorders. However, this cellular point of bifurcation between adipocyte and osteoblast differentiation pathways remains unknown. Here, we identify a new cell type defined by co-expression of skeletal stem cell and adipocyte precursor markers, 9-13 (CD24+CD29+ skeletal stem cells (SSCs)), that serves as a key cellular point of bifurcation between the osteoblast and adipocyte differentiation pathways, giving rise to closely related osteoblast and adipocyte lineage-restricted precursors. CD24+CD29+SSCs comprise a small fraction of SSCs, and only this fraction displays full stemness features, including the ability to undergo serial transplantation. In line with serving as the osteoblast/adipocyte bipotent cell, the "bone to fat" tissue remodeling occurring in models of postmenopausal osteoporosis or after high fat diet exposure occur in part by reprogramming these CD24+CD29+SSCs to change their output of lineage-restricted precursors. Lastly, as subcutaneous white adipose tissue displays a similar set of CD24+CD29+ stem cells and related lineage-restricted progenitors, these findings provide a new schema explaining the stem cell basis of bone versus adipose tissue production that unifies multiple mesenchymal tissues.

Mapping human natural killer cell development in pediatric tonsil by imaging mass cytometry and high-resolution microscopy.

Natural killer (NK) cells develop from CD34+ progenitors in a stage-specific manner defined by changes in cell surface receptor expression and function. Secondary lymphoid tissues, including tonsil, are sites of human NK cell development. Here we present new insights into human NK cell development in pediatric tonsil using cyclic immunofluorescence and imaging mass cytometry. We show that NK cell subset localization and interactions are dependent on NK cell developmental stage and tissue residency. NK cell progenitors are found in the interfollicular domain in proximity to cytokine-expressing stromal cells that promote proliferation and maturation. Mature NK cells are primarily found in the T-cell rich parafollicular domain engaging in cell-cell interactions that differ depending on their stage and tissue residency. The presence of local inflammation results in changes in NK cell interactions, abundance, and localization. This study provides the first comprehensive atlas of human NK cell development in secondary lymphoid tissue.

Insulin feedback is a targetable resistance mechanism of PI3K inhibition in glioblastoma.

BACKGROUND: Insulin feedback is a critical mechanism responsible for the poor clinical efficacy of phosphatidylinositol 3-kinase (PI3K) inhibition in cancer, and hyperglycemia is an independent factor associated with poor prognosis in glioblastoma (GBM). We investigated combination anti-hyperglycemic therapy in a mouse model of GBM and evaluated the association of glycemic control in clinical trial data from patients with GBM. METHODS: The effect of the anti-hyperglycemic regimens, metformin and the ketogenic diet, was evaluated in combination with PI3K inhibition in patient-derived GBM cells and in an orthotopic GBM mouse model. Insulin feedback and the immune microenvironment were retrospectively evaluated in blood and tumor tissue from a Phase 2 clinical trial of buparlisib in patients with recurrent GBM. RESULTS: We found that PI3K inhibition induces hyperglycemia and hyperinsulinemia in mice and that combining metformin with PI3K inhibition improves the treatment efficacy in an orthotopic GBM xenograft model. Through examination of clinical trial data, we found that hyperglycemia was an independent factor associated with poor progression-free survival in patients with GBM. We also found that PI3K inhibition increased insulin receptor activation and T-cell and microglia abundance in tumor tissue from these patients. CONCLUSION: Reducing insulin feedback improves the efficacy of PI3K inhibition in GBM in mice, and hyperglycemia worsens progression-free survival in patients with GBM treated with PI3K inhibition. These findings indicate that hyperglycemia is a critical resistance mechanism associated with PI3K inhibition in GBM and that anti-hyperglycemic therapy may enhance PI3K inhibitor efficacy in GBM patients.

Persistent alveolar type 2 dysfunction and lung structural derangement in post-acute COVID-19.

SARS-CoV-2 infection can manifest as a wide range of respiratory and systemic symptoms well after the acute phase of infection in over 50% of patients. Key questions remain on the long-term effects of infection on tissue pathology in recovered COVID-19 patients. To address these questions we performed multiplexed imaging of post-mortem lung tissue from 12 individuals who died post-acute COVID-19 (PC) and compare them to lung tissue from patients who died during the acute phase of COVID-19, or patients who died with idiopathic pulmonary fibrosis (IPF), and otherwise healthy lung tissue. We find evidence of viral presence in the lung up to 359 days after the acute phase of disease, including in patients with negative nasopharyngeal swab tests. The lung of PC patients are characterized by the accumulation of senescent alveolar type 2 cells, fibrosis with hypervascularization of peribronchial areas and alveolar septa, as the most pronounced pathophysiological features. At the cellular level, lung disease of PC patients, while distinct, shares pathological features with the chronic pulmonary disease of IPF. which may help rationalize interventions for PC patients. Altogether, this study provides an important foundation for the understanding of the long-term effects of SARS-CoV-2 pulmonary infection at the microanatomical, cellular, and molecular level.

Tumor-immune microenvironment revealed by Imaging Mass Cytometry in a metastatic sarcomatoid urothelial carcinoma with a prolonged response to pembrolizumab.

Sarcomatoid urothelial carcinoma (SUC) is a rare subtype of urothelial carcinoma (UC) that typically presents at an advanced stage compared to more common variants of UC. Locally advanced and metastatic UC have a poor long-term survival following progression on first-line platinum-based chemotherapy. Antibodies directed against the programmed cell death 1 protein (PD-1) or its ligand (PD-L1) are now approved to be used in these scenarios. The need for reliable biomarkers for treatment stratification is still under research. Here, we present a novel case report of the first Imaging Mass Cytometry (IMC) analysis done in SUC to investigate the immune cell repertoire and PD-L1 expression in a patient who presented with metastatic SUC and experienced a prolonged response to the anti-PD1 immune checkpoint inhibitor pembrolizumab after progression on first-line chemotherapy. This case report provides an important platform for translating these findings to a larger cohort of UC and UC variants.

The spatio-temporal landscape of lung pathology in SARS-CoV-2 infection.

Recent studies have provided insights into the pathology and immune response to coronavirus disease 2019 (COVID-19)1-8. However thorough interrogation of the interplay between infected cells and the immune system at sites of infection is lacking. We use high parameter imaging mass cytometry9 targeting the expression of 36 proteins, to investigate at single cell resolution, the cellular composition and spatial architecture of human acute lung injury including SARS-CoV-2. This spatially resolved, single-cell data unravels 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 hyper-inflammatory cell state associated with lung damage. By leveraging the temporal range of COVID-19 severe fatal disease in relation to the time of symptom onset, we observe increased macrophage extravasation, mesenchymal cells, and fibroblasts abundance concomitant with increased proximity between these cell types as the disease progresses, possibly as an attempt to repair the damaged lung tissue. This spatially resolved single-cell data allowed us to develop a biologically interpretable landscape of lung pathology from a structural, immunological and clinical standpoint. This spatial single-cell landscape enabled the pathophysiological characterization of the human lung from its macroscopic presentation to the single-cell, providing an important basis for the understanding of COVID-19, and lung pathology in general.

Single-Cell Analysis of the Muscle Stem Cell Hierarchy Identifies Heterotypic Communication Signals Involved in Skeletal Muscle Regeneration.

Muscle regeneration relies on the regulation of muscle stem cells (MuSCs) through paracrine signaling interactions. We analyzed muscle regeneration in mice using single-cell RNA sequencing (scRNA-seq) and generated over 34,000 single-cell transcriptomes spanning four time-points. We identified 15 distinct cell types including heterogenous populations of muscle stem and progenitor cells. We resolved a hierarchical map of these myogenic cells by trajectory inference and observed stage-specific regulatory programs within this continuum. Through ligand-receptor interaction analysis, we identified over 100 candidate regeneration-associated paracrine communication pairs between MuSCs and non-myogenic cells. We show that myogenic stem/progenitor cells exhibit heterogeneous expression of multiple Syndecan proteins in cycling myogenic cells, suggesting that Syndecans may coordinate myogenic fate regulation. We performed ligand stimulation in vitro and confirmed that three paracrine factors (FGF2, TGFß1, and RSPO3) regulate myogenic cell proliferation in a Syndecan-dependent manner. Our study provides a scRNA-seq reference resource to investigate cell communication interactions in muscle regeneration.

Aneuploidy and a deregulated DNA damage response suggest haploinsufficiency in breast tissues of BRCA2 mutation carriers.

Women harboring heterozygous germline mutations of BRCA2 have a 50 to 80% risk of developing breast cancer, yet the pathogenesis of these cancers is poorly understood. To reveal early steps in BRCA2-associated carcinogenesis, we analyzed sorted cell populations from freshly-isolated, non-cancerous breast tissues of BRCA2 mutation carriers and matched controls. Single-cell whole-genome sequencing demonstrates that >25% of BRCA2 carrier (BRCA2mut/+ ) luminal progenitor (LP) cells exhibit sub-chromosomal copy number variations, which are rarely observed in non-carriers. Correspondingly, primary BRCA2mut/+ breast epithelia exhibit DNA damage together with attenuated replication checkpoint and apoptotic responses, and an age-associated expansion of the LP compartment. We provide evidence that these phenotypes do not require loss of the wild-type BRCA2 allele. Collectively, our findings suggest that BRCA2 haploinsufficiency and associated DNA damage precede histologic abnormalities in vivo. Using these hallmarks of cancer predisposition will yield unanticipated opportunities for improved risk assessment and prevention strategies in high-risk patients.

Decoupling genetics, lineages, and microenvironment in IDH-mutant gliomas by single-cell RNA-seq.

Tumor subclasses differ according to the genotypes and phenotypes of malignant cells as well as the composition of the tumor microenvironment (TME). We dissected these influences in isocitrate dehydrogenase (IDH)-mutant gliomas by combining 14,226 single-cell RNA sequencing (RNA-seq) profiles from 16 patient samples with bulk RNA-seq profiles from 165 patient samples. Differences in bulk profiles between IDH-mutant astrocytoma and oligodendroglioma can be primarily explained by distinct TME and signature genetic events, whereas both tumor types share similar developmental hierarchies and lineages of glial differentiation. As tumor grade increases, we find enhanced proliferation of malignant cells, larger pools of undifferentiated glioma cells, and an increase in macrophage over microglia expression programs in TME. Our work provides a unifying model for IDH-mutant gliomas and a general framework for dissecting the differences among human tumor subclasses.