B-Myb Mediates Proliferation and Migration of Non-Small-Cell Lung Cancer via Suppressing IGFBP3

B-Myb has been shown to play an important oncogenic role in several types of human cancers, including non-small-cell lung cancer (NSCLC). We previously found that B-Myb is aberrantly upregulated in NSCLC, and overexpression of B-Myb can significantly promote NSCLC cell growth and motility. In the present study, we have further investigated the therapeutic potential of B-Myb in NSCLC. Kaplan–Meier and Cox proportional hazards analysis indicated that high expression of B-Myb is significantly associated with poor prognosis in NSCLC patients. A loss-of-function study demonstrated that depletion of B-Myb resulted in significant inhibition of cell growth and delayed cell cycle progression in NSCLC cells. Notably, B-Myb depletion also decreased NSCLC cell migration and invasion ability as well as colony-forming ability. Moreover, an in vivo study demonstrated that B-Myb depletion caused significant inhibition of tumor growth in a NSCLC xenograft nude mouse model. A molecular mechanistic study by RNA-seq analysis revealed that B-Myb depletion led to deregulation of various downstream genes, including insulin-like growth factor binding protein 3 (IGFBP3). Overexpression of IGFBP3 suppressed the B-Myb-induced proliferation and migration, whereas knockdown of IGFBP3 significantly rescued the inhibited cell proliferation and motility caused by B-Myb siRNA (small interfering RNA). Expression and luciferase reporter assays revealed that B-Myb could directly suppress the expression of IGFBP3. Taken together, our results suggest that B-Myb functions as a tumor-promoting gene via suppressing IGFBP3 and could serve as a novel therapeutic target in NSCLC.


Introduction
Lung cancer is the most frequent cause of cancer-related deaths for patients in Western society; the five-year relative survival rate for lung cancer is around 18% [1]. Moreover, many studies have reported that lung cancer has become a major public health challenge as the incidence and mortality of lung cancer rapidly increases in developed countries, especially in China [2]. Non-small-cell lung cancer (NSCLC) comprises the majority of diagnosed lung cancer types, which is associated with a relatively poor 15% overall five-year survival rate. Despite there being many reports describing the therapy of NSCLC, the prognosis of the disease remains unfavorable. Thus, developing novel therapeutic strategies and innovative insights for NSCLC is necessary.

Prognostic Significance of B-Myb in NSCLC (Non-Small-Cell Lung Cancer)
We previously found that B-Myb was upregulated in non-small-cell lung cancer (NSCLC). Moreover, the expression level of B-Myb was positively associated with tumor pathological grade, clinical stage, tumor, node, metastasis (TNM) classification, and lymph node metastasis [20]. In this study, we aimed at evaluating the prognostic value of B-Myb expression in silico using published microarray data, which are from primary tumors of lung cancer patients with follow-up information. As shown in Figure 1A, high expression of B-Myb was strongly associated with poor overall survival of lung cancer patients in the Nagoya University lung adenocarcinoma (ADC) cohort and Michigan University lung squamous cell carcinoma (SQCC) cohort (p < 0.05). In consistency with this observation, as shown in Figure 1B,C, online Kaplan-Meier plotter analysis [29] also revealed that B-Myb overexpression was negatively related to significant improvement in patient survival rates in lung ADC and SQCC. Moreover, univariate analyses revealed that high B-Myb expression was significantly associated with poorer survival in both cohorts (hazard ratio (HR) = 1.870, 95% confidence interval (CI) = 1.024-3.416, p = 0.042). Multivariate Cox regression analysis displayed that B-Myb expression was an independent prognostic factor for the Nagoya University cohort (HR = 1.789, 95% CI = 0.974-3.286, p = 0.043). In addition, lymph node metastasis was significantly related to poorer survival (p = 0.003) and the independent prognostic factor (p = 0.002) for the Nagoya University cohort (Table 1).

B-Myb Depletion Delays the Cell Cycle Progression and Inhibits Proliferation in Adenocarcinoma Cells (ADC)
To investigate the therapeutic potential of B-Myb in NSCLC, we depleted the B-Myb expression via small interfering RNA (siRNA)-mediated silencing in A549 lung cancer cell lines, and cell proliferation and cell cycle assays were subsequently performed. Quantitative RT-PCR and Western blot analysis showed that the B-Myb expression was significantly suppressed at both the mRNA and protein levels in A549 lung cancer cell lines (Figure 2A). B-Myb depletion resulted in a significant growth retardation compared with control siRNA from a later time point (96 h) in A549 cells ( Figure 2B). Cell cycle analysis revealed that silencing B-Myb expression caused a remarkable G1 arrest in A549 cells ( Figure 2C). Moreover, our previous study [20] showed that B-Myb depletion affects the cell cycle progression and inhibits proliferation in H1299 cells. These results suggested that B-Myb depletion mainly delays cell cycle progression and significantly inhibits proliferation in both A549 and H1299 cells.

B-Myb Depletion Delays the Cell Cycle Progression and Inhibits Proliferation in Adenocarcinoma Cells (ADC)
To investigate the therapeutic potential of B-Myb in NSCLC, we depleted the B-Myb expression via small interfering RNA (siRNA)-mediated silencing in A549 lung cancer cell lines, and cell proliferation and cell cycle assays were subsequently performed. Quantitative RT-PCR and Western blot analysis showed that the B-Myb expression was significantly suppressed at both the mRNA and protein levels in A549 lung cancer cell lines (Figure 2A). B-Myb depletion resulted in a significant growth retardation compared with control siRNA from a later time point (96 h) in A549 cells ( Figure 2B). Cell cycle analysis revealed that silencing B-Myb expression caused a remarkable G1 arrest in A549 cells ( Figure 2C). Moreover, our previous study [20] showed that B-Myb depletion affects the cell cycle progression and inhibits proliferation in H1299 cells. These results suggested that B-Myb depletion mainly delays cell cycle progression and significantly inhibits proliferation in both A549 and H1299 cells. Figure 2. B-Myb depletion affects cell cycle progression and inhibits proliferation in A549 lung cancer cells. (A) A549 cells of small interfering RNA (siRNA)-mediated B-Myb silencing were transiently transfected with the negative control (NCsi) and B-Myb siRNA (B-Mybsi), respectively. Forty-eight and seventy-two hours after transfection, total RNA and whole cell lysates were respectively prepared and subjected to qRT-PCR and Western blot, and glyceraldehyde-phosphate dehydrogenase GAPDH as control proteins. (B) B-Myb depletion reduced cell proliferation. A549 cells were transiently transfected with negative control or B-Myb siRNA, and cell proliferation was then determined by cell counting kit-8 assay kits (CCK8) at the indicated time points. (C) B-Myb depletion delays G1-S phase transition. A549 cells were seeded on six-well plates and transfected with the indicated siRNAs, and twenty-four hours later, cells were collected and subjected to cell cycle analysis. All experiments were performed in triplicates. Data represent the mean ± standard deviation (SD). * p < 0.05, ** p < 0.01, *** p < 0.001.

B-Myb Depletion Reduces Motility in A549 Lung Cancer Cells
Furthermore, we asked whether depletion of B-Myb expression regulates cell migration ability in A549 lung cancer cells. We performed wound healing and transwell migration assays to evaluate the cell migration ability. The wound healing assay suggested that siRNA-mediated silencing of B-Myb significantly decreased the rate of lateral migration of cells compared with the control cells ( Figure 3A). In addition, the transwell migration assay showed that B-Myb knockdown also significantly suppressed the cell transmigration ability compared with the control cells ( Figure 3B). A549 cells were seeded on six-well plates and transfected with the indicated siRNAs, and twenty-four hours later, cells were collected and subjected to cell cycle analysis. All experiments were performed in triplicates. Data represent the mean ± standard deviation (SD). * p < 0.05, ** p < 0.01, *** p < 0.001.

B-Myb Depletion Reduces Motility in A549 Lung Cancer Cells
Furthermore, we asked whether depletion of B-Myb expression regulates cell migration ability in A549 lung cancer cells. We performed wound healing and transwell migration assays to evaluate the cell migration ability. The wound healing assay suggested that siRNA-mediated silencing of B-Myb significantly decreased the rate of lateral migration of cells compared with the control cells ( Figure 3A). In addition, the transwell migration assay showed that B-Myb knockdown also significantly suppressed the cell transmigration ability compared with the control cells ( Figure 3B). Thus, these results strongly suggested that B-Myb might also regulate the cell motility of A549 lung cancer cells. Thus, these results strongly suggested that B-Myb might also regulate the cell motility of A549 lung cancer cells. Representative images (×200) (left) and quantification results (right) were shown for each assay. All experiments were performed in triplicates. Data represent the mean ± standard deviation (SD). *** p < 0.001.

B-Myb Depletion Inhibits Lung Tumorigenesis In Vitro and In Vivo
To further validate the function of B-Myb in the tumorigenesis of lung cancer, we established stable B-Myb knockdown cells in H1299 cells. Quantitative RT-PCR and Western blot results showed that the B-Myb expression was significantly decreased at both the mRNA and protein levels in the stable B-Myb knockdown cells ( Figure 4A). Next, we analyzed the effect of B-Myb knockdown on the colony formation of H1299 lung cancer cells. As shown in Figure 4B,C, colony formation assays on plastic and soft agar revealed that B-Myb knockdown in H1299 cells significantly inhibited anchorage-dependent and -independent colony-forming ability compared with the control cells.
To examine whether silencing of B-Myb expression could inhibit tumorigenesis in vivo, we used A549 cells to establish stable cells with lentiviral particle B-Myb knockdown. Quantitative RT-PCR and Western blot analyses confirmed that the B-Myb expression was significantly decreased at both the mRNA and protein levels in the stable B-Myb knockdown cells ( Figure 5A). In consistency with the results observed in the siRNA-mediated B-Myb depletion system, stable B-Myb knockdown also significantly decreased cell growth and motility in the lentivirus-mediated B-Myb depletion system, as evidenced by proliferation ( Figure 5B) and transwell migration assays ( Figure 5C). Subsequently, the stable B-Myb knockdown and the negative control cells were injected subcutaneously into the dorsal flank of nude mice, respectively. As shown in Figure 5D-F, stable knockdown of B-Myb expression remarkably suppressed the tumor growth in both tumor volume and tumor weight in the nude mice. Representative images (×200) (left) and quantification results (right) were shown for each assay. All experiments were performed in triplicates. Data represent the mean ± standard deviation (SD). *** p < 0.001.

B-Myb Depletion Inhibits Lung Tumorigenesis in Vitro and in Vivo
To further validate the function of B-Myb in the tumorigenesis of lung cancer, we established stable B-Myb knockdown cells in H1299 cells. Quantitative RT-PCR and Western blot results showed that the B-Myb expression was significantly decreased at both the mRNA and protein levels in the stable B-Myb knockdown cells ( Figure 4A). Next, we analyzed the effect of B-Myb knockdown on the colony formation of H1299 lung cancer cells. As shown in Figure 4B,C, colony formation assays on plastic and soft agar revealed that B-Myb knockdown in H1299 cells significantly inhibited anchorage-dependent and -independent colony-forming ability compared with the control cells.
To examine whether silencing of B-Myb expression could inhibit tumorigenesis in vivo, we used A549 cells to establish stable cells with lentiviral particle B-Myb knockdown. Quantitative RT-PCR and Western blot analyses confirmed that the B-Myb expression was significantly decreased at both the mRNA and protein levels in the stable B-Myb knockdown cells ( Figure 5A). In consistency with the results observed in the siRNA-mediated B-Myb depletion system, stable B-Myb knockdown also significantly decreased cell growth and motility in the lentivirus-mediated B-Myb depletion system, as evidenced by proliferation ( Figure 5B) and transwell migration assays ( Figure 5C). Subsequently, the stable B-Myb knockdown and the negative control cells were injected subcutaneously into the dorsal flank of nude mice, respectively. As shown in Figure 5D-F, stable knockdown of B-Myb expression remarkably suppressed the tumor growth in both tumor volume and tumor weight in the nude mice. cells were seeded on plastic plates for anchorage-dependent colony-formation assay (B) and soft agar for anchorage-independent colony-formation assay (C), respectively. Representative images (×200) (left) and quantification results (right) were shown for assay of B and C. All experiments were performed in triplicates. Data represent the mean ± SD. * p < 0.05, ** p < 0.01. The stable negative control (pNC-sh) and B-Myb knockdown (pB-Myb-sh) H1299 cells were seeded on plastic plates for anchorage-dependent colony-formation assay (B) and soft agar for anchorage-independent colony-formation assay (C), respectively. Representative images (×200) (left) and quantification results (right) were shown for assay of B and C. All experiments were performed in triplicates. Data represent the mean ± SD. * p < 0.05, ** p < 0.01.

B-Myb Promotes Cell Proliferation and Migration by Targeting IGFBP3 in NSCLC Cells
IGFBP3 has been shown to play an important role in proliferation, carcinogenesis, and ERK and Akt signaling pathways [23,24]. Thus, we decided to next investigate whether IGFBP3 acts as an important downstream target gene to mediate the oncogenic role of B-Myb in NSCLC cells. Firstly, we examined whether IGFBP3 overexpression could decrease the oncogenic effects of B-Myb on NSCLC cell proliferation. To this end, IGFBP3 expression plasmid and empty vector were transiently transfected to H1299 cells overexpressing B-Myb. Quantitative RT-PCR analysis confirmed the successful overexpression of IGFBP3 in cells ( Figure 7A). IGFBP3 overexpression significantly suppressed the proliferation of H1299 cells induced by B-Myb ( Figure 7B). Secondly, IGFBP3-siRNA and NC-siRNA were transiently transfected to H1299 cells that had been transfected with the B-Myb-siRNA. Quantitative RT-PCR analysis validated the successful knockdown of IGFBP3 in cells ( Figure 7C). After transfection with IGFBP3-siRNA, proliferation of H1299 cells inhibited by B-Myb knockdown was markedly increased ( Figure 7D). Moreover, in the wound healing assay, IGFBP3 overexpression inhibited the migratory capacity of H1299 cells induced by the B-Myb overexpression ( Figure 7E). These results suggest that B-Myb promotes cell proliferation and migration by targeting IGFBP3 in H1299 cells. IGFBP3 has been shown to play an important role in proliferation, carcinogenesis, and ERK and Akt signaling pathways [23,24]. Thus, we decided to next investigate whether IGFBP3 acts as an important downstream target gene to mediate the oncogenic role of B-Myb in NSCLC cells. Firstly, we examined whether IGFBP3 overexpression could decrease the oncogenic effects of B-Myb on NSCLC cell proliferation. To this end, IGFBP3 expression plasmid and empty vector were transiently transfected to H1299 cells overexpressing B-Myb. Quantitative RT-PCR analysis confirmed the successful overexpression of IGFBP3 in cells ( Figure 7A). IGFBP3 overexpression significantly suppressed the proliferation of H1299 cells induced by B-Myb ( Figure 7B). Secondly, IGFBP3-siRNA and NC-siRNA were transiently transfected to H1299 cells that had been transfected with the B-Myb-siRNA. Quantitative RT-PCR analysis validated the successful knockdown of IGFBP3 in cells ( Figure 7C). After transfection with IGFBP3-siRNA, proliferation of H1299 cells inhibited by B-Myb knockdown was markedly increased ( Figure 7D). Moreover, in the wound healing assay, IGFBP3 overexpression inhibited the migratory capacity of H1299 cells induced by the B-Myb overexpression ( Figure 7E). These results suggest that B-Myb promotes cell proliferation and migration by targeting IGFBP3 in H1299 cells.

B-Myb Activates ERK and Akt Pathways via Targeting IGFBP3 in NSCLC Cells
Our previous report showed that B-Myb overexpression could activate ERK and Akt pathways in NSCLC cells [20]. Therefore, we further investigated whether B-Myb regulates ERK and Akt pathways via IGFBP3. As shown in Figure 8, IGFBP3 overexpression significantly decreased the levels of phosphorylated ERK and Akt in the stable B-Myb expression cells. On the other hand, IGFBP3 knockdown increased the levels of phosphorylated ERK and Akt, whereas B-Myb knockdown decreased the phosphorylated ERK and Akt. Taken together, these results strongly suggest that B-Myb activates ERK and Akt signaling at least partially through inhibiting IGFBP3.
in NSCLC cells [20]. Therefore, we further investigated whether B-Myb regulates ERK and Akt pathways via IGFBP3. As shown in Figure 8, IGFBP3 overexpression significantly decreased the levels of phosphorylated ERK and Akt in the stable B-Myb expression cells. On the other hand, IGFBP3 knockdown increased the levels of phosphorylated ERK and Akt, whereas B-Myb knockdown decreased the phosphorylated ERK and Akt. Taken together, these results strongly suggest that B-Myb activates ERK and AKT signaling at least partially through inhibiting IGFBP3.

B-Myb Inhibits IGFBP3 Promoter Activity
Finally, we asked whether B-Myb could directly regulate IGFBP3 in NSCLC cells. To this end, we firstly analyzed the putative transcription factor binding sites in the IGFBP3 promoter region, and found the existence of potential MYB-binding sites (t/cAACt/gG) in the IGFBP3 promoter ( Figure 9A). To further verify whether B-Myb could directly suppress IGFBP3 gene transcription, H1299 cells were cotransfected with the B-Myb expression plasmid together with the IGFBP3-P2822 luciferase reporter plasmid. The luciferase reporter assay revealed that overexpression of B-Myb inhibits luciferase activities driven by the IGFBP3-P2822 promoter ( Figure 9B). In a sharp contrast to B-Myb overexpression, B-Myb depletion enhances luciferase activities driven by the IGFBP3-P2822 promoter ( Figure 9C). Collectively, our results strongly suggest that B-Myb directly suppresses IGFBP3 gene transcription.
cotransfected with the B-Myb expression plasmid together with the IGFBP3-P2822 luciferase reporter plasmid. The luciferase reporter assay revealed that overexpression of B-Myb inhibits luciferase activities driven by the IGFBP3-P2822 promoter ( Figure 9B). In a sharp contrast to B-Myb overexpression, B-Myb depletion enhances luciferase activities driven by the IGFBP3-P2822 promoter ( Figure 9C). Collectively, our results strongly suggest that B-Myb directly suppresses IGFBP3 gene transcription.

Discussion
It has been well established that B-Myb expression is highly expressed in many types of cancers [13,15,30,31]. Additionally, B-Myb overexpression shows poor prognosis in patients with neuroblastoma [18], acute myeloid leukaemia [19], and breast cancer [15]. Tao et al. verified the B-Myb overexpression in breast cancer, particularly in advanced cancer, and demonstrated that B-Myb overexpression also predicts shorter overall survival of breast cancer patients [32]. Raschella et al. reported that B-Myb is overexpressed in neuroblastoma tumors, and that high expression level of B-Myb is an independent prognostic factor for poor neuroblastoma tumor patient survival [18]. Ren et al. confirmed that B-Myb expression is significantly increased in colorectal cancer (CRC) tissues compared to adjacent noncancerous tissues; that B-Myb overexpression is an independent prognostic factor for worse survival in colorectal cancer patients, as evidenced by Cox multivariate analysis; and that functional analysis showed that downregulation of B-Myb expression suppressed cell proliferation and motility, restricted cell cycle progression, and induced apoptosis [13]. Our previous studies verified that B-Myb is upregulated in non-small-cell lung cancer (NSCLC), and a gain-of-function study has verified that overexpression of B-Myb significantly promotes lung cancer cell growth, migration, and invasion [20]. In the present study, we firstly found prognostic significance of B-Myb in NSCLC, and a loss-of-function study demonstrated that downregulation of B-Myb significantly suppressed NSCLC cell growth, cell cycle progression, colony-forming ability, motility, and tumorigenesis in vivo. These results suggested that B-Myb is required for the accelerated proliferation and motility in NSCLC. Thus, the results of our study in NSCLC are consistent with the aforementioned reports obtained from other types of cancers. Moreover, the prognostic value of B-Myb expression was verified in NSCLC.
Previous studies have also indicated that the product of the B-Myb gene is a transcription factor that regulates the proliferation and progression of malignancies through suppression of the cell cycle and activation of genes and pathways corelated to tumorigenesis [3,4,8,33]. Aberrant B-Myb expression or amplification has been confirmed in different types of human cancers, verifying a role in tumorigenesis [31,34,35]. According to our recent observations [20], exogenous B-Myb overexpression in NSCLC promoted G1-S phase transition progression; caused transactivation of multiple downstream genes such as CCNA1, COL11A1, COL6A1, MMP2, NID1, and FLT4; and activated ERK and Akt signaling pathways. Consistently, our results in this study showed that knockdown of B-Myb caused a significant decrease in the number of S-phase cells. Moreover, RNA-seq analysis showed that knockdown of B-Myb leads to upregulation of different downstream genes such as IGFBP3 and downregulation of various downstream genes such as COL11A1, FLT4, SPARC, IDH2, and PDK3. Previous reports showed that SPARC promotes cancer development in some cancers with highly metastatic characteristics, including melanoma and breast cancer [36,37], and was uncovered as the potential prognostic marker in several cancer types [38][39][40]. IDH2 and PDK3 have been shown to positively regulate the proliferation and tumorigenesis of various cancer cells [41][42][43]. It is worthwhile to mention that IGFBP3 overexpression suppresses the growth of NSCLC cells in vitro and in vivo, and leads to apoptosis in NSCLC cells by suppression of the signal transduction mechanism participating in cellular proliferation and survival [23]. It has been revealed that Akt and MAPK have significant inhibitory effects on apoptosis [44], and the function of Akt in survival has been confirmed in a study in which cells were exposed to different apoptotic stimuli [45][46][47]. Lee et al. reported that IGFBP3 overexpression inhibits the MAPK and Akt pathways, which have crucial roles in cell proliferation and survival [23]. Moreover, we also found that the differentially expressed genes (DEGs) affected by B-Myb depletion could be enriched in MAPK signaling pathways. Of note, our previous studies confirmed that B-Myb overexpression could enhance the activation of ERK and Akt signaling pathways in NSCLC [20]. Furthermore, our studies demonstrated that the expression level of IGFBP3 mRNA was negatively regulated by B-Myb expression, and that IGFBP3 negatively regulated the proliferation and migration of NSCLC induced by B-Myb differential expression. Furthermore, B-Myb depletion could suppress the activation of ERK and Akt signaling pathways though upregulation of IGFBP3 in NSCLC, and B-Myb overexpression could promote the activation of ERK and Akt signaling pathways through inhibition of IGFBP3. In conclusion, our present study demonstrated for the first time that B-Myb is an independent prognostic marker and serves as a potential target in the diagnosis and/or treatment of NSCLC, and that B-Myb functions as a tumor-promoting gene by targeting IGFBP3 in NSCLC cells.

siRNA Synthesis and Transfection
Small interfering RNAs against B-Myb (B-Mybsi) as well as the negative control (NCsi) were chemically synthesized by GenePharma (Shanghai, China). The sequences of the B-Myb siRNA are as follows: 5 -CAGACAAUGCUGUGAAGAATT-3 (sense) and 5 -UUCUUCACAGCAUUGUCUGTT-3 (antisense). The sequences of NC siRNA are as follows: 5 -UUCUCCGAACGUGUCACGUTT-3 (sense) and 5 -ACGUGACACGUUCGGAGAATT-3 (antisense). Knockdown of B-Myb expression was conducted using Lipofectamine RNAiMAX reagent (Invitrogen) according to the manufacturer's instructions. Cells were collected and subjected to subsequent analysis 48 to 72 h after transfection.

RNA Isolation and qRT-PCR
Total RNA was extracted using TRIzol reagent (Invitrogen) according to the manufacturer's instructions. Reverse transcription of total RNA (500 ng) was performed for obtaining cDNA using PrimeScript ® RT reagent Kit (TaKaRa, Otsu, Japan) according to the manufacturer's instructions. qRT-PCR was performed using SYBR Premix Ex Taq™ (Perfect Real Time, TaKaRa). The primers sequences were provided in Table S2. Each experiment was repeated three times.

Western Blot Analysis
Cells were lysed in RIPA buffer including 1% protease inhibitor PMSF (Beyotime, Jiangsu, China). The proteins concentration was quantified by the bicinchoninic acid (BCA) protein assay kit (Thermo Scientific, Waltham, MA, USA). 40 µg of each protein was subjected to 8% SDS-PAGE, following transferral to the polyvinylidene fluoride (PVDF) membrane (Millipore, Burlington, MA, USA). The information of the antibodies used for Western blotting was provided in Table S3. The blots were visualized by the enhanced chemiluminescence (ECL; Bio-Rad Laboratories, Hercules, CA, USA).

Cell Proliferation Assay
For the cell proliferation assay, 10 µL of cell counting kit-8 (CCK8; Dojindo, Tokyo, Japan) was added to the cells seeded in 96-well plates at a density of 0.2 × 10 4 cells per well. The absorbance value (OD) of each well was measured at 450 nm at 0, 24, 48, 72, and 96 h after applying the CCK8 reagent according to the manufacturer's instructions. The experiments were repeated three times.

Cell Cycle Assay
For the cell cycle assay, cells were collected by trypsin digestion and cell cycle assays were performed using the cell cycle kit (Key GEN, Nanjing, China) according to the manufacturer's instructions. Cell cycle distribution was analyzed on a FACScan flow cytometer (BD Biosciences, San Jose, CA, USA).

Colony-Formation Assay
The anchorage-independent soft agar assay was conducted as described previously with minor modifications [33]. Anchorage-dependent colony-formation assay was performed in 6-well plates with a cell density of 2000 cells per well, and continuously cultured for ten days. Colonies were then fixed with 4% paraformaldehyde for 15 min and visualized by staining for 10 min with 0.5% crystal violet.

Wound Healing Assay
The cells were seeded in 6-well plates and cultured. When the cells reached 80% confluence, artificial wounds were firmly scraped with a 10-µL pipette tip. Subsequently floating cells and debris were removed by washing the cells with PBS. The cells were then incubated in medium supplemented with 3% fetal bovine serum (Invitrogen), penicillin (100 IU/mL), and streptomycin (100 mg/mL). Initial gap widths (0 h) and residual gap widths at 24 h and 48 h after wounding were captured in the same wounded region under an inverted microscope (Olympus, Hamburg, Germany).

Transwell Cell Migration Assay
For the transwell cell migration assay, A549 cells were seeded in the upper chamber of a transwell device (8-µm, Merck Millipore, Burlington, MA, USA) at a density of 2 × 10 4 cells per well in serum-free medium, and 700 µL of medium containing 10% FBS was added to the bottom chamber. Thirty-six hours later, cells migrating into the lower surface of the chambers were fixed with methyl alcohol for 20 min and stained with 0.5% crystal violet for 10 min. Then, the cells on the upper surface of the filter were removed using a cotton swab. Five fields were imaged per transwell insert and were reviewed and photographed in five random fields, and the number of cells was counted using the particle-counting module of Photoshop CS5 (Adobe, San Jose, CA, USA).

Transfection
Empty vector plasmids (pSuper-retro-sh) had been preserved in our laboratories, and short hairpin RNA (shRNA) was synthesized by Sangon Biotech (Shanghai, China) according to the sequence of control-siRNA or B-Myb-siRNA. The sequences of the siRNA used were listed in Materials and Methods, Section 4.2. Annealing shRNA was inserted into pSuper-retro-sh via enzyme digestion. Through transformation, selection, and plasmid extraction, the plasmids of negative control and silencing B-Myb gene were constructed well, and were named pControl-sh and pB-Myb-sh. H1299 cells seeded into 6-well plates were transiently transfected with 2 µg of these plasmids using Lipofectamine 2000 reagent (Invitrogen) according to the manufacturer's instructions. After 6 h, the cells were seeded into 96-well plates at a density of about one cell per well and selected in the presence of 1.3 µg/mL puromycin for generating the control and stably silencing B-Myb clones. The survived puromycin-resistant cell clones were then continuously expanded and confirmed to construct the corresponding stable cell lines.
The IGFBP3 overexpression plasmid was a gift of Dr. Hiroshi Nakagawa (PA, USA) [48], and its control vector plasmid was preserved in our laboratories. Short hair IGFBP3 shRNA synthesized by Sangon Biotech (Shanghai, China) was inserted into the psi-LVRH1GP via enzyme digestion, and the IGFBP3 shRNA plasmid was constructed through transformation, selection, and plasmid extraction. The target sequence for IGFBP3 shRNA is as follows: 5 -GCACAGATACCCAGAACTT-3 . H1299 cells of B-Myb overexpression or B-Myb knockdown were transfected with corresponding plasmids to overexpress or knock down IGFBP3 using Lipofectamine 2000 reagent according to the manufacturer's instructions, and incubated at 37 • C for 48 h.

Luciferase Reporter Constructs and Reporter Assays
The human IGFBP3 and (−2282 nucleotide (nt) to +56 nt)-luciferase promoters were gifted by Dr. T. Hanafusa (Okayama University, Okayama, Japan) [49]. Luciferase reporter constructs containing the potential MYB-binding sites were also generated using DNAMAN. All of the reporter constructs were validated by direct sequencing. For luciferase reporter analysis, cells were seeded in triplicate into 12-well plates and cotransfected with the indicated reporter vectors, pRL-TK vector (Promega) encoding Renilla luciferase together with the other vector, using Lipofectamine 2000 reagent (Invitrogen). Forty-eight hours later, cells were lysed, and their luciferase activities were measured using the dual luciferase assay system (Promega) as described previously.

Lentiviral Infection and Construction of Stably Silencing B-Myb Gene A549 Cell Lines
Lentiviral vectors (LV-B-Myb-RNAi and LV-control-RNAi) were obtained from GeneChem (Shanghai, China) and respectively infected into A549 cells according to the manufacturer's instructions. Infected derivative cells stably expressing shRNA were selected in the presence of 0.5 µg/mL puromycin for 10 days, then named B-Myb-sh and NC-sh.

Tumor Xenografts
All procedures conformed to the legal mandates and guidelines of the Laboratory Animal Center of Chongqing Medical University for the care and maintenance of laboratory animals. 3 × 10 7 stable control or B-Myb silencing cells were subcutaneously injected into the dorsal flank of fourto six-week old nude mice (BALB/c). Six mice were used for each group. The xenograft size was measured every four days by a Vernier caliper along two perpendicular axes. The volume of the xenograft was calculated following the formula: volume = 1/2 × length × width 2 . Eight weeks after injection, the mice were killed and the tumor specimens were weighed. All experimental procedures involving animals were conducted based on laboratory animal protocols approved at 3 May 2016 by the Laboratory Animal Center of Chongqing Medical University (No. 20160503).

RNA-Seq Analysis
RNA-seq analysis was performed as described previously [20].

Prognostic Analysis of B-Myb in Lung Cancer Patients
Prognostic analysis of B-Myb in lung cancer patients were performed as described previously [49]. The Nagoya lung adenocarcinoma (ADC) and Michigan lung squamous cell carcinoma (SQCC) datasets were downloaded from the public gene expression omnibus (GEO) database (GSE4573 and GSE13213 [29]. In addition, the prognostic value of B-Myb expression in lung ADC and SQCC was also evaluated by an online tool: Kaplan-Meier plotter analysis [29]].

Lentivirus-Mediated Establishment of Stable B-Myb Overexpression Cell Line
According to the method described previously [20], the lentivirus-mediated establishment of a stable B-Myb overexpression cell line was accomplished.

Statistical Analysis
All of the statistical analyses were performed using the SPSS16.0 software (SPSS INC, Chicago, IL, USA). Statistical significance among different groups in the experiments was tested by ANOVA and unpaired t-test. Overall survival curves were estimated by the Kaplan-Meier method and compared among groups by log-rank test. The Cox proportional hazards regression model was applied to perform univariate and multivariate regression analyses, which calculate the hazard ratios (HRs) and their 95% confidence intervals (95% CIs) for each variable and analyze independent factors affecting prognosis. All comparisons were carried out in two-tailed t-tests, and p values of <0.05 were taken to be statistically significant.