NotI microarrays: novel epigenetic markers for early detection and prognosis of high grade serous ovarian cancer

Construction and analysis of competitive endogenous RNA networks and prognostic models associated with ovarian cancer based on the exoRBase database

OBJECTIVE

To construct a competitive endogenous RNA (ceRNA) regulatory network in blood exosomes of patients with ovarian cancer (OC) using bioinformatics and explore its pathogenesis.

METHODS

The exoRbase2.0 database was used to download blood exosome gene sequencing data from patients OC and normal controls and the expression profiles of exosomal mRNA, long non-coding RNA (lncRNA), and circular RNA (circRNA) were detected independently using R language for differential expression analysis. TargetScan and miRanda databases were combined for the prediction and differential expression of mRNA-binding microRNAs (miRNA). The miRcode and starBase databases were used to predict miRNAs that bind to differentially expressed lncRNAs and circRNAs repectively. The relevant mRNA, circRNA, lncRNA and their corresponding miRNA prediction data were imported into Cytoscape software for visualization of the ceRNA network. The R language and KEGG Orthology-based Annotation System (KOBAS) were used to execute and illustrate the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses. Hub genes were identified using The CytoHubba plugin.

RESULTS

Thirty-one differentially expressed mRNAs, 17 differentially expressed lncRNAs, and 24 differentially expressed circRNAs were screened. Cytoscape software was used to construct the ceRNA network with nine mRNA nodes, two lncRNA nodes, eight circRNA nodes, and 51 miRNA nodes. Both GO and KEGG were focused on the Spliceosome pathway, indicating that spliceosomes are closely linked with the development of OC, while heterogenous nuclear ribonucleoprotein K and RNA binding motif protein X-linked genes were the top 10 score Hub genes screened by Cytoscape software, including two lncRNAs, four mRNAs, and four circRNAs. In patients with OC, the expression of eukaryotic translation initiation factor 4 gamma 2 (EIF4G2), SERPINE 1 mRNA binding protein 1 (SERBP1), ribosomal protein L15 (RPL15) and human leukocyte antigen complex P5 (HCP5) was significantly higher whereas that of testis expressed transcript, Y-linked 15 and DEAD-box helicase 3 Y-linked genes was lower compared to normal controls Immunocorrelation scores revealed that SERBP1 was significantly and negatively correlated with endothelial cells and CD4+ T cells and positively correlated with natural killer (NK) cells and macrophages, respectively; RPL15 was significantly positively correlated with macrophages and endothelial cells and negatively correlated with CD8+ T cells and uncharacterized cells, respectively. EIF4G2 was significantly and negatively correlated with endothelial cells and CD4+ T cells, and positively correlated with uncharacterized cells, respectively. Based on the survival data and the significant correlation characteristics derived from the multifactorial Cox analysis (P < 0.05), the survival prediction curves demonstrated that the prognostic factors associated with 3-year survival in patients with OC were The prognostic factors associated with survival were Macrophage, RPL15.

CONCLUSION

This study successfully constructs a ceRNA regulatory network in blood exosomes of OV patients, which provides the specific targets for diagnosis and treatment of OC.

Tumor Suppressor Properties of Small C-Terminal Domain Phosphatases in Clear Cell Renal Cell Carcinoma

Clear cell renal cell carcinoma (ccRCC) accounts for 80-90% of kidney cancers worldwide. Small C-terminal domain phosphatases CTDSP1, CTDSP2, and CTDSPL (also known as SCP1, 2, 3) are involved in the regulation of several important pathways associated with carcinogenesis. In various cancer types, these phosphatases may demonstrate either antitumor or oncogenic activity. Tumor-suppressive activity of these phosphatases in kidney cancer has been shown previously, but in general case, the antitumor activity may be dependent on the choice of cell line. In the present work, transfection of the Caki-1 cell line (ccRCC morphologic phenotype) with expression constructs containing the coding regions of these genes resulted in inhibition of cell growth in vitro in the case of CTDSP1 (p < 0.001) and CTDSPL (p < 0.05) but not CTDSP2. The analysis of The Cancer Genome Atlas (TCGA) data showed differential expression of some of CTDSP genes and of their target, RB1. These results were confirmed by quantitative RT-PCR using an independent sample of primary ccRCC tumors (n = 52). We observed CTDSPL downregulation and found a positive correlation of expression for two gene pairs: CTDSP1 and CTDSP2 (rs = 0.76; p < 0.001) and CTDSPL and RB1 (rs = 0.38; p < 0.05). Survival analysis based on TCGA data demonstrated a strong association of lower expression of CTDSP1, CTDSP2, CTDSPL, and RB1 with poor survival of ccRCC patients (p < 0.001). In addition, according to TCGA, CTDSP1, CTDSP2, and RB1 were differently expressed in two subtypes of ccRCC-ccA and ccB, characterized by different survival rates. These results confirm that CTDSP1 and CTDSPL have tumor suppressor properties in ccRCC and reflect their association with the more aggressive ccRCC phenotype.

Analysis of the effect of NEKs on the prognosis of patients with non-small-cell lung carcinoma based on bioinformatics

NEKs are proteins that are involved in various cell processes and play important roles in the formation and development of cancer. However, few studies have examined the role of NEKs in the development of non-small-cell lung carcinoma (NSCLC). To address this problem, the Oncomine, UALCAN, and the Human Protein Atlas databases were used to analyze differential NEK expression and its clinicopathological parameters, while the Kaplan-Meier, cBioPortal, GEPIA, and DAVID databases were used to analyze survival, gene mutations, similar genes, and biological enrichments. The rate of NEK family gene mutation was high (> 50%) in patients with NSCLC, in which NEK2/4/6/8/ was overexpressed and significantly correlated with tumor stage and nodal metastasis status. In addition, the high expression of NEK2/3mRNA was significantly associated with poor prognosis in patients with NSCLC, while high expression of NEK1/4/6/7/8/9/10/11mRNA was associated with good prognosis. In summary, these results suggest that NEK2/4/6/8 may be a potential prognostic biomarker for the survival of patients with NSCLC.

Population structure and genetic diversity in yellow catfish (Pelteobagrus fulvidraco) assessed with microsatellites

Yellow catfish (Pelteobagrus fulvidraco) is an important aquaculture species which is widely distributed, especially in the Yangtze River of China. To facilitate its conservation and stock improvement, 273 yellow catfish samples from the Yangtze River (seven populations) and Baiyangdian (BYD) Lake were genotyped using eight microsatellites in combination with capillary electrophoresis. A total of 250 alleles were detected at eight loci in eight populations showing high allelic (N_a= 31.25 ± 7.38) and genetic diversity (H_e = 0.888-0.944). Both F_ST and clustering analyses revealed the presence of subtle population differences between the species of Yangtze River and the BYD lake. Mantel tests suggest that genetic distance is significantly correlated with geographical distance (R = 0.9294 and P < 0.05). The results of genetic diversity and population structure will help in conservation and improvement of yellow catfish.

Regulatory network involving miRNAs and genes in serous ovarian carcinoma

Serous ovarian carcinoma (SOC) is one of the most life-threatening types of gynecological malignancy, but the pathogenesis of SOC remains unknown. Previous studies have indicated that differentially expressed genes and microRNAs (miRNAs) serve important functions in SOC. However, genes and miRNAs are identified in a disperse form, and limited information is known about the regulatory association between miRNAs and genes in SOC. In the present study, three regulatory networks were hierarchically constructed, including a differentially-expressed network, a related network and a global network to reveal associations between each factor. In each network, there were three types of factors, which were genes, miRNAs and transcription factors that interact with each other. Focus was placed on the differentially-expressed network, in which all genes and miRNAs were differentially expressed and therefore may have affected the development of SOC. Following the comparison and analysis between the three networks, a number of signaling pathways which demonstrated differentially expressed elements were highlighted. Subsequently, the upstream and downstream elements of differentially expressed miRNAs and genes were listed, and a number of key elements (differentially expressed miRNAs, genes and TFs predicted using the P-match method) were analyzed. The differentially expressed network partially illuminated the pathogenesis of SOC. It was hypothesized that if there was no differential expression of miRNAs and genes, SOC may be prevented and treatment may be identified. The present study provided a theoretical foundation for gene therapy for SOC.

The Role of Epigenomics in the Study of Cancer Biomarkers and in the Development of Diagnostic Tools

Epigenetics plays a key role in cancer development. Genetics alone cannot explain sporadic cancer and cancer development in individuals with no family history or a weak family history of cancer. Epigenetics provides a mechanism to explain the development of cancer in such situations. Alterations in epigenetic profiling may provide important insights into the etiology and natural history of cancer. Because several epigenetic changes occur before histopathological changes, they can serve as biomarkers for cancer diagnosis and risk assessment. Many cancers may remain asymptomatic until relatively late stages; in managing the disease, efforts should be focused on early detection, accurate prediction of disease progression, and frequent monitoring. This chapter describes epigenetic biomarkers as they are expressed during cancer development and their potential use in cancer diagnosis and prognosis. Based on epigenomic information, biomarkers have been identified that may serve as diagnostic tools; some such biomarkers also may be useful in identifying individuals who will respond to therapy and survive longer. The importance of analytical and clinical validation of biomarkers is discussed, along with challenges and opportunities in this field.

A double-negative feedback loop between EZH2 and miR-26a regulates tumor cell growth in hepatocellular carcinoma

Accumulating evidence demonstrates the important roles of microRNAs (miRNAs) in tumor development and progression. miR-26a has been reported to be downregulated in several types of cancers including hepatocellular carcinoma, but the underlying mechanism of how miR-26a is repressed remains largely unknown. In the present study, we performed western blot analysis, qRT-PCR, luciferase reporter assay and chromatin immunoprecipitation assay to investigate the relationship between miR-26a and the enhancer of zest homologue 2 (EZH2). CCK-8 assay and colony formation assay were carried out to explore the effect of miR-26a on HCC cells proliferation. We demonstrated that miR-26a was epigenetically repressed by EZH2-mediated H3K27 trimethylation within the miR-26a promoter. Moreover, we confirmed that EZH2 was also a direct target of miR-26a in HCC cells, thus, creating a double-negative feedback loop. Furthermore, miR-26a restoration increased the expressions of its host genes (CTDSPL and CTDSP2). Overexpression of EZH2 abrogated miR-26a induction of CTDSPL and CTDSP2. Restoring the balance of the double-negative feedback loop by miR-26a overpression or EZH2 silence significantly inhibited HCC cell growth. Overexpression of EZH2 rescued the growth inhibition effect of miR-26a. These findings suggest that an imbalanced double-negative feedback loop between EZH2 and miR-26a exists in HCC cells, which contributes to miR-26a deregulation and regulates tumor cells proliferation.

Expression and DNA methylation alterations of seven cancer-associated 3p genes and their predicted regulator miRNAs (miR-129-2, miR-9-1) in breast and ovarian cancers

The methylation of promoter CpG islands and interactions between microRNAs (miRNAs) and messenger RNAs (mRNAs) of target genes are considered two crucial epigenetic mechanisms for inducing gene and pathway deregulation in tumors. Here, the expression levels of seven cancer-associated 3p genes (RASSF1(isoform A), RARB(isoform 2), SEMA3B, RHOA, GPX1, NKIRAS1, and CHL1) and their predicted regulator miRNAs (miR-129-2, miR-9-1) were analyzed in breast (BC, 40 samples) and ovarian (OC, 14 samples) cancers using RT-PCR and qPCR. We first revealed a negative correlation between the level of the miR-129-2 precursor and RASSF1(A) and GPX1 mRNA levels in BC (Spearman's correlation coefficient (rs) was − 0.26 in both cases). Similar results were observed for the miR-129-2 precursor and the RASSF1(A), GPX1, RARB(2), and CHL1 genes in OC (rs was in the range − 0.48 to − 0.54). Using methylation-specific PCR, a significant correlation was shown between promoter hypermethylation and the down-regulation of the RASSF1(A), GPX1, RARB(2), SEMA3B, MIR-129-2, and MIR-9-1 genes in BC (rs = 0.41 to 0.75) and of the RASSF1(A) gene in OC (rs = 0.67). We first demonstrated a high hypermethylation frequency of MIR-129-2 and SEMA3B (up to 45 to 48%) in both BC (69 samples) and OC (41 samples). Moreover, we observed a positive correlation between the hypermethylation of MIR-129-2 and the up-regulation of the RASSF1(A) and GPX1 genes in BC (rs = 0.38 and 0.42, respectively). QPCR analysis of the expression of RASSF1(A) and mature miR-129-2 in additional BC sample set (24 samples) revealed a negative correlation between them (rs = − 0.41) that strengthened the results obtained during the analysis of miR-129-2 precursor level. In summary, the obtained data indicate the involvement of methylation in the down-regulation of the studied coding and miRNA genes and suggest the involvement of miR-129-2 in the deregulation of RASSF1(A) via a direct interaction or/and mediators in common pathways (according to KEGG, Gene Ontology (FDR < 0.01), and GeneCards data) in the examined gynecological tumors.

Identification of Novel Epigenetic Markers of Prostate Cancer by NotI-Microarray Analysis

A significant need for reliable and accurate cancer diagnostics and prognosis compels the search for novel biomarkers that would be able to discriminate between indolent and aggressive tumors at the early stages of disease. The aim of this work was identification of potential diagnostic biomarkers for characterization of different types of prostate tumors. NotI-microarrays with 180 clones associated with chromosome 3 genes/loci were applied to determine genetic and epigenetic alterations in 33 prostate tumors. For 88 clones, aberrations were detected in more than 10% of tumors. The major types of alterations were DNA methylation and/or deletions. Frequent methylation of the discovered loci was confirmed by bisulfite sequencing on selective sampling of genes: FGF12, GATA2, and LMCD1. Three genes (BHLHE40, BCL6, and ITGA9) were tested for expression level alterations using qPCR, and downregulation associated with hypermethylation was shown in the majority of tumors. Based on these data, we proposed the set of potential biomarkers for detection of prostate cancer and discrimination between prostate tumors with different malignancy and aggressiveness: BHLHE40, FOXP1, LOC285205, ITGA9, CTDSPL, FGF12, LOC440944/SETD5, VHL, CLCN2, OSBPL10/ZNF860, LMCD1, FAM19A4, CAND2, MAP4, KY, and LRRC58. Moreover, we probabilistically estimated putative functional relations between the genes within each set using the network enrichment analysis.

Downregulation of NEK11 is associated with drug resistance in ovarian cancer

NEKs [NIMA (never in mitosis gene A)-related expressed kinase] are involved in ovarian cancer development and progression, while their association with drug resistance is limited, especially NEK11, and its relationship with drug resistance has never been reported. In this study, on the basis of comprehensive bioinformatic analyses, including mRNA expression according to microarray data, protein/gene interaction, protein-small molecule interaction, annotation of biological process and microRNA-mRNA interaction analysis, we revealed that the NEK11 mRNA was significantly downregulated in 586 cases of ovarian serous cystadenocarcinomas and cisplatin-resistant A2780 ovarian cancer cells, and interacted with 22 proteins and 4 small molecules which all were contributed to drug resistance in ovarian cancer. Furthermore, seven cell cycle-related biological processes were annotated with NEK11, drug resistance and ovarian cancer, suggesting that NEK11 potentially was involved in the drug resistance in ovarian cancer via its regulatory roles in the cell cycle. In addition, among the eight microRNAs predicted to be most strongly targeting NEK11, the majority were involved in drug resistance in ovarian and other cancers. All those results provide a very strong possibility that the notable downregulation of NEK11 in cisplatin-resistant ovarian cancer cells was involved in drug resistance, via its interactions with drug resistance-related genes, proteins, small molecules, microRNAs and biological processes, particularly the cell cycle-related processes. To our knowledge, this is the first report of the association of NEK11 with drug resistance in cancer, and it would pave the way for further investigation of the drug resistance-related functions of this gene.

Epigenetic alterations of chromosome 3 revealed by NotI-microarrays in clear cell renal cell carcinoma

This study aimed to clarify epigenetic and genetic alterations that occur during renal carcinogenesis. The original method includes chromosome 3 specific NotI-microarrays containing 180 NotI-clones associated with 188 genes for hybridization with 23 paired normal/tumor DNA samples of primary clear cell renal cell carcinomas (ccRCC). Twenty-two genes showed methylation and/or deletion in 17–57% of tumors. These genes include tumor suppressors or candidates (VHL, CTDSPL, LRRC3B, ALDH1L1, and EPHB1) and genes that were not previously considered as cancer-associated (e.g., LRRN1, GORASP1, FGD5, and PLCL2). Bisulfite sequencing analysis confirmed methylation as a frequent event in ccRCC. A set of six markers (NKIRAS1/RPL15, LRRN1, LRRC3B, CTDSPL, GORASP1/TTC21A, and VHL) was suggested for ccRCC detection in renal biopsies. The mRNA level decrease was shown for 6 NotI-associated genes in ccRCC using quantitative PCR: LRRN1, GORASP1, FOXP1, FGD5, PLCL2, and ALDH1L1. The majority of examined genes showed distinct expression profiles in ccRCC and papillary RCC. The strongest extent and frequency of downregulation were shown for ALDH1L1 gene both in ccRCC and papillary RCC. Moreover, the extent of ALDH1L1 mRNA level decrease was more pronounced in both histological types of RCC stage III compared with stages I and II (P = 0.03). The same was observed for FGD5 gene in ccRCC (P < 0.06).

Dedicated to thememory of Eugene R. Zabarovsky

BRCA1 promoter hypermethylation in sporadic epithelial ovarian carcinoma: Association with low expression of BRCA1, improved survival and co-expression of DNA methyltransferases

The breast cancer susceptibility gene 1 (BRCA1) inactivation in sporadic epithelial ovarian carcinoma (EOC) is common and low BRCA1 expression is associated with promoter hypermethylation. The clinical validation of BRCA1 methylation as a prognostic marker in EOC remains unresolved. The aim of the present study was to determine the aberrant promoter methylation of BRCA1 in benign and malignant ovarian tumor tissues, to establish the association with the clinicopathological significance and the prognostic value. Additionally, the contribution of DNA methyltransferase (DNMT) expression to BRCA1 promoter hypermethylation was determined. The rate of BRCA1 methylation was observed to be 35.2% (50/142) in the EOCs; however, no methylation (0/32) was observed in the benign tumors. BRCA1 methylation was significantly associated with the downregulation of BRCA1 expression (P<0.001) and the frequency of BRCA1 methylation was greater in the carcinomas of patients whose tumor was bilateral than that of patients with a unilateral carcinoma (P=0.015). BRCA1 methylation was significantly associated with the preoperative serum carbohydrate antigen-125 level (P=0.013), improved overall survival (P=0.005) and disease-free survival (P=0.007). In addition, a significant correlation was observed between the co-expression of DNMTs and the methylation status of BRCA1. Thus, the present study provided support for BRCA1 promoter hypermethylation as a prognostic marker for survival in sporadic EOC, and co-expression of DNMTs was observed to contribute to BRCA1 promoter hypermethylation.

Novel tumor suppressor candidates on chromosome 3 revealed by NotI-microarrays in cervical cancer

Genetic and epigenetic alterations in cervical carcinomas were investigated using NotI-microarrays containing 180 cloned sequences flanking all NotI-sites associated with genes on chromosome 3. In total, 48 paired normal/tumor DNA samples, specifically enriched in NotI-sites, were hybridized to NotI-microarrays. Thirty genes, including tumor suppressors or candidates (for example, VHL, RBSP3/CTDSPL, ITGA9, LRRC3B, ALDH1L1, EPHB1) and genes previously unknown as cancer-associated (ABHD5, C3orf77, PRL32, LOC285375, FGD5 and others), showed methylation/deletion in 21-44% of tumors. The genes were more frequently altered in squamous cell carcinomas (SCC) than in adenocarcinomas (ADC, p<0.01). A set of seven potential markers (LRRN1, PRICKLE2, VHL, BHLHE40, RBSP3, CGGBP1 and SOX14) is promising for discrimination of ADC and SCC. Alterations of more than 20 genes simultaneously were revealed in 23% of SCC. Bisulfite sequencing analysis confirmed methylation as a frequent event in SCC. High down-regulation frequency was shown for RBSP3, ITGA9, VILL, APRG1/C3orf35 and RASSF1 (isoform A) genes (3p21.3 locus) in SCC. Both frequency and extent of RASSF1A and RBSP3 mRNA level decrease were more pronounced in tumors with lymph node metastases compared with non-metastatic ones (p ≤ 0.05). We confirmed by bisulfite sequencing that RASSF1 promoter methylation was a rare event in SCC and, for the first time, demonstrated RASSF1A down-regulation at both the mRNA and protein levels without promoter methylation in tumors of this histological type. Thus, our data revealed novel tumor suppressor candidates located on chromosome 3 and a frequent loss of epigenetic stability of 3p21.3 locus in combination with down-regulation of genes in cervical cancer.