Identification of differentially-expressed genes by DNA methylation in cervical cancer

Non-Coding RNAs of Mitochondrial Origin: Roles in Cell Division and Implications in Cancer

Non-coding RNAs (ncRNAs) are a heterogeneous group, in terms of structure and sequence length, consisting of RNA molecules that do not code for proteins. These ncRNAs have a central role in the regulation of gene expression and are virtually involved in every process analyzed, ensuring cellular homeostasis. Although, over the years, much research has focused on the characterization of non-coding transcripts of nuclear origin, improved bioinformatic tools and next-generation sequencing (NGS) platforms have allowed the identification of hundreds of ncRNAs transcribed from the mitochondrial genome (mt-ncRNA), including long non-coding RNA (lncRNA), circular RNA (circRNA), and microRNA (miR). Mt-ncRNAs have been described in diverse cellular processes such as mitochondrial proteome homeostasis and retrograde signaling; however, the function of the majority of mt-ncRNAs remains unknown. This review focuses on a subgroup of human mt-ncRNAs whose dysfunction is associated with both failures in cell cycle regulation, leading to defects in cell growth, cell proliferation, and apoptosis, and the development of tumor hallmarks, such as cell migration and metastasis formation, thus contributing to carcinogenesis and tumor development. Here we provide an overview of the mt-ncRNAs/cancer relationship that could help the future development of new biomedical applications in the field of oncology.

A New Vista of Aldehyde Dehydrogenase 1A3 (ALDH1A3): New Specific Inhibitors and Activity-Based Probes Targeting ALDH1A3 Dependent Pathways in Glioblastoma, Mesothelioma and Other Cancers

Aldehyde dehydrogenases of the subfamily 1A (ALDH1A) are enzymes necessary for the oxidation of all-trans or 9-cis retinal to retinoic acid (RA). Retinoic acid and its derivatives are important for normal development and maintenance of epithelia, reproduction, memory, and immune function in adults. Moreover, in recent years, it has been demonstrated that ALDH1A members are also expressed and functional in several human cancers where their role is not limited to the synthesis of RA. Here, we review the current knowledge about ALDH1A3, one of the 1A isoforms, in cancers with an emphasis on two of the deadliest tumors that affect humans: glioblastoma multiforme and mesothelioma. In both tumors, ALDH1A3 is considered a negative prognostic factor, and its level correlates with excessive proliferation, chemoresistance, and invasiveness. We also review the recent attempts to develop both ALDH1A3-selective inhibitors for cancer therapy and ALDH1A3-specific fluorescent substrates for fluorescence-guided tumor resection.

AhRR and PPP1R3C: Potential Prognostic Biomarkers for Serous Ovarian Cancer

The lack of effective screening and successful treatment contributes to high ovarian cancer mortality, making it the second most common cause of gynecologic cancer death. Development of chemoresistance in up to 75% of patients is the cause of a poor treatment response and reduced survival. Therefore, identifying potential and effective biomarkers for its diagnosis and prognosis is a strong critical need. Copy number alterations are frequent in cancer, and relevant for molecular tumor stratification and patients' prognoses. In this study, array-CGH analysis was performed in three cell lines and derived cancer stem cells (CSCs) to identify genes potentially predictive for ovarian cancer patients' prognoses. Bioinformatic analyses of genes involved in copy number gains revealed that AhRR and PPP1R3C expression negatively correlated with ovarian cancer patients' overall and progression-free survival. These results, together with a significant association between AhRR and PPP1R3C expression and ovarian cancer stemness markers, suggested their potential role in CSCs. Furthermore, AhRR and PPP1R3C's increased expression was maintained in some CSC subpopulations, reinforcing their potential role in ovarian cancer. In conclusion, we reported for the first time, to the best of our knowledge, a prognostic role of AhRR and PPP1R3C expression in serous ovarian cancer.

Downregulation of CAMK2N1 due to DNA Hypermethylation Mediated by DNMT1 that Promotes the Progression of Prostate Cancer

Calcium/calmodulin-dependentprotein kinase II inhibitor I (CAMK2N1) as one of the tumor suppressor genes is significantly downregulated in prostate cancer (PCa). Reduced expression of CAMK2N1 is positively correlated with PCa progression. However, the mechanisms of CAMK2N1 downregulation in PCa are still unclear. The promoter region of CAMK2N1 contains a large number of CG loci, providing the possibility for DNA methylation. Consequently, we hypothesized that DNA methylation can result in the reduced expression of CAMK2N1 in PCa. In the presented study, the DNA methylation level of CAMK2N1 in prostate cells and clinical specimens was determined by bisulfite sequencing (BS), pyrosequencing, and in silico analysis. Results showed that CAMK2N1 was highly methylated in PCa cells and tissues compared to normal prostate epithelial cells and nonmalignant prostate tissues, which was associated with the clinicopathological characteristics in PCa patients. Afterwards, we explored the expression of CAMK2N1 and its DNA methylation level by qRT-PCR, western blot, BS, and methylation-specific PCR in PCa cells after 5-Aza-CdR treatment or DNMT1 genetic modification, which demonstrated that the reduced expression of CAMK2N1 can be restored by 5-Aza-CdR treatment via demethylation. Moreover, DNMT1 formed a positive feedback loop with CAMK2N1 in PCa cells. The expression of CAMK2N1 was downregulated by DNMT1-mediated DNA methylation, which reversely induced DNMT1 expression through activating AKT or ERK signaling pathway. Finally, functional assays including wound healing, invasion, and migration assay, as well as the xenograft model in nude mice indicated that CAMK2N1 inhibited the invasion, migration, and proliferation of PCa cells and these effects were reversed by DNMT1 overexpression. In conclusion, DNMT1-mediated hypermethylation of CAMK2N1 not only downregulates the gene expression but also promotes the progression of PCa.

MiR-4461 Inhibits Tumorigenesis of Renal Cell Carcinoma by Targeting PPP1R3C

Background: Renal cell carcinoma (RCC) is one of the most common and malignant tumors in the urinary system. The aim of this research was to investigate the mechanism and clinical significance of miR-4461 in the RCC progression. Materials and Methods: Twenty-eight (28) paired RCC tissue samples and adjacent nontumor tissue samples, as well as RCC cell lines were used to measure the expression of miR-4461 and protein phosphatase 1 regulatory subunit 3C (PPP1R3C) transcript by real-time quantitative PCR. The target relationship between miR-4461 and PPP1R3C was predicted by TargetScan and further verified by dual-luciferase reporter gene assay and RNA pull-down assay. Cell Counting Kit-8 (CCK-8) assay and BrdU ELISA assay were performed to measure RCC cell viability and proliferation. In addition, caspase-3 activity assay and cell adhesion assay were implemented to measure RCC cell apoptosis and adhesion. Results: MiR-4461 was lowly expressed both in RCC tissues and cells, while upregulated PPP1R3C was tested in RCC tissues and cells. In addition, miR-4461 was validated to directly target PPP1R3C, thereby negatively regulating PPP1R3C. Particularly, miR-4461 exerted a clear inhibitory effect on the malignant phenotypes of RCC cells by binding and inhibiting PPP1R3C. Conclusion: MiR-4461, which served as a tumor suppressor, inhibited RCC progression by targeting and downregulating PPP1R3C.

Cervical cancer subtypes harbouring integrated and/or episomal HPV16 portray distinct molecular phenotypes based on transcriptome profiling of mRNAs and miRNAs

Heterogeneity in cervical cancers (CaCx) in terms of HPV16 physical status prompted us to investigate the mRNA and miRNA signatures among the different categories of CaCx samples. We performed microarray-based mRNA expression profiling and quantitative real-time PCR-based expression analysis of some prioritised miRNAs implicated in cancer-related pathways among various categories of cervical samples. Such samples included HPV16-positive CaCx cases that harboured either purely integrated HPV16 genomes (integrated) and those that harboured episomal viral genomes, either pure or concomitant with integrated viral genomes (episomal), which were compared with normal cervical samples that were either HPV negative or positive for HPV16. The mRNA expression profile differed characteristically between integrated and episomal CaCx cases for enriched biological pathways. miRNA expression profiles also differed among CaCx cases compared with controls (upregulation—miR-21, miR-16, miR-205, miR-323; downregulation—miR-143, miR-196b, miR-203, miR-34a; progressive upregulation—miR-21 and progressive downregulation—miR-143, miR-34a, miR-196b and miR-203) in the order of HPV-negative controls, HPV16-positive non-malignant samples and HPV16-positive CaCx cases. miR-200a was upregulated in HPV16-positive cervical tissues irrespective of histopathological status. Expression of majority of the predicted target genes was negatively correlated with their corresponding miRNAs, irrespective of the CaCx subtypes. E7 mRNA expression correlated positively with miR-323 expression among episomal cases and miR-203, among integrated cases. miR-181c expression was downregulated only among the episomal CaCx cases and negatively correlated with protein coding transcript of the proliferative target gene, CKS1B of the significantly enriched “G2/M DNA Damage Checkpoint Regulation” pathway among CaCx cases. Thus, the two CaCx subtypes are distinct entities at the molecular level, which could be differentially targeted for therapy. In fact, availability of a small molecule inhibitor of CKS1B, suggests that drugging CKS1B could be a potential avenue of treating the large majority of CaCx cases harbouring episomal HPV16.

Methylation and expression of PTPN22 in esophageal squamous cell carcinoma

Esophageal squamous cell carcinoma (ESCC) is a fatal disease contributed by both genetic and epigenetic factors. The epigenetic alteration of protein tyrosine phosphatase non-receptor type 22 (PTPN22) and its clinical significance in ESCC were still not yet clarified. A quantitative methylation study of PTPN22 and its expression were conducted in 121 and 31 paired tumor and adjacent normal tissue (ANT), respectively. Moreover, the association between PTPN22 methylation and clinicopathological parameters was evaluated. We found that the methylation level of PTPN22 was significantly elevated in tumor tissues (66.3%) relative to ANT (62.1%) (p=0.005). The methylation level of non-smoking ANT (59.1%) was significant lower than smoking ESCC tissue (65.8%) (p=0.03); similarly, the methylation levels in ANT with no lymph node invasion (57.6%) were significant lower than tumor tissues with lymph node invasion (67.5%) (p=0.001). PTPN22 expression in ESCC was lower than normal tissues, however the difference was not statistically significant (p=0.55). Lower expression was more frequently occurred in N_1-3 and III stage patients, while higher expression was more likely to occur in N₀ and I-II stage patients. Lower expression of PTPN22 was associated with poor overall survival (p=0.04). Taken together, PTPN22 was hypermethylationed in ESCC. Hypermethylation was associated with lymph node invasion. The PTPN22 expression may act as a prognostic biomarker to identify patients at risk of high grade.

pETM: a penalized Exponential Tilt Model for analysis of correlated high-dimensional DNA methylation data

Motivation

DNA methylation plays an important role in many biological processes and cancer progression. Recent studies have found that there are also differences in methylation variations in different groups other than differences in methylation means. Several methods have been developed that consider both mean and variance signals in order to improve statistical power of detecting differentially methylated loci. Moreover, as methylation levels of neighboring CpG sites are known to be strongly correlated, methods that incorporate correlations have also been developed. We previously developed a network-based penalized logistic regression for correlated methylation data, but only focusing on mean signals. We have also developed a generalized exponential tilt model that captures both mean and variance signals but only examining one CpG site at a time.

Results

In this article, we proposed a penalized Exponential Tilt Model (pETM) using network-based regularization that captures both mean and variance signals in DNA methylation data and takes into account the correlations among nearby CpG sites. By combining the strength of the two models we previously developed, we demonstrated the superior power and better performance of the pETM method through simulations and the applications to the 450K DNA methylation array data of the four breast invasive carcinoma cancer subtypes from The Cancer Genome Atlas (TCGA) project. The developed pETM method identifies many cancer-related methylation loci that were missed by our previously developed method that considers correlations among nearby methylation loci but not variance signals.

Availability and Implementation

The R package 'pETM' is publicly available through CRAN: http://cran.r-project.org .

Contact

sw2206@columbia.edu.

Supplementary information

Supplementary data are available at Bioinformatics online.

Investigation of differentially-expressed microRNAs and genes in cervical cancer using an integrated bioinformatics analysis

Cervical cancer is one of the most common types of cancer among women worldwide. In order to identify the microRNAs (miRNAs/miRs) and mRNAs associated with the carcinogenesis of cervical cancer, and to investigate the molecular mechanisms of cervical cancer, an miRNA microarray, GSE30656, and 3 mRNA microarrays, GSE63514, GSE39001 and GSE9750, for cervical cancer were retrieved from Gene Expression Omnibus. These datasets were analyzed in order to obtain differentially-expressed genes (DEGs) and miRNAs using the GEO2R tool. Gene Ontology (GO) and pathway enrichment analysis for DEGs were performed using the Database for Annotation, Visualization and Integrated Discovery. Protein-protein interaction (PPI) analysis for DEGs was conducted using The Search Tool for the Retrieval of Interacting Genes software and visualized using Cytoscape, followed by hub gene identification, and biological process and pathway enrichment analysis of the module selected from the PPI network using the Molecular Complex Detection plugin. In addition, miRecords was applied to predict the targets of differentially-expressed miRNAs. A total of 44 DEGs and 15 differentially-expressed miRNAs were identified. These DEGs were mainly enriched in GO terms associated with the cell cycle. In the PPI network, cyclin-dependent kinase 1, topoisomerase DNA IIα, aurora kinase A (AURKA) and minichromosome maintenance complex component 2 (MCM2) had higher degrees of connectivity. A significant module was detected from the PPI network. AURKA, MCM2 and kinesin family member 20A exhibited higher degrees in this module, while the genes in the module were mainly involved in the cell cycle and the DNA replication pathway. In addition, estrogen receptor 1 was predicted as the potential target of 13 miRNAs. A total of 10 DEGs were identified as potential targets of miR-203. In conclusion, the results indicated that microarray dataset analysis may provide a useful method for the identification of key genes and patterns to successfully identify determinants of the carcinogenesis of cervical cancer. The functional studies of candidate genes and miRNAs from these databases may lead to an increased understanding of the development of cervical cancer.

JunD/AP-1 Antagonizes the Induction of DAPK1 To Promote the Survival of v-Src-Transformed Cells

The increase in AP-1 activity is a hallmark of cell transformation by tyrosine kinases. Previously, we reported that blocking AP-1 using the c-Jun dominant negative mutant TAM67 induced senescence, adipogenesis, or apoptosis in v-Src-transformed chicken embryo fibroblasts (CEFs) whereas inhibition of JunD by short hairpin RNA (shRNA) specifically induced apoptosis. To investigate the role of AP-1 in Src-mediated transformation, we undertook a gene profiling study to characterize the transcriptomes of v-Src-transformed CEFs expressing either TAM67 or the JunD shRNA. Our study revealed a cluster of 18 probe sets upregulated exclusively in response to AP-1/JunD impairment and v-Src transformation. Four of these probe sets correspond to genes involved in the interferon pathway. One gene in particular, death-associated protein kinase 1 (DAPK1), is a C/EBPβ-regulated mediator of apoptosis in gamma interferon (IFN-γ)-induced cell death. Here, we show that inhibition of DAPK1 abrogates cell death in v-Src-transformed cells expressing the JunD shRNA. Chromatin immunoprecipitation data indicated that C/EBPβ was recruited to the DAPK1 promoter while the expression of a dominant negative mutant of C/EBPβ abrogated the induction of DAPK1 in response to the inhibition of AP-1. In contrast, as determined by chromatin immunoprecipitation (ChIP) assays, JunD was not detected on the DAPK1 promoter under any conditions, suggesting that JunD promotes survival by indirectly antagonizing the expression of DAPK1 in v-Src transformed cells.

IMPORTANCE

Transformation by the v-Src oncoprotein causes extensive changes in gene expression in primary cells such as chicken embryo fibroblasts. These changes, determining the properties of transformed cells, are controlled in part at the transcriptional level. Much attention has been devoted to transcription factors such as AP-1 and NF-κB and the control of genes associated with a more aggressive phenotype. In this report, we describe a novel mechanism of action determined by the JunD component of AP-1, a factor enhancing cell survival in v-Src-transformed cells. We show that the loss of JunD results in the aberrant activation of a genetic program leading to cell death. This program requires the activation of the tumor suppressor death-associated protein kinase 1 (DAPK1). Since DAPK1 is phosphorylated and inhibited by v-Src, these results highlight the importance of this kinase and the multiple mechanisms controlled by v-Src to antagonize the tumor suppressor function of DAPK1.

Breast cancer upregulating genes in stromal cells by LINE-1 hypermethylation and micrometastatic detection

Changes in the methylation level of genes containing LINE-1 alter host gene regulation.

AIM

This study demonstrates that paracrine signaling of breast cancer influences the epigenetic regulation of stromal cells.

METHODS

We proved in vitro and in vivo breast cancer promoted LINE-1 methylation exists exclusively in female stromal cells.

RESULTS

Genes containing LINE-1 of breast cancer stromal cells were upregulated. Furthermore, one of the genes, MUC-1, was demonstrated to have expression in plasma cells from the lymph nodes of patients with lymph node metastasis or micrometastasis.

CONCLUSION

Breast cancer sends a paracrine signal to stroma cells causing LINE-1 epigenetic regulation. Moreover, the regulated genes in stroma cells are potential biomarkers for detecting breast cancer micrometastasis.

BLCAP arrests G₁/S checkpoint and induces apoptosis through downregulation of pRb1 in HeLa cells

BLCAP (bladder cancer-associated protein) gene exhibited tumor suppressor function in different tumors and is regarded as a candidate tumor suppressor gene; however, the mechanism by which BLCAP exerts its function remains elusive. This study investigated the functional association between BLCAP and proliferation or apoptosis in cervical cancer cells, to identify the functional motifs of BLCAP. The BLCAP-shRNA expression vector based on pRNA-U6.1/Hygro plasmid was used to specifically inhibit BLCAP activity in HeLa cells. The optimal shRNA plasmid was selected to knock down BLCAP expression and the biological effects were investigated. The effects on cell cycle and apoptosis were detected by flow cytometric or Annexin V-FITC staining analysis. The gene expression profiles of HeLa cells transfected with blcap-wt and BLCAP-shRNA were analyzed using human signal pathway gene Oligochips. The levels of protein expression and interaction of BLCAP with Rb1 proteins were determined by western blotting and Co-IP assays. The site-specific mutagenesis assay was used to identify amino acid residues important for BLCAP. Significantly differentially expressed genes were found by gene Oligo chips analysis. These genes were all correlated with proliferation, cell cycle and apoptosis. The results of western blotting and Co-IP assays confirmed that overexpression of BLCAP could interact with Rb1 and inhibit Rb1 phosphorylation. Further investigation revealed that SAXX mutation in the key regions of BLCAP suppressed the function of BLCAP and significantly increased the level of phosphorylated Rb1 protein. Here our findings suggested that the functional association of BLCAP and Rb1 might play important roles in proliferation and apoptosis of HeLa cells. It suggested that BLCAP could be a novel therapeutic target for cervical cancer.

Aberrant Hypermethylation of SALL3 with HPV Involvement Contributes to the Carcinogenesis of Cervical Cancer

OBJECTIVE

This study aimed to investigate the methylation status of the promoter region of spalt-like transcription factor 3 (SALL3) and the expression of SALL3 in cervical cancer to explore the function of this gene in cervical cancer carcinogenesis.

METHODS

The methylation status of SALL3 was detected by methylation-specific PCR, and SALL3 gene expression was assessed by real-time quantitative PCR in the cervical cancer cell lines, SiHa, HeLa and C33A, as well as in cervical cancer tissue samples (n = 23), matched pericarcinomatous tissue samples (n = 23) and normal cervix tissue samples (n = 17). MTT was used to measure the cell viability and proliferation capacity of SiHa and HeLa cells.

RESULTS

The SALL3 promoter was completely methylated in SiHa cells, unmethylated in C33A cells and partially methylated in HeLa cells. After treatment of SiHa and HeLa cells with 5 μM and 10 μM of 5-Azacytidine (5-Aza), respectively, the methylation level of the SALL3 promoter decreased and observed increase in the degree of unmethylation in a dose-dependent manner. Moreover, the relative expression of SALL3 mRNA increased as the concentration of 5-Aza increased in SiHa (p<0.05) and HeLa (p<0.05) cells. This above-mentioned increase in SALL3 mRNA in SiHa cells was more remarkable than that observed in HeLa cells. Cell proliferation capacity also decreased after administration of 5-Aza to SiHa and HeLa cells (p<0.05). Methylation of the SALL3 promoter was observed in 15 of 23 (65.21%) cervical cancer tissue samples, 15 of 23 (65.21%) matched pericarcinomatous tissue samples and 5 of 17 (29.41%) normal cervical tissue samples (p<0.05). SALL3 mRNA expression was significantly lower in cervical cancer and pericarcinomatous tissues compared with normal cervical tissues (p<0.05). In all cervix tissue samples, HPV infection was positively associated with hypermethylation of the promoter region of SALL3 (p<0.05, r = 0.408), and the expression of SALL3 mRNA in HPV-positive tissues was lower than that in HPV-negative tissues (p<0.05).

CONCLUSION

The aberrant hypermethylation of SALL3 together with HPV involvement inactivated its function as a tumor suppressor and contributed to carcinogenesis in cervical cancer.