Genomic aberration patterns and expression profiles of squamous cell carcinomas of the vulva

DNA methylation and copy number alterations in the progression of HPV-associated high-grade vulvar intraepithelial lesion

Human papillomavirus (HPV)-associated high-grade vulvar intraepithelial lesion (HSIL) is a precursor of vulvar squamous cell carcinoma (VSCC). Because of the 8% cancer risk, many vulvar HSIL patients undergo aggressive and mutilating treatments. Characterizing HSIL by their progression risk can help individualize treatment strategies. Accordingly, copy number alterations (CNAs) and DNA methylation have been identified as biomarkers for cancer risk stratification of HSIL. Here, we assessed their potential correlation, and relation to HPV16 (sub)lineages and progression to vulvar cancer. Eighty-two vulvar formalin-fixed paraffin-embedded (FFPE) samples, including controls, HSIL, HSIL adjacent to VSCC and VSCC, with previously determined DNA methylation profiles, were analysed for CNAs using mFAST-SeqS. Genome-wide z-scores were calculated to determine overall aneuploidy (aneuploidy scores), and compared to the methylation levels and status of marker panel ZNF582/SST/miR124-2. For 52 HPV16-positive cases, HPV (sub)lineages were determined by Sanger sequencing. HPV16 lineage A was predominant (86.4%), followed equally by lineages B, C, and D. Frequent chromosomal alterations included chr1pq, chr3q, chr9q gains, and chr2q, chr4q losses. Median aneuploidy scores increased across disease categories, from 0 in controls, to 3 in HSIL, 16 in HSIL adjacent to VSCC and 29 in VSCC. A positive relationship between aneuploidy scores and DNA methylation levels was found (ρ = 0.61, Spearman's rank correlation test). Aneuploidy scores were significantly higher in methylation-positive samples (p < .001). In conclusion, we showed that DNA methylation and CNAs both rise with increasing severity of disease, indicating their prognostic value for cancer risk stratification of HSIL, while no relation to HPV16 (sub)lineages was found.

Expression Profiles Reveal Involvement of VEGF, IGF1, BIRC5, and MMP1 in Vulvar Carcinogenesis

OBJECTIVE

The objective of this study was to identify key genes and shed light on the underlying molecular mechanisms of vulvar squamous cell carcinoma (VSCC).

METHODS

Bioinformatic software was utilized for the identification and characterization of key differentially expressed genes (DEGs) from microarrays GSE63678 and GSE38228, which contain VSCC and normal vulvar tissue data. These microarrays were obtained from Gene Expression Omnibus (GEO). Immunohistochemical assays (55 VSCC and 50 normal vulvar tissues) were utilized to validate the expression of VEGF, IGF1, BIRC5, and MMP1 screened from the identified DEGs. SPSS 18.0 software was used for statistical analyses of the relationships between IGF1, BIRC5, VEGF, MMP1 expression levels and patient clinicopathological characteristics.

RESULTS

A total of 141 DEGs were identified, among which 18 genes were closely correlated with the biological characteristics of VSCC. Four of the 18 genes (VEGF, IGF1, BIRC5, and MMP1) screened from the GEO database were markedly enriched in pathways in cancer (P < 0.05), and could be considered key genes in VSCC based on KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway analysis in DAVID (Database for Annotation, Visualization and Integrated Discovery).The expression levels of these 4 hub genes, determined by immunohistochemical assays, were consistent with the bioinformatics results. Higher expression of IGF1 showed significant association with well-differentiated carcinomas (P = 0.017).BIRC5 expression levels showed a positive correlation with clinical stage (P = 0.039); compared with those in menopause for over 10 years, patients in menopause for less than 10 years at the time of diagnosis tended to have significantly higher expression of BIRC5 (P = 0.003). VEGF and MMP1 expression levels were not correlated with any of the tested clinicopathological characteristics.

CONCLUSION

VEGF, IGF1, BIRC5, and MMP1 were identified as being associated with VSCC using integrated bioinformatic methods, which may provide important insights into the pathogenesis of this disease and help to identify new biomarkers.

Nuclear factor IB is downregulated in vulvar squamous cell carcinoma (VSCC): Unravelling differentially expressed genes in VSCC through gene expression dataset analysis

Vulvar squamous cell carcinoma (VSCC) comprises two distinct etiopathological subtypes: i) Human papilloma virus (HPV)-related VSCC, which arises via the precursor high grade squamous intraepithelial lesion (HSIL); and ii) HPV-independent VSCC, which arises via precursor, differentiated vulvar intraepithelial neoplasia (dVIN), driven by TP53 mutations. However, the mechanism of carcinogenesis of VSCC is poorly understood. The current study aimed to gain insight into VSCC carcinogenesis by identifying differentially expressed genes (DEGs) for each VSCC subtype. The expression of certain DEGs was then further assessed by performing immunohistochemistry (IHC) on whole tissue sections of VSCC and its precursors. Statistical analysis of microarrays was performed on two independent gene expression datasets (GSE38228 and a study from Erasmus MC) on VSCC and normal vulva. DEGs were identified that were similarly (up/down) regulated with statistical significance in both datasets. For HPV-related VSCCs, this constituted 88 DEGs, and for HPV-independent VSCCs, this comprised 46 DEGs. IHC was performed on VSCC (n=11), dVIN (n=6), HSIL (n=6) and normal vulvar tissue (n=7) with i) signal transducer and activator of transcription 1 (STAT1; an upregulated DEGs); ii) nuclear factor IB (NFIB; a downregulated DEG); iii) p16 (to determine the HPV status of tissues); and iv) p53 (to confirm the histological diagnoses). Strong and diffuse NFIB expression was observed in the basal and para-basal layers of normal vulvar tissue, whereas NFIB expression was minimal or completely negative in dVIN and in both subtypes of VSCC. In contrast, no discernable difference was observed in STAT1 expression among normal vulvar tissue, dVIN, HSIL or VSCC. By leveraging bioinformatics, the current study identified DEGs that can facilitate research into VSCC carcinogenesis. The results suggested that NFIB is downregulated in VSCC and its relevance as a diagnostic/prognostic biomarker deserves further exploration.

EZH2 facilitates BMI1-dependent hepatocarcinogenesis through epigenetically silencing microRNA-200c

EZH2, a histone methyltransferase, has been shown to involve in cancer development and progression via epigenetic regulation of tumor suppressor microRNAs, whereas BMI1, a driver of hepatocellular carcinoma (HCC), is a downstream target of these microRNAs. However, it remains unclear whether EZH2 can epigenetically regulate microRNA expression to modulate BMI1-dependent hepatocarcinogenesis. Here, we established that high EZH2 expression correlated with enhanced tumor size, elevated metastasis, increased relapse, and poor prognosis in HCC patients. Further clinical studies revealed that EZH2 overexpression was positively correlated to its gene copy number gain/amplification in HCC. Mechanistically, EZH2 epigenetically suppressed miR-200c expression both in vitro and in vivo, and more importantly, miR-200c post-transcriptionally regulated BMI1 expression by binding to the 3'-UTR region of its mRNA. Furthermore, miR-200c overexpression inhibits the growth of HCC cells in vivo. Silencing miR-200c rescued the tumorigenicity of EZH2-depleted HCC cells, whereas knocking down BMI1 reduced the promoting effect of miR-200c depletion on HCC cell migration. Finally, combination treatment of EZH2 and BMI1 inhibitors further inhibited the viability of HCC cells compared with the cells treated with EZH2 or BMI1 inhibitor alone. Our findings demonstrated that alteration of EZH2 gene copy number status induced BMI1-mediated hepatocarcinogenesis via epigenetically silencing miR-200c, providing novel therapeutic targets for HCC treatment.

Precursor lesions of vulvar squamous cell carcinoma - histology and biomarkers: A systematic review

The precursor lesion of vulvar squamous cell carcinoma (VSCC), namely vulvar intraepithelial neoplasia (VIN), is classified as: human papillomavirus (HPV)-related high grade squamous intraepithelial lesion (HSIL), and HPV-independent differentiated VIN (dVIN). Traditionally, histology and immunohistochemistry (IHC) have been the basis of diagnosis and classification of VIN. HSIL shows conspicuous histological atypia, and positivity on p16-IHC, whereas dVIN shows less obvious histological atypia, and overexpression or null-pattern on p53-IHC. For both types of VIN, other diagnostic immunohistochemical markers have also been evaluated. Molecular characterization of VIN has been attempted in few recent studies, and novel genotypic subtypes of HPV-independent VSCC and VIN have been identified. This systematic review appraises the VSCC precursors identified so far, focusing on histology and biomarkers (immunohistochemical and molecular). To gain further insights into the carcinogenesis and to identify additional potential biomarkers, gene expression omnibus (GEO) datasets on VSCC were analyzed; the results are presented.

A meta-analysis of multiple matched aCGH/expression cancer datasets reveals regulatory relationships and pathway enrichment of potential oncogenes

The copy numbers of genes in cancer samples are often highly disrupted and form a natural amplification/deletion experiment encompassing multiple genes. Matched array comparative genomics and transcriptomics datasets from such samples can be used to predict inter-chromosomal gene regulatory relationships. Previously we published the database METAMATCHED, comprising the results from such an analysis of a large number of publically available cancer datasets. Here we investigate genes in the database which are unusual in that their copy number exhibits consistent heterogeneous disruption in a high proportion of the cancer datasets. We assess the potential relevance of these genes to the pathology of the cancer samples, in light of their predicted regulatory relationships and enriched biological pathways. A network-based method was used to identify enriched pathways from the genes’ inferred targets. The analysis predicts both known and new regulator-target interactions and pathway memberships. We examine examples in detail, in particular the gene POGZ, which is disrupted in many of the cancer datasets and has an unusually large number of predicted targets, from which the network analysis predicts membership of cancer related pathways. The results suggest close involvement in known cancer pathways of genes exhibiting consistent heterogeneous copy number disruption. Further experimental work would clarify their relevance to tumor biology. The results of the analysis presented in the database METAMATCHED, and included here as an R archive file, constitute a large number of predicted regulatory relationships and pathway memberships which we anticipate will be useful in informing such experiments.

Olfactomedin 4 downregulation is associated with tumor initiation, growth and progression in human prostate cancer

The olfactomedin 4 (OLFM4) gene has been analyzed as a tumor‐suppressor gene and a putative biomarker in many cancers. In our study, we analyzed the relationship of OLFM4 expression with clinicopathological features and with CpG site methylation in the OLFM4 gene promoter region in human primary prostate adenocarcinoma. OLFM4 protein expression was significantly reduced in prostate cancer tissue compared to adjacent normal tissue and was further significantly reduced in more advanced cancers. Bioinformatic studies with clinical datasets revealed that primary prostate adenocarcinoma patients with reduced OLFM4 mRNA expression exhibited higher Gleason scores and higher preoperative serum prostate‐specific antigen levels, as well as lower recurrence‐free survival. Three of the eight CpG sites in the OLFM4 gene promoter region were hypermethylated in cancerous prostate cells compared to adjacent normal cells, and reduced methylation of eight CpG sites was associated with increased OLFM4 mRNA expression in RWPE1 and PC‐3 cells. Furthermore, knockdown of OLFM4 gene expression was associated with enhanced epithelial–mesenchymal transition (EMT)‐marker expression in RWPE immortalized normal prostate cells. In contrast, restoration of OLFM4 expression in PC‐3 and DU145 prostate cancer cells lacking OLFM4 significantly inhibited both EMT‐marker expression and tumor cell growth in in vitro and in vivo models, indicating that OLFM4 may play a tumor‐suppressor role in inhibiting the EMT program, as well as tumor initiation and growth, in prostate cells. Taken together, these findings suggest that OLFM4 plays an important tumor‐suppressor role in prostate cancer progression and might be useful as a novel candidate biomarker for prostate cancer.

What's new?

Altered expression of the OLFM4 gene appears to be involved in many cancers. In this study of prostate cancers, the authors found that OLFM4 can suppress tumor initiation, growth and progression. Downregulation of OLFM4 was associated with higher serum PSA levels, higher Gleason scores, and lower recurrence‐free survival in prostate cancer patients. These results indicate that OLFM4 may play an important tumor‐suppressor role in the progression of prostate cancer, and may provide a novel prognostic biomarker for prostate cancer treatment.

Molecular heterogeneity in human papillomavirus-dependent and -independent vulvar carcinogenesis

Vulvar squamous cell carcinoma (VSCC) and precancerous vulvar intraepithelial neoplasia (VIN) can develop through human papillomavirus (HPV)-dependent and -independent pathways, indicating a heterogeneous disease. Only a minority of VIN progress, but current clinicopathological classifications are insufficient to predict the cancer risk. Here we analyzed copy number alterations (CNA) to assess the molecular heterogeneity of vulvar lesions in relation to HPV and cancer risk. HPV-status and CNA by means of whole-genome next-generation shallow-sequencing were assessed in VSCC and VIN. The latter included VIN of women with associated VSCC (VINVSCC ) and women who did not develop VSCC during follow-up (VINnoVSCC ). HPV-testing resulted in 41 HPV-positive (16 VINVSCC , 14 VINnoVSCC , and 11 VSCC) and 24 HPV-negative (11 VINVSCC and 13 VSCC) lesions. HPV-positive and -negative VSCC showed a partially overlapping pattern of recurrent CNA, including frequent gains of 3q and 8q. In contrast, amplification of 11q13/cyclinD1 was exclusively found in HPV-negative lesions. HPV-negative VINVSCC had less CNA than HPV-negative VSCC (P = .009), but shared chromosome 8 alterations. HPV-positive VINnoVSCC had less CNA than VINVSCC (P = .022). Interestingly, 1pq gain was detected in 81% of HPV-positive VINVSCC and only in 21% of VINnoVSCC (P = .001). In conclusion, HPV-dependent and -independent vulvar carcinogenesis is characterized by distinct CNA patterns at the VIN stage, while more comparable patterns are present at the cancer stage. Cancer risk in VIN seems to be reflected by the extent of CNA, in particular chromosome 1 gain in HPV-positive cases.

Expressions of miR-30c and let-7a are inversely correlated with HMGA2 expression in squamous cell carcinoma of the vulva

Malignant tumors of the vulva, most of them squamous cell carcinomas, account for only 5% of cancers of the female genital tract. Though little is known about the genetic features of these tumors, the Fragile Histidine Triad (FHIT) and High Mobility Group AT-hook 2 (HMGA2) genes were found deregulated. We wanted to gain more knowledge about the expression of HMGA2-related miRNAs such as miR-30c and let-7a, and whether a correlation exists between the expression of FHIT and HMGA2, in this tumor type. An inverse correlation was found in-as-much as HMGA2 was highly expressed (mean fold change 8.8) whereas miR30c and let-7a were both downregulated (mean fold change -3.9 and -2.3, respectively). The consistent overexpression of HMGA2 found in all tumors adds to the likelihood that this gene is of importance in SCC pathogenesis. Moreover, we came to the conclusion that miRNAs may be the cause of the deregulation of HMGA2. Our results also show that SCC of the vulva presents a characteristic molecular pattern with FHIT being downregulated whereas HMGA2 is upregulated.

Anti-oncostatin M antibody inhibits the pro-malignant effects of oncostatin M receptor overexpression in squamous cell carcinoma

The oncostatin M (OSM) receptor (OSMR) shows frequent gene copy number gains and overexpression in cervical squamous cell carcinomas (SCCs), associated with adverse clinical outcomes. In SCC cells that overexpress OSMR, the major ligand OSM induces multiple pro-malignant effects, including invasion, secretion of angiogenic factors, and metastasis. Here, we demonstrate, for the first time, that OSMR overexpression in SCC cells activates cell-autonomous feed-forward signalling, via further expression of OSMR and OSM and sustained STAT3 activation, despite expression of the negative regulator suppressor of cytokine signalling 3 (SOCS3). The pro-malignant effects associated with OSMR overexpression are critically mediated by JAK-STAT3 activation, which is induced by exogenous OSM and also by autocrine OSM-OSMR interactions. Importantly, specific inhibition of OSM-OSMR interactions by neutralizing antibodies significantly inhibits STAT3 activation and feed-forward signalling, leading to reduced invasion, angiogenesis, and metastasis. Our findings are supported by data from 1254 clinical SCC samples, in which OSMR levels correlated with multiple cognate genes, including OSM, STAT3, and downstream targets. These data strongly support the development of OSM-OSMR-blocking antibodies as biologically targeted therapies against SCCs of the cervix and other anatomical sites. Copyright © 2017 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

Recurrent fusion transcripts in squamous cell carcinomas of the vulva

Juxtaposition of two different genes or gene parts due to chromosomal rearrangement is a well-known neoplasia-associated pathogenetic mechanism. The detection and characterization of such tumorigenic fusions is of great importance both research-wise, diagnostically because they may be specific for distinct tumor entities, and because they may serve as therapeutic targets for antioncogenic drugs that interact directly with the molecular changes responsible for neoplastic transformation.

At present, more than 10,000 fusion transcripts have been reported in different types of neoplasia, with one tenth of them being identified in squamous cell carcinomas (SCC) of different locations. No recurrent fusion gene has to date been identified in SCC of the vulva.

We performed high-throughput paired-end RNA-sequencing of 12 vulvar SCC and found two recurrent fusions with the STIP1-CREB3L1 and ZDHHC5-GPR137 being present in two tumors each. The transcripts were detected only in the tumor samples, not in normal vulvar tissue from healthy donors used as control. The CREB3L1 and ZDHHC5 genes encode proteins involved in transcription suggesting that the chimeras may alter downstream events in their respective pathways. Expression analysis of the CREB3L1 gene showed the presence of two distinct groups of tumors, one having fusion and downregulation of the gene and the other showing upregulation of CREB3L1.

Lower DSC1 expression is related to the poor differentiation and prognosis of head and neck squamous cell carcinoma (HNSCC)

PURPOSE

Although desmocollins have an important position in cancer-related research, there are little reports about the relations between cancers and desmocollin 1 (DSC1). The present study was designed to investigate the correlations between DSC1 and head and neck squamous cell carcinoma (HNSCC).

METHODS

First we analyzed the GEO database; then, HNSCC and pericarcinous tissues were collected to verify the results. DSC1 expression was detected by western blot and real-time PCR. The co-expression genes of DSC1 were extracted from Cancer Cell Line Encyclopedia database (CCLE database), and their correlation was analyzed in The Cancer Genome Atlas HNSCC database (TCGA HNSCC database). Next the gene ontology analysis (GO) was carried out. Moreover, we suppressed DSC1 in FaDu cell to investigate the internal mechanism.

RESULTS

GEO database showed that DSC1 was higher in HNSCC and patients with higher DSC1 had unfavorable prognosis. The results of the samples showed that DSC1 was significantly higher in HNSCC than in normal tissue, which was consistent with the results of GEO database. The co-expression genes of DSC1 were extracted from CCLE database and verified in TCGA HNSCC database. It revealed that DSC1 was related to cell signal transduction. In FaDu/siDSC1 cells, the proliferation and migration were decreased compared to FaDu cells, and the expression levels of β-catenin, c-myc and cyclin D1 down-regulated significantly.

CONCLUSIONS

The increased expression of DSC1 can promote the occurrence of HNSCC and is associated with tumor. The increased expression of DSC1 also indicates a poor prognosis of the patients with HNSCC.

Cytogenetic and molecular profile of endometrial stromal sarcoma

Recent cytogenetic and molecular investigations have improved our understanding of endometrial stromal tumors, including sarcomas (ESS), and helped redefine their classification into more pathogenetically meaningful categories. Because much more can be gained through such studies, we add information on another 22 ESS examined by karyotyping, PCR analysis, expression array analysis, and transcriptome sequencing. In spite of the known preference for certain pathogenetic pathways, we found considerable genetic heterogeneity in high‐grade (HG) as well as in low‐grade (LG) ESS. Not all HG tumors showed a YWHAE‐NUTM chimeric transcript and as many as six LGESS showed no hitherto known ESS‐related fusions. Among the transcripts identified by transcriptome sequencing and verified by Sanger sequencing, new variants of ZC3H7‐BCOR and its reciprocal BCOR‐ZC3H7 were identified as was involvement of the CREBBP and MLLT4 genes (both well known leukemia‐related genes) in two new fusions. FISH analysis identified a known EPC1‐PHF1 fusion which led to the identification of a new variant at the molecular level. The fact that around 70 genes were found differentially expressed, by microarray analysis, when comparing LGESS showing ESS‐related fusions with LGESS without such transcripts, underscores the biochemical importance of the observed genetic heterogeneity and hints that new subgroups/entities in LGESS still remain undiscovered. © 2016 The Authors. Genes, Chromosomes & Cancer Published by Wiley Periodicals, Inc.

Overexpression of the oncostatin-M receptor in cervical squamous cell carcinoma is associated with epithelial-mesenchymal transition and poor overall survival

Background:

Copy-number gain of the oncostatin-M receptor (OSMR) occurs frequently in cervical squamous cell carcinoma (SCC) and is associated with adverse clinical outcome. We previously showed that OSMR overexpression renders cervical SCC cells more sensitive to the major ligand oncostatin-M (OSM), which increases migration and invasion in vitro. We hypothesised that a major contribution to this phenotype would come from epithelial–mesenchymal transition (EMT).

Methods:

We performed a comprehensive integrated study, involving in vitro cell line studies, in vivo animal models and numerous clinical samples from a variety of anatomical sites.

Results:

In independent sets of cervical, head/neck and lung SCC tissues, OSMR expression levels correlated with multiple EMT-associated phenotypic markers and transcription factors. OSM treatment of OSMR overexpressing cervical SCC cells produced consistent EMT changes and increased tumour sphere formation in suspension culture. In a mouse model, OSMR overexpressing SCC cells treated with OSM showed significant increases in lung colonisation. The biological effects of exogenous OSM were mirrored by highly significant adverse overall survival in cervical SCCs with OSMR overexpression (N=251).

Conclusions:

OSM:OSMR interactions are able to induce EMT, increased cancer stem cell-like properties and enhanced lung colonisation in SCC cells. These changes are likely to contribute to the highly significant adverse outcome associated with OSMR overexpression in cervical SCCs.

An Integrative Approach Uncovers Biomarkers that Associate with Clinically Relevant Disease Outcomes in Vulvar Carcinoma

Vulvar squamous cell carcinoma (VSCC) is a rare disease that has a high mortality rate (∼40%). However, little is known about its molecular signature. Therefore, an integrated genomics approach, based on comparative genome hybridization (aCGH) and genome-wide expression (GWE) array, was performed to identify driver genes in VSCC. To achieve that, DNA and RNA were extracted from frozen VSCC clinical specimens and examined by aCGH and GWE array, respectively. On the basis of the integration of data using the CONEXIC algorithm, PLXDC2 and GNB3 were validated by RT-qPCR. The expression of these genes was then analyzed by IHC in a large set of formalin-fixed paraffin-embedded specimens. These analyses identified 47 putative drivers, 46 of which were characterized by copy number gains that were concomitant with overexpression and one with a copy number loss and downregulation. Two of these genes, PLXDC2 and GNB3, were selected for further validation: PLXDC2 was downregulated and GNB3 was overexpressed compared with non-neoplastic tissue. By IHC, both proteins were ubiquitously expressed throughout vulvar tissue. High expression of GNB3 and low PLXDC2 immunostaining in the same sample was significantly associated with less lymph node metastasis and greater disease-free survival. On the basis of a robust methodology never used before for VSCC evaluation, two novel prognostic markers in vulvar cancer are identified: one with favorable prognosis (GNB3) and the other with unfavorable prognosis (PLXDC2).

IMPLICATIONS

This genomics study reveals markers that associate with prognosis and may provide guidance for better treatment in vulvar cancer. Mol Cancer Res; 14(8); 720-9. ©2016 AACR.

Olfactomedin 4 deficiency promotes prostate neoplastic progression and is associated with upregulation of the hedgehog-signaling pathway

Loss of olfactomedin 4 (OLFM4) gene expression is associated with the progression of human prostate cancer, but its role and the molecular mechanisms involved in this process have not been completely understood. In this study, we found that Olfm4-knockout mice developed prostatic intraepithelial neoplasia and prostatic adenocarcinoma. Importantly, we found that the hedgehog-signaling pathway was significantly upregulated in the Olfm4-knockout mouse model. We also found that restoration of OLFM4 in human prostate-cancer cells that lack OLFM4 expression significantly downregulated hedgehog signaling-pathway component expression. Furthermore, we demonstrated that the OLFM4 protein interacts with sonic hedgehog protein, as well as significantly inhibits GLI-reporter activity. Bioinformatic and immunohistochemistry analyses revealed that decreased OLFM4 and increased SHH expression was significantly associated with advanced human prostate cancer. Thus, olfactomedin 4 appears to play a critical role in regulating progression of prostate cancer, and has potential as a new biomarker for prostate cancer.

Investigating inter-chromosomal regulatory relationships through a comprehensive meta-analysis of matched copy number and transcriptomics data sets

Background

Gene regulatory relationships can be inferred using matched array comparative genomics and transcriptomics data sets from cancer samples. The way in which copy numbers of genes in cancer samples are often greatly disrupted works like a natural gene amplification/deletion experiment. There are now a large number of such data sets publicly available making a meta-analysis of the data possible.

Results

We infer inter-chromosomal acting gene regulatory relationships from a meta-analysis of 31 publicly available matched array comparative genomics and transcriptomics data sets in humans. We obtained statistically significant predictions of target genes for 1430 potential regulatory genes. The regulatory relationships being inferred are either direct relationships, of a transcription factor on its target, or indirect ones, through pathways containing intermediate steps. We analyse the predictions in terms of cocitations, both publications which cite a regulator with any of its inferred targets and cocitations of any genes in a target list.

Conclusions

The most striking observation from the results is the greater number of inter-chromosomal regulatory relationships involving repression compared to those involving activation. The complete results of the meta-analysis are presented in the database METAMATCHED. We anticipate that the predictions contained in the database will be useful in informing experiments and in helping to construct networks of regulatory relationships.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2100-5) contains supplementary material, which is available to authorized users.

Loss-of-Function PTPRD Mutations Lead to Increased STAT3 Activation and Sensitivity to STAT3 Inhibition in Head and Neck Cancer

BACKGROUND

Protein tyrosine phosphatase receptor type D (PTPRD) is a putative tumor suppressor in several cancers including head and neck squamous cell carcinoma (HNSCC). STAT3 is a frequently hyperactivated oncogene in HNSCC. As STAT3 is a direct substrate of PTPRD, we sought to determine the genetic or epigenetic alterations of PTPRD that contribute to overactive STAT3 in HNSCC.

METHODS

We analyzed data from The Cancer Genome Atlas (TCGA) and our previous whole-exome sequencing study and summarized the mutation, methylation, and copy number status of PTPRD in HNSCC and other cancers. In vitro studies involved standard transfection and MTT protocols, as well as methylation-specific PCR.

RESULTS

Our findings indicate that PTPRD mutation, rather than methylation or copy number alteration, is the primary mechanism by which PTPRD function is lost in HNSCC. We demonstrate that overexpression of wild-type PTPRD in HNSCC cells significantly inhibits growth and STAT3 activation while PTPRD mutants do not, suggesting that mutation may lead to loss of function and subsequent hyper-phosphorylation of PTPRD substrates, especially STAT3. Importantly, we determined that HNSCC cells harboring an endogenous PTPRD mutation are more sensitive to STAT3 blockade than PTPRD wild-type cells. We additionally found that PTPRD mRNA expression does not correlate with pSTAT3 expression, suggesting that alterations that manifest through altered mRNA expression, including hypermethylation and gene copy number alterations, do not significantly contribute to STAT3 overactivation in HNSCC.

CONCLUSION

PTPRD mutation, but not methylation or copy number loss, may serve as a predictive biomarker of sensitivity to STAT3 inhibitors in HNSCC.

Deletion of Ptprd and Cdkn2a cooperate to accelerate tumorigenesis

PTPRD encodes the protein tyrosine phosphatase receptor type D and is frequently inactivated across many human cancers. Despite its frequent inactivation, it is unknown whether loss of PTPRD promotes tumorigenesis in vivo. PTPRD is located on chromosome 9p, as is CDKN2A, and the two loci are frequently deleted together. Here, we show that co-deletion of Ptprd and Cdkn2a cooperate to accelerate tumorigenesis. Interestingly, heterozygous loss of Ptprd was sufficient to promote tumorigenesis in our model, suggesting that Ptprd may be a haploinsufficient tumor suppressor. The loss of Ptprd resulted in changes to the tumor spectrum in mice and increased the frequency of lymphomas. In total, we reveal that Ptprd is a tumor suppressor that can promote tumorigenesis in concert with Cdkn2a loss.

HMGA2 expression pattern and TERT mutations in tumors of the vulva

Malignant tumors of the vulva account for only 5% of cancers of the female genital tract in the USA. The most frequent cancers of the vulva are squamous cell carcinoma (SCC) and malignant melanoma (MM). Little is known about the genetic aberrations carried by these tumors. We report a detailed study of 25 vulva tumors [22 SCC, 2 MM, 1 atypical squamous cell hyperplasia (AH)] analyzed for expression of the high-mobility group AT-hook family member genes HMGA2 and HMGA1, for mutations in the IDH1, IDH2 and TERT genes, and for methylation of the MGMT promoter. The RT-PCR and immunohistochemistry analyses showed that HMGA2 was expressed in the great majority of analyzed samples (20 out of 24; SCC as well as MM) but not in the normal controls. HMGA1, on the other hand, was expressed in both tumors and normal tissues. Five of the 24 tumors (all SCC) showed the C228T mutation in the TERT promoter. Our results showed that HMGA2 and TERT may be of importance in the genesis and/or the progression of tumors of the vulva.

Genetic and epigenetic changes in vulvar squamous cell carcinoma and its precursor lesions: a review of the current literature

Vulvar cancer is a relatively rare gynecologic malignancy with an annual incidence in developed countries of approximately 2 per 100,000 women. Vulvar squamous cell carcinoma (VSCC) has two etiological pathways: a high risk human papillomavirus (HPV)-dependent route, which has usual vulvar intraepithelial neoplasia (uVIN) as a precursor lesion, and an HPV-independent route, which is associated with differentiated VIN (dVIN), lichen sclerosus, and genetic alterations, such as TP53 mutations. Research on the molecular etiology of vulvar cancer has increased in the past years, not only regarding genetic alterations, but also epigenetic changes. In genetic alterations, a mutation irreversibly changes the nucleotide sequence of the DNA, or the number of copies of chromosomes per cell is altered. In epigenetics, the nucleotide sequence remains the same but genes can be 'switched' on or off by, for example, DNA methylation or histone modification. We searched the current literature on genetic and epigenetic alterations in VSCC and its precursor lesions. Many studies have reported a higher incidence of somatic mutations in HPV-negative tumors compared to HPV-positive tumors, with TP53 mutations being the most frequent. Allelic imbalances or loss of heterozygosity are more frequently found in higher stages of dysplasia and in invasive carcinomas, but it is not exclusive to HPV-negative tumors. A limited number of studies are available on epigenetic changes in vulvar lesions, with hypermethylation of CDKN2A being the most frequently investigated change. For most genes, hypermethylation occurs more frequently in vulvar squamous cell carcinomas than in precursor lesions. As most studies have focused on HPV infection and TP53 mutations, we suggest that more research should be performed using whole genome or next generation sequencing to determine the true landscape of genetic and epigenetic alterations in vulvar squamous cell carcinoma.

Receptor-type protein tyrosine phosphatases in cancer

Protein tyrosine phosphatases (PTPs) play an important role in regulating cell signaling events in coordination with tyrosine kinases to control cell proliferation, apoptosis, survival, migration, and invasion. Receptor-type protein tyrosine phosphatases (PTPRs) are a subgroup of PTPs that share a transmembrane domain with resulting similarities in function and target specificity. In this review, we summarize genetic and epigenetic alterations including mutation, deletion, amplification, and promoter methylation of PTPRs in cancer and consider the consequences of PTPR alterations in different types of cancers. We also summarize recent developments using PTPRs as prognostic or predictive biomarkers and/or direct targets. Increased understanding of the role of PTPRs in cancer may provide opportunities to improve therapeutic approaches.

A meta-analysis of multiple matched copy number and transcriptomics data sets for inferring gene regulatory relationships

Inferring gene regulatory relationships from observational data is challenging. Manipulation and intervention is often required to unravel causal relationships unambiguously. However, gene copy number changes, as they frequently occur in cancer cells, might be considered natural manipulation experiments on gene expression. An increasing number of data sets on matched array comparative genomic hybridisation and transcriptomics experiments from a variety of cancer pathologies are becoming publicly available. Here we explore the potential of a meta-analysis of thirty such data sets. The aim of our analysis was to assess the potential of in silico inference of trans-acting gene regulatory relationships from this type of data. We found sufficient correlation signal in the data to infer gene regulatory relationships, with interesting similarities between data sets. A number of genes had highly correlated copy number and expression changes in many of the data sets and we present predicted potential trans-acted regulatory relationships for each of these genes. The study also investigates to what extent heterogeneity between cell types and between pathologies determines the number of statistically significant predictions available from a meta-analysis of experiments.

Chromosomal gains and losses in human papillomavirus-associated neoplasia of the lower genital tract - a systematic review and meta-analysis

BACKGROUND

Overexpression of the human papillomavirus (HPV) oncogenes E6 and E7 is necessary for the development of distinct lower genital tract cancers. However, secondary cellular genomic alterations are mandatory to promote progression of HPV-induced premalignant stages. We aimed at identifying the chromosomal regions most frequently gained and lost and the disease stage at which the latter occurs. These regions might be relevant for carcinogenesis and could serve as diagnostic markers to identify premalignant lesions with high progression risk towards invasive cancer.

METHODS

We performed a systematic literature review and meta-analysis of studies listed in PubMed that analysed chromosomal copy number alterations by comparative genomic hybridisation (CGH) in HPV-positive and -negative cancers or premalignant lesions of the anogenital tract (cervix, anus, vagina, penis and vulva).

FINDINGS

Data were extracted and analysed from 32 studies. The most common alterations in cervical squamous cell carcinoma (SCC) (12 studies, 293 samples) were gains at 3q with a rate of 0.55 (95% confidence interval (CI) 0.43-0.70), losses at 3p (0.36, 95%CI 0.27-0.48) and losses at 11q (0.33, 95%CI 0.26-0.43). Gains at 3q were particularly frequent in HPV16-positive cervical SCC (0.84, 95%CI 0.78-0.90). Also more than one quarter of high grade cervical intraepithelial neoplasia (CIN) harboured gains of 3q (0.27, 95%CI 0.20-0.36), but the rate in low grade CIN was low (0.02, 95%CI 0.00-0.09). For HPV-associated vulvar SCC (four studies, 30 samples) the same common alterations as in cervical SCC were reported. Studies on non-cervical and non-vulvar SCC and premalignant lesions of the lower genital tract are scarce.

INTERPRETATION

3q gains were most frequently found in HPV16-positive cervical SCC. The results suggest the selection of HPV-transformed cell clones harbouring 3q gains in high grade premalignant lesions, while alterations in low grade lesions are rare.