Gene Set-Based Integrative Analysis Revealing Two Distinct Functional Regulation Patterns in Four Common Subtypes of Epithelial Ovarian Cancer

Systemic Inflammatory Patterns in Ovarian Cancer Patients: Analysis of Cytokines, Chemokines, and Microparticles

OBJECTIVE

 To compare the patterns of systemic inflammatory response in women with epithelial ovarian cancer (EOC) or no evidence of malignant disease, as well as to evaluate the profile of systemic inflammatory responses in type-1 and type-2 tumors. This is a non-invasive and indirect way to assess both tumor activity and the role of the inflammatory pattern during pro- and antitumor responses.

MATERIALS AND METHODS

 We performed a prospective evaluation of 56 patients: 30 women without evidence of malignant disease and 26 women with EOC. The plasma quantification of cytokines, chemokines, and microparticles (MPs) was performed using flow cytometry.

RESULTS

 Plasma levels of proinflammatory cytokines interleukin-12 (IL12), interleukin-6 (IL-6), tumor necrosis factor alpha (TNF-α) interleukin-1 beta (IL-1β), and interleukin-10 (IL-10), and C-X-C motif chemokine ligand 9 (CXCL-9) and C-X-C motif chemokine ligand 10 (CXCL-10) were significantly higher in patients with EOC than in those in the control group. Plasma levels of cytokine interleukin-17A (IL-17A) and MPs derived from endothelial cells were lower in patients with EOC than in the control group. The frequency of leukocytes and MPs derived from endothelial cells was higher in type-2 tumors than in those without malignancy. We observed an expressive number of inflammatory/regulatory cytokines and chemokines in the cases of EOC, as well as negative and positive correlations involving them, which leads to a higher complexity of these networks.

CONCLUSION

 The present study showed that, through the development of networks consisting of cytokines, chemokines, and MPs, there is a greater systemic inflammatory response in patients with EOC and a more complex correlation of these biomarkers in type-2 tumors.

The role of the inflammasome and its related pathways in ovarian cancer

Ovarian cancer (OC) is the most lethal tumor of the female reproductive tract and one of the most prevalent causes of death among female cancer patients. The absence of suitable procedures for early diagnosis, chemoresistance, and limited surgical debulking are all contributing to poor survival in patients. Despite aggressive treatments, the majority of patients have a recurrence within 16-22 months. Inflammasomes are multimeric protein complexes that play a major role in the innate immune system and inflammation. The overexpression of inflammasome-related pathways, including NOD-, LRR- and pyrin domain-containing protein 3 (NLRP3), Absent in melanoma 2 (AIM2), caspase-1, and Interleukin (IL)-1 have been reported in OC patients and in vitro cell lines. Therefore, inflammasome-related genes and protein might have a role in OC pathogenesis. Considering the potential relationship between inflammasome and OC, this study aimed to provide a literature-based review to explain the role of inflammasome and inflammation in cancer progression in OC.

Synergistic AHR Binding Pathway with EMT Effects on Serous Ovarian Tumors Recognized by Multidisciplinary Integrated Analysis

Epithelial ovarian cancers (EOCs) are fatal and obstinate among gynecological malignancies in advanced stage or relapsed status, with serous carcinomas accounting for the vast majority. Unlike EOCs, borderline ovarian tumors (BOTs), including serous BOTs, maintain a semimalignant appearance. Using gene ontology (GO)-based integrative analysis, we analyzed gene set databases of serous BOTs and serous ovarian carcinomas for dysregulated GO terms and pathways and identified multiple differentially expressed genes (DEGs) in various aspects. The SRC (SRC proto-oncogene, non-receptor tyrosine kinase) gene and dysfunctional aryl hydrocarbon receptor (AHR) binding pathway consistently influenced progression-free survival and overall survival, and immunohistochemical staining revealed elevated expression of related biomarkers (SRC, ARNT, and TBP) in serous BOT and ovarian carcinoma samples. Epithelial-mesenchymal transition (EMT) is important during tumorigenesis, and we confirmed the SNAI2 (Snail family transcriptional repressor 2, SLUG) gene showing significantly high performance by immunohistochemistry. During serous ovarian tumor formation, activated AHR in the cytoplasm could cooperate with SRC, enter cell nuclei, bind to AHR nuclear translocator (ARNT) together with TATA-Box Binding Protein (TBP), and act on DNA to initiate AHR-responsive genes to cause tumor or cancer initiation. Additionally, SNAI2 in the tumor microenvironment can facilitate EMT accompanied by tumorigenesis. Although it has not been possible to classify serous BOTs and serous ovarian carcinomas as the same EOC subtype, the key determinants of relevant DEGs (SRC, ARNT, TBP, and SNAI2) found here had a crucial role in the pathogenetic mechanism of both tumor types, implying gradual evolutionary tendencies from serous BOTs to ovarian carcinomas. In the future, targeted therapy could focus on these revealed targets together with precise detection to improve therapeutic effects and patient survival rates.

Dysregulated Immunological Functionome and Dysfunctional Metabolic Pathway Recognized for the Pathogenesis of Borderline Ovarian Tumors by Integrative Polygenic Analytics

The pathogenesis and molecular mechanisms of ovarian low malignant potential (LMP) tumors or borderline ovarian tumors (BOTs) have not been fully elucidated to date. Surgery remains the cornerstone of treatment for this disease, and diagnosis is mainly made by histopathology to date. However, there is no integrated analysis investigating the tumorigenesis of BOTs with open experimental data. Therefore, we first utilized a functionome-based speculative model from the aggregated obtainable datasets to explore the expression profiling data among all BOTs and two major subtypes of BOTs, serous BOTs (SBOTs) and mucinous BOTs (MBOTs), by analyzing the functional regularity patterns and clustering the separate gene sets. We next prospected and assembled the association between these targeted biomolecular functions and their related genes. Our research found that BOTs can be accurately recognized by gene expression profiles by means of integrative polygenic analytics among all BOTs, SBOTs, and MBOTs; the results exhibited the top 41 common dysregulated biomolecular functions, which were sorted into four major categories: immune and inflammatory response-related functions, cell membrane- and transporter-related functions, cell cycle- and signaling-related functions, and cell metabolism-related functions, which were the key elements involved in its pathogenesis. In contrast to previous research, we identified 19 representative genes from the above classified categories (IL6, CCR2 for immune and inflammatory response-related functions; IFNG, ATP1B1, GAS6, and PSEN1 for cell membrane- and transporter-related functions; CTNNB1, GATA3, and IL1B for cell cycle- and signaling-related functions; and AKT1, SIRT1, IL4, PDGFB, MAPK3, SRC, TWIST1, TGFB1, ADIPOQ, and PPARGC1A for cell metabolism-related functions) that were relevant in the cause and development of BOTs. We also noticed that a dysfunctional pathway of galactose catabolism had taken place among all BOTs, SBOTs, and MBOTs from the analyzed gene set databases of canonical pathways. With the help of immunostaining, we verified significantly higher performance of interleukin 6 (IL6) and galactose-1-phosphate uridylyltransferase (GALT) among BOTs than the controls. In conclusion, a bioinformatic platform of gene-set integrative molecular functionomes and biophysiological pathways was constructed in this study to interpret the complicated pathogenic pathways of BOTs, and these important findings demonstrated the dysregulated immunological functionome and dysfunctional metabolic pathway as potential roles during the tumorigenesis of BOTs and may be helpful for the diagnosis and therapy of BOTs in the future.

Immunology and ovarian cancers

The current standard therapy of epithelial ovarian cancer (EOC) is the combination of surgery and multiagent chemotherapy with/without adding targeted therapy. After treatment, response rate is high and nearly all patients can achieve complete remission, even though they are advanced diseases; however, the majority of patients will relapse and subsequently die of diseases within several years after initial treatment. When treatment options are limited, there is the urgent need for new novel therapeutic approaches for precise cancer control. The development of chemoresistance and evading of the anticancer immune response may be one of the important causes contributing to the therapeutic failure, and therefore, it represents a paradigm shift in cancer research. An individual's immune response and interaction with EOC cells might be one of the key factors for cancer treatment. There are many interventions, including targeting certain type immunogenic EOC-associated antigens, immune checkpoint blockade, and adoptive cellular therapy, which present a profound opportunity to revolutionize EOC treatment. This review will encompass the interaction between EOC and immune system and highlight recent data regarding the research of immunotherapy in EOC.

The Potential Role of Complement System in the Progression of Ovarian Clear Cell Carcinoma Inferred from the Gene Ontology-Based Immunofunctionome Analysis

Ovarian clear cell carcinoma (OCCC) is the second most common epithelial ovarian carcinoma (EOC). It is refractory to chemotherapy with a worse prognosis after the preliminary optimal debulking operation, such that the treatment of OCCC remains a challenge. OCCC is believed to evolve from endometriosis, a chronic immune/inflammation-related disease, so that immunotherapy may be a potential alternative treatment. Here, gene set-based analysis was used to investigate the immunofunctionomes of OCCC in early and advanced stages. Quantified biological functions defined by 5917 Gene Ontology (GO) terms downloaded from the Gene Expression Omnibus (GEO) database were used. DNA microarray gene expression profiles were used to convert 85 OCCCs and 136 normal controls into to the functionome. Relevant offspring were as extracted and the immunofunctionomes were rebuilt at different stages by machine learning. Several dysregulated pathogenic functions were found to coexist in the immunopathogenesis of early and advanced OCCC, wherein the complement-activation-alternative-pathway may be the headmost dysfunctional immunological pathway in duality for carcinogenesis at all OCCC stages. Several immunological genes involved in the complement system had dual influences on patients' survival, and immunohistochemistrical analysis implied the higher expression of C3a receptor (C3aR) and C5a receptor (C5aR) levels in OCCC than in controls.

The recent progress and therapy in endometriosis-associated ovarian cancer

Endometriosis-associated ovarian cancers (EAOCs) including endometrioid and clear cell ovarian carcinoma are subgroups of epithelial ovarian carcinomas (EOCs), which is generally acknowledged as the most lethal gynecological malignancy. Endometriosis (ES), a common clinical disease among women, presents with clinical symptoms of pelvic pain, infertility, or adnexal masses with the formation of endometrioma. It has long been considered to be a potential risk factor for developing EOCs, mainly of endometrioid and clear cell subtypes. Here, we compiled data from previous researches on deregulated molecular functions among ES and EOCs using gene set-based integrative analysis to decipher molecular and genetic relationships between ovarian ES and EOCs, especially EAOCs. We conclude that epidermal growth factor receptor (ERBB) and Phosphoinositide 3-kinases (PI3K)-related pathways are important in the carcinogenesis of type I EOCs, including clear cell, endometrioid, and mucinous ovarian carcinoma. Dysfunctional molecular pathways, such as deregulated oxidoreductase activity, metabolism, hormone activity, inflammatory response, innate immune response, and cell-cell signaling, played key roles in the malignant transformation of EAOCs. Nine genes related to inflammasome complex and inflammasome-related pathway were identified, indicating the importance of inflammation/immunity in EAOC transformation. We also collect progressive treatments of EAOC focused on targeted therapies and immunotherapy so far. This summarized information can contribute toward effective detection and treatment of EAOCs in the future.

Key Immunological Functions Involved in the Progression of Epithelial Ovarian Serous Carcinoma Discovered by the Gene Ontology-Based Immunofunctionome Analysis

Serous carcinoma (SC) is the most common and lethal subtype of epithelial ovarian carcinoma; immunotherapy is a potential treatment for SC, however, the global immunological functions of SC as well as their change during the progression of SC have not been investigated in detail till now. We conducted a genome-wide integrative analysis to investigate the immunofunctionomes of SC at four tumor stages by quantifying the immunological functions defined by the Gene Ontology gene sets. DNA microarray gene expression profiles of 1100 SCs and 136 normal ovarian tissue controls were downloaded from the Gene Expression Omnibus database and converted to the functionome. Then the immunofunctionomes were reconstructed by extracting the offspring from the functionome for the four SC staging groups. The key immunological functions extracted from immunofunctionomes with a series of filters revealed that the immunopathy of SC consisted of a group of deregulated functions with the core members including B cell activation and differentiation, regulation of leukocyte chemotaxis/cellular extravasation, antigen receptor mediated signaling pathway, T helper mediated immunity and macrophage activation; and the auxiliary elements included leukocyte mediated immunity, regulation of inflammatory response, T cell differentiation, mononuclear cell migration, megakaryocyte differentiation, complement activation and cytokine production. These deregulated immunological functions reveal the candidates to target in the immunotherapy.

Integrating the dysregulated inflammasome-based molecular functionome in the malignant transformation of endometriosis-associated ovarian carcinoma

The coexistence of endometriosis (ES) with ovarian clear cell carcinoma (CCC) or endometrioid carcinoma (EC) suggested that malignant transformation of ES leads to endometriosis associated ovarian carcinoma (EAOC). However, there is still lack of an integrating data analysis of the accumulated experimental data to provide the evidence supporting the hypothesis of EAOC transformation. Herein we used a function-based analytic model with the publicly available microarray datasets to investigate the expression profiling between ES, CCC, and EC. We analyzed the functional regularity pattern of the three type of samples and hierarchically clustered the gene sets to identify key mechanisms regulating the malignant transformation of EAOC. We identified a list of 18 genes (NLRP3, AIM2, PYCARD, NAIP, Caspase-4, Caspase-7, Caspase-8, TLR1, TLR7, TOLLIP, NFKBIA, TNF, TNFAIP3, INFGR2, P2RX7, IL-1B, IL1RL1, IL-18) closely related to inflammasome complex, indicating an important role of inflammation/immunity in EAOC transformation. We next explore the association between these target genes and patient survival using Gene Expression Omnibus (GEO), and found significant correlation between the expression levels of the target genes and the progression-free survival. Interestingly, high expression levels of AIM2 and NLRP3, initiating proteins of inflammasomes, were significantly correlated with poor progression-free survival. Immunohistochemistry staining confirmed a correlation between high AIM2 and high Ki-67 in clinical EAOC samples, supporting its role in disease progression. Collectively, we established a bioinformatic platform of gene-set integrative molecular functionome to dissect the pathogenic pathways of EAOC, and demonstrated a key role of dysregulated inflammasome in modulating the malignant transformation of EAOC.

Gene set-based analysis of mucinous ovarian carcinoma

OBJECTIVE

Mucinous ovarian carcinoma (MOC) is an uncommon subtype of epithelial ovarian cancers, and the pathogenesis is still poorly understood because of its rarity. We conducted a gene set-based analysis to investigate the pathogenesis of MOC by integrating microarray gene expression datasets based on the regularity of functions defined by gene ontology or canonical pathway databases.

MATERIALS AND METHODS

Forty-five pairs of MOC and normal ovarian tissue sample gene expression profiles were downloaded from the National Center for Biotechnology Information Gene Expression Omnibus database. The gene expression profiles were converted to the gene set regularity indexes by measuring the change of gene expression ordering in a gene set. Then the pathogenesis of MOC was investigated with the differences of function regularity with the gene set regularity indexes between the MOC and normal control samples.

RESULTS

The informativeness of the gene set regularity indexes was sufficient for machine learning to accurately recognize and classify the functional regulation patterns with an accuracy of 99.44%. The statistical analysis revealed that the GTPase regulators and receptor tyrosine kinase erbB-2 (ERBB2) were the most important aberrations; the exploratory factor analysis revealed phosphoinositide 3-kinase-activating kinase, G-protein coupled receptor pathway, oxidoreductase activity, immune response, peptidase activity, regulation of translation, and transport and channel activity were also involved in the pathogenesis of MOC.

CONCLUSION

Investigating the pathogenesis of MOC with the functionome provided a comprehensive view of the deregulated functions of this disease. In addition to GTPase regulators and ERBB2, a plenty of deregulated functions such as phosphoinositide 3-kinase, G-protein coupled receptor pathway, and immune response also participated in the interaction network of MOC pathogenesis.

Discovering the Deregulated Molecular Functions Involved in Malignant Transformation of Endometriosis to Endometriosis-Associated Ovarian Carcinoma Using a Data-Driven, Function-Based Analysis

The clinical characteristics of clear cell carcinoma (CCC) and endometrioid carcinoma EC) are concomitant with endometriosis (ES), which leads to the postulation of malignant transformation of ES to endometriosis-associated ovarian carcinoma (EAOC). Different deregulated functional areas were proposed accounting for the pathogenesis of EAOC transformation, and there is still a lack of a data-driven analysis with the accumulated experimental data in publicly-available databases to incorporate the deregulated functions involved in the malignant transformation of EOAC. We used the microarray gene expression datasets of ES, CCC and EC downloaded from the National Center for Biotechnology Information Gene Expression Omnibus (NCBI GEO) database. Then, we investigated the pathogenesis of EAOC by a data-driven, function-based analytic model with the quantified molecular functions defined by 1454 Gene Ontology (GO) term gene sets. This model converts the gene expression profiles to the functionome consisting of 1454 quantified GO functions, and then, the key functions involving the malignant transformation of EOAC can be extracted by a series of filters. Our results demonstrate that the deregulated oxidoreductase activity, metabolism, hormone activity, inflammatory response, innate immune response and cell-cell signaling play the key roles in the malignant transformation of EAOC. These results provide the evidence supporting the specific molecular pathways involved in the malignant transformation of EAOC.

Involvement of Chromatin Remodeling Genes and the Rho GTPases RhoB and CDC42 in Ovarian Clear Cell Carcinoma

Objective

Ovarian clear cell carcinomas (OCCCs) constitute a rare ovarian cancer subtype with distinct clinical features, but may nonetheless be difficult to distinguish morphologically from other subtypes. There is limited knowledge of genetic events driving OCCC tumorigenesis beyond ARID1A, which is reportedly mutated in 30–50% of OCCCs. We aimed to further characterize OCCCs by combined global transcriptional profiling and targeted deep sequencing of a panel of well-established cancer genes. Increased knowledge of OCCC-specific genetic aberrations may help in guiding development of targeted treatments and ultimately improve patient outcome.

Methods

Gene expression profiling of formalin-fixed, paraffin-embedded (FFPE) tissue from a cohort of the major ovarian cancer subtypes (cohort 1; n = 67) was performed using whole-genome cDNA-mediated Annealing, Selection, extension and Ligation (WG-DASL) bead arrays, followed by pathway, gene module score, and gene ontology analyses, respectively. A second FFPE cohort of 10 primary OCCCs was analyzed by targeted DNA sequencing of a panel of 60 cancer-related genes (cohort 2). Non-synonymous and non-sense variants affecting single-nucleotide variations and insertions or deletions were further analyzed. A tissue microarray of 43 OCCCs (cohort 3) was used for validation by immunohistochemistry and chromogenic in situ hybridization.

Results

Gene expression analyses revealed a distinct OCCC profile compared to other histological subtypes, with, e.g., ERBB2, TFAP2A, and genes related to cytoskeletal actin regulation being overexpressed in OCCC. ERBB2 was, however, not overexpressed on the protein level and ERBB2 amplification was rare in the validation cohort. Targeted deep sequencing revealed non-synonymous variants or insertions/deletions in 11/60 cancer-related genes. Genes involved in chromatin remodeling, including ARID1A, SPOP, and KMT2D were frequently mutated across OCCC tumors.

Conclusion

OCCCs appear genetically heterogeneous, but harbor frequent alterations in chromatin remodeling genes. Overexpression of TFAP2A and ERBB2 was observed on the mRNA level in relation to other ovarian cancer subtypes. However, overexpression of ERBB2 was not reflected by HER2 amplification or protein overexpression in the OCCC validation cohort. In addition, Rho GTPase-dependent actin organization may also play a role in OCCC pathogenesis and warrants further investigation. The distinct biological features of OCCC discovered here may provide a basis for novel targeted treatment strategies.

Potential Coagulation Factor-Driven Pro-Inflammatory Responses in Ovarian Cancer Tissues Associated with Insufficient O₂ and Plasma Supply

Tissue factor (TF) is a cell surface receptor for coagulation factor VII (fVII). The TF-activated fVII (fVIIa) complex is an essential initiator of the extrinsic blood coagulation process. Interactions between cancer cells and immune cells via coagulation factors and adhesion molecules can promote progression of cancer, including epithelial ovarian cancer (EOC). This process is not necessarily advantageous, as tumor tissues generally undergo hypoxia due to aberrant vasculature, followed by reduced access to plasma components such as coagulation factors. However, hypoxia can activate TF expression. Expression of fVII, intercellular adhesion molecule-1 (ICAM-1), and multiple pro-inflammatory cytokines can be synergistically induced in EOC cells in response to hypoxia along with serum deprivation. Thus, pro-inflammatory responses associated with the TF-fVIIa-ICAM-1 interaction are expected within hypoxic tissues. Tumor tissue consists of multiple components such as stromal cells, interstitial fluid, albumin, and other micro-factors such as proton and metal ions. These factors, together with metabolism reprogramming in response to hypoxia and followed by functional modification of TF, may contribute to coagulation factor-driven inflammatory responses in EOC tissues. The aim of this review was to describe potential coagulation factor-driven inflammatory responses in hypoxic EOC tissues. Arguments were extended to clinical issues targeting this characteristic tumor environment.