Deletion of Chromosomes 13q and 14q Is a Common Feature of Tumors with BRCA2 Mutations

Introduction

Germline BRCA1 or BRCA2 mutations account for 20–30% of familial clustering of breast cancer. The main indication for BRCA2 screening is currently the family history but the yield of mutations identified in patients selected this way is low.

Methods

To develop more efficient approaches to screening we have compared the gene expression and genomic profiles of BRCA2-mutant breast tumors with those of breast tumors lacking BRCA1 or BRCA2 mutations.

Results

We identified a group of 66 genes showing differential expression in our training set of 7 BRCA2-mutant tumors and in an independent validation set of 19 BRCA2-mutant tumors. The differentially expressed genes include a prominent cluster of genes from chromosomes 13 and 14 whose expression is reduced. Gene set enrichment analysis confirmed that genes in specific bands on 13q and 14q showed significantly reduced expression, suggesting that the affected bands may be preferentially deleted in BRCA2-mutant tumors. Genomic profiling showed that the BRCA2-mutant tumors indeed harbor deletions on chromosomes 13q and 14q. To exploit this information we have created a simple fluorescence in situ hybridization (FISH) test and shown that it detects tumors with deletions on chromosomes 13q and 14q.

Conclusion

Together with previous reports, this establishes that deletions on chromosomes 13q and 14q are a hallmark of BRCA2-mutant tumors. We propose that FISH to detect these deletions would be an efficient and cost-effective first screening step to identify potential BRCA2-mutation carriers among breast cancer patients without a family history of breast cancer.     

A transcriptomic computational analysis of mastic oil-treated Lewis lung carcinomas reveals molecular mechanisms targeting tumor cell growth and survival

Background

DNA copy number alterations are frequently observed in ovarian cancer, but it remains a challenge to identify the most relevant alterations and the specific causal genes in those regions.

Methods

We obtained high-resolution 500K SNP array data for 52 ovarian tumors and identified the most statistically significant minimal genomic regions with the most prevalent and highest-level copy number alterations (recurrent CNAs). Within a region of recurrent CNA, comparison of expression levels in tumors with a given CNA to tumors lacking that CNA and to whole normal ovary samples was used to select genes with CNA-specific expression patterns. A public expression array data set of laser capture micro-dissected (LCM) non-malignant fallopian tube epithelia and LCM ovarian serous adenocarcinoma was used to evaluate the effect of cell-type mixture biases.

Results

Fourteen recurrent deletions were detected on chromosomes 4, 6, 9, 12, 13, 15, 16, 17, 18, 22 and most prevalently on X and 8. Copy number and expression data suggest several apoptosis mediators as candidate drivers of the 8p deletions. Sixteen recurrent gains were identified on chromosomes 1, 2, 3, 5, 8, 10, 12, 15, 17, 19, and 20, with the most prevalent gains localized to 8q and 3q. Within the 8q amplicon, PVT1, but not MYC, was strongly over-expressed relative to tumors lacking this CNA and showed over-expression relative to normal ovary. Likewise, the cell polarity regulators PRKCI and ECT2 were identified as putative drivers of two distinct amplicons on 3q. Co-occurrence analyses suggested potential synergistic or antagonistic relationships between recurrent CNAs. Genes within regions of recurrent CNA showed an enrichment of Cancer Census genes, particularly when filtered for CNA-specific expression.

Conclusion

These analyses provide detailed views of ovarian cancer genomic changes and highlight the benefits of using multiple reference sample types for the evaluation of CNA-specific expression changes.     

Distinct Genetic Alterations in Colorectal Cancer

Background

Colon cancer (CRC) development often includes chromosomal instability (CIN) leading to amplifications and deletions of large DNA segments. Epidemiological, clinical, and cytogenetic studies showed that there are considerable differences between CRC tumors from African Americans (AAs) and Caucasian patients. In this study, we determined genomic copy number aberrations in sporadic CRC tumors from AAs, in order to investigate possible explanations for the observed disparities.

Methodology/Principal Findings

We applied genome-wide array comparative genome hybridization (aCGH) using a 105k chip to identify copy number aberrations in samples from 15 AAs. In addition, we did a population comparative analysis with aCGH data in Caucasians as well as with a widely publicized list of colon cancer genes (CAN genes). There was an average of 20 aberrations per patient with more amplifications than deletions. Analysis of DNA copy number of frequently altered chromosomes revealed that deletions occurred primarily in chromosomes 4, 8 and 18. Chromosomal duplications occurred in more than 50% of cases on chromosomes 7, 8, 13, 20 and X. The CIN profile showed some differences when compared to Caucasian alterations.

Conclusions/Significance

Chromosome X amplification in male patients and chromosomes 4, 8 and 18 deletions were prominent aberrations in AAs. Some CAN genes were altered at high frequencies in AAs with EXOC4, EPHB6, GNAS, MLL3 and TBX22 as the most frequently deleted genes and HAPLN1, ADAM29, SMAD2 and SMAD4 as the most frequently amplified genes. The observed CIN may play a distinctive role in CRC in AAs.     

Highly Conserved Non-Coding Sequences and the 18q Critical Region for Short Stature: A Common Mechanism of Disease?

Background

Isolated growth hormone deficiency (IGHD) and multiple pituitary hormone deficiency (MPHD) are heterogeneous disorders with several different etiologies and they are responsible for most cases of short stature. Mutations in different genes have been identified but still many patients did not present mutations in any of the known genes. Chromosomal rearrangements may also be involved in short stature and, among others, deletions of 18q23 defined a critical region for the disorder. No gene was yet identified.

Methodology/Principal Findings

We now report a balanced translocation X;18 in a patient presenting a breakpoint in 18q23 that was surprisingly mapped about 500 Kb distal from the short stature critical region. It separated from the flanking SALL3 gene a region enriched in highly conserved non-coding elements (HCNE) that appeared to be regulatory sequences, active as enhancers or silencers during embryonic development.

Conclusion

We propose that, during pituitary development, the 18q rearrangement may alter expression of 18q genes or of X chromosome genes that are translocated next to the HCNEs. Alteration of expression of developmentally regulated genes by translocation of HCNEs may represent a common mechanism for disorders associated to isolated chromosomal rearrangements.     

Identification of Genomic Alterations in Pancreatic Cancer Using Array-Based Comparative Genomic Hybridization

Background

Genomic aberration is a common feature of human cancers and also is one of the basic mechanisms that lead to overexpression of oncogenes and underexpression of tumor suppressor genes. Our study aims to identify frequent genomic changes in pancreatic cancer.

Materials and Methods

We used array comparative genomic hybridization (array CGH) to identify recurrent genomic alterations and validated the protein expression of selected genes by immunohistochemistry.

Results

Sixteen gains and thirty-two losses occurred in more than 30% and 60% of the tumors, respectively. High-level amplifications at 7q21.3–q22.1 and 19q13.2 and homozygous deletions at 1p33–p32.3, 1p22.1, 1q22, 3q27.2, 6p22.3, 6p21.31, 12q13.2, 17p13.2, 17q21.31 and 22q13.1 were identified. Especially, amplification of AKT2 was detected in two carcinomas and homozygous deletion of CDKN2C in other two cases. In 15 independent validation samples, we found that AKT2 (19q13.2) and MCM7 (7q22.1) were amplified in 6 and 9 cases, and CAMTA2 (17p13.2) and PFN1 (17p13.2) were homozygously deleted in 3 and 1 cases. AKT2 and MCM7 were overexpressed, and CAMTA2 and PFN1 were underexpressed in pancreatic cancer tissues than in morphologically normal operative margin tissues. Both GISTIC and Genomic Workbench software identified 22q13.1 containing APOBEC3A and APOBEC3B as the only homozygous deletion region. And the expression levels of APOBEC3A and APOBEC3B were significantly lower in tumor tissues than in morphologically normal operative margin tissues. Further validation showed that overexpression of PSCA was significantly associated with lymph node metastasis, and overexpression of HMGA2 was significantly associated with invasive depth of pancreatic cancer.

Conclusion

These recurrent genomic changes may be useful for revealing the mechanism of pancreatic carcinogenesis and providing candidate biomarkers.     

Radiogenomic mapping of edema/cellular invasion MRI-phenotypes in glioblastoma multiforme

Background

Despite recent discoveries of new molecular targets and pathways, the search for an effective therapy for Glioblastoma Multiforme (GBM) continues. A newly emerged field, radiogenomics, links gene expression profiles with MRI phenotypes. MRI-FLAIR is a noninvasive diagnostic modality and was previously found to correlate with cellular invasion in GBM. Thus, our radiogenomic screen has the potential to reveal novel molecular determinants of invasion. Here, we present the first comprehensive radiogenomic analysis using quantitative MRI volumetrics and large-scale gene- and microRNA expression profiling in GBM.

Methods

Based on The Cancer Genome Atlas (TCGA), discovery and validation sets with gene, microRNA, and quantitative MR-imaging data were created. Top concordant genes and microRNAs correlated with high FLAIR volumes from both sets were further characterized by Kaplan Meier survival statistics, microRNA-gene correlation analyses, and GBM molecular subtype-specific distribution.

Results

The top upregulated gene in both the discovery (4 fold) and validation (11 fold) sets was PERIOSTIN (POSTN). The top downregulated microRNA in both sets was miR-219, which is predicted to bind to POSTN. Kaplan Meier analysis demonstrated that above median expression of POSTN resulted in significantly decreased survival and shorter time to disease progression (P<0.001). High POSTN and low miR-219 expression were significantly associated with the mesenchymal GBM subtype (P<0.0001).

Conclusion

Here, we propose a novel diagnostic method to screen for molecular cancer subtypes and genomic correlates of cellular invasion. Our findings also have potential therapeutic significance since successful molecular inhibition of invasion will improve therapy and patient survival in GBM.     

The Degree of Segmental Aneuploidy Measured by Total Copy Number Abnormalities Predicts Survival and Recurrence in Superficial Gastroesophageal Adenocarcinoma

Background

Prognostic biomarkers are needed for superficial gastroesophageal adenocarcinoma (EAC) to predict clinical outcomes and select therapy. Although recurrent mutations have been characterized in EAC, little is known about their clinical and prognostic significance. Aneuploidy is predictive of clinical outcome in many malignancies but has not been evaluated in superficial EAC.

Methods

We quantified copy number changes in 41 superficial EAC using Affymetrix SNP 6.0 arrays. We identified recurrent chromosomal gains and losses and calculated the total copy number abnormality (CNA) count for each tumor as a measure of aneuploidy. We correlated CNA count with overall survival and time to first recurrence in univariate and multivariate analyses.

Results

Recurrent segmental gains and losses involved multiple genes, including: HER2, EGFR, MET, CDK6, KRAS (recurrent gains); and FHIT, WWOX, CDKN2A/B, SMAD4, RUNX1 (recurrent losses). There was a 40-fold variation in CNA count across all cases. Tumors with the lowest and highest quartile CNA count had significantly better overall survival (p = 0.032) and time to first recurrence (p = 0.010) compared to those with intermediate CNA counts. These associations persisted when controlling for other prognostic variables.

Significance

SNP arrays facilitate the assessment of recurrent chromosomal gain and loss and allow high resolution, quantitative assessment of segmental aneuploidy (total CNA count). The non-monotonic association of segmental aneuploidy with survival has been described in other tumors. The degree of aneuploidy is a promising prognostic biomarker in a potentially curable form of EAC.     

Mapping of genetic abnormalities of primary tumours from metastatic CRC by high-resolution SNP arrays

Background

For years, the genetics of metastatic colorectal cancer (CRC) have been studied using a variety of techniques. However, most of the approaches employed so far have a relatively limited resolution which hampers detailed characterization of the common recurrent chromosomal breakpoints as well as the identification of small regions carrying genetic changes and the genes involved in them.

Methodology/Principal Findings

Here we applied 500K SNP arrays to map the most common chromosomal lesions present at diagnosis in a series of 23 primary tumours from sporadic CRC patients who had developed liver metastasis. Overall our results confirm that the genetic profile of metastatic CRC is defined by imbalanced gains of chromosomes 7, 8q, 11q, 13q, 20q and X together with losses of the 1p, 8p, 17p and 18q chromosome regions. In addition, SNP-array studies allowed the identification of small (<1.3 Mb) and extensive/large (>1.5 Mb) altered DNA sequences, many of which contain cancer genes known to be involved in CRC and the metastatic process. Detailed characterization of the breakpoint regions for the altered chromosomes showed four recurrent breakpoints at chromosomes 1p12, 8p12, 17p11.2 and 20p12.1; interestingly, the most frequently observed recurrent chromosomal breakpoint was localized at 17p11.2 and systematically targeted the FAM27L gene, whose role in CRC deserves further investigations.

Conclusions/Significance

In summary, in the present study we provide a detailed map of the genetic abnormalities of primary tumours from metastatic CRC patients, which confirm and extend on previous observations as regards the identification of genes potentially involved in development of CRC and the metastatic process.     

The notch ligand delta-like 4 regulates multiple stages of early hemato-vascular development

Background

In mouse embryos, homozygous or heterozygous deletions of the gene encoding the Notch ligand Dll4 result in early embryonic death due to major defects in endothelial remodeling in the yolk sac and embryo. Considering the close developmental relationship between endothelial and hematopoietic cell lineages, which share a common mesoderm-derived precursor, the hemangioblast, and many key regulatory molecules, we investigated whether Dll4 is also involved in the regulation of early embryonic hematopoiesis.

Methodology/Principal Findings

Using Embryoid Bodies (EBs) derived from embryonic stem cells harboring hetero- or homozygous Dll4 deletions, we observed that EBs from both genotypes exhibit an abnormal endothelial remodeling in the vascular sprouts that arise late during EB differentiation, indicating that this in vitro system recapitulates the angiogenic phenotype of Dll4 mutant embryos. However, analysis of EB development at early time points revealed that the absence of Dll4 delays the emergence of mesoderm and severely reduces the number of blast-colony forming cells (BL-CFCs), the in vitro counterpart of the hemangioblast, and of endothelial cells. Analysis of colony forming units (CFU) in EBs and yolk sacs from Dll4^+/− and Dll4^−/− embryos, showed that primitive erythropoiesis is specifically affected by Dll4 insufficiency. In Dll4 mutant EBs, smooth muscle cells (SMCs) were seemingly unaffected and cardiomyocyte differentiation was increased, indicating that SMC specification is Dll4-independent while a normal dose of this Notch ligand is essential for the quantitative regulation of cardiomyogenesis.

Conclusions/Significance

This study highlights a previously unnoticed role for Dll4 in the quantitative regulation of early hemato-vascular precursors, further indicating that it is also involved on the timely emergence of mesoderm in early embryogenesis.     

Relocation of genes generates non-conserved chromosomal segments in Fusarium graminearum that show distinct and co-regulated gene expression patterns

Background

Genome comparisons between closely related species often show non-conserved regions across chromosomes. Some of them are located in specific regions of chromosomes and some are even confined to one or more entire chromosomes. The origin and biological relevance of these non-conserved regions are still largely unknown. Here we used the genome of Fusarium graminearum to elucidate the significance of non-conserved regions.

Results

The genome of F. graminearum harbours thirteen non-conserved regions dispersed over all of the four chromosomes. Using RNA-Seq data from the mycelium of F. graminearum, we found weakly expressed regions on all of the four chromosomes that exactly matched with non-conserved regions. Comparison of gene expression between two different developmental stages (conidia and mycelium) showed that the expression of genes in conserved regions is stable, while gene expression in non-conserved regions is much more influenced by developmental stage. In addition, genes involved in the production of secondary metabolites and secreted proteins are enriched in non-conserved regions, suggesting that these regions could also be important for adaptations to new environments, including adaptation to new hosts. Finally, we found evidence that non-conserved regions are generated by sequestration of genes from multiple locations. Gene relocations may lead to clustering of genes with similar expression patterns or similar biological functions, which was clearly exemplified by the PKS2 gene cluster.

Conclusions

Our results showed that chromosomes can be functionally divided into conserved and non-conserved regions, and both could have specific and distinct roles in genome evolution and regulation of gene expression.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-191) contains supplementary material, which is available to authorized users.     

Gene expression profiling and molecular characterization of antimony resistance in Leishmania amazonensis

Background

Drug resistance is a major problem in leishmaniasis chemotherapy. RNA expression profiling using DNA microarrays is a suitable approach to study simultaneous events leading to a drug-resistance phenotype. Genomic analysis has been performed primarily with Old World Leishmania species and here we investigate molecular alterations in antimony resistance in the New World species L. amazonensis.

Methods/Principal Findings

We selected populations of L. amazonensis promastigotes for resistance to antimony by step-wise drug pressure. Gene expression of highly resistant mutants was studied using DNA microarrays. RNA expression profiling of antimony-resistant L. amazonensis revealed the overexpression of genes involved in drug resistance including the ABC transporter MRPA and several genes related to thiol metabolism. The MRPA overexpression was validated by quantitative real-time RT-PCR and further analysis revealed that this increased expression was correlated to gene amplification as part of extrachromosomal linear amplicons in some mutants and as part of supernumerary chromosomes in other mutants. The expression of several other genes encoding hypothetical proteins but also nucleobase and glucose transporter encoding genes were found to be modulated.

Conclusions/Significance

Mechanisms classically found in Old World antimony resistant Leishmania were also highlighted in New World antimony-resistant L. amazonensis. These studies were useful to the identification of resistance molecular markers.     

A Genotype-First Approach for the Molecular and Clinical Characterization of Uncommon De Novo Microdeletion of 20q13.33

Background

Subtelomeric deletions of the long arm of chromosome 20 are rare, with only 11 described in the literature. Clinical features of individuals with these microdeletions include severe limb malformations, skeletal abnormalities, growth retardation, developmental and speech delay, mental retardation, seizures and mild, non-specific dysmorphic features.

Methodology/Principal Findings

We characterized microdeletions at 20q13.33 in six individuals referred for genetic evaluation of developmental delay, mental retardation, and/or congenital anomalies. A comparison to previously reported cases of 20q13.33 microdeletion shows phenotypic overlap, with clinical features that include mental retardation, developmental delay, speech and language deficits, seizures, and behavior problems such as autistic spectrum disorder. There does not appear to be a clinically recognizable constellation of dysmorphic features among individuals with subtelomeric 20q microdeletions.

Conclusions/Significance

Based on genotype-phenotype correlation among individuals in this and previous studies, we discuss several possible candidate genes for specific clinical features, including ARFGAP1, CHRNA4 and KCNQ2 and neurodevelopmental deficits. Deletion of this region may play an important role in cognitive development.     

Mapping Insertions,Deletions and SNPs on Venter's Chromosomes

Background

The very recent availability of fully sequenced individual human genomes is a major revolution in biology which is certainly going to provide new insights into genetic diseases and genomic rearrangements.

Results

We mapped the insertions, deletions and SNPs (single nucleotide polymorphisms) that are present in Craig Venter''s genome, more precisely on chromosomes 17 to 22, and compared them with the human reference genome hg17. Our results show that insertions and deletions are almost absent in L1 and generally scarce in L2 isochore families (GC-poor L1+L2 isochores represent slightly over half of the human genome), whereas they increase in GC-rich isochores, largely paralleling the densities of genes, retroviral integrations and Alu sequences. The distributions of insertions/deletions are in striking contrast with those of SNPs which exhibit almost the same density across all isochore families with, however, a trend for lower concentrations in gene-rich regions.

Conclusions

Our study strongly suggests that the distribution of insertions/deletions is due to the structure of chromatin which is mostly open in gene-rich, GC-rich isochores, and largely closed in gene-poor, GC-poor isochores. The different distributions of insertions/deletions and SNPs are clearly related to the two different responsible mechanisms, namely recombination and point mutations.     

Genome-wide linkage scan of a pedigree with familial hypercholesterolemia suggests susceptibility loci on chromosomes 3q25-26 and 21q22

Background

Familial hypercholesterolemia (FH) is a heritable disorder that can increase the risk of premature coronary heart disease. Studies suggest there are substantial genetic heterogeneities for different populations. Here we tried to identify novel susceptibility loci for FH in a Chinese pedigree.

Methodology/Principal Findings

We performed a SNP-based genome-wide linkage scan with the Chinese FH pedigree. Two suggestive linkage loci not previously reported were identified on chromosomes 3q25.1-26.1 (NPL = 9.01, nominal P<0.00001, and simulated occurrence per genome scan = 1.08) and 21q22.3 (NPL = 8.95, nominal P<0.00001, and simulated occurrence per genome scan = 1.26). In the interaction analysis with a trimmed version of the pedigree, we obtained a significantly increased joint LOD score (2.70) compared with that obtained when assuming the two loci uncorrelated, suggesting that more than one locus was involved in this pedigree. Exon screening of two candidate genes ABCG1 and LSS from one of the suggestive region 21q22 didn''t report any causative mutations.

Conclusions/Significances

These results confirm complex etiologies and suggest new genetic casual factors for the FH disorder. Further study of the two candidate regions is advocated.     

Association between acquired uniparental disomy and homozygous mutations and HER2/ER/PR status in breast cancer

Background

Genetic alterations in cellular signaling networks are a hallmark of cancer, however, effective methods to discover them are lacking. A novel form of abnormality called acquired uniparental disomy (aUPD) was recently found to pinpoint the region of mutated genes in various cancers, thereby identifying the region for next-generation sequencing.

Methods/Principal Findings

We retrieved large genomic data sets from the Gene Expression Omnibus database to perform genome-wide analysis of aUPD in breast tumor samples and cell lines using approaches that can reliably detect aUPD. aUPD was identified in 52.29% of the tumor samples. The most frequent aUPD regions were located at chromosomes 2q, 3p, 5q, 9p, 9q, 10q, 11q, 13q, 14q and 17q. We evaluated the data for any correlation between the most frequent aUPD regions and HER2/neu, ER, and PR status, and found a statistically significant correlation between the recurrent regions of aUPD and triple negative (TN) breast cancers. aUPD at chromosome 17q (VEZF1, WNT3), 3p (SUMF1, GRM7), 9p (MTAP, NFIB) and 11q (CASP1, CASP4, CASP5) are predictors for TN. The frequency of aUPD was found to be significantly higher in TN breast cancer cases compared to HER2/neu-positive and/or ER or PR-positive cases. Furthermore, using previously published mutation data, we found TP53 homozygously mutated in cell lines having aUPD in that locus.

Conclusions/Significance

We conclude that aUPD is a common and non-random molecular feature of breast cancer that is most prominent in triple negative cases. As aUPD regions are different among the main pathological subtypes, specific aUPD regions may aid the sub-classification of breast cancer. In addition, we provide statistical support using TP53 as an example that identifying aUPD regions can be an effective approach in finding aberrant genes. We thus conclude that a genome-wide scale analysis of aUPD regions for homozygous sequence alterations can provide valuable insights into breast tumorigenesis.     

Molecular characterization of acquired tolerance of tumor cells to picropodophyllin (PPP)

Background

Picropodophyllin (PPP) is a promising novel anti-neoplastic agent that efficiently kills tumor cells in vitro and causes tumor regression and increased survival in vivo. We have previously reported that PPP treatment induced moderate tolerance in two out of 10 cell lines only, and here report the acquired genomic and expression alterations associated with PPP selection over 1.5 years of treatment.

Methodology/Principal Findings

Copy number alterations monitored using metaphase and array-based comparative genomic hybridization analyses revealed largely overlapping alterations in parental and maximally tolerant cells. Gain/ amplification of the MYC and PVT1 loci in 8q24.21 were verified on the chromosome level. Abnormalities observed in connection to PPP treatment included regular gains and losses, as well as homozygous losses in 10q24.1-q24.2 and 12p12.3-p13.2 in one of the lines and amplification at 5q11.2 in the other. Abnormalities observed in both tolerant derivatives include amplification/gain of 5q11.2, gain of 11q12.1-q14.3 and gain of 13q33.3-qter. Using Nexus software analysis we combined the array-CGH data with data from gene expression profilings and identified genes that were altered in both inputs. A subset of genes identified as downregulated (ALDH1A3, ANXA1, TLR4 and RAB5A) or upregulated (COX6A1, NFIX, ME1, MAPK and TAP2) were validated by siRNA in the tolerant or parental cells to alter sensitivity to PPP and confirmed to alter sensitivity to PPP in further cell lines.

Conclusions

Long-term PPP selection lead to altered gene expression in PPP tolerant cells with increase as well as decrease of genes involved in cell death such as PTEN and BCL2. In addition, acquired genomic copy number alterations were observed that were often reflected by altered mRNA expression levels for genes in the same regions.     

OPCML is a broad tumor suppressor for multiple carcinomas and lymphomas with frequently epigenetic inactivation

Background

Identification of tumor suppressor genes (TSGs) silenced by CpG methylation uncovers the molecular mechanism of tumorigenesis and potential tumor biomarkers. Loss of heterozygosity at 11q25 is common in multiple tumors including nasopharyngeal carcinoma (NPC). OPCML, located at 11q25, is one of the downregulated genes we identified through digital expression subtraction.

Methodology/Principal Findings

Semi-quantitative RT-PCR showed frequent OPCML silencing in NPC and other common tumors, with no homozygous deletion detected by multiplex differential DNA-PCR. Instead, promoter methylation of OPCML was frequently detected in multiple carcinoma cell lines (nasopharyngeal, esophageal, lung, gastric, colon, liver, breast, cervix, prostate), lymphoma cell lines (non-Hodgkin and Hodgkin lymphoma, nasal NK/T-cell lymphoma) and primary tumors, but not in any non-tumor cell line and seldom weakly methylated in normal epithelial tissues. Pharmacological and genetic demethylation restored OPCML expression, indicating a direct epigenetic silencing. We further found that OPCML is stress-responsive, but this response is epigenetically impaired when its promoter becomes methylated. Ecotopic expression of OPCML led to significant inhibition of both anchorage-dependent and -independent growth of carcinoma cells with endogenous silencing.

Conclusions/Significance

Thus, through functional epigenetics, we identified OPCML as a broad tumor suppressor, which is frequently inactivated by methylation in multiple malignancies.