Title: | Gene expression-based subtype classification for high-grade serous ovarian cancer |
---|---|
Description: | This package implements four major subtype classifiers for high-grade serous (HGS) ovarian cancer as described by Helland et al. (PLoS One, 2011), Bentink et al. (PLoS One, 2012), Verhaak et al. (J Clin Invest, 2013), and Konecny et al. (J Natl Cancer Inst, 2014). In addition, the package implements a consensus classifier, which consolidates and improves on the robustness of the proposed subtype classifiers, thereby providing reliable stratification of patients with HGS ovarian tumors of clearly defined subtype. |
Authors: | Gregory M Chen [aut], Lavanya Kannan [aut], Ludwig Geistlinger [aut], Victor Kofia [aut], Levi Waldron [aut], Christopher Eeles [ctb], Benjamin Haibe-Kains [aut, cre] |
Maintainer: | Benjamin Haibe-Kains <[email protected]> |
License: | Artistic-2.0 |
Version: | 1.29.0 |
Built: | 2024-11-08 06:18:44 UTC |
Source: | https://github.com/bioc/consensusOV |
Merging all individual esets and merging them into a big eset
dataset.merging( esets, method = c("union", "intersect"), standardization = c("quantile", "robust.scaling", "scaling", "none"), nthread = 1 )
dataset.merging( esets, method = c("union", "intersect"), standardization = c("quantile", "robust.scaling", "scaling", "none"), nthread = 1 )
esets |
The list containing all GSE file that need to be merged. |
method |
either "unique" or "intersect" is use to for selecting geneid |
standardization |
choose between "quantile", "robust.scaling", "scaling" or "none" |
nthread |
number of threads (1 by default) |
The merging eset
Get ovarian cancer subtypes as defined by Bentink et al., 2012
get.bentink.subtypes(expression.matrix, entrez.ids)
get.bentink.subtypes(expression.matrix, entrez.ids)
expression.matrix |
A matrix of gene expression values with rows as genes, columns as samples. |
entrez.ids |
A vector of Entrez Gene IDs, corresponding to the rows of
|
A list with first value Bentink.subtypes
containing a factor
of subtype names; and second value angio
containing the ouput of
genefu::ovcAngiogenic
Bentink et al. Angiogenic mRNA and microRNA gene expression signature predicts a novel subtype of serous ovarian cancer. PloS one (2012).
library(Biobase) library(genefu) data(GSE14764.eset) expression.matrix <- exprs(GSE14764.eset) entrez.ids <- as.character(fData(GSE14764.eset)$EntrezGene.ID) get.bentink.subtypes(expression.matrix, entrez.ids)
library(Biobase) library(genefu) data(GSE14764.eset) expression.matrix <- exprs(GSE14764.eset) entrez.ids <- as.character(fData(GSE14764.eset)$EntrezGene.ID) get.bentink.subtypes(expression.matrix, entrez.ids)
Get consensusOV ovarian cancer subtypes
get.consensus.subtypes( expression.matrix, entrez.ids, concordant.tumors.only = TRUE, remove.using.cutoff = FALSE, percentage.dataset.removed = 0.75, .training.dataset = consensus.training.dataset.full, .dataset.names.to.keep = names(esets.rescaled.classified.filteredgenes) ) margin(rf.probs)
get.consensus.subtypes( expression.matrix, entrez.ids, concordant.tumors.only = TRUE, remove.using.cutoff = FALSE, percentage.dataset.removed = 0.75, .training.dataset = consensus.training.dataset.full, .dataset.names.to.keep = names(esets.rescaled.classified.filteredgenes) ) margin(rf.probs)
expression.matrix |
A matrix of gene expression values with rows as genes, columns as samples. |
entrez.ids |
A vector of Entrez Gene IDs, corresponding to the rows of
|
concordant.tumors.only |
Logical. Should the classifier trained only on tumors that are concordantly classified by Helland, Konecny, and Verhaak? Defaults to TRUE. |
remove.using.cutoff |
Specify whether to classify NA for samples that do not meet a margin cutoff |
percentage.dataset.removed |
If remove.using.cutoff is TRUE, then classify this percentage of samples to NA based on margin values |
.training.dataset |
ExpressionSet containing the training data. Defaults to the pooled dataset across selected MetaGxOvarian datasets. |
.dataset.names.to.keep |
Names of MetaGxOvarian datasets to use for training |
rf.probs |
random forest probabilities for each subtype as returned
by |
get.consensus.subtypes
returns a list with
first value consensusOV.subtypes
containing a
factor of subtype labels; and second value rf.probs
containing a matrix
of subtype probabilities.
margin
returns a numeric vector containing the classification margin
scores, i.e. the difference between the top two subtype scores for each tumor.
library(Biobase) data(GSE14764.eset) expression.matrix <- exprs(GSE14764.eset) entrez.ids <- as.character(fData(GSE14764.eset)$EntrezGene.ID) sts <- get.consensus.subtypes(expression.matrix, entrez.ids) margins <- margin(sts$rf.probs)
library(Biobase) data(GSE14764.eset) expression.matrix <- exprs(GSE14764.eset) entrez.ids <- as.character(fData(GSE14764.eset)$EntrezGene.ID) sts <- get.consensus.subtypes(expression.matrix, entrez.ids) margins <- margin(sts$rf.probs)
Get ovarian cancer subtypes as defined by Hao et al., 2017
get.hao.subtypes(expression.matrix, entrez.ids)
get.hao.subtypes(expression.matrix, entrez.ids)
expression.matrix |
A matrix of gene expression values with genes as rows, samples as columns. |
entrez.ids |
A vector of Entrez Gene IDs, corresponding to the rows of
|
Hao et al., 2017 derived a gene signature to predict the tissue of origin of ovarian tumors as either fallopian tube (FT) or ovarian surface epithelium (OSE).
The authors found that expression patterns of tissue-specific genes, prognostic genes, and molecular markers support a dualistic tissue origin of ovarian cancer, from either FT or OSE.
The subtype classifier considers 112 signature genes including 37 genes upregulated in FT and 75 genes upregulated in OSE. A score is computed that is designed to range from 0 to 1 for FT tumors, while OSE tumors have a score ranging from -1 to 0.
A list with first value tissue
containing a factor
of subtype names (tissue of origin); and second value score
containing
the tissue-of-origin score.
Ludwig Geistlinger
Hao et al. (2017) Integrated analysis reveals tubal- and ovarian-originated serous ovarian cancer and predicts differential therapeutic responses. Clinical Cancer Research, 23:7400-11.
library(Biobase) data(GSE14764.eset) expression.matrix <- exprs(GSE14764.eset) entrez.ids <- as.character(fData(GSE14764.eset)$EntrezGene.ID) get.hao.subtypes(expression.matrix, entrez.ids)
library(Biobase) data(GSE14764.eset) expression.matrix <- exprs(GSE14764.eset) entrez.ids <- as.character(fData(GSE14764.eset)$EntrezGene.ID) get.hao.subtypes(expression.matrix, entrez.ids)
Get ovarian cancer subtypes as defined by Helland et al., 2011
get.helland.subtypes(expression.matrix, entrez.ids)
get.helland.subtypes(expression.matrix, entrez.ids)
expression.matrix |
A matrix of gene expression values with rows as genes, columns as samples. |
entrez.ids |
A vector of Entrez Gene IDs, corresponding to the rows of
|
A list with first value Helland.subtypes
containing a factor
of subtype names; and second value subtype.scores
containing a matrix
of subtype scores
Helland et al. Deregulation of MYCN, LIN28B and LET7 in a molecular subtype of aggressive high-grade serous ovarian cancers. PloS one (2011).
library(Biobase) data(GSE14764.eset) expression.matrix <- exprs(GSE14764.eset) entrez.ids <- as.character(fData(GSE14764.eset)$EntrezGene.ID) get.helland.subtypes(expression.matrix, entrez.ids)
library(Biobase) data(GSE14764.eset) expression.matrix <- exprs(GSE14764.eset) entrez.ids <- as.character(fData(GSE14764.eset)$EntrezGene.ID) get.helland.subtypes(expression.matrix, entrez.ids)
Get ovarian cancer subtypes as defined by Konecny et al., 2014
get.konecny.subtypes(expression.matrix, entrez.ids)
get.konecny.subtypes(expression.matrix, entrez.ids)
expression.matrix |
A matrix of gene expression values with rows as genes, columns as samples. |
entrez.ids |
A vector of Entrez Gene IDs, corresponding to the rows of
|
A list with first value Konecny.subtypes
containing a factor
of subtype names; and second value spearman.cc.vals
containing the
Spearman correlation values per subtype
Konecny et al. Prognostic and therapeutic relevance of molecular subtypes in high-grade serous ovarian cancer. Journal of the National Cancer Institute (2014).
library(Biobase) data(GSE14764.eset) expression.matrix <- exprs(GSE14764.eset) entrez.ids <- as.character(fData(GSE14764.eset)$EntrezGene.ID) get.konecny.subtypes(expression.matrix, entrez.ids)
library(Biobase) data(GSE14764.eset) expression.matrix <- exprs(GSE14764.eset) entrez.ids <- as.character(fData(GSE14764.eset)$EntrezGene.ID) get.konecny.subtypes(expression.matrix, entrez.ids)
Get ovarian cancer subtypes
get.subtypes( expression.dataset, entrez.ids = NULL, method = c("consensusOV", "Helland", "Verhaak", "Konecny", "Bentink"), ... )
get.subtypes( expression.dataset, entrez.ids = NULL, method = c("consensusOV", "Helland", "Verhaak", "Konecny", "Bentink"), ... )
expression.dataset |
Either a matrix of gene expression values with rows
as genes, columns as samples;
or a BioBase::ExpressionSet object from MetaGxOvarian.
If |
entrez.ids |
A vector of Entrez Gene IDs, corresponding to the rows of
|
method |
The subtyping method to use |
... |
Optional parameters to be passed to the low level function |
A list with first value Konecny.subtypes
containing a factor
of subtype names; and second value spearman.cc.vals
containing the
Spearman correlation values per subtype
library(Biobase) data(GSE14764.eset) expression.matrix <- exprs(GSE14764.eset) entrez.ids <- as.character(fData(GSE14764.eset)$EntrezGene.ID) get.subtypes(expression.matrix, entrez.ids, method="Konecny")
library(Biobase) data(GSE14764.eset) expression.matrix <- exprs(GSE14764.eset) entrez.ids <- as.character(fData(GSE14764.eset)$EntrezGene.ID) get.subtypes(expression.matrix, entrez.ids, method="Konecny")
Get ovarian cancer subtypes as defined by Verhaak et al., 2013
get.verhaak.subtypes(expression.matrix, entrez.ids)
get.verhaak.subtypes(expression.matrix, entrez.ids)
expression.matrix |
A matrix of gene expression values with rows as genes, columns as samples. |
entrez.ids |
A vector of Entrez Gene IDs, corresponding to the rows of
|
A list with first value Verhaak.subtypes
containing a factor
of subtype names; and second value gsva
containing the GSVA subtype
scores
Verhaak et al. Prognostically relevant gene signatures of high-grade serous ovarian carcinoma. The Journal of Clinical Investigation (2013)
library(Biobase) data(GSE14764.eset) expression.matrix <- exprs(GSE14764.eset) entrez.ids <- as.character(fData(GSE14764.eset)$EntrezGene.ID) get.konecny.subtypes(expression.matrix, entrez.ids)
library(Biobase) data(GSE14764.eset) expression.matrix <- exprs(GSE14764.eset) entrez.ids <- as.character(fData(GSE14764.eset)$EntrezGene.ID) get.konecny.subtypes(expression.matrix, entrez.ids)
A Biobase::ExpressionSet from package MetaGxOvarian for the dataset GSE14764
GSE14764.eset
GSE14764.eset
A Biobase::ExpressionSet object
http://biorxiv.org/content/biorxiv/early/2016/05/12/052910.full.pdf