RESUMEN
Metastasis is the leading cause of cancer-related deaths. It is unclear how intratumor heterogeneity (ITH) contributes to metastasis and how metastatic cells adapt to distant tissue environments. The study of these adaptations is challenged by the limited access to patient material and a lack of experimental models that appropriately recapitulate ITH. To investigate metastatic cell adaptations and the contribution of ITH to metastasis, we analyzed single-cell transcriptomes of matched primary tumors and metastases from patient-derived xenograft models of breast cancer. We found profound transcriptional differences between the primary tumor and metastatic cells. Primary tumors upregulated several metabolic genes, whereas motility pathway genes were upregulated in micrometastases, and stress response signaling was upregulated during progression. Additionally, we identified primary tumor gene signatures that were associated with increased metastatic potential and correlated with patient outcomes. Immune-regulatory control pathways were enriched in poorly metastatic primary tumors, whereas genes involved in epithelial-mesenchymal transition were upregulated in highly metastatic tumors. We found that ITH was dominated by epithelial-mesenchymal plasticity (EMP), which presented as a dynamic continuum with intermediate EMP cell states characterized by specific genes such as CRYAB and S100A2. Elevated expression of an intermediate EMP signature correlated with worse patient outcomes. Our findings identified inhibition of the intermediate EMP cell state as a potential therapeutic target to block metastasis.
Asunto(s)
Neoplasias de la Mama , Transición Epitelial-Mesenquimal , Metástasis de la Neoplasia , Análisis de la Célula Individual , Humanos , Neoplasias de la Mama/patología , Neoplasias de la Mama/genética , Neoplasias de la Mama/metabolismo , Femenino , Animales , Ratones , Regulación Neoplásica de la Expresión Génica , Transcriptoma , Proteínas de Neoplasias/genética , Proteínas de Neoplasias/metabolismo , Línea Celular TumoralRESUMEN
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.
RESUMEN
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.
Asunto(s)
Linfocitos T CD8-positivos , Colitis , Humanos , Células Endoteliales , Inhibidores del Factor de Necrosis Tumoral , Colitis/inducido químicamente , Colitis/tratamiento farmacológico , Linfocitos T CD4-Positivos , Esteroides/farmacología , Esteroides/uso terapéutico , Células del EstromaRESUMEN
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.
RESUMEN
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.
Asunto(s)
Colitis Ulcerosa , Humanos , Colitis Ulcerosa/tratamiento farmacológico , Colitis Ulcerosa/genética , Integrinas/genética , Multiómica , Proteómica , Fármacos Gastrointestinales/uso terapéutico , Resultado del Tratamiento , Estudios RetrospectivosRESUMEN
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.
RESUMEN
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 ).
Asunto(s)
Ascomicetos , Comunicación Celular , Femenino , Masculino , Animales , Ratones , Envejecimiento/genética , Senescencia Celular , Regulación hacia AbajoRESUMEN
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.
Asunto(s)
COVID-19 , Humanos , COVID-19/genética , SARS-CoV-2/genética , Células Endoteliales , ARN Viral , PulmónRESUMEN
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.
RESUMEN
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.
Asunto(s)
Envejecimiento , Senescencia Celular , Drosophila melanogaster , Animales , Envejecimiento/genética , Perfilación de la Expresión Génica , Transcriptoma , Drosophila melanogaster/citología , Drosophila melanogaster/genética , Drosophila melanogaster/fisiología , Atlas como AsuntoRESUMEN
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.
Asunto(s)
Proteínas Hedgehog , Insuficiencia Renal Crónica , Animales , Humanos , Ratones , Fibrosis , Proteínas Hedgehog/metabolismo , Inflamación , FN-kappa B , Factores de Necrosis Tumoral , Proteína con Dedos de Zinc GLI1RESUMEN
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.
Asunto(s)
Hipocampo , Células-Madre Neurales , Ratones , Animales , Células-Madre Neurales/metabolismo , Neurogénesis/fisiología , Envejecimiento/fisiología , Respuesta de Proteína Desplegada , Proliferación CelularRESUMEN
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.
Asunto(s)
Análisis de la Célula Individual , Transcriptoma , Ratones , Animales , Hibridación Fluorescente in Situ/métodos , Análisis de la Célula Individual/métodos , Transcriptoma/genética , Perfilación de la Expresión Génica/métodos , RNA-SeqRESUMEN
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.
Asunto(s)
COVID-19 , Humanos , COVID-19/diagnóstico , SARS-CoV-2/genética , Prueba de COVID-19 , Pandemias , Sensibilidad y EspecificidadRESUMEN
Nonalcoholic fatty liver disease (NAFLD) can be ameliorated by calorie restriction, which leads to the suppressed somatotroph axis. Paradoxically, the suppressed somatotroph axis is associated with patients with NAFLD and is correlated with the severity of fibrosis. How the somatotroph axis becomes dysregulated and whether the repressed somatotroph axis impacts liver damage during the progression of NAFLD are unclear. Here, we identify a regulatory branch of the hepatic integrated stress response (ISR), which represses the somatotroph axis in hepatocytes through ATF3, resulting in enhanced cell survival and reduced cell proliferation. In mouse models of NAFLD, the ISR represses the somatotroph axis, leading to reduced apoptosis and inflammation but decreased hepatocyte proliferation and exacerbated fibrosis in the liver. NAD+ repletion reduces the ISR, rescues the dysregulated somatotroph axis, and alleviates NAFLD. These results establish that the hepatic ISR suppresses the somatotroph axis to control cell fate decisions and liver damage in NAFLD.
Asunto(s)
Enfermedad del Hígado Graso no Alcohólico , Somatotrofos , Ratones , Animales , Enfermedad del Hígado Graso no Alcohólico/patología , Hígado/patología , Hepatocitos/patología , Cirrosis Hepática/patologíaRESUMEN
We carried out integrated host and pathogen metagenomic RNA and DNA next generation sequencing (mNGS) of whole blood (n = 221) and plasma (n = 138) from critically ill patients following hospital admission. We assigned patients into sepsis groups on the basis of clinical and microbiological criteria. From whole-blood gene expression data, we distinguished patients with sepsis from patients with non-infectious systemic inflammatory conditions using a trained bagged support vector machine (bSVM) classifier (area under the receiver operating characteristic curve (AUC) = 0.81 in the training set; AUC = 0.82 in a held-out validation set). Plasma RNA also yielded a transcriptional signature of sepsis with several genes previously reported as sepsis biomarkers, and a bSVM sepsis diagnostic classifier (AUC = 0.97 training set; AUC = 0.77 validation set). Pathogen detection performance of plasma mNGS varied on the basis of pathogen and site of infection. To improve detection of virus, we developed a secondary transcriptomic classifier (AUC = 0.94 training set; AUC = 0.96 validation set). We combined host and microbial features to develop an integrated sepsis diagnostic model that identified 99% of microbiologically confirmed sepsis cases, and predicted sepsis in 74% of suspected and 89% of indeterminate sepsis cases. In summary, we suggest that integrating host transcriptional profiling and broad-range metagenomic pathogen detection from nucleic acid is a promising tool for sepsis diagnosis.
Asunto(s)
Enfermedad Crítica , Sepsis , Adulto , Humanos , Estudios Prospectivos , Sepsis/diagnóstico , Estudios de Cohortes , ARNRESUMEN
Single-cell RNA sequencing (scRNA-seq) provides unique insights into the pathology and cellular origin of disease. We introduce single-cell disease relevance score (scDRS), an approach that links scRNA-seq with polygenic disease risk at single-cell resolution, independent of annotated cell types. scDRS identifies cells exhibiting excess expression across disease-associated genes implicated by genome-wide association studies (GWASs). We applied scDRS to 74 diseases/traits and 1.3 million single-cell gene-expression profiles across 31 tissues/organs. Cell-type-level results broadly recapitulated known cell-type-disease associations. Individual-cell-level results identified subpopulations of disease-associated cells not captured by existing cell-type labels, including T cell subpopulations associated with inflammatory bowel disease, partially characterized by their effector-like states; neuron subpopulations associated with schizophrenia, partially characterized by their spatial locations; and hepatocyte subpopulations associated with triglyceride levels, partially characterized by their higher ploidy levels. Genes whose expression was correlated with the scDRS score across cells (reflecting coexpression with GWAS disease-associated genes) were strongly enriched for gold-standard drug target and Mendelian disease genes.
Asunto(s)
Estudio de Asociación del Genoma Completo , Análisis de la Célula Individual , Perfilación de la Expresión Génica/métodos , Herencia Multifactorial/genética , RNA-Seq , Análisis de la Célula Individual/métodos , TriglicéridosRESUMEN
Unlike other respiratory viruses, SARS-CoV-2 disproportionately causes severe disease in older adults whereas disease burden in children is lower. To investigate whether differences in the upper airway immune response may contribute to this disparity, we compare nasopharyngeal gene expression in 83 children (<19-years-old; 38 with SARS-CoV-2, 11 with other respiratory viruses, 34 with no virus) and 154 older adults (>40-years-old; 45 with SARS-CoV-2, 28 with other respiratory viruses, 81 with no virus). Expression of interferon-stimulated genes is robustly activated in both children and adults with SARS-CoV-2 infection compared to the respective non-viral groups, with only subtle distinctions. Children, however, demonstrate markedly greater upregulation of pathways related to B cell and T cell activation and proinflammatory cytokine signaling, including response to TNF and production of IFNγ, IL-2 and IL-4. Cell type deconvolution confirms greater recruitment of B cells, and to a lesser degree macrophages, to the upper airway of children. Only children exhibit a decrease in proportions of ciliated cells, among the primary targets of SARS-CoV-2, upon infection. These findings demonstrate that children elicit a more robust innate and especially adaptive immune response to SARS-CoV-2 in the upper airway that likely contributes to their protection from severe disease in the lower airway.
Asunto(s)
COVID-19 , SARS-CoV-2 , Inmunidad Adaptativa/genética , Adulto , Anciano , COVID-19/genética , Niño , Expresión Génica , Humanos , Nasofaringe , Adulto JovenRESUMEN
Molecular characterization of cell types using single-cell transcriptome sequencing is revolutionizing cell biology and enabling new insights into the physiology of human organs. We created a human reference atlas comprising nearly 500,000 cells from 24 different tissues and organs, many from the same donor. This atlas enabled molecular characterization of more than 400 cell types, their distribution across tissues, and tissue-specific variation in gene expression. Using multiple tissues from a single donor enabled identification of the clonal distribution of T cells between tissues, identification of the tissue-specific mutation rate in B cells, and analysis of the cell cycle state and proliferative potential of shared cell types across tissues. Cell type-specific RNA splicing was discovered and analyzed across tissues within an individual.