- Primary research
- Open Access
PITX2 enhances progression of lung adenocarcinoma by transcriptionally regulating WNT3A and activating Wnt/β-catenin signaling pathway
Cancer Cell International volume 19, Article number: 96 (2019)
The homeodomain transcription factor, PITX2 is associated with tumorigenesis of multiple cancers. In this research, we aimed to study the expression, function and mechanism of PITX2 in lung adenocarcinoma (LUAD).
The TCGA dataset was used to analyze the expression and clinical significance of PITX2 in LUAD. The expression of PITX2 in tumor samples and LUAD cell lines was examined by quantitative real-time PCR (qRT-PCR) and western blotting. Small interfering RNAs (siRNAs) were constructed to knockdown PITX2 and to determine the physiological function of PITX2 in vitro. Xenograft model was used to confirm the role of PITX2 in vivo.
PITX2 was overexpressed in LUAD and patients with high level of PITX2 had a worse overall survival and an advanced clinical stage. Knockdown of PITX2 inhibited cell proliferation, migration and invasion of LUAD cells. Further study revealed that the oncogenic role of PITX2 was dependent on activating Wnt/β-catenin signaling pathway, especially by transcriptionally regulating the Wnt gene family member, WNT3A. Lastly, we identified miR-140-5p as a negative mediator of PITX2 by binding its 3′UTR and ectopic expression of miR-140-5p inhibited progression of LUAD cells via suppressing the expression of PITX2.
Up-regulation of PITX2 acts as an oncogene in LUAD by activating Wnt/β-catenin signaling pathway, suggesting that PITX2 may serve as a novel diagnostic and prognostic biomarker in LUAD.
Lung adenocarcinoma (LUAD) is the most common histologic subtype of lung cancer, with an increasing morbidity and mortality rate over the past decades [1, 2]. Multiple genomic alterations are discovered to be associated with the progression of LUAD and significant advances have been achieved in molecularly targeted therapy . However, drug resistance and side effects weaken the therapeutic effects and the overall 5-year survival rate of LUAD is still less than 20% . Thus, a further understanding of the cellular mechanisms that regulate LUAD progression may provide alternative diagnostic markers and therapeutic targets for LUAD patients.
The homeobox genes are a group of genes containing homeobox (a 180 base pairs long DNA sequence) and encode homeodomain proteins that act as transcription factors . Many homeobox genes have been proved to be implicated in various human cancers, acting as oncogenes or tumor suppressors [6,7,8]. Paired-like homeodomain transcription factor 2 (PITX2), a member of the bicoid/paired-like homeobox gene family, is located on chromosome 4q25 and encodes a transcription factor involved in the development of anterior structures [9, 10]. Recent studies have highlighted the association of PITX2 with progression of many human cancers, including gonadotroph tumor, colorectal cancer, thyroid cancer and esophageal squamous cell carcinoma [11,12,13,14]. However, the expression of PITX2 in LUAD and its effect on biological behavior of LUAD have not been elucidated.
Earlier researches have discovered that homeobox genes are actively involved in the regulation of the Wnt pathway [15,16,17]. Particularly, PITX2 is reported to regulate the Wnt/β-catenin pathway by activating Wnt ligand genes and thus contributing to cancer progression . The Wnt/β-catenin signaling pathway was first identified for its role in carcinogenesis and meanwhile it was confirmed to be essential for embryonic development [19, 20]. Herein we investigated the influence of PITX2 on activity of Wnt/β-catenin pathway in LUAD.
MiRNAs, a class of 18–25 nucleotide non-coding RNAs, have emerged as key regulators of gene expression at the post-translational level by combining with the 3′UTRs (untranslated regions) of their target mRNAs . The miRNA-mediated gene post-transcriptional regulation is involved in various aspects of physiological processes . And studies are increasingly implicating miRNAs in the progression of many types of cancer, including LUAD [23,24,25]. Since the regulatory miRNAs could suppress target mRNAs and play a critical role in pathogenesis, we performed experiments to determine whether PITX2 is targeted by specific miRNA.
Materials and methods
To determine the expression of PITX2 in LUAD and normal lung tissues, the TCGA dataset named TCGA_LUAD_exp_Hi-SeqV2-2015-02-24, were downloaded from the UCSC Cancer Browser (https://genomecancer.ucsc.edu/) . The “genomicMatrix” file was browsed to get the normalized expression of PITX2 in LUAD and the clinicopathological information of patients was acquired from “clinical_data” file. The cBioPortal website (http://www.cbioportal.org/)  was used to obtain highly co-expressed genes with PITX2. Totally 317 genes (Additional file 1: Table S1) were submitted to DAVID Bioinformatics Resources 6.8 (http://david.abcc.ncifcrf.gov/)  and Reactome Pathway Database (https://reactome.org/)  for Gene_Ontology pathway and Reactome pathway analysis, respectively. The protein association networks of PITX2 were derived from The String Database (https://string-db.org/) . The transcription factor binding site (TFBS) of PITX2 was exhibited in Human Transcription Factor Database Version3.0 (http://bioinfo.life.hust.edu.cn/AnimalTFDB/#!/)  and hTFtarget database (http://bioinfo.life.hust.edu.cn/hTFtarget#!/) was utilized to analyze the potential binding site of PITX2 with target genes. To search for miRNAs that directly target PITX2, CoMeTa dataset (http://cometa.tigem.it/index.php) , DIANA lab tools webpage (http://diana.imis.athena-innovation.gr/DianaTools/index.php?r=site/index)  and the microRNA database (http://www.microrna.org/microrna/home.do)  were used.
Patient tissue samples
60 pairs of LUAD tissues were obtained from patients who had undergone radical surgical resection in the department of thoracic surgery, Nanjing Jinling Hospital. All tissues were examined by two experienced pathologists who confirmed the TNM stage of each patient. The clinical information of patients was consulted from their medical records. And the study was approved by Ethics Committee of Nanjing Jinling Hospital and informed consent was obtained from all patients.
RNA extraction, reverse transcription, and qRT-PCR
Total RNA of tissue or cell samples were extracted with TRIzol Reagent (Invitrogen). The extracted RNA was dissolved with RNAase-free water, and 1.5 μg RNA from each sample was used for cDNA synthesis using a Reverse Transcription Kit (Takara, cat: RR036A). SYBR Select Master Mix (Applied Biosystems, Cat: 4472908) was used for qRT-PCR and primers are shown in Additional file 2: Table S2. The comparative CT method (ΔΔCT) was used to measure the relative gene expression.
Paraffin-embedded sections were deparaffinized and incubated in retrieval buffer solution for antigen retrieval. The sections were incubated with corresponding antibodies: anti-PITX2 (Sigma-Aldrich, AV35634, 1:100), anti-WNT3A (Abcam, ab28472, 1:100), anti-CyclinD1 (Cell Signaling Technology, 2978, 1:100), anti-C-Myc (R&D Systems, AF3696, 1:100).
Cell lines and cell culture
LUAD cell lines (NCIH1975, NCIH1650, PC9, A549, NCIH1299) and human bronchial epithelial cell (HBE) were purchased from Shanghai Institutes for Biological Science, China. NCIH1975, NCIH1650 and HBE cells were cultured in DMEM medium (KeyGene), and PC9, A549 and NCI-H1299 cells were cultured in RPMI 1640 medium (KeyGene), supplemented with 10% fetal bovine serum. All cells were cultured at 37 °C in a humidified incubator containing 5% CO2.
Western blotting was performed according to standard protocols as previously described . Briefly, cultured cells were lysed with lysis buffer (RIPA, KeyGEN) containing protease inhibitors (PMSF, KeyGEN) on ice, and protein concentration was determined using a BCA Kit (KeyGEN). Western blotting were obtained utilizing 30 μg of lysate protein. The antibodies used were as follows: anti-β-actin (Cell Signaling Technology, 3700, 1:1000), anti-PITX2 (Abcam, ab55599, 1:1000), anti-CyclinD1 (Cell Signaling Technology, 2978, 1:1000), anti-CD44 (Cell Signaling Technology, 3570, 1:1000), anti-C-Myc (R&D Systems, AF3696, 1:1000), anti-C-Jun (Cell Signaling Technology, 9165, 1:1000), anti-WNT3A (Abcam, ab28472, 1:100).
siRNAs against PITX2 (Santa Cruz Biotechnology, sc-44016) were used at 20 nM/well using 2 μL Lipofectamine RNAiMAX (Invitrogen) in the cells seeded in 6-well plates. The full-length cDNA of human PITX2 were PCR-amplified and cloned into the expression vector pENTER (Vigene Biosciences). shRNA and microRNA mimics were conducted and purchased from Vigene Biosciences. Transfection of plasmids and miRNA mimics was performed according to the Lipofectamine 3000 Reagent (Invitrogen) protocol. Totally 10 μg plasmid vectors were transfected into cells seeded in 6-well plates for each group. Nonsense shRNA (sh-nc) and negative control mimics were used as the respective controls. The sequences of shRNAs and microRNA mimics are shown in Additional file 3: Table S3.
Cell proliferation assay
For the Cell Counting Kit-8 (CCK-8) assay, cells were plated in 96-well plates at a density of 4000 cells/100 μL, and the absorbance was measured at 450 nm in an ELx-800 Universal Microplate Reader. Real time xCELLigence system (RTCA) assay was performed according to the instructions provided by the manufacturer . And cells were placed in 6-well plates at 400 cells per well for colony formation assays. After 10–14 days, the cells were fixed using 4% paraformaldehyde and stained with 0.1% crystal violet solution.
Migration and invasion assay
Transwell assay inserts (8 μM PET, 24-well Millicell) or Matrigel-coated membranes (BD Biosciences) were add with 200 μL serum-free medium in the upper chamber and 800 μL serum-containing medium in the lower chamber. Totally 50,000 cells were added to the upper chamber and cells in lower chamber were fixed and stained after 24 h (migration) or 48 h (invasion).
Flow cytometry assay
Flow-cytometry analysis was performed to detect cell cycle distribution. Cells were harvested after transfection with siRNA for 48 h and fixed in 1 mL 70% ethanol at − 20 °C. The ethanol-suspended cells were centrifuged and stained with PI staining solution for 10 min in the dark at 37 °C. A FACSCalibur flow cytometer was used to detect cell cycle distribution. The percentage of cells in G1, S, and G2-M phases were counted and compared.
TUNEL (TdT-mediated dUTP nick-end labeling) assay
TUNEL Kits (KeyGEN BioTECH, KGA7062) were used to detect DNA fragmentation of apoptotic cells. Cells were seeded in 24-well plates and fixed with 4% paraformaldehyde. After washing with PBS and 0.5% Triton X-100, cells were incubated at 37 °C for 1 h from light with terminal deoxynucleotidyl transferase enzyme and cell nuclei were counterstained with DAPI. Images were obtained from a fluorescence microscope and apoptosis ratio was calculated.
Luciferase reporter assay
The WNT3A promoter region (− 2000 bp) or PITX2 3′UTR was amplified and cloned into luciferase reporter plasmid (pGL3-basic). The luciferase reporter plasmids were co-transfected with CMV-PITX2 expression plasmids or miRNA mimics. And CMV-empty vector or miRNA control mimics were used as negative controls. CMV-Renilla was used as an internal control and Firefly/Renilla value was used to measure luciferase activity. After transfection for 48 h, cells were harvested and assessed for luciferase activity using the Dual Luciferase Reporter Assay System (Promega).
ChIP (chromatin immunoprecipitation) assay
Cells were cross-linked in 4% paraformaldehyde and the reaction was quenched with 10× glycine. Then cells were added with pre-cooling PBS (containing 5 μL protease inhibitor, cocktail) and scraped into a centrifuge tube, followed by centrifuging for 10 min at 800×g at 4 °C. Cell precipitates were resuspended in 500 μL cell lysis buffer (containing 2.5 μL cocktail) and 500 μL nucleus lysis buffer (containing 2.5 μL cocktail) successively. After that, the cells were sonicated (amplitude 30%) on ice to break DNA into fragments and were centrifuged for 10 min at 12,000g at 4 °C. For immunoprecipitation, the above products were added with 5 μg anti-PITX2 antibody (Abcam, ab192495) or 5 μg mouse IgG (Santa Cruz Biotechnology, sc-2025), all of which were mixed with 100 μL Dilution buffer and 20 μL magnetic beads (Resin M2, Sigma), and then incubated at 4 °C overnight. For DNA purification, immunoprecipitation products were added with ChIP elution buffer (containing proteinase K) and then was heated at 62 °C for 2 h and 95 °C for 10 min. Finally, the DNA was eluted with elution buffer in adsorption column and used for PCR reaction. The sequence of primers used are shown in Additional file 2: Table S2 and GAPDH was used as a control.
LUAD tumor cells were transfected with shRNAs or control plasmids before transplantation. Six female nude mice (4–6 weeks old) per group were purchased from Nanjing Medical University School of Medicine’s accredited animal facility. 2 × 106 LUAD tumor cells in 100 μL sterile PBS were injected subcutaneously to each mouse and tumors were harvested at 6 weeks after injection. The weight of tumor was measured on the scale and tumor volume was estimated using calipers ([length * width2]/2). And the tumor tissues were immunohistochemically stained with corresponding antibodies. All animal studies were conducted in accordance with NIH animal use guidelines and protocols were approved by Nanjing Medical University Animal Care Committee.
Generally experiments were repeated for three times and data were presented as the mean ± SD. Student’s t test was used to assess the differences between two groups. Chi square test was performed to analyze the differences between variables. Bivariate correlations between study variables were calculated by Pearson’s correlation coefficients. Kaplan–Meier method was used to plot survival curves. Statistical analysis was performed using SPSS Statistics software (version 20.0). p < 0.05 was considered statistically significant. Graphs were made using the GraphPad Prism 6.0 software package.
PITX2 was overexpressed in LUAD and correlated with worse prognosis
By analyzing the TCGA dataset (TCGA_LUAD_exp_Hi-SeqV2-2015-02-24), the expression of PITX2 was upregulated in LUAD tissues (n = 511) compared with normal tissues (n = 58) (Fig. 1a). And for the 511 LUAD tissues included in TCGA dataset, patients with higher PITX2 expression (the upper quartile, n = 128) had a more dismal overall survival than those with lower PITX2 expression (the lower quartile, n = 128) (Fig. 1b). The Chi squared test suggested that the PITX2 level is also related to T stage (p = 0.00311), M stage (p = 0.0367) and TNM stage (p = 0.010652) (Table 1). For further verification, we examined the levels of PITX2 in 60 pairs of LUAD and non-tumor tissues collected from Department of Cardiothoracic Surgery, Jinling Hospital. The results showed that PITX2 was upregulated in 76.67% of LUAD tissues (Fig. 1c). The clinical data of these patients was shown in Additional file 4: Table S4. Unpaired t test was used to validate the correlation between PITX2 and clinical pathological characteristics. The expression of PITX2 in patients of stage T2 (7.116 ± 1.130, N = 27), T3 (15.51 ± 3.166, N = 6) and T4 (20.54 ± 6.943, N = 4) was higher than patients of stage T1 (3.430 ± 0.8230, N = 23) (Fig. 1d). And patients with lymph node metastasis (14.05 ± 3.054, N = 12) had a higher PITX2 expression than patients without lymph node metastasis (5.783 ± 0.8629, N = 48) (Fig. 1e). Moreover, patients of TNM stage II (11.63 ± 1.179, N = 25) and III (16.55 ± 6.693, N = 5) had a higher PITX2 expression than patients of stage I (2.427 ± 0.5142, N = 30) (Fig. 1f). In addition, the mRNA and protein level of PITX2 were elevated in LUAD cell lines (NCIH1975, NCIH1650, PC9, A549, NCIH1299) compared with Human Bronchial Epithelial (HBE) Cells (Fig. 1g, h). The data suggested that PITX2 was a potential upregulated oncogene in LUAD.
Knockdown of PITX2 inhibited progression of LUAD in vitro
To examine the biological function of PITX2 in LUAD, two siRNAs against PITX2 were transfected into NCIH1299 cell lines. Confirmed by qRT-PCR and western blotting (Fig. 2a), both siRNAs effectively silenced the expression of PITX2 (interfering efficiency, si1 75.3%; si2 75.29%). Furtherly we found that knockdown of PITX2 inhibited proliferation of NCIH1299 cells, validated by CCK-8 assay (Fig. 2b), RTCA assay (Fig. 2c) and colony formation assay (Fig. 2d). Transwell (Fig. 2e) and matrigel (Fig. 2f) assay revealed that knockdown of PITX2 suppressed the migration and invasion ability of LUAD cells. In addition, flow cytometry assay suggested that knockdown of PITX2 induced G1/S cell cycle arrest (Fig. 2g) and TUNEL staining proved that silencing PITX2 led to higher apoptosis ratio of NCIH1299 cells (Fig. 2h).
PITX2 displayed its biological function via activating Wnt/β-catenin signaling pathway
To analyze the signaling pathways associated with PITX2, 317 genes highly co-expressed with PITX2 were used for Gene_Ontology Pathway (Fig. 3a) and Reactome Pathway analysis (Fig. 3b). Results showed that PITX2 might participate in “signal transduction”, “apoptotic process”, “regulation of proliferation” and “gene transcription”. In String database, PITX2 was predicted to interact with LEF1, CTNNB1 and AXIN, all of which are the key factors of Wnt/β-catenin signaling pathway (Fig. 3c). As previous studies have reported that Wnt/β-catenin pathway is regulated by PITX2 in many diseases [18, 37, 38], we performed experiments to test the influence of PITX2 on Wnt/β-catenin downstream genes. It turns out that knockdown of PITX2 significantly downregulated the expression of CCND1, CD44, MYC and JUN on both mRNA and protein levels (Fig. 3d, e). The TOP/FOP flash luciferase assay, a method to examine the activity of Wnt/β-catenin pathway, revealed that reporter expression was suppressed by knockdown of PITX2 (Fig. 3f). And Wnt agonist 1, an activator of Wnt/β-catenin pathway, rescued the effect of PITX2 on proliferation (Fig. 3g) and apoptosis (Fig. 3h) in LUAD cells. These results suggested that the oncogenic role of PITX2 was dependent on Wnt/β-catenin pathway.
PITX2 regulated the transcription of WNT3A in LUAD cells
It has been confirmed that the encoded protein of PITX2 acts as a transcription factor to regulate gene expression. As shown by Animal TF DataBase, the TF binding site of PITX2 was presented (Fig. 4a). WNT3A, one of the most important ligands in canonical Wnt/β-catenin signaling pathway, has been well established in the progression of LUAD [39,40,41]. Predicted by hTFtarget website, PITX2 potentially targets the promoter region (1092–1101 bp) of WNT3A (Fig. 4b). Then we constructed luciferase reporter plasmids carrying WNT3A promoter (wild type and mutant type), and luciferase reporter assay revealed that PITX2 increased wild type but not mutant WNT3A promoter luciferase activity (Fig. 4c). Moreover, the results of ChIP assay suggested that PITX2 directly targeted the promoter region of WNT3A (Fig. 4d, e). In NCIH1299 cell lines, knockdown of PITX2 suppressed the expression of WNT3A (Fig. 4f). And in TCGA dataset, PITX2 had a positive correlation with WNT3A (R = 0.22, p < 0.001) (Fig. 4g). These findings indicated that PITX2 controlled the activity of Wnt/β-catenin signaling pathway via transcriptional regulating WNT3A.
Knockdown of PITX2 inhibited tumor growth of LUAD cells in vivo
To investigate the oncogenic role of PITX2 in vivo, xenograft tumor models (six nude mice per group) were used and tumors were harvested at 6 weeks after transplantation (Fig. 5a). The volume and weight of tumors were measured and results showed that tumor growth was significantly suppressed by knockdown of PITX2 (Fig. 5b, c). And immunohistochemistry (IHC) staining with corresponding antibodies revealed that knockdown of PITX2 downregulated the expression of WNT3A, CyclinD1 and C-Myc (Fig. 5d).
MiR-140-5p was a direct target of PITX2 and suppressed progression of LUAD
To determine the regulatory miRNAs of PITX2, three miRNA prediction datasets were analyzed and the results showed that miR-140-5p potentially bound to PITX2 (Fig. 6a). According to the binding site, luciferase reporter plasmids carrying PITX2 3′UTR (wild type and mutant type) were constructed (Fig. 6b) and luciferase reporter assay showed that ectopic expression of miR-140-5p decreased the luciferase activity of wild type PITX2 3′UTR but not mutant type (Fig. 6c). Furtherly confirmed by qRT-PCR and western blotting, the expression of both PITX2 and WNT3A were down-regulated by ectopic expression of miR-140-5p (Fig. 6d, e). And by transfecting miR-140-5p mimics into NCIH1299 cells, we found that miR-140-5p inhibited tumor proliferation (Fig. 6f), reduced the transcriptional activity of Wnt/β-catenin (Fig. 6g) and induced apoptosis (Fig. 6h) in LUAD cells. Meanwhile, these biological functions of miR-140-5p could be restored by overexpression of PITX2 (Fig. 6f–h). As miR-140-5p played vital roles in LUAD progression, we examined the expression of miR-140-5p in 60 pairs of LUAD tissues. Results revealed that miR-140-5p was significantly down-regulated in LUAD tissues (Fig. 6i) and had a negative correlation with PITX2 (p < 0.0001, R = − 0.5866) (Fig. 6j).
In this study, we found that PITX2 was significantly upregulated in LUAD tissues and correlated with dismal overall survival and advanced clinical stages of patients. Knockdown of PITX2 inhibited tumor growth of LUAD both in vitro and in vivo. And PITX2 exerted its oncogenic role via activating Wnt/β-catenin pathway, especially by regulating the transcription of WNT3A. Furthermore, we found that miR-140-5p was a target of PITX2 and suppressed the progression of LUAD (Fig. 7).
The progression of cancer is closely linked to activation of oncogenes and inactivation of tumor suppressor genes [42,43,44]. LUAD is the most predominant histological subtype of lung cancer and multiple genomic alterations function in cancer growth [45, 46]. The best known is that approximately 15–20% of LUAD harbor epidermal growth factor receptor (EGFR) mutations . However, the functional impact of most genomic alterations of LUAD remains largely unknown. To better understand novel functional genomic alterations in LUAD is hopefully to provide new treatment targets for LUAD patients. As a member of homeobox gene, PITX2 plays a role in the terminal differentiation of somatotroph and lactotroph cell phenotypes , is involved in the development of the eye , tooth and abdominal organs , and acts as a transcriptional regulator involved in basal and hormone-regulated activity of prolactin . Recently, more and more studies show the critical role of PITX2 in many kinds of cancers. In breast cancer cells, upregulation of PITX2 promotes letrozole resistance via transcriptional activation of IFITM1 signaling . The PITX2 DNA methylation was confirmed as a biomarker for the risk and prognosis of colorectal carcinoma , head and neck squamous cell carcinoma , urothelial bladder cancer  and prostate cancer . Increasing researches are showing the significance of development-related genes in cancer progression. KDM5B, a critical regulator during early embryonic development, plays a critical role in controlling stem cell function and tumorigenesis . GREM1, a member of antagonist family that has been shown to relay the sonic hedgehog (SHH) signal from the polarizing region to the apical ectodermal ridge during limb bud outgrowth in mouse, is involved progression of various cancers [58,59,60]. Similarly, we identified the oncogenic role of PITX2 in LUAD.
In aspect of mechanism, we found that PITX2 promoted LUAD progression as least partially by activating Wnt/β-catenin pathway. Wnt/β-catenin pathway causes an accumulation of β-catenin in the cytoplasm and its eventual translocation into the nucleus to act as a transcriptional coactivator of transcription factors that belong to the TCF/LEF family . Wnt/β-catenin pathway is involved in virtually every aspect of embryonic development and also controls homeostatic self-renewal in a number of adult tissues . Epigenetic dysregulation of Wnt/β-catenin pathway plays a significant role in the pathogenesis of several human cancers [63,64,65]. The targeting of Wnt protein to the plasma membrane and the binding to its receptor are the initiation process of activating Wnt/β-catenin pathway . As an important member of WNT gene family, WNT3A enhances β-catenin-dependent transcription through GSK3β inhibition or direct β-catenin phosphorylation at Ser552 . WNT3A has been reported to promote the migration and invasion process of LUAD and impact prognosis of LUAD patients . And our results suggested that PITX2 directly bound to the promoter region of WNT3A and enhanced its transcription. The oncogenic role of PITX2 is relied on its transcription of WNT3A and thus activating Wnt/β-catenin pathway.
MicroRNAs regulate post-transcriptional gene expression via binding to the 3′UTR of mRNA, and either inhibit protein translation or destabilize target mRNA . Technologies aimed at replacing tumor suppressor miRNA that are lost during cancer progression have emerged as promising cancer therapies . As we have clarified the role of PITX2 in LUAD, to determine the regulatory miRNA of PITX2 is hopefully to provide alternative approach for managing LUAD. Our results revealed that miR-140-5p directly targeted the 3′UTR of PITX2 and suppressed the progression of LUAD. MiR-140-5p has been reported to inhibit progression of many kinds of cancers and regulate the Wnt/β-catenin pathway [69,70,71]. The miR-140-5p-PITX2-WNT3A axis participates in the regulation of LUAD progression and Wnt/β-catenin pathway, which provides a new sight into the pathogenesis of LUAD.
To sum up, our study identified PITX2 as a novel oncogene in LUAD. PITX2 promoted progression of LUAD partially by activating Wnt/β-catenin pathway, especially by transcriptional regulating the crucial member of Wnt ligands family, WNT3A. Further research revealed that miR-140-5p directly targeted the 3′UTR of PITX2 and suppressed the progression of LUAD. The miR-140-5p-PITX2-WNT3A axis provides a new sight into the tumorigenesis of LUAD and may represent a potential biomarker and a promising therapeutic target for LUAD patients.
The Cancer Genome Atlas
Siegel RL, Miller KD, Jemal A. Cancer statistics, 2017. CA Cancer J Clin. 2017;67:7–30.
Jemal A, Bray F, Center MM, Ferlay J, Ward E, Forman D. Global cancer statistics. CA Cancer J Clin. 2011;61:69–90.
Hirsch FR, Scagliotti GV, Mulshine JL, Kwon R, Curran WJ Jr, Wu YL, Paz-Ares L. Lung cancer: current therapies and new targeted treatments. Lancet. 2016;389:299–311.
Chen Z, Fillmore CM, Hammerman PS, Kim CF, Wong KK. Non-small-cell lung cancers: a heterogeneous set of diseases. Nat Rev Cancer. 2014;14:535–46.
Gehring WJ. The homeobox in perspective. Trends Biochem Sci. 1992;17:277–80.
Dalerba P, Sahoo D, Clarke MF. CDX2 as a prognostic biomarker in colon cancer. N Engl J Med. 2016;374:2184.
Liu Z, Yanagisawa K, Griesing S, Iwai M, Kano K, Hotta N, Kajino T, Suzuki M, Takahashi T. TTF-1/NKX2-1 binds to DDB1 and confers replication stress resistance to lung adenocarcinomas. Oncogene. 2017;36:3740–8.
Terrinoni A, Pagani IS, Zucchi I, Chiaravalli AM, Serra V, Rovera F, Sirchia S, Dionigi G, Miozzo M, Frattini A, Ferrari A, Capella C, Pasquali F, Curto FL, Albertini A, Melino G, Porta G. OTX1 expression in breast cancer is regulated by p53. Oncogene. 2011;30:3096–103.
Logan M, Pagan-Westphal SM, Smith DM, Paganessi L, Tabin CJ. The transcription factor Pitx2 mediates situs-specific morphogenesis in response to left–right asymmetric signals. Cell. 1998;94:307–17.
Campione M, Steinbeisser H, Schweickert A, Deissler K, van Bebber F, Lowe LA, Nowotschin S, Viebahn C, Haffter P, Kuehn MR, Blum M. The homeobox gene Pitx2: mediator of asymmetric left-right signaling in vertebrate heart and gut looping. Development. 1999;126:1225–34.
Acunzo J, Roche C, Defilles C, Thirion S, Quentien MH, Figarella-Branger D, Graillon T, Dufour H, Brue T, Pellegrini I, Enjalbert A, Barlier A. Inactivation of PITX2 transcription factor induced apoptosis of gonadotroph tumoral cells. Endocrinology. 2011;152:3884–92.
Hirose H, Ishii H, Mimori K, Tanaka F, Takemasa I, Mizushima T, Ikeda M, Yamamoto H, Sekimoto M, Doki Y, Mori M. The significance of PITX2 overexpression in human colorectal cancer. Ann Surg Oncol. 2011;18:3005–12.
Huang Y, Guigon CJ, Fan J, Cheng SY, Zhu GZ. Pituitary homeobox 2 (PITX2) promotes thyroid carcinogenesis by activation of cyclin D2. Cell Cycle. 2010;9:1333–41.
Zhang JX, Chen ZH, Xu Y, Chen JW, Weng HW, Yun M, Zheng ZS, Chen C, Wu BL, Li EM, Fu JH, Ye S, Xie D. Downregulation of MicroRNA-644a promotes esophageal squamous cell carcinoma aggressiveness and stem cell-like phenotype via dysregulation of PITX2. Clin Cancer Res. 2017;23:298–310.
McWhirter JR, Neuteboom ST, Wancewicz EV, Monia BP, Downing JR, Murre C. Oncogenic homeodomain transcription factor E2A-Pbx1 activates a novel WNT gene in pre-B acute lymphoblastoid leukemia. Proc Natl Acad Sci USA. 1999;96:11464–9.
Riddle RD, Ensini M, Nelson C, Tsuchida T, Jessell TM, Tabin C. Induction of the LIM homeobox gene Lmx1 by WNT7a establishes dorsoventral pattern in the vertebrate limb. Cell. 1995;83:631–40.
Lagutin OV, Zhu CC, Kobayashi D, Topczewski J, Shimamura K, Puelles L, Russell HR, McKinnon PJ, Solnica-Krezel L, Oliver G. Six3 repression of Wnt signaling in the anterior neuroectoderm is essential for vertebrate forebrain development. Genes Dev. 2003;17:368–79.
Basu M, Roy SS. Wnt/beta-catenin pathway is regulated by PITX2 homeodomain protein and thus contributes to the proliferation of human ovarian adenocarcinoma cell, SKOV-3. J Biol Chem. 2013;288:4355–67.
Nusse R, Varmus HE. Wnt genes. Cell. 1992;69:1073–87.
Nusse R. Wnt signaling in disease and in development. Cell Res. 2005;15:28–32.
Bartel DP. MicroRNAs: target recognition and regulatory functions. Cell. 2009;136:215–33.
Bartel DP. MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004;116:281–97.
Zhang C, Wang H, Liu X, Hu Y, Ding L, Zhang X, Sun Q, Li Y. Oncogenic microRNA-411 promotes lung carcinogenesis by directly targeting suppressor genes SPRY4 and TXNIP. Oncogene. 2018;38:1892–904.
Chen S, Wu J, Jiao K, Wu Q, Ma J, Chen D, Kang J, Zhao G, Shi Y, Fan D, Zhao G. MicroRNA-495-3p inhibits multidrug resistance by modulating autophagy through GRP78/mTOR axis in gastric cancer. Cell Death Dis. 2018;9:1070.
Gao W, Zhang C, Li W, Li H, Sang J, Zhao Q, Bo Y, Luo H, Zheng X, Lu Y, Shi Y, Yang D, Zhang R, Li Z, Cui J, Zhang Y, Niu M, Li J, Wu Z, Guo H, Xiang C, Wang J, Hou J, Zhang L, Thorne RF, Cui Y, Wu Y, Wen S, Wang B. Promoter methylation-regulated miR-145-5p inhibits laryngeal squamous cell carcinoma progression by targeting FSCN1. Mol Ther. 2018;27:365–79.
Goldman M, Craft B, Swatloski T, Cline M, Morozova O, Diekhans M, Haussler D, Zhu J. The UCSC cancer genomics browser: update 2015. Nucleic Acids Res. 2015;43:D812–7.
Gao J, Aksoy BA, Dogrusoz U, Dresdner G, Gross B, Sumer SO, Sun Y, Jacobsen A, Sinha R, Larsson E, Cerami E, Sander C, Schultz N. Integrative analysis of complex cancer genomics and clinical profiles using the cBioPortal. Sci Signal. 2013;6:pl1.
Huang DW, Sherman BT, Tan Q, Kir J, Liu D, Bryant D, Guo Y, Stephens R, Baseler MW, Lane HC, Lempicki RA. DAVID bioinformatics resources: expanded annotation database and novel algorithms to better extract biology from large gene lists. Nucleic Acids Res. 2007;35:W169–75.
Vastrik I, D’Eustachio P, Schmidt E, Gopinath G, Croft D, de Bono B, Gillespie M, Jassal B, Lewis S, Matthews L, Wu G, Birney E, Stein L. Reactome: a knowledge base of biologic pathways and processes. Genome Biol. 2007;8:R39.
von Mering C, Huynen M, Jaeggi D, Schmidt S, Bork P, Snel B. STRING: a database of predicted functional associations between proteins. Nucleic Acids Res. 2003;31:258–61.
Zhang HM, Liu T, Liu CJ, Song S, Zhang X, Liu W, Jia H, Xue Y, Guo AY. AnimalTFDB 2.0: a resource for expression, prediction and functional study of animal transcription factors. Nucleic Acids Res. 2015;43:D76–81.
Gennarino VA, D’Angelo G, Dharmalingam G, Fernandez S, Russolillo G, Sanges R, Mutarelli M, Belcastro V, Ballabio A, Verde P, Sardiello M, Banfi S. Identification of microRNA-regulated gene networks by expression analysis of target genes. Genome Res. 2012;22:1163–72.
Vlachos IS, Kostoulas N, Vergoulis T, Georgakilas G, Reczko M, Maragkakis M, Paraskevopoulou MD, Prionidis K, Dalamagas T, Hatzigeorgiou AG. DIANA miRPath v.2.0: investigating the combinatorial effect of microRNAs in pathways. Nucleic Acids Res. 2012;40:W498–504.
Betel D, Koppal A, Agius P, Sander C, Leslie C. Comprehensive modeling of microRNA targets predicts functional non-conserved and non-canonical sites. Genome Biol. 2010;11:R90.
Luo J, Wang Z, Huang J, Yao Y, Sun Q, Wang J, Shen Y, Xu L, Ren B. HOXC13 promotes proliferation of esophageal squamous cell carcinoma via repressing transcription of CASP3. Cancer Sci. 2018;109:317–29.
Kho D, MacDonald C, Johnson R, Unsworth CP, O’Carroll SJ, du Mez E, Angel CE, Graham ES. Application of xCELLigence RTCA biosensor technology for revealing the profile and window of drug responsiveness in real time. Biosensors (Basel). 2015;5:199–222.
Vela I, Morrissey C, Zhang X, Chen S, Corey E, Strutton GM, Nelson CC, Nicol DL, Clements JA, Gardiner EM. PITX2 and non-canonical Wnt pathway interaction in metastatic prostate cancer. Clin Exp Metastasis. 2014;31:199–211.
Lozano-Velasco E, Hernandez-Torres F, Daimi H, Serra SA, Herraiz A, Hove-Madsen L, Aranega A, Franco D. Pitx2 impairs calcium handling in a dose-dependent manner by modulating Wnt signalling. Cardiovasc Res. 2016;109:55–66.
Benham-Pyle BW, Sim JY, Hart KC, Pruitt BL, Nelson WJ. Increasing beta-catenin/Wnt3A activity levels drive mechanical strain-induced cell cycle progression through mitosis. Elife. 2016;5:e19799.
Schinzari V, Timperi E, Pecora G, Palmucci F, Gallerano D, Grimaldi A, Covino DA, Guglielmo N, Melandro F, Manzi E, Sagnotta A, Lancellotti F, Sacco L, Chirletti P, Grazi GL, Rossi M, Barnaba V. Wnt3a/beta-catenin signaling conditions differentiation of partially exhausted T-effector cells in human cancers. Cancer Immunol Res. 2018;6:941–52.
Xu J, Lv W, Hu Y, Wang L, Wang Y, Cao J, Hu J. Wnt3a expression is associated with epithelial-mesenchymal transition and impacts prognosis of lung adenocarcinoma patients. J Cancer. 2017;8:2523–31.
Rizzolio S, Battistini C, Cagnoni G, Apicella M, Vella V, Giordano S, Tamagnone L. Downregulating neuropilin-2 triggers a novel mechanism enabling EGFR-dependent resistance to oncogene-targeted therapies. Cancer Res. 2018;78:1058–68.
Comoglio PM, Trusolino L, Boccaccio C. Known and novel roles of the MET oncogene in cancer: a coherent approach to targeted therapy. Nat Rev Cancer. 2018;18:341–58.
Meng J, Chen S, Han JX, Tan Q, Wang XR, Wang HZ, Zhong WL, Qin Y, Qiao KL, Zhang C, Gao WF, Lei YY, Liu HJ, Liu YR, Zhou HG, Sun T, Yang C. Derepression of co-silenced tumor suppressor genes by nanoparticle-loaded circular ssDNA reduces tumor malignancy. Sci Transl Med. 2018;10:6321.
Imielinski M, Berger AH, Hammerman PS, Hernandez B, Pugh TJ, Hodis E, Cho J, Suh J, Capelletti M, Sivachenko A, Sougnez C, Auclair D, Lawrence MS, Stojanov P, Cibulskis K, Choi K, de Waal L, Sharifnia T, Brooks A, Greulich H, Banerji S, Zander T, Seidel D, Leenders F, Ansen S, Ludwig C, Engel-Riedel W, Stoelben E, Wolf J, Goparju C, Thompson K, Winckler W, Kwiatkowski D, Johnson BE, Janne PA, Miller VA, Pao W, Travis WD, Pass HI, Gabriel SB, Lander ES, Thomas RK, Garraway LA, Getz G, Meyerson M. Mapping the hallmarks of lung adenocarcinoma with massively parallel sequencing. Cell. 2012;150:1107–20.
Weir BA, Woo MS, Getz G, Perner S, Ding L, Beroukhim R, Lin WM, Province MA, Kraja A, Johnson LA, Shah K, Sato M, Thomas RK, Barletta JA, Borecki IB, Broderick S, Chang AC, Chiang DY, Chirieac LR, Cho J, Fujii Y, Gazdar AF, Giordano T, Greulich H, Hanna M, Johnson BE, Kris MG, Lash A, Lin L, Lindeman N, Mardis ER, McPherson JD, Minna JD, Morgan MB, Nadel M, Orringer MB, Osborne JR, Ozenberger B, Ramos AH, Robinson J, Roth JA, Rusch V, Sasaki H, Shepherd F, Sougnez C, Spitz MR, Tsao MS, Twomey D, Verhaak RG, Weinstock GM, Wheeler DA, Winckler W, Yoshizawa A, Yu S, Zakowski MF, Zhang Q, Beer DG, Wistuba II, Watson MA, Garraway LA, Ladanyi M, Travis WD, Pao W, Rubin MA, Gabriel SB, Gibbs RA, Varmus HE, Wilson RK, Lander ES, Meyerson M. Characterizing the cancer genome in lung adenocarcinoma. Nature. 2007;450:893–8.
Seo JS, Ju YS, Lee WC, Shin JY, Lee JK, Bleazard T, Lee J, Jung YJ, Kim JO, Shin JY, Yu SB, Kim J, Lee ER, Kang CH, Park IK, Rhee H, Lee SH, Kim JI, Kang JH, Kim YT. The transcriptional landscape and mutational profile of lung adenocarcinoma. Genome Res. 2012;22:2109–19.
Quentien MH, Barlier A, Franc JL, Pellegrini I, Brue T, Enjalbert A. Pituitary transcription factors: from congenital deficiencies to gene therapy. J Neuroendocrinol. 2006;18:633–42.
Gage PJ, Suh H, Camper SA. Dosage requirement of Pitx2 for development of multiple organs. Development. 1999;126:4643–51.
Lin CR, Kioussi C, O’Connell S, Briata P, Szeto D, Liu F, Izpisua-Belmonte JC, Rosenfeld MG. Pitx2 regulates lung asymmetry, cardiac positioning and pituitary and tooth morphogenesis. Nature. 1999;401:279–82.
Quentien MH, Manfroid I, Moncet D, Gunz G, Muller M, Grino M, Enjalbert A, Pellegrini I. Pitx factors are involved in basal and hormone-regulated activity of the human prolactin promoter. J Biol Chem. 2002;277:44408–16.
Xu YY, Yu HR, Sun JY, Zhao Z, Li S, Zhang XF, Liao ZX, Cui MK, Li J, Li C, Zhang Q. Upregulation of PITX2 promotes letrozole resistance via transcriptional activation of IFITM1 signaling in breast cancer cells. Cancer Res Treat. 2018. https://doi.org/10.4143/crt.2018.100.
Semaan A, Uhl B, Branchi V, Lingohr P, Bootz F, Kristiansen G, Kalff JC, Matthaei H, Pantelis D, Dietrich D. Significance of PITX2 promoter methylation in colorectal carcinoma prognosis. Clin Colorectal Cancer. 2018;17:e385–93.
Sailer V, Gevensleben H, Dietrich J, Goltz D, Kristiansen G, Bootz F, Dietrich D. Clinical performance validation of PITX2 DNA methylation as prognostic biomarker in patients with head and neck squamous cell carcinoma. PLoS ONE. 2017;12:e0179412.
Lopez JI, Angulo JC, Martin A, Sanchez-Chapado M, Gonzalez-Corpas A, Colas B, Ropero S. A DNA hypermethylation profile reveals new potential biomarkers for the evaluation of prognosis in urothelial bladder cancer. APMIS. 2017;125:787–96.
Uhl B, Gevensleben H, Tolkach Y, Sailer V, Majores M, Jung M, Meller S, Stein J, Ellinger J, Dietrich D, Kristiansen G. PITX2 DNA methylation as biomarker for individualized risk assessment of prostate cancer in core biopsies. J Mol Diagn. 2017;19:107–14.
Xhabija B, Kidder BL. KDM5B is a master regulator of the H3K4-methylome in stem cells, development and cancer. Semin Cancer Biol. 2018;1:1. https://doi.org/10.1016/j.semcancer.2018.11.001.
Liu Y, Li Y, Hou R, Shu Z. Knockdown GREM1 suppresses cell growth, angiogenesis, and epithelial–mesenchymal transition in colon cancer. J Cell Biochem. 2018;120:5583.
Hong D, Liu T, Huang W, Liao Y, Wang L, Zhang Z, Chen H, Zhang X, Xiang Q. Gremlin1 delivered by mesenchymal stromal cells promoted epithelial–mesenchymal transition in human esophageal squamous cell carcinoma. Cell Physiol Biochem. 2018;47:1785–99.
Honma R, Sakamoto N, Ishikawa A, Taniyama D, Fukada K, Hattori T, Sentani K, Oue N, Tanabe K, Ohdan H, Yasui W. Clinicopathological and prognostic significance of epithelial gremlin1 expression in gastric cancer. Anticancer Res. 2018;38:1419–25.
Spiegelman VS, Slaga TJ, Pagano M, Minamoto T, Ronai Z, Fuchs SY. Wnt/beta-catenin signaling induces the expression and activity of betaTrCP ubiquitin ligase receptor. Mol Cell. 2000;5:877–82.
Clevers H. Wnt/beta-catenin signaling in development and disease. Cell. 2006;127:469–80.
Tenbaum SP, Ordonez-Moran P, Puig I, Chicote I, Arques O, Landolfi S, Fernandez Y, Herance JR, Gispert JD, Mendizabal L, Aguilar S, Ramon y Cajal S, Schwartz S Jr, Vivancos A, Espin E, Rojas S, Baselga J, Tabernero J, Munoz A, Palmer HG. Beta-catenin confers resistance to PI3K and AKT inhibitors and subverts FOXO3a to promote metastasis in colon cancer. Nat Med. 2012;18:892–901.
Urakami S, Shiina H, Enokida H, Kawakami T, Tokizane T, Ogishima T, Tanaka Y, Li LC, Ribeiro-Filho LA, Terashima M, Kikuno N, Adachi H, Yoneda T, Kishi H, Shigeno K, Konety BR, Igawa M, Dahiya R. Epigenetic inactivation of Wnt inhibitory factor-1 plays an important role in bladder cancer through aberrant canonical Wnt/beta-catenin signaling pathway. Clin Cancer Res. 2006;12:383–91.
Chim CS, Pang R, Fung TK, Choi CL, Liang R. Epigenetic dysregulation of Wnt signaling pathway in multiple myeloma. Leukemia. 2007;21:2527–36.
Hosseini V, Dani C, Geranmayeh MH, Mohammadzadeh F, Nazari S, Ahmad S, Darabi M. Wnt lipidation: roles in trafficking, modulation, and function. J Cell Physiol. 2018;234:8040–54.
Zhang S, Chen X, Hu Y, Wu J, Cao Q, Chen S, Gao Y. All-trans retinoic acid modulates Wnt3A-induced osteogenic differentiation of mesenchymal stem cells via activating the PI3K/AKT/GSK3beta signalling pathway. Mol Cell Endocrinol. 2016;422:243–53.
Garzon R, Marcucci G, Croce CM. Targeting microRNAs in cancer: rationale, strategies and challenges. Nat Rev Drug Discov. 2010;9:775–89.
Cha Y, He Y, Ouyang K, Xiong H, Li J, Yuan X. MicroRNA-140-5p suppresses cell proliferation and invasion in gastric cancer by targeting WNT1 in the WNT/beta-catenin signaling pathway. Oncol Lett. 2018;16:6369–76.
Wu D, Zhang J, Lu Y, Bo S, Li L, Wang L, Zhang Q, Mao J. miR-140-5p inhibits the proliferation and enhances the efficacy of doxorubicin to breast cancer stem cells by targeting Wnt1. Cancer Gene Ther. 2018;26:74–82.
Han XR, Wen X, Wang YJ, Wang S, Shen M, Zhang ZF, Fan SH, Shan Q, Wang L, Li MQ, Hu B, Sun CH, Wu DM, Lu J, Zheng YL. MicroRNA-140-5p elevates cerebral protection of dexmedetomidine against hypoxic-ischaemic brain damage via the Wnt/beta-catenin signalling pathway. J Cell Mol Med. 2018;22:3167–82.
YS, YQ and QD designed and supervised the study. JL performed the experiments in vitro and wrote the manuscript. YY performed the experiments in vivo. SJ was responsible for statistical analysis. QS and KL provided technical support and helped to refine the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
This work was supported by the National Natural Science Foundation of China (No. 81172032).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.