A 5'-region polymorphism modulates promoter activity of the tumor suppressor gene MFSD2A
- Francesca Colombo†1,
- Felicia S Falvella†1,
- Antonella Galvan1,
- Elisa Frullanti1,
- Hideo Kunitoh2,
- Toshikazu Ushijima2 and
- Tommaso A Dragani1Email author
© Colombo et al; licensee BioMed Central Ltd. 2011
Received: 29 September 2010
Accepted: 7 July 2011
Published: 7 July 2011
The MFSD2A gene maps within a linkage disequilibrium block containing the MYCL1-EcoRI polymorphism associated with prognosis and survival in lung cancer patients. Survival discrepancies between Asians and Caucasians point to ethnic differences in allelic frequencies of the functional genetic variations.
Analysis of three single-nucleotide polymorphisms (SNPs) mapping in the MFSD2A 5'-regulatory region using a luciferase reporter system showed that SNP rs12072037, in linkage disequilibrium with the MYCL1-EcoRI polymorphism and polymorphic in Asians but not in Caucasians, modulated transcriptional activity of the MFSD2A promoter in cell lines expressing AHR and ARNT transcription factors, which potentially bind to the SNP site.
SNP rs12072037 modulates MFSD2A promoter activity and thus might affect MFSD2A levels in normal lung and in lung tumors, representing a candidate ethnically specific genetic factor underlying the association between the MYCL1 locus and lung cancer patients' survival.
A 106-kb linkage disequilibrium (LD) block on chromosome 1p34, which includes the TRIT1, MYCL1, and MFSD2A genes, is associated with lung cancer prognosis and survival , although conflicting results of the association of this region with prognosis, in particular of the MYCL1-Eco RI polymorphism, have been reported [2, 3]. Indeed, association of MYCL1-EcoRI with lung cancer patients' survival was observed in all of 4 studies of Asians, but in none of 3 studies on Caucasians . The discrepancies might reflect ethnic differences in allelic frequencies of the functional genetic variants mapping in this locus, as suggested by the significant difference between Caucasian and Asian subjects in the frequencies of several SNPs located in the TRIT1, MYCL1, and MFSD2A gene .
Modulation of expression of a gene mapping in the MYCL1 region may represent a mechanism underlying the association of this region with cancer patients' survival. Indeed, MYCL1 expression is not detected in normal or tumor tissue. Both the TRIT1 and MFSD2A genes are downregulated in lung adenocarcinomas (ADCA), whereas overexpression of either gene has tumor-suppressor effects [1, 4, 5]. The MFSD2A gene was also strongly downregulated in a panel of non-small cell lung cancer (NSCLC) cell lines, where it inhibits cell adhesion and migration when overexpressed . Thus, available data suggest that downregulation of MFSD2A plays a role in lung tumor progression.
Since functional polymorphisms in the promoter region may affect mRNA levels of target genes by altering transcription factor (TF) binding sites [6, 7], we analyzed three single-nucleotide polymorphisms (SNPs) (rs3131703, rs12072037, and rs3738668) mapping in the MFSD2A 5' region for a potential role in altering MFSD2A promoter activity.
SNP rs3131703 in the MFSD2A 5' regulatory region has no functional effects on transcriptional activity
Allele frequency of the three validated SNPs in MFSD2A 5' region by ethnicity
Chromosome position a
Frequency of the minor frequency allele (MAF) in HapMap populations a
Analysis by quantitative real-time (qRT)-PCR to test for an association between genotype and MFSD2A mRNA levels revealed no significant difference between two genotype groups of normal lung tissue samples from 20 subjects in our series selected for homozygosity at either allele (10 GG versus 10 AA samples) (data not shown).
SNP rs12072037 is the major polymorphism involved in modulation of MFSD2A promoter activity
The MFSD2A 5' region rs12072037 (C/A alleles) and rs3738668 (C/A alleles) SNPs, mapping at -759 bp and -246 bp from ATG, respectively, are not polymorphic in Caucasians, whereas in Asian populations, the frequencies of the minor allele of these polymorphisms are quite high (MAF = 0.47-0.48; Table 1). Accordingly, sequencing analysis in a small number of Italian (n = 15) and Japanese (n = 15) lung cancer patients showed that SNPs rs12072037 and rs3738668 defined two haplotypes in Asians, AA and CC, whereas only the CC haplotype was detected in Caucasians.
Similar analyses using a 699-bp fragment containing only SNP rs3738668 (A and C alleles; Figure 1C) revealed no overall modulation of MFSD2A promoter activity, as indicated by luciferase values, except for a marginally statistically significant effect only in Hek293T cells (P = 0.022). As reported for the other constructs, a > 100-fold increase in luciferase activity in Hek293T cells was observed in the presence of either alleles of rs3738668 SNP.
These findings point to SNP rs12072037 as the main modulator of MFSD2A promoter activity.
SNP rs12072037 alleles create different putative transcription factor binding sites
To test the binding of AHR/ARNT heterodimer to the construct containing the A allele, a preliminary chromatin immunoprecipitation (ChIP) assay was carried out using chromatin from cells transfected with either of the two constructs. Real-time PCR analysis revealed an approximately 2-fold increase of immunoprecipitated DNA in the presence of the Asian allele (not shown).
We investigated the possible functional role of three SNPs, located in the MFSD2A 5'-regulatory region and showing population-specific allele frequencies, in modulating MFSD2A promoter activity. Luciferase reporter assay indicated no such modulation by SNP rs3131703, polymorphic only in the Caucasian population, and MFSD2A mRNA levels in normal lung tissue of lung cancer patients showed no statistically significant association with rs3131703, thus excluding the candidacy of this SNP as a functional genetic element regulating MFSD2A promoter activity. By contrast, analyses of SNPs rs12072037 and rs3738668, polymorphic in Asians, pointed to a functional role for rs12072037 in modulating MFSD2A transcriptional activity, showing allele-specific promoter activity in 3 of 7 cell lines.
Since promoter activity depends on the TF expression profile , a SNP in a regulatory region could modulate gene expression by altering TF binding sites, eliminating the natural site or generating a novel one and thus altering binding strength .
To investigate the basis of the observed cell type-specific differential effects on transcriptional activity, we measured expression of in silico- predicted TFs in the cell lines used for promoter activity analysis. Transcriptional levels of both AHR and ARNT were high in the same three cell lines (HepG2, HT-29, and IGROV1) showing higher MFSD2A promoter activity associated with the A allele of rs12072037.
Accordingly, AHR and ARNT specifically bind to the MFSD2A promoter containing the A allele of the rs12072037 polymorphism. Moreover, NCI-H520 and NCI-H596 cells, in which no modulation of MFSD2A promoter activity by genetic polymorphisms was detected, showed high levels of ARNT, but not of AHR, suggesting that the AHR/ARNT heterodimer rather than the single elements is necessary to increase MFSD2A promoter activity, as reported for other target genes . Thus, the cell type-specific modulation of reporter mRNA levels by the MFSD2A promoter might be associated to the joint expression of AHR and ARNT, which may form a dimer upregulating the MFSD2A promoter variant containing the A allele at rs12072037. Our preliminary ChIP experiment supported this hypothesis since differences in immunoprecipitated DNA levels were measured between the A and C alleles of rs12072037. However, transfection of recombinant plasmids in cancer cells where the endogenous MFSD2A gene is strongly downregulated may not be adequate to test the TFs modulating the natural MFSD2A promoter in normal cells. To this aim, comparison of AHR and ARNT binding to the MFSD2A promoter between normal cells carrying the different rs12072037 genotypes would be required. Also, the MFSD2A promoter fragment that we have transfected contains an additional AHR/ARNT binding site downstream to the assayed site, causing a potential technical bias in the ChIP assay. Therefore, we cannot exclude a role for other TFs in affecting MFSD2A promoter activity.
SNP rs12072037, which modulated MFSD2A transcriptional activity, showed the most statistically significant differences of allele frequencies between Caucasian and Asian subjects among 12 SNPs mapping in the LD block ; indeed, SNP rs12072037 is almost invariant in Caucasians (Table 1). Furthermore, in Asians, the rs12072037 functional promoter polymorphism shows a significant LD with the MYCL1 SNP rs3134613 (EcoRI) (D' = 0.75, r^2 = 0.48, n = 88; genotypes of JPT samples downloaded from HapMap, accessed on August 24, 2009; analysis carried out using the JLIN program ), a SNP reported to be associated with prognostic factors and survival of lung cancer patients in Asian patients [3, 13, 14]. Together, these data suggest the candidacy of rs12072037 as the functional variation in the MYCL1 locus responsible for modulating nodal status, metastasis occurrence, and survival of lung cancer patients of Asian ethnicity.
We identified SNP rs12072037 as the major polymorphism modulating MFSD2A promoter activity. Further studies in Asian lung cancer patients are needed to test the association of rs12072037 genotypes with MFSD2A mRNA levels in normal lung tissue and to clarify the role of this polymorphism in lung cancer progression and prognosis.
Materials and Methods
Human lung cancer (A549, NCI-H520, NCI-H596), hepatocellular carcinoma (HepG2), colorectal adenocarcinoma (HT-29), ovarian adenocarcinoma (IGROV1) and embryonal kidney (Hek293T) cells were used. All represent adherent epithelial cells. Each cell line was propagated in the appropriate culture medium as recommended by American Type Culture Collection, except for IGROV1 cells, which were cultured in RPMI 1640 medium with 2 mM L-glutamine, 1.5 g/l sodium bicarbonate, 10 mM HEPES, 1.0 mM sodium pyruvate, 4.5 g/l glucose and 10% FBS.
Search and genotyping of MFSD2A SNPs in 5'-region
Genomic DNAs were obtained from lung ADCA patients (n = 151) enrolled at Istituto Nazionale Tumori, Milan, Italy, and from lung cancer cases (n = 15) at the National Cancer Center, Tokyo, Japan. Patients gave written permission to use their biological material for research purposes, and study protocols were approved by the Institute committees for ethics.
SNP rs3131703 was genotyped in 151 Caucasian samples using pyrosequencing analysis on a PSQ96MA system (Biotage AB, Uppsala, Sweden), according to the manufacturer's instructions using specific primers reported in Additional file 1.
To identify haplotypes defined by SNPs rs12072037 and rs3738668, genomic DNAs of 15 Caucasian and 15 Asian lung cancer samples were PCR-amplified for 1160 bp spanning from the 5'-region to an initial part of intron 1-2 of MFSD2A using primer pairs reported in Additional file 1. Sequences were determined using an automatic sequencer (Applied Biosystems, Foster City, CA) and aligned and compared using Genomatix Dialign software (http://www.genomatix.de).
Quantitative real-time PCR (qRT-PCR)
Total RNAs were extracted with TRIZOL® Plus RNA Purification Kit (Invitrogen, Carlsbad, CA) from cell lines and 20 normal lung tissue samples excised during thoracic surgery for lung cancer at Istituto Nazionale Tumori, Milan, Italy. Total RNAs were reverse-transcribed using a 1:1 mix of oligo(dT)18 and random hexamer primers, according to the protocol recommended for the Transcriptor First-Strand cDNA Synthesis Kit (Roche). MFSD2A gene expression levels in lung samples were measured by qRT-PCR as described . Transcription factor mRNA levels in cell lines were assayed by qRT-PCR using primer pairs listed in Additional file 1. Amplification mixtures contained cDNA template diluted in RNase-free water, 12.5 μl 2 × Power SYBR® Green PCR Master Mix (Applied Biosystems) and 0.3 μM specific PCR primers in a final volume of 25 μl. The human hypoxanthine phosphoribosyltransferase 1 (HPRT1) gene was used as housekeeping control. Reactions were run in duplicate on the 7900HT System (Applied Biosystems). Each experiment was repeated three times. Relative changes in mRNA levels were assessed using the comparative cycle threshold (Ct) method, and relative quantities were calculated using the A549 sample as calibrator.
Subcloning of MFSD2A 5'-region allelic variants into luciferase reporter vector
Genomic DNAs of previously genotyped Caucasian individuals were selected for the two different alleles (G/A) of rs3131703 (Figure 1A) to be amplified for the 5' region (1499 bp from ATG). Genomic DNAs of previously genotyped Caucasian or Asian individuals were selected for the CC or AA haplotype (rs12072037 and rs3738668), respectively, to be amplified for 933 bp upstream of ATG (Figure 1B). Finally, the shorter fragments (699 bp from ATG), containing either of the two alleles (A/C) of rs3738668 (Figure 1C), were PCR-amplified from the previously subcloned longer fragments (Figure 1B). PCR amplifications were carried out using BIO-X-ACT Short DNA Polymerase kit (Bioline, London, UK) and primer pairs listed in Additional file 1. Fragments were cloned into pGL3-Basic vector (Promega, Madison, WI) upstream of the firefly luciferase reporter gene using Rapid DNA Dephos & Ligation Kit (Roche, Basel, Switzerland), according to the manufacturer's instructions.
Luciferase reporter assays
Seven human cell lines derived from different tissues were co-transfected with pGL3-Basic recombinant vectors containing the different fragments of the MFSD2A 5' region and the pRL-TK plasmid containing the Renilla gene (Promega), using FuGene HD transfection reagent (Roche) according to the manufacturer's instructions. Firefly luciferase was used as primary reporter to monitor mRNA levels, whereas Renilla luciferase served as control reporter for normalization. Cells were assayed for luciferase activity 48 h after transfection using the Dual-Luciferase Reporter Assay System kit (Promega), according to the manufacturer's instructions. Normalized signals were expressed as relative luciferase units. Data from at least 6 replicas were evaluated.
Chromatin immunoprecipitation (ChIP) assay
HT-29 cells (2 × 106 cells/immunoprecipitation) were transiently transfected with the pGL3-Basic recombinant vectors containing the intermediate length fragment (Figure 1B) of the MFSD2A 5' region and after 24 hours crosslinking with 1% formaldehyde was carried out. Cell lysis, micrococcal nuclease digestion of chromatin, immunoprecipitation (IP) and DNA recovery were carried out following the protocol of the Pierce Agarose ChIP kit (Pierce, Thermo Scientific, Rockford, IL). IP was carried out using anti-AHR and anti-ARNT transcription factors (MA1-513 and PA3-16508 respectively, Thermo Scientific) or normal rabbit IgG, as negative control. Recovered and purified DNA was quantified with Nanodrop Spectrophotometer ND-1000 (NanoDrop products, Wilmington, DE, USA) and quantitative real-time PCR was carried out with primers flanking the transcription factor binding site (Additional File 1). DNA levels were calculated by subtracting the Ct of the negative control sample from the Ct of the immunoprecipitated sample, to correct for the plasmidic DNA carry-over.
Associations of quantitative measurements with genotypes or MFSD2A promoter variants were assessed by analysis of variance, using the Rcmdr package in R .
- Spinola M, Falvella FS, Galvan A, Pignatiello C, Leoni VP, Pastorino U, Paroni R, Chen S, Skaug V, Haugen A, Dragani TA: Ethnic differences in frequencies of gene polymorphisms in the MYCL1 region and modulation of lung cancer patients' survival. Lung Cancer. 2007, 55: 271-277. 10.1016/j.lungcan.2006.10.023View ArticlePubMedGoogle Scholar
- Shih CM, Kuo YY, Wang YC, Jian SL, Hsu YT, Wu HY, Guo MW, Wang YC: Association of L-myc polymorphism with lung cancer susceptibility and prognosis in relation to age-selected controls and stratified cases. Lung Cancer. 2002, 36: 125-132. 10.1016/S0169-5002(01)00467-6View ArticlePubMedGoogle Scholar
- Spinola M, Pedotti P, Dragani TA, Taioli E: Meta-analysis suggests association of L-myc EcoRI polymorphism with cancer prognosis. Clin Cancer Res. 2004, 10: 4769-4775. 10.1158/1078-0432.CCR-04-0055View ArticlePubMedGoogle Scholar
- Spinola M, Galvan A, Pignatiello C, Conti B, Pastorino U, Nicander B, Paroni R, Dragani TA: Identification and functional characterization of the candidate tumor suppressor gene TRIT1 in human lung cancer. Oncogene. 2005, 24: 5502-5509. 10.1038/sj.onc.1208687View ArticlePubMedGoogle Scholar
- Spinola M, Falvella FS, Colombo F, Sullivan JP, Shames DS, Girard L, Spessotto P, Minna JD, Dragani TA: MFSD2A is a novel lung tumor suppressor gene modulating cell cycle and matrix attachment. Mol Cancer. 2010, 9: 62- 10.1186/1476-4598-9-62PubMed CentralView ArticlePubMedGoogle Scholar
- Mohrenweiser H: Survey of polymorphic sequence variation in the immediate 5' region of human DNA repair genes. Mutat Res. 2007, 616 (1-2): 221-226.View ArticlePubMedGoogle Scholar
- Guo Y, Jamison DC: The distribution of SNPs in human gene regulatory regions. BMC Genomics. 2005, 6: 140- 10.1186/1471-2164-6-140PubMed CentralView ArticlePubMedGoogle Scholar
- Cartharius K, Frech K, Grote K, Klocke B, Haltmeier M, Klingenhoff A, Frisch M, Bayerlein M, Werner T: MatInspector and beyond: promoter analysis based on transcription factor binding sites. Bioinformatics. 2005, 21 (13): 2933-2942. 10.1093/bioinformatics/bti473View ArticlePubMedGoogle Scholar
- Mottagui-Tabar S, Faghihi MA, Mizuno Y, Engstrom PG, Lenhard B, Wasserman WW, Wahlestedt C: Identification of functional SNPs in the 5-prime flanking sequences of human genes. BMC Genomics. 2005, 6 (1): 18- 10.1186/1471-2164-6-18PubMed CentralView ArticlePubMedGoogle Scholar
- Chorley BN, Wang X, Campbell MR, Pittman GS, Noureddine MA, Bell DA: Discovery and verification of functional single nucleotide polymorphisms in regulatory genomic regions: current and developing technologies. Mutat Res. 2008, 659 (1-2): 147-157. 10.1016/j.mrrev.2008.05.001PubMed CentralView ArticlePubMedGoogle Scholar
- Swanson HI: DNA binding and protein interactions of the AHR/ARNT heterodimer that facilitate gene activation. Chem Biol Interact. 2002, 141 (1-2): 63-76. 10.1016/S0009-2797(02)00066-2View ArticlePubMedGoogle Scholar
- Carter KW, McCaskie PA, Palmer LJ: JLIN: a java based linkage disequilibrium plotter. BMC Bioinformatics. 2006, 7: 60- 10.1186/1471-2105-7-60PubMed CentralView ArticlePubMedGoogle Scholar
- Kawashima K, Nomura S, Hirai H, Fukushi S, Karube T, Takeuchi K, Naruke T, Nishimura S: Correlation of L-myc RFLP with metastasis, prognosis and multiple cancer in lung-cancer patients. Int J Cancer. 1992, 50: 557-561. 10.1002/ijc.2910500411View ArticlePubMedGoogle Scholar
- Tamai S, Sugimura H, Caporaso NE, Resau JH, Trump BF, Weston A, Harris CC: Restriction fragment length polymorphism analysis of the L-myc gene locus in a case-control study of lung cancer. Int J Cancer. 1990, 46: 411-415. 10.1002/ijc.2910460314View ArticlePubMedGoogle Scholar
- Fox J: The R Commander: A basic-statistics graphical user interface to R. JStat Software. 2005, 19 (9): 1-42.Google Scholar
This article is published under license to 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.