1
|
Missarova A, Dann E, Rosen L, Satija R, Marioni J. Leveraging neighborhood representations of single-cell data to achieve sensitive DE testing with miloDE. Genome Biol 2024; 25:189. [PMID: 39026254 PMCID: PMC11256449 DOI: 10.1186/s13059-024-03334-3] [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: 09/01/2023] [Accepted: 07/10/2024] [Indexed: 07/20/2024] Open
Abstract
Single-cell RNA-sequencing enables testing for differential expression (DE) between conditions at a cell type level. While powerful, one of the limitations of such approaches is that the sensitivity of DE testing is dictated by the sensitivity of clustering, which is often suboptimal. To overcome this, we present miloDE-a cluster-free framework for DE testing (available as an open-source R package). We illustrate the performance of miloDE on both simulated and real data. Using miloDE, we identify a transient hemogenic endothelia-like state in mouse embryos lacking Tal1 and detect distinct programs during macrophage activation in idiopathic pulmonary fibrosis.
Collapse
Affiliation(s)
- Alsu Missarova
- European Molecular Biology Laboratory, European Bioinformatics Institute, Wellcome Genome Campus, Cambridge, UK
- Cancer Research UK Cambridge Institute, University of Cambridge, Cambridge, UK
| | - Emma Dann
- Wellcome Sanger Institute, Wellcome Genome Campus, Hinxton, Cambridge, UK
| | - Leah Rosen
- European Molecular Biology Laboratory, European Bioinformatics Institute, Wellcome Genome Campus, Cambridge, UK
- Cancer Research UK Cambridge Institute, University of Cambridge, Cambridge, UK
| | - Rahul Satija
- Center for Genomics and Systems Biology, NYU, New York, USA.
- New York Genome Center, New York, USA.
| | - John Marioni
- European Molecular Biology Laboratory, European Bioinformatics Institute, Wellcome Genome Campus, Cambridge, UK.
- Cancer Research UK Cambridge Institute, University of Cambridge, Cambridge, UK.
| |
Collapse
|
2
|
Russell D, Pisu D, Mattila J, Johnston L. CD38+ Alveolar macrophages mediate early control of M. tuberculosis proliferation in the lung. RESEARCH SQUARE 2024:rs.3.rs-3934768. [PMID: 39070650 PMCID: PMC11275981 DOI: 10.21203/rs.3.rs-3934768/v1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 07/30/2024]
Abstract
Tuberculosis, caused by M.tuberculosis (Mtb), remains an enduring global health challenge, especially given the limited efficacy of current therapeutic interventions. Much of existing research has focused on immune failure as a driver of tuberculosis. However, the crucial role of host macrophage biology in controlling the disease remains underappreciated. While we have gained deeper insights into how alveolar macrophages (AMs) interact with Mtb, the precise AM subsets that mediate protection and potentially prevent tuberculosis progression have yet to be identified. In this study, we employed multi-modal scRNA-seq analyses to evaluate the functional roles of diverse macrophage subpopulations across different infection timepoints, allowing us to delineate the dynamic landscape of controller and permissive AM populations during the course of infection. Our analyses at specific time-intervals post-Mtb challenge revealed macrophage populations transitioning between distinct anti- and pro-inflammatory states. Notably, early in Mtb infection, CD38- AMs showed a muted response. As infection progressed, we observed a phenotypic shift in AMs, with CD38+ monocyte-derived AMs (moAMs) and a subset of tissue-resident AMs (TR-AMs) emerging as significant controllers of bacterial growth. Furthermore, scATAC-seq analysis of naïve lungs demonstrated that CD38+ TR-AMs possessed a distinct chromatin signature prior to infection, indicative of epigenetic priming and predisposition to a pro-inflammatory response. BCG intranasal immunization increased the numbers of CD38+ macrophages, substantially enhancing their capability to restrict Mtb growth. Collectively, our findings emphasize the pivotal, dynamic roles of different macrophage subsets in TB infection and reveal rational pathways for the development of improved vaccines and immunotherapeutic strategies.
Collapse
|
3
|
Amancherla K, Schlendorf KH, Chow N, Sheng Q, Freedman JE, Rathmell JC. Single-cell RNA-sequencing identifies unique cell-specific gene expression profiles in high-grade cardiac allograft vasculopathy. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.07.10.602989. [PMID: 39026730 PMCID: PMC11257508 DOI: 10.1101/2024.07.10.602989] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 07/20/2024]
Abstract
Background Cardiac allograft vasculopathy (CAV), a diffuse thickening of the intima of the coronary arteries and microvasculature, is the leading cause of late graft failure and mortality after heart transplantation (HT). Diagnosis involves invasive coronary angiography, which carries substantial risk, and minimally-invasive approaches to CAV diagnosis are urgently needed. Using single-cell RNA-sequencing in peripheral blood mononuclear cells (PBMCs), we sought to identify cell-specific gene expression profiles in CAV. Methods Whole blood was collected from 22 HT recipients with angiographically-confirmed CAV and 18 HT recipients without CAV. PBMCs were isolated and subjected to single-cell RNA-sequencing using a 10X Genomics microfluidic platform. Downstream analyses focused on differential expression of genes, cell compositional changes, and T cell receptor repertoire analyses. Results Across 40 PBMC samples, we isolated 134,984 cells spanning 8 major clusters and 31 subclusters of cell types. Compositional analyses showed subtle, but significant increases in CD4+ T central memory cells, and CD14+ and CD16+ monocytes in high-grade CAV (CAV-2 and CAV-3) as compared to low-grade or absent CAV. After adjusting for age, gender, and prednisone use, 745 genes were differentially expressed in a cell-specific manner in high-grade CAV. Weighted gene co-expression network analyses showed enrichment for putative pathways involved in inflammation and angiogenesis. There were no significant differences in T cell clonality or diversity with increasing CAV severity. Conclusions Unbiased whole transcriptomic analyses at single-cell resolution identify unique, cell-specific gene expression patterns in CAV, suggesting the potential utility of peripheral gene expression biomarkers in diagnosing CAV.
Collapse
|
4
|
Bilous M, Hérault L, Gabriel AA, Teleman M, Gfeller D. Building and analyzing metacells in single-cell genomics data. Mol Syst Biol 2024; 20:744-766. [PMID: 38811801 PMCID: PMC11220014 DOI: 10.1038/s44320-024-00045-6] [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: 02/04/2024] [Revised: 05/03/2024] [Accepted: 05/08/2024] [Indexed: 05/31/2024] Open
Abstract
The advent of high-throughput single-cell genomics technologies has fundamentally transformed biological sciences. Currently, millions of cells from complex biological tissues can be phenotypically profiled across multiple modalities. The scaling of computational methods to analyze and visualize such data is a constant challenge, and tools need to be regularly updated, if not redesigned, to cope with ever-growing numbers of cells. Over the last few years, metacells have been introduced to reduce the size and complexity of single-cell genomics data while preserving biologically relevant information and improving interpretability. Here, we review recent studies that capitalize on the concept of metacells-and the many variants in nomenclature that have been used. We further outline how and when metacells should (or should not) be used to analyze single-cell genomics data and what should be considered when analyzing such data at the metacell level. To facilitate the exploration of metacells, we provide a comprehensive tutorial on the construction and analysis of metacells from single-cell RNA-seq data ( https://github.com/GfellerLab/MetacellAnalysisTutorial ) as well as a fully integrated pipeline to rapidly build, visualize and evaluate metacells with different methods ( https://github.com/GfellerLab/MetacellAnalysisToolkit ).
Collapse
Affiliation(s)
- Mariia Bilous
- Department of Oncology, Ludwig Institute for Cancer Research Lausanne, University of Lausanne, 1011, Lausanne, Switzerland
- Agora Cancer Research Centre, 1011, Lausanne, Switzerland
- Swiss Cancer Center Leman (SCCL), Lausanne, Switzerland
- Swiss Institute of Bioinformatics (SIB), 1015, Lausanne, Switzerland
| | - Léonard Hérault
- Department of Oncology, Ludwig Institute for Cancer Research Lausanne, University of Lausanne, 1011, Lausanne, Switzerland
- Agora Cancer Research Centre, 1011, Lausanne, Switzerland
- Swiss Cancer Center Leman (SCCL), Lausanne, Switzerland
- Swiss Institute of Bioinformatics (SIB), 1015, Lausanne, Switzerland
| | - Aurélie Ag Gabriel
- Department of Oncology, Ludwig Institute for Cancer Research Lausanne, University of Lausanne, 1011, Lausanne, Switzerland
- Agora Cancer Research Centre, 1011, Lausanne, Switzerland
- Swiss Cancer Center Leman (SCCL), Lausanne, Switzerland
- Swiss Institute of Bioinformatics (SIB), 1015, Lausanne, Switzerland
| | - Matei Teleman
- Department of Oncology, Ludwig Institute for Cancer Research Lausanne, University of Lausanne, 1011, Lausanne, Switzerland
- Agora Cancer Research Centre, 1011, Lausanne, Switzerland
- Swiss Cancer Center Leman (SCCL), Lausanne, Switzerland
- Swiss Institute of Bioinformatics (SIB), 1015, Lausanne, Switzerland
| | - David Gfeller
- Department of Oncology, Ludwig Institute for Cancer Research Lausanne, University of Lausanne, 1011, Lausanne, Switzerland.
- Agora Cancer Research Centre, 1011, Lausanne, Switzerland.
- Swiss Cancer Center Leman (SCCL), Lausanne, Switzerland.
- Swiss Institute of Bioinformatics (SIB), 1015, Lausanne, Switzerland.
| |
Collapse
|
5
|
Wang W, Cen Y, Lu Z, Xu Y, Sun T, Xiao Y, Liu W, Li JJ, Wang C. scCDC: a computational method for gene-specific contamination detection and correction in single-cell and single-nucleus RNA-seq data. Genome Biol 2024; 25:136. [PMID: 38783325 PMCID: PMC11112958 DOI: 10.1186/s13059-024-03284-w] [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/30/2023] [Accepted: 05/16/2024] [Indexed: 05/25/2024] Open
Abstract
In droplet-based single-cell and single-nucleus RNA-seq assays, systematic contamination of ambient RNA molecules biases the quantification of gene expression levels. Existing methods correct the contamination for all genes globally. However, there lacks specific evaluation of correction efficacy for varying contamination levels. Here, we show that DecontX and CellBender under-correct highly contaminating genes, while SoupX and scAR over-correct lowly/non-contaminating genes. Here, we develop scCDC as the first method to detect the contamination-causing genes and only correct expression levels of these genes, some of which are cell-type markers. Compared with existing decontamination methods, scCDC excels in decontaminating highly contaminating genes while avoiding over-correction of other genes.
Collapse
Affiliation(s)
- Weijian Wang
- Centre of Biomedical Systems and Informatics, International Campus, ZJU-UoE Institute, Zhejiang University School of Medicine, Zhejiang University, Haining, Zhejiang, 314400, China
| | - Yihui Cen
- Centre of Biomedical Systems and Informatics, International Campus, ZJU-UoE Institute, Zhejiang University School of Medicine, Zhejiang University, Haining, Zhejiang, 314400, China
| | - Zezhen Lu
- Centre of Biomedical Systems and Informatics, International Campus, ZJU-UoE Institute, Zhejiang University School of Medicine, Zhejiang University, Haining, Zhejiang, 314400, China
| | - Yueqing Xu
- Centre of Biomedical Systems and Informatics, International Campus, ZJU-UoE Institute, Zhejiang University School of Medicine, Zhejiang University, Haining, Zhejiang, 314400, China
| | - Tianyi Sun
- Department of Statistics and Data Science, University of California, Los Angeles, CA, 90095, USA
| | - Ying Xiao
- Sir Run Run Shaw Hospital, School of Medicine, Zhejiang University, Hangzhou, Zhejiang, 310020, China
| | - Wanlu Liu
- Centre of Biomedical Systems and Informatics, International Campus, ZJU-UoE Institute, Zhejiang University School of Medicine, Zhejiang University, Haining, Zhejiang, 314400, China
| | - Jingyi Jessica Li
- Department of Statistics and Data Science, University of California, Los Angeles, CA, 90095, USA.
| | - Chaochen Wang
- Centre of Biomedical Systems and Informatics, International Campus, ZJU-UoE Institute, Zhejiang University School of Medicine, Zhejiang University, Haining, Zhejiang, 314400, China.
- Department of Gynecology, The Second Affiliated Hospital, Zhejiang University School of Medicine, Zhejiang University, Hangzhou, Zhejiang, 310020, China.
- Biomedical and Health Translational Research Centre, Zhejiang University, Haining, Zhejiang, 314400, China.
| |
Collapse
|
6
|
Konopka G, Bhaduri A. Functional genomics and systems biology in human neuroscience. Nature 2023; 623:274-282. [PMID: 37938705 DOI: 10.1038/s41586-023-06686-1] [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: 02/23/2023] [Accepted: 09/27/2023] [Indexed: 11/09/2023]
Abstract
Neuroscience research has entered a phase of key discoveries in the realm of neurogenomics owing to strong financial and intellectual support for resource building and tool development. The previous challenge of tissue heterogeneity has been met with the application of techniques that can profile individual cells at scale. Moreover, the ability to perturb genes, gene regulatory elements and neuronal activity in a cell-type-specific manner has been integrated with gene expression studies to uncover the functional underpinnings of the genome at a systems level. Although these insights have necessarily been grounded in model systems, we now have the opportunity to apply these approaches in humans and in human tissue, thanks to advances in human genetics, brain imaging and tissue collection. We acknowledge that there will probably always be limits to the extent to which we can apply the genomic tools developed in model systems to human neuroscience; however, as we describe in this Perspective, the neuroscience field is now primed with an optimal foundation for tackling this ambitious challenge. The application of systems-level network analyses to these datasets will facilitate a deeper appreciation of human neurogenomics that cannot otherwise be achieved from directly observable phenomena.
Collapse
Affiliation(s)
- Genevieve Konopka
- Department of Neuroscience, UT Southwestern Medical Center, Dallas, TX, USA.
- Peter O'Donnell Jr Brain Institute, UT Southwestern Medical Center, Dallas, TX, USA.
| | - Aparna Bhaduri
- Department of Biological Chemistry, University of California, Los Angeles, CA, USA.
| |
Collapse
|
7
|
Rao SB, Brundu F, Chen Y, Sun Y, Zhu H, Shprintzen RJ, Tomer R, Rabadan R, Leong KW, Markx S, Xu B, Gogos JA. Aberrant pace of cortical neuron development in brain organoids from patients with 22q11.2 deletion syndrome and schizophrenia. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.10.04.557612. [PMID: 37873382 PMCID: PMC10592956 DOI: 10.1101/2023.10.04.557612] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/25/2023]
Abstract
Adults and children afflicted with the 22q11.2 deletion syndrome (22q11.2DS) exhibit cognitive, social, and emotional impairments, and are at significantly heightened risk for schizophrenia (SCZ). The impact of this deletion on early human brain development, however, has remained unclear. Here we harness organoid models of the developing human cerebral cortex, cultivated from subjects with 22q11.2DS and SCZ, as well as unaffected control samples, to identify cell-type-specific developmental abnormalities arising from this genomic lesion. Leveraging single-cell RNA-sequencing in conjunction with experimental validation, we find that the loss of genes within the 22q11.2 locus leads to a delayed development of cortical neurons. This compromised development was reflected in an elevated proportion of actively proliferating neural progenitor cells, coupled with a decreased fraction of more mature neurons. Furthermore, we identify perturbed molecular imprints linked to neuronal maturation, observe the presence of sparser neurites, and note a blunted amplitude in glutamate-induced Ca2+ transients. The aberrant transcription program underlying impaired development contains molecular signatures significantly enriched in neuropsychiatric genetic liability. MicroRNA profiling and target gene investigation suggest that microRNA dysregulation may drive perturbations of genes governing the pace at which maturation unfolds. Using protein-protein interaction network analysis we define complementary effects stemming from additional genes residing within the deleted locus. Our study uncovers reproducible neurodevelopmental and molecular alterations due to 22q11.2 deletions. These findings have the potential to facilitate disease modeling and promote the pursuit of therapeutic interventions.
Collapse
|
8
|
Zhang H, Lu X, Lu B, Chen L. scGEM: Unveiling the Nested Tree-Structured Gene Co-Expressing Modules in Single Cell Transcriptome Data. Cancers (Basel) 2023; 15:4277. [PMID: 37686554 PMCID: PMC10486867 DOI: 10.3390/cancers15174277] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/20/2023] [Revised: 08/22/2023] [Accepted: 08/25/2023] [Indexed: 09/10/2023] Open
Abstract
BACKGROUND Single-cell transcriptome analysis has fundamentally changed biological research by allowing higher-resolution computational analysis of individual cells and subsets of cell types. However, few methods have met the need to recognize and quantify the underlying cellular programs that determine the specialization and differentiation of the cell types. METHODS In this study, we present scGEM, a nested tree-structured nonparametric Bayesian model, to reveal the gene co-expression modules (GEMs) reflecting transcriptome processes in single cells. RESULTS We show that scGEM can discover shared and specialized transcriptome signals across different cell types using peripheral blood mononuclear single cells and early brain development single cells. scGEM outperformed other methods in perplexity and topic coherence (p < 0.001) on our simulation data. Larger datasets, deeper trees and pre-trained models are shown to be positively associated with better scGEM performance. The GEMs obtained from triple-negative breast cancer single cells exhibited better correlations with lymphocyte infiltration (p = 0.009) and the cell cycle (p < 0.001) than other methods in additional validation on the bulk RNAseq dataset. CONCLUSIONS Altogether, we demonstrate that scGEM can be used to model the hidden cellular functions of single cells, thereby unveiling the specialization and generalization of transcriptomic programs across different types of cells.
Collapse
Affiliation(s)
- Han Zhang
- Department of Biomedical Informatics, University of Pittsburgh, Pittsburgh, PA 15206, USA; (H.Z.)
| | - Xinghua Lu
- Department of Biomedical Informatics, University of Pittsburgh, Pittsburgh, PA 15206, USA; (H.Z.)
- UPMC Hillman Cancer Center, Pittsburgh, PA 15232, USA
| | - Binfeng Lu
- Center for Discovery and Innovation, Hackensack Meridian Health, Nutley, NJ 07110, USA
| | - Lujia Chen
- Department of Biomedical Informatics, University of Pittsburgh, Pittsburgh, PA 15206, USA; (H.Z.)
| |
Collapse
|
9
|
Zeng H, Huang J, Zhou H, Meilandt WJ, Dejanovic B, Zhou Y, Bohlen CJ, Lee SH, Ren J, Liu A, Tang Z, Sheng H, Liu J, Sheng M, Wang X. Integrative in situ mapping of single-cell transcriptional states and tissue histopathology in a mouse model of Alzheimer's disease. Nat Neurosci 2023; 26:430-446. [PMID: 36732642 PMCID: PMC11332722 DOI: 10.1038/s41593-022-01251-x] [Citation(s) in RCA: 15] [Impact Index Per Article: 15.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/29/2021] [Accepted: 12/20/2022] [Indexed: 02/04/2023]
Abstract
Complex diseases are characterized by spatiotemporal cellular and molecular changes that may be difficult to comprehensively capture. However, understanding the spatiotemporal dynamics underlying pathology can shed light on disease mechanisms and progression. Here we introduce STARmap PLUS, a method that combines high-resolution spatial transcriptomics with protein detection in the same tissue section. As proof of principle, we analyze brain tissues of a mouse model of Alzheimer's disease at 8 and 13 months of age. Our approach provides a comprehensive cellular map of disease progression. It reveals a core-shell structure where disease-associated microglia (DAM) closely contact amyloid-β plaques, whereas disease-associated astrocyte-like (DAA-like) cells and oligodendrocyte precursor cells (OPCs) are enriched in the outer shells surrounding the plaque-DAM complex. Hyperphosphorylated tau emerges mainly in excitatory neurons in the CA1 region and correlates with the local enrichment of oligodendrocyte subtypes. The STARmap PLUS method bridges single-cell gene expression profiles with tissue histopathology at subcellular resolution, providing a tool to pinpoint the molecular and cellular changes underlying pathology.
Collapse
Affiliation(s)
- Hu Zeng
- Broad Institute of MIT and Harvard, Cambridge, MA, USA
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, MA, USA
| | - Jiahao Huang
- Broad Institute of MIT and Harvard, Cambridge, MA, USA
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, MA, USA
| | - Haowen Zhou
- Broad Institute of MIT and Harvard, Cambridge, MA, USA
| | - William J Meilandt
- Department of Neuroscience, Genentech, Inc., South San Francisco, CA, USA
| | - Borislav Dejanovic
- Broad Institute of MIT and Harvard, Cambridge, MA, USA
- Stanley Center for Psychiatric Research, Broad Institute of MIT and Harvard, Cambridge, MA, USA
| | - Yiming Zhou
- Broad Institute of MIT and Harvard, Cambridge, MA, USA
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, MA, USA
| | | | - Seung-Hye Lee
- Department of Neuroscience, Genentech, Inc., South San Francisco, CA, USA
| | - Jingyi Ren
- Broad Institute of MIT and Harvard, Cambridge, MA, USA
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, MA, USA
| | - Albert Liu
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, MA, USA
| | - Zefang Tang
- Broad Institute of MIT and Harvard, Cambridge, MA, USA
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, MA, USA
| | - Hao Sheng
- John A. Paulson School of Engineering and Applied Sciences, Harvard University, Boston, MA, USA
| | - Jia Liu
- John A. Paulson School of Engineering and Applied Sciences, Harvard University, Boston, MA, USA
| | - Morgan Sheng
- Broad Institute of MIT and Harvard, Cambridge, MA, USA.
- Stanley Center for Psychiatric Research, Broad Institute of MIT and Harvard, Cambridge, MA, USA.
- Department of Brain and Cognitive Sciences, Massachusetts Institute of Technology, Cambridge, MA, USA.
| | - Xiao Wang
- Broad Institute of MIT and Harvard, Cambridge, MA, USA.
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, MA, USA.
- Stanley Center for Psychiatric Research, Broad Institute of MIT and Harvard, Cambridge, MA, USA.
| |
Collapse
|
10
|
Identification of Transcription Factor Networks during Mouse Hindlimb Development. Cells 2022; 12:cells12010028. [PMID: 36611822 PMCID: PMC9818828 DOI: 10.3390/cells12010028] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/16/2022] [Revised: 12/17/2022] [Accepted: 12/20/2022] [Indexed: 12/24/2022] Open
Abstract
Mammalian hindlimb development involves a variety of cells and the regulation of spatiotemporal molecular events, but regulatory networks of transcription factors contributing to hindlimb morphogenesis are not well understood. Here, we identified transcription factor networks during mouse hindlimb morphology establishment through transcriptome analysis. We used four stages of embryonic hindlimb transcription profiles acquired from the Gene Expression Omnibus database (GSE30138), including E10.5, E11.5, E12.5 and E13.5, to construct a gene network using Weighted Gene Co-expression Network Analysis (WGCNA), and defined seven stage-associated modules. After filtering 7625 hub genes, we further prioritized 555 transcription factors with AnimalTFDB3.0. Gene ontology enrichment showed that transcription factors of different modules were enriched in muscle tissue development, connective tissue development, embryonic organ development, skeletal system morphogenesis, pattern specification process and urogenital system development separately. Six regulatory networks were constructed with key transcription factors, which contribute to the development of different tissues. Knockdown of four transcription factors from regulatory networks, including Sox9, Twist1, Snai2 and Klf4, showed that the expression of limb-development-related genes was also inhibited, which indicated the crucial role of transcription factor networks in hindlimb development.
Collapse
|
11
|
Zhou D, Zhou X, Sun C, Tang G, Liu L, Chen L, He H, Xiong Q. Transcriptome and Metabolome Analysis Provides Insights into the Heterosis of Yield and Quality Traits in Two Hybrid Rice Varieties (Oryza sativa L.). Int J Mol Sci 2022; 23:ijms232112934. [PMID: 36361748 PMCID: PMC9654843 DOI: 10.3390/ijms232112934] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/27/2022] [Revised: 10/18/2022] [Accepted: 10/20/2022] [Indexed: 11/24/2022] Open
Abstract
Heterosis is a common biological phenomenon that is useful for breeding superior lines. Using heterosis to increase the yield and quality of crops is one of the main achievements of modern agricultural science. In this study, we analysed the transcriptome and metabolome of two three-line hybrid rice varieties, Taiyou 871 (TY871), and Taiyou 398 (TY398) and the parental grain endosperm using RNA-seq (three biological repeats per variety) and untargeted metabolomic (six biological repeats per variety) methods. TY871 and TY398 showed specific heterosis in yield and quality. Transcriptome analysis of the hybrids revealed 638 to 4059 differentially expressed genes in the grain when compared to the parents. Metabolome analysis of the hybrids revealed 657 to 3714 differential grain metabolites when compared to the parents. The honeydew1 and grey60 module core genes Os04g0350700 and Os05g0154700 are involved in the regulation of awn development, grain size, and grain number, as well as the regulation of grain length and plant height, respectively. Rice grain length may be an important indicator for improving the quality of three-line hybrid rice. In addition, the rice quality-related metabolite NEG_M341T662 was highly connected to the module core genes Os06g0254300 and Os03g0168100. The functions of Os06g0254300 and Os03g0168100 are EF-hand calcium binding protein and late embroideries absolute protein repeat containing protein, respectively. These genes may play a role in the formation of rice quality. We constructed a gene and metabolite coexpression network, which provides a scientific basis for the utilization of heterosis in producing high-yield and high-quality hybrid rice.
Collapse
Affiliation(s)
- Dahu Zhou
- Key Laboratory of Crop Physiology, Ecology and Genetic Breeding, Ministry of Education, College of Agronomy, Jiangxi Agricultural University, Nanchang 330045, China
| | - Xinyi Zhou
- Key Laboratory of Crop Physiology, Ecology and Genetic Breeding, Ministry of Education, College of Agronomy, Jiangxi Agricultural University, Nanchang 330045, China
| | - Changhui Sun
- Jiangsu Key Laboratory of Crop Genetics and Physiology/Jiangsu Key Laboratory of Crop Cultivation and Physiology, Agricultural College of Yangzhou University, Yangzhou 225009, China
| | - Guoping Tang
- Jiangxi Academy of Agricultural Sciences Rice Research Institute, Nanchang 330200, China
| | - Lin Liu
- Key Laboratory of Crop Physiology, Ecology and Genetic Breeding, Ministry of Education, College of Agronomy, Jiangxi Agricultural University, Nanchang 330045, China
| | - Le Chen
- Jiangsu Key Laboratory of Crop Genetics and Physiology/Jiangsu Key Laboratory of Crop Cultivation and Physiology, Agricultural College of Yangzhou University, Yangzhou 225009, China
| | - Haohua He
- Key Laboratory of Crop Physiology, Ecology and Genetic Breeding, Ministry of Education, College of Agronomy, Jiangxi Agricultural University, Nanchang 330045, China
- Correspondence: (H.H.); (Q.X.)
| | - Qiangqiang Xiong
- Jiangsu Key Laboratory of Crop Genetics and Physiology/Jiangsu Key Laboratory of Crop Cultivation and Physiology, Agricultural College of Yangzhou University, Yangzhou 225009, China
- Correspondence: (H.H.); (Q.X.)
| |
Collapse
|
12
|
Wang M, Bi C, Li H, Lu L, Gao T, Huang P, Liu C, Wang B. The emerging double-edged sword role of Sirtuins in the gastric inflammation-carcinoma sequence revealed by bulk and single-cell transcriptomes. Front Oncol 2022; 12:1004726. [PMID: 36324577 PMCID: PMC9619065 DOI: 10.3389/fonc.2022.1004726] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/27/2022] [Accepted: 09/12/2022] [Indexed: 12/03/2022] Open
Abstract
Histone modification and the inflammation-carcinoma sequence (ICS) have been acknowledgedly implicated in gastric carcinogenesis. However, the extremum expression of some histone modification genes (HMGs) in intestinal metaplasia (IM) rather than GC obscures the roles of HMGs in ICS. In this study, we assumed an explanation that the roles of HMGs in ICS were stage specific. Bulk RNA-seq on endoscopy biopsy samples from a total of 50 patients was accompanied by reanalysis of a set of published single-cell transcriptomes, which cross-sectionally profiled the transcriptomic features of chronic superficial gastritis (SG), atrophy gastritis (AG), IM, and early gastric cancer (GC). Differential analysis observed significantly peaked expression of SIRT6 and SIRT7 at IM. Weighted correlation network analysis on bulk transcriptome recognized significant correlations between SIRT1/6 and IM. The single-cell atlas identified one subgroup of B cells expressing high level of TFF1 (TFF1hi naive B cell) that theoretically played important roles in defending microbial infection, while SIRT6 displayed a positive correlation with TFF1low naive B cells. Moreover, gene set enrichment analysis at different lesions (SG-AG, AG-IM, and IM-GC) highlighted that gene sets contributing to IM, e.g., Brush Border, were largely enriched from co-expressing genes of Sirtuins (SIRTs) in AG-IM. Surveys of the genes negatively correlated with SIRT6 in public databases considered SIRT6 as tumor suppressors, which was confirmed by the cell proliferation and migration assays after transient transfection of SIRT6 overexpression vector into AGS cells. All the above observations were then confirmed by serial section-based immunohistochemistry against Ki-67, MUC2, MUC5AC, p53, and SIRT6 on the endoscopic submucosal dissection tissue. By contrast, the expression of the other HMGs varied even opposite within same family. Taken together, this study preliminarily demonstrated the two-edged sword role of SIRTs in ICS and, by extension, showed that the roles of HMGs in ICS were probably stage specific. Our study may provide new insights into and attract attention on gastric prevention and therapy targeting HMGs.
Collapse
Affiliation(s)
- Mengyang Wang
- Department of Immunology, Binzhou Medical University, Yantai, China
- Department of Gastroenterology, Binzhou Medical University Hospital, Binzhou, China
| | - Chenxiao Bi
- Department of Immunology, Binzhou Medical University, Yantai, China
- Department of Gastroenterology, Binzhou Medical University Hospital, Binzhou, China
| | - Hong Li
- Department of Pathology, Binzhou Medical University Hospital, Binzhou, China
| | - Lizhen Lu
- Department of Pathology, Binzhou Medical University Hospital, Binzhou, China
| | - Tao Gao
- Department of Gastroenterology, Binzhou Medical University Hospital, Binzhou, China
| | - Panpan Huang
- Department of Immunology, Binzhou Medical University, Yantai, China
| | - Chengxia Liu
- Department of Gastroenterology, Binzhou Medical University Hospital, Binzhou, China
- *Correspondence: Chengxia Liu, ; Bin Wang,
| | - Bin Wang
- Department of Immunology, Binzhou Medical University, Yantai, China
- *Correspondence: Chengxia Liu, ; Bin Wang,
| |
Collapse
|
13
|
Gros J, Davey M, Sears KE. Limb development, evolution, and regeneration and repair: Part two. Dev Dyn 2022; 251:1387-1388. [PMID: 36052833 DOI: 10.1002/dvdy.521] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 08/12/2022] [Indexed: 11/08/2022] Open
Affiliation(s)
- Jérôme Gros
- Université Paris Cité, CNRS UMR3738, Developmental and Stem Cell Biology Department, Institut Pasteur, Paris, France
| | - Megan Davey
- Functional Genetics and Development, The Roslin Institute, The University of Edinburgh, Edinburgh, Scotland
| | - Karen E Sears
- Departments of Ecology and Evolutionary Biology and Molecular, Cell, and Developmental Biology, UCLA, Los Angeles, California, USA
| |
Collapse
|
14
|
Zhang T, Wong G. Gene expression data analysis using Hellinger correlation in weighted gene co-expression networks (WGCNA). Comput Struct Biotechnol J 2022; 20:3851-3863. [PMID: 35891798 PMCID: PMC9307959 DOI: 10.1016/j.csbj.2022.07.018] [Citation(s) in RCA: 18] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/13/2022] [Revised: 07/09/2022] [Accepted: 07/09/2022] [Indexed: 12/24/2022] Open
Abstract
Weighted gene co-expression network analysis (WGCNA) is used to detect clusters with highly correlated genes. Measurements of correlation most typically rely on linear relationships. However, a linear relationship does not always model pairwise functional-related dependence between genes. In this paper, we first compared 6 different correlation methods in their ability to capture complex dependence between genes in three different tissues. Next, we compared their gene-pairwise coefficient results and corresponding WGCNA results. Finally, we applied a recently proposed correlation method, Hellinger correlation, as a more sensitive correlation measurement in WGCNA. To test this method, we constructed gene networks containing co-expression gene modules from RNA-seq data of human frontal cortex from Alzheimer's disease patients. To test the generality, we also used a microarray data set from human frontal cortex, single cell RNA-seq data from human prefrontal cortex, RNA-seq data from human temporal cortex, and GTEx data from heart. The Hellinger correlation method captures essentially similar results as other linear correlations in WGCNA, but provides additional new functional relationships as exemplified by uncovering a link between inflammation and mitochondria function. We validated the network constructed with the microarray and single cell sequencing data sets and a RNA-seq dataset of temporal cortex. We observed that this new correlation method enables the detection of non-linear biologically meaningful relationships among genes robustly and provides a complementary new approach to WGCNA. Thus, the application of Hellinger correlation to WGCNA provides a more flexible correlation approach to modelling networks in gene expression analysis that uncovers novel network relationships.
Collapse
Affiliation(s)
- Tianjiao Zhang
- Cancer Centre, Centre for Reproduction, Development and Aging, Department of Public Health and Medicinal Administration, Faculty of Health Sciences, University of Macau, Taipa 999078, Macau Special Administrative Region
| | - Garry Wong
- Cancer Centre, Centre for Reproduction, Development and Aging, Department of Public Health and Medicinal Administration, Faculty of Health Sciences, University of Macau, Taipa 999078, Macau Special Administrative Region
| |
Collapse
|
15
|
Li A, Yang J, Qian J, Shao X, Liao J, Lu X, Fan X. Tracing the cell-type-specific modules of immune responses during COVID-19 progression using scDisProcema. Comput Struct Biotechnol J 2022; 20:3545-3555. [PMID: 35811838 PMCID: PMC9250167 DOI: 10.1016/j.csbj.2022.06.066] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/14/2022] [Revised: 06/30/2022] [Accepted: 06/30/2022] [Indexed: 12/16/2022] Open
Abstract
COVID-19 has caused severe threats to lives and damage to property worldwide. The immunopathology of the disease is of particular concern. Currently, researchers have used gene co-expression networks (GCNs) to deepen the study of molecular mechanisms of immune responses to COVID-19. However, most efforts have not fully explored dynamic changes of cell-type-specific molecular networks in the disease process. This study proposes a GCN construction pipeline named single-cell Disease Progression cellular module analysis (scDisProcema), which can trace dynamic changes of immune system response during disease progression using single-cell data. Here, scDisProcema considers changes in cell fate and expression patterns during disease development, identifying gene modules responsible for different immune cells. The hub genes are screened for each module by the specific expression level and the intercellular connectivity of modules. Based on functional items enriched by each gene module, we elucidate the biological processes of different cells involved in disease development and explain the molecular mechanisms underlying the process of cell depletion or proliferation caused by disease. Compared with traditional WGCNA methods, scDisProcema can make more convenient use of the heterogeneity information provided by scRNA-seq data and has great potential in exploring molecular changes during disease progression and organ development.
Collapse
Affiliation(s)
- Anyao Li
- College of Pharmaceutical Sciences, Zhejiang University, Hangzhou 310058, China
- Innovation Center of Yangtze River Delta, Zhejiang University, Jiaxing 314102, China
| | - Jihong Yang
- College of Pharmaceutical Sciences, Zhejiang University, Hangzhou 310058, China
- Innovation Center in Zhejiang University, State Key Laboratory of Component-Based Chinese Medicine, Hangzhou 310058, China
- Zhang Boli Intelligent Health Innovation Lab, Hangzhou 311100, China
| | - Jingyang Qian
- College of Pharmaceutical Sciences, Zhejiang University, Hangzhou 310058, China
- Innovation Center of Yangtze River Delta, Zhejiang University, Jiaxing 314102, China
| | - Xin Shao
- College of Pharmaceutical Sciences, Zhejiang University, Hangzhou 310058, China
- Innovation Center of Yangtze River Delta, Zhejiang University, Jiaxing 314102, China
| | - Jie Liao
- College of Pharmaceutical Sciences, Zhejiang University, Hangzhou 310058, China
- Innovation Center of Yangtze River Delta, Zhejiang University, Jiaxing 314102, China
| | - Xiaoyan Lu
- College of Pharmaceutical Sciences, Zhejiang University, Hangzhou 310058, China
- Innovation Center of Yangtze River Delta, Zhejiang University, Jiaxing 314102, China
- Innovation Center in Zhejiang University, State Key Laboratory of Component-Based Chinese Medicine, Hangzhou 310058, China
| | - Xiaohui Fan
- College of Pharmaceutical Sciences, Zhejiang University, Hangzhou 310058, China
- Innovation Center of Yangtze River Delta, Zhejiang University, Jiaxing 314102, China
- Innovation Center in Zhejiang University, State Key Laboratory of Component-Based Chinese Medicine, Hangzhou 310058, China
- Corresponding author at: College of Pharmaceutical Sciences, Zhejiang University, Hangzhou 310058, China.
| |
Collapse
|
16
|
Newton AH. Marsupials and Multi-Omics: Establishing New Comparative Models of Neural Crest Patterning and Craniofacial Development. Front Cell Dev Biol 2022; 10:941168. [PMID: 35813210 PMCID: PMC9260703 DOI: 10.3389/fcell.2022.941168] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/11/2022] [Accepted: 06/06/2022] [Indexed: 11/15/2022] Open
Abstract
Studies across vertebrates have revealed significant insights into the processes that drive craniofacial morphogenesis, yet we still know little about how distinct facial morphologies are patterned during development. Studies largely point to evolution in GRNs of cranial progenitor cell types such as neural crest cells, as the major driver underlying adaptive cranial shapes. However, this hypothesis requires further validation, particularly within suitable models amenable to manipulation. By utilizing comparative models between related species, we can begin to disentangle complex developmental systems and identify the origin of species-specific patterning. Mammals present excellent evolutionary examples to scrutinize how these differences arise, as sister clades of eutherians and marsupials possess suitable divergence times, conserved cranial anatomies, modular evolutionary patterns, and distinct developmental heterochrony in their NCC behaviours and craniofacial patterning. In this review, I lend perspectives into the current state of mammalian craniofacial biology and discuss the importance of establishing a new marsupial model, the fat-tailed dunnart, for comparative research. Through detailed comparisons with the mouse, we can begin to decipher mammalian conserved, and species-specific processes and their contribution to craniofacial patterning and shape disparity. Recent advances in single-cell multi-omics allow high-resolution investigations into the cellular and molecular basis of key developmental processes. As such, I discuss how comparative evolutionary application of these tools can provide detailed insights into complex cellular behaviours and expression dynamics underlying adaptive craniofacial evolution. Though in its infancy, the field of "comparative evo-devo-omics" presents unparalleled opportunities to precisely uncover how phenotypic differences arise during development.
Collapse
|