MFSD2A is a novel lung tumor suppressor gene modulating cell cycle and matrix attachment
© Spinola et al. 2010
Received: 1 October 2009
Accepted: 17 March 2010
Published: 17 March 2010
MFSD2A (major facilitator superfamily domain containing 2) gene maps on chromosome 1p34 within a linkage disequilibrium block containing genetic elements associated with progression of lung cancer.
Here we show that MFSD2A expression is strongly downregulated in non-small cell lung cancer cell lines of different histotypes and in primary lung adenocarcinomas. Experimental modulation of MFSD2A in lung cancer cells is associated with alteration of mRNA levels of genes involved in cell cycle control and interaction with the extracellular matrix. Exogenous expression of MFSD2A in lung cancer cells induced a G1 block, impaired adhesion and migration in vitro, and significantly reduced tumor colony number in vitro (4- to 27-fold, P < 0.0001) and tumor volume in vivo (~3-fold, P < 0.0001). siRNA knockdown studies in normal human bronchial epithelial cells confirmed the role of MFSD2A in G1 regulation.
Together these data suggest that MFSD2A is a novel lung cancer tumor suppressor gene that regulates cell cycle progression and matrix attachment.
Cancer progression is defined as the stepwise process through which cells evolve towards a more malignant and aggressive phenotype . This process results from the accumulation of somatic genetic and epigenetic changes occurring within neoplastic cells . However, a growing body of evidence also points to the role of genetic background in cancer susceptibility, progression, and prognosis [3–5]. We previously identified a 106 kb linkage disequilibrium block containing genetic elements associated with survival in lung adenocarcinoma (ADCA) patients . The refined region maps to chromosome 1p34 and includes MYCL1, TRIT1 (tRNA isopentenyltransferase 1), and MFSD2A (major facilitator superfamily domain containing 2). While the role of MYCL1 and TRIT1 in lung tumor growth and development has been studied [6, 7], no information is available on MFSD2A. Thus, we addressed the functional role of MFSD2A in lung tumorigenesis.
Downregulation of MFSD2A in lung cancer
Measurement of MFSD2A mRNA levels in NSCLC cell lines and normal human bronchial epithelial cell (HBEC) lines (Additional file 1), normalizing the data to the average expression of HBECs, revealed downregulation of MFSD2A (2- to 44-fold) in 33/47 (70%) NSCLC cell lines but only in 4/20 (20%) HBEC lines (Fig. 1B).
A specific transcriptional profile is associated with MFSD2A over-expression
Changes in gene expression dependent on or correlated to MFSD2A were studied using two different approaches both based on microarray analysis.
Class comparison analysis of the expression profile of lung tumor cells transiently transfected with MFSD2A identified 460 genes whose expression differed significantly after MFSD2A exogenous expression (P < 0.001, FDR < 0.03) (Additional file 2). Microarray results were validated on a subset of 15 genes by real-time PCR (rho = 0.83, P < 0.0001). Analysis with the DAVID (Database for Annotation, Visualization and Integrated Discovery) Functional Annotation Tool  pointed to four main gene ontology categories: regulation of transcription, mitosis, apoptosis, and cell cycle.
Comparison of microarray data for NSCLC cell lines showing the highest and lowest MFSD2A mRNA levels identified 200 genes displaying at least a 4-fold difference (P < 0.005; Additional file 3). Most of the MFSD2A-correlated genes control developmental processes, neurodevelopment, cell motility, and adhesion.
MFSD2A controls tumor growth and G0/G1 phase of lung cells
Analysis of tumor growth in nude mice injected subcutaneously with control or MFSD2A-transfected A549 cells revealed xenograft tumors at 2 weeks after injection in all control cell-treated mice, whereas tumors were measurable in all mice treated with MFSD2A-transfected cells only at 5 weeks post-injection. At the end of the observation period, tumors grown from MFSD2A-expressing cells were, on average, ~3-fold smaller than those from control cells (P < 0.0001) (Fig. 3C).
To confirm this finding we studied changes in the cell cycle profile of normal human bronchial epithelial cells (HBECs) transfected with siRNAs targeting MFSD2A. Quantitative real-time PCR showed that 3 of 5 oligos (5'-catggagagtaacctcatcat-3', 5'-gagtgtcactgggcatttcta-3', 5'-ccactgtgaatatgccaagga-3') induced >3-fold reduction in MFSD2A mRNA levels. Low MFSD2A expression resulted in fewer cells accumulated in G1 as compared to control cells (from 59% to 41%) (Fig. 4, right panel), whereas distribution of cells in S and G2/M phases varied among the different samples (three independent replicates). Overall statistical analysis showed a significant effect of MFSD2A category of expression (over-expression or silencing) on modulation of cells accumulating in G1 fraction (P < 0.0001), with significant differences between the two cell lines (P = 0.003), and no significant interaction between MFSD2A expression and cell lines.
MFSD2A over-expression modulates cell adhesion and cell migration
The effects of MFSD2A expression on cancer cell migration were evaluated in a time-course experiment with the same ECM substrates used in the adhesion assay. No difference was seen in migration on collagen I-coated membranes, which also displayed the smallest effects in adhesion experiments (not shown). In contrast, the different migratory behavior of MFSD2A-transfected cells on collagen IV- and fibronectin-coated membranes was already apparent at early time points and was maintained throughout the whole experiment (Fig. 5B). At the final time point of the experiment, migration was reduced by ~60% on collagen IV and by ~80% on fibronectin. Statistical analysis indicated a significant difference between control and MFSD2A cells (P < 0.0001) as well as between substrates (collagen-IV and fibronectin) (P < 0.001).
List of extracellular matrix genes >3-fold downregulated by MFSD2A over-expression in stable-transfected A549 cells.
ADAM metallopeptidase with thrombospondin type 1 motif, 1
cadherin 1, type 1, E-cadherin (epithelial)
collagen, type XI, alpha 1
collagen, type XII, alpha 1
collagen, type IV, alpha 2
collagen, type V, alpha 1
collagen, type VI, alpha 1
collagen, type VII, alpha 1
connective tissue growth factor
catenin (cadherin-associated protein), delta 1
intercellular adhesion molecule 1
integrin, alpha 2 (CD49B, alpha 2 subunit of VLA-2 receptor)
integrin, alpha 7
laminin, beta 3
matrix metallopeptidase 1 (interstitial collagenase)
matrix metallopeptidase 10 (stromelysin 2)
matrix metallopeptidase 14 (membrane-inserted)
matrix metallopeptidase 15 (membrane-inserted)
secreted protein, acidic, cysteine-rich (osteonectin)
transforming growth factor, beta-induced, 68 kDa
TIMP metallopeptidase inhibitor 3
In this study, most of the primary lung tumors and NSCLC cell lines examined showed significantly lower MFSD2A expression as compared to their normal counterparts. The mechanisms underlying the inhibition of MFSD2A expression are unclear, although our preliminary experiments indicating restoration of MFSD2A expression in NSCLC cell lines upon treatment with 5'-azacytidine, a methyltransferase inhibitor and demethylating agent, suggest a role for methylation. Other possible mechanisms of gene inactivation, such as the presence of somatic mutations or the influence of microRNAs, remain to be investigated.
Over-expression of MFSD2A in transfected lung cancer cell lines was associated with reduced clonogenicity in vitro and diminished tumorigenicity in vivo, an effect due likely to the ability of MFSD2A to block the cell cycle in the G1 phase and to impair adhesive and migratory properties. Indeed, microarray data indicated that MFSD2A regulates expression of genes controlling cell cycle progression and correlates with expression of genes affecting adhesion and motility. Although microarray analysis also pointed to regulation of apoptotic genes, the FACS profile of MFSD2A-transfected cells revealed no increase in the sub-G1 phase fraction, no change was observed after annexin V staining, and neither PARP nor caspase cleavage was detected (data not shown).
MFSD2A has recently been described as the human receptor for syncytin-2, a retrovirus-derived protein mediating fusion of placental trophoblasts; however, in silico analysis predicts that the primary function of the gene is in transport of carbohydrates . The presence of this gene within the MYCL1 linkage disequilibrium block associated with differences in the survival of lung cancer patients  suggests a role for MFSD2A in controlling predisposition to lung cancer progression, although sequencing of coding regions identified no functional polymorphism in MFSD2A that could account for this effect . Studies addressing the possible existence of genetic variations in the promoter region and their role in populations of different ethnicities will shed light on the potential role of MFSD2A in modifying lung cancer progression.
Based on the present findings we can conclude that MFSD2A is a novel suppressor gene in lung cancer acting on tumor growth and development through control of cell cycle profile, matrix attachment, and cell motility.
Lung Primary Samples and Cell lines
Matched specimens of normal lung parenchyma and lung adenocarcinoma tissue were obtained from patients who underwent lobectomy at Istituto Nazionale Tumori (Milan, Italy). Normal human bronchial epithelial cells (HBECs) and lung tumor cell lines were available at the tissue culture repository of the Hamon Center for Therapeutic Oncology Research, UTSW Medical Center (Dallas, TX) . Cell lines that have been used are listed in Additional file 1.
Quantitative real-time PCR
Total RNA was extracted with the RNeasy Midi kit (Qiagen, Valencia, CA) and reverse-transcribed with either the SuperScript First-Strand Synthesis System (Invitrogen, Carlsbad, CA). MFSD2A mRNA expression was analyzed using TaqMan gene expression assays (Applied Biosystems, Foster City, CA): MFSD2A (Hs00293017_m1), HPRT1 (Hs99999909_m1), and GAPDH (4352934E).
Intron-spanning primers were designed to validate microarray results of NCI-H520 transfected cells for the following genes: DDIT3, DNAJB9, ELMO3, HRK, IRF1, LAMP3, MX1, NUPR1, PARL, PPP1R15A, RYR1, S100P, TCTA, UHMK1, WIPI1, and HPRT1 (housekeeping control). Amplification mixtures contained cDNA template, Power SYBR® Green PCR Master Mix (Applied Biosystems), and gene-specific PCR primers (sequences of oligonucleotide primers are available upon request).
For the expression levels analysis of genes involved in cell adhesion and cell migration pathways, we have used the TaqMan Array Human Extracellular Matrix & Adhesion Molecules 96-well plate (Applied Biosystems), according to the manufacturer's instructions.
Relative expression values were calculated using the comparative Ct method.
Gene expression profile of NCI-H520 cells transfected with recombinant MFSD2A or empty control vector (4 replicas/each) was analyzed using the Human-8 v3 Expression BeadChips (Illumina Inc., San Diego, CA, USA). Intensity values of each hybridization were quality-checked and the data set was normalized using a cubic spline algorithm, with BeadStudio Version 3 software. A P-value < 0.05 was set as a cutoff to filter reliably detected genes.
Gene expression profile of the 47 NSCLC cell lines was analyzed on the Affymetrix GeneChips HG-U133A and HG-U133B (together 44,928 elements; 23,583 unique genes) according to the manufacturer's protocol. Intensity values were quality-checked and the data set was quantile-normalized.
MFSD2A over-expression and silencing
MFSD2A full-length transcript (NM_032793) was amplified from a pool of human lung cDNAs (primers 5'-ggtcatggccaaaggagaa-3' and 5'-gaggatgctagccagctctgtg-3'), cloned in pEF6/V5-His TOPO vector (Invitrogen), and sequenced. Cell lines were transiently transfected using FuGENE HD Transfection Reagent (Roche, Basel, Switzerland). Stable transfectant clones were obtained after selection with 5 (A549) or 2 (NCI-H596, NCI-H520) μg/ml of blasticidin (Invitrogen) for 2 weeks.
siRNAs targeted against MFSD2A (5'-gcttcacaaagtgccaaccat-3', 5'-catggagagtaacctcatcat-3', 5'-gagtgtcactgggcatttcta-3', 5'-cacggcccatacatcaaactt-3', 5'-ccactgtgaatatgccaagga-3') were purchased from Qiagen. HBECs cells, which express endogenous MFSD2A, were transfected with 25 nM siRNA oligonucleotides using Oligofectamine (Invitrogen) and harvested after 96 hours for quantitative real-time PCR or cell cycle analysis.
Western blots and immunohistochemistry
MFSD2A expression in stably transfected cells was detected by Western blotting with anti-V5-HRP antibody (Invitrogen).
Samples of paraffin-embedded sections of MFSD2A- or empty vector-transfected HEK-293T cells, lung ADCA tissue, and surrounding normal lung tissue were used to prepare histological sections that were immunostained using standard methods after antigen retrieval performed in 0.07 M citrate buffer (pH = 6) at 95°C for 10 min. Mouse polyclonal anti-MFSD2A antibody (H00084879-B01P; Abnova, Taipei City, Taiwan) was used at a 1:170 dilution. Immunoreactive signals were detected with ChemMate DAB (Dako, Glostrup, Denmark).
In vitro assays
Colony formation assay was carried out in stably transfected A549, NCI-H520, and NCI-H596 cells (6 replicas). After 2 weeks of selection, colonies were methanol-fixed, stained with 10% Giemsa, and manually counted.
Cell adhesion was measured using CAFCA (Centrifugal Assay for Fluorescence based Cell Adhesion) assay . Briefly, six well strips were coated with different substrates from Becton-Dickinson (Falcon, Milan, Italy). Cells were labeled with the vital fluorochrome calcein AM (Invitrogen) for 15 minutes at 37°C and aliquoted into the bottom CAFCA miniplates, which were centrifuged to synchronize the contact of the cells with the substrate. The miniplates were incubated for 20 minutes at 37°C and mounted together with a similar CAFCA miniplate to create communicating chambers for reverse centrifugation. The relative number of cells bound to the substrate (i.e. remaining in the wells of the bottom miniplates) and cells that failed to bind to the substrate (i.e. remaining in the wells of the top miniplates) was estimated by top/bottom fluorescence detection with GENios Plus microplate fluorometer (TECAN, Italy). Percentage of adherent cells was determined 8 hours after plating.
Cell migration in response to extracellular matrix substrates was assessed by FATIMA (Fluorescence-Assisted Transmigration Invasion and Motility) assay as described . Briefly, membranes of HTS FluoroBlokTM transwell inserts with 8 μm pores (Becton-Dickinson, Falcon) were coated on the underside with various ECM molecules at 4°C and blocked with 1% BSA. Cells were fluorescently tagged with 5 μg/ml DiI lipophilic dye (Invitrogen) and plated in the upper chamber. Migration was monitored at every hour by independent fluorescence detection from the top (corresponding to non-transmigrated cells) and bottom (corresponding to transmigrated cells) side of the membrane with GENios Plus microplate fluorometer (TECAN).
Cell cycle analysis of stably transfected cells was performed on a FACScan (BD Biosciences, San Jose, CA) or FACSCalibur flow cytometer (BD Biosciences). MFSD2A- or empty vector-transfected cells were harvested and fixed in 70% EtOH over night at 4°C. Cells were then incubated in 500 μl of buffered propidium iodide (PI) staining solution containing 0.05% Triton X-100 (Sigma-Aldrich), 0.1 mg/ml RNase A (Millipore), and 50 g/ml PI (Sigma-Aldrich) in PBS for 30 min at 37°C. Cells were briefly spun down and resuspended in PBS. Flow cytometry data were analyzed using FlowJo software (TreeStar, Ashland, OR) with the aid of the Watson modeling algorithm.
In vivo tumor growth assay
Adult female CD-1 nude mice (purchased from Charles River, Calco, Italy) were injected subcutaneously with 4 × 106 A549 control or MFSD2A-transfected cells (2 independent clones). Tumor size was measured weekly by calipers and animals were sacrificed 8 weeks after injection.
Genes differentially expressed in the two classes of either vector- or MFSD2A-transfected cells were identified using random variance t-statistics . Microarray data were analyzed using BRB ArrayTools developed by Dr. Richard Simon and Amy Peng Lam http://linus.nci.nih.gov/BRB-ArrayTools.html (Illumina experiment) or using in-house Visual Basic software MATRIX 1.4 and the Bioconductor R package affy (Affymetrix experiment). Functional annotation in gene ontology (GO) categories was carried out with the DAVID Functional Annotation Tools . Differences in quantitative measures were assessed by analysis of variance. Correlations between expression levels were expressed by the Spearman's correlation coefficient. Cell cycle parameters by MFSD2A categories of expression (high, low) were analyzed using the generalized linear model, family = quasibinomial, procedure.
We thank Dr. A. Gazdar for discussion and cell lines and we thank Ms. Lucia Gioiosa for technical assistance in immunohistochemistry. MS was supported by a fellowship from Associazione Italiana Ricerca Cancro (AIRC). This work was funded in part by grants from Associazione and Fondazione Italiana Ricerca Cancro (AIRC and FIRC), Fondo Investimenti Ricerca di Base (FIRB), Italy, NCI Lung Cancer SPORE (P50CA70907), and Gillson Longenbaugh Foundation.
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.