Inference of Molecular Subtypes of Uterine Corpus Endometrioid Carcinoma with Different Survival Based on Cancer Hallmarks-associated Long Non-Coding RNAs and Gene Modules

Inference of Molecular Subtypes of Uterine Corpus Endometrioid Carcinoma with Different Survival Based on Cancer Hallmarks-associated Long Non-Coding RNAs and Gene Modules | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Research Inference of Molecular Subtypes of Uterine Corpus Endometrioid Carcinoma with Different Survival Based on Cancer Hallmarks-associated Long Non-Coding RNAs and Gene Modules Nan Lv, Shuyu Zhao, Yan Li, Tianyi Liu, Xiaodan Chu This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-37407/v1 This work is licensed under a CC BY 4.0 License Status: Posted Version 1 posted You are reading this latest preprint version Abstract Background: Uterine corpus endometrioid carcinoma (UCEC), a common gynecological malignancy with high incidence, affects the mental and physical health of women. It is increasingly evident that long non-coding RNAs (lncRNAs) have causative roles in cancers including UCEC. However, very little is known about the cancer hallmark-related risks for UCEC based on lncRNAs and genes. Methods: In this study, a computational integrated pipeline was development to evaluate cancer hallmark-related risk for each UCEC patient based on gene and lncRNA expression. Results: Some UCEC-specific cancer hallmark-related genes and lncRNAs were identified. Core modules were extracted from co-expressed UCEC-specific lncRNAs and genes networks for each cancer hallmark. Some core modules showed specific features and functions. Follow a multi-dimensional rank approach and core modules, each UCEC sample was given an integrated cancer hallmark-related risk score. We divided all the UCEC patients to diverse hallmark risk groups including multi-hallmark, media-hallmark and single-hallmark groups and they showed different prognosis. We also identified some key lncRNAs which participated in multiple kinds of cancer hallmarks. These key lncRNAs were associated with some essential pathways such as estrogen receptor signaling. Conclusions: In conclusion, the present study provide novel insights to the classification and treatment of UCEC. Endocrinology & Metabolism Long non-coding RNA cancer hallmark uterine corpus endometrioid carcinoma survival molecular subtypes Figures Figure 1 Figure 2 Figure 3 Figure 4 Figure 5 Figure 6 Background Uterine corpus endometrioid carcinoma (UCEC) is the sixth-most-common leading cause of cancer-related death among women in the United States[ 1 , 2 ], with an estimated 65,620 new cases and 12,590 deaths in 2020 [ 3 ]. Over time, the number of newly diagnosed UCEC in the United States has steadily increased and this trend is expected to continue [ 4 ]. Most UCEC women of early stages diagnosed show favorable outcomes. However, there are some low-grade, early, well differentiated UCEC women in which unexpected recurrence and adverse outcomes may occur. UCEC patients usually have diverse survival and drug response. Thus, UCEC is one of the few human malignant tumors for which mortality is increasing, which underscores the urgent need to construct novel and effective strategies for distinguishing subtypes with diverse survival and drug response of UCEC. Recent years, many kinds of non-coding RNAs have been discovered and researched in multiple kinds of diseases including cancer [ 5 , 6 ]. Long non-coding RNA (lncRNA) was a major kind of non-coding RNAs which >200 nt in length [ 7 , 8 ]. Accumulating evidence has suggested that lncRNAs serve as important and essential roles in physiological and pathological processes of various cancers [ 9 , 10 ]. Numerous studies have reported that the occurrence, development, prognosis and treatment of UCEC were closely related to the abnormal expression of a large number of lncRNAs. For example, lncRNA SNHG16 induced by TFAP2A modulates glycolysis and proliferation of UCEC and is associated with poor survival rate [ 11 ]. Yao et al. reported that lncRNA LSINCT5 is up-regulated and significantly inhibited cell proliferation, cell cycle progression, and induced apoptosis in UCEC. Thus, comprehensive and systematic exploration of the characteristics about lncRNA in UCEC could provide assistance for identifying novel molecular associations of potential mechanistic significance in the development of UCEC. It is generally known that the biology of cancer is extremely complex, individualized and various. However, some key traits have been revealed to reduce cancer complexity during the past decade. These traits could be represented by a few distinctive and complementary capabilities (which are considered as “cancer hallmarks”) that could promote tumor growth and metastasis. These cancer hallmarks could provide a logical framework for understanding the significant diversity of multiple kinds of cancers [ 12 ]. Hanahan and Weinberg proposed six hallmarks and ten hallmarks of cancer in 2000 and 2011 [ 13 , 14 ]. These two researches both suggest that multiple kinds of cancers share some common cancer hallmarks which dominate the convert from normal cells to cancer cells. These ten cancer hallmarks are effective principles to depict the characteristics of cancers. However, we have little insight into how to use these hallmarks to depict the roles of lncRNAs and distinguish UCEC patients. In present study, we developed an integrated algorithm for evaluating hallmark-related risks of UCEC patients and distinguishing them to diverse groups based on gene and lncRNA expression profiles (Figure 1). Some UCEC-specific lncRNAs and genes in ten kinds of hallmarks were identified. Co-expressed UCEC-specific lncRNAs and genes networks in diverse cancer hallmark were constructed. Some core modules were extracted from these co-expressed UCEC-specific lncRNAs and genes networks. Each UCEC patient was given a comprehensive score for a specific core module in each cancer hallmark based on multi-dimensional rank approach. All the UCEC patients were divided to diverse hallmark risk groups including multi-hallmark, media-hallmark and single-hallmark groups. These diverse hallmark risk groups were associated with different survival. Some key lncRNAs which were related to multiple cancer hallmarks were also discovered and showed specific functions. Collectively, this study clarified the roles of cancer hallmark-related lncRNAs in UCSC and distinguished them to diverse risk groups. Methods Collection of gene and lncRNA expression profiles of UCEC patients The lncRNA expression and gene expression (level 3) data, as well as clinical data of UCEC patients were obtained from The Cancer Genome Atlas (TCGA, Release: 2019-07-21). The download link of these data was https://gdc.xenahubs.net/download/TCGA-UCEC.htseq_fpkm.tsv.gz. We also obtained genome annotation data, including genome sites and symbols of genes and lncRNAs from GENCODE 31 (19.06.19). The dataset includes 545 UCEC tumor and 38 control samples. All the lncRNAs and genes which were not expressed in all samples would be excluded. The minimum value of all samples were given to any remaining expression values of 0. All the expression values were transformed follow log2(value+1). Gathering cancer hallmark-related GO terms and genes All the cancer hallmark-related GO terms were obtained from a previous study [ 15 ]. We downloaded all the genes of these cancer hallmark-related GO terms from Gene Ontology using AmiGo (version: 2.5.12; http://amigo.geneontology.org/amigo) [ 16 ]. Thus, all the cancer hallmark-related genes were got for follow analysis. Identification of UCEC-specific cancer hallmark-related genes and lncRNAs test was applied to identify differential expressed cancer hallmark-related genes and lncRNAs between UCEC and control normal samples based on expression profiles. lncRNAs and cancer hallmark-related genes were considered as UCEC-specific lncRNAs and cancer hallmark-related genes if they were differentially expressed (P < 0.01). All the UCEC-specific lncRNAs and cancer hallmark-related genes were divided to up- and down-regulated lncRNAs and genes based on fold values. Construction of UCEC-specific cancer hallmark-related genes and lncRNAs co-expressed networks and identification of core modules We calculated Pearson’s correlation coefficients (PCCs) for each UCEC-specific cancer hallmark-related gene and lncRNA pair in all cancer hallmarks. The UCEC-specific gene and lncRNA pairs which their PCCs were higher than 0.3 or smaller than -0.3 and the p-values were smaller than 0.01 were extracted and considered as co-expressed UCEC-specific gene and lncRNA pairs. UCEC-specific cancer hallmark-related genes and lncRNAs co-expressed networks were constructed by Cytoscape 3.3.0 (http://www.cytoscape.org/) based on co-expressed pairs which their PCCs were higher than 0.5 or smaller than -0.5. In order to more accurately identify cancer hallmark-related genes and lncRNAs pairs in UCEC, we extracted core modules from co-expressed networks using MCODE module in cytoscape. Thus, certain numbers of core modules were identified for each cancer hallmark in UCEC. Evaluating hallmark-related risk for UCEC patients and dividing them to diverse risk groups based on core modules For each UCEC patient, risk score in each cancer hallmark was calculated and given. Multidimensional rank approach was applied to obtain risk score for each cancer hallmark based on genes and lncRNAs expression values in core modules. Up- and down-regulated genes and lncRNAs ranked in positive and negative order. 1000 permutation test was performed follow randomly perturbing cancer samples. The UCEC patient was considered as a risk sample for a specific cancer hallmark if permutation P values was smaller than 0.05 in any core module. Then, all the UCEC patients were divided to non-hallmark (0 risk hallmark), media-hallmark (1-5 risk hallmarks) and multi-hallmark groups (6-8 risk hallmarks). Survival analysis for diverse risk cancer hallmark groups and functional analysis for UCEC-specific genes and lncRNAs Kaplan-Meier survival analysis was performed for diverse risk cancer hallmark groups. Log-rank test was using to asses statistical significance (P< 0.05). R 3.6.2 framework was performed for all analyses. Enrichr tools (http://amp.pharm.mssm.edu/Enrichr) with default parameters was used for functional analysis based on UCEC-specific genes in each cancer hallmark [ 17 ]. Significantly enriched pathways (P< 0.05) were obtained for each cancer hallmark in UCEC. Results Some lncRNAs and cancer hallmark-related genes were specific and co-expressed in UCEC In order to depict the roles of cancer hallmark in UCEC, we identified UCEC-specific lncRNAs and cancer hallmark-related genes by differential expression. 1880 (12.43%) UCEC-specific lncRNAs were identified between UCEC and control samples (Figure 2A). These UCEC-specific lncRNAs included 891 and 989 up- and down-regulated lncRNAs (Figure 2B). UCEC-specific genes were also identified for each cancer hallmark (Figure 2C). In all kinds of cancer hallmarks, there were more than 50% UCEC-specific genes. Specially, there was 70.57% UCEC-specific genes in cancer hallmark genome instability and mutation. The results indicated that these cancer hallmark-related genes maybe serve as essential roles in UCEC. We assumed that lncRNAs and cancer hallmark-related genes could function by cooperating in UCEC. Thus, some co-expressed lncRNAs and cancer hallmark-related genes pairs were identified in UCEC for each kind of cancer hallmark. Most PCCs of co-expressed lncRNAs and cancer hallmark-related genes pairs were concentrated between 0.3 and 0.5 (Figure 2D). There were 8631 co-expressed lncRNAs and cancer hallmark-related genes pairs which their PCCs were higher than 0.5 (Figure 2E). In each kind of cancer hallmark, the numbers of pairs, lncRNAs and genes were diverse (Figure 2F). For example, there were 2525, 365 and 384 pairs, genes and lncRNAs in cancer hallmark self sufficiency in growth signals. However, only there were 84, 47 and 12 pairs, genes and lncRNAs in cancer hallmark reprogramming energy metabolism. The results indicated that diverse cancer hallmarks play different roles in UCEC. All above results indicated that lncRNAs and cancer hallmark-related genes cooperative pairs were important in UCEC. Some core modules were extracted from co-expressed lncRNAs and cancer hallmark-related genes networks in each cancer hallmark For each cancer hallmark, co-expressed lncRNAs and cancer hallmark-related genes which their PCCs were higher than 0.5 were extracted for constructing co-expressed networks. In cancer hallmark evading apoptosis, lncRNAs and cancer hallmark-related genes co-expressed network was constructed (Figure 3A). This co-expressed network contained 380 nodes (234 UCEC-specific lncRNAs and 147 cancer hallmark-related genes) and 906 edges. We found that some cancer hallmark-related genes and lncRNAs such as SLC25A27, AC005288 and GD5-AS1 played core roles in this co-expressed network. Most of cancer hallmark-related genes and lncRNAs showed positive correlations in UCEC. In eight cancer hallmarks, there were diverse numbers of core modules were identified (Figure 3B). Cancer hallmarks insensitivity to antigrowth signals, self sufficiency in growth signals and tissue invasion and metastasis had most core modules. The numbers of cancer hallmark-related genes and lncRNAs were also different (Figure 3C). For example, there were more lncRNA in core module 1 in cancer hallmark insensitivity to antigrowth signals. These core modules maybe show specific functions. For example, a core module in cancer hallmark evading apoptosis contained three lncRNAs and three genes (Figure 3D). PSMB8, PSMB10, PSME2 and PSMB8-AS1 were all proteasome-related genes or lncRNAs. The proteasome is a multicatalytic proteinase complex which is characterized by its ability to cleave peptides. Specially, cancer hallmark-related gene PSMB8 and lncRNA PSMB8-AS1 showed strong positive correlation (P< 0.001, PCC=0.79). These genes and lncRNAs in this core module showed close interactions. Another core module also showed close interactions (Figure 3E). All the results explained these core modules in co-expressed networks could function and serve as specific biomarks in UCEC. Specific cancer hallmark-related risk were evaluated for each UCEC patient based on core modules We inferred that each UCEC patient maybe have diverse hallmark-related risk. Thus, we calculated risk scores for each UCEC patient based on core modules in each cancer hallmark. Only eight cancer hallmarks were extracted for calculating risk scores due to core modules. The density distribution of risk scores in all core modules of each cancer hallmark were similar (Figure 4A). Only a small number of UCEC patients showed higher risk scores. The differences of average risk scores for diverse top 20 core modules were also present (Figure 4B). These 20 core modules were significantly associated with more UCEC patients (Figure 4C). For example, there were 32.03% samples were significant in core module 1 in cancer hallmark insensitivity to antigrowth signals. This core module was a key module which had most related samples and also been explained in above results. We also discovered the percent of significant risk-related samples in 0%corresponding top ranked samples. 70% UCEC samples ranked before corresponding orders in most core modules (Figure 4D). These results indicated that UCEC patients showed differences of cancer hallmark-related risk scores. UCEC groups with diverse cancer hallmark risk showed specific features Each UCEC patient could be associated with diverse cancer hallmark follow above pipeline. The numbers of UCEC in each cancer hallmark were different (Figure 5A). For example, there were more than 350 UCEC samples were associated with cancer hallmark tissue invasion and metastasis. The cancer hallmark reprogramming energy metabolism was related to 300 UCEC samples. Specially, some UCEC patients were associated with multiple cancer hallmarks. 13.59% UCEC patients had some relationships with any kinds of cancer hallmarks (Figure 5B). 11.55% UCEC patients were only related to one kind of cancer hallmark. Thus, we could divide all the UCEC patients to diverse groups with different numbers of cancer hallmarks (Figure 5C). The three diverse groups contained non-, media- and multi-hallmarks. 70.76% samples belonged to media-hallmarks groups. The UCEC patients with more cancer hallmarks usually had better survival days (Figure 5D). In addition, we also divided all the UCEC patients to another three cancer hallmark-related risk groups based on hierarchical clustering (Figure 5E). We also discovered that these diverse cancer hallmark-related risk groups showed different prognosis (Figure 5F, G). Group 2 significantly had better survival than group 1 and 3 (P= 0.014 and 0.023). All the results suggested that diverse hallmark-related risk groups had respective features and prognosis. Some key lncRNAs could participate in multiple kinds of cancer hallmarks and showed specific functions In order to further depict the roles of lncRNAs in cancer hallmarks for UCEC, we extracted some key lncRNAs which participate in multiple kinds of cancer hallmarks. There were 11 lncRNAs were associated with more than three kinds of cancer hallmarks (Figure 6A). lncRNAs AL590764.1, ANKRD10-IT1, NORD and AP000766.1 could participated in four kinds of cancer hallmarks (Figure 6B). However, the classes of these four kinds of cancer hallmarks were diverse. We inferred that these lncRNAs may serve as essential roles in UCEC. Thus, we further performed functional analysis for these lncRNAs in each kind of cancer hallmark. The lncRNAs were enrichment in some essential functions for cancer development (Figure 6C). For example, lncRNAs associated with evading apoptosis were enrichment in some hormone-related pathways such as negative regulation of intracellular estrogen receptor signaling, negative regulation of intracellular steroid hormone receptor signaling and regulation of intracellular estrogen receptor signaling. The estrogen receptor status is reported to be an important marker of UCEC [ 18 , 19 ]. In addition, we found most lncRNAs were associated with gap junction assembly pathway. Nishimura M et al. reported that gap junctional intercellular communication was suppressed via 5' CpG island methylation in promoter region of E-cadherin gene in UCEC cells [ 20 ]. These results indicated that cancer hallmark-related lncRNAs could serve as essential roles in UCEC. Discussion In this study, a calculated integrated approach which evaluate cancer hallmark-risk for each UCEC patient based on gene and lncRNA expression was developed. Some core modules were extracted from co-expressed UCEC-specific lncRNAs and genes networks in diverse cancer hallmark. All the UCEC patients were divided to diverse hallmark risk groups based on respective comprehensive risk scores. These diverse hallmark risk groups showed different survival. Some key lncRNAs were related to multiple cancer hallmarks and showed specific functions. Dividing cancer samples to groups with diverse molecular characteristics could provide assistance for cancer diagnosis and individualized treatment. Wright GW et al. described an algorithm that determines the probability that a patient's lymphoma belongs to one of seven genetic subtypes based on its genetic features [ 21 ]. Lee S et al. presented an explainable deep learning model with attention mechanism and network propagation for cancer subtype classification [ 22 ]. Most of these developed clarification methods only focused on coding genes. In our work, lncRNAs were considered as essential and key factors for constructing the computational approach about cancer patients classification. The computational approach not only divided UCEC patients to diverse hallmark risk groups but also discovered some key lncRNAs which serve as important roles in UCEC. These lncRNAs were significantly enriched in some essential cancer development functions. Analysis of cancer hallmarks could greatly improve our understanding of the occurrence, development and metastasis of many cancer types. In present study, we applied cancer hallmarks to evaluate and distinguish UCEC patients. We discovered that most cancer hallmark-related genes were differentially expressed in UCEC. It indicated that cancer hallmark-related genes maybe had essential functions in UCEC. Specially, we found that each UCEC patient showed obvious differences for risk of cancer hallmarks. Although cancer hallmarks were considered as common properties for cancer, each UCEC patient still showed personalized features. This result supported that there were great differences among different cancer patients and the application of personalized medicine in cancer treatment is very important. Thus, our study could help to evaluate risk for UCEC patient and establishes personalized therapy. Conclusions In summary, the present study developed a standardized computational procedures to evaluate cancer patient hallmark-related risk based on genes and lncRNAs expression. All the UCEC patients were divided to diverse cancer hallmark-related risk groups and showed different prognosis. Some key cancer hallmark-related lncRNAs were significantly associated with UCEC development. Collectively, our study leads to a novel starting point for future functional explorations, the identification of biomarkers, and lncRNA-based targeted therapy for UCEC. Abbreviations UCEC: uterine corpus endometrioid carcinoma; lncRNA: long non-coding RNA; TCGA: The Cancer Genome Atlas; GO: gene ontology; PCC: Pearson correlation coefficients. Declarations Funding Not applicable. Disclosure of interest The authors declare that they have no competing interest. Acknowledgments Not applicable. Availability of data and material The data that support the findings of this study are available from the corresponding author upon reasonable request. Ethics approval and consent to participate Not applicable. Patient consent for publication Not applicable. Competing interests The authors declare that they have no competing interests. Consent for publication Not applicable. Authors' contributions LN conceived and designed the experiments, ZSY, LY, LTY and CXD analysed the data, and LN and ZSY wrote the manuscript. References Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A: Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin 2018, 68: 394-424. Gaber C, Meza R, Ruterbusch JJ, Cote ML: Endometrial Cancer Trends by Race and Histology in the USA: Projecting the Number of New Cases from 2015 to 2040. J Racial Ethn Health Disparities 2016. Siegel RL, Miller KD, Jemal A: Cancer statistics, 2020. CA Cancer J Clin 2020, 70: 7-30. Bell DW, Ellenson LH: Molecular Genetics of Endometrial Carcinoma. Annu Rev Pathol 2019, 14: 339-367. Anastasiadou E, Jacob LS, Slack FJ: Non-coding RNA networks in cancer. Nat Rev Cancer 2018, 18: 5-18. Kaikkonen MU, Adelman K: Emerging Roles of Non-Coding RNA Transcription. Trends Biochem Sci 2018, 43: 654-667. Ransohoff JD, Wei Y, Khavari PA: The functions and unique features of long intergenic non-coding RNA. Nat Rev Mol Cell Biol 2018, 19: 143-157. Uszczynska-Ratajczak B, Lagarde J, Frankish A, Guigo R, Johnson R: Towards a complete map of the human long non-coding RNA transcriptome. Nat Rev Genet 2018, 19: 535-548. Wang Y, Fang Z, Hong M, Yang D, Xie W: Long-noncoding RNAs (lncRNAs) in drug metabolism and disposition, implications in cancer chemo-resistance. Acta Pharm Sin B 2020, 10: 105-112. Wu P, Mo Y, Peng M, Tang T, Zhong Y, Deng X, Xiong F, Guo C, Wu X, Li Y, et al: Emerging role of tumor-related functional peptides encoded by lncRNA and circRNA. Mol Cancer 2020, 19: 22. Zhang G, Ma A, Jin Y, Pan G, Wang C: LncRNA SNHG16 induced by TFAP2A modulates glycolysis and proliferation of endometrial carcinoma through miR-490-3p/HK2 axis. Am J Transl Res 2019, 11: 7137-7145. Wang E, Zaman N, McGee S, Milanese JS, Masoudi-Nejad A, O'Connor-McCourt M: Predictive genomics: a cancer hallmark network framework for predicting tumor clinical phenotypes using genome sequencing data. Semin Cancer Biol 2015, 30: 4-12. Hanahan D, Weinberg RA: The hallmarks of cancer. Cell 2000, 100: 57-70. Hanahan D, Weinberg RA: Hallmarks of cancer: the next generation. Cell 2011, 144: 646-674. Plaisier CL, Pan M, Baliga NS: A miRNA-regulatory network explains how dysregulated miRNAs perturb oncogenic processes across diverse cancers. Genome Res 2012, 22: 2302-2314. Carbon S, Ireland A, Mungall CJ, Shu S, Marshall B, Lewis S, Ami GOH, Web Presence Working G: AmiGO: online access to ontology and annotation data. Bioinformatics 2009, 25: 288-289. Kuleshov MV, Jones MR, Rouillard AD, Fernandez NF, Duan Q, Wang Z, Koplev S, Jenkins SL, Jagodnik KM, Lachmann A, et al: Enrichr: a comprehensive gene set enrichment analysis web server 2016 update. Nucleic Acids Res 2016, 44: W90-97. Sho T, Hachisuga T, Nguyen TT, Urabe R, Kurita T, Kagami S, Kawagoe T, Matsuura Y, Shimajiri S: Expression of estrogen receptor-alpha as a prognostic factor in patients with uterine serous carcinoma. Int J Gynecol Cancer 2014, 24: 102-106. Togami S, Sasajima Y, Oi T, Ishikawa M, Onda T, Ikeda S, Kato T, Tsuda H, Kasamatsu T: Clinicopathological and prognostic impact of human epidermal growth factor receptor type 2 (HER2) and hormone receptor expression in uterine papillary serous carcinoma. Cancer Sci 2012, 103: 926-932. Nishimura M, Saito T, Yamasaki H, Kudo R: Suppression of gap junctional intercellular communication via 5' CpG island methylation in promoter region of E-cadherin gene in endometrial cancer cells. Carcinogenesis 2003, 24: 1615-1623. Wright GW, Huang DW, Phelan JD, Coulibaly ZA, Roulland S, Young RM, Wang JQ, Schmitz R, Morin RD, Tang J, et al: A Probabilistic Classification Tool for Genetic Subtypes of Diffuse Large B Cell Lymphoma with Therapeutic Implications. Cancer Cell 2020, 37: 551-568 e514. Lee S, Lim S, Lee T, Sung I, Kim S: Cancer subtype classification and modeling by pathway attention and propagation. Bioinformatics 2020. Cite Share Download PDF Status: Posted Version 1 posted You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-37407","acceptedTermsAndConditions":true,"allowDirectSubmit":true,"archivedVersions":[],"articleType":"Research","associatedPublications":[],"authors":[{"id":725907,"identity":"cc904ee6-c057-4042-bbbd-45ecf9c3c1de","order_by":0,"name":"Nan Lv","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAAA5klEQVRIiWNgGAWjYBACNv7mww8+/qmR42dvbHyQUFFDWAufxLE0w5kNx4wlew43Gzw4c4ywFjmGHANpzgbmxA030tskH7YwE+EwhjMGxow72IwZDiS2VSQ2sDHwt3cn4NfC3FbwuPCMjBxjw8G2G4k7ZBgkzpzdQMCWwxuMZ7CxGTMzNgK1nGFjMJDIJaQlwUCah405sY2Zsa0ASBKjJcVAmreNObGHjbGNgTgtoECeceaYsQQPY7NEwpljPAT9It8PjMoPFTVy9vefP/z4A8jgb+/FrwUD8JCmfBSMglEwCkYBVgAAje5Mu/KV1pwAAAAASUVORK5CYII=","orcid":"","institution":"Second Affiliated Hospital of Harbin Medical University","correspondingAuthor":true,"prefix":"","firstName":"Nan","middleName":"","lastName":"Lv","suffix":""},{"id":725908,"identity":"62145c2c-8bad-44a7-aeaa-79b01c21dc49","order_by":1,"name":"Shuyu Zhao","email":"","orcid":"","institution":"Department of Endocrinology, Second Affiliated Hospital of Harbin Medical University","correspondingAuthor":false,"prefix":"","firstName":"Shuyu","middleName":"","lastName":"Zhao","suffix":""},{"id":725909,"identity":"ce1e755a-259e-481b-bc35-d922c33ab965","order_by":2,"name":"Yan Li","email":"","orcid":"","institution":"Department of Endocrinology, Second Affiliated Hospital of Harbin Medical University","correspondingAuthor":false,"prefix":"","firstName":"Yan","middleName":"","lastName":"Li","suffix":""},{"id":725910,"identity":"b645c3dd-93e3-4218-b915-340eef99654c","order_by":3,"name":"Tianyi Liu","email":"","orcid":"","institution":"Department of Endocrinology, Second Affiliated Hospital of Harbin Medical University","correspondingAuthor":false,"prefix":"","firstName":"Tianyi","middleName":"","lastName":"Liu","suffix":""},{"id":725911,"identity":"a49ac575-421c-4f7c-9d8f-22af2fce8c04","order_by":4,"name":"Xiaodan Chu","email":"","orcid":"","institution":"Department of Endocrinology, Second Affiliated Hospital of Harbin Medical University","correspondingAuthor":false,"prefix":"","firstName":"Xiaodan","middleName":"","lastName":"Chu","suffix":""}],"badges":[],"createdAt":"2020-06-22 06:00:33","currentVersionCode":1,"declarations":"","doi":"10.21203/rs.3.rs-37407/v1","doiUrl":"https://doi.org/10.21203/rs.3.rs-37407/v1","draftVersion":[],"editorialEvents":[],"editorialNote":"","failedWorkflow":false,"files":[{"id":1462196,"identity":"ec6c1097-abe1-438d-8012-cbf6debf6680","added_by":"auto","created_at":"2020-06-30 19:41:58","extension":"jpg","order_by":1,"title":"Figure 1","display":"","copyAsset":false,"role":"figure","size":1663760,"visible":true,"origin":"","legend":"The workflow of evaluating cancer hallmark-related risk for UCEC patients based on genes and lncRNAs expression. Step 1. Identification of UCEC-specific cancer hallmark-related genes and lncRNAs and calculating co-expression pairs based on genes and lncRNAs expression profiles. Step 2. Construction of UCEC-specific cancer hallmark-related genes and lncRNAs co-expressed networks and extraction of core modules. Step 2. Evaluating cancer hallmark-related risks for UCEC patients and dividing them to diverse groups.","description":"","filename":"Fig1.jpg","url":"https://assets-eu.researchsquare.com/files/rs-37407/v1/Fig1.jpg"},{"id":1462197,"identity":"268780f9-e498-472d-9608-0bdbb549fa7f","added_by":"auto","created_at":"2020-06-30 19:41:58","extension":"jpg","order_by":2,"title":"Figure 2","display":"","copyAsset":false,"role":"figure","size":1299885,"visible":true,"origin":"","legend":"Some lncRNAs and cancer hallmark-related genes were specific and co-expressed in UCEC. (A) The pie chart shows percent of differential expressed lncRNAs in UCEC. (B) The bar plots show numbers of up- (red) and down-regulated (blue) lncRNAs. (C) Pie charts show percents of differential expressed genes in UCEC for each kind of cancer hallmark. (D) The violin plots show PCCs between UCEC-specific lncRNAs and genes in each cancer hallmark. (E) The bar charts show numbers of co-expressed UCEC-specific lncRNAs and genes pairs. (F) The size of the circle represent the number of co-expressed UCEC-specific lncRNAs and genes pairs, genes and lncRNAs in each cancer hallmark.","description":"","filename":"Fig2.jpg","url":"https://assets-eu.researchsquare.com/files/rs-37407/v1/Fig2.jpg"},{"id":1462198,"identity":"8a59b2ef-6c4b-4939-bdc3-4c8d16c842e2","added_by":"auto","created_at":"2020-06-30 19:41:58","extension":"jpg","order_by":3,"title":"Figure 3","display":"","copyAsset":false,"role":"figure","size":2963384,"visible":true,"origin":"","legend":"Extraction of core modules from co-expressed lncRNAs and genes networks in each cancer hallmark. (A) Co-expressed UCEC-specific lncRNAs and genes networks in cancer hallmark evading apoptosis. Green and yellow nodes represent lncRNAs and genes, respectively. (B) The radar chart shows numbers of core modules in each cancer hallmark. (C) The rose diagram shows numbers of lncRNAs and genes in each core module. (D, E) Two core modules as examples in cancer hallmark evading apoptosis.","description":"","filename":"Fig3.jpg","url":"https://assets-eu.researchsquare.com/files/rs-37407/v1/Fig3.jpg"},{"id":1462199,"identity":"0d83d454-0f36-4e86-b342-355b75c73053","added_by":"auto","created_at":"2020-06-30 19:41:58","extension":"jpg","order_by":4,"title":"Figure 4","display":"","copyAsset":false,"role":"figure","size":2444078,"visible":true,"origin":"","legend":"Evaluating specific cancer hallmark-related risk for each UCEC patient based on core modules. (A) Density curves of risk scores for all core modules in each cancer hallmark. (B) The bar charts show average risk scores for top 20 core modules with significant cancer hallmark-related UCEC samples. (C) The diagram shows percents of significant risk UCEC samples using diverse core modules. (D) The point plot shows percents of UCEC samples ranked before corresponding orders in all core modules.","description":"","filename":"Fig4.jpg","url":"https://assets-eu.researchsquare.com/files/rs-37407/v1/Fig4.jpg"},{"id":1462200,"identity":"6fbdbfbb-c17d-410c-b464-224532810e31","added_by":"auto","created_at":"2020-06-30 19:41:59","extension":"jpg","order_by":5,"title":"Figure 5","display":"","copyAsset":false,"role":"figure","size":1701579,"visible":true,"origin":"","legend":"Specific features of UCEC groups with diverse cancer hallmark risk. (A) The bar plot shows numbers of samples related to diverse cancer hallmark. (B) The pie chart shows percents of UCEC samples with different numbers of cancer hallmarks. (C) The numbers of UCEC samples in diverse cancer hallmark-related risk groups. (D) The box plots show overall survival days in diverse cancer hallmark-related risk groups. (E) Clustering trees for UCEC samples. (F, G) The Kaplan-Meier curve for the overall survival of two cancer hallmark-related risk groups. The difference between the two curves was evaluated by a two-sided log-rank test.","description":"","filename":"Fig5.jpg","url":"https://assets-eu.researchsquare.com/files/rs-37407/v1/Fig5.jpg"},{"id":1462201,"identity":"cc9bb4bc-55e7-4f9e-89a8-2a491a312ca0","added_by":"auto","created_at":"2020-06-30 19:41:59","extension":"jpg","order_by":6,"title":"Figure 6","display":"","copyAsset":false,"role":"figure","size":3217098,"visible":true,"origin":"","legend":"Some key lncRNAs could participate in multiple kinds of cancer hallmarks and showed specific functions. (A) The heatmap shows lncRNAs if are associated with cancer hallmark. (B) The bar plots show numbers of cance hallmarks which lncRNAs participate in. (C) Enrichment pathways of lncRNAs in each cancer hallmark. The bar plots represent -log (P).","description":"","filename":"Fig6.jpg","url":"https://assets-eu.researchsquare.com/files/rs-37407/v1/Fig6.jpg"},{"id":13548593,"identity":"27d976c9-96c2-48d4-8ba5-5eac55affa79","added_by":"auto","created_at":"2021-09-17 02:18:06","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":1934773,"visible":true,"origin":"","legend":"","description":"","filename":"manuscript.pdf","url":"https://assets-eu.researchsquare.com/files/rs-37407/v1/dbeada26-6b99-4940-8212-e45e187a127c.pdf"}],"financialInterests":"","formattedTitle":"\u003cp\u003eInference of Molecular Subtypes of Uterine Corpus Endometrioid Carcinoma with Different Survival Based on Cancer Hallmarks-associated Long Non-Coding RNAs and Gene Modules\u003c/p\u003e","fulltext":[{"header":"Background ","content":"\u003cp\u003eUterine corpus endometrioid carcinoma (UCEC) is the sixth-most-common leading cause of cancer-related death among women in the United States[\u003ca href=\"#_ENREF_1\"\u003e1\u003c/a\u003e, \u003ca href=\"#_ENREF_2\"\u003e2\u003c/a\u003e], with an estimated 65,620 new cases and 12,590 deaths in 2020 [\u003ca href=\"#_ENREF_3\"\u003e3\u003c/a\u003e]. Over time, the number of newly diagnosed UCEC in the United States has steadily increased and this trend is expected to continue [\u003ca href=\"#_ENREF_4\"\u003e4\u003c/a\u003e]. Most UCEC women of early stages diagnosed show favorable outcomes. However, there are some low-grade, early, well differentiated UCEC women in which unexpected recurrence and adverse outcomes may occur. UCEC patients usually have diverse survival and drug response. Thus, UCEC is one of the few human malignant tumors for which mortality is increasing, which underscores the urgent need to construct novel and effective strategies for distinguishing subtypes with diverse survival and drug response of UCEC.\u003c/p\u003e\n\u003cp\u003eRecent years, many kinds of non-coding RNAs have been discovered and researched in multiple kinds of diseases including cancer [\u003ca href=\"#_ENREF_5\"\u003e5\u003c/a\u003e, \u003ca href=\"#_ENREF_6\"\u003e6\u003c/a\u003e]. Long non-coding RNA (lncRNA) was a major kind of non-coding RNAs which \u0026gt;200 nt in length [\u003ca href=\"#_ENREF_7\"\u003e7\u003c/a\u003e, \u003ca href=\"#_ENREF_8\"\u003e8\u003c/a\u003e]. Accumulating evidence has suggested that lncRNAs serve as important and essential roles in physiological and pathological processes of various cancers [\u003ca href=\"#_ENREF_9\"\u003e9\u003c/a\u003e, \u003ca href=\"#_ENREF_10\"\u003e10\u003c/a\u003e]. Numerous studies have reported that the occurrence, development, prognosis and treatment of UCEC were closely related to the abnormal expression of a large number of lncRNAs. For example, lncRNA SNHG16 induced by TFAP2A modulates glycolysis and proliferation of UCEC and is associated with poor survival rate [\u003ca href=\"#_ENREF_11\"\u003e11\u003c/a\u003e]. Yao et al. reported that lncRNA LSINCT5 is up-regulated and significantly inhibited cell proliferation, cell cycle progression, and induced apoptosis in UCEC. Thus, comprehensive and systematic exploration of the characteristics about lncRNA in UCEC could provide assistance for identifying novel molecular associations of potential mechanistic significance in the development of UCEC.\u003c/p\u003e\n\u003cp\u003eIt is generally known that the biology of cancer is extremely complex, individualized and various. However, some key traits have been revealed to reduce cancer complexity during the past decade. These traits could be represented by a few distinctive and complementary capabilities (which are considered as \u0026ldquo;cancer hallmarks\u0026rdquo;) that could promote tumor growth and metastasis. These cancer hallmarks could provide a logical framework for understanding the significant diversity of multiple kinds of cancers [\u003ca href=\"#_ENREF_12\"\u003e12\u003c/a\u003e]. Hanahan and Weinberg proposed six hallmarks and ten hallmarks of cancer in 2000 and 2011 [\u003ca href=\"#_ENREF_13\"\u003e13\u003c/a\u003e, \u003ca href=\"#_ENREF_14\"\u003e14\u003c/a\u003e]. These two researches both suggest that multiple kinds of cancers share some common cancer hallmarks which dominate the convert from normal cells to cancer cells. These ten cancer hallmarks are effective principles to depict the characteristics of cancers. However, we have little insight into how to use these hallmarks to depict the roles of lncRNAs and distinguish UCEC patients.\u003c/p\u003e\n\u003cp\u003eIn present study, we developed an integrated algorithm for evaluating hallmark-related risks of UCEC patients and distinguishing them to diverse groups based on gene and lncRNA expression profiles (Figure 1). Some UCEC-specific lncRNAs and genes in ten kinds of hallmarks were identified. Co-expressed\u0026nbsp;\u0026nbsp; UCEC-specific lncRNAs and genes networks in diverse cancer hallmark were constructed. Some core modules were extracted from these co-expressed\u0026nbsp;\u0026nbsp; UCEC-specific lncRNAs and genes networks. Each UCEC patient was given a comprehensive score for a specific core module in each cancer hallmark based on multi-dimensional rank approach. All the UCEC patients were divided to diverse hallmark risk groups including multi-hallmark, media-hallmark and single-hallmark groups. These diverse hallmark risk groups were associated with different survival. Some key lncRNAs which were related to multiple cancer hallmarks were also discovered and showed specific functions. Collectively, this study clarified the roles of cancer hallmark-related lncRNAs in UCSC and distinguished them to diverse risk groups.\u003c/p\u003e"},{"header":"Methods","content":"\u003cp\u003e\u003cstrong\u003eCollection of gene and lncRNA expression profiles of UCEC patients\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe lncRNA expression and gene expression (level 3) data, as well as clinical data of UCEC patients were obtained from The Cancer Genome Atlas (TCGA, Release: 2019-07-21). The download link of these data was https://gdc.xenahubs.net/download/TCGA-UCEC.htseq_fpkm.tsv.gz. We also obtained genome annotation data, including genome sites and symbols of genes and lncRNAs from GENCODE 31 (19.06.19). The dataset includes 545 UCEC tumor and 38 control samples. All the lncRNAs and genes which were not expressed in all samples would be excluded. The minimum value of all samples were given to any remaining expression values of 0. All the expression values were transformed follow log2(value+1).\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eGathering cancer hallmark-related GO terms and genes\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eAll the cancer hallmark-related GO terms were obtained from a previous study [\u003ca href=\"#_ENREF_15\"\u003e15\u003c/a\u003e]. We downloaded all the genes of these cancer hallmark-related GO terms from Gene Ontology using AmiGo (version: 2.5.12; http://amigo.geneontology.org/amigo) [\u003ca href=\"#_ENREF_16\"\u003e16\u003c/a\u003e]. Thus, all the cancer hallmark-related genes were got for follow analysis.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eIdentification of UCEC-specific cancer hallmark-related genes and lncRNAs\u003c/strong\u003e\u003c/p\u003e\n\u003cul\u003e\n\u003cli\u003etest was applied to identify differential expressed cancer hallmark-related genes and lncRNAs between UCEC and control normal samples based on expression profiles. lncRNAs and cancer hallmark-related genes were considered as UCEC-specific lncRNAs and cancer hallmark-related genes if they were differentially expressed (P \u0026lt; 0.01). All the UCEC-specific lncRNAs and cancer hallmark-related genes were divided to up- and down-regulated lncRNAs and genes based on fold values.\u003c/li\u003e\n\u003c/ul\u003e\n\u003cp\u003e\u003cstrong\u003eConstruction of UCEC-specific cancer hallmark-related genes and lncRNAs co-expressed networks and identification of core modules\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eWe calculated Pearson\u0026rsquo;s correlation coefficients (PCCs) for each UCEC-specific cancer hallmark-related gene and lncRNA pair in all cancer hallmarks. The UCEC-specific gene and lncRNA pairs which their PCCs were higher than 0.3 or smaller than -0.3 and the p-values were smaller than 0.01 were extracted and considered as co-expressed UCEC-specific gene and lncRNA pairs. UCEC-specific cancer hallmark-related genes and lncRNAs co-expressed networks were constructed by Cytoscape 3.3.0 (http://www.cytoscape.org/) based on co-expressed pairs which their PCCs were higher than 0.5 or smaller than -0.5. In order to more accurately identify cancer hallmark-related genes and lncRNAs pairs in UCEC, we extracted core modules from co-expressed networks using MCODE module in cytoscape. Thus, certain numbers of core modules were identified for each cancer hallmark in UCEC.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eEvaluating hallmark-related risk for UCEC patients and dividing them to diverse risk groups based on core modules\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eFor each UCEC patient, risk score in each cancer hallmark was calculated and given. Multidimensional rank approach was applied to obtain risk score for each cancer hallmark based on genes and lncRNAs expression values in core modules. Up- and down-regulated genes and lncRNAs ranked in positive and negative order. 1000 permutation test was performed follow randomly perturbing cancer samples. The UCEC patient was considered as a risk sample for a specific cancer hallmark if permutation P values was smaller than 0.05 in any core module. Then, all the UCEC patients were divided to non-hallmark (0 risk hallmark), media-hallmark (1-5 risk hallmarks) and multi-hallmark groups (6-8 risk hallmarks).\u0026nbsp;\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eSurvival analysis for diverse risk cancer hallmark groups and functional analysis for UCEC-specific genes and lncRNAs\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eKaplan-Meier survival analysis was performed for diverse risk cancer hallmark groups. Log-rank test was using to asses statistical significance (P\u0026lt; 0.05). R 3.6.2 framework was performed for all analyses. Enrichr tools (http://amp.pharm.mssm.edu/Enrichr) with default parameters was used for functional analysis based on UCEC-specific genes in each cancer hallmark [\u003ca href=\"#_ENREF_17\"\u003e17\u003c/a\u003e]. Significantly enriched pathways (P\u0026lt; 0.05) were obtained for each cancer hallmark in UCEC.\u003c/p\u003e"},{"header":"Results","content":"\u003cp\u003e\u003cstrong\u003eSome lncRNAs and cancer hallmark-related genes were specific and co-expressed in UCEC\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eIn order to depict the roles of cancer hallmark in UCEC, we identified UCEC-specific lncRNAs and cancer hallmark-related genes by differential expression. 1880 (12.43%) UCEC-specific lncRNAs were identified between UCEC and control samples (Figure 2A). These UCEC-specific lncRNAs included 891 and 989 up- and down-regulated lncRNAs (Figure 2B). UCEC-specific \u0026nbsp;genes were also identified for each cancer hallmark (Figure 2C). In all kinds of cancer hallmarks, there were more than 50% UCEC-specific genes. Specially, there was 70.57% UCEC-specific genes in cancer hallmark genome instability and mutation. The results indicated that these cancer hallmark-related genes maybe serve as essential roles in UCEC. We assumed that lncRNAs and cancer hallmark-related genes could function by cooperating in UCEC. Thus, some co-expressed lncRNAs and cancer hallmark-related genes pairs were identified in UCEC for each kind of cancer hallmark. Most PCCs of co-expressed lncRNAs and cancer hallmark-related genes pairs were concentrated between 0.3 and 0.5 (Figure 2D). There were 8631 co-expressed lncRNAs and cancer hallmark-related genes pairs which their PCCs were higher than 0.5 (Figure 2E). In each kind of cancer hallmark, the numbers of pairs, lncRNAs and genes were diverse (Figure 2F). For example, there were 2525, 365 and 384 pairs, genes and lncRNAs in cancer hallmark self sufficiency in growth signals. However, only there were 84, 47 and 12 pairs, genes and lncRNAs in cancer hallmark reprogramming energy metabolism. The results indicated that diverse cancer hallmarks play different roles in UCEC. All above results indicated that lncRNAs and cancer hallmark-related genes cooperative pairs were important in UCEC.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eSome core modules were extracted from co-expressed lncRNAs and cancer hallmark-related genes networks in each cancer hallmark\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eFor each cancer hallmark, co-expressed lncRNAs and cancer hallmark-related genes which their PCCs were higher than 0.5 were extracted for constructing co-expressed networks. In cancer hallmark evading apoptosis, lncRNAs and cancer hallmark-related genes co-expressed network was constructed (Figure 3A). This co-expressed network contained 380 nodes (234 UCEC-specific lncRNAs and 147 cancer hallmark-related genes) and 906 edges. We found that some cancer hallmark-related genes and lncRNAs such as SLC25A27, AC005288 and GD5-AS1 played core roles in this co-expressed network. Most of cancer hallmark-related genes and lncRNAs showed positive correlations in UCEC. In eight cancer hallmarks, there were diverse numbers of core modules were identified (Figure 3B). Cancer hallmarks insensitivity to antigrowth signals, self sufficiency in growth signals and tissue invasion and metastasis had most core modules. The numbers of cancer hallmark-related genes and lncRNAs were also different (Figure 3C). For example, there were more lncRNA in core module 1 in cancer hallmark insensitivity to antigrowth signals. These core modules maybe show specific functions. For example, a core module in cancer hallmark evading apoptosis contained three lncRNAs and three genes (Figure 3D). PSMB8, PSMB10, PSME2 and PSMB8-AS1 were all proteasome-related genes or lncRNAs. The proteasome is a multicatalytic proteinase complex which is characterized by its ability to cleave peptides. Specially, cancer hallmark-related gene PSMB8 and lncRNA PSMB8-AS1 showed strong positive correlation (P\u0026lt; 0.001, PCC=0.79). These genes and lncRNAs in this core module showed close interactions. Another core module also showed close interactions (Figure 3E). All the results explained these core modules in co-expressed networks could function and serve as specific biomarks in UCEC.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eSpecific cancer hallmark-related risk were evaluated for each UCEC patient based on core modules\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eWe inferred that each UCEC patient maybe have diverse hallmark-related risk. Thus, we calculated risk scores for each UCEC patient based on core modules in each cancer hallmark. Only eight cancer hallmarks were extracted for calculating risk scores due to core modules. The density distribution of risk scores in all core modules of each cancer hallmark were similar (Figure 4A). Only a small number of UCEC patients showed higher risk scores. The differences of average risk scores for diverse top 20 core modules were also present (Figure 4B). These 20 core modules were significantly associated with more UCEC patients (Figure 4C). For example, there were 32.03% samples were significant in core module 1 in cancer hallmark insensitivity to antigrowth signals. This core module was a key module which had most related samples and also been explained in above results. We also discovered the percent of significant risk-related samples in 0%corresponding top ranked samples. 70% UCEC samples ranked before corresponding orders in most core modules (Figure 4D). These results indicated that UCEC patients showed differences of cancer hallmark-related risk scores.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eUCEC groups with diverse cancer hallmark risk showed specific features\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eEach UCEC patient could be associated with diverse cancer hallmark follow above pipeline. The numbers of UCEC in each cancer hallmark were different (Figure 5A). For example, there were more than 350 UCEC samples were associated with cancer hallmark tissue invasion and metastasis. The cancer hallmark reprogramming energy metabolism was related to 300 UCEC samples. Specially, some UCEC patients were associated with multiple cancer hallmarks. 13.59% UCEC patients had some relationships with any kinds of cancer hallmarks (Figure 5B). 11.55% UCEC patients were only related to one kind of cancer hallmark. Thus, we could divide all the UCEC patients to diverse groups with different numbers of cancer hallmarks (Figure 5C). The three diverse groups contained non-, media- and multi-hallmarks. 70.76% samples belonged to media-hallmarks groups. The UCEC patients with more cancer hallmarks usually had better survival days (Figure 5D). In addition, we also divided all the UCEC patients to another three cancer hallmark-related risk groups based on hierarchical clustering (Figure 5E). We also discovered that these diverse cancer hallmark-related risk groups showed different prognosis (Figure 5F, G). Group 2 significantly had better survival than group 1 and 3 (P= 0.014 and 0.023). All the results suggested that diverse hallmark-related risk groups had respective features and prognosis.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eSome key lncRNAs could participate in multiple kinds of cancer hallmarks and showed specific functions\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eIn order to further depict the roles of lncRNAs in cancer hallmarks for UCEC, we extracted some key lncRNAs which participate in multiple kinds of cancer hallmarks. There were 11 lncRNAs were associated with more than three kinds of cancer hallmarks (Figure 6A). lncRNAs AL590764.1, ANKRD10-IT1, NORD and AP000766.1 could participated in four kinds of cancer hallmarks (Figure 6B). However, the classes of these four kinds of cancer hallmarks were diverse. We inferred that these lncRNAs may serve as essential roles in UCEC. Thus, we further performed functional analysis for these lncRNAs in each kind of cancer hallmark. The lncRNAs were enrichment in some essential functions for cancer development (Figure 6C). For example, lncRNAs associated with evading apoptosis were enrichment in some hormone-related pathways such as negative regulation of intracellular estrogen receptor signaling, negative regulation of intracellular steroid hormone receptor signaling and regulation of intracellular estrogen receptor signaling. The estrogen receptor status is reported to be an important marker of UCEC [\u003ca href=\"#_ENREF_18\"\u003e18\u003c/a\u003e, \u003ca href=\"#_ENREF_19\"\u003e19\u003c/a\u003e]. In addition, we found most lncRNAs were associated with gap junction assembly pathway. Nishimura M et al. reported that gap junctional intercellular communication was suppressed via 5' CpG island methylation in promoter region of E-cadherin gene in UCEC cells [\u003ca href=\"#_ENREF_20\"\u003e20\u003c/a\u003e]. These results indicated that cancer hallmark-related lncRNAs could serve as essential roles in UCEC.\u003c/p\u003e"},{"header":"Discussion","content":"\u003cp\u003eIn this study, a calculated integrated approach which evaluate cancer hallmark-risk for each UCEC patient based on gene and lncRNA expression was developed. Some core modules were extracted from co-expressed\u0026nbsp;\u0026nbsp; UCEC-specific lncRNAs and genes networks in diverse cancer hallmark. All the UCEC patients were divided to diverse hallmark risk groups based on respective comprehensive risk scores. These diverse hallmark risk groups showed different survival. Some key lncRNAs were related to multiple cancer hallmarks and showed specific functions.\u003c/p\u003e\n\u003cp\u003e\u0026nbsp; Dividing cancer samples to groups with diverse molecular characteristics could provide assistance for cancer diagnosis and individualized treatment. Wright GW et al. described an algorithm that determines the probability that a patient's lymphoma belongs to one of seven genetic subtypes based on its genetic features [\u003ca href=\"#_ENREF_21\"\u003e21\u003c/a\u003e]. Lee S et al. presented an explainable deep learning model with attention mechanism and network propagation for cancer subtype classification [\u003ca href=\"#_ENREF_22\"\u003e22\u003c/a\u003e]. Most of these developed clarification methods only focused on coding genes. In our work, lncRNAs were considered as essential and key factors for constructing the computational approach about cancer patients classification. The computational approach not only divided UCEC patients to diverse hallmark risk groups but also discovered some key lncRNAs which serve as important roles in UCEC. These lncRNAs were significantly enriched in some essential cancer development functions.\u003c/p\u003e\n\u003cp\u003e\u0026nbsp; Analysis of cancer hallmarks could greatly improve our understanding of the occurrence, development and metastasis of many cancer types. In present study, we applied cancer hallmarks to evaluate and distinguish UCEC patients. We discovered that most cancer hallmark-related genes were differentially expressed in UCEC. It indicated that cancer hallmark-related genes maybe had essential functions in UCEC. Specially, we found that each UCEC patient showed obvious differences for risk of cancer hallmarks. Although cancer hallmarks were considered as common properties for cancer, each UCEC patient still showed personalized features. This result supported that there were great differences among different cancer patients and the application of personalized medicine in cancer treatment is very important. Thus, our study could help to evaluate risk for UCEC patient and establishes personalized therapy.\u003c/p\u003e"},{"header":"Conclusions","content":"\u003cp\u003eIn summary, the present study developed a standardized computational procedures to evaluate cancer patient hallmark-related risk based on genes and lncRNAs expression. All the UCEC patients were divided to diverse cancer hallmark-related risk groups and showed different prognosis. Some key cancer hallmark-related lncRNAs were significantly associated with UCEC development. Collectively, our study leads to a novel starting point for future functional explorations, the identification of biomarkers, and lncRNA-based targeted therapy for UCEC.\u003c/p\u003e"},{"header":"Abbreviations","content":"\u003cp\u003eUCEC: uterine corpus endometrioid carcinoma; lncRNA: long non-coding RNA; TCGA: The Cancer Genome Atlas; GO: gene ontology; PCC: Pearson correlation coefficients.\u003c/p\u003e"},{"header":"Declarations","content":"\u003cp\u003e\u003cstrong\u003eFunding\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eNot applicable.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eDisclosure of interest\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe authors declare that they have no competing interest.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAcknowledgments\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eNot applicable.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAvailability of data and material\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe data that support the findings of this study are available from the corresponding author upon reasonable request.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eEthics approval and consent to participate\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eNot applicable.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003ePatient consent for publication\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eNot applicable.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eCompeting interests\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe authors declare that they have no competing interests.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eConsent for publication\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eNot applicable.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAuthors' contributions\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eLN conceived and designed the experiments, ZSY, LY, LTY and CXD analysed the data, and LN and ZSY wrote the manuscript.\u003c/p\u003e"},{"header":"References","content":"\u003col\u003e\n\u003cli\u003eBray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A: \u003cstrong\u003eGlobal cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries.\u003c/strong\u003e \u003cem\u003eCA Cancer J Clin \u003c/em\u003e2018, \u003cstrong\u003e68:\u003c/strong\u003e394-424.\u003c/li\u003e\n\u003cli\u003eGaber C, Meza R, Ruterbusch JJ, Cote ML: \u003cstrong\u003eEndometrial Cancer Trends by Race and Histology in the USA: Projecting the Number of New Cases from 2015 to 2040.\u003c/strong\u003e \u003cem\u003eJ Racial Ethn Health Disparities \u003c/em\u003e2016.\u003c/li\u003e\n\u003cli\u003eSiegel RL, Miller KD, Jemal A: \u003cstrong\u003eCancer statistics, 2020.\u003c/strong\u003e \u003cem\u003eCA Cancer J Clin \u003c/em\u003e2020, \u003cstrong\u003e70:\u003c/strong\u003e7-30.\u003c/li\u003e\n\u003cli\u003eBell DW, Ellenson LH: \u003cstrong\u003eMolecular Genetics of Endometrial Carcinoma.\u003c/strong\u003e \u003cem\u003eAnnu Rev Pathol \u003c/em\u003e2019, \u003cstrong\u003e14:\u003c/strong\u003e339-367.\u003c/li\u003e\n\u003cli\u003eAnastasiadou E, Jacob LS, Slack FJ: \u003cstrong\u003eNon-coding RNA networks in cancer.\u003c/strong\u003e \u003cem\u003eNat Rev Cancer \u003c/em\u003e2018, \u003cstrong\u003e18:\u003c/strong\u003e5-18.\u003c/li\u003e\n\u003cli\u003eKaikkonen MU, Adelman K: \u003cstrong\u003eEmerging Roles of Non-Coding RNA Transcription.\u003c/strong\u003e \u003cem\u003eTrends Biochem Sci \u003c/em\u003e2018, \u003cstrong\u003e43:\u003c/strong\u003e654-667.\u003c/li\u003e\n\u003cli\u003eRansohoff JD, Wei Y, Khavari PA: \u003cstrong\u003eThe functions and unique features of long intergenic non-coding RNA.\u003c/strong\u003e \u003cem\u003eNat Rev Mol Cell Biol \u003c/em\u003e2018, \u003cstrong\u003e19:\u003c/strong\u003e143-157.\u003c/li\u003e\n\u003cli\u003eUszczynska-Ratajczak B, Lagarde J, Frankish A, Guigo R, Johnson R: \u003cstrong\u003eTowards a complete map of the human long non-coding RNA transcriptome.\u003c/strong\u003e \u003cem\u003eNat Rev Genet \u003c/em\u003e2018, \u003cstrong\u003e19:\u003c/strong\u003e535-548.\u003c/li\u003e\n\u003cli\u003eWang Y, Fang Z, Hong M, Yang D, Xie W: \u003cstrong\u003eLong-noncoding RNAs (lncRNAs) in drug metabolism and disposition, implications in cancer chemo-resistance.\u003c/strong\u003e \u003cem\u003eActa Pharm Sin B \u003c/em\u003e2020, \u003cstrong\u003e10:\u003c/strong\u003e105-112.\u003c/li\u003e\n\u003cli\u003eWu P, Mo Y, Peng M, Tang T, Zhong Y, Deng X, Xiong F, Guo C, Wu X, Li Y, et al: \u003cstrong\u003eEmerging role of tumor-related functional peptides encoded by lncRNA and circRNA.\u003c/strong\u003e \u003cem\u003eMol Cancer \u003c/em\u003e2020, \u003cstrong\u003e19:\u003c/strong\u003e22.\u003c/li\u003e\n\u003cli\u003eZhang G, Ma A, Jin Y, Pan G, Wang C: \u003cstrong\u003eLncRNA SNHG16 induced by TFAP2A modulates glycolysis and proliferation of endometrial carcinoma through miR-490-3p/HK2 axis.\u003c/strong\u003e \u003cem\u003eAm J Transl Res \u003c/em\u003e2019, \u003cstrong\u003e11:\u003c/strong\u003e7137-7145.\u003c/li\u003e\n\u003cli\u003eWang E, Zaman N, McGee S, Milanese JS, Masoudi-Nejad A, O'Connor-McCourt M: \u003cstrong\u003ePredictive genomics: a cancer hallmark network framework for predicting tumor clinical phenotypes using genome sequencing data.\u003c/strong\u003e \u003cem\u003eSemin Cancer Biol \u003c/em\u003e2015, \u003cstrong\u003e30:\u003c/strong\u003e4-12.\u003c/li\u003e\n\u003cli\u003eHanahan D, Weinberg RA: \u003cstrong\u003eThe hallmarks of cancer.\u003c/strong\u003e \u003cem\u003eCell \u003c/em\u003e2000, \u003cstrong\u003e100:\u003c/strong\u003e57-70.\u003c/li\u003e\n\u003cli\u003eHanahan D, Weinberg RA: \u003cstrong\u003eHallmarks of cancer: the next generation.\u003c/strong\u003e \u003cem\u003eCell \u003c/em\u003e2011, \u003cstrong\u003e144:\u003c/strong\u003e646-674.\u003c/li\u003e\n\u003cli\u003ePlaisier CL, Pan M, Baliga NS: \u003cstrong\u003eA miRNA-regulatory network explains how dysregulated miRNAs perturb oncogenic processes across diverse cancers.\u003c/strong\u003e \u003cem\u003eGenome Res \u003c/em\u003e2012, \u003cstrong\u003e22:\u003c/strong\u003e2302-2314.\u003c/li\u003e\n\u003cli\u003eCarbon S, Ireland A, Mungall CJ, Shu S, Marshall B, Lewis S, Ami GOH, Web Presence Working G: \u003cstrong\u003eAmiGO: online access to ontology and annotation data.\u003c/strong\u003e \u003cem\u003eBioinformatics \u003c/em\u003e2009, \u003cstrong\u003e25:\u003c/strong\u003e288-289.\u003c/li\u003e\n\u003cli\u003eKuleshov MV, Jones MR, Rouillard AD, Fernandez NF, Duan Q, Wang Z, Koplev S, Jenkins SL, Jagodnik KM, Lachmann A, et al: \u003cstrong\u003eEnrichr: a comprehensive gene set enrichment analysis web server 2016 update.\u003c/strong\u003e \u003cem\u003eNucleic Acids Res \u003c/em\u003e2016, \u003cstrong\u003e44:\u003c/strong\u003eW90-97.\u003c/li\u003e\n\u003cli\u003eSho T, Hachisuga T, Nguyen TT, Urabe R, Kurita T, Kagami S, Kawagoe T, Matsuura Y, Shimajiri S: \u003cstrong\u003eExpression of estrogen receptor-alpha as a prognostic factor in patients with uterine serous carcinoma.\u003c/strong\u003e \u003cem\u003eInt J Gynecol Cancer \u003c/em\u003e2014, \u003cstrong\u003e24:\u003c/strong\u003e102-106.\u003c/li\u003e\n\u003cli\u003eTogami S, Sasajima Y, Oi T, Ishikawa M, Onda T, Ikeda S, Kato T, Tsuda H, Kasamatsu T: \u003cstrong\u003eClinicopathological and prognostic impact of human epidermal growth factor receptor type 2 (HER2) and hormone receptor expression in uterine papillary serous carcinoma.\u003c/strong\u003e \u003cem\u003eCancer Sci \u003c/em\u003e2012, \u003cstrong\u003e103:\u003c/strong\u003e926-932.\u003c/li\u003e\n\u003cli\u003eNishimura M, Saito T, Yamasaki H, Kudo R: \u003cstrong\u003eSuppression of gap junctional intercellular communication via 5' CpG island methylation in promoter region of E-cadherin gene in endometrial cancer cells.\u003c/strong\u003e \u003cem\u003eCarcinogenesis \u003c/em\u003e2003, \u003cstrong\u003e24:\u003c/strong\u003e1615-1623.\u003c/li\u003e\n\u003cli\u003eWright GW, Huang DW, Phelan JD, Coulibaly ZA, Roulland S, Young RM, Wang JQ, Schmitz R, Morin RD, Tang J, et al: \u003cstrong\u003eA Probabilistic Classification Tool for Genetic Subtypes of Diffuse Large B Cell Lymphoma with Therapeutic Implications.\u003c/strong\u003e \u003cem\u003eCancer Cell \u003c/em\u003e2020, \u003cstrong\u003e37:\u003c/strong\u003e551-568 e514.\u003c/li\u003e\n\u003cli\u003eLee S, Lim S, Lee T, Sung I, Kim S: \u003cstrong\u003eCancer subtype classification and modeling by pathway attention and propagation.\u003c/strong\u003e \u003cem\u003eBioinformatics \u003c/em\u003e2020.\u003c/li\u003e\n\u003c/ol\u003e"}],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":true,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":false,"hideJournal":true,"highlight":"","institution":"","isAcceptedByJournal":false,"isAuthorSuppliedPdf":false,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":false,"isPdf":false,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"[email protected]","identity":"researchsquare","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":true,"externalIdentity":"","sideBox":"","snPcode":"","submissionUrl":"/submission","title":"Research Square","twitterHandle":"researchsquare","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"","reportingPortfolio":"","inReviewEnabled":false,"inReviewRevisionsEnabled":true},"keywords":"Long non-coding RNA, cancer hallmark, uterine corpus endometrioid carcinoma, survival, molecular subtypes","lastPublishedDoi":"10.21203/rs.3.rs-37407/v1","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-37407/v1","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"\u003cp\u003e\u003cstrong\u003eBackground: \u003c/strong\u003eUterine corpus endometrioid carcinoma (UCEC), a common gynecological malignancy with high incidence, affects the mental and physical health of women. It is increasingly evident that long non-coding RNAs (lncRNAs) have causative roles in cancers including UCEC. However, very little is known about the cancer hallmark-related risks for UCEC based on lncRNAs and genes. \u003c/p\u003e\u003cp\u003e\u003cstrong\u003eMethods: \u003c/strong\u003eIn this study, a computational integrated pipeline was development to evaluate cancer hallmark-related risk for each UCEC patient based on gene and lncRNA expression. \u003c/p\u003e\u003cp\u003e\u003cstrong\u003eResults:\u003c/strong\u003e Some UCEC-specific cancer hallmark-related genes and lncRNAs were identified. Core modules were extracted from co-expressed UCEC-specific lncRNAs and genes networks for each cancer hallmark. Some core modules showed specific features and functions. Follow a multi-dimensional rank approach and core modules, each UCEC sample was given an integrated cancer hallmark-related risk score. We divided all the UCEC patients to diverse hallmark risk groups including multi-hallmark, media-hallmark and single-hallmark groups and they showed different prognosis. We also identified some key lncRNAs which participated in multiple kinds of cancer hallmarks. These key lncRNAs were associated with some essential pathways such as estrogen receptor signaling. \u003c/p\u003e\u003cp\u003e\u003cstrong\u003eConclusions:\u003c/strong\u003e In conclusion, the present study provide novel insights to the classification and treatment of UCEC.\u0026nbsp;\u003c/p\u003e","manuscriptTitle":"Inference of Molecular Subtypes of Uterine Corpus Endometrioid Carcinoma with Different Survival Based on Cancer Hallmarks-associated Long Non-Coding RNAs and Gene Modules","msid":"","msnumber":"","nonDraftVersions":[{"code":1,"date":"2020-06-30 19:41:57","doi":"10.21203/rs.3.rs-37407/v1","editorialEvents":[{"type":"communityComments","content":0}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"researchsquare","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":true,"externalIdentity":"","sideBox":"","snPcode":"","submissionUrl":"/submission","title":"Research Square","twitterHandle":"researchsquare","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"","reportingPortfolio":"","inReviewEnabled":false,"inReviewRevisionsEnabled":true}}],"origin":"","ownerIdentity":"549b73e2-efc6-4c5a-96a3-1b96f95f5550","owner":[],"postedDate":"June 30th, 2020","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"posted","subjectAreas":[{"id":136673,"name":"Endocrinology \u0026 Metabolism"}],"tags":[],"updatedAt":"2020-07-20T19:14:31+00:00","versionOfRecord":[],"versionCreatedAt":"2020-06-30 19:41:57","video":"","vorDoi":"","vorDoiUrl":"","workflowStages":[]},"version":"v1","identity":"rs-37407","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-37407","identity":"rs-37407","version":["v1"]},"buildId":"J0_U0BvcaRcwD8yVFaRlm","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}