Network-Based Integration of GWAS and Gene Expression Identifies a HOX-Centric Network Associated with Serous Ovarian Cancer Risk

High expression of homeobox B2 predicts poor survival of colon adenocarcinoma by enhancing tumor proliferation and invasion

OBJECTIVE

Homeobox B2 (HOXB2) is known to be correlated with the development and prognosis of various cancers. However, its role in colon cancer remains unclear.

AIMS

In this study, we explored the prognostic value of HOXB2 in colon adenocarcinoma (COAD).

METHODS

A total of 264 colon adenocarcinoma cases were retrospectively enrolled to evaluate HOXB2 expression and clinical significance. Chi-square test was applied to identify relationship between clinical features and HOXB2 expression. The effect of HOXB2 expression and clinical features on the survival of COAD patients was evaluated using Kaplan-Meier and Cox regression analyses. Cellular assays and mice models were conducted to validate the tumor-related role of HOXB2 in COAD.

RESULTS

Higher expression of HOXB2 in COAD tissues was significantly associated with tumor size, invasion depth, and lymph node metastasis (all P < 0.05). Univariate and multivariate analyses showed that high expression of HOXB2 was significantly correlated with a poor overall survival. In vitro cellular assays combined with knockdown strategies demonstrated that HOXB2 can promote tumor proliferation and invasion of COAD, which was further confirmed by in vivo xenograft experiments.

CONCLUSIONS

HOXB2 may be a valuable biomarker and potential therapeutic target for the treatment of COAD.

DNA methylation QTL mapping across diverse human tissues provides molecular links between genetic variation and complex traits

Studies of DNA methylation (DNAm) in solid human tissues are relatively scarce; tissue-specific characterization of DNAm is needed to understand its role in gene regulation and its relevance to complex traits. We generated array-based DNAm profiles for 987 human samples from the Genotype-Tissue Expression (GTEx) project, representing 9 tissue types and 424 subjects. We characterized methylome and transcriptome correlations (eQTMs), genetic regulation in cis (mQTLs and eQTLs) across tissues and e/mQTLs links to complex traits. We identified mQTLs for 286,152 CpG sites, many of which (>5%) show tissue specificity, and mQTL colocalizations with 2,254 distinct GWAS hits across 83 traits. For 91% of these loci, a candidate gene link was identified by integration of functional maps, including eQTMs, and/or eQTL colocalization, but only 33% of loci involved an eQTL and mQTL present in the same tissue type. With this DNAm-focused integrative analysis, we contribute to the understanding of molecular regulatory mechanisms in human tissues and their impact on complex traits.

Dysregulated Immune and Metabolic Microenvironment Is Associated with the Post-Operative Relapse in Stage I Non-Small Cell Lung Cancer

Simple Summary

The underlying mechanism of post-operative relapse of non-small cell lung cancer (NSCLC) remained poorly understood. This study highlights that both tumors and adjacent tissues from stage I NSCLC with relapse show a dysregulated immune and metabolic environment. Immune response shifts from an active state in primary tumors to a suppressive state in recurrent tumors. A model based on the enriched biological features in the primary tumors with relapse could effectively predict recurrence for stage I NSCLC. These results provide insights into the underpinning of the post-operative relapse and suggest that identifying NSCLC patients with a high risk of relapse could help the clinical decision of applying appropriate therapeutic interventions.

Abstract

The underlying mechanism of post-operative relapse of non-small cell lung cancer (NSCLC) remains poorly understood. We enrolled 57 stage I NSCLC patients with or without relapse and performed whole-exome sequencing (WES) and RNA sequencing (RNA-seq) on available primary and recurrent tumors, as well as on matched tumor-adjacent tissues (TATs). The WES analysis revealed that primary tumors from patients with relapse were enriched with USH2A mutation and 2q31.1 amplification. RNA-seq data showed that the relapse risk was associated with aberrant immune response and metabolism in the microenvironment of primary lesions. TATs from the patients with relapse showed an immunosuppression state. Moreover, recurrent lesions exhibited downregulated immune response compared with their paired primary tumors. Genomic and transcriptomic features were further subjected to build a prediction model classifying patients into groups with different relapse risks. We show that the recurrence risk of stage I NSCLC could be ascribed to the altered immune and metabolic microenvironment. TATs might be affected by cancer cells and facilitate the invasion of tumors. The immune microenvironment in the recurrent lesions is suppressed. Patients with a high risk of relapse need active post-operative intervention.

HOX cluster and their cofactors showed an altered expression pattern in eutopic and ectopic endometriosis tissues

Endometriosis is major gynecological disease that affects over 10% of women worldwide and 30%-50% of these women have pelvic pain, abnormal uterine bleeding and infertility. The cause of endometriosis is unknown and there is no definite cure mainly because of our limited knowledge about its pathophysiology at the cellular and molecular levels. Therefore, demystifying the molecular mechanisms that underlie endometriosis is essential to develop advanced therapies for this disease. In this regard, HOX genes are remarkable because of their critical role in endometrial development and receptivity during implantation, which is attributed to their ability to mediate some of the sex steroid functions during the reproductive period. Access to the expression profiles of these genes would provide the necessary information to uncover new genes for endometriosis and assist with disease diagnosis and treatment. In this study we demonstrate an altered expression pattern for the HOX clusters (A-D) and their cofactors in both eutopic and ectopic conditions compared to control tissue biopsies. Remarkably, most of the intensive changes occurred in eutopic samples from endometriosis patients compared to control tissue biopsies. Pathway analysis revealed the involvement of differentially expressed genes in cancer that correlate with an association between endometriosis and cancer. Our results suggest critical roles for the HOX cluster and their cofactors in endometriosis pathophysiology.

Induction of HOX Genes by Hepatitis C Virus Infection via Impairment of Histone H2A Monoubiquitination

Hepatitis C virus (HCV) infection causes liver pathologies, including hepatocellular carcinoma (HCC). Homeobox (HOX) gene products regulate embryonic development and are associated with tumorigenesis, although the regulation of HOX genes by HCV infection has not been clarified in detail. We examined the effect of HCV infection on HOX gene expression. In this study, HCV infection induced more than half of the HOX genes and reduced the level of histone H2A monoubiquitination on lysine 119 (K119) (H2Aub), which represses HOX gene promoter activity. HCV infection also promoted proteasome-dependent degradation of RNF2, which is an E3 ligase mediating H2A monoubiquitination as a component of polycomb repressive complex 1. Since full-genomic replicon cells but not subgenomic replicon cells exhibited reduced RNF2 and H2Aub levels and induction of HOX genes, we focused on the core protein. Expression of the core protein reduced the amounts of RNF2 and H2Aub and induced HOX genes. Treatment with LY-411575, which can reduce HCV core protein expression via signal peptide peptidase (SPP) inhibition without affecting other viral proteins, dose-dependently restored the amounts of RNF2 and H2Aub in HCV-infected cells and impaired the induction of HOX genes and production of viral particles but not viral replication. The chromatin immunoprecipitation assay results also indicated infection- and proteasome-dependent reductions in H2Aub located in HOX gene promoters. These results suggest that HCV infection or core protein induces HOX genes by impairing histone H2A monoubiquitination via a reduction in the RNF2 level.IMPORTANCE Recently sustained virologic response can be achieved by direct-acting antiviral (DAA) therapy in most hepatitis C patients. Unfortunately, DAA therapy does not completely eliminate a risk of hepatocellular carcinoma (HCC). Several epigenetic factors, including histone modifications, are well known to contribute to hepatitis C virus (HCV)-associated HCC. However, the regulation of histone modifications by HCV infection has not been clarified in detail. In this study, our data suggest that HCV infection or HCV core protein expression impairs monoubiquitination of histone H2A K119 in the homeobox (HOX) gene promoter via destabilization of RNF2 and then induces HOX genes. Several lines of evidence suggest that the expression of several HOX genes is dysregulated in certain types of tumors. These findings reveal a novel mechanism of HCV-related histone modification and may provide information about new targets for diagnosis and prevention of HCC occurrence.

Loss of HOXB3 correlates with the development of hormone receptor negative breast cancer

Background

The homeobox gene family, encoding a specific nuclear protein, is essential for embryonic development, differentiation, and homeostasis. The role of the HOXB3 protein varies in different tumors. This study aims to explore the role of the HOXB3 gene in breast cancer.

Method

Differentially expressed genes were screened by analyzing metastatic breast cancer gene chip data from TCGA and GEO databases. The function of the selected HOXB3 gene was also analyzed in different databases and through molecular biology methods, such as qRT-PCR, western blot and IF to verify bioinformatics findings.

Results

Both bioinformatics analyses and western blot showed that HOXB3 was lost in breast cancer compared to normal breast tissue. Survival analysis also showed that lower expression of HOXB3 was associated with poor prognosis. Bioinformatics analyses further showed that HOXB3 was positively correlated with hormone receptors. Metascape for GO analysis of GEO data provided possible mechanisms that HOXB3 could positively regulate cell adhesion, inhibit cell proliferation and activate immune response in breast cancer; moreover, GSEA included several cancer-associated pathways.

Conclusion

In summary, HOXB3 expression was decreased in breast cancer, and it was associated with poor prognosis. It might become a new biomarker to predict prognosis of breast cancer.

Genome-wide Network-assisted Association and Enrichment Study of Amyloid Imaging Phenotype in Alzheimer's Disease

BACKGROUND

The etiology of Alzheimer's disease remains poorly understood at the mechanistic level, and genome-wide network-based genetics have the potential to provide new insights into the disease mechanisms.

OBJECTIVE

The study aimed to explore the collective effects of multiple genetic association signals on an AV-45 PET measure, which is a well-known Alzheimer's disease biomarker, by employing a network assisted strategy.

METHODS

First, we took advantage of a dense module search algorithm to identify modules enriched by genetic association signals in a protein-protein interaction network. Next, we performed statistical evaluation to the modules identified by dense module search, including a normalization process to adjust the topological bias in the network, a replication test to ensure the modules were not found randomly , and a permutation test to evaluate unbiased associations between the modules and amyloid imaging phenotype. Finally, topological analysis, module similarity tests and functional enrichment analysis were performed for the identified modules.

RESULTS

We identified 24 consensus modules enriched by robust genetic signals in a genome-wide association analysis. The results not only validated several previously reported AD genes (APOE, APP, TOMM40, DDAH1, PARK2, ATP5C1, PVRL2, ELAVL1, ACTN1 and NRF1), but also nominated a few novel genes (ABL1, ABLIM2) that have not been studied in Alzheimer's disease but have shown associations with other neurodegenerative diseases.

CONCLUSION

The identified genes, consensus modules and enriched pathways may provide important clues to future research on the neurobiology of Alzheimer's disease and suggest potential therapeutic targets.

HOXB5 promotes the progression of breast cancer through wnt/beta-catenin pathway

OBJECTIVE

The present study was designed to explore the function of HOXB5 in breast cancer and related signaling pathway.

METHODS

Breast cancer tissues and non-cancerous tissues were collected from 82 cases who were pathologically diagnosed with breast cancer. The mRNA level of HOXB5 was detected via quantitative real-time polymerase chain reaction (qRT-PCR). Chi-square test was adopted to analyze the association of HOXB5 with clinical features. The viability, migration and invasion of breast cancer cells were detected through MTT and Transwell assays, respectively. Protein analysis was performed adopting western blot analysis.

RESULTS

HOXB5 expression was increased in breast cancer tissues and cells, and showed positive correlation with tumor size (P = 0.028), TNM stage (P = 0.048), and lymph node metastasis (P = 0.002). Losing HOXB5 expression suppressed clone formation, proliferation, migration and invasion of breast cancer cells. The knockdown of HOXB5 significantly inactivated wnt/β-catenin pathway. Furthermore, wnt/β-catenin pathway had the potential to neutralize the oncogenic function of HOXB5 in breast cancer.

CONCLUSION

HOXB5 may be involved in the invasive progression of breast cancer. The function of HOXB5 in breast cancer was mediated by wnt/β-catenin pathway.

HOXB5 acts as an oncogenic driver in head and neck squamous cell carcinoma via EGFR/Akt/Wnt/β-catenin signaling axis

INTRODUCTION

Identification of therapeutic targets in head and neck squamous cell carcinoma (HNSCC) is essential because most of the patients with advanced HNSCC have a poor prognosis. Homeobox genes constitute a large cluster of transcription factors with important regulatory roles in mammalian embryonic development and cell differentiation. The oncogenic role of homeobox B5 (HOXB5) in HNSCC has not been investigated.

MATERIALS AND METHODS

We used The Cancer Genome Atlas (TCGA) data to evaluate the correlations between HOXB5 expression and various HNSCC clinicopathological factors. We knocked down HOXB5 expression in HNSCC cell lines and explored the in vitro and in vivo effects on cell proliferation and motility, and HOXB5 signaling.

RESULTS

The Cancer Genome Atlas (TCGA) data shows that HOXB5 is overexpressed in HNSCC compared to normal tissues and significantly associates with tumor stage (P = 0.003), lymph node metastasis (P = 0.031), disease stage (P = 0.002), and angiolymphatic invasion (P = 0.004). Our results also show that HOXB5 expression is up-regulated in HNSCC cell lines, and HOXB5 knockdown significantly reduced cell proliferation and tumor growth in vitro and in vivo. Inhibition of HOXB5 decreases cell migration and invasion via suppression of epithelial-to-mesenchymal transition (EMT)-associated proteins expression. Moreover, HOXB5 directly binds to the promoter region of EGFR and consequently regulates the activity of the Akt/Wnt/β-catenin signaling axis.

CONCLUSION

HOXB5 may be a novel therapeutic target as an oncogenic driver by regulating EGFR transcription in HNSCC.

The prospect of discovering new biomarkers for ovarian cancer based on current knowledge of susceptibility loci and genetic variation (Review)

Ovarian cancer is the most lethal gynaecological malignancy. The cancer initially presents with non‑specific symptoms; thus, it is typically not discovered until the patient has reached the late, considerably more lethal, stages of the disease. Research focus is currently on finding novel biomarkers, especially for early detection and stratification of the disease. One promising approach has been to focus on mutations or variations in the genetic code that are associated with the risk of developing ovarian cancer. A certain heritable component is already known regarding genes such as BRCA1/2, TP53, MSH6, BRIP1 and RAD51C, yet these are estimated to only account for ~3.1% of the total risk. Recent advances in sequencing technologies have enabled the investigation of hundreds of thousands of genetic variants in genome‑wide association studies in tens of thousands of patients, which has led to the discovery of 108 (39 loci with P<5.0x10‑8) novel susceptibility loci for ovarian cancer, presented in this review. Using the published variants in a patient cohort screening, together with variants identified in our ongoing whole exome sequencing project, future aims are to ascertain whether certain of the novel variants could be used as biomarkers for early diagnosis and/or treatment decisions.

HOX Genes in High Grade Ovarian Cancer

HOX genes are highly conserved members of the homeobox superfamily that have a crucial role in determining cellular identity. High grade ovarian cancer is the most lethal gynaecological malignancy. Our understanding of the role of HOX genes in the oncogenesis of ovarian cancer is evolving, and here we review their dysregulated expression patterns, their function in cell survival and invasion, their potential uses as biomarkers, and ways in which HOX genes are being targeted with new and existing drugs.

Genetic Data from Nearly 63,000 Women of European Descent Predicts DNA Methylation Biomarkers and Epithelial Ovarian Cancer Risk

DNA methylation is instrumental for gene regulation. Global changes in the epigenetic landscape have been recognized as a hallmark of cancer. However, the role of DNA methylation in epithelial ovarian cancer (EOC) remains unclear. In this study, high-density genetic and DNA methylation data in white blood cells from the Framingham Heart Study (N = 1,595) were used to build genetic models to predict DNA methylation levels. These prediction models were then applied to the summary statistics of a genome-wide association study (GWAS) of ovarian cancer including 22,406 EOC cases and 40,941 controls to investigate genetically predicted DNA methylation levels in association with EOC risk. Among 62,938 CpG sites investigated, genetically predicted methylation levels at 89 CpG were significantly associated with EOC risk at a Bonferroni-corrected threshold of P < 7.94 × 10-7. Of them, 87 were located at GWAS-identified EOC susceptibility regions and two resided in a genomic region not previously reported to be associated with EOC risk. Integrative analyses of genetic, methylation, and gene expression data identified consistent directions of associations across 12 CpG, five genes, and EOC risk, suggesting that methylation at these 12 CpG may influence EOC risk by regulating expression of these five genes, namely MAPT, HOXB3, ABHD8, ARHGAP27, and SKAP1. We identified novel DNA methylation markers associated with EOC risk and propose that methylation at multiple CpG may affect EOC risk via regulation of gene expression. SIGNIFICANCE: Identification of novel DNA methylation markers associated with EOC risk suggests that methylation at multiple CpG may affect EOC risk through regulation of gene expression.

Functional Analysis and Fine Mapping of the 9p22.2 Ovarian Cancer Susceptibility Locus

Genome-wide association studies have identified 40 ovarian cancer risk loci. However, the mechanisms underlying these associations remain elusive. In this study, we conducted a two-pronged approach to identify candidate causal SNPs and assess underlying biological mechanisms at chromosome 9p22.2, the first and most statistically significant associated locus for ovarian cancer susceptibility. Three transcriptional regulatory elements with allele-specific effects and a scaffold/matrix attachment region were characterized and, through physical DNA interactions, BNC2 was established as the most likely target gene. We determined the consensus binding sequence for BNC2 in vitro, verified its enrichment in BNC2 ChIP-seq regions, and validated a set of its downstream target genes. Fine-mapping by dense regional genotyping in over 15,000 ovarian cancer cases and 30,000 controls identified SNPs in the scaffold/matrix attachment region as among the most likely causal variants. This study reveals a comprehensive regulatory landscape at 9p22.2 and proposes a likely mechanism of susceptibility to ovarian cancer. SIGNIFICANCE: Mapping the 9p22.2 ovarian cancer risk locus identifies BNC2 as an ovarian cancer risk gene.See related commentary by Choi and Brown, p. 439.

A Transcriptome-Wide Association Study Among 97,898 Women to Identify Candidate Susceptibility Genes for Epithelial Ovarian Cancer Risk

Large-scale genome-wide association studies (GWAS) have identified approximately 35 loci associated with epithelial ovarian cancer (EOC) risk. The majority of GWAS-identified disease susceptibility variants are located in noncoding regions, and causal genes underlying these associations remain largely unknown. Here, we performed a transcriptome-wide association study to search for novel genetic loci and plausible causal genes at known GWAS loci. We used RNA sequencing data (68 normal ovarian tissue samples from 68 individuals and 6,124 cross-tissue samples from 369 individuals) and high-density genotyping data from European descendants of the Genotype-Tissue Expression (GTEx V6) project to build ovarian and cross-tissue models of genetically regulated expression using elastic net methods. We evaluated 17,121 genes for their cis-predicted gene expression in relation to EOC risk using summary statistics data from GWAS of 97,898 women, including 29,396 EOC cases. With a Bonferroni-corrected significance level of P < 2.2 × 10-6, we identified 35 genes, including FZD4 at 11q14.2 (Z = 5.08, P = 3.83 × 10-7, the cross-tissue model; 1 Mb away from any GWAS-identified EOC risk variant), a potential novel locus for EOC risk. All other 34 significantly associated genes were located within 1 Mb of known GWAS-identified loci, including 23 genes at 6 loci not previously linked to EOC risk. Upon conditioning on nearby known EOC GWAS-identified variants, the associations for 31 genes disappeared and three genes remained (P < 1.47 × 10-3). These data identify one novel locus (FZD4) and 34 genes at 13 known EOC risk loci associated with EOC risk, providing new insights into EOC carcinogenesis.Significance: Transcriptomic analysis of a large cohort confirms earlier GWAS loci and reveals FZD4 as a novel locus associated with EOC risk. Cancer Res; 78(18); 5419-30. ©2018 AACR.

HOXA4/HOXB3 gene expression signature as a biomarker of recurrence in patients with high-grade serous ovarian cancer following primary cytoreductive surgery and first-line adjuvant chemotherapy

OBJECTIVES

Aberrant homeobox (HOX) gene expression is reported in high-grade serous ovarian carcinoma (HGSOC), however, its prognostic significance remains unclear.

METHODS

HOX genes associated with progression-free survival (PFS) in a discovery cohort of primary HGSOC samples with RNA sequencing data, and those previously reported to be associated with clinical outcomes, were selected for qPCR testing in an independent training cohort of primary HGSOC samples (n=71). A prognostic model for PFS was developed using univariate and multivariate Cox regression. Patients were stratified into risk groups that optimized the test statistic. The model was tested in an independent HGSOC cohort from The Cancer Genome Atlas (TCGA) (n=320). The effect of selected HOX genes on drug sensitivity and reactive oxygen species (ROS) accumulation was examined in vitro.

RESULTS

Of 23 HOX genes tested in the training cohort, HOXA4 (HR=1.20, 95% CI=1.07-1.34, P=0.002) and HOXB3 (HR=1.09, 95% CI=1.01-1.17, P=0.027) overexpression were significantly associated with shorter PFS in multivariate analysis. Based on the optimal cutoff of the HOXA4/HOXB3 risk score, median PFS was 16.9months (95% CI=14.6-21.2months) and not reached (>80months) for patients with high and low risk scores, respectively (HR=8.89, 95% CI=2.09-37.74, P<0.001). In TCGA, the HOXA4/HOXB3 risk score was significantly associated with disease-free survival (HR=1.44, 95% CI=1.00-2.09, P=0.048). HOXA4 or HOXB3 overexpression in ovarian cancer cells decreased sensitivity to cisplatin and attenuated the generation of cisplatin-induced ROS (P<0.05).

CONCLUSIONS

HOXA4/HOXB3 gene expression-based risk score may be useful for prognostic risk stratification and warrants prospective validation in HGSOC patients.

HOXB5 promotes retinoblastoma cell migration and invasion via ERK1/2 pathway-mediated MMPs production

Homeobox genes (HOX genes) have been implicated in many tumors. As a member of HOX genes, HOXB5 is overexpressed in bladder cancer and contributes to the proliferation and invasion of breast cancer cells. However, functions of HOXB5 in retinoblastoma remain elusive. In this study, we found that HOXB5 expression is upregulated in retinoblastoma cell lines and tissues. Overexpression of HOXB5 promoted retinoblastoma cell migration and invasion, but knockdown of HOXB5 suppressed the migration and invasion. Moreover, HOXB5 induced the activation of ERK1/2 and upregulated the production of MMP-3 and MMP-13. In addition, ERK1/2 pathway was required for HOXB5-mediated retinoblastoma cell migration and invasion. Taken together, our study suggests that HOXB5 promotes the migration and invasion of retinoblastoma cells by inducing the activation of ERK1/2 and increasing the production of MMP-3 and MMP-13. Therefore, HOXB5 may represent an effective target for treatment of retinoblastoma.

Common Genetic Variation and Susceptibility to Ovarian Cancer: Current Insights and Future Directions

In this review, we summarize current progress in the genetic epidemiology of epithelial ovarian cancer (EOC), focusing exclusively on elucidating the role of common germline genetic variation in conferring susceptibility to EOC. We provide an overview of the more than 30 EOC risk loci identified to date by genome-wide association studies (GWAS) and describe the contribution of large-scale, cross-cancer type, custom genotyping projects, such as the OncoArray and the Collaborative Oncological Gene-Environment Study, to locus discovery and replication. We discuss the histotype-specific nature of these EOC risk loci, pleiotropy, or overlapping genetic effects between EOC and other hormone-related cancer types, and the application of findings to polygenic risk prediction for EOC. The second part of the article offers a concise review of primarily laboratory-based studies that have led to the identification of several putative EOC susceptibility genes using common variants at the known EOC risk loci as starting points. More global biological insights emerging from network- and pathway-based analyses of GWAS for EOC susceptibility are also highlighted. Finally, we delve into potential future directions, including the need to identify EOC risk loci in non-European populations and the next generation of GWAS functional studies that are likely to involve genome editing to establish the cell type-specific carcinogenic effects of EOC risk variants Cancer Epidemiol Biomarkers Prev; 27(4); 395-404. ©2018 AACRSee all articles in this CEBP Focus section, "Genome-Wide Association Studies in Cancer."

Dysregulated expression of homeobox family genes may influence survival outcomes of patients with epithelial ovarian cancer: analysis of data from The Cancer Genome Atlas

Homeobox (HOX) family genes encode key transcription factors for embryogenesis and may be correlated with carcinogenesis. The aim of this study was to elucidate whether aberrant expression of HOX genes influences outcomes in epithelial ovarian cancer (EOC). Gene expression data and clinicopathologic information from 630 patients with EOC were downloaded from The Cancer Genome Atlas database. We explored correlations between expression levels of HOX gene family members and clinicopathological variables. Higher expression of HOXA1, A4, A5, A7, A10, A11, B13, C13, D1, and D3 was associated with advanced FIGO stage. Suboptimal residual disease after debulking surgery was significantly correlated with higher expression of HOXB9, B13, and C13. Additionally, patients with high expression of HOXC6 and C11 were significantly more likely to have poor performance status. Overall survival was significantly shorter in patients with high, rather than low, expression of two HOX genes (HOXA10 and B3), and significantly longer in patients with high rather than low HOXC5 expression. Dysregulated expression of the HOXA10, B3, and C5 was significantly correlated with overall survival in EOC patients. HOX gene expression levels are potentially useful as a prognostic indicator in EOC, and HOX genes may represent a novel and promising target for anticancer therapeutics.

Enrichment of putative PAX8 target genes at serous epithelial ovarian cancer susceptibility loci

BACKGROUND

Genome-wide association studies (GWAS) have identified 18 loci associated with serous ovarian cancer (SOC) susceptibility but the biological mechanisms driving these findings remain poorly characterised. Germline cancer risk loci may be enriched for target genes of transcription factors (TFs) critical to somatic tumorigenesis.

METHODS

All 615 TF-target sets from the Molecular Signatures Database were evaluated using gene set enrichment analysis (GSEA) and three GWAS for SOC risk: discovery (2196 cases/4396 controls), replication (7035 cases/21 693 controls; independent from discovery), and combined (9627 cases/30 845 controls; including additional individuals).

RESULTS

The PAX8-target gene set was ranked 1/615 in the discovery (PGSEA<0.001; FDR=0.21), 7/615 in the replication (PGSEA=0.004; FDR=0.37), and 1/615 in the combined (PGSEA<0.001; FDR=0.21) studies. Adding other genes reported to interact with PAX8 in the literature to the PAX8-target set and applying an alternative to GSEA, interval enrichment, further confirmed this association (P=0.006). Fifteen of the 157 genes from this expanded PAX8 pathway were near eight loci associated with SOC risk at P<10-5 (including six with P<5 × 10-8). The pathway was also associated with differential gene expression after shRNA-mediated silencing of PAX8 in HeyA8 (PGSEA=0.025) and IGROV1 (PGSEA=0.004) SOC cells and several PAX8 targets near SOC risk loci demonstrated in vitro transcriptomic perturbation.

CONCLUSIONS

Putative PAX8 target genes are enriched for common SOC risk variants. This finding from our agnostic evaluation is of particular interest given that PAX8 is well-established as a specific marker for the cell of origin of SOC.

HOXC10 up-regulation promotes gastric cancer cell proliferation and metastasis through MAPK pathway

Objective

As an important regulator of embryonic morphogenesis, homeodomain-containing gene 10 (HOXC10) has been found to promote progression of human cancers and its expression indicates poor survival outcome. However, very few studies are available on the role of HOXC10 in gastric carcinoma. Therefore, the aim of this study was to determine the role of HOXC10 in gastric cancer and the potential mechanism underlying its function for cancer biology.

Methods

A primary gastric cancer mouse model was obtained via intra-gastric wall injection of gastric cancer cells and was used to evaluate the function of HOXC10 during gastric cancer progression in vivo. Immunohistochemistry was performed to visualize and measure HOXC10 protein expression in gastric cancer tissue. Cells were transfected with plasmids to increase the expression of HOXC10, and siRNA transfection was performed to suppress HOXC10 expression. Reverse transcription polymerase chain reaction (RT-PCR) and western blotting were utilized to measure mRNA and protein expression, respectively. Proliferation, migration, and invasion were investigated using 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT) assay, wound healing assay, and matrigel invasion assay in vitro, respectively.

Results

HOXC10 expression was significantly increased in gastric cancer tissues compared to matched normal tissues. HOXC10 up-regulation significantly increased tumor volumes in nude mice. Plasmid transfection significantly increased HOXC10 protein and mRNA expressions and effectively promoted cell proliferation. Moreover, HOXC10 up-regulation significantly promoted migration and invasion of gastric cancer cells. Mechanistic investigation showed that HOXC10 up-regulation significantly increased mRNA and protein expression of mitogen-activated protein kinase (MAPK) signaling related genes, including c-myc, c-jun and p53, while also modulating the phosphorylation of c-Jun N-terminal kinase (JNK), extracellular signal-regulated kinase (ERK) and P38 but not their total protein levels.

Conclusions

This study demonstrated the tight link between HOXC10 and gastric cancer cell proliferation and metastasis via involvement of the MAPK pathway.

Dysregulated expression of homebox gene HOXA13 is correlated with the poor prognosis in bladder cancer

BACKGROUND

The homeobox (HOX) genes have been implicated playing important roles in many human cancers. HOXA13 is a member of HOX genes that encode transcription factors regulating embryonic development and cell fate. In the present study, we aimed to investigate the expression and prognostic significance of HOXA13 in bladder cancer.

METHODS

Immunohistochemical staining was initially performed to screen the differentially expressed HOXA13 between bladder cancer tissues and paired adjacent non-cancerous tissues. Subsequent Western blotting analysis validation was conducted using tissue samples from patients with bladder cancer.

RESULTS

The expression level of HOXA13 was significantly higher in bladder cancer tissues compared to that in adjacent non-cancerous tissues (P < 0.001). The χ²-test showed that expression of HOXA13 was positively correlated with lymphatic metastasis (P = 0.013), bladder tumor TNM stage (P = 0.002) and pathological grade (P < 0.001). Kaplan-Meier survival analysis revealed that patients with bladder cancer with high expression of HOXA13 had shorter overall survival time (P = 0.001) and disease-free survival time (P = 0.001) compared to patients with low expression of HOXA13. Multivariate analysis indicated that HOXA13 was an independent prognostic factor for overall survival for bladder cancer patients.

CONCLUSIONS

The results of our study show that high expression of HOXA13 is associated with the progression of bladder cancer and that HOXA13 might serve as a biomarker for prognosis of bladder cancer.

Functional mechanisms underlying pleiotropic risk alleles at the 19p13.1 breast-ovarian cancer susceptibility locus

A locus at 19p13 is associated with breast cancer (BC) and ovarian cancer (OC) risk. Here we analyse 438 SNPs in this region in 46,451 BC and 15,438 OC cases, 15,252 BRCA1 mutation carriers and 73,444 controls and identify 13 candidate causal SNPs associated with serous OC (P=9.2 × 10(-20)), ER-negative BC (P=1.1 × 10(-13)), BRCA1-associated BC (P=7.7 × 10(-16)) and triple negative BC (P-diff=2 × 10(-5)). Genotype-gene expression associations are identified for candidate target genes ANKLE1 (P=2 × 10(-3)) and ABHD8 (P<2 × 10(-3)). Chromosome conformation capture identifies interactions between four candidate SNPs and ABHD8, and luciferase assays indicate six risk alleles increased transactivation of the ADHD8 promoter. Targeted deletion of a region containing risk SNP rs56069439 in a putative enhancer induces ANKLE1 downregulation; and mRNA stability assays indicate functional effects for an ANKLE1 3'-UTR SNP. Altogether, these data suggest that multiple SNPs at 19p13 regulate ABHD8 and perhaps ANKLE1 expression, and indicate common mechanisms underlying breast and ovarian cancer risk.

Genome-wide association studies and epigenome-wide association studies go together in cancer control

Completion of the human genome a decade ago laid the foundation for: using genetic information in assessing risk to identify individuals and populations that are likely to develop cancer, and designing treatments based on a person's genetic profiling (precision medicine). Genome-wide association studies (GWAS) completed during the past few years have identified risk-associated single nucleotide polymorphisms that can be used as screening tools in epidemiologic studies of a variety of tumor types. This led to the conduct of epigenome-wide association studies (EWAS). This article discusses the current status, challenges and research opportunities in GWAS and EWAS. Information gained from GWAS and EWAS has potential applications in cancer control and treatment.

Expression profiles of long noncoding RNAs in cutaneous squamous cell carcinoma

BACKGROUND

Despite there being over 35,000 different long noncoding RNA (lncRNA) sequences described little is known regarding their molecular-pathological role in cutaneous squamous cell carcinoma (cSCC).

MATERIALS & METHODS

In this pilot study, lncRNA and mRNA expression profiles were determined in cSCC and control (n = 6) by an Arraystar human lncRNA Microarray. Kyoto Encyclopedia of Genes and Genomes pathway enrichment and gene ontology analysis of mRNAs was performed.

RESULTS

Analysis of differential expression revealed 1516 upregulated lncRNAs and 2586 downregulated lncRNAs in cSCC compared with controls. Data analysis identified known oncogenic lncRNAs, such as the HOX transcript antisense RNA HOTAIR, among the differentially expressed lncRNA sequences. Kyoto Encyclopedia of Genes and Genomes pathway analysis showed that focal adhesion, extracellular matrix and the oncogenic phosphatidylinositol 3'-kinase-Akt signaling pathway had the highest enrichment scores.

CONCLUSION

This study provides the first evidence for differential expression of lncRNA in cSCC and serves as a template for further, larger functional in-depth analyses regarding cSCC molecular lncRNAs.