Design and experimental validation of OPERA_MET-A panel for deep methylation analysis by next generation sequencing.

DNA methylation is the most recognized epigenetic mark that leads to a massive distortion in cancer cells. It has been observed that a large number of DNA aberrant methylation events occur simultaneously in a group of genes, thus providing a growth advantage to the cell in promoting cell differentiation and neoplastic transformation. Due to this reason, methylation profiles have been suggested as promising cancer biomarkers. Here, we designed and performed a first step of validation of a novel targeted next generation sequencing (NGS) panel for methylation analysis, which can simultaneously evaluate the methylation levels at CpG sites of multiple cancer-related genes. The OPERA_MET-A methylation panel was designed using the Ion AmpliSeq™ technology to amplify 155 regions with 125-175 bp mean length and covers a total of 1107 CpGs of 18 cancer-related genes. The performance of the panel was assessed by running commercially available fully methylated and unmethylated control human genomic DNA (gDNA) samples and a variable mixture of them. The libraries were run on Ion Torrent platform and the sequencing output was analyzed using the "methylation_analysis" plugin. DNA methylation calls on both Watson (W) and Crick (C) strands and methylated:unmethylated ratio for each CpG site were obtained. Cell lines, fresh frozen and formalin-fixed paraffin-embedded (FFPE) lung cancer tissues were tested. The OPERA_MET-A panel allows to run a minimum of 6 samples/530 chip to reach an observed mean target depth ≥2,500X (W and C strands) and an average number of mapped reads >750,000/sample. The conversion efficiency, determined by spiking-in unmethylated Lambda DNA into each sample before the bisulfite conversion process, was >97% for all samples. The observed percentage of global methylation for all CpGs was >95% and <5% for fully methylated and unmethylated gDNA samples, respectively, and the observed results for the variable mixtures were in agreement with what was expected. Methylation-specific NGS analysis represents a feasible method for a fast and multiplexed screening of cancer patients by a high-throughput approach. Moreover, it offers the opportunity to construct a more robust algorithm for disease prediction in cancer patients having a low quantity of biological material available.

Novel association of MEN1 gene mutations with parathyroid carcinoma.

Inactivating mutations of the multiple endocrine neoplasia 1 (MEN1) gene cause MEN1 syndrome, characterized by primary hyperparathyroidism (pHPT), and parathyroid and gastro-entero-pancreatic pituitary tumors. At present, only 14 cases of malignant parathyroid tumor have been associated with the syndrome, with 6 cases carrying an inactivating mutation of the MEN1 gene. The present study presents the case of a 48-year-old female who presented with multigland pHPT and multiple pancreatic lesions. The patient underwent surgery several times for the excision of parathyroid hyperplasia, carcinoma and adenoma. The MEN1 gene was screened, revealing three variants (in cis) at the intron/exon 3 boundary (IVS2-3G>C, c.497A>T and c.499G>T) detected on the DNA of the proband, not shared by her relatives. RNA sequencing revealed that the IVS2-3C>G variant caused the skipping of the exon 3. Therefore, the present study reports on a novel rare association of MEN1 syndrome and parathyroid carcinoma. The reported splicing mutation was previously identified in subjects who always developed malignant lesions; thus, a possible genotype-phenotype association may be considered.

Stepwise analysis of MIR9 loci identifies miR-9-5p to be involved in Oestrogen regulated pathways in breast cancer patients.

miR-9 was initially identified as an epigenetically regulated miRNA in tumours, but inconsistent findings have been reported so far. We analysed the expression of miR-9-5p, miR-9-3p, pri-miRs and MIR9 promoters methylation status in 131 breast cancer cases and 12 normal breast tissues (NBTs). The expression of both mature miRs was increased in tumours as compared to NBTs (P < 0.001) and negatively correlated with ER protein expression (P = 0.005 and P = 0.003, for miR-9-3p and miR-9-5p respectively). In addition, miR-9-5p showed a significant negative correlation with PgR (P = 0.002). Consistently, miR-9-5p and miR-9 3p were differentially expressed in the breast cancer subgroups identified by ER and PgR expression and HER2 amplification. No significant correlation between promoter methylation and pri-miRNAs expressions was found either in tumours or in NBTs. In the Luminal breast cancer subtype the expression of miR-9-5p was associated with a worse prognosis in both univariable and multivariable analyses. Ingenuity Pathway Analysis exploring the putative interactions among miR-9-5p/miR-9-3p, ER and PgR upstream and downstream regulators suggested a regulatory loop by which miR-9-5p but not miR-9-3p is induced by steroid hormone receptor and acts within hormone-receptor regulated pathways.

Competitive allele-specific TaqMan PCR (Cast-PCR) is a sensitive, specific and fast method for BRAF V600 mutation detection in Melanoma patients.

BRAF codon 600 mutation testing of melanoma patients is mandatory for the choice of the most appropriate therapy in the clinical setting. Competitive allele specific TaqMan PCR (Cast-PCR) technology allows not only the selective amplification of minor alleles, but it also blocks the amplification of non-mutant allele. We genotyped codon 600 of the BRAF gene in 54 patients' samples by Cast-PCR and bidirectional direct sequence analysis. All the mutations detected by sequencing were also identified by Cast-PCR. In addition, Cast-PCR assay detected four samples carrying mutations and was able to clearly identify two mutations of uncertain interpretation by Sanger sequencing. The limit of detection of Cast-PCR was evaluated by constructing dilution curves of BRAF(V600E) and BRAF(V600K) mutated clinical samples mixed with a not-mutated specimens. Both mutations could be detected until a 1:100 mutated/not mutated ratio. Cloning and sequencing of the clones was used to confirm mutations on representative discrepant cases. Cast PCR performances were not affected by intratumour heterogeneity, and less affected by melanin content. Our results indicate that Cast-PCR is a reliable diagnostic tool for the identification of melanoma patients as eligible to be treated with TKIs and might be implemented in the clinical setting as elective screening method.

A miRNA signature for defining aggressive phenotype and prognosis in gliomas.

Gliomas represent a disparate group of tumours for which there are to date no cure. Thus, there is a recognized need for new diagnostic and therapeutic approaches based on increased understanding of their molecular nature. We performed the comparison of the microRNA (miRNA) profile of 8 WHO grade II gliomas and 24 higher grade tumours (2 WHO grade III and 22 glioblastomas) by using the Affymetrix GeneChip miRNA Array v. 1.0. A relative quantification method (RT-qPCR) with standard curve was used to confirm the 22 miRNA signature resulted by array analysis. The prognostic performances of the confirmed miRNAs were estimated on the Tumor Cancer Genome Atlas (TCGA) datasets. We identified 22 miRNAs distinguishing grade II gliomas from higher grade tumours. RT-qPCR confirmed the differential expression in the two patients' groups for 13 out of the 22 miRNAs. The analysis of the Glioblastoma Multiforme (GBM) and Lower Grade Glioma (LGG) datasets from TCGA demonstrated the association with prognosis for 6 of those miRNAs. Moreover, in the GBM dataset miR-21 and miR-210 were predictors of worse prognosis in both univariable and multivariable Cox regression analyses (HR 1.19, p = 0.04, and HR 1.18, p = 0.029 respectively). Our results support a direct contribution of miRNAs to glioma cancerogenesis and suggest that miR-21 and miR-210 may play a role in the aggressive clinical behaviour of glioblastomas.

FOXP1 and TP63 involvement in the progression of myelodysplastic syndrome with 5q- and additional cytogenetic abnormalities.

BACKGROUND: The progression of low-risk del(5q) myelodysplastic syndrome to acute myeloid leukemia is increased when associated with mutations of TP53, or with additional chromosomal abnormalities. However, to date the prognostic impact and molecular consequences of these rearrangements were poorly investigated. Single additional alterations to del(5q) by balanced chromosome rearrangements were rarely found in myelodysplasia. In particular, balanced alterations involving TP63 and FOXP1 genes were never reported in the literature. CASE PRESENTATION: Here we report on a 79-year woman with an aggressive form of myelodysplastic syndrome with del(5q), no TP53 mutation, and a novel complex rearrangement of chromosome 3 in bone marrow cells. Our results revealed that the FOXP1 and TP63 genes were both relocated along chromosome 3. Strikingly, immunohistochemistry analysis showed altered protein levels, disclosing that this rearrangement triggered the expression of FOXP1 and TP63 genes. FOXP1 was also found activated in other patients with myelodysplasia and acute myeloid leukemia, showing that it is an important, recurrent event. CONCLUSIONS: We document an apparent role of FOXP1 and TP63, up to now poorly documented, in the progression of MDS in our patient who is lacking mutations in the TP53 tumor suppressor gene normally associated with poor outcome in myelodysplastic syndrome with 5q-. Finally, our results may suggest a possible broader role of FOXP1 in the pathogenesis and progression of myelodysplasia and acute myeloid leukemia.

Evaluation of microRNA-10b prognostic significance in a prospective cohort of breast cancer patients.

BACKGROUND: MicroRNA-10b (miR-10b) has a prominent role in regulating tumor invasion and metastasis by targeting the HOXD10 transcriptional repressor and has been found up-regulated in several tumor types. METHODS: We evaluated the expression of miR-10b in paired tumor and normal specimens obtained from a prospective cohort of breast cancer patients with at least 36 months follow-up enrolled according to the REMARK guidelines (n = 150). RNA quality was measured and only samples with RNA Integrity Number (RIN) ≥7.0 were analyzed. RESULTS: The relative expression of miR-10b in tumor as compared to its normal counterpart (RER) was determined by RT-qPCR. miR-10b RERs were higher in the subgroup of patients with synchronous metastases (n = 11, Median 0.25; IQR 0.11-1.02) as compared with patients without metastases (n = 90, Median 0.09; IQR 0.04-0.29) (p = 0.028). In the subgroup of patients without synchronous metastases (n = 90), higher miR-10b RERs were associated with increased risk of disease progression and death in both univariable (HR 1.16, p = 0.021 and HR 1.20, p = 0.015 respectively for 0.10 unitary increase of miR-10b RERs levels) and multivariable (HR1.30, p < 0.001, and HR 1.31, p = 0.003 respectively for 0.10 unitary increase of miR-10b RERs levels) Cox regression models. The addition of miR-10b RERs to the Nottingham Prognostic Index (NPI) provided an improvement in discrimination power and risk reclassification abilities for the clinical outcomes at 36 months. Survival C-indices significantly increased from 0.849 to 0.889 (p = 0.009) for OS and from 0.735 to 0.767 (p = 0.050) for DFS. CONCLUSIONS: Our results provide evidences that the addition of miR-10b RERs to the prognostic factors used in clinical routine could improve the prediction abilities for both overall mortality and disease progression in breast cancer patients.

Molecular dissection of the VHL gene in solitary capillary hemangioblastoma of the central nervous system.

Capillary hemangioblastomas (HGBs) of the CNS occur either sporadically or as part of the von Hippel-Lindau (VHL) syndrome. Molecular characterizations of the VHL gene in sporadic HGBs at the somatic level have been limited to date. We investigated the VHL gene in 57 patients most of whom (55 [96%] of 57) had a solitary CNS HGB at the time of surgery. Tissues from 23 HGBs of these patients (2 VHL related and 21 unrelated) were also investigated at genetic and epigenetic levels. Two of the 51 patients with apparently sporadic HGBs and no additional evidence of VHL (∼4%) were found to have a germline VHL gene mutation; both of these patients subsequently developed evidence of VHL syndrome. Somatic VHL gene mutations were found in 11 (52%) of the 21 non-VHL-related cases. A germline mutation was identified in 5 (84%) of 6 VHL-associated HGBs; double gene inactivation was observed in tumor tissue from VHL syndrome patients. Seven different previously unreported VHL gene alterations (6 somatic and 1 germline) were identified; double hits were identified in 7 (12%) of 57 cases. Our findings confirm the usefulness of VHL gene analysis at the germline level in patients who present with apparently solitary HGB. Moreover, the genetic and epigenetic VHL gene investigations performed support a key role for functional alterations of the VHL gene in sporadic neuraxial HGB.

VHL gene alterations in Italian patients with isolated renal cell carcinomas.

Clear cell renal cell carcinoma (ccRCC) is the most common malignant neoplasm of the kidney and belongs to the few human tumors known to develop from mutations of the VHL tumor suppressor gene. VHL germline mutations are associated with hereditary ccRCCs in VHL disease. However, somatic VHL gene defects may also occur in sporadic ccRCCs. In this study, we analyzed the frequency and the spectrum of VHL gene alterations in 35 Italian patients with sporadic renal cell carcinoma (RCC). Tumor-specific intragenic VHL pathogenic mutations were detected in 38% (11/29) of the ccRCC patients and 33% (2/6) of the patients with other types of RCC. One novel 18-bp in-tandem duplication and 4 previously unreported nucleotide changes in the VHL gene were described. Microsatellite analysis showed loss of heterozygosity for at least 1 informative marker in 43% (9/21) of the ccRCCs and 50% (3/6) of the non-ccRCCs; 5 of the 13 tumors (38%) harboring VHL gene alterations also had loss of heterozygosity for at least 1 microsatellite marker. Our results confirm that somatic inactivation of the VHL gene may play a pivotal role in the tumorigenesis of sporadic ccRCCs in Italian patients and suggests that mutation analysis of the VHL gene may be helpful for discriminating sporadic, VHL-gene-related ccRCCs from those related to VHL disease.

20 novel point mutations and one large deletion in EXT1 and EXT2 genes: report of diagnostic screening in a large Italian cohort of patients affected by hereditary multiple exostosis.

BACKGROUND: Hereditary multiple exostosis represents the most frequent bone tumor disease in humans. It consists of cartilage deformities affecting the juxta-ephyseal region of long bones. Usually benign, exostosis could degenerate in malignant chondrosarcoma form in less than 5% of the cases. Being caused by mutations in the predicted tumor suppressor genes, EXT1 (chr 8q23-q24) and EXT2 (chr 11p11-p12) genes, HMEs are usually inherited with an autosomal dominant pattern, although "de novo" cases are not infrequent. AIM: Here we present our genetic diagnostic report on the largest Southern Italy cohort of HME patients consisting of 90 subjects recruited over the last 5years. RESULTS: Molecular screening performed by direct sequencing of both EXT1 and EXT2 genes, by MLPA and Array CGH analyses led to the identification of 66 mutations (56 different occurrences) and one large EXT2 deletion out of 90 patients (74.4%). The total of 21 mutations (20 different occurrences, 33.3%) and the EXT2 gene deletion were novel. In agreement with literature data, EXT1 gene mutations were scattered along all the protein sequence, while EXT2 lesions fell in the first part of the protein. Conservation, damaging prediction and 3-D modeling, in-silico, analyses, performed on three novel missense variants, confirmed that at least in two cases the novel aminoacidic changes could alter the structure stability causing a strong protein misfolding. CONCLUSIONS: Here we present 20 novel EXT1/EXT2 mutations and one large EXT2 deletion identified in the largest Southern Italy cohort of patients affected by hereditary multiple exostosis.

Aberrant Keap1 methylation in breast cancer and association with clinicopathological features.

Keap1 (Kelch-like ECH-associated protein 1) is an adaptor protein that mediates the ubiquitination/degradation of genes regulating cell survival and apoptosis under oxidative stress conditions. We determined methylation status of the KEAP1 promoter in 102 primary breast cancers, 14 pre-invasive lesions, 38 paired normal breast tissues and 6 normal breast from reductive mammoplasty by quantitative methylation specific PCR (QMSP). Aberrant promoter methylation was detected in 52 out of the 102 primary breast cancer cases (51%) and 10 out of 14 pre-invasive lesions (71%). No mutations of the KEAP1 gene were identified in the 20 breast cancer cases analyzed by fluorescence based direct sequencing. Methylation was more frequent in the subgroup of patients identified as ER positive-HER2 negative tumors (66.7%) as compared with triple-negative breast cancers (35%) (p = 0.05, Chi-square test). The impact of the interactions between Er, PgR, Her2 expression and KEAP1 methylation on mortality was investigated by RECPAM multivariable statistical analysis, identifying four prognostic classes at different mortality risks. Triple-negative breast cancer patients with KEAP1 methylation had higher mortality risk than patients without triple-negative breast cancer (HR = 14.73, 95%CI: 3.65-59.37). Both univariable and multivariable COX regressions analyses showed that KEAP1 methylation was associated with a better progression free survival in patients treated with epirubicin/cyclophosfamide and docetaxel as sequential chemotherapy (HR = 0.082; 95%CI: 0.007-0.934). These results indicate that aberrant promoter methylation of the KEAP1 gene is involved in breast cancerogenesis. In addition, identifying patients with KEAP1 epigenetic abnormalities may contribute to disease progression prediction in breast cancer patients.

Hypermethylated levels of E-cadherin promoter in Huh-7 cells expressing the HCV core protein.

BACKGROUND AND AIM: The mechanisms of hepatocarcinogenesis induced by hepatitis C virus remain unclear. Our aim was to investigate the effect of the HCV core protein on the promoter methylation status of selected genes potentially involved in the hepatocellular carcinoma (HCC). MATERIALS AND METHODS: We evaluated the promoter methylation levels of the E-cadherin (CDH1), the glutathione S-transferase p1 (GSTP1), adenomatosis polyposis coli (APC), tissue inhibitor of metalloproteinase 3 (TIMP3), catenin (cadherin-associated protein) beta 1 (CNNTB1) genes by a quantitative methylation-specific polymerase chain reaction (QMSP) in the in vitro model of Huh-7 cells expressing the HCV core protein of genotype 1b. RESULTS: We found that CDH1 promoter was hypermethylated in genotype 1b HCV core protein-positive cells as compared to control cells expressing the GFP protein alone (HCV core 1b vs GFP p=0.00; HCV core 1b vs Huh-7 p=0.03). This resulted in reduced levels of CDH1 protein as evaluated by immunoblot and by immunofluorescence. On the other hand no significant changes were observed for the other genes investigated. Furthermore, we present evidence that genotype 1b HCV core protein expression induces SIRT1 upregulation and that treatment with SIRT1 inhibitor sirtinol decreases the methylation levels of CDH1 promoter (1b+sirtinol vs 1b p=0.05; 1b+sirtinol vs GFP+sirtinol p=NS) resulting in 1.7-fold increased CDH1 mRNA expression (1b+sirtinol vs 1b p=0.05). CONCLUSIONS: Our findings suggest that HCV core protein could play a role in HCC at least in part by altering the methylation status of CDH1 promoter. These findings could also suggest a novel therapeutic approach for HCC.

Frequent epigenetics inactivation of KEAP1 gene in non-small cell lung cancer.

The KEAP1/Nrf2 pathway is a master regulator of several redox-sensitive genes implicated in resistance of tumor cells against chemotherapeutic drugs. Recent data suggest that epigenetic mechanisms may play a pivotal role in the regulation of KEAP1 expression. We performed a comprehensive genetic and epigenetic analysis of the KEAP1 gene in 47 non-small cell lung cancer tissues and normal specimens. Promoter methylation analysis was performed using a quantitative methylation specific PCR assay in real time. Methylation at the KEAP1 promoter region was detected in 22 out of the 47 NSCLCs (47%) and in none of the normal tissues analyzed. Somatic mutations were detected in 7 out of the 47 tumors (15%) and loss of heterozygosity (LOH) in 10 out of the 47 (21%) of the cases. Overall, we found at least one molecular alteration in 57% of the cases. Approximately one third of the tumors had two alterations and this feature was associated with higher risk of disease progression in univariate COX regression analysis (HR = 3.62; 95% CI 1.24-10.65, p = 0.02). This result was confirmed by Kaplan-Meier analysis, which demonstrated an association between worst outcome and KEAP1 double alterations (p = 0.01, Log rank test). Our results further suggest that deregulation of the NRF2/KEAP1 system could play a pivotal role in the cancerogenesis of NSCLC. In addition identifying patients with KEAP1 genetic and epigenetic abnormalities may contribute to disease progression prediction and response to therapy in lung cancer patients.

Regulation of KEAP1 expression by promoter methylation in malignant gliomas and association with patient's outcome.

In light with the view that KEAP1 loss of function may impact tumour behavior and modify response to chemotherapeutical agents, we sought to determine whether KEAP1 gene is epigenetically regulated in malignant gliomas. We developed a Quantitative Methylation Specific PCR (QMSP) assay to analyze 86 malignant gliomas and 20 normal brain tissues. The discriminatory power of the assay was assessed by Receiving Operating Characteristics (ROC) curve analysis. The AUC value of the curve was 0.823 (95%CI: 0.764-0.883) with an optimal cut off value of 0.133 yielding a 74% sensitivity (95%CI: 63%-82%) and an 85% specificity (95%CI: 64%-95%). Bisulfite sequencing analysis confirmed QMSP results and demonstrated a direct correlation between percentage of methylated CpGs and methylation levels (Spearman's Rho 0.929, P=0.003). Remarkably, a strong inverse correlation was observed between methylation levels and KEAP1 mRNA transcript in tumour tissue (Spearman's Rho -0.656 P=0.0001) and in a cell line before and after treatment with 5-azacytidine (P=0.003). RECPAM multivariate statistical analysis studying the interaction between MGMT and KEAP1 methylation in subjects treated with radiotherapy and temozolomide (n=70), identified three prognostic classes of glioma patients at different risk to progress. While simultaneous methylation of MGMT and KEAP1 promoters was associated with the lowest risk to progress, patients showing only MGMT methylation were the subgroup at the higher risk (HR 5.54, 95% CI 1.35-22.74). Our results further suggest that KEAP1 expression is epigenetically regulated. In addition we demonstrated that KEAP1 is frequently methylated in malignant gliomas and a predictor of patient's outcome.

Small deletion at the 7q21.2 locus in a CCM family detected by real-time quantitative PCR.

Cerebral cavernous malformations (CCMs) represent a common autosomal dominant disorder that predisposes patients to haemorrhagic strokes and focal neurological signs. About 56% of the hereditary forms of CCMs have been so far associated with mutations in the KRIT1 (Krev Interaction Trapped 1) gene, located at 7q21.2 (CCM1 locus). We described the complete loss of 7q21.2 locus encompassing the KRIT1 gene and 4 flanking genes in a CCM family by using a dense set of 12 microsatellite markers. The complete loss of the maternal copy of KRIT1 gene region was confirmed by Real-Time Quantitative Polymerase Chain Reaction (RT-QPCR) and the same approach was used for expression analysis. Additional RT-QPCR analysis showed the extension of the deletion, for a total of 700 kb, to the adjacent downstream and upstream-located genes, MTERF, AKAP9, CYP51A1, as well as a partial loss of the ANKIB1 gene. Here we report the molecular characterization of an interstitial small genomic deletion of the 7q21.2 region in a CCMs affected family, encompassing the KRIT1 gene. Our findings confirm the loss of function mechanism for the already known CCM1 locus, without any evident involvement of the other deleted genes. Moreover, our investigations highlight the usefulness of the RT-QPCR to the molecular characterization of the breakpoints genomic deletions and to the identification of internal deleted genes involved in the human genetic diseases.

Identification of two novel mutations and of a novel critical region in the KRIT1 gene.

Cerebral cavernous malformations (CCMs) represent a common autosomal dominant disorder that predisposes patients to hemorrhagic strokes and focal neurological signs. Mutations in three genes (KRIT1, MGC4607, and PDCD10) have been associated with CCMs. We investigated the role of two new mutations in the KRIT1 gene in two Italian families affected by CCMs. Whole blood DNA was extracted and the mutations were detected after polymerase chain reaction (PCR), denaturing high-performance liquid chromatography screening, and sequencing of the coding regions of the three CCMs-associated genes. Total RNA was extracted, and the KRIT1 cDNA was sequenced and subsequently subjected to real-time quantitative PCR in order to examine the translational outcome of each genomic mutation. A novel splicing acceptor site deletion of the exon 14 in one family and an intronic nucleotide change close to the exon 19 in the other one were identified, both in the KRIT1 gene. These mutations were proven to alter the correct splicing mechanism, resulting, respectively, in a truncated protein of 432 amino acids and in a protein lacking an internal segment. We report two novel cases of splicing affecting genomic variants, suggesting a careful reanalysis of previously identified splice site variations in KRIT1 to look for their possible causative roles of similar missplicing events and their consequent involvement in the pathogenesis of CCMs. Moreover, our genotype-phenotype functional correlation suggests that the C-terminal portion of the KRIT1 protein is likely to contain a short, previously unrecognized segment necessary for its activity.