Dataset Information

De novo Gene Signature Identification from Single-Cell RNA-Seq with Hierarchical Poisson Factorization

ABSTRACT: Common approaches to gene signature discovery in single cell RNA-sequencing depend upon predefined structures like clustering or pseudo-temporal orderings, do not account for the sparsity of single cell data, or require prior normalization. We present single cell Hierarchical Poisson Factorization (scHPF), a Bayesian factorization method that adapts Hierarchical Poisson Factorization for de novo discovery of both continuous and discrete expression patterns in complex tissues. scHPF does not require prior normalization and outperforms other methods in benchmark datasets. Applied to single cell RNA-sequencing of the core and margin of a high-grade glioma, scHPF uncovers subtle regional expression biases within glioma subpopulations and an expression signature associated with inferior survival in glioblastoma.

ORGANISM(S): Homo sapiens

PROVIDER: GSE116621 | GEO | 2019/02/21

REPOSITORIES: GEO

ACCESS DATA

Similar Datasets

Project description:Microarray data analysis Intensity of 21,939 gene features per array was extracted from scanned microarray images using Feature Extraction 5.1.1 software (Agilent Technologies), which performs background subtractions and dye normalization. This normalization method is targeted at detecting changes in relative expression of individual genes rather than global expression. Global expression change would require external normalization controls (van de Peppel et al., 2003). Text output was processed using an application developed in-house to perform ANOVA analysis (http://lgsun.grc.nia.nih.gov/ANOVA/). Intensity of features measured with a >50% error were replaced with missing values except features with very low intensity. Surrogate values equal to mean error were inserted for values that were negative or less than the probe error. Data were analyzed using ANOVA with embryonic stage as a factor. The small number of biological replications typical in expression profiling experiments results in a highly variable error variance, and this problem is usually addressed by log-ratio thresholds (Schena et al., 1995) that require subjective decisions about biological significance, or by Bayesian adjustment of error variance (Baldi and Long, 2001), which may still underestimate error variance and result in false positive results. To reduce false-positives, we opted for a very conservative error model in which error variance that is used for estimating F-statistics is the maximum of the actual error variance for this gene and the average error variance in 500 genes with similar average intensity. Statistical significance was determined using the False Discovery Rate (FDR = 10%) method (Benjamini and Hochberg, 1995). Pair-wise mean comparison was done with t-statistics and FDR=10%. Further data processing including scatter plots, hierarchical clustering, and principal component analysis (PCA) were also performed through NIA microarray analysis tool (http://lgsun.grc.nia.nih.gov/ANOVA/). The input file for NIA microarray analysis tool is available at http://lgsun.grc.nia.nih.gov/microarray/data.html

Project description:Feature Extraction 5.1.1 software (Agilent Technologies) performs background subtractions and dye normalization. This normalization method is targeted at detecting changes in relative expression of individual genes rather than global expression. Global expression change would require external normalization controls (van de Peppel et al., 2003). Text output was processed using an application developed in-house to perform ANOVA analysis (http://lgsun.grc.nia.nih.gov/ANOVA/). Intensity of features measured with a >50% error were replaced with missing values except features with very low intensity. Surrogate values equal to mean error were inserted for values that were negative or less than the probe error. Data were analyzed using ANOVA with embryonic stage as a factor. The small number of biological replications typical in expression profiling experiments results in a highly variable error variance, and this problem is usually addressed by log-ratio thresholds (Schena et al., 1995) that require subjective decisions about biological significance, or by Bayesian adjustment of error variance (Baldi and Long, 2001), which may still underestimate error variance and result in false positive results. To reduce false-positives, we opted for a very conservative error model in which error variance that is used for estimating F-statistics is the maximum of the actual error variance for this gene and the average error variance in 500 genes with similar average intensity. Statistical significance was determined using the False Discovery Rate (FDR = 10%) method (Benjamini and Hochberg, 1995). Pair-wise mean comparison was done with t-statistics and FDR=10%. Further data processing including scatter plots, hierarchical clustering, and principal component analysis (PCA) were also performed through NIA microarray analysis tool (http://lgsun.grc.nia.nih.gov/ANOVA/). Details are described in the publication (Hamatani T. et al., Developmental Cell. Published online Dec. 18, 2003).

Project description:Objective. Microarray analysis was used to determine whether children with recent onset polyarticular juvenile idiopathic arthritis (JIA) exhibit biologically or clinically informative gene expression signatures in peripheral blood mononuclear cells (PBMC). Methods. Peripheral blood samples were obtained from 59 healthy children and 61 children with polyarticular JIA prior to treatment with second-line medications, such as methotrexate or biological agents. RNA was purified from Ficoll-isolated mononuclear cells, fluorescently labeled and then hybridized to Affymetrix U133 Plus 2.0 GeneChips. Data were analyzed using ANOVA at a 5% false discovery rate threshold after Robust Multi-Array Average pre-processing and Distance Weighted Discrimination normalization. Results. Initial analysis revealed 873 probe sets for genes that were differentially expressed between polyarticular JIA and controls. Hierarchical clustering of these probe sets distinguished three subgroups within polyarticular JIA. Prototypical subjects within each subgroup were identified and used to define subgroup-specific gene expression signatures. One of these signatures was associated with monocyte markers, another with transforming growth factor-beta-inducible genes, and a third with immediate-early genes. Correlation of these gene expression signatures with clinical and biological features of JIA subgroups suggests direct relevance to aspects of disease activity and supports the division of polyarticular JIA into distinct subsets. Conclusions. PBMC gene expression signatures in recent onset polyarticular JIA reflect discrete disease processes and offer a molecular classification of disease. Keywords: Patient vs. control, reassessment of phenotype PBMC samples were obtained from 59 healthy children and 61 children with polyarticular JIA prior to treatment with second-line medications, such as methotrexate or biological agents. RNA was purified from Ficoll-isolated mononuclear cells, fluorescently labeled and then hybridized to Affymetrix U133 Plus 2.0 GeneChips. Data were analyzed using ANOVA at a 5% false discovery rate threshold after Robust Multi-Array Average pre-processing and Distance Weighted Discrimination normalization.

Dataset Information

De novo Gene Signature Identification from Single-Cell RNA-Seq with Hierarchical Poisson Factorization

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure