Loss, Mutation and Deregulation of L3MBTL4 in Breast Cancers
Total Page:16
File Type:pdf, Size:1020Kb
Addou-Klouche et al. Molecular Cancer 2010, 9:213 http://www.molecular-cancer.com/content/9/1/213 RESEARCH Open Access Loss, mutation and deregulation of L3MBTL4 in breast cancers Lynda Addou-Klouche1, José Adélaïde1, Pascal Finetti1, Nathalie Cervera1, Anthony Ferrari1, Ismahane Bekhouche1, Fabrice Sircoulomb1, Christos Sotiriou2, Patrice Viens3,4, Soraya Moulessehoul5, François Bertucci1,3,4, Daniel Birnbaum1, Max Chaffanet1* Abstract Background: Many alterations are involved in mammary oncogenesis, including amplifications of oncogenes and losses of tumor suppressor genes (TSG). Losses may affect almost all chromosome arms and many TSGs remain to be identified. Results: We studied 307 primary breast tumors and 47 breast cancer cell lines by high resolution array comparative genomic hybridization (aCGH). We identified a region on 18p11.31 lost in about 20% of the tumors and 40% of the cell lines. The minimal common region of loss (Chr18:6,366,938-6,375,929 bp) targeted the L3MBTL4 gene. This gene was also targeted by breakage in one tumor and in two cell lines. We studied the exon sequence of L3MBTL4 in 180 primary tumor samples and 47 cell lines and found six missense and one nonsense heterozygous mutations. Compared with normal breast tissue, L3MBTL4 mRNA expression was downregulated in 73% of the tumors notably in luminal, ERBB2 and normal-like subtypes. Losses of the 18p11 region were associated with low L3MBTL4 expression level. Integrated analysis combining genome and gene expression profiles of the same tumors pointed to 14 other potential 18p TSG candidates. Downregulated expression of ZFP161, PPP4R1 and YES1 was correlated with luminal B molecular subtype. Low ZFP161 gene expression was associated with adverse clinical outcome. Conclusion: We have identified L3MBTL4 as a potential TSG of chromosome arm 18p. The gene is targeted by deletion, breakage and mutations and its mRNA is downregulated in breast tumors. Additional 18p TSG candidates might explain the aggressive phenotype associated with the loss of 18p in breast tumors. Background and 40% of the cell lines. We delineated a minimal com- The development and progression of breast cancer is the monregionofdeletionthattargetedtheL3MBTL4 result of the accumulation of genetic alterations such as gene. The L3MBTL4 gene is one of the four human amplification of oncogenes and deletions of tumor sup- orthologs of Drosophila lethal (3) malignant brain tumor pressor genes (TSG) in the epithelial cells of the mam- (l(3)mbt). L(3)mbt is a bona fide TSG in the fly [4,5]. mary gland. Frequent deletions have been reported on The L3MBTL4 protein contains three “malignant brain chromosome arms 1p, 3p, 7q, 8p 9p, 16q and 17p but tumor” (MBT) domains. This domain of about 100 only few reports describe such deletions on chromo- amino acid residues is conserved in protostomians and some arm 18p [1-3]. deuterostomiansandoftenexistsasrepeats[6].The We profiled a series of 307 primary breast tumors and MBT domain binds methylated histone residues. The 47 breast cancer cell lines by using high resolution array human genome contains several MBT-containing pro- comparative genomic hybridization (aCGH). We identi- teins, some of which have been linked to gene regula- fied a region on 18p11.31 deleted in 25% of the tumors tory pathways and polycomb-mediated repression, and to cancer [7]. L3MBTL1 is a TSG implicated in myeloid * Correspondence: [email protected] malignancies [8], and L3MBTL3 deregulation is asso- 1 Marseille Cancer Research Center; Department of Molecular Oncology, ciated with neuroblastoma [9]. UMR891 Inserm; Institut Paoli-Calmettes; Marseille, France Full list of author information is available at the end of the article © 2010 Addou-Klouche et al; licensee BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. Addou-Klouche et al. Molecular Cancer 2010, 9:213 Page 2 of 13 http://www.molecular-cancer.com/content/9/1/213 To document the involvement of L3MBTL4 in breast extracted from a normal area of paraffin-embedded cancer we searched for mutations by sequence analysis T8584, T8847, T9193 and T8525 tissues. of 180 primary tumor samples and 47 cell lines. We found that, in addition to deletions, L3MBTL4 is tar- Genome analysis by array-comparative genomic geted by mutations. Finally, we found that L3MBTL4 hybridization (aCGH) deletions correlate with low mRNA expression and with Genome profiles were established on 307 tumors and 47 the presence of lymph node metastasis, high Scarf- breast cancer cell lines by aCGH using 244K CGH Bloom-Richardson (SBR) grade and luminal B molecular Microarrays (Hu-244A, Agilent Technologies, Massy, subtype. France) as previously described [16]. A pool of 13 nor- Our study is the first to identify a region on chromo- mal male DNAs was used as reference. Scanning was some arm 18p likely to contain a putative TSG involved done with Agilent Autofocus Dynamic Scanner in breast oncogenesis, and to show that this TSG may (G2565BA, Agilent Technologies). Data analysis was be L3MBTL4. We did not exclude the existence of other done and visualized with CGH Analytics 3.4 software potential 18p TSGs, some of which could explain the (Agilent Technologies). Extraction of data (log2 ratio) aggressive phenotype of breast tumors with 18p loss. was done from CGH analytics, while normalized and fil- tered log2 ratio were obtained from “Feature extraction” Samples and Methods software (Agilent Technologies). Data generated by Breast tumors probes mapped to × and Y chromosomes were elimi- Tumor tissues were collected from 307 patients with nated. The final dataset contained 225,388 unique primary adenocarcinoma who underwent initial surgery probes covering 22,509 genes and intergenic regions at the Institut Paoli-Calmettes (Marseille, France) according to the hg17/NCBI human genome mapping between 1992 and 2004. Immediately after macroscopic database (build 35). Data were analyzed using circular examination of the surgery specimen by two patholo- binary segmentation (CBS) [17] as implemented in the gists, samples containing more than 60% of tumor cells DNA copy R/Bioconductor package [18] with default were obtained, frozen in liquid nitrogen and stored at parameters to translate intensity measurements in -80°C until nucleic acids extraction. The main histoclini- regions of equal copy number, each region being defined cal characteristics of tumors are listed in the Additional by at least five consecutive probes. Thus, each probe file 1, Table S1. Each patient gave written informed con- was assigned a segment value referred to as its sent and the study was approved by our institutional “smoothed” value. review committee. We used a threshold value of |0.33| (log2 ratio) to define a copy number aberration (CNA) [16]. Identifica- Cell lines tion of copy number variations (CNV) was done using The 47 breast cell lines used in this study were BT-20, the regions published by McCarroll et al [19] which are BT-483, HCC1937, CAMA-1, HCC38, HCC1500 [10], stored in the Database of Genomic Variants (release v8) Hs 578T [11], BT-474, BT-549, HCC202, HCC1395, [20]. To determine altered regions, we used the GISTIC HCC1569, HCC1806, HCC1954, HCC2218, HME-1, algorithm [21], which computes for each genomic seg- MCF7, MDA-MB-134, MDA-MB-157, MDA-MB-175, ment through the whole genome a score based on the MDA-MB-231, MDA-MB-361, MDA-MB-415, MDA- frequency of CNA combined with its amplitude, with MB-436, MDA-MB-453, MDA-MB-468, SK-BR-3, bootstrapping to calculate the significance level (p < SK-BR-7, T47 D, UACC-812, UACC-893, ZR-75-1, ZR- 0.001). To establish significant association between CNA 75-30 (ATTC, Manassas, VA), BrCa-Mz-01, BrCa-Mz- and categorical variables, Fisher’sexacttestwasused. 02 [12], SUM44, SUM52, SUM102, SUM149, SUM159, Gains and losses were handled separately. SUM185, SUM190, SUM206, SUM225, SUM229 [13], S68 (V. Catros, Rennes, France), and CAL51 [14]. All Mutation analysis cell lines were grown as recommended by their Sequence analysis was done on 180 out of the 307 suppliers. breast tumors and on the 47 breast cell lines after amplification of genomic DNA. Seventeen primers pairs DNA and RNA extraction were designed to amplify L3MBTL4 DNA by polymerase DNA and RNA were extracted from frozen samples by chain reaction (PCR) (Additional file 1, Table S2). PCR using guanidium isothiocynanate and cesium chloride amplifications were done in a total volume of 25 μl PCR gradient, as previously described [15]. DNA quality and mix containing at least 10 ng template DNA, Taq buffer, RNA integrity were respectively controled on polyacryla- 500 μmol of each deoxynucleotide triphosphate, 100 mide gel electrophoresis and on Agilent Bioanalyzer μmol of each primer and 1 unit of Hot Star Taq DNA (AgilentTechnologies,Massy,France).DNAwasalso polymerase (Qiagen, Courtaboeuf, France). PCR Addou-Klouche et al. Molecular Cancer 2010, 9:213 Page 3 of 13 http://www.molecular-cancer.com/content/9/1/213 amplification conditions were as follows: 95°C 10 min All probes were mapped based on their EntrezGeneID. and 95°C 30 sec, annealing T° 55°C 30 sec, 72°C 45 sec When multiple probes were mapped to the same gene, to to 1 min depending on PCR product length, 72°C 10 the probe sets with an extension « at », next « s_at », min, 35 cycles. PCR products were purified using Milli- and followed by all other extensions were preferentially pore plate MSNU030 (Millipore SAS, Molsheim, kept. When several probe sets with the best extension France). Aliquots (1 μl) of the purified PCR products were available, the one with the highest median value were sequenced using Big Dye terminator v1.1 Cycling was retained. Sequencing Kit (Applied Biosystems, Courtaboeuf, The five molecular subtypes related to the intrinsic France) including the forward or reverse primer.