Next Article in Journal
Implementation of External Beam Five-Fraction Adjuvant Breast Irradiation in a US Center
Previous Article in Journal
CDK6 Degradation Is Counteracted by p16INK4A and p18INK4C in AML
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Cancers
Volume 14
Issue 6
10.3390/cancers14061557
cancers-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Long Non-Coding RNA-Based Functional Prediction Reveals Novel Targets in Notch-Upregulated Ovarian Cancer

by
Seonhyang Jeong
1,†,
Sunmi Park
1,†,
Young Suk Jo
1,
Moon Jung Choi
2,
Gibbeum Lee
2,
Seul Gi Lee
3,
Min Chul Choi
4,
Hyun Park
4,
Won Duk Joo
4,
Sang Geun Jung
4,* and
Jandee Lee
2,*
1
Department of Internal Medicine, Yonsei University College of Medicine, 50-1 Yonsei-ro, Seodaemun-gu, Seoul 03722, Korea
2
Department of Surgery, Open NBI Convergence Technology Research Laboratory, Severance Hospital, Yonsei Cancer Center, Yonsei University College of Medicine, 50-1 Yonsei-ro, Seodaemun-gu, Seoul 03722, Korea
3
Department of Surgery, Eulji University School of Medicine, 95 Dunsanseo-ro, Seo-gu, Daejeon 35233, Korea
4
Department of Gynecological Oncology, Comprehensive Gynecologic Cancer Center, CHA University, 59 Yatap-ro, Bundang-gu, Seongnam 13496, Korea
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Cancers 2022, 14(6), 1557; https://doi.org/10.3390/cancers14061557
Submission received: 22 February 2022 / Revised: 14 March 2022 / Accepted: 15 March 2022 / Published: 18 March 2022
(This article belongs to the Topic Advances in Ovarian Cancer Research: From Biology to Therapeutics)
Browse Figures
Review Reports Versions Notes

Abstract

:

Simple Summary

The study contributes to our understanding of the role of lncRNAs in the regulation of Notch signaling and their target genes in ovarian cancer. These lncRNAs were identified with in silico analyses using The Cancer Genome Atlas and the results were validated using our transcriptome data from the in vitro NOTCH1/3 silencing experiments and the qRT-PCR analysis, using tissue samples from patients with ovarian tumors. Finally, to improve our understanding of the complexity of Notch signaling, we identified master transcriptional factors which might orchestrate the upregulation of Notch-related lncRNAs and coding-genes in ovarian cancer.

Abstract

Notch signaling is a druggable target in high-grade serous ovarian cancers; however, its complexity is not clearly understood. Recent revelations of the biological roles of lncRNAs have led to an increased interest in the oncogenic action of lncRNAs in various cancers. In this study, we performed in silico analyses using The Cancer Genome Atlas data to discover novel Notch-related lncRNAs and validated our transcriptome data via NOTCH1/3 silencing in serous ovarian cancer cells. The expression of novel Notch-related lncRNAs was down-regulated by a Notch inhibitor and was upregulated in high-grade serous ovarian cancers, compared to benign or borderline ovarian tumors. Functionally, Notch-related lncRNAs were tightly linked to Notch-related changes in diverse gene expressions. Notably, genes related to DNA repair and spermatogenesis showed specific correlations with Notch-related lncRNAs. Master transcription factors, including EGR1, CTCF, GABPα, and E2F4 might orchestrate the upregulation of Notch-related lncRNAs, along with the associated genes. The discovery of Notch-related lncRNAs significantly contributes to our understanding of the complex crosstalk of Notch signaling with other oncogenic pathways at the transcriptional level.

1. Introduction

Ovarian cancer is one of the common types of malignancies in women and is often fatal due to the unavailability of precise early-detection tests, as well as its aggressive features, such as an early progression to metastatic disease [1,2,3,4,5]. More than 60% of patients first diagnosed with ovarian cancer have stage-III or stage-IV disease, indicating that it has spread into the abdominal cavity and other peritoneal structures [6,7,8]. Histologic subtypes of ovarian cancer can also affect clinical outcomes, along with the stage. In this regard, a high-grade serous carcinoma (HGSC) is considered, due to structural and ultrastructural evidence, as the most aggressive subtype of ovarian cancer, with a debated origin, a frequent resistance to chemotherapy, and a poor prognosis [9,10,11,12]. To overcome these diagnostic and therapeutic limitations, clinicians and researchers have investigated the biological and clinical features of ovarian cancer, especially HGSC [10,13].
About 98% of the human genome is composed of protein non-coding regions, being transcribed into non-coding RNAs [14]. Among the non-coding RNAs, long non-coding RNAs (lncRNAs) are defined as transcripts longer than 200 nucleotides. According to the Functional Annotation of the Mouse/Mammalian Genomes (FANTOM), 27,919 lncRNAs have been identified in humans [15,16]. Moreover, recent transcriptomic analyses by next generation sequencing have indicated that the number of lncRNAs in mammals could be in the order of tens of thousands [17]. Due to the rapid development in next generation sequencing methods, various lncRNAs, such as ANRIL, CCAT1, FAL1, H19, MALAT1, and HOTAIR, have been discovered in ovarian cancer [18]. HOTAIR has been recently discovered to enhance the invasion and migration of ovarian cancer progression by regulating EZH2, as well as regulating the expression of miR-193a and ODK2, via introducing H3K27me3 [19]. However, despite accumulating evidence suggesting that the majority of lncRNAs are biologically relevant, the insight into their functions in cell biology is still limited.
Over the past two decades, the major components of Notch signaling pathways, including Notch receptors 1 and 3 (NOTCH1 and NOTCH3); their ligands, such as jagged canonical Notch ligand 1 (JAG1) and delta like canonical Notch ligand 4 (DLL4); and downstream proteins, such as the Hes family bHLH transcription factor 1 (HES1) and DLG associated protein 5 (DLGAP5), have been intensively studied in ovarian cancers [20,21,22]. The Notch signaling pathway has important roles in angiogenesis and chemo-resistance. Additionally, Notch signaling inhibitors, such as γ-secretase inhibitors (GSIs) and DLL4 inhibitors, as well as Notch antibodies, have been examined in preclinical and clinical settings [23,24]. However, more work is needed to understand the diverse cell signaling pathways and cell biological processes that are affected by Notch activation in order to develop more specific and effective means of inhibiting Notch pathway components [25]. In this regard, Notch-related lncRNAs have been described in other cancers [26].
In the present study, we identified potential Notch-related lncRNA clusters in ovarian cancer (OVCA) using data from The Cancer Genome Atlas (TCGA) and our transcriptome analysis [27,28,29]. We also compared the cell biological functions of Notch-related mRNAs proposed by a gene set enrichment analysis (GSEA) and that of Notch-related lncRNAs identified by the in silico prediction of lncRNA function (FuncPred) [17,30]. In addition, to understand the regulatory mechanism of Notch-related lncRNAs, we examined master transcriptional factors (TFs) using the Encyclopedia of DNA Elements (ENCODE) [31]. This study provides a fresh insight into how Notch-related lncRNAs could cooperate with Notch-related gene transcripts through novel signaling mechanisms.

2. Materials and Methods

2.1. Public Databases

Two representative public datasets were utilized in this study. First, the mRNA expression dataset from TCGA (https://portal.gdc.cancer.gov/ (access on 1 November 2021)) OVCA RNA-seq database (http://cancergenome.nih.gov/ (access on 1 November 2021), version 1.0.6) was collected and analyzed using GSEA, with hallmark gene sets of the Molecular Signatures Database (MSigDB), gene sets derived from the Kyoto Encyclopedia of Genes and Genomes (KEGG), and gene sets related to TF targets based on sharing upstream cis-regulatory motifs as potential TF binding sites [30,32]. For a large number of gene sets, nominal p-values were considered to indicate a significant enrichment of a set if its normalized enrichment scores (NES) had a False Discovery Rate (FDR) q-value below 0.25. Second, to identify Notch-related lncRNA, its expression dataset was collected from The Atlas of Non-coding RNAs in Cancer (TANRIC) and was analyzed for its functional aspects using FuncPred (http://www.funcpred.com/ (access on 17 August 2021) [17,29]. The information on gene regulation, including the binding site and TFs of Notch-related lncRNA, was investigated with ENCODE data using the Ensembl Genome Browser. The ENCODE analysis, based on chromatin immunoprecipitation sequencing data, was used to identify the regulatory elements for each gene; scores were presented in the ENCODE table. The higher the score, the stronger the binding of the relevant motif, and the greater likelihood of it being identified in all cell lines [31]. The overall schematic flow of our research strategy is presented in Figure S1.

2.2. Patients and Specimens

Human tissue samples were obtained from 100 female patients (mean age = 54.4 ± 9.1) following a unilateral or bilateral salpingo-oophorectomy for benign or borderline tumors (n = 16), clear cell carcinoma (n = 5), or HGSC (n = 79) between October 2018 and May 2020 at Bundang CHA Medical Center (CHA University, Gyeonggi-do, Korea). All samples were frozen in liquid nitrogen and stored at −80 °C until use. This study was approved by the institutional review board of Bundang CHA Medical Center (CHAMC 2018-10-008, 10 October 2018) and was conducted in accordance with the recommendations of the institutional review board. Informed consent was obtained from all individual participants included in this study.

2.3. RNA Isolation and Real-Time PCR

Total RNA was isolated from frozen tissues using a Trizol reagent (Invitrogen, Carlsbad, CA, USA), and the RNA quality was assessed using a 2100 Bioanalyzer System (Agilent Technologies, Santa Clara, CA, USA). Complementary DNA (cDNA) was prepared from the total RNA using the SuperScript™ III First-Strand Synthesis System (Invitrogen, 18080051). Quantitative real-time PCR (qRT-PCR) was performed using a QuantiTect SYBR Green RT-PCR kit (Qiagen, Valencia, CA, USA). The relative expression was calculated using the Power SYBR™ Green PCR Master Mix (Applied Biosystems, Foster City, CA, USA). Primers used in qRT-PCR are listed in Table S1. All experiments were repeated at least three times.

2.4. Cell Culture and Reagents

The normal ovarian surface epithelium, the HOSEpic cell line (cat#7310), was purchased from ScienCell (Carlsbad, CA, USA). The clear cell line, ES-2 (cat#CRL-1978) and TOV21G (cat#CRL-11730) cell lines, were purchased from ATCC (Manassas, VA, USA). High-grade serous carcinoma (HGSC), Caov3 (cat#30075), OVCAR3 (cat#30161), and SKOV3 (cat#30077) cell lines were purchased from the Korean Cell Line Bank (Seoul, Korea). The HOSEpic cell line was cultured in an ovarian epithelial cell medium (cat#7311, ScienCell) supplied with an ovarian epithelial cell growth supplement (cat#7352, ScienCell) and an antibiotic solution (cat#0503, ScienCell). All cell lines have been authenticated using STR profiling. The ES-2, TOV21G, OVCAR3, and SKOV3 cell lines were cultured in RPMI (cat# 11875093, Thermo Fisher Scientific, Waltham, MA, USA), supplemented with 10% fetal bovine serum (cat# 35015CV, Thermo Fisher Scientific), 1% penicillin, and streptomycin. The GSI, DAPT (N-[N-(3,5-difluorophenacetyl)-L-alanyl]-S-phenylglycine t-butyl ester), was purchased from Sigma-Aldrich (cat# D5942, Burlington, MA, USA). For DAPT treatment, cells were cultured for 48 h at the concentration indicated in the figures.

2.5. siRNA Transfection and Western Blot Analysis

Three different small interfering RNAs (siRNAs) targeting NOTCH1 (cat#4851-1, 4851-2, and 4851-3) and NOTCH3 (cat#4854-1, 4854-2, and 4854-3) were purchased from Bioneer (Deajeon, Korea). Non-targeting siRNA was used as a negative control. The RNAi oligonucleotide or RNAi negative control was transfected into the cells using a Lipofectamine RNAiMAX Transfection Reagent (Thermo Fisher Scientific, Waltham, MA, USA) according to the manufacturer’s instructions. Cells were lysed with a cell extraction buffer (Thermo Fisher Scientific) containing protease inhibitors. Protein concentrations were determined using the BCA protein assay (Thermo Fisher Scientific). Protein samples were separated using SDS-PAGE and transferred to polyvinylidene difluoride membranes. The membranes were incubated with the indicated primary antibodies overnight at 4 °C. The following antibodies were used: NOTCH1 (cat#4380, Cell Signaling Technology, Danvers, MA, USA), NOTCH3 (cat#5276, Cell Signaling Technology), and α-tubulin (cat#sc-8035, Santa Cruz Biotechnology, Inc., Dallas, TX, USA). The blots were visualized using the Pierce™ ECL Western Blotting Substrate (Thermo Fisher Scientific, Waltham, MA, USA). The original WB can be found at Figure S5.

2.6. Total RNA Sequencing

Total RNA was isolated from OVCAR3 cells using Trizol (Thermo Fisher Scientific). After confirming RNA quality, the TruSeq Stranded mRNA LT Sample Prep Kit (Illumina, San Diego, CA, USA) was used to prepare a cDNA library according to the Sample Preparation Guide (Part #15031048 Rev E) by Macrogen Inc (Seoul, Korea). The 100-nt paired-end sequencing was conducted using NovaSeq 6000 (Illumina, San Diego, CA, USA). To remove any bias in the analysis, low-quality RNA or artifacts, such as adaptor sequences, contaminant DNA, and PCR duplicates were trimmed from raw data. Trimmed data was mapped to the reference genome using HISAT2 (version 2.1.0) and the aligned reads were generated. StringTie (version 2.1.3b) was used to produce reconstructions of genes and estimated expression levels using the reference-based aligned reads information. In all, 60,606 transcripts were classified based on their overlap with the GENCODE annotation (version 38).

2.7. Statistical Analysis

Heatmaps and correlation graphs were drawn using R studio (version 3.3.3). Continuous variables were compared using Student’s t-tests or Mann–Whitney U tests, accordingly. Categorical variables were compared using either χ2 tests or Fisher’s exact tests, accordingly. Pearson’s correlation coefficient was used to examine the relationship between lncRNAs with the genes of interest. Statistical analyses were performed using Microsoft Excel (ver. 2016, Microsoft, Redmond, WA, USA), SPSS version 25 for Windows (IBM Corporation, Armonk, NY, USA) and GraphPad Prism (GraphPad Software, Inc., San Diego, CA, USA). Data are presented as the means ± SD, and all p-values are two-sided. p < 0.05 was considered significant.

3. Results

3.1. NOTCH1/3 and Their Target Gene Expression in TCGA HGS-OVCA

To discover potential Notch-related lncRNAs, we used the public data from TCGA HGS-OVCA (n = 303) and analyzed it according to the NOTCH1/3 expression status. First, we classified the dataset into a high group with a high expression of both NOTCH1 and NOTCH3 (n = 97), and a low group with a low expression of both NOTCH1 and NOTCH3 (n = 96). The criteria for high and low expression levels were based on median transcripts per million (TPM)-mapped read values. To validate whether these group classifications reflected the activity of Notch signaling, we generated a Notch target gene set (n = 36) based on a literature review and performed GSEA using our homemade gene set [24]. As shown in Figure S2A, the high NOTCH1/3 group showed a coordinated enrichment of Notch target genes (p < 0.001). Following the validation of our group classification, we analyzed the expression of the entire gene set and generated a heatmap using 100 representative genes showing the greatest difference in expression between the two groups (Figure S2B). This heatmap clearly suggests that Notch signaling contributes to the up- or down-regulation of diverse genes. To understand the cell biological changes by Notch signaling, we performed GSEA using KEGG and Hallmark gene sets (Figure S2C and Figure 1A). As expected, a gene set related to Notch signaling was coordinately enriched in the high NOTCH1/3 group, as shown by the analyses of both gene sets (Figure S2C and Figure 1A). Besides Notch signaling, gene sets related to WNT signaling, epithelial mesenchymal transition (EMT), and Hedgehog signaling were highly enriched in the high NOTCH1/3 group (Figure 1A). The results suggested that Notch signaling could affect invasive features of HGSC by orchestrating EMT-related gene expression; our classification appropriately reflected the activity of Notch-signaling.

3.2. Discovery of Notch-Related lncRNAs Using TANRIC Database

Ten of the available TCGA HGS-OVCA cases (n = 303) did not have data related to the lncRNA expression status, and the remaining 293 cases were analyzed to match mRNA and lncRNA expression levels. First, we divided the HGS-OVCA cases into four groups: I (n = 92), the low expression of both NOTCH-1 and -3; II (n = 54), the high expression of NOTCH3 and the low expression of NOTCH1; III (n = 54), the high expression of NOTCH1 and the low expression of NOTCH3; IV (n = 93), the high expressions of both NOTCH1 and NOTCH3. Next, based on this grouping, we analyzed the expression status of lncRNAs and generated a heatmap using representative lncRNAs (n = 1737), showing a differential gene expression pattern among the four groups (Figure 1B). To investigate the biological impact of the discovered lncRNAs in HGS-OVCA, we used FuncPred, a web-based in silico predictor of lncRNA functions. Applying the top 10% highly expressed lncRNAs in Group IV to FuncPred, we found 46 related Hallmark gene sets (Figure 1C). Interestingly, every 13th gene set significantly enriched in the high NOTCH1/3 expression group of our GSEA overlapped with the lncRNA-related gene sets in FuncPred (Figure 1D). This data suggested that Notch-related lncRNA may cooperate with Notch-related mRNA to generate the invasiveness of HGSC; again, our grouping and analytic method appropriately reflected Notch signaling activation. Among the 33 gene sets discovered only in FuncPred, the top-ranked gene sets related to DNA repair and spermatogenesis were not significantly enriched in the high NOTCH1/3 expression group by GSEA using mRNA data from TCGA OVCA. Taken together, these data suggested that Notch-related lncRNAs could be involved in complex biological functions in Notch up-regulated ovarian cancer, along with Notch-related mRNA, and the functional prediction of lncRNAs might be a powerful tool to investigate NOTCH1/3 functions in HGSC.

3.3. Validation of Notch-Related lncRNAs by NOTCH1/3 Silencing and Analysis of Clinical Tissue Samples

To validate our FuncPred analysis results, we investigated the effect of NOTCH1/3 silencing using siRNAs. First, we performed a Western blot analysis to check the NOTCH1/3 expression status in various normal ovarian cell lines and ovarian cancer cell lines (Figure S3A). Among the cell lines tested, TOV21G and OVCAR3 showed simultaneous NOTCH1 and NOTCH3 expressions. Because TCGA OVCA includes data from patients with HGSC but not from clear cell carcinomas, we selected the HGSC cell line, OVCAR3, for this study. We tested three types of siNOTCH1 and siNOTCH3 (Figure S3B,C, respectively). These efficiently reduced transmembrane NOTCH1 or NOTCH3, as analyzed by Western blotting (Figure S3B,C). The total RNAs extracted from OVCAR3 treated with two kinds of siRNA combinations were then subjected to total RNA sequencing. In all, 550 differentially expressed lncRNAs were identified (Figure 2A). From these, we selected 216 lncRNAs that were downregulated by NOTCH1/3 silencing to investigate the related gene sets using FuncPred. The lncRNA-related gene sets obtained from the FuncPred analysis (Figure 2B) are similar to those from our TANRIC analysis (Figure 1C). We found 43 Notch-related lncRNAs in common between the 1670 up-regulated lncRNAs in high NOTCH1/3 OVCA and the 216 lncRNAs that were down-regulated by NOTCH1/3 silencing (Figure 2C,D). The FuncPred analysis, using these 43 lncRNAs, revealed that genes related to spermatogenesis and DNA repair were consistently top ranked among the related genes (Figure 2E).
To validate our transcriptomic data, we silenced NOTCH1/3 in OVCAR3 cells using different combinations of siRNAs specific for NOTCH1 and NOTCH3 (Figure S4A). Using total RNAs from these cells, we detected a decreased NOTCH1/3 expression, along with their well-known target genes, such as the hairy and enhancer-of-split-1 (HES1) and the hairy/enhancer-of-split related to YRPW motif protein 1 (HEY1) (Figure S4A,B). The levels of representative common lncRNAs, such as AP000525.1, ASB16-AS1, PRKCZ-AS, GUSBP11, and AL355488.1 were also decreased by different combinations of siRNAs (Figure 3A). In addition, genes related to DNA repair and spermatogenesis were also decreased by these siRNAs (Figure 3B,C). Mature Notch receptors are processed and assembled as heterodimeric proteins [33]. The binding of Notch receptors to ligands triggers the dissociation of heterodimeric proteins, which causes the cleavage of the receptors by γ-secretase. These result in the release of the Notch intracellular domain (NICD), which, in turn, translocates to the nucleus and activates the transcription of Notch targets [34,35]. DAPT is a potent γ-secretase inhibitor (GSI) that targets Notch signaling. DAPT diminishes the breakdown product transmembrane/intracellular region which reflects the activity of Notch in a dose-dependent manner (Figure S4C). Treating OVCAR3 with DAPT, we observed a decreased expression of NOTCH1/3, along with their target genes HES1 and HEY1 (Figure S4D,E). Representative Notch-related lncRNAs and genes related to DNA repair and spermatogenesis were also down-regulated in OVCAR3 following DAPT treatment (Figure 3D–F).
To validate these results in clinical tissue samples, we extracted RNAs from HGSC (n = 79) and divided them into two groups: the high NOTCH1/3 expression group (n = 21) and the low NOTCH1/3 expression group (n = 19) (Figure 4A). We observed the increased expression of representative Notch-related lncRNAs and genes related to DNA repair and spermatogenesis (Figure 4B,C). On the other hand, the analysis of benign or borderline ovarian tumors (n = 16) and clear cell ovarian cancers (n = 5) identified low NOTCH1/3 expressions and also presented low expressions of Notch-related lncRNAs and their related gene sets, compared to those in HGSC (n = 79) (Figure 4D–G). In summary, all our validation experiments and the FuncPred analysis clearly supported our analysis using TANRIC database, showing that Notch-related lncRNAs and genes were frequently upregulated in HGSC.

3.4. Simultaneous Expression of Common lncRNAs and Genes Related to DNA Repair and Spermatogenesis

Because the transcription of lncRNA is simultaneously accompanied by the expression of subsets of genes, we performed correlation analyses between the 43 common lncRNAs and genes related to DNA repair (n = 150) and spermatogenesis (n = 135) using TANRIC and our total transcriptome data from siRNA experiments. Whereas a large number of genes and lncRNAs showed significant positive correlations, the analysis of our total transcriptome data from the siRNA experiments presented higher correlation coefficients (Figure 5A,B), suggesting a direct relationship in their expression mechanisms. Based on these findings, we sought to reveal the TFs commonly utilized by the 43 lncRNAs and genes related to DNA repair and spermatogenesis, using ENCODE (Figure 5C). Although only 21 out of 43 lncRNAs could be evaluated, we were able to discover common TFs for them, including Early Growth Response 1 (EGR1), CCCTC-binding factor (CTCF), GA-binding protein α (GABPA), and E2F Transcription Factor 4 (E2F4) (Figure 5D). The genes related to DNA repair and spermatogenesis also utilize these TFs (Figure 5E). GSEA also showed that transcriptional targets of EGR1 and E2F4 were highly enriched in the high NOTCH1/3 group in TCGA OVCA and in the siNOTCH1/3 group in our transcriptome data. However, the target gene sets of CTCF and GABPα could not be analyzed as these TFs were not provided by MSigDB (Figure 5F). Taken together, our correlation analyses and the identification of common TFs using ENCODE suggested that Notch-related lncRNAs may share master TFs and work cooperatively with their related gene sets.

4. Discussion

Ovarian cancer, especially high-grade serous ovarian carcinoma, is one of the most lethal human malignancies [2,10,24]. Patients with ovarian cancer usually present with advanced disease stages and show frequent recurrences, that are refractory against platinum-based conventional chemotherapy [8]. To overcome these therapeutic limitations, novel targets, such as Notch signaling, have been rigorously investigated [24,36]. Several Notch-targeting drugs, including GSIs, have been tested in phase I and phase II trials [23]. Transcription factors of the Hes and Hes-related (Hey) families are the most well-known targets in the Notch pathway [37]. Cell cycle regulators, cyclin D1 and p21, and NF-κB family members, c-Myc and Deltex, are also included in Notch target genes [38,39,40,41]. As shown by our results from GSEA, the Notch pathway has a cooperative relationship with various signaling pathways, including WNT/β-catenin, TGFβ, and Hedgehog signaling [42,43,44]. This complex crosstalk with other signaling pathways have been reported to highly depend on the cellular context and environment [45]. However, the exact mechanisms of this crosstalk remain to be elucidated.
In addition to its interaction with other signaling pathways, Notch has been reported to induce lncRNA expressions to promote tumor progression in colorectal cancer [26]. The expression of lncRNAs in ovarian cancer has been investigated, suggesting a possible role in tumor progression, estrogen responses, and resistance to conventional chemotherapy [46]. Among the previously reported lncRNAs, some lncRNAs are Notch-related, for example, HOTAIR and MALAT1 [47,48,49,50,51,52,53,54]. However, to the best of our knowledge, this study is the first to investigate Notch-related lncRNAs in ovarian cancer using a genome-wide approach. Using TANRIC database, we identified 1757 differentially expressed lncRNAs based on NOTCH1/3 expressions; 1670 were upregulated, and 87 downregulated. Because data acquired from tumor tissue can be distorted by the complexity of multiple signaling pathways, we tried to confirm our results from TANRIC with siRNA experiments, which allowed the observation of relatively short time changes. After the transfection with NOTCH1/3-target siRNAs, we obtained 550 differentially expressed lncRNAs and selected 216 down-regulated lncRNAs. Finally, a comparative analysis of these downregulated lncRNAs with our data from TANRIC resulted in the identification of 43 common lncRNAs.
An interesting aspect of our functional analysis was that the genes related to specific functions were consistently associated with Notch-related lncRNAs. GSEA using TCGA OVCA revealed the signaling complexity related to Notch signaling. However, our functional analysis using FuncPred suggested more diverse gene sets related to Notch-related lncRNAs. Furthermore, as the number of representative lncRNAs decreased as we proceeded from TANRIC to siRNA experiments, the complexity of these signal changes was reproducible. Even a reduction in the number of representative lncRNAs from 1670 to 43 did not significantly change the results of these functional analyses.
The diverse roles of lncRNA in cell biology have been intensively investigated [55]. However, the regulatory mechanism of lncRNA expression itself has not been well understood. Independent transcriptional regulatory mechanisms might exist, or their expression might be related to DNA methylation [56]. In addition, the regulation of lncRNA expression might occur simultaneously with the process of regulating the coding genes [18,57]. Our correlation analysis supported the simultaneous expression of lncRNAs with coding genes. As we performed the correlation analysis using total transcriptome data from siRNA experiments, more than 50% of genes related to DNA repair and spermatogenesis presented high correlation coefficients, indicating the simultaneous expression of lncRNAs and coding genes. In fact, the FuncPred algorithm for predicting the lncRNA functions is based on the proximity of the lncRNAs and genes in the chromosome. Our analysis using ENCODE also supported this hypothesis. The 21 lncRNAs out of 43 investigated in ENCODE showed putative common TFs, such as EGR1, CTCF, GABPα, and E2F4, which were likewise found in novel gene sets related to DNA repair and spermatogenesis. In addition, our GSEA showed that the transcriptional targets of these common TFs were coordinately enriched in the high NOTCH1/3 and siNOTCH1/3 groups, indicating that the molecular targets in Notch upregulated ovarian cancer might be orchestrated by these master TFs [58,59,60,61].

5. Conclusions

In this study, we investigated the Notch-related lncRNAs using an in silico analysis and in vitro experiments, followed by validation using clinical samples in human ovarian cancers. Through this, we identified novel lncRNAs related to Notch signaling and analyzed their close relationship with gene sets involved in functional aspects, such as DNA repair and spermatogenesis, as well as regulatory mechanisms of expression. The results from this study enhance our understanding of the role of lncRNAs in the regulation of Notch signaling and their target genes in ovarian cancer. Although our data did not reveal the regulatory mechanism of master TFs, the precise regulatory mechanisms of Notch-related lncRNAs should be elucidated in the future with further research.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/cancers14061557/s1, Figure S1: Scheme of our analytic strategy. Figure S2: Grouping of ovarian cancer according to NOTCH1/3 mRNA expression status derived from TCGA HGS-OVCA (n = 303). Figure S3: Validation of the silencing effect of siNOTCH1 and siNOTCH3. Figure S4: NOTCH1 and NOTCH3 knock-down by siRNA and NOTCH inhibitor treatment (DAPT). Figure S5: Original uncropped western blot images related to Figures. Table S1: Primer sequences used for the RT-PCR.

Author Contributions

S.G.J. and J.L. designed the project and supervised the research; S.J., S.P., Y.S.J. and M.J.C. performed the experiments and the statistical analyses; M.J.C., G.L., S.G.L., M.C.C., H.P. and W.D.J. collected and provided data; S.J., S.P., Y.S.J., S.G.J. and J.L. wrote the manuscript; all authors contributed to discussions on the research and reviewed the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

S.P. was supported by the National Research Foundation of Korea (NRF) grant funded by the Korean government (MSIT) (NRF-2020R1I1A1A01069524). J.L. was supported by the National Research Foundation of Korea (NRF) grant funded by the Korean government (MSIT) (NRF-2020R1A2C1006047). S.G.J. was supported by the National Research Foundation of Korea (NRF) grant funded by the Korean government (MSIT) (NRF-2018R1D1A1B07045227).

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki, and approved by the Institutional Review Board of Bundang CHA Medical Center (IRB No. 2018-10-008, 10 October 2018).

Informed Consent Statement

Written informed consent has been obtained from the patients to publish this paper.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

The authors thank Ji Young Kim (Severance Hospital, Seoul, South Korea), Hwa Young Lee (Severance Hospital), Hee Chang Yu (Severance Hospital), and Hoyoung Kim (Severance Hospital) for technical support.

Conflicts of Interest

The authors declare that there is no conflict of interest that could be perceived as prejudicing the impartiality of this study.

References

  1. Shih, I.M.; Kurman, R.J. Ovarian tumorigenesis: A proposed model based on morphological and molecular genetic analysis. Am. J. Pathol. 2004, 164, 1511–1518. [Google Scholar] [CrossRef]
  2. Hennessy, B.T.; Coleman, R.L.; Markman, M. Ovarian cancer. Lancet 2009, 374, 1371–1382. [Google Scholar] [CrossRef]
  3. Siegel, R.; Ma, J.; Zou, Z.; Jemal, A. Cancer statistics, 2014. CA Cancer J. Clin. 2014, 64, 9–29. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Jones, M.R.; Kamara, D.; Karlan, B.Y.; Pharoah, P.D.P.; Gayther, S.A. Genetic epidemiology of ovarian cancer and prospects for polygenic risk prediction. Gynecol. Oncol. 2017, 147, 705–713. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Zhang, Y.; Luo, G.; Li, M.; Guo, P.; Xiao, Y.; Ji, H.; Hao, Y. Global patterns and trends in ovarian cancer incidence: Age, period and birth cohort analysis. BMC Cancer 2019, 19, 984. [Google Scholar] [CrossRef] [PubMed]
  6. Olivier, R.I.; van Beurden, M.; Lubsen, M.A.; Rookus, M.A.; Mooij, T.M.; van de Vijver, M.J.; van’t Veer, L.J. Clinical outcome of prophylactic oophorectomy in BRCA1/BRCA2 mutation carriers and events during follow-up. Br. J. Cancer 2004, 90, 1492–1497. [Google Scholar] [CrossRef] [PubMed]
  7. Armstrong, D.K.; Bundy, B.; Wenzel, L.; Huang, H.Q.; Baergen, R.; Lele, S.; Copeland, L.J.; Walker, J.L.; Burger, R.A.; Gynecologic Oncology, G. Intraperitoneal cisplatin and paclitaxel in ovarian cancer. N. Engl. J. Med. 2006, 354, 34–43. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Manzanedo, I.; Pereira, F.; Serrano, A.; Perez-Viejo, E.; Martinez-Torres, B.; Carrion, L.; Calzas, J. The use of cisplatin plus doxorubicin or paclitaxel in hyperthermic intraperitoneal chemotherapy (HIPEC) for stage IIIC or IV epithelial ovarian cancer: A comparative study. Clin. Transl. Oncol. 2019, 21, 1357–1363. [Google Scholar] [CrossRef]
  9. Schmeler, K.M.; Sun, C.C.; Bodurka, D.C.; Deavers, M.T.; Malpica, A.; Coleman, R.L.; Ramirez, P.T.; Gershenson, D.M. Neoadjuvant chemotherapy for low-grade serous carcinoma of the ovary or peritoneum. Gynecol. Oncol. 2008, 108, 510–514. [Google Scholar] [CrossRef]
  10. Kurman, R.J.; Shih, I.M. Molecular pathogenesis and extraovarian origin of epithelial ovarian cancer--shifting the paradigm. Hum. Pathol. 2011, 42, 918–931. [Google Scholar] [CrossRef] [Green Version]
  11. Cobb, L.P.; Sun, C.C.; Iyer, R.; Nick, A.M.; Fleming, N.D.; Westin, S.N.; Sood, A.K.; Wong, K.K.; Silva, E.G.; Gershenson, D.M. The role of neoadjuvant chemotherapy in the management of low-grade serous carcinoma of the ovary and peritoneum: Further evidence of relative chemoresistance. Gynecol. Oncol. 2020, 158, 653–658. [Google Scholar] [CrossRef]
  12. Giusti, I.; Bianchi, S.; Nottola, S.A.; Macchiarelli, G.; Dolo, V. Clinical electron microscopy in the study of human ovarian tissues. EuroMediterr. Biomed. J. 2019, 14, 145–151. [Google Scholar]
  13. Pinato, D.J.; Graham, J.; Gabra, H.; Sharma, R. Evolving concepts in the management of drug resistant ovarian cancer: Dose dense chemotherapy and the reversal of clinical platinum resistance. Cancer Treat. Rev. 2013, 39, 153–160. [Google Scholar] [CrossRef] [PubMed]
  14. Fatica, A.; Bozzoni, I. Long non-coding RNAs: New players in cell differentiation and development. Nat. Rev. Genet. 2014, 15, 7–21. [Google Scholar] [CrossRef]
  15. Esteller, M. Non-coding RNAs in human disease. Nat. Rev. Genet. 2011, 12, 861–874. [Google Scholar] [CrossRef]
  16. Chen, X.; Sun, Y.-Z.; Guan, N.-N.; Qu, J.; Huang, Z.-A.; Zhu, Z.-X.; Li, J.-Q. Computational models for lncRNA function prediction and functional similarity calculation. Brief. Funct. Genom. 2019, 18, 58–82. [Google Scholar] [CrossRef]
  17. Perron, U.; Provero, P.; Molineris, I. In silico prediction of lncRNA function using tissue specific and evolutionary conserved expression. BMC Bioinform. 2017, 18, 144. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Akrami, R.; Jacobsen, A.; Hoell, J.; Schultz, N.; Sander, C.; Larsson, E. Comprehensive analysis of long non-coding RNAs in ovarian cancer reveals global patterns and targeted DNA amplification. PLoS ONE 2013, 8, e80306. [Google Scholar] [CrossRef] [Green Version]
  19. Dai, Z.Y.; Jin, S.M.; Luo, H.Q.; Leng, H.L.; Fang, J.D. LncRNA HOTAIR regulates anoikis-resistance capacity and spheroid formation of ovarian cancer cells by recruiting EZH2 and influencing H3K27 methylation. Neoplasma 2021, 68, 509–518. [Google Scholar] [CrossRef] [PubMed]
  20. Bray, S.J. Notch signalling: A simple pathway becomes complex. Nat. Rev. Mol. Cell Biol. 2006, 7, 678–689. [Google Scholar] [CrossRef]
  21. Shawber, C.; Boulter, J.; Lindsell, C.E.; Weinmaster, G. Jagged2: A serrate-like gene expressed during rat embryogenesis. Dev. Biol. 1996, 180, 370–376. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Ikeuchi, T.; Sisodia, S.S. The Notch ligands, Delta1 and Jagged2, are substrates for presenilin-dependent “gamma-secretase” cleavage. J. Biol. Chem. 2003, 278, 7751–7754. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Espinoza, I.; Miele, L. Notch inhibitors for cancer treatment. Pharmacol. Ther. 2013, 139, 95–110. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Groeneweg, J.W.; Foster, R.; Growdon, W.B.; Verheijen, R.H.; Rueda, B.R. Notch signaling in serous ovarian cancer. J. Ovarian Res. 2014, 7, 95. [Google Scholar] [CrossRef]
  25. Locatelli, M.; Curigliano, G. Notch inhibitors and their role in the treatment of triple negative breast cancer: Promises and failures. Curr. Opin. Oncol. 2017, 29, 411–427. [Google Scholar] [CrossRef]
  26. Zhang, Z.; Li, G.; Qiu, H.; Yang, J.; Bu, X.; Zhu, S.; Zheng, J.; Dang, C.; Wang, W.; Chu, D. The Novel Notch-induced Long Noncoding RNA LUNAR1 Determines the Proliferation and Prognosis of Colorectal Cancer. Sci. Rep. 2019, 9, 19915. [Google Scholar] [CrossRef] [PubMed]
  27. The Cancer Genome Atlas Research Network. Integrated genomic analyses of ovarian carcinoma. Nature 2011, 474, 609–615. [Google Scholar] [CrossRef]
  28. Zhang, Q.; Burdette, J.E.; Wang, J.P. Integrative network analysis of TCGA data for ovarian cancer. BMC Syst. Biol. 2014, 8, 1338. [Google Scholar] [CrossRef] [Green Version]
  29. Li, J.; Han, L.; Roebuck, P.; Diao, L.; Liu, L.; Yuan, Y.; Weinstein, J.N.; Liang, H. TANRIC: An Interactive Open Platform to Explore the Function of lncRNAs in Cancer. Cancer Res. 2015, 75, 3728–3737. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Subramanian, A.; Tamayo, P.; Mootha, V.K.; Mukherjee, S.; Ebert, B.L.; Gillette, M.A.; Paulovich, A.; Pomeroy, S.L.; Golub, T.R.; Lander, E.S.; et al. Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles. Proc. Natl. Acad. Sci. USA 2005, 102, 15545–15550. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Guigo, R.; Flicek, P.; Abril, J.F.; Reymond, A.; Lagarde, J.; Denoeud, F.; Antonarakis, S.; Ashburner, M.; Bajic, V.B.; Birney, E.; et al. EGASP: The human ENCODE Genome Annotation Assessment Project. Genome Biol. 2006, 7, S2. [Google Scholar] [CrossRef] [Green Version]
  32. Xie, X.; Lu, J.; Kulbokas, E.J.; Golub, T.R.; Mootha, V.; Lindblad-Toh, K.; Lander, E.S.; Kellis, M. Systematic discovery of regulatory motifs in human promoters and 3’ UTRs by comparison of several mammals. Nature 2005, 434, 338–345. [Google Scholar] [CrossRef] [PubMed]
  33. Chan, Y.M.; Jan, Y.N. Roles for proteolysis and trafficking in notch maturation and signal transduction. Cell 1998, 94, 423–426. [Google Scholar] [CrossRef] [Green Version]
  34. Schroeter, E.H.; Kisslinger, J.A.; Kopan, R. Notch-1 signalling requires ligand-induced proteolytic release of intracellular domain. Nature 1998, 393, 382–386. [Google Scholar] [CrossRef] [PubMed]
  35. Rand, M.D.; Grimm, L.M.; Artavanis-Tsakonas, S.; Patriub, V.; Blacklow, S.C.; Sklar, J.; Aster, J.C. Calcium depletion dissociates and activates heterodimeric notch receptors. Mol. Cell Biol. 2000, 20, 1825–1835. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Jung, S.G.; Kwon, Y.D.; Song, J.A.; Back, M.J.; Lee, S.Y.; Lee, C.; Hwang, Y.Y.; An, H.J. Prognostic significance of Notch 3 gene expression in ovarian serous carcinoma. Cancer Sci. 2010, 101, 1977–1983. [Google Scholar] [CrossRef] [PubMed]
  37. Iso, T.; Kedes, L.; Hamamori, Y. HES and HERP families: Multiple effectors of the Notch signaling pathway. J. Cell Physiol. 2003, 194, 237–255. [Google Scholar] [CrossRef] [PubMed]
  38. Oswald, F.; Liptay, S.; Adler, G.; Schmid, R.M. NF-kappaB2 is a putative target gene of activated Notch-1 via RBP-Jkappa. Mol. Cell Biol. 1998, 18, 2077–2088. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Ronchini, C.; Capobianco, A.J. Induction of cyclin D1 transcription and CDK2 activity by Notch(ic): Implication for cell cycle disruption in transformation by Notch(ic). Mol. Cell Biol. 2001, 21, 5925–5934. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Choi, J.W.; Pampeno, C.; Vukmanovic, S.; Meruelo, D. Characterization of the transcriptional expression of Notch-1 signaling pathway members, Deltex and HES-1, in developing mouse thymocytes. Dev. Comp. Immunol. 2002, 26, 575–588. [Google Scholar] [CrossRef]
  41. Weng, A.P.; Millholland, J.M.; Yashiro-Ohtani, Y.; Arcangeli, M.L.; Lau, A.; Wai, C.; Del Bianco, C.; Rodriguez, C.G.; Sai, H.; Tobias, J.; et al. c-Myc is an important direct target of Notch1 in T-cell acute lymphoblastic leukemia/lymphoma. Genes Dev. 2006, 20, 2096–2109. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Hurlbut, G.D.; Kankel, M.W.; Lake, R.J.; Artavanis-Tsakonas, S. Crossing paths with Notch in the hyper-network. Curr. Opin. Cell Biol. 2007, 19, 166–175. [Google Scholar] [CrossRef] [PubMed]
  43. Brechbiel, J.; Miller-Moslin, K.; Adjei, A.A. Crosstalk between hedgehog and other signaling pathways as a basis for combination therapies in cancer. Cancer Treat. Rev. 2014, 40, 750–759. [Google Scholar] [CrossRef]
  44. Duncan, A.W.; Rattis, F.M.; DiMascio, L.N.; Congdon, K.L.; Pazianos, G.; Zhao, C.; Yoon, K.; Cook, J.M.; Willert, K.; Gaiano, N.; et al. Integration of Notch and Wnt signaling in hematopoietic stem cell maintenance. Nat. Immunol. 2005, 6, 314–322. [Google Scholar] [CrossRef] [PubMed]
  45. D’Souza, B.; Miyamoto, A.; Weinmaster, G. The many facets of Notch ligands. Oncogene 2008, 27, 5148–5167. [Google Scholar] [CrossRef] [Green Version]
  46. Nikpayam, E.; Tasharrofi, B.; Sarrafzadeh, S.; Ghafouri-Fard, S. The Role of Long Non-Coding RNAs in Ovarian Cancer. Iran. Biomed. J. 2017, 21, 3–15. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Qiu, J.J.; Lin, Y.Y.; Ding, J.X.; Feng, W.W.; Jin, H.Y.; Hua, K.Q. Long non-coding RNA ANRIL predicts poor prognosis and promotes invasion/metastasis in serous ovarian cancer. Int. J. Oncol. 2015, 46, 2497–2505. [Google Scholar] [CrossRef] [Green Version]
  48. Hu, X.; Feng, Y.; Zhang, D.; Zhao, S.D.; Hu, Z.; Greshock, J.; Zhang, Y.; Yang, L.; Zhong, X.; Wang, L.P.; et al. A functional genomic approach identifies FAL1 as an oncogenic long noncoding RNA that associates with BMI1 and represses p21 expression in cancer. Cancer Cell 2014, 26, 344–357. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Cheng, Z.; Guo, J.; Chen, L.; Luo, N.; Yang, W.; Qu, X. A long noncoding RNA AB073614 promotes tumorigenesis and predicts poor prognosis in ovarian cancer. Oncotarget 2015, 6, 25381–25389. [Google Scholar] [CrossRef] [Green Version]
  50. Yin, G.; Chen, R.; Alvero, A.B.; Fu, H.H.; Holmberg, J.; Glackin, C.; Rutherford, T.; Mor, G. TWISTing stemness, inflammation and proliferation of epithelial ovarian cancer cells through MIR199A2/214. Oncogene 2010, 29, 3545–3553. [Google Scholar] [CrossRef] [Green Version]
  51. Tanos, V.; Prus, D.; Ayesh, S.; Weinstein, D.; Tykocinski, M.L.; De-Groot, N.; Hochberg, A.; Ariel, I. Expression of the imprinted H19 oncofetal RNA in epithelial ovarian cancer. Eur. J. Obstet. Gynecol. Reprod. Biol. 1999, 85, 7–11. [Google Scholar] [CrossRef]
  52. Gao, Y.; Meng, H.; Liu, S.; Hu, J.; Zhang, Y.; Jiao, T.; Liu, Y.; Ou, J.; Wang, D.; Yao, L.; et al. LncRNA-HOST2 regulates cell biological behaviors in epithelial ovarian cancer through a mechanism involving microRNA let-7b. Hum. Mol. Genet. 2015, 24, 841–852. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Qiu, J.J.; Wang, Y.; Ding, J.X.; Jin, H.Y.; Yang, G.; Hua, K.Q. The long non-coding RNA HOTAIR promotes the proliferation of serous ovarian cancer cells through the regulation of cell cycle arrest and apoptosis. Exp. Cell Res. 2015, 333, 238–248. [Google Scholar] [CrossRef]
  54. Silva, J.M.; Boczek, N.J.; Berres, M.W.; Ma, X.; Smith, D.I. LSINCT5 is over expressed in breast and ovarian cancer and affects cellular proliferation. RNA Biol. 2011, 8, 496–505. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Klinge, C.M. Non-coding RNAs: Long non-coding RNAs and microRNAs in endocrine-related cancers. Endocr. Relat. Cancer 2018, 25, R259–R282. [Google Scholar] [CrossRef] [PubMed]
  56. Sheng, X.; Li, J.; Yang, L.; Chen, Z.; Zhao, Q.; Tan, L.; Zhou, Y.; Li, J. Promoter hypermethylation influences the suppressive role of maternally expressed 3, a long non-coding RNA, in the development of epithelial ovarian cancer. Oncol. Rep. 2014, 32, 277–285. [Google Scholar] [CrossRef] [Green Version]
  57. Jeong, S.; Lee, S.G.; Kim, H.; Lee, G.; Park, S.; Kim, I.K.; Lee, J.; Jo, Y.S. Simultaneous Expression of Long Non-Coding RNA FAL1 and Extracellular Matrix Protein 1 Defines Tumour Behaviour in Young Patients with Papillary Thyroid Cancer. Cancers 2021, 13, 3223. [Google Scholar] [CrossRef] [PubMed]
  58. Helleday, T.; Petermann, E.; Lundin, C.; Hodgson, B.; Sharma, R.A. DNA repair pathways as targets for cancer therapy. Nat. Rev. Cancer 2008, 8, 193–204. [Google Scholar] [CrossRef] [PubMed]
  59. Old, L.J. Cancer/testis (CT) antigens—A new link between gametogenesis and cancer. Cancer Immun. 2001, 1, 1. [Google Scholar] [PubMed]
  60. Simpson, A.J.; Caballero, O.L.; Jungbluth, A.; Chen, Y.T.; Old, L.J. Cancer/testis antigens, gametogenesis and cancer. Nat. Rev. Cancer 2005, 5, 615–625. [Google Scholar] [CrossRef] [PubMed]
  61. Whyte, W.A.; Orlando, D.A.; Hnisz, D.; Abraham, B.J.; Lin, C.Y.; Kagey, M.H.; Rahl, P.B.; Lee, T.I.; Young, R.A. Master transcription factors and mediator establish super-enhancers at key cell identity genes. Cell 2013, 153, 307–319. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Cancers 14 01557 g001 550
Figure 1. Discovery of Notch-related lncRNAs in HGSC. (A) Results of GSEA using Hallmark gene sets in the TCGA HGS-OVCA database (N = 303) to validate our grouping according to NOTCH1/3 expression. “High” indicates a group with high expressions of both NOTCH1 and NOTCH3 (n = 97). “Low” indicates a group with low expressions of both NOTCH1 and NOTHC3 (n = 96). The classification into high and low was based on median TPM values. High NOTCH1 with low NOTCH3 (n = 55) and high NOTCH3 with low NOTCH1 (n = 55) were excluded from this analysis. Red coloured dots indicated statistically significant p-value (<0.05) and FDR q-value (<0.025). (B) Identification of differentially expressed lncRNAs using TANRIC database (n = 293). Ten of the TCGA HGS-OVCA cases (n = 303) did not show lncRNA expression. Heatmap showing representative lncRNAs (n = 1757) according to NOTCH1/3 expression status. Based on the expression status of NOTCH1/3, TCGA HGS-OVCA dataset was divided into four groups: I (n = 92), low expression of both NOTCH1 and 3; II (n = 54), high expression of NOTCH3 and low expression of NOTCH1; III (n = 54), high expression of NOTCH1 and low expression of NOTCH3; and IV (n = 93), high expression of both NOTCH1 and 3. (C) Potential function of Notch-related lncRNAs using FuncPred. Red coloured dots indicate unique gene sets highly ranked in this analysis but not in GSEA depicted in (A). (D) Venn diagram presenting the logical relation between FuncPred data with lncRNAs and GSEA data with mRNAs. Abbreviations: NOTCH1/3, Notch receptor 1 and 3; TCGA, The Cancer Genome Atlas; HGS-OVCA, high-grade serous ovarian cancer; lncRNA, long non-coding RNA; GSEA, gene set enrichment analysis; FuncPred, in silico prediction of gene function for protein-coding and lncRNA human genes.
Figure 1. Discovery of Notch-related lncRNAs in HGSC. (A) Results of GSEA using Hallmark gene sets in the TCGA HGS-OVCA database (N = 303) to validate our grouping according to NOTCH1/3 expression. “High” indicates a group with high expressions of both NOTCH1 and NOTCH3 (n = 97). “Low” indicates a group with low expressions of both NOTCH1 and NOTHC3 (n = 96). The classification into high and low was based on median TPM values. High NOTCH1 with low NOTCH3 (n = 55) and high NOTCH3 with low NOTCH1 (n = 55) were excluded from this analysis. Red coloured dots indicated statistically significant p-value (<0.05) and FDR q-value (<0.025). (B) Identification of differentially expressed lncRNAs using TANRIC database (n = 293). Ten of the TCGA HGS-OVCA cases (n = 303) did not show lncRNA expression. Heatmap showing representative lncRNAs (n = 1757) according to NOTCH1/3 expression status. Based on the expression status of NOTCH1/3, TCGA HGS-OVCA dataset was divided into four groups: I (n = 92), low expression of both NOTCH1 and 3; II (n = 54), high expression of NOTCH3 and low expression of NOTCH1; III (n = 54), high expression of NOTCH1 and low expression of NOTCH3; and IV (n = 93), high expression of both NOTCH1 and 3. (C) Potential function of Notch-related lncRNAs using FuncPred. Red coloured dots indicate unique gene sets highly ranked in this analysis but not in GSEA depicted in (A). (D) Venn diagram presenting the logical relation between FuncPred data with lncRNAs and GSEA data with mRNAs. Abbreviations: NOTCH1/3, Notch receptor 1 and 3; TCGA, The Cancer Genome Atlas; HGS-OVCA, high-grade serous ovarian cancer; lncRNA, long non-coding RNA; GSEA, gene set enrichment analysis; FuncPred, in silico prediction of gene function for protein-coding and lncRNA human genes.
Cancers 14 01557 g001
Cancers 14 01557 g002 550
Figure 2. Discovery of Notch-related lncRNAs using transcriptome data from OVCAR3 cells treated with siNOTCH1/3. (A) Heatmap showing differentially expressed lncRNAs (n = 550) according to siNOTCH1/3 treatment. (B) Predicted function of Notch-related lncRNAs based on our siRNA experiments using FuncPred. Red coloured dots indicated unique gene sets highly ranked in this analysis, but not in GSEA, as depicted in Figure 1A. (C) Venn diagram presenting the logical relation between upregulated lncRNAs in high NOTCH1/3 group in TCGA OVCA and downregulated lncRNAs in our siRNA experiments. (D) List of common lncRNAs, as depicted in (C). (E) Potential function of the common lncRNAs listed in (D). Red coloured dots indicate unique gene sets highly ranked in this analysis, but not in GSEA, as depicted in Figure 1A. Abbreviations: TANRIC, The Atlas of Noncoding RNAs in Cancer.
Figure 2. Discovery of Notch-related lncRNAs using transcriptome data from OVCAR3 cells treated with siNOTCH1/3. (A) Heatmap showing differentially expressed lncRNAs (n = 550) according to siNOTCH1/3 treatment. (B) Predicted function of Notch-related lncRNAs based on our siRNA experiments using FuncPred. Red coloured dots indicated unique gene sets highly ranked in this analysis, but not in GSEA, as depicted in Figure 1A. (C) Venn diagram presenting the logical relation between upregulated lncRNAs in high NOTCH1/3 group in TCGA OVCA and downregulated lncRNAs in our siRNA experiments. (D) List of common lncRNAs, as depicted in (C). (E) Potential function of the common lncRNAs listed in (D). Red coloured dots indicate unique gene sets highly ranked in this analysis, but not in GSEA, as depicted in Figure 1A. Abbreviations: TANRIC, The Atlas of Noncoding RNAs in Cancer.
Cancers 14 01557 g002
Cancers 14 01557 g003 550
Figure 3. Validation of in silico analysis and transcriptome data in HGSC cell line. Results of qRT-PCR for mRNAs of representative Notch-related lncRNAs (A), DNA repair-related targets (B), and Spermatogenesis-related targets (C) in OVCAR3 cells transfected with control siRNA (siControl) or NOTCH 1/3 siRNAs as indicated for 48 h. Comparisons of the relative expression of representative Notch-related lncRNAs (D), genes related to DNA repair (E), and spermatogenesis (F) in OVCAR3 cells following treatment with DAPT for 48 h. Data is shown as mean ± SD and p-values were calculated by two-tailed Mann–Whitney U-test. All experiments were repeated three times, and each experiment was performed in triplicate. * p < 0.05, ** p < 0.01, *** p < 0.001. Abbreviations: CTL, Control; TAF1C, TATA-Box Binding Protein Associated Factor; ADCY6, Adenylate Cyclase Type 6; DGCR8, DiGeorge Syndrome Critical Region Gene 8; RAD52, DNA Repair Protein RAD52 Homolog; ACRV1, Acrosomal Vesicle Protein 1; NPY5R, Neuropeptide Y Receptor Y5; DMC1, DNA Meiotic Recombinase 1; TCP11, T-Complex 11, DAPT: N-[N-(3,5-difluorophenacetyl)-L-alanyl]-S-phenylglycine t-butyl ester.
Figure 3. Validation of in silico analysis and transcriptome data in HGSC cell line. Results of qRT-PCR for mRNAs of representative Notch-related lncRNAs (A), DNA repair-related targets (B), and Spermatogenesis-related targets (C) in OVCAR3 cells transfected with control siRNA (siControl) or NOTCH 1/3 siRNAs as indicated for 48 h. Comparisons of the relative expression of representative Notch-related lncRNAs (D), genes related to DNA repair (E), and spermatogenesis (F) in OVCAR3 cells following treatment with DAPT for 48 h. Data is shown as mean ± SD and p-values were calculated by two-tailed Mann–Whitney U-test. All experiments were repeated three times, and each experiment was performed in triplicate. * p < 0.05, ** p < 0.01, *** p < 0.001. Abbreviations: CTL, Control; TAF1C, TATA-Box Binding Protein Associated Factor; ADCY6, Adenylate Cyclase Type 6; DGCR8, DiGeorge Syndrome Critical Region Gene 8; RAD52, DNA Repair Protein RAD52 Homolog; ACRV1, Acrosomal Vesicle Protein 1; NPY5R, Neuropeptide Y Receptor Y5; DMC1, DNA Meiotic Recombinase 1; TCP11, T-Complex 11, DAPT: N-[N-(3,5-difluorophenacetyl)-L-alanyl]-S-phenylglycine t-butyl ester.
Cancers 14 01557 g003
Cancers 14 01557 g004 550
Figure 4. Validation of in silico analysis and transcriptome data in human ovarian tissue samples. (A–C) Comparisons of the relative expression of NOTCH1/NOTCH3 (A), representative Notch-related lncRNAs (B), genes related to DNA repair (TAF1C and ADCY6) and spermatogenesis (C). (ACRV1) using our own HGSC tissue samples according to NOTCH1/3 expression (high vs. low, n = 21 vs. n = 19, respectively). (D–G). Relative expression values of NOTC1/3 (D), notch-related lncRNAs (E), and representative genes related to DNA repair (F), and spermatogenesis (G), in benign or borderline ovarian tumors (n = 16), clear cell ovarian carcinomas (n = 5), and high-grade serous ovarian cancers (HGSC, n = 79). Data was shown as mean ± SD and p-values were calculated by two-tailed Mann–Whitney U-test. All experiments were repeated three times, and each experiment was performed in triplicate. * p < 0.05, ** p < 0.01, *** p < 0.001. Abbreviations: TAF1C, TATA-Box Binding Protein Associated Factor; ADCY6, Adenylate Cyclase Type 6; DGCR8, DiGeorge Syndrome Critical Region Gene 8; RAD52, DNA Repair Protein RAD52 Homolog; ACRV1, Acrosomal Vesicle Protein 1; NPY5R, Neuropeptide Y Receptor Y5; DMC1, DNA Meiotic Recombinase 1; TCP11, T-Complex 11.
Figure 4. Validation of in silico analysis and transcriptome data in human ovarian tissue samples. (A–C) Comparisons of the relative expression of NOTCH1/NOTCH3 (A), representative Notch-related lncRNAs (B), genes related to DNA repair (TAF1C and ADCY6) and spermatogenesis (C). (ACRV1) using our own HGSC tissue samples according to NOTCH1/3 expression (high vs. low, n = 21 vs. n = 19, respectively). (D–G). Relative expression values of NOTC1/3 (D), notch-related lncRNAs (E), and representative genes related to DNA repair (F), and spermatogenesis (G), in benign or borderline ovarian tumors (n = 16), clear cell ovarian carcinomas (n = 5), and high-grade serous ovarian cancers (HGSC, n = 79). Data was shown as mean ± SD and p-values were calculated by two-tailed Mann–Whitney U-test. All experiments were repeated three times, and each experiment was performed in triplicate. * p < 0.05, ** p < 0.01, *** p < 0.001. Abbreviations: TAF1C, TATA-Box Binding Protein Associated Factor; ADCY6, Adenylate Cyclase Type 6; DGCR8, DiGeorge Syndrome Critical Region Gene 8; RAD52, DNA Repair Protein RAD52 Homolog; ACRV1, Acrosomal Vesicle Protein 1; NPY5R, Neuropeptide Y Receptor Y5; DMC1, DNA Meiotic Recombinase 1; TCP11, T-Complex 11.
Cancers 14 01557 g004
Cancers 14 01557 g005 550
Figure 5. Identification of common transcription factors for Notch-related lncRNAs and genes related to DNA repair and spermatogenesis. (A) Correlation analysis of common lncRNA clusters with genes related to DNA repair and spermatogenesis using data from TCGA OVCA and (B) our siRNA experiments. (C) Representative result of ENCODE predicting putative transcription factors. (D) Results from ENCODE using 21 out of 43 common lncRNAs presenting number of lncRNAs associated with putative transcription factors. (E) Results from ENCODE using genes related to DNA repair and spermatogenesis presenting number of genes associated with putative transcription factors. (F) Transcription factors that are up- or down-regulated according to NOTCH1/3 status. GSEA was performed using legacy subset of transcription target gene sets as described in the Methods section. Upper and lower panels presented the results of GSEA using TCGA OVCA and our transcriptome data, respectively. The legacy subset of transcription target gene sets originally supported 610 gene sets, but we excluded unknown transcription factor and applied gene set size filter (min = 15, max = 500), and finally used 466 gene sets. Each dot represents one transcription target gene set which showed adjusted p-value < 0.05 and FDR q-value < 0.25. The dotted line indicates FDR = 0.25. Abbreviations: EGR1, Early Growth Response 1; CCTF, CCCTC-Binding Factor; GABPA, GA Binding Protein Transcription Factor Subunit Alpha; E2F4, E2F Transcription Factor 4; FDR, False discovery rate.
Figure 5. Identification of common transcription factors for Notch-related lncRNAs and genes related to DNA repair and spermatogenesis. (A) Correlation analysis of common lncRNA clusters with genes related to DNA repair and spermatogenesis using data from TCGA OVCA and (B) our siRNA experiments. (C) Representative result of ENCODE predicting putative transcription factors. (D) Results from ENCODE using 21 out of 43 common lncRNAs presenting number of lncRNAs associated with putative transcription factors. (E) Results from ENCODE using genes related to DNA repair and spermatogenesis presenting number of genes associated with putative transcription factors. (F) Transcription factors that are up- or down-regulated according to NOTCH1/3 status. GSEA was performed using legacy subset of transcription target gene sets as described in the Methods section. Upper and lower panels presented the results of GSEA using TCGA OVCA and our transcriptome data, respectively. The legacy subset of transcription target gene sets originally supported 610 gene sets, but we excluded unknown transcription factor and applied gene set size filter (min = 15, max = 500), and finally used 466 gene sets. Each dot represents one transcription target gene set which showed adjusted p-value < 0.05 and FDR q-value < 0.25. The dotted line indicates FDR = 0.25. Abbreviations: EGR1, Early Growth Response 1; CCTF, CCCTC-Binding Factor; GABPA, GA Binding Protein Transcription Factor Subunit Alpha; E2F4, E2F Transcription Factor 4; FDR, False discovery rate.
Cancers 14 01557 g005
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Jeong, S.; Park, S.; Jo, Y.S.; Choi, M.J.; Lee, G.; Lee, S.G.; Choi, M.C.; Park, H.; Joo, W.D.; Jung, S.G.; et al. Long Non-Coding RNA-Based Functional Prediction Reveals Novel Targets in Notch-Upregulated Ovarian Cancer. Cancers 2022, 14, 1557. https://doi.org/10.3390/cancers14061557

AMA Style

Jeong S, Park S, Jo YS, Choi MJ, Lee G, Lee SG, Choi MC, Park H, Joo WD, Jung SG, et al. Long Non-Coding RNA-Based Functional Prediction Reveals Novel Targets in Notch-Upregulated Ovarian Cancer. Cancers. 2022; 14(6):1557. https://doi.org/10.3390/cancers14061557

Chicago/Turabian Style

Jeong, Seonhyang, Sunmi Park, Young Suk Jo, Moon Jung Choi, Gibbeum Lee, Seul Gi Lee, Min Chul Choi, Hyun Park, Won Duk Joo, Sang Geun Jung, and et al. 2022. "Long Non-Coding RNA-Based Functional Prediction Reveals Novel Targets in Notch-Upregulated Ovarian Cancer" Cancers 14, no. 6: 1557. https://doi.org/10.3390/cancers14061557

APA Style

Jeong, S., Park, S., Jo, Y. S., Choi, M. J., Lee, G., Lee, S. G., Choi, M. C., Park, H., Joo, W. D., Jung, S. G., & Lee, J. (2022). Long Non-Coding RNA-Based Functional Prediction Reveals Novel Targets in Notch-Upregulated Ovarian Cancer. Cancers, 14(6), 1557. https://doi.org/10.3390/cancers14061557

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop