CCPLS reveals cell-type-specific spatial dependence of transcriptomes in single cells.
Ontology highlight
ABSTRACT:
SUBMITTER:
PROVIDER: S-EPMC9933841 | biostudies-literature | 2023 Feb
REPOSITORIES: biostudies-literature
Publications
CCPLS reveals cell-type-specific spatial dependence of transcriptomes in single cells.
Bioinformatics (Oxford, England) 20230201 2
Similar Datasets
Project description:MotivationCell-cell communications regulate internal cellular states, e.g. gene expression and cell functions, and play pivotal roles in normal development and disease states. Furthermore, single-cell RNA sequencing methods have revealed cell-to-cell expression variability of highly variable genes (HVGs), which is also crucial. Nevertheless, the regulation of cell-to-cell expression variability of HVGs via cell-cell communications is still largely unexplored. The recent advent of spatial transcriptome methods has linked gene expression profiles to the spatial context of single cells, which has provided opportunities to reveal those regulations. The existing computational methods extract genes with expression levels influenced by neighboring cell types. However, limitations remain in the quantitativeness and interpretability: they neither focus on HVGs nor consider the effects of multiple neighboring cell types.ResultsHere, we propose CCPLS (Cell-Cell communications analysis by Partial Least Square regression modeling), which is a statistical framework for identifying cell-cell communications as the effects of multiple neighboring cell types on cell-to-cell expression variability of HVGs, based on the spatial transcriptome data. For each cell type, CCPLS performs PLS regression modeling and reports coefficients as the quantitative index of the cell-cell communications. Evaluation using simulated data showed our method accurately estimated the effects of multiple neighboring cell types on HVGs. Furthermore, applications to the two real datasets demonstrate that CCPLS can extract biologically interpretable insights from the inferred cell-cell communications.Availability and implementationThe R package is available at https://github.com/bioinfo-tsukuba/CCPLS. The data are available at https://github.com/bioinfo-tsukuba/CCPLS_paper.Supplementary informationSupplementary data are available at Bioinformatics online.
| S-EPMC9620831 | biostudies-literature
Project description:Accurate annotation of gene function in individual samples and even in each cell type is essential for understanding the pathogenesis of cancers. Single-cell RNA-sequencing (scRNA-seq) provides unprecedented resolution to decipher gene function. In order to explore how scRNA-seq contributes to the understanding of gene function in cancers, we constructed an assessment framework based on co-expression network and neighbor-voting method using 116,814 cells. Compared with bulk transcriptome, scRNA-seq recalled more experimentally verified gene functions. Surprisingly, scRNA-seq revealed cell-type-specific functions, especially in immune cells, whose expression profile recalled immune-related functions that were not discovered in cancer cells. Furthermore, scRNA-seq discovered sample-specific functions, highlighting that it provided sample-specific information. We also explored factors affecting the performance of gene function prediction. We found that 500 or more cells should be considered in the prediction with scRNA-seq, and that scRNA-seq datasets generated from 10x Genomics platform had a better performance than those from Smart-seq2. Collectively, we compared the prediction performance of bulk data and scRNA-seq data from multiple perspectives, revealing the irreplaceable role of single-cell sequencing in decoding the biological progresses in which the gene involved.
| S-EPMC11830296 | biostudies-literature
Project description:Cell type-specific regulatory programs that drive type 1 diabetes (T1D) in the pancreas are poorly understood. Here we performed single nucleus multiomics and spatial transcriptomics in up to 32 non-diabetic (ND), autoantibody-positive (AAB+), and T1D pancreas donors. Genomic profiles from 853,005 cells mapped to 12 pancreatic cell types, including multiple exocrine sub-types. Beta, acinar, and other cell types, and related cellular niches, had altered abundance and gene activity in T1D progression, including distinct pathways altered in AAB+ compared to T1D. We identified epigenomic drivers of gene activity in T1D and AAB+ which, combined with genetic association, revealed causal pathways of T1D risk including antigen presentation in beta cells. Finally, single cell and spatial profiles together revealed widespread changes in cell-cell signaling in T1D including signals affecting beta cell regulation. Overall, these results revealed drivers of T1D progression in the pancreas, which form the basis for therapeutic targets for disease prevention.
| S-EPMC11870426 | biostudies-literature
Project description:Single-cell RNA sequencing methods can profile the transcriptomes of single cells but cannot preserve spatial information. Conversely, spatial transcriptomics assays can profile spatial regions in tissue sections, but do not have single-cell resolution. Here, we developed a computational method called CellTrek that combines these two datasets to achieve single-cell spatial mapping through coembedding and metric learning approaches. We benchmarked CellTrek using simulation and in situ hybridization datasets, which demonstrated its accuracy and robustness. We then applied CellTrek to existing mouse brain and kidney datasets and showed that CellTrek can detect topological patterns of different cell types and cell states. We performed single-cell RNA sequencing and spatial transcriptomics experiments on two ductal carcinoma in situ tissues and applied CellTrek to identify tumor subclones that were restricted to different ducts, and specific T cell states adjacent to the tumor areas. Our data show that CellTrek can accurately map single cells in diverse tissue types to resolve their spatial organization.
| S-EPMC9673606 | biostudies-literature
Project description:Blood glucose levels are tightly controlled by the coordinated action of at least four cell types constituting pancreatic islets. Changes in the proportion and/or function of these cells are associated with genetic and molecular pathophysiology of monogenic, type 1, and type 2 (T2D) diabetes. Cellular heterogeneity impedes precise understanding of the molecular components of each islet cell type that govern islet (dys)function, particularly the less abundant delta and gamma/pancreatic polypeptide (PP) cells. Here, we report single-cell transcriptomes for 638 cells from nondiabetic (ND) and T2D human islet samples. Analyses of ND single-cell transcriptomes identified distinct alpha, beta, delta, and PP/gamma cell-type signatures. Genes linked to rare and common forms of islet dysfunction and diabetes were expressed in the delta and PP/gamma cell types. Moreover, this study revealed that delta cells specifically express receptors that receive and coordinate systemic cues from the leptin, ghrelin, and dopamine signaling pathways implicating them as integrators of central and peripheral metabolic signals into the pancreatic islet. Finally, single-cell transcriptome profiling revealed genes differentially regulated between T2D and ND alpha, beta, and delta cells that were undetectable in paired whole islet analyses. This study thus identifies fundamental cell-type-specific features of pancreatic islet (dys)function and provides a critical resource for comprehensive understanding of islet biology and diabetes pathogenesis.
| S-EPMC5287227 | biostudies-literature
Project description:BackgroundSingle-cell sequencing approaches allow gene expression to be measured at the single-cell level, providing opportunities and challenges to study the aetiology of complex diseases, including cancer.MethodsBased on single-cell gene and lncRNA expression levels, we proposed a computational framework for cell type identification that fully considers cell dropout characteristics. First, we defined the dropout features of the cells and identified the dropout clusters. Second, we constructed a differential co-expression network and identified differential modules. Finally, we identified cell types based on the differential modules.ResultsThe method was applied to single-cell melanoma data, and eight cell types were identified. Enrichment analysis of the candidate cell marker genes for the two key cell types showed that both key cell types were closely related to the physiological activities of the major histocompatibility complex (MHC); one key cell type was associated with mitosis-related activities, and the other with pathways related to ten diseases.ConclusionsThrough identification and analysis of key melanoma-related cell types, we explored the molecular mechanism of melanoma, providing insight into melanoma research. Moreover, the candidate cell markers for the two key cell types are potential therapeutic targets for melanoma.
| S-EPMC8596920 | biostudies-literature
Project description:The mammalian nucleus is compartmentalized by diverse subnuclear structures. These subnuclear structures, marked by nuclear bodies and histone modifications, are often cell-type specific and affect gene regulation and 3D genome organization1-3. Understanding nuclear organization requires identifying the molecular constituents of subnuclear structures and mapping their associations with specific genomic loci in individual cells, within complex tissues. Here, we introduce two-layer DNA seqFISH+, which allows simultaneous mapping of 100,049 genomic loci, together with nascent transcriptome for 17,856 genes and a diverse set of immunofluorescently labeled subnuclear structures all in single cells in cell lines and adult mouse cerebellum. Using these multi-omics datasets, we showed that repressive chromatin compartments are more variable by cell type than active compartments. We also discovered a single exception to this rule: an RNA polymerase II (RNAPII)-enriched compartment was associated with long, cell-type specific genes (> 200kb), in a manner distinct from nuclear speckles. Further, our analysis revealed that cell-type specific facultative and constitutive heterochromatin compartments marked by H3K27me3 and H4K20me3 are enriched at specific genes and gene clusters, respectively, and shape radial chromosomal positioning and inter-chromosomal interactions in neurons and glial cells. Together, our results provide a single-cell high-resolution multi-omics view of subnuclear compartments, associated genomic loci, and their impacts on gene regulation, directly within complex tissues.
| S-EPMC10197539 | biostudies-literature
Project description:Atherosclerotic cardiovascular diseases are the major cause of death worldwide. CD4 T cells responding to Apolipoprotein B (ApoB), the core protein of most lipoproteins, have been identified as critical disease modulators. In healthy individuals, ApoB-reactive (ApoB+) CD4 T cells are mostly regulatory T cells (Tregs), which exert anti-inflammatory effects. Yet, they may obtain pro-inflammatory features and thus become proatherogenic. Evidence from animal studies suggests that vaccination against certain major histocompatibility complex (MHC) II-binding ApoB peptides induces an expansion of ApoB+ Tregs and thus confers atheroprotection. To date, in-depth phenotyping of vaccine-expanded ApoB+ T cells has not yet been performed. To this end, we vaccinated C57BL/6J mice with the ApoB-peptide P6 (ApoB978-993 TGAYSNASSTESASY) and performed single-cell RNA sequencing of tetramer-sorted P6+ T cells. P6+ cells were clonally expanded (one major, two minor clones) and formed a transcriptional cluster distinct from clusters mainly containing non-expanded P6+ and P6- cells. Transcriptomic profiling revealed that most expanded P6+ cells had a strong Treg signature and highly expressed genes mediating suppressive functions. Yet, some expanded P6+ cells only had a residual Treg signature and expressed genes related to T helper 1 (TH1) cells, which are proatherogenic. Modeling the T cell receptor (TCR) and P6:MHC-II interaction showed that only three amino acid residues in the α and β chain contact the P6 peptide in the MHC-II groove and thus determine the specificity of this TCR to P6. Our data begin to reveal the vaccination-induced response to an ApoB epitope.
| S-EPMC9849899 | biostudies-literature
Project description:Analysis of cell-type-specific transcriptomes is vital for understanding the biology of tissues and organs in the context of multicellular organisms. In this Protocol Extension, we combine a previously developed cell-type-specific metabolic RNA labeling method (thiouracil (TU) tagging) and a pipeline to detect the labeled transcripts by a novel RNA sequencing (RNA-seq) method, SLAMseq (thiol (SH)-linked alkylation for the metabolic sequencing of RNA). By injecting a uracil analog, 4-thiouracil, into transgenic mice that express cell-type-specific uracil phosphoribosyltransferase (UPRT), an enzyme required for 4-thiouracil incorporation into newly synthesized RNA, only cells expressing UPRT synthesize thiol-containing RNA. Total RNA isolated from a tissue of interest is then sequenced with SLAMseq, which introduces thymine to cytosine (T>C) conversions at the sites of the incorporated 4-thiouracil. The resulting sequencing reads are then mapped with the T>C-aware alignment software, SLAM-DUNK, which allows mapping of reads containing T>C mismatches. The number of T>C conversions per transcript is further analyzed to identify which transcripts are synthesized in the UPRT-expressing cells. Thus, our method, SLAM-ITseq (SLAMseq in tissue), enables cell-specific transcriptomics without laborious FACS-based cell sorting or biochemical isolation of the labeled transcripts used in TU tagging. In the murine tissues we assessed previously, this method identified ~5,000 genes that are expressed in a cell type of interest from the total RNA pool from the tissue. Any laboratory with access to a high-throughput sequencer and high-power computing can adapt this protocol with ease, and the entire pipeline can be completed in <5 d.
| S-EPMC6997029 | biostudies-literature
Project description:Genome-wide transcriptomic analyses in whole tissues reflect the aggregate gene expression in heterogeneous cell populations comprising resident and migratory cells, and they are unable to identify cell type-specific information. We used a computational method (population-specific expression analysis [PSEA]) to decompose gene expression in gingival tissues into cell type-specific signatures for 8 cell types (epithelial cells, fibroblasts, endothelial cells, neutrophils, monocytes/macrophages, plasma cells, T cells, and B cells). We used a gene expression data set generated using microarrays from 120 persons (310 tissue samples; 241 periodontitis affected and 69 healthy). Decomposition of the whole-tissue transcriptomes identified differentially expressed genes in each of the cell types, which mapped to biologically relevant pathways, including dysregulation of Th17 cell differentiation, AGE-RAGE signaling, and epithelial-mesenchymal transition in epithelial cells. We validated selected PSEA-predicted, differentially expressed genes in purified gingival epithelial cells and B cells from an unrelated cohort (n = 15 persons), each of whom contributed with 1 periodontitis-affected and 1 healthy gingival tissue sample. Differential expression of these genes by quantitative reverse transcription polymerase chain reaction corroborated the PSEA predictions and pointed to dysregulation of biologically important pathways in periodontitis. Collectively, our results demonstrate the robustness of the PSEA in the decomposition of gingival tissue transcriptomes and its ability to identify differentially regulated transcripts in particular cellular constituents. These genes may serve as candidates for further investigation with respect to their roles in the pathogenesis of periodontitis.
| S-EPMC8060774 | biostudies-literature