1
|
Liang Z, Anderson HD, Locher V, O'Leary C, Riesenfeld SJ, Jabri B, McDonald BD, Bendelac A. Eomes expression identifies the early bone marrow precursor to classical NK cells. Nat Immunol 2024; 25:1172-1182. [PMID: 38871999 DOI: 10.1038/s41590-024-01861-6] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/03/2023] [Accepted: 05/01/2024] [Indexed: 06/15/2024]
Abstract
Natural killer (NK) cells traffic through the blood and mount cytolytic and interferon-γ (IFNγ)-focused responses to intracellular pathogens and tumors. Type 1 innate lymphoid cells (ILC1s) also produce type 1 cytokines but reside in tissues and are not cytotoxic. Whether these differences reflect discrete lineages or distinct states of a common cell type is not understood. Using single-cell RNA sequencing and flow cytometry, we focused on populations of TCF7+ cells that contained precursors for NK cells and ILC1s and identified a subset of bone marrow lineage-negative NK receptor-negative cells that expressed the transcription factor Eomes, termed EomeshiNKneg cells. Transfer of EomeshiNKneg cells into Rag2-/-Il2rg-/- recipients generated functional NK cells capable of preventing metastatic disease. By contrast, transfer of PLZF+ ILC precursors generated a mixture of ILC1s, ILC2s and ILC3s that lacked cytotoxic potential. These findings identified EomeshiNKneg cells as the bone marrow precursor to classical NK cells and demonstrated that the NK and ILC1 lineages diverged early during development.
Collapse
Affiliation(s)
- Zhitao Liang
- Committee on Immunology, University of Chicago, Chicago, IL, USA
- Department of Pathology, University of Chicago, Chicago, IL, USA
| | - Hope D Anderson
- Biophysical Sciences Graduate Program, University of Chicago, Chicago, IL, USA
| | - Veronica Locher
- Committee on Immunology, University of Chicago, Chicago, IL, USA
| | - Crystal O'Leary
- Department of Pathology, University of Chicago, Chicago, IL, USA
| | - Samantha J Riesenfeld
- Committee on Immunology, University of Chicago, Chicago, IL, USA
- Pritzker School of Molecular Engineering, University of Chicago, Chicago, IL, USA
- Department of Medicine, University of Chicago, Chicago, IL, USA
- Institute for Biophysical Dynamics, University of Chicago, Chicago, IL, USA
| | - Bana Jabri
- Committee on Immunology, University of Chicago, Chicago, IL, USA
- Department of Medicine, University of Chicago, Chicago, IL, USA
| | | | - Albert Bendelac
- Committee on Immunology, University of Chicago, Chicago, IL, USA
- Department of Pathology, University of Chicago, Chicago, IL, USA
| |
Collapse
|
2
|
Lazarov O, Disouky A, Sanborn M, Mostafa M, Sabitha K, Schantz A, Kim N, Pawlowski S, Honer W, Bennett D, Zhou Y, Keene C, Maienschein-Cline M, Rehman J. A roadmap to human hippocampal neurogenesis in adulthood, aging and AD. RESEARCH SQUARE 2024:rs.3.rs-4469965. [PMID: 38854131 PMCID: PMC11160907 DOI: 10.21203/rs.3.rs-4469965/v1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/11/2024]
Abstract
In the rodent, hippocampal neurogenesis plays critical roles in learning and memory1,2, is tightly regulated by inhibitory neurons3-7 and contributes to memory dysfunction in Alzheimer's disease (AD) mouse models8-10. In contrast, the mechanisms regulating neurogenesis in the adult human hippocampus, the dynamic shifts in the transcriptomic and epigenomic profiles in aging and AD and putative niche interactions within the cellular environment, remain largely unknown. Using single nuclei multi-omics of postmortem human hippocampi we map the molecular mechanisms of hippocampal neurogenesis across aging, cognitive decline, and AD neuropathology. Transcriptomic and epigenetic profiling of neural stem cells (NSCs), neuroblasts and immature neurons suggests that the earliest shift in the characteristics of neurogenesis takes place in NSCs in aging. Cognitive impairment was associated with changes in neuroblast profile. In AD, there was a widespread cessation of the transcription machinery in immature neurons, with robust downregulation of genes regulating ribosomal and mitochondrial function. Further, there was substantial loss of parvalbumin+ inhibitory neurons in the hippocampus in aging. The number of the rest of inhibitory neurons were reduced as a function of age and diagnosis. Notably, a similar system-level effect was observed between immature and inhibitory neurons in the transition from aging to AD, manifested by common molecular pathways that were ultimately lost in AD. The numbers of neuroblasts, immature and GABAergic neurons inversely correlated with extent of neuropathology. Using CellChat and NeuronChat, we inferred the ligands and receptors by which neurogenic cells communicate with their cellular environment. Loss of synaptic adhesion molecules and neurotransmitters, either sent or received by neurogenic cells, was observed in AD. Together, this study delineates the molecular mechanisms and dynamics of human neurogenesis, functional association with inhibitory neurons and a mechanism of hippocampal hyperexcitability in AD.
Collapse
Affiliation(s)
| | | | | | | | - K Sabitha
- The University of Illinois at Chicago
| | | | | | | | | | | | - Yi Zhou
- Institute of Neuroscience, Chinese Academy of Sciences
| | | | | | | |
Collapse
|
3
|
Wang L, Jin G, Zhou Q, Liu Y, Zhao X, Li Z, Yin N, Peng M. Induction of immortal-like and functional CAR T cells by defined factors. J Exp Med 2024; 221:e20232368. [PMID: 38530240 DOI: 10.1084/jem.20232368] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/22/2023] [Revised: 01/10/2024] [Accepted: 03/04/2024] [Indexed: 03/27/2024] Open
Abstract
Long-term antitumor efficacy of chimeric antigen receptor (CAR) T cells depends on their functional persistence in vivo. T cells with stem-like properties show better persistence, but factors conferring bona fide stemness to T cells remain to be determined. Here, we demonstrate the induction of CAR T cells into an immortal-like and functional state, termed TIF. The induction of CARTIF cells depends on the repression of two factors, BCOR and ZC3H12A, and requires antigen or CAR tonic signaling. Reprogrammed CARTIF cells possess almost infinite stemness, similar to induced pluripotent stem cells while retaining the functionality of mature T cells, resulting in superior antitumor effects. Following the elimination of target cells, CARTIF cells enter a metabolically dormant state, persisting in vivo with a saturable niche and providing memory protection. TIF represents a novel state of T cells with unprecedented stemness, which confers long-term functional persistence of CAR T cells in vivo and holds broad potential in T cell therapies.
Collapse
Affiliation(s)
- Lixia Wang
- State Key Laboratory of Molecular Oncology, Beijing Key Laboratory for Immunological Research on Chronic Diseases, School of Medicine, Institute for Immunology, Tsinghua University, Beijing, China
- SXMU-Tsinghua Collaborative Innovation Center for Frontier Medicine , Taiyuan, China
- Tsinghua-Peking Center for Life Sciences , Beijing, China
| | - Gang Jin
- State Key Laboratory of Molecular Oncology, Beijing Key Laboratory for Immunological Research on Chronic Diseases, School of Medicine, Institute for Immunology, Tsinghua University, Beijing, China
- SXMU-Tsinghua Collaborative Innovation Center for Frontier Medicine , Taiyuan, China
- Tsinghua-Peking Center for Life Sciences , Beijing, China
| | - Qiuping Zhou
- State Key Laboratory of Molecular Oncology, Beijing Key Laboratory for Immunological Research on Chronic Diseases, School of Medicine, Institute for Immunology, Tsinghua University, Beijing, China
- SXMU-Tsinghua Collaborative Innovation Center for Frontier Medicine , Taiyuan, China
- Tsinghua-Peking Center for Life Sciences , Beijing, China
| | - Yanyan Liu
- State Key Laboratory of Molecular Oncology, Beijing Key Laboratory for Immunological Research on Chronic Diseases, School of Medicine, Institute for Immunology, Tsinghua University, Beijing, China
- SXMU-Tsinghua Collaborative Innovation Center for Frontier Medicine , Taiyuan, China
- Tsinghua-Peking Center for Life Sciences , Beijing, China
| | - Xiaocui Zhao
- State Key Laboratory of Molecular Oncology, Beijing Key Laboratory for Immunological Research on Chronic Diseases, School of Medicine, Institute for Immunology, Tsinghua University, Beijing, China
- SXMU-Tsinghua Collaborative Innovation Center for Frontier Medicine , Taiyuan, China
- Tsinghua-Peking Center for Life Sciences , Beijing, China
| | - Zhuoyang Li
- State Key Laboratory of Molecular Oncology, Beijing Key Laboratory for Immunological Research on Chronic Diseases, School of Medicine, Institute for Immunology, Tsinghua University, Beijing, China
- SXMU-Tsinghua Collaborative Innovation Center for Frontier Medicine , Taiyuan, China
- Tsinghua-Peking Center for Life Sciences , Beijing, China
| | - Na Yin
- State Key Laboratory of Molecular Oncology, Beijing Key Laboratory for Immunological Research on Chronic Diseases, School of Medicine, Institute for Immunology, Tsinghua University, Beijing, China
- SXMU-Tsinghua Collaborative Innovation Center for Frontier Medicine , Taiyuan, China
- Tsinghua-Peking Center for Life Sciences , Beijing, China
| | - Min Peng
- State Key Laboratory of Molecular Oncology, Beijing Key Laboratory for Immunological Research on Chronic Diseases, School of Medicine, Institute for Immunology, Tsinghua University, Beijing, China
- SXMU-Tsinghua Collaborative Innovation Center for Frontier Medicine , Taiyuan, China
- Tsinghua-Peking Center for Life Sciences , Beijing, China
| |
Collapse
|
4
|
Gao CF, Vaikuntanathan S, Riesenfeld SJ. Dissection and integration of bursty transcriptional dynamics for complex systems. Proc Natl Acad Sci U S A 2024; 121:e2306901121. [PMID: 38669186 PMCID: PMC11067469 DOI: 10.1073/pnas.2306901121] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/26/2023] [Accepted: 03/06/2024] [Indexed: 04/28/2024] Open
Abstract
RNA velocity estimation is a potentially powerful tool to reveal the directionality of transcriptional changes in single-cell RNA-sequencing data, but it lacks accuracy, absent advanced metabolic labeling techniques. We developed an approach, TopicVelo, that disentangles simultaneous, yet distinct, dynamics by using a probabilistic topic model, a highly interpretable form of latent space factorization, to infer cells and genes associated with individual processes, thereby capturing cellular pluripotency or multifaceted functionality. Focusing on process-associated cells and genes enables accurate estimation of process-specific velocities via a master equation for a transcriptional burst model accounting for intrinsic stochasticity. The method obtains a global transition matrix by leveraging cell topic weights to integrate process-specific signals. In challenging systems, this method accurately recovers complex transitions and terminal states, while our use of first-passage time analysis provides insights into transient transitions. These results expand the limits of RNA velocity, empowering future studies of cell fate and functional responses.
Collapse
Affiliation(s)
- Cheng Frank Gao
- Department of Chemistry, University of Chicago, Chicago, IL60637
| | - Suriyanarayanan Vaikuntanathan
- Department of Chemistry, University of Chicago, Chicago, IL60637
- Institute for Biophysical Dynamics, University of Chicago, Chicago, IL60637
| | - Samantha J. Riesenfeld
- Institute for Biophysical Dynamics, University of Chicago, Chicago, IL60637
- Pritzker School of Molecular Engineering, University of Chicago, Chicago, IL60637
- Department of Medicine, University of Chicago, Chicago, IL60637
- Committee on Immunology, Biological Sciences Division, University of Chicago, Chicago, IL60637
| |
Collapse
|
5
|
Chen F, Zou G, Wu Y, Ou-Yang L. Clustering single-cell multi-omics data via graph regularized multi-view ensemble learning. Bioinformatics 2024; 40:btae169. [PMID: 38547401 PMCID: PMC11015955 DOI: 10.1093/bioinformatics/btae169] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/11/2024] [Revised: 02/21/2024] [Accepted: 03/26/2024] [Indexed: 04/15/2024] Open
Abstract
MOTIVATION Single-cell clustering plays a crucial role in distinguishing between cell types, facilitating the analysis of cell heterogeneity mechanisms. While many existing clustering methods rely solely on gene expression data obtained from single-cell RNA sequencing techniques to identify cell clusters, the information contained in mono-omic data is often limited, leading to suboptimal clustering performance. The emergence of single-cell multi-omics sequencing technologies enables the integration of multiple omics data for identifying cell clusters, but how to integrate different omics data effectively remains challenging. In addition, designing a clustering method that performs well across various types of multi-omics data poses a persistent challenge due to the data's inherent characteristics. RESULTS In this paper, we propose a graph-regularized multi-view ensemble clustering (GRMEC-SC) model for single-cell clustering. Our proposed approach can adaptively integrate multiple omics data and leverage insights from multiple base clustering results. We extensively evaluate our method on five multi-omics datasets through a series of rigorous experiments. The results of these experiments demonstrate that our GRMEC-SC model achieves competitive performance across diverse multi-omics datasets with varying characteristics. AVAILABILITY AND IMPLEMENTATION Implementation of GRMEC-SC, along with examples, can be found on the GitHub repository: https://github.com/polarisChen/GRMEC-SC.
Collapse
Affiliation(s)
- Fuqun Chen
- College of Electronic and Information Engineering, Shenzhen University, Shenzhen 518060, Guangdong, China
- Guangdong Key Laboratory of Intelligent Information Processing, Shenzhen University, Shenzhen 518060, Guangdong, China
- Shenzhen Key Laboratory of Media Security and Guangdong Laboratory of Artificial Intelligence and Digital Economy (SZ), Shenzhen University, Shenzhen 518060, Guangdong, China
| | - Guanhua Zou
- College of Electronic and Information Engineering, Shenzhen University, Shenzhen 518060, Guangdong, China
- Guangdong Key Laboratory of Intelligent Information Processing, Shenzhen University, Shenzhen 518060, Guangdong, China
- Shenzhen Key Laboratory of Media Security and Guangdong Laboratory of Artificial Intelligence and Digital Economy (SZ), Shenzhen University, Shenzhen 518060, Guangdong, China
| | - Yongxian Wu
- College of Electronic and Information Engineering, Shenzhen University, Shenzhen 518060, Guangdong, China
- Guangdong Key Laboratory of Intelligent Information Processing, Shenzhen University, Shenzhen 518060, Guangdong, China
- Shenzhen Key Laboratory of Media Security and Guangdong Laboratory of Artificial Intelligence and Digital Economy (SZ), Shenzhen University, Shenzhen 518060, Guangdong, China
| | - Le Ou-Yang
- College of Electronic and Information Engineering, Shenzhen University, Shenzhen 518060, Guangdong, China
- Guangdong Key Laboratory of Intelligent Information Processing, Shenzhen University, Shenzhen 518060, Guangdong, China
- Shenzhen Key Laboratory of Media Security and Guangdong Laboratory of Artificial Intelligence and Digital Economy (SZ), Shenzhen University, Shenzhen 518060, Guangdong, China
| |
Collapse
|
6
|
Rezaie N, Rebboah E, Williams BA, Liang HY, Reese F, Balderrama-Gutierrez G, Dionne LA, Reinholdt L, Trout D, Wold BJ, Mortazavi A. Identification of robust cellular programs using reproducible LDA that impact sex-specific disease progression in different genotypes of a mouse model of AD. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.02.26.582178. [PMID: 38464087 PMCID: PMC10925135 DOI: 10.1101/2024.02.26.582178] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 03/12/2024]
Abstract
The gene expression profiles of distinct cell types reflect complex genomic interactions among multiple simultaneous biological processes within each cell that can be altered by disease progression as well as genetic background. The identification of these active cellular programs is an open challenge in the analysis of single-cell RNA-seq data. Latent Dirichlet Allocation (LDA) is a generative method used to identify recurring patterns in counts data, commonly referred to as topics that can be used to interpret the state of each cell. However, LDA's interpretability is hindered by several key factors including the hyperparameter selection of the number of topics as well as the variability in topic definitions due to random initialization. We developed Topyfic, a Reproducible LDA (rLDA) package, to accurately infer the identity and activity of cellular programs in single-cell data, providing insights into the relative contributions of each program in individual cells. We apply Topyfic to brain single-cell and single-nucleus datasets of two 5xFAD mouse models of Alzheimer's disease crossed with C57BL6/J or CAST/EiJ mice to identify distinct cell types and states in different cell types such as microglia. We find that 8-month 5xFAD/Cast F1 males show higher level of microglial activation than matching 5xFAD/BL6 F1 males, whereas female mice show similar levels of microglial activation. We show that regulatory genes such as TFs, microRNA host genes, and chromatin regulatory genes alone capture cell types and cell states. Our study highlights how topic modeling with a limited vocabulary of regulatory genes can identify gene expression programs in single-cell data in order to quantify similar and divergent cell states in distinct genotypes.
Collapse
Affiliation(s)
- Narges Rezaie
- Department of Developmental and Cell Biology, University of California, Irvine, CA, USA
- Center for Complex Biological Systems, University of California, Irvine, CA, USA
| | - Elisabeth Rebboah
- Department of Developmental and Cell Biology, University of California, Irvine, CA, USA
- Center for Complex Biological Systems, University of California, Irvine, CA, USA
| | - Brian A Williams
- Division of Biology, California Institute of Technology, Pasadena, CA, USA
| | - Heidi Yahan Liang
- Department of Developmental and Cell Biology, University of California, Irvine, CA, USA
- Center for Complex Biological Systems, University of California, Irvine, CA, USA
| | - Fairlie Reese
- Department of Developmental and Cell Biology, University of California, Irvine, CA, USA
- Center for Complex Biological Systems, University of California, Irvine, CA, USA
| | - Gabriela Balderrama-Gutierrez
- Department of Developmental and Cell Biology, University of California, Irvine, CA, USA
- Center for Complex Biological Systems, University of California, Irvine, CA, USA
| | | | | | - Diane Trout
- Division of Biology, California Institute of Technology, Pasadena, CA, USA
| | - Barbara J Wold
- Division of Biology, California Institute of Technology, Pasadena, CA, USA
| | - Ali Mortazavi
- Department of Developmental and Cell Biology, University of California, Irvine, CA, USA
- Center for Complex Biological Systems, University of California, Irvine, CA, USA
| |
Collapse
|
7
|
Aldana JA, Moa B, Mattsson J, Russell JH, Hawkins BJ. Histological, chemical and gene expression differences between western redcedar seedlings resistant and susceptible to cedar leaf blight. FRONTIERS IN PLANT SCIENCE 2024; 15:1309762. [PMID: 38379949 PMCID: PMC10878471 DOI: 10.3389/fpls.2024.1309762] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 10/08/2023] [Accepted: 01/08/2024] [Indexed: 02/22/2024]
Abstract
Introduction Western redcedar (Thuja plicata) is an important species in the Cupressaceae both at economic and cultural levels in the Pacific Northwest of North America. In adult trees, the species produces one of the most weathering-resistant heartwoods among conifers, making it one of the preferred species for outdoor applications. However, young T. plicata plants are susceptible to infection with cedar leaf blight (Didymascella thujina), an important foliar pathogen that can be devastating in nurseries and small-spaced plantations. Despite that, variability in the resistance against D. thujina in T. plicata has been documented, and such variability can be used to breed T. plicata for resistance against the pathogen. Objective This investigation aimed to discern the phenotypic and gene expression differences between resistant and susceptible T. plicata seedlings to shed light on the potential constitutive resistance mechanisms against cedar leaf blight in western redcedar. Methods The study consisted of two parts. First, the histological differences between four resistant and four susceptible families that were never infected with the pathogen were investigated. And second, the differences between one resistant and one susceptible family that were infected and not infected with the pathogen were analyzed at the chemical (C, N, mineral nutrients, lignin, fiber, starch, and terpenes) and gene expression (RNA-Seq) levels. Results The histological part showed that T. plicata seedlings resistant to D. thujina had constitutively thicker cuticles and lower stomatal densities than susceptible plants. The chemical analyses revealed that, regardless of their infection status, resistant plants had higher foliar concentrations of sabinene and α-thujene, and higher levels of expression of transcripts that code for leucine-rich repeat receptor-like protein kinases and for bark storage proteins. Conclusion The data collected in this study shows that constitutive differences at the phenotypic (histological and chemical) and gene expression level exist between T. plicata seedlings susceptible and resistant to D. thujina. Such differences have potential use for marker-assisted selection and breeding for resistance against cedar leaf blight in western redcedar in the future.
Collapse
Affiliation(s)
- Juan A. Aldana
- School of Arts, Science, and Education, Medicine Hat College, Medicine Hat, AB, Canada
| | - Belaid Moa
- Electrical and Computer Engineering Department, University of Victoria, Victoria, BC, Canada
| | - Jim Mattsson
- Department of Biological Sciences, Simon Fraser University, Burnaby, BC, Canada
| | - John H. Russell
- British Columbia Ministry of Forests, Mesachie Lake, BC, Canada
| | | |
Collapse
|
8
|
Takahama M, Patil A, Richey G, Cipurko D, Johnson K, Carbonetto P, Plaster M, Pandey S, Cheronis K, Ueda T, Gruenbaum A, Kawamoto T, Stephens M, Chevrier N. A pairwise cytokine code explains the organism-wide response to sepsis. Nat Immunol 2024; 25:226-239. [PMID: 38191855 PMCID: PMC10834370 DOI: 10.1038/s41590-023-01722-8] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/20/2023] [Accepted: 11/29/2023] [Indexed: 01/10/2024]
Abstract
Sepsis is a systemic response to infection with life-threatening consequences. Our understanding of the molecular and cellular impact of sepsis across organs remains rudimentary. Here, we characterize the pathogenesis of sepsis by measuring dynamic changes in gene expression across organs. To pinpoint molecules controlling organ states in sepsis, we compare the effects of sepsis on organ gene expression to those of 6 singles and 15 pairs of recombinant cytokines. Strikingly, we find that the pairwise effects of tumor necrosis factor plus interleukin (IL)-18, interferon-gamma or IL-1β suffice to mirror the impact of sepsis across tissues. Mechanistically, we map the cellular effects of sepsis and cytokines by computing changes in the abundance of 195 cell types across 9 organs, which we validate by whole-mouse spatial profiling. Our work decodes the cytokine cacophony in sepsis into a pairwise cytokine message capturing the gene, cell and tissue responses of the host to the disease.
Collapse
Affiliation(s)
- Michihiro Takahama
- Pritzker School of Molecular Engineering, University of Chicago, Chicago, IL, USA
- Laboratory of Bioresponse Regulation, Graduate School of Pharmaceutical Sciences, Osaka University, Osaka, Japan
| | | | - Gabriella Richey
- Pritzker School of Molecular Engineering, University of Chicago, Chicago, IL, USA
| | - Denis Cipurko
- Pritzker School of Molecular Engineering, University of Chicago, Chicago, IL, USA
| | - Katherine Johnson
- Pritzker School of Molecular Engineering, University of Chicago, Chicago, IL, USA
| | - Peter Carbonetto
- Department of Human Genetics, University of Chicago, Chicago, IL, USA
- Research Computing Center, University of Chicago, Chicago, IL, USA
| | - Madison Plaster
- Pritzker School of Molecular Engineering, University of Chicago, Chicago, IL, USA
| | - Surya Pandey
- Pritzker School of Molecular Engineering, University of Chicago, Chicago, IL, USA
| | - Katerina Cheronis
- Pritzker School of Molecular Engineering, University of Chicago, Chicago, IL, USA
| | - Tatsuki Ueda
- Pritzker School of Molecular Engineering, University of Chicago, Chicago, IL, USA
| | - Adam Gruenbaum
- Pritzker School of Molecular Engineering, University of Chicago, Chicago, IL, USA
| | | | - Matthew Stephens
- Department of Human Genetics, University of Chicago, Chicago, IL, USA
- Department of Statistics, University of Chicago, Chicago, IL, USA
| | - Nicolas Chevrier
- Pritzker School of Molecular Engineering, University of Chicago, Chicago, IL, USA.
| |
Collapse
|
9
|
Mosteiro L, Nguyen TTT, Hankeova S, Alvarez-Sierra D, Reichelt M, Vandriel SM, Lai Z, Choudhury FK, Sangaraju D, Kamath BM, Scherl A, Pujol-Borrell R, Piskol R, Siebel CW. Notch signaling in thyrocytes is essential for adult thyroid function and mammalian homeostasis. Nat Metab 2023; 5:2094-2110. [PMID: 38123718 DOI: 10.1038/s42255-023-00937-1] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 09/19/2022] [Accepted: 10/31/2023] [Indexed: 12/23/2023]
Abstract
The thyroid functions as an apex endocrine organ that controls growth, differentiation and metabolism1, and thyroid diseases comprise the most common endocrine disorders2. Nevertheless, high-resolution views of the cellular composition and signals that govern the thyroid have been lacking3,4. Here, we show that Notch signalling controls homeostasis and thermoregulation in adult mammals through a mitochondria-based mechanism in a subset of thyrocytes. We discover two thyrocyte subtypes in mouse and human thyroids, identified in single-cell analyses by different levels of metabolic activity and Notch signalling. Therapeutic antibody blockade of Notch in adult mice inhibits a thyrocyte-specific transcriptional program and induces thyrocyte defects due to decreased mitochondrial activity and ROS production. Thus, disrupting Notch signalling in adult mice causes hypothyroidism, characterized by reduced levels of circulating thyroid hormone and dysregulation of whole-body thermoregulation. Inducible genetic deletion of Notch1 and 2 in thyrocytes phenocopies this antibody-induced hypothyroidism, establishing a direct role for Notch in adult murine thyrocytes. We confirm that hypothyroidism is enriched in children with Alagille syndrome, a genetic disorder marked by Notch mutations, suggesting that these findings translate to humans.
Collapse
Grants
- NA Genentech (Genentech, Inc.)
- NA Genentech (Genentech, Inc.)
- NA Genentech (Genentech, Inc.)
- NA Genentech (Genentech, Inc.)
- NA Genentech (Genentech, Inc.)
- NA Genentech (Genentech, Inc.)
- NA Genentech (Genentech, Inc.)
- NA Genentech (Genentech, Inc.)
- NA Genentech (Genentech, Inc.)
- NA Genentech (Genentech, Inc.)
Collapse
Affiliation(s)
- Lluc Mosteiro
- Department of Discovery Oncology, Genentech, South San Francisco, CA, USA.
| | - Thi Thu Thao Nguyen
- Department of Oncology Bioinformatics, Genentech, South San Francisco, CA, USA
| | - Simona Hankeova
- Department of Discovery Oncology, Genentech, South San Francisco, CA, USA
| | - Daniel Alvarez-Sierra
- Translational Immunology Group, Vall d'Hebron Institut de Recerca (VHIR), Campus Vall Hebron, Barcelona, Spain
- Department of Cell Biology, Physiology and Immunology, Universitat Autònoma de Barcelona, Bellaterra, Spain
| | - Mike Reichelt
- Department of Research Pathology, Genentech, South San Francisco, CA, USA
| | - Shannon M Vandriel
- Division of Gastroenterology, Hepatology and Nutrition, Hospital for Sick Children, University of Toronto, Toronto, Ontario, Canada
| | - Zijuan Lai
- Department of Drug Metabolism and Pharmacokinetics, Genentech, South San Francisco, CA, USA
| | - Feroza K Choudhury
- Department of Drug Metabolism and Pharmacokinetics, Genentech, South San Francisco, CA, USA
| | - Dewakar Sangaraju
- Department of Drug Metabolism and Pharmacokinetics, Genentech, South San Francisco, CA, USA
| | - Binita M Kamath
- Division of Gastroenterology, Hepatology and Nutrition, Hospital for Sick Children, University of Toronto, Toronto, Ontario, Canada
| | - Alexis Scherl
- Department of Research Pathology, Genentech, South San Francisco, CA, USA
| | - Ricardo Pujol-Borrell
- Translational Immunology Group, Vall d'Hebron Institut de Recerca (VHIR), Campus Vall Hebron, Barcelona, Spain
- Department of Cell Biology, Physiology and Immunology, Universitat Autònoma de Barcelona, Bellaterra, Spain
- Vall Hebron Institute of Oncology (VHIO), Campus Vall Hebron, Barcelona, Spain
| | - Robert Piskol
- Department of Oncology Bioinformatics, Genentech, South San Francisco, CA, USA
| | - Christian W Siebel
- Department of Discovery Oncology, Genentech, South San Francisco, CA, USA.
| |
Collapse
|
10
|
Uttley K, Papanastasiou AS, Lahne M, Brisbane JM, MacDonald RB, Bickmore WA, Bhatia S. Unique activities of two overlapping PAX6 retinal enhancers. Life Sci Alliance 2023; 6:e202302126. [PMID: 37643867 PMCID: PMC10465922 DOI: 10.26508/lsa.202302126] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/01/2023] [Revised: 08/16/2023] [Accepted: 08/17/2023] [Indexed: 08/31/2023] Open
Abstract
Enhancers play a critical role in development by precisely modulating spatial, temporal, and cell type-specific gene expression. Sequence variants in enhancers have been implicated in diseases; however, establishing the functional consequences of these variants is challenging because of a lack of understanding of precise cell types and developmental stages where the enhancers are normally active. PAX6 is the master regulator of eye development, with a regulatory landscape containing multiple enhancers driving the expression in the eye. Whether these enhancers perform additive, redundant or distinct functions is unknown. Here, we describe the precise cell types and regulatory activity of two PAX6 retinal enhancers, HS5 and NRE. Using a unique combination of live imaging and single-cell RNA sequencing in dual enhancer-reporter zebrafish embryos, we uncover differences in the spatiotemporal activity of these enhancers. Our results show that although overlapping, these enhancers have distinct activities in different cell types and therefore likely nonredundant functions. This work demonstrates that unique cell type-specific activities can be uncovered for apparently similar enhancers when investigated at high resolution in vivo.
Collapse
Affiliation(s)
- Kirsty Uttley
- https://ror.org/011jsc803 MRC Human Genetics Unithttps://ror.org/01nrxwf90 , Institute of Genetics and Cancer, The University of Edinburgh, Edinburgh, UK
| | - Andrew S Papanastasiou
- https://ror.org/011jsc803 MRC Human Genetics Unithttps://ror.org/01nrxwf90 , Institute of Genetics and Cancer, The University of Edinburgh, Edinburgh, UK
| | - Manuela Lahne
- https://ror.org/02jx3x895 UCL Institute of Ophthalmology, University College London, Greater London, UK
| | - Jennifer M Brisbane
- https://ror.org/011jsc803 MRC Human Genetics Unithttps://ror.org/01nrxwf90 , Institute of Genetics and Cancer, The University of Edinburgh, Edinburgh, UK
| | - Ryan B MacDonald
- https://ror.org/02jx3x895 UCL Institute of Ophthalmology, University College London, Greater London, UK
| | - Wendy A Bickmore
- https://ror.org/011jsc803 MRC Human Genetics Unithttps://ror.org/01nrxwf90 , Institute of Genetics and Cancer, The University of Edinburgh, Edinburgh, UK
| | - Shipra Bhatia
- https://ror.org/011jsc803 MRC Human Genetics Unithttps://ror.org/01nrxwf90 , Institute of Genetics and Cancer, The University of Edinburgh, Edinburgh, UK
| |
Collapse
|
11
|
Carbonetto P, Luo K, Sarkar A, Hung A, Tayeb K, Pott S, Stephens M. GoM DE: interpreting structure in sequence count data with differential expression analysis allowing for grades of membership. Genome Biol 2023; 24:236. [PMID: 37858253 PMCID: PMC10588049 DOI: 10.1186/s13059-023-03067-9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/03/2023] [Accepted: 09/20/2023] [Indexed: 10/21/2023] Open
Abstract
Parts-based representations, such as non-negative matrix factorization and topic modeling, have been used to identify structure from single-cell sequencing data sets, in particular structure that is not as well captured by clustering or other dimensionality reduction methods. However, interpreting the individual parts remains a challenge. To address this challenge, we extend methods for differential expression analysis by allowing cells to have partial membership to multiple groups. We call this grade of membership differential expression (GoM DE). We illustrate the benefits of GoM DE for annotating topics identified in several single-cell RNA-seq and ATAC-seq data sets.
Collapse
Affiliation(s)
- Peter Carbonetto
- Department of Human Genetics, University of Chicago, Chicago, IL, USA
- Research Computing Center, University of Chicago, Chicago, IL, USA
| | - Kaixuan Luo
- Department of Human Genetics, University of Chicago, Chicago, IL, USA
| | - Abhishek Sarkar
- Department of Human Genetics, University of Chicago, Chicago, IL, USA
- Vesalius Therapeutics, Cambridge, MA, USA
| | - Anthony Hung
- Department of Human Genetics, University of Chicago, Chicago, IL, USA
- Section of Genetic Medicine, University of Chicago, Chicago, IL, USA
| | - Karl Tayeb
- Department of Human Genetics, University of Chicago, Chicago, IL, USA
- Committee on Genetics, Genomics and Systems Biology, University of Chicago, Chicago, IL, USA
| | - Sebastian Pott
- Department of Human Genetics, University of Chicago, Chicago, IL, USA
- Section of Genetic Medicine, University of Chicago, Chicago, IL, USA
| | - Matthew Stephens
- Department of Human Genetics, University of Chicago, Chicago, IL, USA.
- Department of Statistics, University of Chicago, Chicago, IL, USA.
| |
Collapse
|
12
|
Fukuyama J, Sankaran K, Symul L. Multiscale analysis of count data through topic alignment. Biostatistics 2023; 24:1045-1065. [PMID: 35657012 DOI: 10.1093/biostatistics/kxac018] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/12/2021] [Revised: 03/10/2022] [Accepted: 03/21/2022] [Indexed: 10/19/2023] Open
Abstract
Topic modeling is a popular method used to describe biological count data. With topic models, the user must specify the number of topics $K$. Since there is no definitive way to choose $K$ and since a true value might not exist, we develop a method, which we call topic alignment, to study the relationships across models with different $K$. In addition, we present three diagnostics based on the alignment. These techniques can show how many topics are consistently present across different models, if a topic is only transiently present, or if a topic splits into more topics when $K$ increases. This strategy gives more insight into the process of generating the data than choosing a single value of $K$ would. We design a visual representation of these cross-model relationships, show the effectiveness of these tools for interpreting the topics on simulated and real data, and release an accompanying R package, alto.
Collapse
Affiliation(s)
- Julia Fukuyama
- Department of Statistics, Indiana University Bloomington, 919 E 10th Street, Bloomington, IN 47408, USA
| | - Kris Sankaran
- Department of Statistics, University of Wisconsin - Madison, 1300 University Ave, Madison, WI 53706, USA
| | - Laura Symul
- Department of Statistics, Stanford University, 390 Jane Stanford Way, Stanford, CA 94305, USA
| |
Collapse
|
13
|
Taraborrelli L, Şenbabaoğlu Y, Wang L, Lim J, Blake K, Kljavin N, Gierke S, Scherl A, Ziai J, McNamara E, Owyong M, Rao S, Calviello AK, Oreper D, Jhunjhunwala S, Argiles G, Bendell J, Kim TW, Ciardiello F, Wongchenko MJ, de Sauvage FJ, de Sousa E Melo F, Yan Y, West NR, Murthy A. Tumor-intrinsic expression of the autophagy gene Atg16l1 suppresses anti-tumor immunity in colorectal cancer. Nat Commun 2023; 14:5945. [PMID: 37741832 PMCID: PMC10517947 DOI: 10.1038/s41467-023-41618-7] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/31/2023] [Accepted: 09/12/2023] [Indexed: 09/25/2023] Open
Abstract
Microsatellite-stable colorectal cancer (MSS-CRC) is highly refractory to immunotherapy. Understanding tumor-intrinsic determinants of immunotherapy resistance is critical to improve MSS-CRC patient outcomes. Here, we demonstrate that high tumor expression of the core autophagy gene ATG16L1 is associated with poor clinical response to anti-PD-L1 therapy in KRAS-mutant tumors from IMblaze370 (NCT02788279), a large phase III clinical trial of atezolizumab (anti-PD-L1) in advanced metastatic MSS-CRC. Deletion of Atg16l1 in engineered murine colon cancer organoids inhibits tumor growth in primary (colon) and metastatic (liver and lung) niches in syngeneic female hosts, primarily due to increased sensitivity to IFN-γ-mediated immune pressure. ATG16L1 deficiency enhances programmed cell death of colon cancer organoids induced by IFN-γ and TNF, thus increasing their sensitivity to host immunity. In parallel, ATG16L1 deficiency reduces tumor stem-like populations in vivo independently of adaptive immune pressure. This work reveals autophagy as a clinically relevant mechanism of immune evasion and tumor fitness in MSS-CRC and provides a rationale for autophagy inhibition to boost immunotherapy responses in the clinic.
Collapse
Affiliation(s)
- Lucia Taraborrelli
- Department of Cancer Immunology, Genentech Inc., South San Francisco, USA
| | - Yasin Şenbabaoğlu
- Department of Oncology Bioinformatics, Genentech Inc., South San Francisco, USA
| | - Lifen Wang
- Department of Cancer Immunology, Genentech Inc., South San Francisco, USA
| | - Junghyun Lim
- Department of Cancer Immunology, Genentech Inc., South San Francisco, USA
| | - Kerrigan Blake
- Department of Cancer Immunology, Genentech Inc., South San Francisco, USA
| | - Noelyn Kljavin
- Department of Molecular Oncology, Genentech Inc., South San Francisco, USA
| | - Sarah Gierke
- Center for Advanced Light Microscopy, Genentech Inc., South San Francisco, USA
- Department of Pathology, Genentech Inc., South San Francisco, USA
| | - Alexis Scherl
- Department of Pathology, Genentech Inc., South San Francisco, USA
| | - James Ziai
- Department of Pathology, Genentech Inc., South San Francisco, USA
| | - Erin McNamara
- Department of In Vivo Pharmacology, Genentech Inc., South San Francisco, USA
| | - Mark Owyong
- Department of In Vivo Pharmacology, Genentech Inc., South San Francisco, USA
| | - Shilpa Rao
- Department of Oncology Bioinformatics, Genentech Inc., South San Francisco, USA
| | | | - Daniel Oreper
- Department of Oncology Bioinformatics, Genentech Inc., South San Francisco, USA
| | - Suchit Jhunjhunwala
- Department of Oncology Bioinformatics, Genentech Inc., South San Francisco, USA
| | - Guillem Argiles
- Vall d'Hebrón Institute of Oncology, Vall d'Hebrón University Hospital, Universitat Autònoma de Barcelona, Barcelona, Spain
| | - Johanna Bendell
- Sarah Cannon Research Institute/Tennessee Oncology, Nashville, TN, USA
| | - Tae Won Kim
- Department of Oncology, Medical Center, University of Ulsan, Seoul, Korea
| | - Fortunato Ciardiello
- Department of Precision Medicine, Università degli Studi della Campania Luigi Vanvitelli, Naples, Italy
| | | | | | | | - Yibing Yan
- Oncology Biomarker Development, Genentech, Inc., South San Francisco, CA, USA
| | - Nathaniel R West
- Department of Cancer Immunology, Genentech Inc., South San Francisco, USA.
| | - Aditya Murthy
- Department of Cancer Immunology, Genentech Inc., South San Francisco, USA.
- Gilead Sciences, Foster City, USA.
| |
Collapse
|
14
|
Karger A, Mansouri S, Leisegang MS, Weigert A, Günther S, Kuenne C, Wittig I, Zukunft S, Klatt S, Aliraj B, Klotz LV, Winter H, Mahavadi P, Fleming I, Ruppert C, Witte B, Alkoudmani I, Gattenlöhner S, Grimminger F, Seeger W, Pullamsetti SS, Savai R. ADPGK-AS1 long noncoding RNA switches macrophage metabolic and phenotypic state to promote lung cancer growth. EMBO J 2023; 42:e111620. [PMID: 37545364 PMCID: PMC10505917 DOI: 10.15252/embj.2022111620] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/09/2022] [Revised: 07/06/2023] [Accepted: 07/08/2023] [Indexed: 08/08/2023] Open
Abstract
Long noncoding RNAs (lncRNAs) influence the transcription of gene networks in many cell types, but their role in tumor-associated macrophages (TAMs) is still largely unknown. We found that the lncRNA ADPGK-AS1 was substantially upregulated in artificially induced M2-like human macrophages, macrophages exposed to lung cancer cells in vitro, and TAMs from human lung cancer tissue. ADPGK-AS1 is partly located within mitochondria and binds to the mitochondrial ribosomal protein MRPL35. Overexpression of ADPGK-AS1 in macrophages upregulates the tricarboxylic acid cycle and promotes mitochondrial fission, suggesting a phenotypic switch toward an M2-like, tumor-promoting cytokine release profile. Macrophage-specific knockdown of ADPGK-AS1 induces a metabolic and phenotypic switch (as judged by cytokine profile and production of reactive oxygen species) to a pro-inflammatory tumor-suppressive M1-like state, inhibiting lung tumor growth in vitro in tumor cell-macrophage cocultures, ex vivo in human tumor precision-cut lung slices, and in vivo in mice. Silencing ADPGK-AS1 in TAMs may thus offer a novel therapeutic strategy for lung cancer.
Collapse
Affiliation(s)
- Annika Karger
- Max Planck Institute for Heart and Lung ResearchMember of the German Center for Lung Research (DZL), Member of the Cardio‐Pulmonary Institute (CPI)Bad NauheimGermany
- Institute for Lung Health (ILH)Justus Liebig UniversityGiessenGermany
| | - Siavash Mansouri
- Max Planck Institute for Heart and Lung ResearchMember of the German Center for Lung Research (DZL), Member of the Cardio‐Pulmonary Institute (CPI)Bad NauheimGermany
- Institute for Lung Health (ILH)Justus Liebig UniversityGiessenGermany
| | - Matthias S Leisegang
- Institute for Cardiovascular Physiology, Medical FacultyGoethe University FrankfurtFrankfurtGermany
| | - Andreas Weigert
- Institute of Biochemistry I, Faculty of MedicineGoethe University FrankfurtFrankfurtGermany
- Frankfurt Cancer Institute (FCI)Goethe University FrankfurtFrankfurtGermany
| | - Stefan Günther
- Max Planck Institute for Heart and Lung ResearchMember of the German Center for Lung Research (DZL), Member of the Cardio‐Pulmonary Institute (CPI)Bad NauheimGermany
| | - Carsten Kuenne
- Max Planck Institute for Heart and Lung ResearchMember of the German Center for Lung Research (DZL), Member of the Cardio‐Pulmonary Institute (CPI)Bad NauheimGermany
| | - Ilka Wittig
- Functional Proteomics, Medical SchoolGoethe University FrankfurtFrankfurtGermany
| | - Sven Zukunft
- Institute for Vascular Signalling, Centre for Molecular MedicineGoethe UniversityFrankfurtGermany
| | - Stephan Klatt
- Institute for Vascular Signalling, Centre for Molecular MedicineGoethe UniversityFrankfurtGermany
| | - Blerina Aliraj
- Institute of Biochemistry I, Faculty of MedicineGoethe University FrankfurtFrankfurtGermany
| | - Laura V Klotz
- Translational Lung Research Center (TLRC), Member of the DZLHeidelbergGermany
- Department of Thoracic SurgeryThoraxklinik at the University Hospital HeidelbergHeidelbergGermany
| | - Hauke Winter
- Translational Lung Research Center (TLRC), Member of the DZLHeidelbergGermany
- Department of Thoracic SurgeryThoraxklinik at the University Hospital HeidelbergHeidelbergGermany
| | - Poornima Mahavadi
- Department of Internal MedicineMember of the DZL, Member of CPI, Justus Liebig UniversityGiessenGermany
| | - Ingrid Fleming
- Institute for Vascular Signalling, Centre for Molecular MedicineGoethe UniversityFrankfurtGermany
| | - Clemens Ruppert
- Department of Internal MedicineMember of the DZL, Member of CPI, Justus Liebig UniversityGiessenGermany
| | - Biruta Witte
- Department of General and Thoracic SurgeryUniversity Hospital GiessenGiessenGermany
| | - Ibrahim Alkoudmani
- Department of General and Thoracic SurgeryUniversity Hospital GiessenGiessenGermany
| | | | - Friedrich Grimminger
- Institute for Lung Health (ILH)Justus Liebig UniversityGiessenGermany
- Department of Internal MedicineMember of the DZL, Member of CPI, Justus Liebig UniversityGiessenGermany
| | - Werner Seeger
- Max Planck Institute for Heart and Lung ResearchMember of the German Center for Lung Research (DZL), Member of the Cardio‐Pulmonary Institute (CPI)Bad NauheimGermany
- Institute for Lung Health (ILH)Justus Liebig UniversityGiessenGermany
- Department of Internal MedicineMember of the DZL, Member of CPI, Justus Liebig UniversityGiessenGermany
| | - Soni Savai Pullamsetti
- Max Planck Institute for Heart and Lung ResearchMember of the German Center for Lung Research (DZL), Member of the Cardio‐Pulmonary Institute (CPI)Bad NauheimGermany
- Institute for Lung Health (ILH)Justus Liebig UniversityGiessenGermany
- Department of Internal MedicineMember of the DZL, Member of CPI, Justus Liebig UniversityGiessenGermany
| | - Rajkumar Savai
- Max Planck Institute for Heart and Lung ResearchMember of the German Center for Lung Research (DZL), Member of the Cardio‐Pulmonary Institute (CPI)Bad NauheimGermany
- Institute for Lung Health (ILH)Justus Liebig UniversityGiessenGermany
- Frankfurt Cancer Institute (FCI)Goethe University FrankfurtFrankfurtGermany
- Department of Internal MedicineMember of the DZL, Member of CPI, Justus Liebig UniversityGiessenGermany
| |
Collapse
|
15
|
Carbonetto P, Luo K, Sarkar A, Hung A, Tayeb K, Pott S, Stephens M. GoM DE: interpreting structure in sequence count data with differential expression analysis allowing for grades of membership. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.03.03.531029. [PMID: 36945441 PMCID: PMC10028846 DOI: 10.1101/2023.03.03.531029] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 03/11/2023]
Abstract
Parts-based representations, such as non-negative matrix factorization and topic modeling, have been used to identify structure from single-cell sequencing data sets, in particular structure that is not as well captured by clustering or other dimensionality reduction methods. However, interpreting the individual parts remains a challenge. To address this challenge, we extend methods for differential expression analysis by allowing cells to have partial membership to multiple groups. We call this grade of membership differential expression (GoM DE). We illustrate the benefits of GoM DE for annotating topics identified in several single-cell RNA-seq and ATAC-seq data sets.
Collapse
Affiliation(s)
- Peter Carbonetto
- Department of Human Genetics, University of Chicago, Chicago, IL, USA
- Research Computing Center, University of Chicago, Chicago, IL, USA
| | - Kaixuan Luo
- Department of Human Genetics, University of Chicago, Chicago, IL, USA
| | - Abhishek Sarkar
- Department of Human Genetics, University of Chicago, Chicago, IL, USA
- Vesalius Therapeutics, Cambridge, MA, USA
| | - Anthony Hung
- Department of Human Genetics, University of Chicago, Chicago, IL, USA
- Section of Genetic Medicine, University of Chicago, Chicago, IL, USA
| | - Karl Tayeb
- Department of Human Genetics, University of Chicago, Chicago, IL, USA
- Committee on Genetics, Genomics and Systems Biology, University of Chicago, Chicago, IL, USA
| | - Sebastian Pott
- Department of Human Genetics, University of Chicago, Chicago, IL, USA
- Section of Genetic Medicine, University of Chicago, Chicago, IL, USA
| | - Matthew Stephens
- Department of Human Genetics, University of Chicago, Chicago, IL, USA
- Department of Statistics, University of Chicago, Chicago, IL, USA
| |
Collapse
|
16
|
Wang P, Kljavin N, Nguyen TTT, Storm EE, Marsh B, Jiang J, Lin W, Menon H, Piskol R, de Sauvage FJ. Adrenergic nerves regulate intestinal regeneration through IL-22 signaling from type 3 innate lymphoid cells. Cell Stem Cell 2023; 30:1166-1178.e8. [PMID: 37597516 DOI: 10.1016/j.stem.2023.07.013] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2022] [Revised: 06/23/2023] [Accepted: 07/25/2023] [Indexed: 08/21/2023]
Abstract
The intestinal epithelium has high intrinsic turnover rate, and the precise renewal of the epithelium is dependent on the microenvironment. The intestine is innervated by a dense network of peripheral nerves that controls various aspects of intestinal physiology. However, the role of neurons in regulating epithelial cell regeneration remains largely unknown. Here, we investigated the effects of gut-innervating adrenergic nerves on epithelial cell repair following irradiation (IR)-induced injury. We observed that adrenergic nerve density in the small intestine increased post IR, while chemical adrenergic denervation impaired epithelial regeneration. Single-cell RNA sequencing experiments revealed a decrease in IL-22 signaling post IR in denervated animals. Combining pharmacologic and genetic tools, we demonstrate that β-adrenergic receptor signaling drives IL-22 production from type 3 innate lymphoid cells (ILC3s) post IR, which in turn promotes epithelial regeneration. These results define an adrenergic-ILC3 axis important for intestinal regeneration.
Collapse
Affiliation(s)
- Putianqi Wang
- Molecular Oncology, Genentech, Inc., 1 DNA Way, South San Francisco, CA 94080, USA
| | - Noelyn Kljavin
- Molecular Oncology, Genentech, Inc., 1 DNA Way, South San Francisco, CA 94080, USA
| | - Thi Thu Thao Nguyen
- Oncology Bioinformatics, Genentech, Inc., 1 DNA Way, South San Francisco, CA 94080, USA
| | - Elaine E Storm
- Molecular Oncology, Genentech, Inc., 1 DNA Way, South San Francisco, CA 94080, USA
| | - Bryan Marsh
- Molecular Oncology, Genentech, Inc., 1 DNA Way, South San Francisco, CA 94080, USA
| | - Jian Jiang
- Research Pathology, Genentech, Inc., 1 DNA Way, South San Francisco, CA 94080, USA
| | - William Lin
- Research Pathology, Genentech, Inc., 1 DNA Way, South San Francisco, CA 94080, USA
| | - Hari Menon
- Microchemistry, Proteomics, Lipidomics and Next Generation Sequencing, Genentech, Inc., 1 DNA Way, South San Francisco, CA 94080, USA
| | - Robert Piskol
- Oncology Bioinformatics, Genentech, Inc., 1 DNA Way, South San Francisco, CA 94080, USA
| | | |
Collapse
|
17
|
Smith MH, Gao VR, Periyakoil PK, Kochen A, DiCarlo EF, Goodman SM, Norman TM, Donlin LT, Leslie CS, Rudensky AY. Drivers of heterogeneity in synovial fibroblasts in rheumatoid arthritis. Nat Immunol 2023; 24:1200-1210. [PMID: 37277655 PMCID: PMC10307631 DOI: 10.1038/s41590-023-01527-9] [Citation(s) in RCA: 19] [Impact Index Per Article: 19.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/26/2022] [Accepted: 05/04/2023] [Indexed: 06/07/2023]
Abstract
Inflammation of non-barrier immunologically quiescent tissues is associated with a massive influx of blood-borne innate and adaptive immune cells. Cues from the latter are likely to alter and expand activated states of the resident cells. However, local communications between immigrant and resident cell types in human inflammatory disease remain poorly understood. Here, we explored drivers of fibroblast-like synoviocyte (FLS) heterogeneity in inflamed joints of patients with rheumatoid arthritis using paired single-cell RNA and ATAC sequencing, multiplexed imaging and spatial transcriptomics along with in vitro modeling of cell-extrinsic factor signaling. These analyses suggest that local exposures to myeloid and T cell-derived cytokines, TNF, IFN-γ, IL-1β or lack thereof, drive four distinct FLS states some of which closely resemble fibroblast states in other disease-affected tissues including skin and colon. Our results highlight a role for concurrent, spatially distributed cytokine signaling within the inflamed synovium.
Collapse
Affiliation(s)
- Melanie H Smith
- Division of Rheumatology, Department of Medicine, Hospital for Special Surgery, New York, NY, USA.
- Howard Hughes Medical Institute and Immunology Program at Sloan Kettering Institute, Ludwig Center for Cancer Immunotherapy, Memorial Sloan Kettering Cancer Center, New York, NY, USA.
| | - Vianne R Gao
- Computational and Systems Biology Program, Memorial Sloan Kettering Cancer Center, New York, NY, USA
- Weill Cornell Medical College and Graduate School, New York, NY, USA
| | - Preethi K Periyakoil
- Computational and Systems Biology Program, Memorial Sloan Kettering Cancer Center, New York, NY, USA
| | - Alejandro Kochen
- Arthritis and Tissue Degeneration Program and the David Z. Rosensweig Genomics Research Center, Hospital for Special Surgery, New York, NY, USA
| | - Edward F DiCarlo
- Department of Pathology and Laboratory Medicine, Hospital for Special Surgery, New York, NY, USA
| | - Susan M Goodman
- Division of Rheumatology, Department of Medicine, Hospital for Special Surgery, New York, NY, USA
- Weill Cornell Medical College and Graduate School, New York, NY, USA
| | - Thomas M Norman
- Computational and Systems Biology Program, Memorial Sloan Kettering Cancer Center, New York, NY, USA
| | - Laura T Donlin
- Weill Cornell Medical College and Graduate School, New York, NY, USA
- Arthritis and Tissue Degeneration Program and the David Z. Rosensweig Genomics Research Center, Hospital for Special Surgery, New York, NY, USA
| | - Christina S Leslie
- Computational and Systems Biology Program, Memorial Sloan Kettering Cancer Center, New York, NY, USA.
| | - Alexander Y Rudensky
- Howard Hughes Medical Institute and Immunology Program at Sloan Kettering Institute, Ludwig Center for Cancer Immunotherapy, Memorial Sloan Kettering Cancer Center, New York, NY, USA.
| |
Collapse
|
18
|
Cain A, Taga M, McCabe C, Green GS, Hekselman I, White CC, Lee DI, Gaur P, Rozenblatt-Rosen O, Zhang F, Yeger-Lotem E, Bennett DA, Yang HS, Regev A, Menon V, Habib N, De Jager PL. Multicellular communities are perturbed in the aging human brain and Alzheimer's disease. Nat Neurosci 2023; 26:1267-1280. [PMID: 37336975 PMCID: PMC10789499 DOI: 10.1038/s41593-023-01356-x] [Citation(s) in RCA: 23] [Impact Index Per Article: 23.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/23/2020] [Accepted: 05/10/2023] [Indexed: 06/21/2023]
Abstract
The role of different cell types and their interactions in Alzheimer's disease (AD) is a complex and open question. Here, we pursued this question by assembling a high-resolution cellular map of the aging frontal cortex using single-nucleus RNA sequencing of 24 individuals with a range of clinicopathologic characteristics. We used this map to infer the neocortical cellular architecture of 638 individuals profiled by bulk RNA sequencing, providing the sample size necessary for identifying statistically robust associations. We uncovered diverse cell populations associated with AD, including a somatostatin inhibitory neuronal subtype and oligodendroglial states. We further identified a network of multicellular communities, each composed of coordinated subpopulations of neuronal, glial and endothelial cells, and we found that two of these communities are altered in AD. Finally, we used mediation analyses to prioritize cellular changes that might contribute to cognitive decline. Thus, our deconstruction of the aging neocortex provides a roadmap for evaluating the cellular microenvironments underlying AD and dementia.
Collapse
Affiliation(s)
- Anael Cain
- Edmond & Lily Safra Center for Brain Sciences, The Hebrew University of Jerusalem, Jerusalem, Israel
| | - Mariko Taga
- Center for Translational & Computational Immunology, Department of Neurology and Taub Institute for Research on Alzheimer's Disease and the Aging Brain, Columbia University Irving Medical Center, New York, NY, USA
| | - Cristin McCabe
- Klarman Cell Observatory, Broad Institute of MIT and Harvard, Cambridge, MA, USA
| | - Gilad S Green
- Edmond & Lily Safra Center for Brain Sciences, The Hebrew University of Jerusalem, Jerusalem, Israel
| | - Idan Hekselman
- Department of Clinical Biochemistry and Pharmacology, Faculty of Health Sciences, Ben-Gurion University of the Negev, Beer-Sheva, Israel
| | | | - Dylan I Lee
- Center for Translational & Computational Immunology, Department of Neurology and Taub Institute for Research on Alzheimer's Disease and the Aging Brain, Columbia University Irving Medical Center, New York, NY, USA
| | - Pallavi Gaur
- Center for Translational & Computational Immunology, Department of Neurology and Taub Institute for Research on Alzheimer's Disease and the Aging Brain, Columbia University Irving Medical Center, New York, NY, USA
| | - Orit Rozenblatt-Rosen
- Klarman Cell Observatory, Broad Institute of MIT and Harvard, Cambridge, MA, USA
- Genentech, South San Francisco, CA, USA
| | - Feng Zhang
- Broad Institute, Cambridge, MA, USA
- Department of Brain and Cognitive Sciences, McGovern Institute for Brain Research, Massachusetts Institute of Technology, Cambridge, MA, USA
| | - Esti Yeger-Lotem
- Department of Clinical Biochemistry and Pharmacology, Faculty of Health Sciences, Ben-Gurion University of the Negev, Beer-Sheva, Israel
- National Institute for Biotechnology in the Negev, Ben-Gurion University of the Negev, Beer-Sheva, Israel
| | - David A Bennett
- Rush Alzheimer's Disease Center, Rush University Medical Center, Chicago, IL, USA
| | - Hyun-Sik Yang
- Broad Institute, Cambridge, MA, USA
- Harvard Medical School, Boston, MA, USA
- Center for Alzheimer Research and Treatment, Department of Neurology, Brigham and Women's Hospital, Boston, MA, USA
| | - Aviv Regev
- Klarman Cell Observatory, Broad Institute of MIT and Harvard, Cambridge, MA, USA
- Department of Biology, Koch Institute of Integrative Cancer Research, Massachusetts Institute of Technology, Cambridge, MA, USA
- Genentech, South San Francisco, CA, USA
| | - Vilas Menon
- Center for Translational & Computational Immunology, Department of Neurology and Taub Institute for Research on Alzheimer's Disease and the Aging Brain, Columbia University Irving Medical Center, New York, NY, USA.
| | - Naomi Habib
- Edmond & Lily Safra Center for Brain Sciences, The Hebrew University of Jerusalem, Jerusalem, Israel.
| | - Philip L De Jager
- Center for Translational & Computational Immunology, Department of Neurology and Taub Institute for Research on Alzheimer's Disease and the Aging Brain, Columbia University Irving Medical Center, New York, NY, USA.
- Broad Institute, Cambridge, MA, USA.
| |
Collapse
|
19
|
Gao CF, Vaikuntanathan S, Riesenfeld SJ. Dissection and Integration of Bursty Transcriptional Dynamics for Complex Systems. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.06.13.544828. [PMID: 37398022 PMCID: PMC10312759 DOI: 10.1101/2023.06.13.544828] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 07/04/2023]
Abstract
RNA velocity estimation is a potentially powerful tool to reveal the directionality of transcriptional changes in single-cell RNA-seq data, but it lacks accuracy, absent advanced metabolic labeling techniques. We developed a novel approach, TopicVelo, that disentangles simultaneous, yet distinct, dynamics by using a probabilistic topic model, a highly interpretable form of latent space factorization, to infer cells and genes associated with individual processes, thereby capturing cellular pluripotency or multifaceted functionality. Focusing on process-associated cells and genes enables accurate estimation of process-specific velocities via a master equation for a transcriptional burst model accounting for intrinsic stochasticity. The method obtains a global transition matrix by leveraging cell topic weights to integrate process-specific signals. In challenging systems, this method accurately recovers complex transitions and terminal states, while our novel use of first-passage time analysis provides insights into transient transitions. These results expand the limits of RNA velocity, empowering future studies of cell fate and functional responses.
Collapse
Affiliation(s)
| | | | - Samantha J Riesenfeld
- Institute for Biophysical Dynamics, University of Chicago, IL
- Pritzker School of Molecular Engineering, University of Chicago, IL
- Department of Medicine, University of Chicago, IL
- Committee on Immunology, University of Chicago, IL
| |
Collapse
|
20
|
Min A, Durham T, Gevirtzman L, Noble WS. Matrix prior for data transfer between single cell data types in latent Dirichlet allocation. PLoS Comput Biol 2023; 19:e1011049. [PMID: 37146053 DOI: 10.1371/journal.pcbi.1011049] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/22/2022] [Revised: 05/17/2023] [Accepted: 03/26/2023] [Indexed: 05/07/2023] Open
Abstract
Single cell ATAC-seq (scATAC-seq) enables the mapping of regulatory elements in fine-grained cell types. Despite this advance, analysis of the resulting data is challenging, and large scale scATAC-seq data are difficult to obtain and expensive to generate. This motivates a method to leverage information from previously generated large scale scATAC-seq or scRNA-seq data to guide our analysis of new scATAC-seq datasets. We analyze scATAC-seq data using latent Dirichlet allocation (LDA), a Bayesian algorithm that was developed to model text corpora, summarizing documents as mixtures of topics defined based on the words that distinguish the documents. When applied to scATAC-seq, LDA treats cells as documents and their accessible sites as words, identifying "topics" based on the cell type-specific accessible sites in those cells. Previous work used uniform symmetric priors in LDA, but we hypothesized that nonuniform matrix priors generated from LDA models trained on existing data sets may enable improved detection of cell types in new data sets, especially if they have relatively few cells. In this work, we test this hypothesis in scATAC-seq data from whole C. elegans nematodes and SHARE-seq data from mouse skin cells. We show that nonsymmetric matrix priors for LDA improve our ability to capture cell type information from small scATAC-seq datasets.
Collapse
Affiliation(s)
- Alan Min
- Department of Statistics, University of Washington, Seattle, Washington, United States of America
| | - Timothy Durham
- Department of Genomics, University of Washington, Seattle, Washington, United States of America
- Broad Institute, Cambridge, Massachusetts, United States of America
| | - Louis Gevirtzman
- Department of Genomics, University of Washington, Seattle, Washington, United States of America
| | - William Stafford Noble
- Department of Genomics, University of Washington, Seattle, Washington, United States of America
- Paul G. Allen School of Computer Science and Engineering, University of Washington, Seattle, Washington, United States of America
| |
Collapse
|
21
|
McCabe SD, Nobel AB, Love MI. ACTOR: a latent Dirichlet model to compare expressed isoform proportions to a reference panel. Biostatistics 2023; 24:388-405. [PMID: 33948626 PMCID: PMC10102900 DOI: 10.1093/biostatistics/kxab013] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/23/2020] [Revised: 03/19/2021] [Accepted: 03/23/2021] [Indexed: 11/13/2022] Open
Abstract
The relative proportion of RNA isoforms expressed for a given gene has been associated with disease states in cancer, retinal diseases, and neurological disorders. Examination of relative isoform proportions can help determine biological mechanisms, but such analyses often require a per-gene investigation of splicing patterns. Leveraging large public data sets produced by genomic consortia as a reference, one can compare splicing patterns in a data set of interest with those of a reference panel in which samples are divided into distinct groups, such as tissue of origin, or disease status. We propose A latent Dirichlet model to Compare expressed isoform proportions TO a Reference panel (ACTOR), a latent Dirichlet model with Dirichlet Multinomial observations to compare expressed isoform proportions in a data set to an independent reference panel. We use a variational Bayes procedure to estimate posterior distributions for the group membership of one or more samples. Using the Genotype-Tissue Expression project as a reference data set, we evaluate ACTOR on simulated and real RNA-seq data sets to determine tissue-type classifications of genes. ACTOR is publicly available as an R package at https://github.com/mccabes292/actor.
Collapse
Affiliation(s)
- Sean D McCabe
- Department of Biostatistics, University of North Carolina at Chapel Hill, 135 Dauer Drive, Chapel Hill, NC 27599-7400, USA
| | - Andrew B Nobel
- Department of Statistics and Operations Research, University of North Carolina at Chapel Hill, 318 Hanes Hall, Chapel Hill, NC 27599-3260, USA and Department of Biostatistics, University of North Carolina at Chapel Hill, 135 Dauer Drive, Chapel Hill, NC 27599-7400, USA
| | - Michael I Love
- Department of Biostatistics, University of North Carolina at Chapel Hill, 135 Dauer Drive, Chapel Hill, NC 27599-7400, USA and Department of Genetics, University of North Carolina at Chapel Hill, 120 Mason Farm Rd, Chapel Hill, NC 27514, USA
| |
Collapse
|
22
|
Himmels P, Nguyen TTT, Mitzner MC, Arrazate A, Yeung S, Burton J, Clark R, Totpal K, Jesudason R, Yang A, Solon M, Eastham J, Modrusan Z, Webster JD, Lo AA, Piskol R, Ye W. T cell-dependent bispecific antibodies alter organ-specific endothelial cell-T cell interaction. EMBO Rep 2023; 24:e55532. [PMID: 36621885 PMCID: PMC9986820 DOI: 10.15252/embr.202255532] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/02/2022] [Revised: 12/08/2022] [Accepted: 12/09/2022] [Indexed: 01/10/2023] Open
Abstract
Preclinical and clinical studies demonstrate that T cell-dependent bispecific antibodies (TDBs) induce systemic changes in addition to tumor killing, leading to adverse events. Here, we report an in-depth characterization of acute responses to TDBs in tumor-bearing mice. Contrary to modest changes in tumors, rapid and substantial lymphocyte accumulation and endothelial cell (EC) activation occur around large blood vessels in normal organs including the liver. We hypothesize that organ-specific ECs may account for the differential responses in normal tissues and tumors, and we identify a list of genes selectively upregulated by TDB in large liver vessels. Using one of the genes as an example, we demonstrate that CD9 facilitates ICAM-1 to support T cell-EC interaction in response to soluble factors released from a TDB-mediated cytotoxic reaction. Our results suggest that multiple factors may cooperatively promote T cell infiltration into normal organs as a secondary response to TDB-mediated tumor killing. These data shed light on how different vascular beds respond to cancer immunotherapy and may help improve their safety and efficacy.
Collapse
Affiliation(s)
- Patricia Himmels
- Department of Molecular OncologyGenentechSouth San FranciscoCAUSA
| | | | - Maresa Caunt Mitzner
- Department of Molecular OncologyGenentechSouth San FranciscoCAUSA
- Product DevelopmentGenentechSouth San FranciscoCAUSA
| | - Alfonso Arrazate
- Department of Translational OncologyGenentechSouth San FranciscoCAUSA
| | - Stacey Yeung
- Department of Molecular OncologyGenentechSouth San FranciscoCAUSA
| | - Jeremy Burton
- Department of Molecular OncologyGenentechSouth San FranciscoCAUSA
| | - Robyn Clark
- Department of Translational OncologyGenentechSouth San FranciscoCAUSA
| | - Klara Totpal
- Department of Translational OncologyGenentechSouth San FranciscoCAUSA
| | - Raj Jesudason
- Department of Research PathologyGenentechSouth San FranciscoCAUSA
| | - Angela Yang
- GSK‐Laboratory for Genomic ResearchSan FranciscoCAUSA
- Department of Microchemistry, Proteomics and Lipidomics, and Next Generation Sequencing (MPL‐NGS)GenentechSouth San FranciscoCAUSA
| | - Margaret Solon
- Department of Research PathologyGenentechSouth San FranciscoCAUSA
| | - Jeffrey Eastham
- Department of Research PathologyGenentechSouth San FranciscoCAUSA
| | - Zora Modrusan
- Department of Microchemistry, Proteomics and Lipidomics, and Next Generation Sequencing (MPL‐NGS)GenentechSouth San FranciscoCAUSA
| | - Joshua D Webster
- Department of Research PathologyGenentechSouth San FranciscoCAUSA
| | - Amy A Lo
- Department of Research PathologyGenentechSouth San FranciscoCAUSA
| | - Robert Piskol
- Department of Oncology BioinformaticsGenentechSouth San FranciscoCAUSA
| | - Weilan Ye
- Department of Molecular OncologyGenentechSouth San FranciscoCAUSA
| |
Collapse
|
23
|
El Hachem EJ, Sokolovska N, Soula H. Latent dirichlet allocation for double clustering (LDA-DC): discovering patients phenotypes and cell populations within a single Bayesian framework. BMC Bioinformatics 2023; 24:61. [PMID: 36823548 PMCID: PMC9948385 DOI: 10.1186/s12859-023-05177-4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/25/2022] [Accepted: 02/08/2023] [Indexed: 02/25/2023] Open
Abstract
BACKGROUND Current clinical routines rely more and more on "omics" data such as flow cytometry data from host and microbiota. Cohorts variability in addition to patients' heterogeneity and huge dimensions make it difficult to understand underlying structure of the data and decipher pathologies. Patients stratification and diagnostics from such complex data are extremely challenging. There is an acute need to develop novel statistical machine learning methods that are robust with respect to the data heterogeneity, efficient from the computational viewpoint, and can be understood by human experts. RESULTS We propose a novel approach to stratify cell-based observations within a single probabilistic framework, i.e., to extract meaningful phenotypes from both patients and cells data simultaneously. We define this problem as a double clustering problem that we tackle with the proposed approach. Our method is a practical extension of the Latent Dirichlet Allocation and is used for the Double Clustering task (LDA-DC). We first validate the method on artificial datasets, then we apply our method to two real problems of patients stratification based on cytometry and microbiota data. We observe that the LDA-DC returns clusters of patients and also clusters of cells related to patients' conditions. We also construct a graphical representation of the results that can be easily understood by humans and are, therefore, of a big help for experts involved in pre-clinical research.
Collapse
Affiliation(s)
- Elie-Julien El Hachem
- Sorbonne University, INSERM, Nutrition and Obesities: Systemic Approaches, NutriOmique, 91 Boulevard de l'hôpital, 75013, Paris, France.
| | - Nataliya Sokolovska
- Sorbonne University, INSERM, Nutrition and Obesities: Systemic Approaches, NutriOmique, 91 Boulevard de l'hôpital, 75013, Paris, France
| | - Hedi Soula
- Sorbonne University, INSERM, Nutrition and Obesities: Systemic Approaches, NutriOmique, 91 Boulevard de l'hôpital, 75013, Paris, France
| |
Collapse
|
24
|
Molecular and cellular evolution of the amygdala across species analyzed by single-nucleus transcriptome profiling. Cell Discov 2023; 9:19. [PMID: 36788214 PMCID: PMC9929086 DOI: 10.1038/s41421-022-00506-y] [Citation(s) in RCA: 13] [Impact Index Per Article: 13.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/18/2022] [Accepted: 11/24/2022] [Indexed: 02/16/2023] Open
Abstract
The amygdala, or an amygdala-like structure, is found in the brains of all vertebrates and plays a critical role in survival and reproduction. However, the cellular architecture of the amygdala and how it has evolved remain elusive. Here, we generated single-nucleus RNA-sequencing data for more than 200,000 cells in the amygdala of humans, macaques, mice, and chickens. Abundant neuronal cell types from different amygdala subnuclei were identified in all datasets. Cross-species analysis revealed that inhibitory neurons and inhibitory neuron-enriched subnuclei of the amygdala were well-conserved in cellular composition and marker gene expression, whereas excitatory neuron-enriched subnuclei were relatively divergent. Furthermore, LAMP5+ interneurons were much more abundant in primates, while DRD2+ inhibitory neurons and LAMP5+SATB2+ excitatory neurons were dominant in the human central amygdalar nucleus (CEA) and basolateral amygdalar complex (BLA), respectively. We also identified CEA-like neurons and their species-specific distribution patterns in chickens. This study highlights the extreme cell-type diversity in the amygdala and reveals the conservation and divergence of cell types and gene expression patterns across species that may contribute to species-specific adaptations.
Collapse
|
25
|
Takahama M, Patil A, Johnson K, Cipurko D, Miki Y, Taketomi Y, Carbonetto P, Plaster M, Richey G, Pandey S, Cheronis K, Ueda T, Gruenbaum A, Dudek SM, Stephens M, Murakami M, Chevrier N. Organism-Wide Analysis of Sepsis Reveals Mechanisms of Systemic Inflammation. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.01.30.526342. [PMID: 36778287 PMCID: PMC9915512 DOI: 10.1101/2023.01.30.526342] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/05/2023]
Abstract
Sepsis is a systemic response to infection with life-threatening consequences. Our understanding of the impact of sepsis across organs of the body is rudimentary. Here, using mouse models of sepsis, we generate a dynamic, organism-wide map of the pathogenesis of the disease, revealing the spatiotemporal patterns of the effects of sepsis across tissues. These data revealed two interorgan mechanisms key in sepsis. First, we discover a simplifying principle in the systemic behavior of the cytokine network during sepsis, whereby a hierarchical cytokine circuit arising from the pairwise effects of TNF plus IL-18, IFN-γ, or IL-1β explains half of all the cellular effects of sepsis on 195 cell types across 9 organs. Second, we find that the secreted phospholipase PLA2G5 mediates hemolysis in blood, contributing to organ failure during sepsis. These results provide fundamental insights to help build a unifying mechanistic framework for the pathophysiological effects of sepsis on the body.
Collapse
|
26
|
Gassler J, Kobayashi W, Gáspár I, Ruangroengkulrith S, Mohanan A, Gómez Hernández L, Kravchenko P, Kümmecke M, Lalic A, Rifel N, Ashburn RJ, Zaczek M, Vallot A, Cuenca Rico L, Ladstätter S, Tachibana K. Zygotic genome activation by the totipotency pioneer factor Nr5a2. Science 2022; 378:1305-1315. [PMID: 36423263 DOI: 10.1126/science.abn7478] [Citation(s) in RCA: 48] [Impact Index Per Article: 24.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
Abstract
Life begins with a switch in genetic control from the maternal to the embryonic genome during zygotic genome activation (ZGA). Despite its importance, the essential regulators of ZGA remain largely unknown in mammals. On the basis of de novo motif searches, we identified the orphan nuclear receptor Nr5a2 as a key activator of major ZGA in mouse two-cell embryos. Nr5a2 is required for progression beyond the two-cell stage. It binds to its motif within SINE B1/Alu retrotransposable elements found in cis-regulatory regions of ZGA genes. Chemical inhibition suggests that 72% of ZGA genes are regulated by Nr5a2 and potentially other orphan nuclear receptors. Nr5a2 promotes chromatin accessibility during ZGA and binds nucleosomal DNA in vitro. We conclude that Nr5a2 is an essential pioneer factor that regulates ZGA.
Collapse
Affiliation(s)
- Johanna Gassler
- Department of Totipotency, Max Planck Institute of Biochemistry (MPIB), Munich, Germany.,Institute of Molecular Biotechnology of the Austrian Academy of Sciences (IMBA), Vienna Biocenter (VBC), Vienna, Austria
| | - Wataru Kobayashi
- Department of Totipotency, Max Planck Institute of Biochemistry (MPIB), Munich, Germany
| | - Imre Gáspár
- Department of Totipotency, Max Planck Institute of Biochemistry (MPIB), Munich, Germany
| | | | - Adarsh Mohanan
- Department of Totipotency, Max Planck Institute of Biochemistry (MPIB), Munich, Germany
| | - Laura Gómez Hernández
- Department of Totipotency, Max Planck Institute of Biochemistry (MPIB), Munich, Germany
| | - Pavel Kravchenko
- Department of Totipotency, Max Planck Institute of Biochemistry (MPIB), Munich, Germany
| | - Maximilian Kümmecke
- Institute of Molecular Biotechnology of the Austrian Academy of Sciences (IMBA), Vienna Biocenter (VBC), Vienna, Austria
| | - Aleksandar Lalic
- Department of Totipotency, Max Planck Institute of Biochemistry (MPIB), Munich, Germany
| | - Nikita Rifel
- Department of Totipotency, Max Planck Institute of Biochemistry (MPIB), Munich, Germany
| | - Robert John Ashburn
- Department of Totipotency, Max Planck Institute of Biochemistry (MPIB), Munich, Germany
| | - Maciej Zaczek
- Institute of Molecular Biotechnology of the Austrian Academy of Sciences (IMBA), Vienna Biocenter (VBC), Vienna, Austria
| | - Antoine Vallot
- Institute of Molecular Biotechnology of the Austrian Academy of Sciences (IMBA), Vienna Biocenter (VBC), Vienna, Austria
| | - Laura Cuenca Rico
- Institute of Molecular Biotechnology of the Austrian Academy of Sciences (IMBA), Vienna Biocenter (VBC), Vienna, Austria
| | - Sabrina Ladstätter
- Institute of Molecular Biotechnology of the Austrian Academy of Sciences (IMBA), Vienna Biocenter (VBC), Vienna, Austria
| | - Kikuë Tachibana
- Department of Totipotency, Max Planck Institute of Biochemistry (MPIB), Munich, Germany.,Institute of Molecular Biotechnology of the Austrian Academy of Sciences (IMBA), Vienna Biocenter (VBC), Vienna, Austria
| |
Collapse
|
27
|
Bastide S, Chomsky E, Saudemont B, Loe-Mie Y, Schmutz S, Novault S, Marlow H, Tanay A, Spitz F. TATTOO-seq delineates spatial and cell type-specific regulatory programs in the developing limb. SCIENCE ADVANCES 2022; 8:eadd0695. [PMID: 36516250 PMCID: PMC9750149 DOI: 10.1126/sciadv.add0695] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 05/18/2022] [Accepted: 11/14/2022] [Indexed: 06/17/2023]
Abstract
The coordinated differentiation of progenitor cells into specialized cell types and their spatial organization into distinct domains is central to embryogenesis. Here, we developed and applied an unbiased spatially resolved single-cell transcriptomics method to identify the genetic programs underlying the emergence of specialized cell types during mouse limb development and their spatial integration. We identify multiple transcription factors whose expression patterns are predominantly associated with cell type specification or spatial position, suggesting two parallel yet highly interconnected regulatory systems. We demonstrate that the embryonic limb undergoes a complex multiscale reorganization upon perturbation of one of its spatial organizing centers, including the loss of specific cell populations, alterations of preexisting cell states' molecular identities, and changes in their relative spatial distribution. Our study shows how multidimensional single-cell, spatially resolved molecular atlases can allow the deconvolution of spatial identity and cell fate and reveal the interconnected genetic networks that regulate organogenesis and its reorganization upon genetic alterations.
Collapse
Affiliation(s)
- Sébastien Bastide
- (Epi)genomics of Animal Development, Department of Developmental and Stem Cell Biology, Institut Pasteur, Paris, France
- École Doctorale “Complexité du Vivant”, Sorbonne Université, 75005 Paris, France
- Department of Human Genetics, The University of Chicago, Chicago, IL, USA
| | - Elad Chomsky
- Department of Computer Science and Applied Mathematics, Weizmann Institute, Rehovot, Israel
- Department of Biological Regulation, Weizmann Institute, Rehovot, Israel
| | - Baptiste Saudemont
- (Epi)genomics of Animal Development, Department of Developmental and Stem Cell Biology, Institut Pasteur, Paris, France
| | - Yann Loe-Mie
- (Epi)genomics of Animal Development, Department of Developmental and Stem Cell Biology, Institut Pasteur, Paris, France
- Hub de Bioinformatique et Biostatistique, Département Biologie Computationnelle, Institut Pasteur, Paris, France
| | - Sandrine Schmutz
- Cytometry and Biomarkers, Center for Technological Resources and Research, Institut Pasteur, Paris, France
| | - Sophie Novault
- Cytometry and Biomarkers, Center for Technological Resources and Research, Institut Pasteur, Paris, France
| | - Heather Marlow
- (Epi)genomics of Animal Development, Department of Developmental and Stem Cell Biology, Institut Pasteur, Paris, France
- Department of Organismal Biology and Anatomy, The University of Chicago, Chicago, IL, USA
| | - Amos Tanay
- Department of Computer Science and Applied Mathematics, Weizmann Institute, Rehovot, Israel
| | - François Spitz
- (Epi)genomics of Animal Development, Department of Developmental and Stem Cell Biology, Institut Pasteur, Paris, France
- Department of Human Genetics, The University of Chicago, Chicago, IL, USA
| |
Collapse
|
28
|
Ors A, Chitsazan AD, Doe AR, Mulqueen RM, Ak C, Wen Y, Haverlack S, Handu M, Naldiga S, Saldivar J, Mohammed H. Estrogen regulates divergent transcriptional and epigenetic cell states in breast cancer. Nucleic Acids Res 2022; 50:11492-11508. [PMID: 36318267 PMCID: PMC9723652 DOI: 10.1093/nar/gkac908] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/24/2022] [Revised: 09/20/2022] [Accepted: 10/20/2022] [Indexed: 12/12/2022] Open
Abstract
Breast cancers are known to be driven by the transcription factor estrogen receptor and its ligand estrogen. While the receptor's cis-binding elements are known to vary between tumors, heterogeneity of hormone signaling at a single-cell level is unknown. In this study, we systematically tracked estrogen response across time at a single-cell level in multiple cell line and organoid models. To accurately model these changes, we developed a computational tool (TITAN) that quantifies signaling gradients in single-cell datasets. Using this approach, we found that gene expression response to estrogen is non-uniform, with distinct cell groups expressing divergent transcriptional networks. Pathway analysis suggested the two most distinct signatures are driven separately by ER and FOXM1. We observed that FOXM1 was indeed activated by phosphorylation upon estrogen stimulation and silencing of FOXM1 attenuated the relevant gene signature. Analysis of scRNA-seq data from patient samples confirmed the existence of these divergent cell groups, with the FOXM1 signature predominantly found in ER negative cells. Further, multi-omic single-cell experiments indicated that the different cell groups have distinct chromatin accessibility states. Our results provide a comprehensive insight into ER biology at the single-cell level and potential therapeutic strategies to mitigate resistance to therapy.
Collapse
Affiliation(s)
| | | | - Aaron Reid Doe
- Cancer Early Detection Advanced Research Center, Knight Cancer Institute, Oregon Health & Science University, Portland, Oregon 97201, USA
| | - Ryan M Mulqueen
- Cancer Early Detection Advanced Research Center, Knight Cancer Institute, Oregon Health & Science University, Portland, Oregon 97201, USA
| | - Cigdem Ak
- Cancer Early Detection Advanced Research Center, Knight Cancer Institute, Oregon Health & Science University, Portland, Oregon 97201, USA
| | - Yahong Wen
- Cancer Early Detection Advanced Research Center, Knight Cancer Institute, Oregon Health & Science University, Portland, Oregon 97201, USA
| | - Syber Haverlack
- Cancer Early Detection Advanced Research Center, Knight Cancer Institute, Oregon Health & Science University, Portland, Oregon 97201, USA
| | - Mithila Handu
- Cancer Early Detection Advanced Research Center, Knight Cancer Institute, Oregon Health & Science University, Portland, Oregon 97201, USA
| | - Spandana Naldiga
- Cancer Early Detection Advanced Research Center, Knight Cancer Institute, Oregon Health & Science University, Portland, Oregon 97201, USA
| | - Joshua C Saldivar
- Cancer Early Detection Advanced Research Center, Knight Cancer Institute, Oregon Health & Science University, Portland, Oregon 97201, USA,Division of Oncological Sciences, Knight Cancer Institute, Oregon Health & Science University, Portland, Oregon 97201, USA
| | | |
Collapse
|
29
|
Quintanal-Villalonga Á, Chan JM, Masilionis I, Gao VR, Xie Y, Allaj V, Chow A, Poirier JT, Pe'er D, Rudin CM, Mazutis L. Protocol to dissociate, process, and analyze the human lung tissue using single-cell RNA-seq. STAR Protoc 2022; 3:101776. [PMID: 36313536 PMCID: PMC9597186 DOI: 10.1016/j.xpro.2022.101776] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022] Open
Abstract
We report a protocol for obtaining high-quality single-cell transcriptomics data from human lung biospecimens acquired from core needle biopsies, fine-needle aspirates, surgical resection, and pleural effusions. The protocol relies upon the brief mechanical and enzymatic disruption of tissue, enrichment of live cells by fluorescence-activated cell sorting (FACS), and droplet-based single-cell RNA sequencing (scRNA-seq). The protocol also details a procedure for analyzing the scRNA-seq data. For complete details on the use and execution of this protocol, please refer to Chan et al. (2021).
Collapse
Affiliation(s)
- Álvaro Quintanal-Villalonga
- Department of Medicine, Thoracic Oncology Service, Memorial Sloan Kettering Cancer Center, New York, NY, USA
| | - Joseph M Chan
- Department of Medicine, Thoracic Oncology Service, Memorial Sloan Kettering Cancer Center, New York, NY, USA; Program for Computational and Systems Biology, Sloan Kettering Institute, Memorial Sloan Kettering Cancer Center, New York, NY, USA
| | - Ignas Masilionis
- Program for Computational and Systems Biology, Sloan Kettering Institute, Memorial Sloan Kettering Cancer Center, New York, NY, USA
| | - Vianne Ran Gao
- Program for Computational and Systems Biology, Sloan Kettering Institute, Memorial Sloan Kettering Cancer Center, New York, NY, USA
| | - Yubin Xie
- Program for Computational and Systems Biology, Sloan Kettering Institute, Memorial Sloan Kettering Cancer Center, New York, NY, USA
| | - Viola Allaj
- Department of Medicine, Thoracic Oncology Service, Memorial Sloan Kettering Cancer Center, New York, NY, USA
| | - Andrew Chow
- Department of Medicine, Thoracic Oncology Service, Memorial Sloan Kettering Cancer Center, New York, NY, USA
| | - John T Poirier
- Perlmutter Cancer Center, New York University Langone Health, New York, NY, USA
| | - Dana Pe'er
- Program for Computational and Systems Biology, Sloan Kettering Institute, Memorial Sloan Kettering Cancer Center, New York, NY, USA; Parker Institute for Cancer Immunotherapy, Memorial Sloan Kettering Cancer Center, New York, NY, USA
| | - Charles M Rudin
- Department of Medicine, Thoracic Oncology Service, Memorial Sloan Kettering Cancer Center, New York, NY, USA
| | - Linas Mazutis
- Institute of Biotechnology, Life Sciences Centre, Vilnius University, Vilnius, Lithuania.
| |
Collapse
|
30
|
Mahalanabis A, Turinsky A, Husic M, Christensen E, Luo P, Naidas A, Brudno M, Pugh T, Ramani A, Shooshtari P. Evaluation of Single-cell RNA-seq Clustering Algorithms on Cancer Tumor Datasets. Comput Struct Biotechnol J 2022; 20:6375-6387. [DOI: 10.1016/j.csbj.2022.10.029] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/29/2022] [Revised: 10/19/2022] [Accepted: 10/20/2022] [Indexed: 11/03/2022] Open
|
31
|
Ke ZT, Wang M. Using SVD for Topic Modeling. J Am Stat Assoc 2022. [DOI: 10.1080/01621459.2022.2123813] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/17/2022]
|
32
|
Gewirtz AD, Townes FW, Engelhardt BE. Telescoping bimodal latent Dirichlet allocation to identify expression QTLs across tissues. Life Sci Alliance 2022; 5:e202101297. [PMID: 35977827 PMCID: PMC9387650 DOI: 10.26508/lsa.202101297] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/11/2021] [Revised: 07/15/2022] [Accepted: 07/18/2022] [Indexed: 11/24/2022] Open
Abstract
Expression quantitative trait loci (eQTLs), or single-nucleotide polymorphisms that affect average gene expression levels, provide important insights into context-specific gene regulation. Classic eQTL analyses use one-to-one association tests, which test gene-variant pairs individually and ignore correlations induced by gene regulatory networks and linkage disequilibrium. Probabilistic topic models, such as latent Dirichlet allocation, estimate latent topics for a collection of count observations. Prior multimodal frameworks that bridge genotype and expression data assume matched sample numbers between modalities. However, many data sets have a nested structure where one individual has several associated gene expression samples and a single germline genotype vector. Here, we build a telescoping bimodal latent Dirichlet allocation (TBLDA) framework to learn shared topics across gene expression and genotype data that allows multiple RNA sequencing samples to correspond to a single individual's genotype. By using raw count data, our model avoids possible adulteration via normalization procedures. Ancestral structure is captured in a genotype-specific latent space, effectively removing it from shared components. Using GTEx v8 expression data across 10 tissues and genotype data, we show that the estimated topics capture meaningful and robust biological signal in both modalities and identify associations within and across tissue types. We identify 4,645 cis-eQTLs and 995 trans-eQTLs by conducting eQTL mapping between the most informative features in each topic. Our TBLDA model is able to identify associations using raw sequencing count data when the samples in two separate data modalities are matched one-to-many, as is often the case in biological data. Our code is freely available at https://github.com/gewirtz/TBLDA.
Collapse
Affiliation(s)
- Ariel Dh Gewirtz
- Lewis-Sigler Institute of Integrative Genomics, Princeton University, Princeton, NJ, USA
| | - F William Townes
- Department of Computer Science, Princeton University, Princeton, NJ, USA
| | - Barbara E Engelhardt
- Department of Computer Science, Princeton University, Princeton, NJ, USA
- Gladstone Institutes, San Francisco, CA, USA
| |
Collapse
|
33
|
Kaur G, Porter CBM, Ashenberg O, Lee J, Riesenfeld SJ, Hofree M, Aggelakopoulou M, Subramanian A, Kuttikkatte SB, Attfield KE, Desel CAE, Davies JL, Evans HG, Avraham-Davidi I, Nguyen LT, Dionne DA, Neumann AE, Jensen LT, Barber TR, Soilleux E, Carrington M, McVean G, Rozenblatt-Rosen O, Regev A, Fugger L. Mouse fetal growth restriction through parental and fetal immune gene variation and intercellular communications cascade. Nat Commun 2022; 13:4398. [PMID: 35906236 PMCID: PMC9338297 DOI: 10.1038/s41467-022-32171-w] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/28/2021] [Accepted: 07/18/2022] [Indexed: 11/08/2022] Open
Abstract
Fetal growth restriction (FGR) affects 5-10% of pregnancies, and can have serious consequences for both mother and child. Prevention and treatment are limited because FGR pathogenesis is poorly understood. Genetic studies implicate KIR and HLA genes in FGR, however, linkage disequilibrium, genetic influence from both parents, and challenges with investigating human pregnancies make the risk alleles and their functional effects difficult to map. Here, we demonstrate that the interaction between the maternal KIR2DL1, expressed on uterine natural killer (NK) cells, and the paternally inherited HLA-C*0501, expressed on fetal trophoblast cells, leads to FGR in a humanized mouse model. We show that the KIR2DL1 and C*0501 interaction leads to pathogenic uterine arterial remodeling and modulation of uterine NK cell function. This initial effect cascades to altered transcriptional expression and intercellular communication at the maternal-fetal interface. These findings provide mechanistic insight into specific FGR risk alleles, and provide avenues of prevention and treatment.
Collapse
Affiliation(s)
- Gurman Kaur
- MRC Human Immunology Unit, MRC Weatherall Institute of Molecular Medicine, John Radcliffe Hospital, University of Oxford, Oxford, UK
- Klarman Cell Observatory, Broad Institute of MIT and Harvard, Cambridge, MA, USA
| | - Caroline B M Porter
- Klarman Cell Observatory, Broad Institute of MIT and Harvard, Cambridge, MA, USA
| | - Orr Ashenberg
- Klarman Cell Observatory, Broad Institute of MIT and Harvard, Cambridge, MA, USA
| | - Jack Lee
- Department of Biomedical Engineering, School of Biomedical Engineering and Imaging Sciences, King's College London, London, UK
| | - Samantha J Riesenfeld
- Klarman Cell Observatory, Broad Institute of MIT and Harvard, Cambridge, MA, USA
- Pritzker School of Molecular Engineering, University of Chicago, Chicago, IL, USA
- Department of Medicine, University of Chicago, Chicago, IL, USA
| | - Matan Hofree
- Klarman Cell Observatory, Broad Institute of MIT and Harvard, Cambridge, MA, USA
| | - Maria Aggelakopoulou
- Oxford Centre for Neuroinflammation, Nuffield Department of Clinical Neurosciences, MRC Weatherall Institute of Molecular Medicine, John Radcliffe Hospital, University of Oxford, Oxford, UK
| | | | - Subita Balaram Kuttikkatte
- MRC Human Immunology Unit, MRC Weatherall Institute of Molecular Medicine, John Radcliffe Hospital, University of Oxford, Oxford, UK
| | - Kathrine E Attfield
- Oxford Centre for Neuroinflammation, Nuffield Department of Clinical Neurosciences, MRC Weatherall Institute of Molecular Medicine, John Radcliffe Hospital, University of Oxford, Oxford, UK
| | - Christiane A E Desel
- Oxford Centre for Neuroinflammation, Nuffield Department of Clinical Neurosciences, MRC Weatherall Institute of Molecular Medicine, John Radcliffe Hospital, University of Oxford, Oxford, UK
- University Department of Neurology, University Hospital Magdeburg, Magdeburg, Germany
| | - Jessica L Davies
- Oxford Centre for Neuroinflammation, Nuffield Department of Clinical Neurosciences, MRC Weatherall Institute of Molecular Medicine, John Radcliffe Hospital, University of Oxford, Oxford, UK
| | - Hayley G Evans
- Oxford Centre for Neuroinflammation, Nuffield Department of Clinical Neurosciences, MRC Weatherall Institute of Molecular Medicine, John Radcliffe Hospital, University of Oxford, Oxford, UK
| | - Inbal Avraham-Davidi
- Klarman Cell Observatory, Broad Institute of MIT and Harvard, Cambridge, MA, USA
| | - Lan T Nguyen
- Klarman Cell Observatory, Broad Institute of MIT and Harvard, Cambridge, MA, USA
| | - Danielle A Dionne
- Klarman Cell Observatory, Broad Institute of MIT and Harvard, Cambridge, MA, USA
| | | | - Lise Torp Jensen
- Department of Clinical Medicine, Aarhus University Hospital, Aarhus, Denmark
| | - Thomas R Barber
- MRC Human Immunology Unit, MRC Weatherall Institute of Molecular Medicine, John Radcliffe Hospital, University of Oxford, Oxford, UK
| | - Elizabeth Soilleux
- Department of Pathology, Tennis Court Rd, University of Cambridge, Cambridge, England
| | - Mary Carrington
- Basic Science Program, Frederick National Laboratory for Cancer Research in the Laboratory of Integrative Cancer Immunology, National Cancer Institute, Bethesda, MD, USA
- Ragon Institute of MGH, MIT, and Harvard, Cambridge, MA, USA
| | - Gil McVean
- Big Data Institute, Li Ka Shing Centre for Health Information and Discovery, University of Oxford, Oxford, UK
| | - Orit Rozenblatt-Rosen
- Klarman Cell Observatory, Broad Institute of MIT and Harvard, Cambridge, MA, USA
- Genentech, 1 DNA Way, South San Francisco, CA, USA
| | - Aviv Regev
- Klarman Cell Observatory, Broad Institute of MIT and Harvard, Cambridge, MA, USA.
- Massachusetts Institute of Technology, Department of Biology, Cambridge, MA, USA.
- Howard Hughes Medical Institute, Chevy Chase, MD, USA.
- Genentech, 1 DNA Way, South San Francisco, CA, USA.
| | - Lars Fugger
- MRC Human Immunology Unit, MRC Weatherall Institute of Molecular Medicine, John Radcliffe Hospital, University of Oxford, Oxford, UK.
- Oxford Centre for Neuroinflammation, Nuffield Department of Clinical Neurosciences, MRC Weatherall Institute of Molecular Medicine, John Radcliffe Hospital, University of Oxford, Oxford, UK.
- Department of Clinical Medicine, Aarhus University Hospital, Aarhus, Denmark.
| |
Collapse
|
34
|
Selective activation and expansion of regulatory T cells using lipid encapsulated mRNA encoding a long-acting IL-2 mutein. Nat Commun 2022; 13:3866. [PMID: 35790728 PMCID: PMC9256694 DOI: 10.1038/s41467-022-31130-9] [Citation(s) in RCA: 20] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/08/2021] [Accepted: 05/31/2022] [Indexed: 11/30/2022] Open
Abstract
Interleukin-2 (IL-2) is critical for regulatory T cell (Treg) function and homeostasis. At low doses, IL-2 can suppress immune pathologies by expanding Tregs that constitutively express the high affinity IL-2Rα subunit. However, even low dose IL-2, signaling through the IL2-Rβ/γ complex, may lead to the activation of proinflammatory, non-Treg T cells, so improving specificity toward Tregs may be desirable. Here we use messenger RNAs (mRNA) to encode a half-life-extended human IL-2 mutein (HSA-IL2m) with mutations promoting reliance on IL-2Rα. Our data show that IL-2 mutein subcutaneous delivery as lipid-encapsulated mRNA nanoparticles selectively activates and expands Tregs in mice and non-human primates, and also reduces disease severity in mouse models of acute graft versus host disease and experimental autoimmune encephalomyelitis. Single cell RNA-sequencing of mouse splenic CD4+ T cells identifies multiple Treg states with distinct response dynamics following IL-2 mutein treatment. Our results thus demonstrate the potential of mRNA-encoded HSA-IL2m immunotherapy to treat autoimmune diseases. IL-2 has been used to expand regulatory T (Treg) cells for treating inflammatory disorders. Here the authors test an engineered IL-2 mutein, delivered subcutaneously as mRNA, to show its increased specificity for activating and expanding Treg cells in both rodents and non-human primates, and to demonstrate its ability to suppress autoimmunity in mouse models.
Collapse
|
35
|
Chen J, Liu W, Luo T, Yu Z, Jiang M, Wen J, Gupta GP, Giusti P, Zhu H, Yang Y, Li Y. A comprehensive comparison on cell-type composition inference for spatial transcriptomics data. Brief Bioinform 2022; 23:6618233. [PMID: 35753702 PMCID: PMC9294426 DOI: 10.1093/bib/bbac245] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/21/2022] [Revised: 05/20/2022] [Accepted: 05/25/2022] [Indexed: 01/11/2023] Open
Abstract
Spatial transcriptomics (ST) technologies allow researchers to examine transcriptional profiles along with maintained positional information. Such spatially resolved transcriptional characterization of intact tissue samples provides an integrated view of gene expression in its natural spatial and functional context. However, high-throughput sequencing-based ST technologies cannot yet reach single cell resolution. Thus, similar to bulk RNA-seq data, gene expression data at ST spot-level reflect transcriptional profiles of multiple cells and entail the inference of cell-type composition within each ST spot for valid and powerful subsequent analyses. Realizing the critical importance of cell-type decomposition, multiple groups have developed ST deconvolution methods. The aim of this work is to review state-of-the-art methods for ST deconvolution, comparing their strengths and weaknesses. In particular, we construct ST spots from single-cell level ST data to assess the performance of 10 methods, with either ideal reference or non-ideal reference. Furthermore, we examine the performance of these methods on spot- and bead-level ST data by comparing estimated cell-type proportions to carefully matched single-cell ST data. In comparing the performance on various tissues and technological platforms, we concluded that RCTD and stereoscope achieve more robust and accurate inferences.
Collapse
Affiliation(s)
| | | | | | - Zhentao Yu
- Department of Biostatistics, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, USA
| | - Minzhi Jiang
- Department of Applied Physical Sciences, University of North Carolina at Chapel Hill, North Carolina, USA
| | - Jia Wen
- Department of Genetics, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, USA
| | - Gaorav P Gupta
- Department of Radiation Oncology, The University of North Carolina at Chapel Hill, Chapel Hill, NC 27514, USA
| | - Paola Giusti
- Department of Psychiatry, University of Florida, Gainesville, Florida, USA
| | - Hongtu Zhu
- Department of Biostatistics, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, USA
| | - Yuchen Yang
- State Key Laboratory of Biocontrol, School of Ecology, Sun Yat-sen University, 510275 Guangzhou, China
| | - Yun Li
- Corresponding author. Yun Li, Department of Genetics, 120 Mason Farm Road, Campus Box 7264, University North Carolina, Chapel Hill, NC 27599, USA. Tel: (919) 843-2832; Fax: (919) 843-4682; E-mail:
| |
Collapse
|
36
|
Khunsriraksakul C, McGuire D, Sauteraud R, Chen F, Yang L, Wang L, Hughey J, Eckert S, Dylan Weissenkampen J, Shenoy G, Marx O, Carrel L, Jiang B, Liu DJ. Integrating 3D genomic and epigenomic data to enhance target gene discovery and drug repurposing in transcriptome-wide association studies. Nat Commun 2022; 13:3258. [PMID: 35672318 PMCID: PMC9171100 DOI: 10.1038/s41467-022-30956-7] [Citation(s) in RCA: 14] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/20/2021] [Accepted: 05/25/2022] [Indexed: 02/08/2023] Open
Abstract
Transcriptome-wide association studies (TWAS) are popular approaches to test for association between imputed gene expression levels and traits of interest. Here, we propose an integrative method PUMICE (Prediction Using Models Informed by Chromatin conformations and Epigenomics) to integrate 3D genomic and epigenomic data with expression quantitative trait loci (eQTL) to more accurately predict gene expressions. PUMICE helps define and prioritize regions that harbor cis-regulatory variants, which outperforms competing methods. We further describe an extension to our method PUMICE +, which jointly combines TWAS results from single- and multi-tissue models. Across 79 traits, PUMICE + identifies 22% more independent novel genes and increases median chi-square statistics values at known loci by 35% compared to the second-best method, as well as achieves the narrowest credible interval size. Lastly, we perform computational drug repurposing and confirm that PUMICE + outperforms other TWAS methods.
Collapse
Affiliation(s)
- Chachrit Khunsriraksakul
- grid.29857.310000 0001 2097 4281Bioinformatics and Genomics Graduate Program, Pennsylvania State University College of Medicine, Hershey, PA 17033 USA ,grid.29857.310000 0001 2097 4281Institute for Personalized Medicine, Pennsylvania State University College of Medicine, Hershey, PA 17033 USA
| | - Daniel McGuire
- grid.29857.310000 0001 2097 4281Institute for Personalized Medicine, Pennsylvania State University College of Medicine, Hershey, PA 17033 USA ,grid.29857.310000 0001 2097 4281Department of Public Health Sciences, Pennsylvania State University College of Medicine, Hershey, PA 17033 USA
| | - Renan Sauteraud
- grid.29857.310000 0001 2097 4281Institute for Personalized Medicine, Pennsylvania State University College of Medicine, Hershey, PA 17033 USA ,grid.29857.310000 0001 2097 4281Department of Public Health Sciences, Pennsylvania State University College of Medicine, Hershey, PA 17033 USA
| | - Fang Chen
- grid.29857.310000 0001 2097 4281Institute for Personalized Medicine, Pennsylvania State University College of Medicine, Hershey, PA 17033 USA ,grid.29857.310000 0001 2097 4281Department of Public Health Sciences, Pennsylvania State University College of Medicine, Hershey, PA 17033 USA
| | - Lina Yang
- grid.29857.310000 0001 2097 4281Institute for Personalized Medicine, Pennsylvania State University College of Medicine, Hershey, PA 17033 USA ,grid.29857.310000 0001 2097 4281Department of Public Health Sciences, Pennsylvania State University College of Medicine, Hershey, PA 17033 USA
| | - Lida Wang
- grid.29857.310000 0001 2097 4281Institute for Personalized Medicine, Pennsylvania State University College of Medicine, Hershey, PA 17033 USA ,grid.29857.310000 0001 2097 4281Department of Public Health Sciences, Pennsylvania State University College of Medicine, Hershey, PA 17033 USA
| | - Jordan Hughey
- grid.29857.310000 0001 2097 4281Bioinformatics and Genomics Graduate Program, Pennsylvania State University College of Medicine, Hershey, PA 17033 USA ,grid.29857.310000 0001 2097 4281Institute for Personalized Medicine, Pennsylvania State University College of Medicine, Hershey, PA 17033 USA
| | - Scott Eckert
- grid.29857.310000 0001 2097 4281Bioinformatics and Genomics Graduate Program, Pennsylvania State University College of Medicine, Hershey, PA 17033 USA ,grid.29857.310000 0001 2097 4281Institute for Personalized Medicine, Pennsylvania State University College of Medicine, Hershey, PA 17033 USA
| | - J. Dylan Weissenkampen
- grid.29857.310000 0001 2097 4281Institute for Personalized Medicine, Pennsylvania State University College of Medicine, Hershey, PA 17033 USA ,grid.29857.310000 0001 2097 4281Department of Public Health Sciences, Pennsylvania State University College of Medicine, Hershey, PA 17033 USA
| | - Ganesh Shenoy
- grid.29857.310000 0001 2097 4281Department of Neurosurgery, Pennsylvania State University College of Medicine, Hershey, PA 17033 USA
| | - Olivia Marx
- grid.29857.310000 0001 2097 4281Biomedical Science Program, Pennsylvania State University College of Medicine, Hershey, PA 17033 USA
| | - Laura Carrel
- grid.29857.310000 0001 2097 4281Department of Biochemistry and Molecular Biology, Pennsylvania State University College of Medicine, Hershey, PA 17033 USA
| | - Bibo Jiang
- grid.29857.310000 0001 2097 4281Department of Public Health Sciences, Pennsylvania State University College of Medicine, Hershey, PA 17033 USA
| | - Dajiang J. Liu
- grid.29857.310000 0001 2097 4281Bioinformatics and Genomics Graduate Program, Pennsylvania State University College of Medicine, Hershey, PA 17033 USA ,grid.29857.310000 0001 2097 4281Institute for Personalized Medicine, Pennsylvania State University College of Medicine, Hershey, PA 17033 USA ,grid.29857.310000 0001 2097 4281Department of Public Health Sciences, Pennsylvania State University College of Medicine, Hershey, PA 17033 USA
| |
Collapse
|
37
|
Miller BF, Huang F, Atta L, Sahoo A, Fan J. Reference-free cell type deconvolution of multi-cellular pixel-resolution spatially resolved transcriptomics data. Nat Commun 2022; 13:2339. [PMID: 35487922 PMCID: PMC9055051 DOI: 10.1038/s41467-022-30033-z] [Citation(s) in RCA: 37] [Impact Index Per Article: 18.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/16/2021] [Accepted: 04/12/2022] [Indexed: 12/12/2022] Open
Abstract
Recent technological advancements have enabled spatially resolved transcriptomic profiling but at multi-cellular pixel resolution, thereby hindering the identification of cell-type-specific spatial patterns and gene expression variation. To address this challenge, we develop STdeconvolve as a reference-free approach to deconvolve underlying cell types comprising such multi-cellular pixel resolution spatial transcriptomics (ST) datasets. Using simulated as well as real ST datasets from diverse spatial transcriptomics technologies comprising a variety of spatial resolutions such as Spatial Transcriptomics, 10X Visium, DBiT-seq, and Slide-seq, we show that STdeconvolve can effectively recover cell-type transcriptional profiles and their proportional representation within pixels without reliance on external single-cell transcriptomics references. STdeconvolve provides comparable performance to existing reference-based methods when suitable single-cell references are available, as well as potentially superior performance when suitable single-cell references are not available. STdeconvolve is available as an open-source R software package with the source code available at https://github.com/JEFworks-Lab/STdeconvolve. Identifying cell-type-specific spatial patterns in ST data is critical for understanding tissue organization but current methods rely on external references. Here the authors develop a reference-free method to effectively recover cell-type transcriptional profiles and proportions.
Collapse
Affiliation(s)
- Brendan F Miller
- Center for Computational Biology, Whiting School of Engineering, Johns Hopkins University, Baltimore, MD, 21211, United States.,Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD, 21218, United States
| | - Feiyang Huang
- Center for Computational Biology, Whiting School of Engineering, Johns Hopkins University, Baltimore, MD, 21211, United States.,Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD, 21218, United States
| | - Lyla Atta
- Center for Computational Biology, Whiting School of Engineering, Johns Hopkins University, Baltimore, MD, 21211, United States.,Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD, 21218, United States
| | - Arpan Sahoo
- Center for Computational Biology, Whiting School of Engineering, Johns Hopkins University, Baltimore, MD, 21211, United States.,Department of Computer Science, Johns Hopkins University, Baltimore, MD, 21218, United States
| | - Jean Fan
- Center for Computational Biology, Whiting School of Engineering, Johns Hopkins University, Baltimore, MD, 21211, United States. .,Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD, 21218, United States. .,Department of Computer Science, Johns Hopkins University, Baltimore, MD, 21218, United States.
| |
Collapse
|
38
|
Ruiz Tejada Segura ML, Abou Moussa E, Garabello E, Nakahara TS, Makhlouf M, Mathew LS, Wang L, Valle F, Huang SSY, Mainland JD, Caselle M, Osella M, Lorenz S, Reisert J, Logan DW, Malnic B, Scialdone A, Saraiva LR. A 3D transcriptomics atlas of the mouse nose sheds light on the anatomical logic of smell. Cell Rep 2022; 38:110547. [PMID: 35320714 PMCID: PMC8995392 DOI: 10.1016/j.celrep.2022.110547] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/13/2021] [Revised: 01/26/2022] [Accepted: 03/01/2022] [Indexed: 12/26/2022] Open
Abstract
The sense of smell helps us navigate the environment, but its molecular architecture and underlying logic remain understudied. The spatial location of odorant receptor genes (Olfrs) in the nose is thought to be independent of the structural diversity of the odorants they detect. Using spatial transcriptomics, we create a genome-wide 3D atlas of the mouse olfactory mucosa (OM). Topographic maps of genes differentially expressed in space reveal that both Olfrs and non-Olfrs are distributed in a continuous and overlapping fashion over at least five broad zones in the OM. The spatial locations of Olfrs correlate with the mucus solubility of the odorants they recognize, providing direct evidence for the chromatographic theory of olfaction. This resource resolves the molecular architecture of the mouse OM and will inform future studies on mechanisms underlying Olfr gene choice, axonal pathfinding, patterning of the nervous system, and basic logic for the peripheral representation of smell.
Collapse
Affiliation(s)
- Mayra L Ruiz Tejada Segura
- Institute of Epigenetics and Stem Cells, Helmholtz Zentrum München, Feodor-Lynen-Strasse 21, 81377 München, Germany; Institute of Functional Epigenetics, Helmholtz Zentrum München, Ingolstädter Landstraße 1, 85764 Neuherberg, Germany; Institute of Computational Biology, Helmholtz Zentrum München, Ingolstädter Landstraße 1, 85764 Neuherberg, Germany
| | | | - Elisa Garabello
- Institute of Epigenetics and Stem Cells, Helmholtz Zentrum München, Feodor-Lynen-Strasse 21, 81377 München, Germany; Physics Department, University of Turin and INFN, Via P. Giuria 1, 10125 Turin, Italy; Department of Civil and Environmental Engineering, Cornell University, Ithaca, NY 14853, USA
| | - Thiago S Nakahara
- Department of Biochemistry, University of São Paulo, São Paulo, Brazil
| | | | | | - Li Wang
- Sidra Medicine, P.O. Box 26999, Doha, Qatar
| | - Filippo Valle
- Physics Department, University of Turin and INFN, Via P. Giuria 1, 10125 Turin, Italy
| | | | - Joel D Mainland
- Monell Chemical Senses Center, 3500 Market Street, Philadelphia, PA 19104, USA; Department of Neuroscience, University of Pennsylvania, Philadelphia, PA 19104, USA
| | - Michele Caselle
- Physics Department, University of Turin and INFN, Via P. Giuria 1, 10125 Turin, Italy
| | - Matteo Osella
- Physics Department, University of Turin and INFN, Via P. Giuria 1, 10125 Turin, Italy
| | - Stephan Lorenz
- Sidra Medicine, P.O. Box 26999, Doha, Qatar; Wellcome Sanger Institute, Wellcome Genome Campus, Hinxton, Cambridge CB10 1SD, UK
| | - Johannes Reisert
- Monell Chemical Senses Center, 3500 Market Street, Philadelphia, PA 19104, USA
| | - Darren W Logan
- Wellcome Sanger Institute, Wellcome Genome Campus, Hinxton, Cambridge CB10 1SD, UK
| | - Bettina Malnic
- Department of Biochemistry, University of São Paulo, São Paulo, Brazil
| | - Antonio Scialdone
- Institute of Epigenetics and Stem Cells, Helmholtz Zentrum München, Feodor-Lynen-Strasse 21, 81377 München, Germany; Institute of Functional Epigenetics, Helmholtz Zentrum München, Ingolstädter Landstraße 1, 85764 Neuherberg, Germany; Institute of Computational Biology, Helmholtz Zentrum München, Ingolstädter Landstraße 1, 85764 Neuherberg, Germany.
| | - Luis R Saraiva
- Sidra Medicine, P.O. Box 26999, Doha, Qatar; Monell Chemical Senses Center, 3500 Market Street, Philadelphia, PA 19104, USA; College of Health and Life Sciences, Hamad Bin Khalifa University, P.O. Box 34110, Doha, Qatar.
| |
Collapse
|
39
|
Adossa NA, Rytkönen KT, Elo LL. Dirichlet process mixture models for single-cell RNA-seq clustering. Biol Open 2022; 11:274586. [PMID: 35237784 PMCID: PMC9002799 DOI: 10.1242/bio.059001] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/21/2021] [Accepted: 02/17/2022] [Indexed: 11/20/2022] Open
Abstract
Clustering of cells based on gene expression is one of the major steps in single-cell RNA-sequencing (scRNA-seq) data analysis. One key challenge in cluster analysis is the unknown number of clusters and, for this issue, there is still no comprehensive solution. To enhance the process of defining meaningful cluster resolution, we compare Bayesian latent Dirichlet allocation (LDA) method to its non-parametric counterpart, hierarchical Dirichlet process (HDP) in the context of clustering scRNA-seq data. A potential main advantage of HDP is that it does not require the number of clusters as an input parameter from the user. While LDA has been used in single-cell data analysis, it has not been compared in detail with HDP. Here, we compare the cell clustering performance of LDA and HDP using four scRNA-seq datasets (immune cells, kidney, pancreas and decidua/placenta), with a specific focus on cluster numbers. Using both intrinsic (DB-index) and extrinsic (ARI) cluster quality measures, we show that the performance of LDA and HDP is dataset dependent. We describe a case where HDP produced a more appropriate clustering compared to the best performer from a series of LDA clusterings with different numbers of clusters. However, we also observed cases where the best performing LDA cluster numbers appropriately capture the main biological features while HDP tended to inflate the number of clusters. Overall, our study highlights the importance of carefully assessing the number of clusters when analyzing scRNA-seq data. Summary: Dirichlet mixture models (LDA and HDP) are applied for clustering cells in scRNA-Seq data. Here we made a comprehensive comparison of LDA and HDP model-based clustering for scRNA-seq data.
Collapse
Affiliation(s)
- Nigatu A Adossa
- Turku Bioscience Centre, University of Turku and Åbo Akademi University, FI-20520, Turku, Finland
| | - Kalle T Rytkönen
- Turku Bioscience Centre, University of Turku and Åbo Akademi University, FI-20520, Turku, Finland.,Institute of Biomedicine, Research Centre for Integrative Physiology and Pharmacology, University of Turku, FI-20014, Finland
| | - Laura L Elo
- Turku Bioscience Centre, University of Turku and Åbo Akademi University, FI-20520, Turku, Finland.,Institute of Biomedicine, University of Turku, FI-20014, Finland
| |
Collapse
|
40
|
Housman G, Briscoe E, Gilad Y. Evolutionary insights into primate skeletal gene regulation using a comparative cell culture model. PLoS Genet 2022; 18:e1010073. [PMID: 35263340 PMCID: PMC8936463 DOI: 10.1371/journal.pgen.1010073] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/02/2021] [Revised: 03/21/2022] [Accepted: 02/02/2022] [Indexed: 01/10/2023] Open
Abstract
The evolution of complex skeletal traits in primates was likely influenced by both genetic and environmental factors. Because skeletal tissues are notoriously challenging to study using functional genomic approaches, they remain poorly characterized even in humans, let alone across multiple species. The challenges involved in obtaining functional genomic data from the skeleton, combined with the difficulty of obtaining such tissues from nonhuman apes, motivated us to consider an alternative in vitro system with which to comparatively study gene regulation in skeletal cell types. Specifically, we differentiated six human (Homo sapiens) and six chimpanzee (Pan troglodytes) induced pluripotent stem cell lines (iPSCs) into mesenchymal stem cells (MSCs) and subsequently into osteogenic cells (bone cells). We validated differentiation using standard methods and collected single-cell RNA sequencing data from over 100,000 cells across multiple samples and replicates at each stage of differentiation. While most genes that we examined display conserved patterns of expression across species, hundreds of genes are differentially expressed (DE) between humans and chimpanzees within and across stages of osteogenic differentiation. Some of these interspecific DE genes show functional enrichments relevant in skeletal tissue trait development. Moreover, topic modeling indicates that interspecific gene programs become more pronounced as cells mature. Overall, we propose that this in vitro model can be used to identify interspecific regulatory differences that may have contributed to skeletal trait differences between species. Primates display a range of skeletal morphologies and susceptibilities to skeletal diseases, but the molecular basis of these phenotypic differences is unclear. Studies of gene expression variation in primate skeletal tissues are extremely restricted due to the ethical and practical challenges associated with collecting samples. Nevertheless, the ability to study gene regulation in primate skeletal tissues is crucial for understanding how the primate skeleton has evolved. We therefore developed a comparative primate skeletal cell culture model that allows us to access a spectrum of human and chimpanzee cell types as they differentiate from stem cells into bone cells. While most gene expression patterns are conserved across species, we also identified hundreds of differentially expressed genes between humans and chimpanzees within and across stages of differentiation. We also classified cells by osteogenic stage and identified additional interspecific differentially expressed genes which may contribute to skeletal trait differences. We anticipate that this model will be extremely useful for exploring questions related to gene regulation variation in primate bone biology and development.
Collapse
Affiliation(s)
- Genevieve Housman
- Section of Genetic Medicine, Department of Medicine, University of Chicago, Chicago, Illinois, United States of America
- * E-mail:
| | - Emilie Briscoe
- Section of Genetic Medicine, Department of Medicine, University of Chicago, Chicago, Illinois, United States of America
| | - Yoav Gilad
- Section of Genetic Medicine, Department of Medicine, University of Chicago, Chicago, Illinois, United States of America
- Department of Human Genetics, University of Chicago, Chicago, Illinois, United States of America
| |
Collapse
|
41
|
Multiomics Topic Modeling for Breast Cancer Classification. Cancers (Basel) 2022; 14:cancers14051150. [PMID: 35267458 PMCID: PMC8909787 DOI: 10.3390/cancers14051150] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/11/2022] [Accepted: 02/18/2022] [Indexed: 12/04/2022] Open
Abstract
The integration of transcriptional data with other layers of information, such as the post-transcriptional regulation mediated by microRNAs, can be crucial to identify the driver genes and the subtypes of complex and heterogeneous diseases such as cancer. This paper presents an approach based on topic modeling to accomplish this integration task. More specifically, we show how an algorithm based on a hierarchical version of stochastic block modeling can be naturally extended to integrate any combination of 'omics data. We test this approach on breast cancer samples from the TCGA database, integrating data on messenger RNA, microRNAs, and copy number variations. We show that the inclusion of the microRNA layer significantly improves the accuracy of subtype classification. Moreover, some of the hidden structures or "topics" that the algorithm extracts actually correspond to genes and microRNAs involved in breast cancer development and are associated to the survival probability.
Collapse
|
42
|
Rhodes K, Barr KA, Popp JM, Strober BJ, Battle A, Gilad Y. Human embryoid bodies as a novel system for genomic studies of functionally diverse cell types. eLife 2022; 11:71361. [PMID: 35142607 PMCID: PMC8830892 DOI: 10.7554/elife.71361] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/17/2021] [Accepted: 12/18/2021] [Indexed: 12/23/2022] Open
Abstract
Practically all studies of gene expression in humans to date have been performed in a relatively small number of adult tissues. Gene regulation is highly dynamic and context-dependent. In order to better understand the connection between gene regulation and complex phenotypes, including disease, we need to be able to study gene expression in more cell types, tissues, and states that are relevant to human phenotypes. In particular, we need to characterize gene expression in early development cell types, as mutations that affect developmental processes may be of particular relevance to complex traits. To address this challenge, we propose to use embryoid bodies (EBs), which are organoids that contain a multitude of cell types in dynamic states. EBs provide a system in which one can study dynamic regulatory processes at an unprecedentedly high resolution. To explore the utility of EBs, we systematically explored cellular and gene expression heterogeneity in EBs from multiple individuals. We characterized the various cell types that arise from EBs, the extent to which they recapitulate gene expression in vivo, and the relative contribution of technical and biological factors to variability in gene expression, cell composition, and differentiation efficiency. Our results highlight the utility of EBs as a new model system for mapping dynamic inter-individual regulatory differences in a large variety of cell types. One major goal of human genetics is to understand how changes in the way genes are regulated affect human traits, including disease susceptibility. To date, most studies of gene regulation have been performed in adult tissues, such as liver or kidney tissue, that were collected at a single time point. Yet, gene regulation is highly dynamic and context-dependent, meaning that it is important to gather data from a greater variety of cell types at different stages of their development. Additionally, observing which genes switch on and off in response to external treatments can shed light on how genetic variation can drive errors in gene regulation and cause diseases. Stem cells can produce more cells like themselves or differentiate – acquire the characteristics – of many cell types. These cells have been used in the laboratory to research gene regulation. Unfortunately, these studies often fail to capture the complex spatial and temporal dynamics of stem cell differentiation; in particular, these studies are unable to observe gene regulation in the transient cell types that appear early in embryonic development. To overcome these limitations, scientists developed systems such as embryoid bodies: three-dimensional aggregates of stem cells that, when grown under certain conditions, spontaneously develop into a variety of cell types. Rhodes, Barr et al. wanted to assess the utility of embryoid bodies as a model to study how genes are dynamically regulated in different cell types, by different individuals who have distinct genetic makeups. To do this, they grew embryoid bodies made from human stem cells from different individuals to examine which genes switched on and off as the stem cells that formed the embryoid bodies differentiated into different types of cells. The results showed that it was possible to grow embryoid bodies derived from genetically distinct individuals that consistently produce diverse cell types, similar to those found during human fetal development. Rhodes, Barr et al.’s findings suggest that embryoid bodies are a useful model to study gene regulation across individuals with different genetic backgrounds. This could accelerate research into how genetics are associated with disease by capturing gene regulatory dynamics at an unprecedentedly high spatial and temporal resolution. Additionally, embryoid bodies could be used to explore how exposure to different environmental factors during early development affect disease-related outcomes in adulthood in different individuals.
Collapse
Affiliation(s)
- Katherine Rhodes
- Department of Medicine, University of Chicago, Chicago, United States
| | - Kenneth A Barr
- Department of Medicine, University of Chicago, Chicago, United States
| | - Joshua M Popp
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, United States
| | - Benjamin J Strober
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, United States
| | - Alexis Battle
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, United States.,Department of Computer Science, Johns Hopkins University, Baltimore, United States.,Department of Genetic Medicine, Johns Hopkins University, Baltimore, Maryland, United States
| | - Yoav Gilad
- Department of Medicine, University of Chicago, Chicago, United States
| |
Collapse
|
43
|
Akagbosu B, Tayyebi Z, Shibu G, Paucar Iza YA, Deep D, Parisotto YF, Fisher L, Pasolli HA, Thevin V, Elmentaite R, Knott M, Hemmers S, Jahn L, Friedrich C, Verter J, Wang ZM, van den Brink M, Gasteiger G, Grünewald TGP, Marie JC, Leslie C, Rudensky AY, Brown CC. Novel antigen-presenting cell imparts T reg-dependent tolerance to gut microbiota. Nature 2022; 610:752-760. [PMID: 36070798 PMCID: PMC9605865 DOI: 10.1038/s41586-022-05309-5] [Citation(s) in RCA: 81] [Impact Index Per Article: 40.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/09/2021] [Accepted: 09/01/2022] [Indexed: 01/21/2023]
Abstract
Establishing and maintaining tolerance to self-antigens or innocuous foreign antigens is vital for the preservation of organismal health. Within the thymus, medullary thymic epithelial cells (mTECs) expressing autoimmune regulator (AIRE) have a critical role in self-tolerance through deletion of autoreactive T cells and promotion of thymic regulatory T (Treg) cell development1-4. Within weeks of birth, a separate wave of Treg cell differentiation occurs in the periphery upon exposure to antigens derived from the diet and commensal microbiota5-8, yet the cell types responsible for the generation of peripheral Treg (pTreg) cells have not been identified. Here we describe the identification of a class of RORγt+ antigen-presenting cells called Thetis cells, with transcriptional features of both mTECs and dendritic cells, comprising four major sub-groups (TC I-TC IV). We uncover a developmental wave of Thetis cells within intestinal lymph nodes during a critical window in early life, coinciding with the wave of pTreg cell differentiation. Whereas TC I and TC III expressed the signature mTEC nuclear factor AIRE, TC IV lacked AIRE expression and was enriched for molecules required for pTreg generation, including the TGF-β-activating integrin αvβ8. Loss of either major histocompatibility complex class II (MHCII) or ITGB8 by Thetis cells led to a profound impairment in intestinal pTreg differentiation, with ensuing colitis. By contrast, MHCII expression by RORγt+ group 3 innate lymphoid cells (ILC3) and classical dendritic cells was neither sufficient nor required for pTreg generation, further implicating TC IV as the tolerogenic RORγt+ antigen-presenting cell with an essential function in early life. Our studies reveal parallel pathways for the establishment of tolerance to self and foreign antigens in the thymus and periphery, respectively, marked by the involvement of shared cellular and transcriptional programmes.
Collapse
Affiliation(s)
- Blossom Akagbosu
- grid.51462.340000 0001 2171 9952Immuno-Oncology, Human Oncology and Pathogenesis Program, Memorial Sloan Kettering Cancer Center, New York, USA
| | - Zakieh Tayyebi
- grid.51462.340000 0001 2171 9952Computational and Systems Biology Program, Memorial Sloan Kettering Cancer Center, New York, NY USA ,grid.5386.8000000041936877XTri-Institutional Program in Computational Biology and Medicine, Weill Cornell Graduate School, New York, NY USA
| | - Gayathri Shibu
- grid.51462.340000 0001 2171 9952Immuno-Oncology, Human Oncology and Pathogenesis Program, Memorial Sloan Kettering Cancer Center, New York, USA ,grid.5386.8000000041936877XImmunology and Microbial Pathogenesis Program, Weill Cornell Medicine Graduate School of Medical Sciences, New York, NY USA
| | - Yoselin A. Paucar Iza
- grid.51462.340000 0001 2171 9952Immuno-Oncology, Human Oncology and Pathogenesis Program, Memorial Sloan Kettering Cancer Center, New York, USA ,grid.5386.8000000041936877XImmunology and Microbial Pathogenesis Program, Weill Cornell Medicine Graduate School of Medical Sciences, New York, NY USA
| | - Deeksha Deep
- grid.5386.8000000041936877XImmunology and Microbial Pathogenesis Program, Weill Cornell Medicine Graduate School of Medical Sciences, New York, NY USA ,grid.51462.340000 0001 2171 9952Howard Hughes Medical Institute and Immunology Program, Sloan Kettering Institute and Ludwig Center at Memorial Sloan Kettering Cancer Center, New York, NY USA ,grid.5386.8000000041936877XTri-Institutional MD-PhD Program, Weill Cornell Medicine, The Rockefeller University and Memorial Sloan Kettering Cancer Center, New York, NY USA
| | - Yollanda Franco Parisotto
- grid.51462.340000 0001 2171 9952Immuno-Oncology, Human Oncology and Pathogenesis Program, Memorial Sloan Kettering Cancer Center, New York, USA
| | - Logan Fisher
- grid.51462.340000 0001 2171 9952Immuno-Oncology, Human Oncology and Pathogenesis Program, Memorial Sloan Kettering Cancer Center, New York, USA ,grid.5386.8000000041936877XImmunology and Microbial Pathogenesis Program, Weill Cornell Medicine Graduate School of Medical Sciences, New York, NY USA
| | - H. Amalia Pasolli
- grid.134907.80000 0001 2166 1519Electron Microscopy Resource Center, The Rockefeller University, New York, NY USA
| | - Valentin Thevin
- grid.462282.80000 0004 0384 0005Tumor Escape Resistance Immunity Department, CRCL, INSERM U1052, CNRS 5286, Centre Léon Bérard, Université de Lyon, Lyon, France ,Equipe Labellisée Ligue Nationale contre le Cancer, Lyon, France
| | - Rasa Elmentaite
- grid.10306.340000 0004 0606 5382Wellcome Sanger Institute, Wellcome Genome Campus, Hinxton UK
| | - Maximilian Knott
- grid.5252.00000 0004 1936 973XInstitute of PathologyFaculty of Medicine, LMU Munich, Munich, Germany
| | - Saskia Hemmers
- grid.5386.8000000041936877XImmunology and Microbial Pathogenesis Program, Weill Cornell Medicine Graduate School of Medical Sciences, New York, NY USA ,grid.51462.340000 0001 2171 9952Howard Hughes Medical Institute and Immunology Program, Sloan Kettering Institute and Ludwig Center at Memorial Sloan Kettering Cancer Center, New York, NY USA ,grid.26009.3d0000 0004 1936 7961Present Address: Department of Immunology, Duke University, Durham, NC USA
| | - Lorenz Jahn
- grid.51462.340000 0001 2171 9952Department of Immunology, Memorial Sloan Kettering Cancer Center, New York, NY USA
| | - Christin Friedrich
- grid.8379.50000 0001 1958 8658Würzburg Institute of Systems Immunology, Max Planck Research Group, Julius-Maximilians-Universität Würzburg, Würzburg, Germany
| | - Jacob Verter
- grid.51462.340000 0001 2171 9952Howard Hughes Medical Institute and Immunology Program, Sloan Kettering Institute and Ludwig Center at Memorial Sloan Kettering Cancer Center, New York, NY USA
| | - Zhong-Min Wang
- grid.51462.340000 0001 2171 9952Howard Hughes Medical Institute and Immunology Program, Sloan Kettering Institute and Ludwig Center at Memorial Sloan Kettering Cancer Center, New York, NY USA
| | - Marcel van den Brink
- grid.5386.8000000041936877XImmunology and Microbial Pathogenesis Program, Weill Cornell Medicine Graduate School of Medical Sciences, New York, NY USA ,grid.51462.340000 0001 2171 9952Department of Immunology, Memorial Sloan Kettering Cancer Center, New York, NY USA ,grid.51462.340000 0001 2171 9952Department of Medicine, Memorial Sloan Kettering Cancer Center, New York, NY USA
| | - Georg Gasteiger
- grid.8379.50000 0001 1958 8658Würzburg Institute of Systems Immunology, Max Planck Research Group, Julius-Maximilians-Universität Würzburg, Würzburg, Germany
| | - Thomas G. P. Grünewald
- grid.510964.fHopp—Children’s Cancer Center Heidelberg (KiTZ), Heidelberg, Germany ,grid.7497.d0000 0004 0492 0584Division of Translational Pediatric Sarcoma Research, German Cancer Research Center (DKFZ), German Cancer Consortium (DKTK), Heidelberg, Germany ,grid.5253.10000 0001 0328 4908Institute of Pathology, Heidelberg University Hospital, Heidelberg, Germany
| | - Julien C. Marie
- grid.462282.80000 0004 0384 0005Tumor Escape Resistance Immunity Department, CRCL, INSERM U1052, CNRS 5286, Centre Léon Bérard, Université de Lyon, Lyon, France ,Equipe Labellisée Ligue Nationale contre le Cancer, Lyon, France
| | - Christina Leslie
- grid.51462.340000 0001 2171 9952Computational and Systems Biology Program, Memorial Sloan Kettering Cancer Center, New York, NY USA
| | - Alexander Y. Rudensky
- grid.5386.8000000041936877XImmunology and Microbial Pathogenesis Program, Weill Cornell Medicine Graduate School of Medical Sciences, New York, NY USA ,grid.51462.340000 0001 2171 9952Howard Hughes Medical Institute and Immunology Program, Sloan Kettering Institute and Ludwig Center at Memorial Sloan Kettering Cancer Center, New York, NY USA
| | - Chrysothemis C. Brown
- grid.51462.340000 0001 2171 9952Immuno-Oncology, Human Oncology and Pathogenesis Program, Memorial Sloan Kettering Cancer Center, New York, USA ,grid.5386.8000000041936877XImmunology and Microbial Pathogenesis Program, Weill Cornell Medicine Graduate School of Medical Sciences, New York, NY USA ,grid.51462.340000 0001 2171 9952Department of Pediatrics, Memorial Sloan Kettering Cancer Center, New York, NY USA ,grid.51462.340000 0001 2171 9952Parker Institute for Cancer Immunotherapy, Memorial Sloan Kettering Cancer Center, New York, NY USA
| |
Collapse
|
44
|
G-protein-coupled receptor P2Y10 facilitates chemokine-induced CD4 T cell migration through autocrine/paracrine mediators. Nat Commun 2021; 12:6798. [PMID: 34815397 PMCID: PMC8611058 DOI: 10.1038/s41467-021-26882-9] [Citation(s) in RCA: 19] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/18/2020] [Accepted: 10/25/2021] [Indexed: 12/12/2022] Open
Abstract
G-protein-coupled receptors (GPCRs), especially chemokine receptors, play a central role in the regulation of T cell migration. Various GPCRs are upregulated in activated CD4 T cells, including P2Y10, a putative lysophospholipid receptor that is officially still considered an orphan GPCR, i.e., a receptor with unknown endogenous ligand. Here we show that in mice lacking P2Y10 in the CD4 T cell compartment, the severity of experimental autoimmune encephalomyelitis and cutaneous contact hypersensitivity is reduced. P2Y10-deficient CD4 T cells show normal activation, proliferation and differentiation, but reduced chemokine-induced migration, polarization, and RhoA activation upon in vitro stimulation. Mechanistically, CD4 T cells release the putative P2Y10 ligands lysophosphatidylserine and ATP upon chemokine exposure, and these mediators induce P2Y10-dependent RhoA activation in an autocrine/paracrine fashion. ATP degradation impairs RhoA activation and migration in control CD4 T cells, but not in P2Y10-deficient CD4 T cells. Importantly, the P2Y10 pathway appears to be conserved in human T cells. Taken together, P2Y10 mediates RhoA activation in CD4 T cells in response to auto-/paracrine-acting mediators such as LysoPS and ATP, thereby facilitating chemokine-induced migration and, consecutively, T cell-mediated diseases.
Collapse
|
45
|
Schenkel JM, Herbst RH, Canner D, Li A, Hillman M, Shanahan SL, Gibbons G, Smith OC, Kim JY, Westcott P, Hwang WL, Freed-Pastor WA, Eng G, Cuoco MS, Rogers P, Park JK, Burger ML, Rozenblatt-Rosen O, Cong L, Pauken KE, Regev A, Jacks T. Conventional type I dendritic cells maintain a reservoir of proliferative tumor-antigen specific TCF-1 + CD8 + T cells in tumor-draining lymph nodes. Immunity 2021; 54:2338-2353.e6. [PMID: 34534439 PMCID: PMC8604155 DOI: 10.1016/j.immuni.2021.08.026] [Citation(s) in RCA: 113] [Impact Index Per Article: 37.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/28/2020] [Revised: 05/22/2021] [Accepted: 08/20/2021] [Indexed: 12/28/2022]
Abstract
In tumors, a subset of CD8+ T cells expressing the transcription factor TCF-1 drives the response to immune checkpoint blockade. We examined the mechanisms that maintain these cells in an autochthonous model of lung adenocarcinoma. Longitudinal sampling and single-cell sequencing of tumor-antigen specific TCF-1+ CD8+ T cells revealed that while intratumoral TCF-1+ CD8+ T cells acquired dysfunctional features and decreased in number as tumors progressed, TCF-1+ CD8+ T cell frequency in the tumor draining LN (dLN) remained stable. Two discrete intratumoral TCF-1+ CD8+ T cell subsets developed over time-a proliferative SlamF6+ subset and a non-cycling SlamF6- subset. Blocking dLN egress decreased the frequency of intratumoral SlamF6+ TCF-1+ CD8+ T cells. Conventional type I dendritic cell (cDC1) in dLN decreased in number with tumor progression, and Flt3L+anti-CD40 treatment recovered SlamF6+ T cell frequencies and decreased tumor burden. Thus, cDC1s in tumor dLN maintain a reservoir of TCF-1+ CD8+ T cells and their decrease contributes to failed anti-tumor immunity.
Collapse
Affiliation(s)
- Jason M Schenkel
- David H. Koch Institute for Integrative Cancer Research, Massachusetts Institute of Technology, 500 Main Street, Cambridge, MA 02139, USA; Department of Pathology, Brigham and Women's Hospital, Boston, MA 02115, USA; Harvard Medical School, 25 Shattuck Street, Boston, MA 02115, USA
| | - Rebecca H Herbst
- Harvard Medical School, 25 Shattuck Street, Boston, MA 02115, USA; Broad Institute of MIT and Harvard, 415 Main Street, Cambridge, MA 02142, USA
| | - David Canner
- David H. Koch Institute for Integrative Cancer Research, Massachusetts Institute of Technology, 500 Main Street, Cambridge, MA 02139, USA; Department of Biology, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, MA 02139, USA
| | - Amy Li
- David H. Koch Institute for Integrative Cancer Research, Massachusetts Institute of Technology, 500 Main Street, Cambridge, MA 02139, USA; Harvard Medical School, 25 Shattuck Street, Boston, MA 02115, USA; Department of Biology, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, MA 02139, USA
| | - Michelle Hillman
- David H. Koch Institute for Integrative Cancer Research, Massachusetts Institute of Technology, 500 Main Street, Cambridge, MA 02139, USA
| | - Sean-Luc Shanahan
- David H. Koch Institute for Integrative Cancer Research, Massachusetts Institute of Technology, 500 Main Street, Cambridge, MA 02139, USA
| | - Grace Gibbons
- David H. Koch Institute for Integrative Cancer Research, Massachusetts Institute of Technology, 500 Main Street, Cambridge, MA 02139, USA
| | - Olivia C Smith
- David H. Koch Institute for Integrative Cancer Research, Massachusetts Institute of Technology, 500 Main Street, Cambridge, MA 02139, USA
| | - Jonathan Y Kim
- David H. Koch Institute for Integrative Cancer Research, Massachusetts Institute of Technology, 500 Main Street, Cambridge, MA 02139, USA
| | - Peter Westcott
- David H. Koch Institute for Integrative Cancer Research, Massachusetts Institute of Technology, 500 Main Street, Cambridge, MA 02139, USA
| | - William L Hwang
- David H. Koch Institute for Integrative Cancer Research, Massachusetts Institute of Technology, 500 Main Street, Cambridge, MA 02139, USA; Broad Institute of MIT and Harvard, 415 Main Street, Cambridge, MA 02142, USA; Department of Radiation Oncology, Massachusetts General Hospital, Boston MA 02114
| | - William A Freed-Pastor
- David H. Koch Institute for Integrative Cancer Research, Massachusetts Institute of Technology, 500 Main Street, Cambridge, MA 02139, USA; Department of Medical Oncology, Dana-Farber Cancer Institute, Boston, MA, USA
| | - George Eng
- David H. Koch Institute for Integrative Cancer Research, Massachusetts Institute of Technology, 500 Main Street, Cambridge, MA 02139, USA; Department of Pathology, Massachusetts General Hospital, Boston, MA 02114, USA
| | - Michael S Cuoco
- Broad Institute of MIT and Harvard, 415 Main Street, Cambridge, MA 02142, USA
| | - Patricia Rogers
- Broad Institute of MIT and Harvard, 415 Main Street, Cambridge, MA 02142, USA
| | - Jin K Park
- David H. Koch Institute for Integrative Cancer Research, Massachusetts Institute of Technology, 500 Main Street, Cambridge, MA 02139, USA; Harvard Medical School, 25 Shattuck Street, Boston, MA 02115, USA
| | - Megan L Burger
- David H. Koch Institute for Integrative Cancer Research, Massachusetts Institute of Technology, 500 Main Street, Cambridge, MA 02139, USA
| | | | - Le Cong
- Departments of Pathology, Stanford University, Stanford, CA 94305, USA
| | - Kristen E Pauken
- Department of Immunology, Blavatnik Institute, Harvard Medical School, Boston, MA 02115, USA; Evergrande Center for Immunologic Diseases, Harvard Medical School and Brigham and Women's Hospital, Boston, MA 02115, USA
| | - Aviv Regev
- David H. Koch Institute for Integrative Cancer Research, Massachusetts Institute of Technology, 500 Main Street, Cambridge, MA 02139, USA; Broad Institute of MIT and Harvard, 415 Main Street, Cambridge, MA 02142, USA; Department of Biology, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, MA 02139, USA.
| | - Tyler Jacks
- David H. Koch Institute for Integrative Cancer Research, Massachusetts Institute of Technology, 500 Main Street, Cambridge, MA 02139, USA; Department of Biology, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, MA 02139, USA.
| |
Collapse
|
46
|
Vaginal microbiome topic modeling of laboring Ugandan women with and without fever. NPJ Biofilms Microbiomes 2021; 7:75. [PMID: 34508087 PMCID: PMC8433417 DOI: 10.1038/s41522-021-00244-1] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/10/2021] [Accepted: 08/13/2021] [Indexed: 12/12/2022] Open
Abstract
The composition of the maternal vaginal microbiome influences the duration of pregnancy, onset of labor, and even neonatal outcomes. Maternal microbiome research in sub-Saharan Africa has focused on non-pregnant and postpartum composition of the vaginal microbiome. Here we aimed to illustrate the relationship between the vaginal microbiome of 99 laboring Ugandan women and intrapartum fever using routine microbiology and 16S ribosomal RNA gene sequencing from two hypervariable regions (V1–V2 and V3–V4). To describe the vaginal microbes associated with vaginal microbial communities, we pursued two approaches: hierarchical clustering methods and a novel Grades of Membership (GoM) modeling approach for vaginal microbiome characterization. Leveraging GoM models, we created a basis composed of a preassigned number of microbial topics whose linear combination optimally represents each patient yielding more comprehensive associations and characterization between maternal clinical features and the microbial communities. Using a random forest model, we showed that by including microbial topic models we improved upon clinical variables to predict maternal fever. Overall, we found a higher prevalence of Granulicatella, Streptococcus, Fusobacterium, Anaerococcus, Sneathia, Clostridium, Gemella, Mobiluncus, and Veillonella genera in febrile mothers, and higher prevalence of Lactobacillus genera (in particular L. crispatus and L. jensenii), Acinobacter, Aerococcus, and Prevotella species in afebrile mothers. By including clinical variables with microbial topics in this model, we observed young maternal age, fever reported earlier in the pregnancy, longer labor duration, and microbial communities with reduced Lactobacillus diversity were associated with intrapartum fever. These results better defined relationships between the presence or absence of intrapartum fever, demographics, peripartum course, and vaginal microbial topics, and expanded our understanding of the impact of the microbiome on maternal and potentially neonatal outcome risk.
Collapse
|
47
|
Li L, Xiong F, Wang Y, Zhang S, Gong Z, Li X, He Y, Shi L, Wang F, Liao Q, Xiang B, Zhou M, Li X, Li Y, Li G, Zeng Z, Xiong W, Guo C. What are the applications of single-cell RNA sequencing in cancer research: a systematic review. JOURNAL OF EXPERIMENTAL & CLINICAL CANCER RESEARCH : CR 2021; 40:163. [PMID: 33975628 PMCID: PMC8111731 DOI: 10.1186/s13046-021-01955-1] [Citation(s) in RCA: 28] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Subscribe] [Scholar Register] [Received: 01/26/2021] [Accepted: 04/20/2021] [Indexed: 12/18/2022]
Abstract
Single-cell RNA sequencing (scRNA-seq) is a tool for studying gene expression at the single-cell level that has been widely used due to its unprecedented high resolution. In the present review, we outline the preparation process and sequencing platforms for the scRNA-seq analysis of solid tumor specimens and discuss the main steps and methods used during data analysis, including quality control, batch-effect correction, normalization, cell cycle phase assignment, clustering, cell trajectory and pseudo-time reconstruction, differential expression analysis and gene set enrichment analysis, as well as gene regulatory network inference. Traditional bulk RNA sequencing does not address the heterogeneity within and between tumors, and since the development of the first scRNA-seq technique, this approach has been widely used in cancer research to better understand cancer cell biology and pathogenetic mechanisms. ScRNA-seq has been of great significance for the development of targeted therapy and immunotherapy. In the second part of this review, we focus on the application of scRNA-seq in solid tumors, and summarize the findings and achievements in tumor research afforded by its use. ScRNA-seq holds promise for improving our understanding of the molecular characteristics of cancer, and potentially contributing to improved diagnosis, prognosis, and therapeutics.
Collapse
Affiliation(s)
- Lvyuan Li
- NHC Key Laboratory of Carcinogenesis and Hunan Key Laboratory of Cancer Metabolism, Hunan Cancer Hospital and the Affiliated Cancer Hospital of Xiangya School of Medicine, Central South University, Changsha, China.,Key Laboratory of Carcinogenesis and Cancer Invasion of the Chinese Ministry of Education, Cancer Research Institute, Central South University, Changsha, China
| | - Fang Xiong
- Department of Stomatology, Xiangya Hospital, Central South University, Changsha, China
| | - Yumin Wang
- Key Laboratory of Carcinogenesis and Cancer Invasion of the Chinese Ministry of Education, Cancer Research Institute, Central South University, Changsha, China.,Department of Stomatology, Xiangya Hospital, Central South University, Changsha, China
| | - Shanshan Zhang
- Department of Stomatology, Xiangya Hospital, Central South University, Changsha, China
| | - Zhaojian Gong
- Department of Oral and Maxillofacial Surgery, The Second Xiangya Hospital, Central South University, Changsha, China
| | - Xiayu Li
- Hunan Key Laboratory of Nonresolving Inflammation and Cancer, Disease Genome Research Center, The Third Xiangya Hospital, Central South University, Changsha, China
| | - Yi He
- NHC Key Laboratory of Carcinogenesis and Hunan Key Laboratory of Cancer Metabolism, Hunan Cancer Hospital and the Affiliated Cancer Hospital of Xiangya School of Medicine, Central South University, Changsha, China
| | - Lei Shi
- Department of Oral and Maxillofacial Surgery, The Second Xiangya Hospital, Central South University, Changsha, China
| | - Fuyan Wang
- Key Laboratory of Carcinogenesis and Cancer Invasion of the Chinese Ministry of Education, Cancer Research Institute, Central South University, Changsha, China
| | - Qianjin Liao
- NHC Key Laboratory of Carcinogenesis and Hunan Key Laboratory of Cancer Metabolism, Hunan Cancer Hospital and the Affiliated Cancer Hospital of Xiangya School of Medicine, Central South University, Changsha, China
| | - Bo Xiang
- NHC Key Laboratory of Carcinogenesis and Hunan Key Laboratory of Cancer Metabolism, Hunan Cancer Hospital and the Affiliated Cancer Hospital of Xiangya School of Medicine, Central South University, Changsha, China.,Key Laboratory of Carcinogenesis and Cancer Invasion of the Chinese Ministry of Education, Cancer Research Institute, Central South University, Changsha, China
| | - Ming Zhou
- NHC Key Laboratory of Carcinogenesis and Hunan Key Laboratory of Cancer Metabolism, Hunan Cancer Hospital and the Affiliated Cancer Hospital of Xiangya School of Medicine, Central South University, Changsha, China.,Key Laboratory of Carcinogenesis and Cancer Invasion of the Chinese Ministry of Education, Cancer Research Institute, Central South University, Changsha, China
| | - Xiaoling Li
- NHC Key Laboratory of Carcinogenesis and Hunan Key Laboratory of Cancer Metabolism, Hunan Cancer Hospital and the Affiliated Cancer Hospital of Xiangya School of Medicine, Central South University, Changsha, China.,Key Laboratory of Carcinogenesis and Cancer Invasion of the Chinese Ministry of Education, Cancer Research Institute, Central South University, Changsha, China
| | - Yong Li
- Department of Medicine, Dan L Duncan Comprehensive Cancer Center, Baylor College of Medicine, Houston, TX, USA
| | - Guiyuan Li
- NHC Key Laboratory of Carcinogenesis and Hunan Key Laboratory of Cancer Metabolism, Hunan Cancer Hospital and the Affiliated Cancer Hospital of Xiangya School of Medicine, Central South University, Changsha, China.,Key Laboratory of Carcinogenesis and Cancer Invasion of the Chinese Ministry of Education, Cancer Research Institute, Central South University, Changsha, China
| | - Zhaoyang Zeng
- NHC Key Laboratory of Carcinogenesis and Hunan Key Laboratory of Cancer Metabolism, Hunan Cancer Hospital and the Affiliated Cancer Hospital of Xiangya School of Medicine, Central South University, Changsha, China.,Key Laboratory of Carcinogenesis and Cancer Invasion of the Chinese Ministry of Education, Cancer Research Institute, Central South University, Changsha, China
| | - Wei Xiong
- NHC Key Laboratory of Carcinogenesis and Hunan Key Laboratory of Cancer Metabolism, Hunan Cancer Hospital and the Affiliated Cancer Hospital of Xiangya School of Medicine, Central South University, Changsha, China. .,Key Laboratory of Carcinogenesis and Cancer Invasion of the Chinese Ministry of Education, Cancer Research Institute, Central South University, Changsha, China.
| | - Can Guo
- NHC Key Laboratory of Carcinogenesis and Hunan Key Laboratory of Cancer Metabolism, Hunan Cancer Hospital and the Affiliated Cancer Hospital of Xiangya School of Medicine, Central South University, Changsha, China. .,Key Laboratory of Carcinogenesis and Cancer Invasion of the Chinese Ministry of Education, Cancer Research Institute, Central South University, Changsha, China.
| |
Collapse
|
48
|
Dani N, Herbst RH, McCabe C, Green GS, Kaiser K, Head JP, Cui J, Shipley FB, Jang A, Dionne D, Nguyen L, Rodman C, Riesenfeld SJ, Prochazka J, Prochazkova M, Sedlacek R, Zhang F, Bryja V, Rozenblatt-Rosen O, Habib N, Regev A, Lehtinen MK. A cellular and spatial map of the choroid plexus across brain ventricles and ages. Cell 2021; 184:3056-3074.e21. [PMID: 33932339 DOI: 10.1016/j.cell.2021.04.003] [Citation(s) in RCA: 142] [Impact Index Per Article: 47.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/15/2019] [Revised: 12/15/2020] [Accepted: 04/02/2021] [Indexed: 02/06/2023]
Abstract
The choroid plexus (ChP) in each brain ventricle produces cerebrospinal fluid (CSF) and forms the blood-CSF barrier. Here, we construct a single-cell and spatial atlas of each ChP in the developing, adult, and aged mouse brain. We delineate diverse cell types, subtypes, cell states, and expression programs in epithelial and mesenchymal cells across ages and ventricles. In the developing ChP, we predict a common progenitor pool for epithelial and neuronal cells, validated by lineage tracing. Epithelial and fibroblast cells show regionalized expression by ventricle, starting at embryonic stages and persisting with age, with a dramatic transcriptional shift with maturation, and a smaller shift in each aged cell type. With aging, epithelial cells upregulate host-defense programs, and resident macrophages upregulate interleukin-1β (IL-1β) signaling genes. Our atlas reveals cellular diversity, architecture and signaling across ventricles during development, maturation, and aging of the ChP-brain barrier.
Collapse
Affiliation(s)
- Neil Dani
- Department of Pathology, Boston Children's Hospital, Boston, MA 02115, USA
| | - Rebecca H Herbst
- Klarman Cell Observatory, Broad Institute of MIT and Harvard, Cambridge, MA 02142, USA; Department of Systems Biology, Harvard Medical School, Boston, MA 02115, USA
| | - Cristin McCabe
- Klarman Cell Observatory, Broad Institute of MIT and Harvard, Cambridge, MA 02142, USA
| | - Gilad S Green
- Edmond and Lily Safra Center for Brain Sciences, Hebrew University of Jerusalem, Jerusalem 9190401, Israel
| | - Karol Kaiser
- Department of Experimental Biology, Faculty of Science, Masaryk University, Brno 611 37, Czech Republic
| | - Joshua P Head
- Department of Pathology, Boston Children's Hospital, Boston, MA 02115, USA
| | - Jin Cui
- Department of Pathology, Boston Children's Hospital, Boston, MA 02115, USA
| | - Frederick B Shipley
- Department of Pathology, Boston Children's Hospital, Boston, MA 02115, USA; Graduate Program in Biophysics, Harvard University, Cambridge, MA 02115, USA
| | - Ahram Jang
- Department of Pathology, Boston Children's Hospital, Boston, MA 02115, USA
| | - Danielle Dionne
- Klarman Cell Observatory, Broad Institute of MIT and Harvard, Cambridge, MA 02142, USA
| | - Lan Nguyen
- Klarman Cell Observatory, Broad Institute of MIT and Harvard, Cambridge, MA 02142, USA
| | - Christopher Rodman
- Klarman Cell Observatory, Broad Institute of MIT and Harvard, Cambridge, MA 02142, USA
| | - Samantha J Riesenfeld
- Klarman Cell Observatory, Broad Institute of MIT and Harvard, Cambridge, MA 02142, USA
| | - Jan Prochazka
- Czech Centre for Phenogenomics and Laboratory of Transgenic Models of Diseases, Institute of Molecular Genetics of the CAS, Prague 142 20, Czech Republic
| | - Michaela Prochazkova
- Czech Centre for Phenogenomics and Laboratory of Transgenic Models of Diseases, Institute of Molecular Genetics of the CAS, Prague 142 20, Czech Republic
| | - Radislav Sedlacek
- Czech Centre for Phenogenomics and Laboratory of Transgenic Models of Diseases, Institute of Molecular Genetics of the CAS, Prague 142 20, Czech Republic
| | - Feng Zhang
- Broad Institute of MIT and Harvard, Cambridge, MA 02142, USA; Howard Hughes Medical Institute, Chevy Chase, MD 20815, USA; McGovern Institute for Brain Research, Massachusetts Institute of Technology, Cambridge, MA 02139, USA; Department of Brain and Cognitive Sciences and Department of Biological Engineering, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
| | - Vitezslav Bryja
- Department of Experimental Biology, Faculty of Science, Masaryk University, Brno 611 37, Czech Republic
| | - Orit Rozenblatt-Rosen
- Klarman Cell Observatory, Broad Institute of MIT and Harvard, Cambridge, MA 02142, USA
| | - Naomi Habib
- Broad Institute of MIT and Harvard, Cambridge, MA 02142, USA; Edmond and Lily Safra Center for Brain Sciences, Hebrew University of Jerusalem, Jerusalem 9190401, Israel.
| | - Aviv Regev
- Klarman Cell Observatory, Broad Institute of MIT and Harvard, Cambridge, MA 02142, USA; Broad Institute of MIT and Harvard, Cambridge, MA 02142, USA; Howard Hughes Medical Institute, Chevy Chase, MD 20815, USA; Koch Institute of Integrative Cancer Research, Department of Biology, Massachusetts Institute of Technology, Cambridge, MA 02142, USA.
| | - Maria K Lehtinen
- Department of Pathology, Boston Children's Hospital, Boston, MA 02115, USA; Graduate Program in Biophysics, Harvard University, Cambridge, MA 02115, USA.
| |
Collapse
|
49
|
Del Giudice M, Peirone S, Perrone S, Priante F, Varese F, Tirtei E, Fagioli F, Cereda M. Artificial Intelligence in Bulk and Single-Cell RNA-Sequencing Data to Foster Precision Oncology. Int J Mol Sci 2021; 22:ijms22094563. [PMID: 33925407 PMCID: PMC8123853 DOI: 10.3390/ijms22094563] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/20/2021] [Revised: 04/21/2021] [Accepted: 04/23/2021] [Indexed: 02/01/2023] Open
Abstract
Artificial intelligence, or the discipline of developing computational algorithms able to perform tasks that requires human intelligence, offers the opportunity to improve our idea and delivery of precision medicine. Here, we provide an overview of artificial intelligence approaches for the analysis of large-scale RNA-sequencing datasets in cancer. We present the major solutions to disentangle inter- and intra-tumor heterogeneity of transcriptome profiles for an effective improvement of patient management. We outline the contributions of learning algorithms to the needs of cancer genomics, from identifying rare cancer subtypes to personalizing therapeutic treatments.
Collapse
Affiliation(s)
- Marco Del Giudice
- Cancer Genomics and Bioinformatics Unit, IIGM—Italian Institute for Genomic Medicine, c/o IRCCS, Str. Prov.le 142, km 3.95, 10060 Candiolo, TO, Italy; (M.D.G.); (S.P.); (S.P.); (F.P.); (F.V.)
- Candiolo Cancer Institute, FPO—IRCCS, Str. Prov.le 142, km 3.95, 10060 Candiolo, TO, Italy
| | - Serena Peirone
- Cancer Genomics and Bioinformatics Unit, IIGM—Italian Institute for Genomic Medicine, c/o IRCCS, Str. Prov.le 142, km 3.95, 10060 Candiolo, TO, Italy; (M.D.G.); (S.P.); (S.P.); (F.P.); (F.V.)
- Department of Physics and INFN, Università degli Studi di Torino, via P.Giuria 1, 10125 Turin, Italy
| | - Sarah Perrone
- Cancer Genomics and Bioinformatics Unit, IIGM—Italian Institute for Genomic Medicine, c/o IRCCS, Str. Prov.le 142, km 3.95, 10060 Candiolo, TO, Italy; (M.D.G.); (S.P.); (S.P.); (F.P.); (F.V.)
- Department of Physics, Università degli Studi di Torino, via P.Giuria 1, 10125 Turin, Italy
| | - Francesca Priante
- Cancer Genomics and Bioinformatics Unit, IIGM—Italian Institute for Genomic Medicine, c/o IRCCS, Str. Prov.le 142, km 3.95, 10060 Candiolo, TO, Italy; (M.D.G.); (S.P.); (S.P.); (F.P.); (F.V.)
- Department of Physics, Università degli Studi di Torino, via P.Giuria 1, 10125 Turin, Italy
| | - Fabiola Varese
- Cancer Genomics and Bioinformatics Unit, IIGM—Italian Institute for Genomic Medicine, c/o IRCCS, Str. Prov.le 142, km 3.95, 10060 Candiolo, TO, Italy; (M.D.G.); (S.P.); (S.P.); (F.P.); (F.V.)
- Department of Life Science and System Biology, Università degli Studi di Torino, via Accademia Albertina 13, 10123 Turin, Italy
| | - Elisa Tirtei
- Paediatric Onco-Haematology Division, Regina Margherita Children’s Hospital, City of Health and Science of Turin, 10126 Turin, Italy; (E.T.); (F.F.)
| | - Franca Fagioli
- Paediatric Onco-Haematology Division, Regina Margherita Children’s Hospital, City of Health and Science of Turin, 10126 Turin, Italy; (E.T.); (F.F.)
- Department of Public Health and Paediatric Sciences, University of Torino, 10124 Turin, Italy
| | - Matteo Cereda
- Cancer Genomics and Bioinformatics Unit, IIGM—Italian Institute for Genomic Medicine, c/o IRCCS, Str. Prov.le 142, km 3.95, 10060 Candiolo, TO, Italy; (M.D.G.); (S.P.); (S.P.); (F.P.); (F.V.)
- Candiolo Cancer Institute, FPO—IRCCS, Str. Prov.le 142, km 3.95, 10060 Candiolo, TO, Italy
- Correspondence: ; Tel.: +39-011-993-3969
| |
Collapse
|
50
|
Durham TJ, Daza RM, Gevirtzman L, Cusanovich DA, Bolonduro O, Noble WS, Shendure J, Waterston RH. Comprehensive characterization of tissue-specific chromatin accessibility in L2 Caenorhabditis elegans nematodes. Genome Res 2021; 31:1952-1969. [PMID: 33888511 DOI: 10.1101/gr.271791.120] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/15/2020] [Accepted: 04/13/2021] [Indexed: 11/24/2022]
Abstract
Recently developed single-cell technologies allow researchers to characterize cell states at ever greater resolution and scale. Caenorhabditis elegans is a particularly tractable system for studying development, and recent single-cell RNA-seq studies characterized the gene expression patterns for nearly every cell type in the embryo and at the second larval stage (L2). Gene expression patterns give insight about gene function and into the biochemical state of different cell types; recent advances in other single-cell genomics technologies can now also characterize the regulatory context of the genome that gives rise to these gene expression levels at a single-cell resolution. To explore the regulatory DNA of individual cell types in C. elegans, we collected single-cell chromatin accessibility data using the sci-ATAC-seq assay in L2 larvae to match the available single-cell RNA-seq data set. By using a novel implementation of the latent Dirichlet allocation algorithm, we identify 37 clusters of cells that correspond to different cell types in the worm, providing new maps of putative cell type-specific gene regulatory sites, with promise for better understanding of cellular differentiation and gene regulation.
Collapse
Affiliation(s)
- Timothy J Durham
- Department of Genome Sciences, University of Washington, Seattle, Washington 98195, USA
| | - Riza M Daza
- Department of Genome Sciences, University of Washington, Seattle, Washington 98195, USA
| | - Louis Gevirtzman
- Department of Genome Sciences, University of Washington, Seattle, Washington 98195, USA
| | - Darren A Cusanovich
- Department of Cellular and Molecular Medicine, University of Arizona, Tucson, Arizona 85721, USA
| | - Olubusayo Bolonduro
- Department of Genome Sciences, University of Washington, Seattle, Washington 98195, USA
| | - William Stafford Noble
- Department of Genome Sciences, University of Washington, Seattle, Washington 98195, USA.,Paul G. Allen School of Computer Science and Engineering, University of Washington, Seattle, Washington 98195, USA
| | - Jay Shendure
- Department of Genome Sciences, University of Washington, Seattle, Washington 98195, USA.,Howard Hughes Medical Institute, Chevy Chase, Maryland 20815, USA.,Brotman Baty Institute for Precision Medicine, Seattle, Washington 98195, USA.,Allen Discovery Center for Cell Lineage Tracing, University of Washington, Seattle, Washington 98195, USA
| | - Robert H Waterston
- Department of Genome Sciences, University of Washington, Seattle, Washington 98195, USA
| |
Collapse
|