ABSTRACT: 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.
