Targeting multiple genes containing long mononucleotide A-T repeats in lung cancer stem cells
Journal of Translational Medicine volume 19, Article number: 231 (2021)
Intratumour heterogeneous gene expression among cancer and cancer stem cells (CSCs) can cause failure of current targeted therapies because each drug aims to target the function of a single gene. Long mononucleotide A-T repeats are cis-regulatory transcriptional elements that control many genes, increasing the expression of numerous genes in various cancers, including lung cancer. Therefore, targeting A-T repeats may dysregulate many genes driving cancer development. Here, we tested a peptide nucleic acid (PNA) oligo containing a long A-repeat sequence [A(15)] to disrupt the transcriptional control of the A-T repeat in lung cancer and CSCs.
First, we separated CSCs from parental lung cancer cell lines. Then, we evaluated the role of A-T repeat gene regulation by counting the number of repeats in differentially regulated genes between CSCs and the parental cells of the CSCs. After testing the dosage and effect of PNA-A15 on normal and cancer cell toxicity and CSC phenotypes, we analysed genome-wide expression to identify dysregulated genes in CSCs.
The number of A-T repeats in genes differentially regulated between CSCs and parental cells differed. PNA-A15 was toxic to lung cancer cells and CSCs but not to noncancer cells. Finally, PNA-A15 dysregulated a number of genes in lung CSCs.
PNA-A15 is a promising novel targeted therapy agent that targets the transcriptional control activity of multiple genes in lung CSCs.
Cancer stem cells (CSCs) are the major reason for the devastating effects of cancer . CSCs are critical factors that drive cancer aggressiveness, metastasis and drug resistance [1, 2]. Thus, CSC-targeted cancer therapy is crucial. Thus, developing new drugs that can effectively kill the cell population is not only a promising solution but also very challenging.
Intratumour heterogeneous gene expression among cancer and CSCs can cause a failure of current targeted therapy because each drug aims to target the function of a single gene . Currently, CSCs can acquire resistance to EGFR (Epidermal growth factor receptor) inhibitors, which is an effective therapy for lung adenocarcinoma patients with activating mutations, making effective therapy more difficult . Therefore, a novel CSC-targeted therapy is required.
Long mononucleotide A-T repeats are cis-regulatory transcriptional controls of many genes, and the element increases the expression of many genes in various cancers, including lung cancer [5, 6]. A-T repeats are cis-regulatory elements that regulate gene expression by binding Argonaute proteins (Agos) . Therefore, we hypothesized that disruption of A-T repeat interactions with Agos would dysregulate many cancer and CSC genes and that targeting A-T repeats may dysregulate many genes driving cancer development and activate transcription in several ways . Thus, the binding of A and T repeats could be a targeted therapy for cancer treatment. The transcriptional control of the A-T repeat can be disrupted by a peptide nucleic acid (PNA) oligo containing a long A-repeat sequence[A(15)] . Here, we demonstrated that PNA-A15 has promising potential in lung cancer and CSC targeted therapy by dysregulating multiple cancer and CSC genes.
Materials and methods
Cell culture and treatment
Human NSCLC-derived cell lines H460, H23, H292, and A549; human kidney cells HK2 and HEK293; and the periodontal ligament (PDL) fibroblast cells were obtained from the American Type Culture Collection (ATCC®, Manassas, VA, USA). H460, H23, and H292 cells were cultivated in Roswell Park Memorial Institute (RPMI) 1640 medium supplemented with 10% fetal bovine serum (FBS), 2 mM L-glutamine, 100 U/mL penicillin and 100 µg/mL streptomycin. A549, PDL fibroblast, HK2, and HEK293 cells were cultured in Dulbecco’s modified Eagle’s medium (DMEM) supplemented with 10% FBS, 2 mM L-glutamine, 100 U/mL penicillin and 100 µg/mL streptomycin. Cells were incubated in a 37 °C humidified incubator with 5% CO2 and were routinely subcultured using a 0.25% trypsin solution with 0.53 mM EDTA. RPMI 1640 medium, DMEM medium, FBS, L-glutamine, penicillin/streptomycin, phosphate-buffered saline (PBS), trypsin and EDTA were obtained from GIBCO (Grand Island, NY). In addition, 3-(4,5-Dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT), DMSO, Hoechst 33342, propidium iodide (PI), and bovine serum albumin (BSA) were purchased from Sigma Chemical, Inc. (St. Louis, MO, USA). Antibodies directed against CD44, ABCG2, ALDH1A1 and β-actin and the respective secondary antibodies were purchased from Cell Signaling (Danvers, MA, USA). Antibodies directed against CD133 were purchased from Cell Applications, Inc. (San Diego, CA, USA). PNA-A15 was synthesized from PANAGENE (PANAGENE, Deajeon, Korea. Cells were treated with either PNA-A15 or scramble (control) for 48 h. The culture medium was replaced by medium containing PNA-A15 or scramble every day. The cells were then collected for microarray, anchorage-independent growth, spheroid formation, flow cytometry, and Western blotting analysis.
Cells were harvested by centrifugation, and resuspended cells were incubated with lysis buffer (20 mM Tris–HCl (pH 7.5), 1% (v/v) Triton X-100, 150 mM sodium chloride, 10% (v/v) glycerol, 1 mM sodium orthovanadate, 50 mM sodium fluoride, protease inhibitor cocktail (Roche Molecular Biochemical) and 100 mM phenylmethylsulfonyl fluoride) for 30 min on ice. The cellular lysates were collected, and their protein contents were determined using a BCA protein assay kit (Pierce Biotechnology, Rockford, IL, USA). Equal amounts of protein from each sample (80 µg for the CSC marker in the total cell population or 30 µg for other proteins and enriched CSC experiments) were separated by SDS-PAGE and transferred to 0.45-μm nitrocellulose membranes (Bio-Rad). The resulting blots were blocked for 1 h with 5% (w/v) nonfat dry milk in TBST (25 mM Tris–HCl (pH 7.5), 125 mM NaCl and 0.1% (v/v) Tween 20) and incubated with the appropriate primary antibodies at 4 °C overnight. After three washes in TBST, the blots were incubated with horseradish peroxidase (HRP)-conjugated secondary antibodies for 2 h at room temperature. Finally, protein bands were detected using an enhanced chemiluminescence substrate (Supersignal West Pico; Pierce, Rockford, IL, USA) and quantified using the analyst densitometry software package (Bio-Rad).
Spheroid formation assay
The formation of spheroids was performed under serum-free conditions in an ultralow attachment plate as previously reported . NSCLC-derived H292 and A549 cells were pretreated with PNA-A15 and scramble (5 μM) for 48 h. Then, cells were detached using 1 mM EDTA and suspended into single cells. These cells were grown in a 24-well ultralow attachment plate at a density of 2.5 × 103 cells/well in stem cell media (SCM) (i.e., 0.8% (w/v) methylcellulose-based serum-free medium (Stem Cell Technologies, Vancouver, BC, Canada) supplemented with 20 ng/ml epidermal growth factor (BD Biosciences, San Jose, CA, USA), 20 ng/ml basic fibroblast growth factor and 4 mg/ml insulin (Sigma) for 7 days to form primary spheroids. These primary spheroids were harvested, resuspended as single cells using 1 mM EDTA and cultured in SCM for 14 days in a 24-well ultralow attachment plate to form secondary spheroids.
RNA extraction and microarray analysis
Total RNA was extracted using TRIzol (Invitrogen Life Technologies, Carlsbad, USA); purity and integrity were evaluated using the ND-1000 Spectrophotometer (NanoDrop, Wilmington, USA) and Agilent 2100 Bioanalyzer (Agilent Technologies, Palo Alto, USA). Total RNA was amplified and purified using the TargetAmp-Nano Labelling Kit for Illumina Expression BeadChip (EPICENTRE, Madison, USA) to yield biotinylated cRNA according to the manufacturer’s instructions. After purification, cRNA was quantified using the ND-1000 Spectrophotometer (NanoDrop, Wilmington, USA). Next, 750 ng of labelled cRNA samples were hybridized to each Human HT-12 v4.0 Expression Beadchip for 18 h at 58 °C according to the manufacturer's instructions (Illumina, Inc., San Diego, USA). Detection of the array signal was performed using Amershamfluorolink streptavidin-Cy3 (GE Healthcare Bio-Sciences, Little Chalfont, UK) following the bead array manual. Arrays were scanned with an Illumina bead array Reader confocal scanner according to the manufacturer's instructions. To investigate whether A and T repeats regulate altered gene expression in spheroid formation, two groups of spheroid H292 and A549 cells were compared with H292 and A549 parental cells. Furthermore, whether PNA-A15 alters the expression of genes containing A and T repeats was also assessed. Two groups of spheroids, H292 and A549 cells, were transfected with PNA-A15 and scramble. Each group consisting of 2 samples was studied. Data were uploaded to the Gene Expression Omnibus (GEO: GSE142616).
The down- and upregulated genes in microarray experiments were identified by CU-DREAM (P-value threshold = 0.05) . Microarray samples are shown in Table 1. We intersected the down- and upregulated genes between the two experiments. Next, we searched for A-T repeats around the transcription start sites (TSSs) of those genes . The human genome (GRCh38) was downloaded via Entrez Direct (EDirect) and its software package ncbi-entrez-direct. The chromosomal locations of genes were downloaded from ftp://ftp.ncbi.nlm.nih.gov/gene/DATA/gene2refseq.gz. The genome sequence around a TSS was divided into 200 bins, and each bin was 100 bp in length. Bins 1 to 100 covered 10,000 bp upstream of a gene, whereas bins 101 to 200 covered a 10,000bp intragenic region. In each bin, intact A repeats (length 13–27 bp) were identified, and their lengths were summed. For instance, a bin consisted of three A repeats that were 12, 13, and 15 in length. The length sum of this bin was 12 + 13 + 15 = 40. A repeat that overlapped between two bins contributed to each bin proportionately.
For the cytotoxicity assay, H460, H23, H292, A549, HK2, HEK293 and PDL fibroblast cells were seeded onto 96-well plates at a density of 1 × 104 cells/well and allowed to adhere by incubation overnight. Cells were then treated with various concentrations of PNA-A15 or scramble plasmid (0–40 μM) for 48 h and then analysed for cell viability using the 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT) assay. Cells were incubated with 500 µg/ml MTT for 4 h at 37 °C, and the intensity of the formazan product after solubilization in 100 µl DMSO was measured at 570 nm using a microplate reader (Anthros, Durham, NC, USA). Relative cell viabilities were calculated by dividing the absorbance of the treated cells by that of the control cells. The half-maximal inhibition concentration (IC50) was determined from four independent experiments using GraphPad Prism 5.0 software (La Jolla, CA).
Anchorage-independent growth assay
After treatment, cells were subjected to an anchorage-independent growth assay based on colony formation in soft agar. The bottom layer was prepared using a 1:1 (v/v) mixture of RPMI-1640 or DMEM medium containing 10% (v/v) FBS and 1% (w/v) agarose. This mixture was allowed to solidify for 30 min, after which the upper cellular layer composed of RPMI-1640 or DMEM medium with 0.3% (w/v) agarose, 10% (v/v) FBS and 1 × 103 cells/ml was prepared and added. Finally, RPMI-1640 or DMEM medium containing 10% (v/v) FBS was added over the upper layer, and the cells were incubated at 37 °C for 14 days. Afterwards, colony formation was observed and imaged using phase-contrast microscopy (Olympus IX51 with DP70). The relative colony numbers and sizes were calculated relative to those for the untreated cells.
H460, H292, H23 and A549 cells were pretreated with PNA-A15 and scramble (5 μM) for 48 h, and then, cells were plated and allowed to form primary and secondary spheroids as detailed above. At day 14 of the secondary spheroid, spheroids were harvested by centrifugation and suspended into single cells using 1 mM EDTA. The resuspended cells were incubated on ice with a rabbit anti-CD133 antibody for 1 h. Next, the primary antibody was removed, and the cells were washed and incubated for 30 min with an Alexa Fluor 488-conjugated goat anti-rabbit IgG (H + L) secondary antibody (Life Technologies, Eugene, OR, USA). After washing, the fluorescence intensity was determined by flow cytometry using a 488-nm excitation beam and a 519-nm bandpass filter (FACSort; Becton Dickinson, Rutherford, NJ, USA). The mean fluorescence intensity was quantified using CellQuest software (Becton Dickinson).
Spheroid viability assay
Spheroid viability was determined using a water-soluble tetrazolium salt (WST) surrogate viability assay according to the manufacturer’s instructions (Roche Diagnostic GmbH, Mannheim, Germany). Briefly, at day 14 of secondary spheroid formation, as detailed in the “Single 3D spheroid formation assay” section, the cells were incubated with 10% (v/v) WST-1 at 37 °C for 30 min. The intensity of the formazan product was determined using a microplate reader (Anthros, Durham, NC, USA), and the relative cell viability (%) was calculated as the absorbance of PNA-A15 or scramble-treated cells relative to untreated cells.
Single 3D spheroid formation assay
For CSC-rich population establishment, the spheroid culture assay used in this study was slightly modified from a previously described method [2, 9, 10]. Cells were seeded onto a 24-well ultralow attachment plate with approximately 2.5 × 103 cells/well in 0.8% methylcellulose-based serum-free medium supplemented with 20 ng/ml epidermal growth factor and 4 mg/ml insulin. The primary spheroids were allowed to form for 7 days. At day 7 of primary spheroid culture, primary spheroids were resuspended into single cells using 1 mM EDTA. Again, 2.5 × 103 cells/well were seeded onto a 24-well ultralow attached plate. Secondary spheroids were allowed to form for 14 days. For the single three-dimensional (3D) spheroid-formation assay, cells were allowed to form primary and secondary spheroids as detailed above. At day 14 of secondary spheroid formation, they were dissociated into a single spheroid of the same size, and each spheroid was then treated with 5 µM PNA-A15 or scramble and allowed to grow in the indicated time. Phase-contrast images of the secondary spheroids were obtained at days 0, 3 and 7 after PNA-A15 or scramble treatment under a phase-contrast microscope (Nikon Eclipse Ts2).
Nuclear staining assay
At day 14 of secondary spheroid formation, these spheroids were dissociated into a single spheroid of the same size, and each spheroid was then treated with 5 µM PNA-A15 or scramble and allowed to grow in the indicated time. Phase-contrast images of the secondary spheroids were obtained at days 0, 3 and 7 after PNA-A15 or scramble treatment under a phase-contrast microscope. On day 7, every individual spheroid was incubated with 10 μg/ml Hoechst 33342 for 30 min followed by 5 μg/ml PI for 5 min. Nuclear condensation and DNA fragmentation of apoptotic cells were visualized using a fluorescence microscope (Olympus IX5; 40 ×) equipped with a DP70 digital camera system (Olympus, Tokyo, Japan).
Data from three or more independent experiments are presented as the mean ± standard deviation (SD). Multiple comparisons for significant differences between multiple groups were performed using analysis of variance (ANOVA) followed by individual comparisons with Scheffe’s post hoc test. Statistical significance was accepted at the P < 0.05 level.
A-T repeats regulate altered gene expression in CSC-rich spheroids
CSCs from H292 and parental cell lines were prepared following previous studies [7, 10] and tested for CD133, CD44, ABCG2, and ALDH1A1, which are well-known CSC markers [10, 11]. We found that the expression levels of CSC markers CD133, CD44, ABCG2 and ALDH1A1 were all dramatically increased in the enriched CSC population compared to their parental cells (Fig. 1a).
A previous study demonstrated that genes containing long mononucleotide A-T repeats are overexpressed in various cancers . To determine whether A-T repeats regulate genes differently between lung cancer and CSCs, we compared A-T repeats around TSS to the spheroid formation population of H292 and A549 lung cancer cells. We divided the expression profiling of genes in the spheroid population compared to the total population of their parental cells into 3 groups, up, down, and neutral, depending on the levels of mRNA. The results showed that the T repeat densities of downregulated genes were significantly increased compared with those of neutral genes in the spheroid population of H292 and A549 cell lines (P-values = 1.18E-04 and 4.65E-10, respectively). These results suggest that genes containing A-T repeats of cancer cells and CSCs are differentially expressed (Fig. 1b, c).
Then, the candidate genes were selected from overlapping T repeat densities of downregulated genes between H292 and A549 cells and were classified by a biological process based on Gene Ontology analysis as shown in Fig. 1d and Additional file 1: Table S1. Figure 1e shows that the 75 overlapping genes were filtered with a P-value ≤ 0.05 in both H292 and A549 cells and were related to genes depending on their functions. We found that dysregulated A and T repeat densities of CSCs in both cell lines were different.
PNA-A15 is not toxic to normal cells but effectively kills cancer cells
Prior to determining the effect of PNA-A15 on CSC properties, the appropriate noncytotoxic concentrations were evaluated. To study the effect of PNA-A15, human NSCLC H460, H292, H23, and A549 cells were treated with various concentrations of compound (0, 2.5, 5, 10, 20, and 40 µM) for 48 h, and then their cell viability was determined by the MTT viability assay. PNA-A15 was found to be nontoxic at concentrations below 5 µM, whereas no cytotoxic effect was found for the scramble control (Fig. 2a–d). Interestingly, PNA-A15 showed no cytotoxic effect on normal cells that contained HK2, HEK293 and PDL fibroblast at all tested concentrations (0–40 µM) (Fig. 2e–g). Thus, nontoxic concentrations of PNA-A15 were further tested for their effects on CSC-like gene expression and phenotypes.
PNA-A15 altered the transcription of CSC genes containing A-T repeats
We determined genome-wide gene expression after spheroid cancer cell transfection with PNA-A15. Cells were pretreated with PNA-A15 or scramble PNA oligo for 2 days and then evaluated for spheroid formation in ultralow attachment plates. After that, pretreated cells were seeded at a low density and allowed to form primary spheroids for 7 days. The primary spheroids were then resuspended into single cells, and secondary spheroids were allowed to grow for 14 days. Then, we analysed the expression profiling of A and T repeat genes in the CSC-rich population treated with PNA-A15 compared to their scramble cells. We found that the A repeat densities of upregulated genes were significantly reduced compared with those of neutral genes (P values = 1.30E-08), and the T repeat densities of downregulated and upregulated genes were also significantly altered compared with neutral (P values = 3.69E-04 and 1.65E-19, respectively) in H292 spheroids. Furthermore, T repeat densities of down- and upregulated genes were also significantly altered compared with neutral genes (P-values = 6.53E-07 and 2.39E-02, respectively) in spheroid A549 cells. These results suggest that PNA-A15 significantly alters the expression of genes containing A- and T- repeats in spheroid H292 and A549 cancer cell lines (Fig. 3a, b).
To determine PNA-A15 dysregulated genes in CSCs, we selected candidate genes from overlapping A and T repeat densities of both down- and upregulated genes between the CSC-rich population treated with PNA-A15 and their scramble cells. All genes were classified by a biological process based on Gene Ontology analysis as shown in Fig. 3c and Additional file 1: Table S2. The validation genes were selected based on function and their biological processes. The data represent 69 overlapping genes filtered with a P value ≤ 0.01, and we found that PNA-A15 dysregulated the A and T repeat densities of downregulated genes that are related to genes acting as proto-oncogenes and inducing tumorigenesis in both cell lines. PNA-A15 dysregulated the A and T repeat densities of upregulated genes, which act as a tumour suppressor and induce apoptotic processes (Fig. 3e, f and Additional file 1: Table S3).
PNA-A15 suppresses CSC phenotypes
The ability of cancer cells to form 3D spheroids and to grow in anchorage-independent conditions is widely used to assess CSC phenotypes [2, 10]. Therefore, we investigated the role of PNA-A15 in CSC behaviours. Accordingly, the effect of PNA-A15 on the growth and survival of NSCLC-derived cells was evaluated under these conditions. First, H460, H292, H23 and A549 cells were exposed to noncytotoxic concentrations of PNA-A15 for 48 h and subsequently analysed for colony formation in soft agar by recording the relative number and area of colonies after 14 days compared to untreated control or scramble-treated cells (Fig. 4). We found that the PNA-A15-treated cells exhibited a significant decrease in the number and size of colonies compared to the untreated or scramble-treated control cells.
Then, we evaluated spheroid formation in ultralow attachment plates. We found that treatment of the cells with PNA-A15 significantly inhibited spheroid formation (Fig. 5). These results suggest that PNA-A15 inhibits the growth and survival characteristics of CSCs in NSCLC-derived cells.
To determine lung CSC markers, immunocytochemistry experiments assessing the expression of CD133 were performed and analysed by flow cytometry. The results showed that PNA-A15 treatment decreased the mean fluorescence intensity of CD133 compared with the untreated or scramble-treated control (Fig. 6).
PNA-A15 suppresses CSC maintenance in CSC-rich spheroids
We further monitored the effect of PNA-A15 on the maintenance of CSCs in the CSC-rich population. The spheroids were either treated with PNA-A15 or scramble plasmid or untreated for 48 h. Cell viability was then evaluated by WST assay. PNA-A15 significantly decreased the viability of CSCs in H292 cells compared to parental and scrambled cells (Fig. 7a). Every individual spheroid was subsequently treated with PNA-A15 at 5 μM for 0–7 days. We found that PNA-A15 treatment significantly decreased the size of the spheroids and induced apoptosis at days 3 and 7 after treatment as determined by increased chromatin condensation and nuclear fluorescence following Hoechst 33342 and PI staining (Fig. 7b). These results support the above observed inhibitory effect of PNA-A15 on CSCs.
Furthermore, CSC markers, including CD133, CD44, ABCG2 and ALDH1A1, were analysed by Western blot analysis. Figure 7c shows that PNA-A15 caused a significant decrease in CD133, CD44, ABCG2 and ALDH1A1 expression compared with the untreated or scramble-treated control level in parental cells. The results were confirmed in the 3D-CSC-enriched population. The enriched CSC population was treated with PNA-A15 at 5 μM for 48 h, and the expression of CD133, CD44, ABCG2 and ALDH1A1 was determined by Western blot analysis. Figure 7d shows that PNA-A15 treatment reduced CD133, CD44, ABCG2 and ALDH1A1 expression compared with the untreated or scramble-treated control level. These results strongly support the role of PNA-A15 as a multiple gene targeting agent for lung CSCs.
Here, we found that A and T repeats regulate gene expression in CSCs; therefore, PNA-A15 represents a promising targeted therapy agent in lung cancer by downregulating multiple oncogenes and upregulating multiple tumour-suppressor genes. We also revealed that A-T repeat-containing genes that are downregulated in CSCs are involved in CSC morphology, such as BRWD1, FLOT2, SKA3, and TGF-β1. BRWD1 is a transcriptional activator that inhibits cell proliferation by coordinately regulating MYC, and this mechanism plays a role in the self-renewal of CSCs . TGF-β1 is a multifunctional cytokine that plays a critical role in the regulation of CD133 expression. Previous work revealed that TGF-β1 upregulates CD133 expression in hepatocellular carcinoma . FLOT2 is a proto-oncogene that regulates cell migration and invasion in cancer. FLOT2 is involved in cell cycle progression and regulates signalling pathways that are responsible for stemness and induce cadherin binding . FLOT2 is necessary for TGF-β1 and induced epithelial-mesenchymal transition (EMT) in gastric cancer. Suppression of FLOT2 results in decreased cell invasion through repressing TGF-β1-mediated EMT and decreased CSCs [14, 15]. Moreover, SKA3 is a proto-oncogene that induces cell cycle progression. Previous work revealed that SKA3 mRNA expression was elevated and correlated with poor survival outcomes in lung cancer patients . SKA3 acts as an oncogene by directly binding to EGFR and promoting cancer metastasis [16, 17]. Moreover, the SKA3 gene acts as a proto-oncogene and induces tumorigenesis by inhibiting p53 in the apoptotic process. Additionally, p53 regulates the expression of CSC markers. Thus, activating p53 not only increases apoptotic induction in tumour cells but also suppresses CSC self-renewal.
After treatment of PNA-A15 into spheroids of H292 and A549 cells, we found that the expression of many A and T repeat-containing genes changed in both the up- and downregulation directions. CSCs are difficult to therapeutically target because they are inherently resistant to cytotoxic and targeted drugs, and they also evade radiotherapy and immune surveillance [1, 18]. Recent data show that CSC populations underlie significant diversification and plasticity and disguise heterogeneity . Our experiments also demonstrated heterogeneity in the gene expression of CSCs. This heterogeneity prevents the effectiveness of therapy targeting a single target. Therefore, PNA-A15 is a promising targeted therapy in lung cancer by downregulating multiple oncogenes and upregulating multiple tumour suppressor genes. PNA-A15 treatment downregulated A-T repeat-containing genes. Interestingly, many of the genes possess stem cell maintenance and functions, such as the ARID1A, MAPK12, PLOD3, and RAC1 genes [19,20,21,22]. Previous studies revealed that ARID1A regulates SOX9 expression, which is a stem cell transcription factor that maintains the self-renewal property of intestinal stem cells . In contrast, the deletion of the ARID1A gene results in the loss of stem cells and increases apoptosis in adult mice [19, 23]. MAPK12 also modulates the stemness of various CSC types by inducing tumorigenesis and cancer aggressiveness . A previous study revealed that MAPK12 significantly increased EMT and promoted CSC regulation in breast cancer . In addition, PNA-A15 interferes with both cancer and CSC phenotypes through the regulation of PLOD3, which is a potent inducer of lung cancer metastasis via the RAS-MAPK pathway in vivo . Moreover, PLOD3 also interacts with STAT3, which is a stem cell transcription factor that maintains the self-renewal property of lung stem cells . PNA-A15 also downregulated the RAC1 gene, which mediated lung tumour growth and increased cell proliferation in vivo. Previous work demonstrated that RAC1 is critically involved in NSCLC migration and lung CSC formation and that RAC1 served as a useful therapeutic target by inhibiting tumour initiation and metastasis of CSCs in lung cancer . Thus, PNA-A15 not only interferes with tumour cells but also suppresses several genes possessing the self-renewal of CSC functions.
PNA-A15 treatment upregulated A-T repeat-containing genes that act as tumour suppressor genes, such as BAHD1, CCAR1, EGR1, and SLC5A11. These genes prevent tumour formation and increase apoptotic processes [24,25,26]. CCAR1 is a tumour suppressor gene related to p53-induced apoptosis . BAHD1 is a tumour suppressor gene, and this anti-proliferative protein inhibits cell cycle progression from G0/G1 to s phase . EGR1 is a tumour suppressor gene that regulates cell survival by activating the expression of p53 to prevent tumour formation . This role is important because p53 regulates the expression of CSC genes. Thus, activating p53 not only increases apoptotic induction in tumour cells but also suppresses self-renewal of CSCs. CCAR1 mediates the expression of multiple cell cycle regulatory genes and plays a role in cell cycle progression as a transducer of Notch signalling . The Notch signalling pathway is involved in the process of stem cell initiation [28, 29]. In some solid tumours, dysregulation of the Notch signalling pathway is correlated with tumour initiation and increased tumour sphere formation . These findings suggest that PNA-A15 may target the Notch signalling pathway in lung CSCs.
PNA-A15 is not toxic to normal cells but effectively kills cancer cells. Previously, we reported that short A and T repeats (2–9 bp) are more abundant in tissue-specific genes, whereas long A and T repeats (10–30 bp) are more abundant in housekeeping genes . For normal cells, genes controlling differentiation are essential [5, 6]. On the other hand, to maintain the cancer and stem cell phenotypes, cancer cells must prevent full differentiation and use housekeeping genes to play a major role in cell survival. Therefore, long A and T repeats may influence cancer and stem cell survival more than completely differentiated cells.
The stemness of a cancer cell can drive drug resistance through certain mechanisms. For instance, ABCG2, an indicator of cancer stemness, is highly involved in acquiring multidrug resistance to cancer chemotherapeutics such as Cisplatin and Etoposide [31, 32]. In accompany with reduction of CSC, we found that PNA-A15 treatment could reduce ABCG2 expression. Future studies to determine the efficacy of PNA-A15 over drug-resistant phenotype are interesting.
In conclusion, A-T repeat transcriptional control is essential for CSC biology. Furthermore, interfering with PNA-A15 activity can disrupt CSC growth by dysregulating many genes. This property of PNA-A15 suggests the potential of PNA-A15 in lung CSC treatment to overcome genetic heterogeneity. Moreover, PNA-A15 might be a potential therapeutic drug for many other CSCs in addition to lung cancer.
While PNA-A15 targets CSC, standard chemotherapeutic agents are effective in cancer cells, future study a combinatorial drug between PNA-A15 and standard chemotherapeutic agents may improve treatment efficacy.
Availability of data and materials
The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.
Cancer stem cells
Peptide nucleic acid
Epidermal growth factor receptor
Non-small cell lung cancer
Fetal bovine serum
Cluster of differentiation 133
Cluster of differentiation 44
ATP-binding cassette super-family G member 2
Aldehyde dehydrogenalse 1
Stem cell media
Gene Expression Omnibus
Connecting up- and down-regulation expression analysis of microarrays program
- IC50 :
The half-maximal inhibition concentration
Water-soluble tetrazolium salt
Bromodomain and WD repeat domain containing 1
Spindle and kinetochore associated complex subunit 3
Transforming growth beta 1
AT-rich interaction domain 1A
Mitogen-activated protein kinase 12
Procollagen-lysine,2-oxoglutarate 5-dioxygenase 3
Ras-related C3 botulinum toxin substrate 1
SRY-box transcription factor 9
Signal transducer and activator of transcription 3
Bromo adjacent homology domain containing 1
Cell division cycle and apoptosis regulator 1
Early growth response protein 1
Solute carrier family 5 member 11
Ayob AZ, Ramasamy TS. Cancer stem cells as key drivers of tumour progression. J Biomed Sci. 2018;25(1):20.
Cao L, et al. Sphere-forming cell subpopulations with cancer stem cell properties in human hepatoma cell lines. BMC Gastroenterol. 2011;11:71.
Vinogradov S, Wei X. Cancer stem cells and drug resistance: the potential of nanomedicine. Nanomedicine. 2012;7(4):597–615.
Nan X, et al. EGFR TKI as first-line treatment for patients with advanced EGFR mutation-positive non-small-cell lung cancer. Oncotarget. 2017;8(43):75712–26.
Aporntewan C, et al. Upstream mononucleotide A-repeats play a cis-regulatory role in mammals through the DICER1 and Ago proteins. Nucleic Acids Res. 2013;41(19):8872–85.
Pin-On P, et al. Targeting high transcriptional control activity of long mononucleotide A-T repeats in cancer by Argonaute 1. Gene. 2019;699:54–61.
Bhummaphan N, Chanvorachote P. Gigantol suppresses cancer stem cell-like phenotypes in lung cancer cells. Evid Based Compl Alt Med. 2015;2015:836564.
Aporntewan C. Connection up-and down-regulation expression analysis of microarrays. Asian Biomed. 2011. https://doi.org/10.5372/1905-7415.0502.034.
Bhummaphan N, et al. Lusianthridin targeting of lung cancer stem cells via Src-STAT3 suppression. Phytomedicine. 2019;62:152932.
Bhummaphan N, et al. cancer stem cell-suppressing activity of chrysotoxine, a Bibenzyl from Dendrobium pulchellum. J Pharmacol Exp Ther. 2018;364(2):332–46.
Chanvorachote P, Luanpitpong S. Iron induces cancer stem cells and aggressive phenotypes in human lung cancer cells. Am J Physiol Cell Physiol. 2016;310(9):C728-39.
Mandal M, et al. BRWD1 orchestrates epigenetic landscape of late B lymphopoiesis. Nat Commun. 2018;9(1):3888.
You H, Ding W, Rountree CB. Epigenetic regulation of cancer stem cell marker CD133 by transforming growth factor-beta. Hepatology. 2010;51(5):1635–44.
Li Q, et al. miR-449a targets Flot2 and inhibits gastric cancer invasion by inhibiting TGF-β-mediated EMT. Diagn Pathol. 2015;10:202.
Papageorgis P. TGF β signaling in tumor initiation, epithelial-to-mesenchymal transition, and metastasis. J Oncol. 2015;2015:587193.
Hu D-D, et al. SKA3 promotes lung adenocarcinoma metastasis through the EGFR-PI3K-Akt axis. Biosci Reports. 2020. https://doi.org/10.1042/BSR20194335.
Hu R, et al. SKA3 promotes cell proliferation and migration in cervical cancer by activating the PI3K/Akt signaling pathway. Cancer Cell Int. 2018;18(1):183.
Boesch M, et al. Heterogeneity of cancer stem cells: rationale for targeting the stem cell niche. Biochem Biophys Acta. 2016;1866(2):276–89.
Hiramatsu Y, et al. Arid1a is essential for intestinal stem cells through Sox9 regulation. Proc Natl Acad Sci USA. 2019;116(5):1704–13.
Myant KB, et al. Rac1 drives intestinal stem cell proliferation and regeneration. Cell Cycle. 2013;12(18):2973–7.
Xu M, et al. (2018) Role of p38γ MAPK in regulation of EMT and cancer stem cells. Biochimica et Biophysica Acta BBA Mol Basis Dis. 1864;11:3605–17.
Baek J-H, et al. PLOD3 promotes lung metastasis via regulation of STAT3. Cell Death Dis. 2018;9(12):1138.
Cui N, et al. Epigenetic regulator ARID1A and stem cell transcription factor SOX9 in the maintenance of pancreatic ductal cell differentiation state and development of intraductal papillary mucinous neoplasia (IPMN) and pancreatic ductal adenocarcinoma (PDAC). Transl Cancer Res. 2018. https://doi.org/10.21037/tcr.2018.07.25.
Lu P, et al. p53-Bad: a novel tumor suppressor/proapoptotic factor hybrid directed to the mitochondria for ovarian cancer gene therapy. Mol Pharm. 2019;16(8):3386–98.
Chang TS, et al. Inhibition of CCAR1, a coactivator of β-Catenin suppresses the proliferation and migration of gastric cancer cells. Int J Mol Sci. 2017. https://doi.org/10.3390/ijms18020460.
Inoue K, Fry EA. Tumor suppression by the EGR1, DMP1, ARF, p53, and PTEN network. Cancer Invest. 2018;36(9–10):520–36.
Muthu M, Cheriyan VT, Rishi AK. CARP-1/CCAR1: a biphasic regulator of cancer cell growth and apoptosis. Oncotarget. 2015;6(9):6499–510.
Pelullo M, et al. Wnt, Notch, and TGF-β pathways impinge on hedgehog signaling complexity: an open window on cancer. Front Genet. 2019;10:711.
Yang L, et al. Targeting cancer stem cell pathways for cancer therapy. Signal Transduct Target Ther. 2020;5(1):8.
Abel EV, et al. The Notch pathway is important in maintaining the cancer stem cell population in pancreatic cancer. PLoS ONE. 2014;9(3):e91983–e91983.
Wang X, Zhang H, Chen X. Drug resistance and combating drug resistance in cancer. Cancer Drug Resist. 2019;2(2):141–60.
Tang Y, et al. ABCG2 regulates the pattern of self-renewing divisions in cisplatin-resistant non-small cell lung cancer cell lines. Oncol Rep. 2014;32(5):2168–74.
This research project was supported by a 2019 Research Chair Grant from the National Science and Technology Development Agency, Thailand; the Second Century Fund (C2F), Chulalongkorn University; and The Chulalongkorn Academic Advancement into Its 2nd Century Project. Piyapat Pin-on is supported by the research fund, King Mongkut's Institute of Technology Ladkrabang (2563–0216002).
Ethics approval and consent to participate
Not applicable. All experiments referred to in this report were conducted with commercially available cell lines and/or primary human cells with all appropriate consents available from the suppliers.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The lists of down-regulated genes between CSCs and parental cells in H292 and A549 cell lines. The data were filtered with P ≤ 0.05. Table S2. The lists of down-regulated genes between PNA-A15 and their scramble in H292 and A549 cell lines. The data were filtered with P ≤ 0.05. Table S3. The lists of up-regulated genes between PNA-A15 and their scramble in H292 and A549 cell lines. The data were filtered with P ≤ 0.05.
About this article
Cite this article
Bhummaphan, N., Pin-on, P., Phiboonchaiyanan, P.P. et al. Targeting multiple genes containing long mononucleotide A-T repeats in lung cancer stem cells. J Transl Med 19, 231 (2021). https://doi.org/10.1186/s12967-021-02902-6
- Cancer stem cells
- Lung cancer
- Mononucleotide A-T repeats
- Hallmark of cancer
- Universal target