Analysis of the mitogen-activated protein kinase kinase 4 (MAP2K4) tumor suppressor gene in ovarian cancer

Genomic analyses of germline and somatic variation in high-grade serous ovarian cancer

BACKGROUND

High-grade serous ovarian cancers (HGSCs) display a high degree of complex genetic alterations. In this study, we identified germline and somatic genetic alterations in HGSC and their association with relapse-free and overall survival. Using a targeted capture of 557 genes involved in DNA damage response and PI3K/AKT/mTOR pathways, we conducted next-generation sequencing of DNA from matched blood and tumor tissue from 71 HGSC participants. In addition, we performed the OncoScan assay on tumor DNA from 61 participants to examine somatic copy number alterations (SCNA).

RESULTS

Approximately one-third of tumors had loss-of-function (LOF) germline (18/71, 25.4%) or somatic (7/71, 9.9%) variants in the DNA homologous recombination repair pathway genes BRCA1, BRCA2, CHEK2, MRE11A, BLM, and PALB2. LOF germline variants also were identified in other Fanconi anemia genes and in MAPK and PI3K/AKT/mTOR pathway genes. Most tumors harbored somatic TP53 variants (65/71, 91.5%). Using the OncoScan assay on tumor DNA from 61 participants, we identified focal homozygous deletions in BRCA1, BRCA2, MAP2K4, PTEN, RB1, SLX4, STK11, CREBBP, and NF1. In total, 38% (27/71) of HGSC patients harbored pathogenic variants in DNA homologous recombination repair genes. For patients with multiple tissues from the primary debulking or from multiple surgeries, the somatic mutations were maintained with few newly acquired point mutations suggesting that tumor evolution was not through somatic mutations. There was a significant association of LOF variants in homologous recombination repair pathway genes and high-amplitude somatic copy number alterations. Using GISTIC analysis, we identified NOTCH3, ZNF536, and PIK3R2 in these regions that were significantly associated with an increase in cancer recurrence and a reduction in overall survival.

CONCLUSIONS

From 71 patients with HGCS, we performed targeted germline and tumor sequencing and provided a comprehensive analysis of these 557 genes. We identified germline and somatic genetic alterations including somatic copy number alterations and analyzed their associations with relapse-free and overall survival. This single-site long-term follow-up study provides additional information on genetic alterations related to occurrence and outcome of HGSC. Our findings suggest that targeted treatments based on both variant and SCNA profile potentially could improve relapse-free and overall survival.

Genomic Analyses of Germline and Somatic Variation in High-Grade Serous Ovarian Cancer

Background High-grade serous ovarian cancers (HGSCs) display a high degree of complex genetic alterations. In this study, we identified germline and somatic genetic alterations in HGSC and their association with relapse-free and overall survival. Using a targeted capture of 577 genes involved in DNA damage response and PI3K/AKT/mTOR pathways, we conducted next-generation sequencing of DNA from matched blood and tumor tissue from 71 HGSC participants. In addition, we performed the OncoScan assay on tumor DNA from 61 participants to examine somatic copy number alterations. Results Approximately one-third of tumors had loss-of-function germline (18/71, 25.4%) or somatic (7/71, 9.9%) variants in the DNA homologous recombination repair pathway genes BRCA1, BRCA2, CHEK2, MRE11A, BLM , and PALB2 . Loss-of-function germline variants also were identified in other Fanconi anemia genes and in MAPK and PI3K/AKT/mTOR pathway genes. Most tumors harbored somatic TP53 variants (65/71, 91.5%). Using the OncoScan assay on tumor DNA from 61 participants, we identified focal homozygous deletions in BRCA1, BRCA2, MAP2K4, PTEN, RB1, SLX4, STK11, CREBBP , and NF1 . In total, 38% (27/71) of HGSC patients harbored pathogenic variants in DNA homologous recombination repair genes. For patients with multiple tissues from the primary debulking or from multiple surgeries, the somatic mutations were maintained with few newly acquired point mutations suggesting that tumor evolution was not through somatic mutations. There was a significant association of loss-of-function variants in homologous recombination repair pathway genes and high-amplitude somatic copy number alterations. Using GISTIC analysis, we identified NOTCH3, ZNF536 , and PIK3R2 in these regions that were significantly associated with an increase in cancer recurrence and a reduction in overall survival. Conclusions From 71 patients with HGCS, we performed targeted germline and tumor sequencing and provided a comprehensive analysis of these 577 genes. We identified germline and somatic genetic alterations including somatic copy number alterations and analyzed their associations with relapse-free and overall survival. This single-site long-term follow-up study provides additional information on genetic alterations related to occurrence and outcome of HGSC. Our findings suggest that targeted treatments based on both variant and SCNA profile potentially could improve relapse-free and overall survival.

The circular RNA expression profile in ovarian serous cystadenocarcinoma reveals a complex circRNA-miRNA regulatory network

Background

Ovarian serous cystadenocarcinoma is one of the most serious gynecological malignancies. Circular RNA (circRNA) is a type of noncoding RNA with a covalently closed continuous loop structure. Abnormal circRNA expression might be associated with tumorigenesis because of its complex biological mechanisms by, for example, functioning as a microRNA (miRNA) sponge. However, the circRNA expression profile in ovarian serous cystadenocarcinoma and their associations with other RNAs have not yet been characterized. The main purpose of this study was to reveal the circRNA expression profile in ovarian serous cystadenocarcinoma.

Methods

We collected six specimens from three patients with ovarian serous cystadenocarcinoma and adjacent normal tissues. After RNA sequencing, we analyzed the expression of circRNAs with relevant mRNAs and miRNAs to characterize potential function.

Results

15,092 unique circRNAs were identified in six specimens. Approximately 46% of these circRNAs were not recorded in public databases. We then reported 353 differentially expressed circRNAs with oncogenes and tumor-suppressor genes. Furthermore, a conjoint analysis with relevant mRNAs revealed consistent changes between circRNAs and their homologous mRNAs. Overall, construction of a circRNA–miRNA network suggested that 4 special circRNAs could be used as potential biomarkers.

Conclusions

Our study revealed the circRNA expression profile in the tissues of patients with ovarian serous cystadenocarcinoma. The differential expression of circRNAs was thought to be associated with ovarian serous cystadenocarcinoma in the enrichment analysis, and co-expression analysis with relevant mRNAs and miRNAs illustrated the latent regulatory network. We also constructed a complex circRNA–miRNA interaction network and then demonstrated the potential function of certain circRNAs to aid future diagnosis and treatment.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12920-021-01132-5.

MAP2K4 interacts with Vimentin to activate the PI3K/AKT pathway and promotes breast cancer pathogenesis

Mitogen-activated protein kinase kinase 4 (MAP2K4) is a member of the mitogen-activated protein kinase (MAPK) activator family. MAPK signaling plays a significant role in cell proliferation, differentiation, transcriptional regulation, and development. However, specific function and mechanism of MAP2K4 in breast cancer have not been clarified. According to our study, overexpressed MAP2K4 in breast cancer cells increased proliferation, migration, and invasion in vivo and in vitro, while MAP2K4 knockdown restored the effects. Subsequent mechanistic analyses demonstrated that MAP2K4 promoted cell proliferation, migration, and invasion by activating phosphoinositide-3-kinase (PI3K)/AKT signaling, the downstream proteins, c-JUN, the G1/S cell cycle, and the epithelial-to-mesenchymal transition (EMT). Meanwhile, MAP2K4 interacted with Vimentin and further propagated the malignant phenotype. Furthermore, patients with high MAP2K4 and Vimentin expression levels had poorer overall survival rates than those with low expression levels of both proteins. Our studies demonstrated that MAP2K4 has the potential to serve as an oncogene in breast cancer and it activates the phosphorylated PI3K/AKT signaling pathway to activate downstream cycle-associated proteins and EMT signals while interacting with Vimentin to promote breast cancer cells proliferation, migration, and invasion. In our study, MAP2K4 and Vimentin co-expression is confirmed to be an unfavorable factor in breast cancer.

Comprehensive network of miRNA-induced intergenic interactions and a biological role of its core in cancer

MicroRNAs (miRNAs) are a family of short noncoding RNAs that posttranscriptionally regulate gene expression and play an important role in multiple cellular processes. A significant percentage of miRNAs are intragenic, which is often functionally related to their host genes playing either antagonistic or synergistic roles. In this study, we constructed and analyzed the entire network of intergenic interactions induced by intragenic miRNAs. We further focused on the core of this network, which was defined as a union of nontrivial strongly connected components, i.e., sets of nodes (genes) mutually connected via directed paths. Both the entire network and its core possessed statistically significant non-random properties. Specifically, genes forming the core had high expression levels and low expression variance. Furthermore, the network core did not split into separate components corresponding to individual signalling or metabolic pathways, but integrated genes involved in key cellular processes, including DNA replication, transcription, protein homeostasis and cell metabolism. We suggest that the network core, consisting of genes mutually regulated by their intragenic miRNAs, could coordinate adjacent pathways or homeostatic control circuits, serving as a horizontal inter-circuit link. Notably, expression patterns of these genes had an efficient prognostic potential for breast and colorectal cancer patients.

Case report: primary acinar cell carcinoma of the liver treated with multimodality therapy

We describe a case of primary acinar cell carcinoma (ACC) originating in the liver in a 54-year-old female, diagnosed following persistent abnormal elevated liver function. Imaging revealed two masses, one dominant lesion in the right hepatic lobe and another in segment IVA. A right hepatectomy was performed to remove the larger lesion, while the mass in segment IVA was unresectable due to its proximity to the left hepatic vein. Immunohistochemical staining showed positivity for trypsin and chymotrypsin. Postoperatively the patient underwent hepatic arterial embolization of the other unresectable lesion followed by FOLFOX chemotherapy. At 20 months from diagnosis the patient is currently under observation with a decreasing necrotic mass and no other disease evident. Based on histology, immunohistochemistry and radiological findings a diagnosis of primary ACC of the liver was made. Genomic assessment of somatic mutations within the patient's tumor was also performed through next generation sequencing and findings were consistent with an acinar malignancy. This case highlights a rare tumor subtype treated with a combination of therapeutic modalities through a multidisciplinary approach.

Exome and transcriptome sequencing identifies loss of PDLIM2 in metastatic colorectal cancers

Background

Understanding the genomic determinants associated with metastasis in colorectal cancers (CRCs) provides crucial clues for improving patient care.

Patients and methods

In this study, we performed whole-exome sequencing as well as RNA sequencing analyses on five pairs of primary and liver metastasized samples from CRC patients together with blood/normal control samples for each pair.

Results

We identified genomic deletions in the region of 8p21-23 (q value <0.01) from analysis of recurrent regions with copy number variations in both primary and matched metastatic lesions. Consistent with this result, we found significantly decreased expression levels of all 12 genes (ADAMDEC1, C8orf80, CLDN23, EPHX2, GFRA2, NEFL, NEFM, PDLIM2, PTK2B, SCARA5, SLC18A1 and STMN4) located within this region (adjusted P<0.01). Notably, the mRNA levels of PDLIM2, a key regulator of well-known cancer-associated genes including the proto-oncogene c-MYC, an early response gene IER3, and regulators of apoptosis such as BCL2, FAS, and FASLG, were highly downregulated in tumors compared to normal tissues.

Conclusion

Taken together, our findings uncovered various genomic alterations potentially leading to metastasis in CRC and provide important insights into the development of potential therapeutic targets for preventing metastatic progression of CRC.

MicroRNA-802 plays a tumour suppressive role in tongue squamous cell carcinoma through directly targeting MAP2K4

OBJECTIVES

Tongue squamous cell carcinoma (TSCC) is the most common oral tumours. MicroRNAs play crucial roles in many cell processes including cell viability, development, apoptosis, migration and invasion. The role of miR-802 in the TSCC is still unknown.

MATERIALS AND METHODS

The miR-802 expression in TSCC tissues and cell lines was determined by quantitative real-time polymerase chain reaction. CCK-8 assay was performed to measure the cell viability, while the cell invasion assay was used to determine the cell invasion. Dual-luciferase reporter and western blot were used to confirm the potential target gene of miR-802.

RESULTS

In our study, we demonstrated that miR-802 expression was downregulated in TSCC tissues and cell lines. Elevated expression of miR-802 suppressed the TSCC cell viability and invasion. Moreover, enforced expression of miR-802 increased the expression of E-cadherin, while suppressed the expression of N-cadherin, Snail and Vimentin in the TSCC cell. In addition, we identified the mitogen-activated protein kinase 4 (MAP2K4) as a direct target gene of miR-802 in the TSCC cell. We also demonstrated that the expression of MAP2K4 was higher in the TSCC tissues than that in the adjacent normal tissues. Furthermore, the expression level of MAP2K4 was inversely associated with the expression of miR-802 in TSCC tissues. We also demonstrated that the MAP2K4 expression was upregulated in TSCC cell lines. Elevated expression of miR-802 inhibited TSCC cell viability and invasion through inhibiting MAP2K4 expression.

CONCLUSIONS

Our data revealed that miR-802 played as a tumour suppressor gene and might act as a therapeutic target in TSCC patients.

miR-141 promotes colon cancer cell proliferation by inhibiting MAP2K4

MicroRNAs (miRNAs or miRs) can function as tumor-suppressor or oncogenic genes. Upregulation of miRNA-141 has been frequently observed in colorectal cancer (CRC) samples. The experimentally observed targets of miR-141 include the tumor-suppressor gene mitogen-activated protein kinase kinase 4 (MAP2K4). The aim of the present study was to investigate the role of miR-141 in the proliferation of colonic cancer. Western blotting, immunohistochemistry and reverse transcription-quantitative polymerase chain reaction were used to detect the expression levels of miR-141 and MAP2K4 in colonic adenocarcinoma (CAC) and adjacent non-cancerous (NC) tissue samples, as well as in human CAC cell lines (HT29, T94 and LS174). MTT assay was used to investigate the proliferation and apoptosis of these three cell lines. The expression levels of miR-141 were significantly upregulated in clinical samples of CAC, compared with adjacent NC tissues. By contrast, MAP2K4 was downregulated in CAC. The in vitro assays demonstrated that overexpression of miR-141 resulted in cell proliferation of CAC by inhibiting MAP2K4 activity. Our study suggests that targeting the miR-141-MAP2K4 signaling pathway may represent a novel approach for the treatment of CRC.

An exact algorithm for finding cancer driver somatic genome alterations: the weighted mutually exclusive maximum set cover problem

Background

The mutual exclusivity of somatic genome alterations (SGAs), such as somatic mutations and copy number alterations, is an important observation of tumors and is widely used to search for cancer signaling pathways or SGAs related to tumor development. However, one problem with current methods that use mutual exclusivity is that they are not signal-based; another problem is that they use heuristic algorithms to handle the NP-hard problems, which cannot guarantee to find the optimal solutions of their models.

Method

In this study, we propose a novel signal-based method that utilizes the intrinsic relationship between SGAs on signaling pathways and expression changes of downstream genes regulated by pathways to identify cancer signaling pathways using the mutually exclusive property. We also present a relatively efficient exact algorithm that can guarantee to obtain the optimal solution of the new computational model.

Results

We have applied our new model and exact algorithm to the breast cancer data. The results reveal that our new approach increases the capability of finding better solutions in the application of cancer research. Our new exact algorithm has a time complexity of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$O^{*}(1.325^{m})$$\end{document}O∗(1.325m)(Note: Following the recent convention, we use a star * to represent that the polynomial part of the time complexity is neglected), which has solved the NP-hard problem of our model efficiently.

Conclusion

Our new method and algorithm can discover the true causes behind the phenotypes, such as what SGA events lead to abnormality of the cell cycle or make the cell metastasis lose control in tumors; thus, it identifies the target candidates for precision (or target) therapeutics.

Electronic supplementary material

The online version of this article (doi:10.1186/s13015-016-0073-9) contains supplementary material, which is available to authorized users.

Registered report: Diverse somatic mutation patterns and pathway alterations in human cancers

The Reproducibility Project: Cancer Biology seeks to address growing concerns about reproducibility in scientific research by conducting replications of selected experiments from a number of high-profile papers in the field of cancer biology. The papers, which were published between 2010 and 2012, were selected on the basis of citations and Altmetric scores (Errington et al., 2014). This Registered Report describes the proposed replication plan of key experiments from "Diverse somatic mutation patterns and pathway alterations in human cancers" by Kan and colleagues published in Nature in 2010 (Kan et al., 2010). The experiments to be replicated are those reported in Figures 3D-F and 4C-F. Kan and colleagues utilized mismatch repair detection (MRD) technology to identify somatic mutations in primary human tumor samples and identified a previously uncharacterized arginine 243 to histidine (R243H) mutation in the G-protein α subunit GNAO1 in breast carcinoma tissue. In Figures 3D-F, Kan and colleagues demonstrated that stable expression of mutant GNAO1(R243D) conferred a significant growth advantage in human mammary epithelial cells, confirming the oncogenic potential of this mutation. Similarly, expression of variants with somatic mutations in MAP2K4, a JNK pathway kinase (shown in Figures 4C-E) resulted in a significant increase in anchorage-independent growth. Interestingly, these mutants exhibited reduced kinase activity compared to wild type MAP2K4, indicating these mutations impose a dominant-negative influence to promote growth (Figure 4F). The Reproducibility Project: Cancer Biology is a collaboration between the Center for Open Science and Science Exchange and the results of the replications will be published in eLife.

Revealing the Complexity of Breast Cancer by Next Generation Sequencing

Over the last few years the increasing usage of "-omic" platforms, supported by next-generation sequencing, in the analysis of breast cancer samples has tremendously advanced our understanding of the disease. New driver and passenger mutations, rare chromosomal rearrangements and other genomic aberrations identified by whole genome and exome sequencing are providing missing pieces of the genomic architecture of breast cancer. High resolution maps of breast cancer methylomes and sequencing of the miRNA microworld are beginning to paint the epigenomic landscape of the disease. Transcriptomic profiling is giving us a glimpse into the gene regulatory networks that govern the fate of the breast cancer cell. At the same time, integrative analysis of sequencing data confirms an extensive intertumor and intratumor heterogeneity and plasticity in breast cancer arguing for a new approach to the problem. In this review, we report on the latest findings on the molecular characterization of breast cancer using NGS technologies, and we discuss their potential implications for the improvement of existing therapies.

Loss of heterozygosity: what is it good for?

BACKGROUND

Loss of heterozygosity (LOH) is a common genetic event in cancer development, and is known to be involved in the somatic loss of wild-type alleles in many inherited cancer syndromes. The wider involvement of LOH in cancer is assumed to relate to unmasking a somatically mutated tumour suppressor gene through loss of the wild type allele.

METHODS

We analysed 86 ovarian carcinomas for mutations in 980 genes selected on the basis of their location in common regions of LOH.

RESULTS

We identified 36 significantly mutated genes, but these could only partly account for the quanta of LOH in the samples. Using our own and TCGA data we then evaluated five possible models to explain the selection for non-random accumulation of LOH in ovarian cancer genomes: 1. Classic two-hit hypothesis: high frequency biallelic genetic inactivation of tumour suppressor genes. 2. Epigenetic two-hit hypothesis: biallelic inactivation through methylation and LOH. 3. Multiple alternate-gene biallelic inactivation: low frequency gene disruption. 4. Haplo-insufficiency: Single copy gene disruption. 5. Modified two-hit hypothesis: reduction to homozygosity of low penetrance germline predisposition alleles. We determined that while high-frequency biallelic gene inactivation under model 1 is rare, regions of LOH (particularly copy-number neutral LOH) are enriched for deleterious mutations and increased promoter methylation, while copy-number loss LOH regions are likely to contain under-expressed genes suggestive of haploinsufficiency. Reduction to homozygosity of cancer predisposition SNPs may also play a minor role.

CONCLUSION

It is likely that selection for regions of LOH depends on its effect on multiple genes. Selection for copy number neutral LOH may better fit the classic two-hit model whereas selection for copy number loss may be attributed to its effect on multi-gene haploinsufficiency. LOH mapping alone is unlikely to be successful in identifying novel tumour suppressor genes; a combined approach may be more effective.

Evaluation and integration of cancer gene classifiers: identification and ranking of plausible drivers

The number of mutated genes in cancer cells is far larger than the number of mutations that drive cancer. The difficulty this creates for identifying relevant alterations has stimulated the development of various computational approaches to distinguishing drivers from bystanders. We develop and apply an ensemble classifier (EC) machine learning method, which integrates 10 classifiers that are publically available, and apply it to breast and ovarian cancer. In particular we find the following: (1) Using both standard and non-standard metrics, EC almost always outperforms single method classifiers, often by wide margins. (2) Of the 50 highest ranked genes for breast (ovarian) cancer, 34 (30) are associated with other cancers in either the OMIM, CGC or NCG database (P < 10⁻²²). (3) Another 10, for both breast and ovarian cancer, have been identified by GWAS studies. (4) Several of the remaining genes--including a protein kinase that regulates the Fra-1 transcription factor which is overexpressed in ER negative breast cancer cells; and Fyn, which is overexpressed in pancreatic and prostate cancer, among others--are biologically plausible. Biological implications are briefly discussed. Source codes and detailed results are available at http://www.visantnet.org/misi/driver_integration.zip.

Mining TCGA data using Boolean implications

Boolean implications (if-then rules) provide a conceptually simple, uniform and highly scalable way to find associations between pairs of random variables. In this paper, we propose to use Boolean implications to find relationships between variables of different data types (mutation, copy number alteration, DNA methylation and gene expression) from the glioblastoma (GBM) and ovarian serous cystadenoma (OV) data sets from The Cancer Genome Atlas (TCGA). We find hundreds of thousands of Boolean implications from these data sets. A direct comparison of the relationships found by Boolean implications and those found by commonly used methods for mining associations show that existing methods would miss relationships found by Boolean implications. Furthermore, many relationships exposed by Boolean implications reflect important aspects of cancer biology. Examples of our findings include cis relationships between copy number alteration, DNA methylation and expression of genes, a new hierarchy of mutations and recurrent copy number alterations, loss-of-heterozygosity of well-known tumor suppressors, and the hypermethylation phenotype associated with IDH1 mutations in GBM. The Boolean implication results used in the paper can be accessed at http://crookneck.stanford.edu/microarray/TCGANetworks/.