Cross-platform comparison and visualisation of gene expression data using co-inertia analysis

BACKGROUND: Rapid development of DNA microarray technology has resulted in different laboratories adopting numerous different protocols and technological platforms, which has severely impacted on the comparability of array data. Current cross-platform comparison of microarray gene expression data are usually based on cross-referencing the annotation of each gene transcript represented on the arrays, extracting a list of genes common to all arrays and comparing expression data of this gene subset. Unfortunately, filtering of genes to a subset represented across all arrays often excludes many thousands of genes, because different subsets of genes from the genome are represented on different arrays. We wish to describe the application of a powerful yet simple method for cross-platform comparison of gene expression data. Co-inertia analysis (CIA) is a multivariate method that identifies trends or co-relationships in multiple datasets which contain the same samples. CIA simultaneously finds ordinations (dimension reduction diagrams) from the datasets that are most similar. It does this by finding successive axes from the two datasets with maximum covariance. CIA can be applied to datasets where the number of variables (genes) far exceeds the number of samples (arrays) such is the case with microarray analyses. RESULTS: We illustrate the power of CIA for cross-platform analysis of gene expression data by using it to identify the main common relationships in expression profiles on a panel of 60 tumour cell lines from the National Cancer Institute (NCI) which have been subjected to microarray studies using both Affymetrix and spotted cDNA array technology. The co-ordinates of the CIA projections of the cell lines from each dataset are graphed in a bi-plot and are connected by a line, the length of which indicates the divergence between the two datasets. Thus, CIA provides graphical representation of consensus and divergence between the gene expression profiles from different microarray platforms. Secondly, the genes that define the main trends in the analysis can be easily identified. CONCLUSIONS: CIA is a robust, efficient approach to coupling of gene expression datasets. CIA provides simple graphical representations of the results making it a particularly attractive method for the identification of relationships between large datasets

iBBiG: iterative binary bi-clustering of gene sets

Motivation: Meta-analysis of genomics data seeks to identify genes associated with a biological phenotype across multiple datasets; however, merging data from different platforms by their features (genes) is challenging. Meta-analysis using functionally or biologically characterized gene sets simplifies data integration is biologically intuitive and is seen as having great potential, but is an emerging field with few established statistical methods. Results: We transform gene expression profiles into binary gene set profiles by discretizing results of gene set enrichment analyses and apply a new iterative bi-clustering algorithm (iBBiG) to identify groups of gene sets that are coordinately associated with groups of phenotypes across multiple studies. iBBiG is optimized for meta-analysis of large numbers of diverse genomics data that may have unmatched samples. It does not require prior knowledge of the number or size of clusters. When applied to simulated data, it outperforms commonly used clustering methods, discovers overlapping clusters of diverse sizes and is robust in the presence of noise. We apply it to meta-analysis of breast cancer studies, where iBBiG extracted novel gene set—phenotype association that predicted tumor metastases within tumor subtypes

A multivariate approach to the integration of multi-omics datasets

Background: To leverage the potential of multi-omics studies, exploratory data analysis methods that provide systematic integration and comparison of multiple layers of omics information are required. We describe multiple co-inertia analysis (MCIA), an exploratory data analysis method that identifies co-relationships between multiple high dimensional datasets. Based on a covariance optimization criterion, MCIA simultaneously projects several datasets into the same dimensional space, transforming diverse sets of features onto the same scale, to extract the most variant from each dataset and facilitate biological interpretation and pathway analysis. Results: We demonstrate integration of multiple layers of information using MCIA, applied to two typical “omics” research scenarios. The integration of transcriptome and proteome profiles of cells in the NCI-60 cancer cell line panel revealed distinct, complementary features, which together increased the coverage and power of pathway analysis. Our analysis highlighted the importance of the leukemia extravasation signaling pathway in leukemia that was not highly ranked in the analysis of any individual dataset. Secondly, we compared transcriptome profiles of high grade serous ovarian tumors that were obtained, on two different microarray platforms and next generation RNA-sequencing, to identify the most informative platform and extract robust biomarkers of molecular subtypes. We discovered that the variance of RNA-sequencing data processed using RPKM had greater variance than that with MapSplice and RSEM. We provided novel markers highly associated to tumor molecular subtype combined from four data platforms. MCIA is implemented and available in the R/Bioconductor “omicade4” package. Conclusion: We believe MCIA is an attractive method for data integration and visualization of several datasets of multi-omics features observed on the same set of individuals. The method is not dependent on feature annotation, and thus it can extract important features even when there are not present across all datasets. MCIA provides simple graphical representations for the identification of relationships between large datasets

Stem Cell-Like Gene Expression in Ovarian Cancer Predicts Type II Subtype and Prognosis

Although ovarian cancer is often initially chemotherapy-sensitive, the vast majority of tumors eventually relapse and patients die of increasingly aggressive disease. Cancer stem cells are believed to have properties that allow them to survive therapy and may drive recurrent tumor growth. Cancer stem cells or cancer-initiating cells are a rare cell population and difficult to isolate experimentally. Genes that are expressed by stem cells may characterize a subset of less differentiated tumors and aid in prognostic classification of ovarian cancer. The purpose of this study was the genomic identification and characterization of a subtype of ovarian cancer that has stem cell-like gene expression. Using human and mouse gene signatures of embryonic, adult, or cancer stem cells, we performed an unsupervised bipartition class discovery on expression profiles from 145 serous ovarian tumors to identify a stem-like and more differentiated subgroup. Subtypes were reproducible and were further characterized in four independent, heterogeneous ovarian cancer datasets. We identified a stem-like subtype characterized by a 51-gene signature, which is significantly enriched in tumors with properties of Type II ovarian cancer; high grade, serous tumors, and poor survival. Conversely, the differentiated tumors share properties with Type I, including lower grade and mixed histological subtypes. The stem cell-like signature was prognostic within high-stage serous ovarian cancer, classifying a small subset of high-stage tumors with better prognosis, in the differentiated subtype. In multivariate models that adjusted for common clinical factors (including grade, stage, age), the subtype classification was still a significant predictor of relapse. The prognostic stem-like gene signature yields new insights into prognostic differences in ovarian cancer, provides a genomic context for defining Type I/II subtypes, and potential gene targets which following further validation may be valuable in the clinical management or treatment of ovarian cancer

GeneSigDB—a curated database of gene expression signatures

The primary objective of most gene expression studies is the identification of one or more gene signatures; lists of genes whose transcriptional levels are uniquely associated with a specific biological phenotype. Whilst thousands of experimentally derived gene signatures are published, their potential value to the community is limited by their computational inaccessibility. Gene signatures are embedded in published article figures, tables or in supplementary materials, and are frequently presented using non-standard gene or probeset nomenclature. We present GeneSigDB (http://compbio.dfci.harvard.edu/genesigdb) a manually curated database of gene expression signatures. GeneSigDB release 1.0 focuses on cancer and stem cells gene signatures and was constructed from more than 850 publications from which we manually transcribed 575 gene signatures. Most gene signatures (n = 560) were successfully mapped to the genome to extract standardized lists of EnsEMBL gene identifiers. GeneSigDB provides the original gene signature, the standardized gene list and a fully traceable gene mapping history for each gene from the original transcribed data table through to the standardized list of genes. The GeneSigDB web portal is easy to search, allows users to compare their own gene list to those in the database, and download gene signatures in most common gene identifier formats

CENP-F expression is associated with poor prognosis and chromosomal instability in patients with primary breast cancer

DNA microarrays have the potential to classify tumors according to their transcriptome. Tissue microarrays (TMAs) facilitate the validation of biomarkers by offering a high-throughput approach to sample analysis. We reanalyzed a high profile breast cancer DNA microarray dataset containing 96 tumor samples using a powerful statistical approach, between group analyses. Among the genes we identified was centromere protein-F (CENP-F), a gene associated with poor prognosis. In a published follow-up breast cancer DNA microarray study, comprising 295 tumour samples, we found that CENP-F upregulation was significantly associated with worse overall survival (p < 0.001) and reduced metastasis-free survival (p < 0.001). To validate and expand upon these findings, we used 2 independent breast cancer patient cohorts represented on TMAs. CENP-F protein expression was evaluated by immunohistochemistry in 91 primary breast cancer samples from cohort I and 289 samples from cohort II. CENP-F correlated with markers of aggressive tumor behavior including ER negativity and high tumor grade. In cohort I, CENP-F was significantly associated with markers of CIN including cyclin E, increased telomerase activity, c-Myc amplification and aneuploidy. In cohort II, CENP-F correlated with VEGFR2, phosphorylated Ets-2 and Ki67, and in multivariate analysis, was an independent predictor of worse breast cancer-specific survival (p = 0.036) and overall survival (p = 0.040). In conclusion, we identified CENP-F as a biomarker associated with poor outcome in breast cancer and showed several novel associations of biological significance

The Stem Cell Discovery Engine: an integrated repository and analysis system for cancer stem cell comparisons

Mounting evidence suggests that malignant tumors are initiated and maintained by a subpopulation of cancerous cells with biological properties similar to those of normal stem cells. However, descriptions of stem-like gene and pathway signatures in cancers are inconsistent across experimental systems. Driven by a need to improve our understanding of molecular processes that are common and unique across cancer stem cells (CSCs), we have developed the Stem Cell Discovery Engine (SCDE)—an online database of curated CSC experiments coupled to the Galaxy analytical framework. The SCDE allows users to consistently describe, share and compare CSC data at the gene and pathway level. Our initial focus has been on carefully curating tissue and cancer stem cell-related experiments from blood, intestine and brain to create a high quality resource containing 53 public studies and 1098 assays. The experimental information is captured and stored in the multi-omics Investigation/Study/Assay (ISA-Tab) format and can be queried in the data repository. A linked Galaxy framework provides a comprehensive, flexible environment populated with novel tools for gene list comparisons against molecular signatures in GeneSigDB and MSigDB, curated experiments in the SCDE and pathways in WikiPathways. The SCDE is available at http://discovery.hsci.harvard.edu

The CIN4 chromosomal instability qPCR classifier defines tumor aneuploidy and stratifies outcome in grade 2 breast cancer.

Purpose: Quantifying chromosomal instability (CIN) has both prognostic and predictive clinical utility in breast cancer. In order to establish a robust and clinically applicable gene expression-based measure of CIN, we assessed the ability of four qPCR quantified genes selected from the 70-gene Chromosomal Instability (CIN70) expression signature to stratify outcome in patients with grade 2 breast cancer. Methods: AURKA, FOXM1, TOP2A and TPX2 (CIN4), were selected from the CIN70 signature due to their high level of correlation with histological grade and mean CIN70 signature expression in silico. We assessed the ability of CIN4 to stratify outcome in an independent cohort of patients diagnosed between 1999 and 2002. 185 formalin-fixed, paraffin-embedded (FFPE) samples were included in the qPCR measurement of CIN4 expression. In parallel, ploidy status of tumors was assessed by flow cytometry. We investigated whether the categorical CIN4 score derived from the CIN4 signature was correlated with recurrence-free survival (RFS) and ploidy status in this cohort. Results: We observed a significant association of tumor proliferation, defined by Ki67 and mitotic index (MI), with both CIN4 expression and aneuploidy. The CIN4 score stratified grade 2 carcinomas into good and poor prognostic cohorts (mean RFS: 83.864.9 and 69.4 +- 8.2 months, respectively, p = 0.016) and its predictive power was confirmed by multivariate analysis outperforming MI and Ki67 expression. Conclusions: The first clinically applicable qPCR derived measure of tumor aneuploidy from FFPE tissue, stratifies grade 2 tumors into good and poor prognosis groups

TGFBR1 Intralocus Epistatic Interaction as a Risk Factor for Colorectal Cancer

In colorectal cancer (CRC), an inherited susceptibility risk affects about 35% of patients, whereas high-penetrance germline mutations account for <6% of cases. A considerable proportion of sporadic tumors could be explained by the coinheritance of multiple low-penetrance variants, some of which are common. We assessed the susceptibility to CRC conferred by genetic variants at the TGFBR1 locus. We analyzed 14 polymorphisms and the allele-specific expression (ASE) of TGFBR1 in 1025 individuals from the Spanish population. A case-control study was undertaken with 504 controls and 521 patients with sporadic CRC. Fourteen polymorphisms located at the TGFBR1 locus were genotyped with the iPLEX Gold (MassARRAY-Sequenom) technology. Descriptive analyses of the polymorphisms and haplotypes and association studies were performed with the SNPator workpackage. No relevant associations were detected between individual polymorphisms or haplotypes and the risk of CRC. The TGFBR1*9A/6A polymorphism was used for the ASE analysis. Heterozygous individuals were analyzed for ASE by fragment analysis using cDNA from normal tissue. The relative level of allelic expression was extrapolated from a standard curve. The cutoff value was calculated with Youden's index. ASE was found in 25.4% of patients and 16.4% of controls. Considering both bimodal and continuous types of distribution, no significant differences between the ASE values of patients and controls were identified. Interestingly, a combined analysis of the polymorphisms and ASE for the association with CRC occurrence revealed that ASE-positive individuals carrying one of the most common haplotypes (H2: 20.7%) showed remarkable susceptibility to CRC (RR: 5.25; 95% CI: 2.547–5.250; p<0.001) with a synergy factor of 3.7. In our study, 54.1% of sporadic CRC cases were attributable to the coinheritance of the H2 haplotype and TGFBR1 ASE. These results support the hypothesis that the allelic architecture of cancer genes, rather than individual polymorphisms, more accurately defines the CRC risk

Characteristics and outcomes of over 300,000 patients with COVID-19 and history of cancer in the United States and Spain

Background: We described the demographics, cancer subtypes, comorbidities, and outcomes of patients with a history of cancer and coronavirus disease 2019 (COVID-19). Second, we compared patients hospitalized with COVID-19 to patients diagnosed with COVID-19 and patients hospitalized with influenza. Methods: We conducted a cohort study using eight routinely collected health care databases from Spain and the United States, standardized to the Observational Medical Outcome Partnership common data model. Three cohorts of patients with a history of cancer were included: (i) diagnosed with COVID-19, (ii) hospitalized with COVID-19, and (iii) hospitalized with influenza in 2017 to 2018. Patients were followed from index date to 30 days or death. We reported demographics, cancer subtypes, comorbidities, and 30-day outcomes. Results: We included 366,050 and 119,597 patients diagnosed and hospitalized with COVID-19, respectively. Prostate and breast cancers were the most frequent cancers (range: 5%–18% and 1%–14% in the diagnosed cohort, respectively). Hematologic malignancies were also frequent, with non-Hodgkin’s lymphoma being among the five most common cancer subtypes in the diagnosed cohort. Overall, patients were aged above 65 years and had multiple comorbidities. Occurrence of death ranged from 2% to 14% and from 6% to 26% in the diagnosed and hospitalized COVID-19 cohorts, respectively. Patients hospitalized with influenza (n ¼ 67,743) had a similar distribution of cancer subtypes, sex, age, and comorbidities but lower occurrence of adverse events. Conclusions: Patients with a history of cancer and COVID-19 had multiple comorbidities and a high occurrence of COVID-19-related events. Hematologic malignancies were frequent. Impact: This study provides epidemiologic characteristics that can inform clinical care and etiologic studies.</p