Quantifying the effect of experimental perturbations in single-cell3RNA-sequencing data using graph signal processing
Ontology highlight
ABSTRACT:
ORGANISM(S): Homo sapiens
PROVIDER: GSE161465 | GEO | 2020/11/14
REPOSITORIES: GEO
Similar Datasets
Project description:Small cell lung cancer (SCLC) tumors comprise heterogeneous mixtures of cell states, categorized into neuroendocrine (NE) and non-neuroendocrine (non-NE) transcriptional subtypes. NE to non-NE state transitions, fueled by plasticity, likely underlie adaptability to treatment and dismal survival rates. Here, we apply an archetypal analysis to model plasticity by recasting SCLC phenotypic heterogeneity through multi-task evolutionary theory. Cell line and tumor transcriptomics data fit well in a five-dimensional convex polytope whose vertices optimize tasks reminiscent of pulmonary NE cells, the SCLC normal counterparts. These tasks, supported by knowledge and experimental data, include proliferation, slithering, metabolism, secretion, and injury repair, reflecting cancer hallmarks. SCLC subtypes, either at the population or single-cell level, can be positioned in archetypal space by bulk or single-cell transcriptomics, respectively, and characterized as task specialists or multi-task generalists by the distance from archetype vertex signatures. In the archetype space, modeling single-cell plasticity as a Markovian process along an underlying state manifold indicates that task trade-offs, in response to microenvironmental perturbations or treatment, may drive cell plasticity. Stifling phenotypic transitions and plasticity may provide new targets for much-needed translational advances in SCLC.
2022-08-22 | GSE211751 | GEO
Project description:A key challenge in single cell RNA-sequencing (scRNA-seq) data analysis are dataset- and batch-specific differences that can obscure the biological signal of interest. While there are various tools and methods to perform data integration and correct for batch effects, their performance can vary between datasets and according to the nature of the bias. Therefore, it is important to understand how batch effects manifest in order to adjust for them in a reliable way. Here, we systematically explore batch effects in scRNA-seq data from a variety of datasets according to magnitude, cell type specificity and complexity. We developed a cell-specific mixing score (\texttt{cms}) that quantifies how well cells from multiple batches are mixed. By considering distance distributions (in a lower dimensional space), the score is able to detect local batch bias and differentiate between unbalanced batches (i.e., when one cell type is more abundant in a batch) and systematic differences between cells of the same cell type. We implemented the \texttt{cms}, as well as related metrics to detect batch effects or measure structure preservation, in the CellMixS R/Bioconductor package. We systematically compare different metrics that have been proposed to quantify batch effects or bias in scRNA-seq data using real datasets with known batch effects and synthetic data that mimic various real data scenarios. While these metrics target the same question and are used interchangeably, we find differences in inter- and intra-dataset scalability, sensitivity and in a metric's ability to handle batch effects with differentially abundant cell types. We find that cell-specific metrics outperform cell type-specific and global metrics and recommend them for both method benchmarks and batch exploration.
2021-01-31 | E-MTAB-9916 | biostudies-arrayexpress
Project description:The paper describes a model of multiple myeloma.
Created by COPASI 4.26 (Build 213)
This model is described in the article:
A mathematical model of cell equilibrium and joint cell formation in multiple myeloma
M.A. Koenders, R. Saso
Journal of Theoretical Biology 390 (2016) 73–79
Abstract:
In Multiple Myeloma Bone Disease healthy bone remodelling is affected by tumour cells by means of paracrine cytokinetic signalling in such a way that osteoclast formation is enhanced and the growth of osteoblast cells inhibited. The participating cytokines are described in the literature. Osteoclast-induced myeloma cell growth is also reported. Based on existing mathematical models for healthy bone remo- delling a three-way equilibrium model is presented for osteoclasts, osteoblasts and myeloma cell populations to describe the progress of the illness in a scenario in which there is a secular increase in the cytokinetic interactive effectiveness of paracrine processes. The equilibrium state for the system is obtained. The paracrine interactive effectiveness is explored by parameter variation and the stable region in the parameter space is identified. Then recently-discovered joint myeloma–osteoclast cells are added to the model to describe the populations inside lytic lesions. It transpires that their presence expands the available parameter space for stable equilibrium, thus permitting a detrimental, larger population of osteoclasts and myeloma cells. A possible relapse mechanism for the illness is explored by letting joint cells dissociate. The mathematics then permits the evaluation of the evolution of the cell populations as a function of time during relapse.
To cite BioModels Database, please use: BioModels Database: An enhanced, curated and annotated resource for published quantitative kinetic models .
To the extent possible under law, all copyright and related or neighbouring rights to this encoded model have been dedicated to the public domain worldwide.
Please refer to CC0 Public Domain Dedication for more information.
2019-08-29 | BIOMD0000000804 | BioModels
Project description:10 pilosebaceous units (PUs) at the anagen phase were isolated from each involved (balding; vertex) and non-involved (control; occipital) scalp region using the follicular unit extraction (FUE) technique.
2020-11-30 | E-MTAB-9557 | biostudies-arrayexpress
Project description:Transcriptomes from multiple pre-meiotic stages of wild type, mac1, and msca1 maize anthers were characterized by microarray hybridization. The goal was to characterize the developmental progression as the anther specifies five cell types and grows rapidly precedeing meiotic entry. The stages characterized were immature anther primordia (0.15 mm long in maize) containing just stem cells, through somatic and germinal cell fate specification (0.20 and 0.25 mm), mitotic proliferation (0.4 mm), and finally the birth of the middle layer and tapetum (0.7 mm). To obtain cell-type specific markers, at 0.7 mm we also compared whole anthers to collections of laser-microdissected anther cell types including the archesporial cells (pre-meiotic germinal cells), nutritive layers (middle layer and tapetum) and structural layers (endothecium and epidemis) of the anther lobe. keyword: anther development, maize, male-sterile Three loop designs covered the early stages (up to 0.7 mm) with two replicates for each comparison. The first loop had 0.2 mm long anthers and compared wild type versus mac1 mutant versus msca1 mutant in a three vertex loop design. The second loop had four vertices and compared 0.15 mm WT anther primordia, 0.25 mm WT anthers, 0.4 mm WT anthers and finally 0.4 mm mac1 mutant anthers. The third had 0.7 mm anthers in a three vertex loop with the nutritive layers (middle layer and tapetum) at one vertex, the germinal pre-meiotic cells at another vertex, and whole anthers at a third vertex. The whole anther samples were also, separately and outside of the loop, compared in four replicates to the structural layers (endothecium and epidermis).
2013-05-01 | E-GEOD-43982 | biostudies-arrayexpress
Project description:In this study, we characterize transciprtional phenotypes of airway macrophagages (AMs) throughout homeostatsis, inflammation, and repair at single cell granularity. We confirm that cell origin is the major determinant of AM programing and describe two previously uncharacterized, transcriptionally distinct subdivisions of AMs based on proliferative capacity and inflammatory programing.
2019-04-07 | GSE120000 | GEO
Project description:The combinatorial complexity of histone samples is manifold higher than what is usually encountered in proteomics. Consequently, a considerably bigger part of the acquired MSMS spectra remains unannotated to date. Adapted search parameters can dig deeper into the dark histone ion space, but the lack of false discovery rate (FDR) control and the high level of ambiguity when searching combinatorial PTMs makes it very hard to assess whether the newly assigned ions are informative. Therefore, we use an easily adoptable time-lapse enzymatic deacetylation (HDAC1) of a commercial histone extract as a quantify-first strategy that allows isolating ion populations of interest, when studying e.g. acetylation on histones, that currently remain in the dark. By adapting search parameters to study potential issues in sample preparation, data acquisition and data analysis, we stepwise managed to double the portion of annotated precursors of interest from 10.5% to 21.6%. This strategy is intended to make up for the lack of validated FDR control and has led to several adaptations of our current workflow that will reduce the portion of the dark histone ion space in the future.
2022-10-04 | PXD009910 | Pride
Project description:To describe the global profile of genes regulated by miR-31, we compared the transcriptomics between miR-31 silenced and antagomir control treated naive CD4+ T cells. Interestingly, the cells could be separated into two subsets according FCS and SSC by flow cytometer after the silence of miR-31, one (~65%) was small sized cells with low granularity, the other (~30%) was large cells with high granularity, whereas the antagoNC treated naïve CD4+ T cells only contain a marginal fraction of large cells. We pooled “antagoNC”, “antagomiR31S” and “antagomiR31L” naïve CD4+ T cells from 5 healthy donors, respectively, and analyzed the differential genes among those three subsets.
2022-12-24 | GSE83853 | GEO
Project description:The aim of this pilot study was to combine standard clinical diagnostics with RNA and protein analyses from non-invasively obtained material. In 6 male volunteers with AGA, Norwood-Hamilton stage IIIv-IV, we performed global photography, phototrichogram, pH and sebum measurements, as well as skin surface cytokine measurements on AGA-affected frontal and vertex scalp skin compared to clinically unaffected occiput.
2017-06-30 | GSE78722 | GEO
Project description:To date, most epigenetic datasets have been generated on DNA samples isolated from bulk tissues. As the proportion of individual cell types within a sample can vary across individuals, systematic differences in cellular proportions that correlate with the phenotype of interest may manifest as differences in the overall epigenetic profile. Deconvolution algorithms calculate a series of continuous variables reflecting the underlying cellular heterogeneity of each sample from the bulk tissue profile that can be used to adjust for confounding in epigenome-wide association analyses. Here we provide novel reference profiles for brain cell types for use with reference based deconvolution algorithms. We used a FANS protocol recently described by our group to purify nuclei populations from prefrontal cortex tissue from 43 adult donors. Our initial gating strategy used an antibody against NeuN to isolate neuronal nuclei in combination with an antibody against SOX10 to separate oligodendrocyte nuclei from other glial nuclei. Subsequently, in a second gating strategy we additionally included an antibody against IRF8 to enrich microglia from the NeuNNeg/SOX10Neg fraction. Our third gating strategy used an antibody against SATB2 in place of NeuN to isolate excitatory neurons. We generated DNA methylation profiles using the Illumina EPIC array for NeuNPos (neuron-enriched; n = 28), NeuNNeg/SOX10Pos (oligodendrocyte-enriched; n = 24), NeuNNeg/SOX10Neg (microglia- and astrocyte-enriched; n = 21), NeuNNeg/SOX10Neg/IRF8Pos (microglia-enriched; n = 17), NeuNNeg/SOX10Neg/IRF8Neg, (astrocyte-enriched; n = 7), SATB2Pos (excitatory neuron-enriched; n = 9), and SATB2Neg (inhibitory neuron- and glial- enriched; n = 6) nuclei populations. We have demonstrated that these are applicable for use with established deconvolution algorithms to quantify the cellular heterogeneity of the cortex and other regions of the human brain from bulk DNAm data. These variables will be critical covariates to include in future epigenetic studies of brain disorders to minimise the risk of false positive associations and improve our understanding of the changes in the brain that underpin the development of psychiatric disorders and neurodegenerative diseases.
2023-11-05 | GSE234520 | GEO