Genome-wide profiling of long non-coding RNA expression patterns in the EGFR-TKI resistance of lung adenocarcinoma by microarray

  • Authors:
    • Ying Wu
    • Dan-Dan Yu
    • Yong Hu
    • Dali Yan
    • Xiu Chen
    • Hai‑Xia Cao
    • Shao-Rong Yu
    • Zhuo Wang
    • Ji-Feng Feng
  • View Affiliations

  • Published online on: April 20, 2016     https://doi.org/10.3892/or.2016.4758
  • Pages: 3371-3386
Metrics: Total Views: 0 (Spandidos Publications: | PMC Statistics: )
Total PDF Downloads: 0 (Spandidos Publications: | PMC Statistics: )


Abstract

Mutations in the epidermal growth factor receptor (EGFR) make lung adenocarcinoma cells sensitive to EGFR tyrosine kinase inhibitors (TKIs). Long-term cancer therapy may cause the occurrence of acquired resistance to EGFR TKIs. Long non-coding RNAs (lncRNAs) play important roles in tumor formation, tumor metastasis and the development of EGFR-TKI resistance in lung cancer. To gain insight into the molecular mechanisms of EGFR-TKI resistance, we generated an EGFR-TKI-resistant HCC827-8-1 cell line and analyzed expression patterns by lncRNA microarray and compared it with its parental HCC827 cell line. A total of 1,476 lncRNA transcripts and 1,026 mRNA transcripts were dysregulated in the HCC827‑8-1 cells. The expression levels of 7 chosen lncRNAs were validated by real-time quantitative PCR. As indicated by functional analysis, several groups of lncRNAs may be involved in the bio-pathways associated with EGFR-TKI resistance through their cis- and/or trans‑regulation of protein-coding genes. Thus, lncRNAs may be used as novel candidate biomarkers and potential targets in EGFR-TKI therapy in the future.

Introduction

Lung cancers with epidermal growth factor receptor (EGFR) gene mutations are a well-understood subgroup of lung adenocarcinomas (LACs). This subgroup is characterized by sensitivity to EGFR tyrosine kinase inhibitors (TKIs; e.g. gefitinib or erlotinib), and has high occurrence rates in females, non-smokers and Asians (1). However, a median therapy (10–14 months) with EGFR TKI may aggravate the disease progression in most patients in this subgroup. The mechanisms of EGFR-TKI resistance are multi-factorial and several of them have been reported, such as T790M mutation in exon 20 of EGFR, MET amplification and depletion of phosphatase and tensin homolog (PTEN) (24). However, nearly 30% of resistance mechanisms are unknown. Thus, elucidating the molecular mechanisms of EGFR-TKI resistance is essential for the identification of key biomarkers.

In the last two decades, non-coding RNAs (ncRNAs) have been gradually accepted as key factors in the process of epigenetic regulation rather than 'transcription noise' (5). As recently reported, ncRNAs take part in the pathogenesis of non-small cell lung cancer (NSCLC), providing new biological insights into this disease (6,7). Long non-coding RNAs (lncRNAs), which are non-protein coding transcripts with length >200 nucleotides (nt), are involved in regulating the occurrence, invasion, metastasis and chemotherapeutic resistance of lung cancer (6,8,9). lncRNAs also take part in the regulation of EGFR-TKI resistance. GAS5, one lncRNA, was found to be overexpressed in EGFR-TKI-sensitive cells compared with its expression in resistant cells, and was found to enhance gefitinib-induced cell death in innate EGFR-TKI-resistant LAC cells with wild-type EGFR via downregulation of IGF-1R expression (10). These findings indicate that lncRNAs may be promising biomarkers as diagnostic and therapeutic targets in the resistance of EGFR-TKIs.

The present study presents the lncRNA expression profiles in 3 replicate gefitinib-sensitive HCC827 and gefitinib-resistant HCC827-8-1 cells in pairs by microarray. Then, 7 of the differentially expressed lncRNAs were validated by real-time quantitative PCR (RT-qPCR) in HCC827 and HCC827-8-1 cells. We also predicted the functions of differentially expressed lncRNAs through their co-expressed protein-coding genes.

Materials and methods

Cell culture

The human NSCLC H827 cell line harboring the EGFR exon 19 deletion (Del E746-A750) was obtained from the Shanghai Institutes for Biological Sciences, Chinese Academy of Cell Resource Center. We generated a gefitinib-resistant cell line by exposing HCC827 cells to increasing concentrations of gefitinib as described in our previous study (11). One individual clone HCC827-8-1 was isolated and independently confirmed to be resistant to gefitinib. HCC827-8-1 cells were 348-fold more resistant to gefitinib than the parental HCC827 cells. The cells were cultured in RPMI-1640 medium supplemented with 10% fetal bovine serum (FBS) (Gibco, Carlsbad, CA, USA) at 37°C, in a 5% CO2 and humidified atmosphere.

RNA extraction

HCC827 and HCC827-8-1 cells were seeded in 6-well plates (1×105 cells/well) for 72 h in 3 replicate wells, and then resuspended in 500 μl lysis buffer. Total RNA was extracted from lysis buffer using the mirVana miRNA Isolation kit procedure (Applied Biosystem, Foster City, CA, USA), according to the manufacturer's specifications, and eluted with 100 ml of nuclease-free water. The yield of RNA was quantified by the NanoDrop ND 2000 (Thermo Scientific, Waltham, MA, USA) and the RNA integrity was assessed using Agilent Bioanalyzer 2100 (Agilent Technologies, Santa Clara, CA, USA).

lncRNA and mRNA microarray expression profiling

Total RNA was transcribed to double-stranded cDNA, and then synthesized into cRNA and labeled with cyanine 3-CTP. The labeled cRNAs were hybridized onto the microarray. The Agilent human lncRNA array (4×180K) contains 32,776 human mRNAs and 78,243 human lncRNAs, which are derived from authoritative databases, including RefSeq, Ensemble, GenBank and the Broad Institute. After washing, the arrays were scanned by the Agilent Scanner G2505C (Agilent Technologies). Raw datum was extracted using Feature Extraction (version 10.7.1.1; Agilent Technologies). The microarray profiling was conducted in the laboratory of the OE Biotechnology Co. (Shanghai, China).

RT-qPCR validation of 7 differentially expressed lncRNAs

Total RNA was extracted from HCC827 and HCC827-8-1 cells using TRIzol reagent (Invitrogen, Carlsbad, CA, USA) according to the manufacturer's instructions. One microgram of total RNA was reverse transcribed in a final volume of 20 μl using PrimerScript RT Master Mix (Takara, Dalian, China). The reverse transcription reaction was carried out under the following conditions: 37°C for 15 min; 85°C for 5 sec, and then hold on 4°C. The RT-qPCR was performed using SYBR-Green PCR Mix (Roche, Mannheim, Germany) on the ABI 7900 system (Applied Biosystems) according to the manufacturer's instructions. The primer sequences are listed in Table I. β-actin was used as an internal control to normalize the amount of total RNA in each sample. Each sample was run in triplicate for analysis. At the end of the PCR cycles, melting curve analysis was performed to validate the specific generation of the expected PCR product. The expression levels of lncRNAs were normalized to internal control gene β-actin and were calculated using the 2−ΔΔCt method (12).

Table I

The primer sequences used in the present study.

Table I

The primer sequences used in the present study.

Target IDForward primerReverse primer
FR165245 GAGGGTTTGGCTGTTTGCTG ACCCCCACTTAGAGACCAGAA
ENST00000464359 GCAACAACCACTTGGCTCAG GCAGAGGACACGAACTCACA
ENST00000602301 GTTACCTCCTCATGCCGGAC AAAAGGGTCAGTAAGCACCCG
NONHSAT107900 GGCTGCATTTGTTTCTCGCA CCCGCCCAGCTATAGTCAAG
NONHSAT082241 TGCCAAAACTCACCAGCTACA GGAGCGGTATGTGCTAGACC
ENST00000434951 TGGGAGTGAATGTTCCGGTG CAAGAGGAGCTGTTGTTTGTCC
NONHSAG031748 GGATGTGCACGCATGAACTG ACTCCAGCCAAGGTGGTTTT
β-actin GATGAGATTGGCATGGCTTT CACCTTCACCGTTCCAGTTT
Differential expression level of mRNAs and lncRNAs from the microarray

After quantile normalization, raw signals from the microarray were log2 transformed. Differential expression of an mRNA or lncRNA was defined by the absolute value of fold-change (FC) >2 (gefitinib-sensitive HCC827 cells=1) and P-value <0.05 (Student's t-test), with the 3 parallel samples in the HCC827 or HCC827-8-1 group having detectable signals compared with the background. The differentially expressed mRNAs were submitted to the NCBI (gene_go_information) and KEGG database analyzed by Python program to be classified into different Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) annotation groups.

Co-expression of lncRNAs with mRNAs and functional prediction

Most of the lncRNAs in the current databases have not yet been functionally annotated. Thus, the prediction of their functions is based on the functional annotations of their co-expressed mRNAs. This method was originally described by Guttman et al (13). In brief, first, for every dysregulated lncRNA, Pearson correlation coefficient (PCC) of its expression with that of each dysregulated mRNA was calculated to find its co-expressed mRNAs, with PCC >0.7 or ≤0.7 and P-value of PCC <0.05 being statistically significant. Then, a functional enrichment analysis of the co-expressed mRNAs was conducted using the hypergeometric cumulative distribution function, and the enriched GO/KEGG pathway annotations were assigned to the lncRNA as its predicted functions. The threshold of statistical significance is set as a P-value <0.05 and false discovery rate (FDR) <0.05 [under the control of the Benjamini and Hochberg procedure (14)].

In addition, on the basis of co-expression, we further explored how these dysregulated lncRNAs may exert their functions through cis- and/or trans-regulating protein-coding genes. We defined cis-regulated genes as protein-coding genes co-expressed with one dysregulated lncRNA and within 300 kb in genomic distance in the same allele.

We potentially defined trans-regulated protein-coding genes as co-expressed and beyond 300 kb in genomic distance from, or, on the other allele of, differentially expressed lncRNAs. In other words, for any protein-coding gene co-expressed with an lncRNA, if it did not fit the criteria of cis-regulated, it was categorized as potentially trans-regulated.

According to Guttman et al (13,15), specific lncRNAs participating in certain biological pathways are transcriptionally regulated by key transcription factors (TFs) that regulate these pathways. Thus, to categorize lncRNAs that possibly have trans-regulating functions, we compared the mRNAs that co-expressed with these lncRNAs with the mRNAs that are regulatory targets of certain TFs. If the intersection of these 2 groups is large enough (P<0.05; calculated by hypergeometric cumulative distribution function and FDR <0.05, under the control of the Benjamini and Hochberg procedure), then, we predict that these lncRNAs possibly participate in pathways regulated by these TFs. The lncRNA-TF network was constructed using hypergeometric cumulative distribution function with the help of Perl. The graph of the lncRNA-TF network was drawn with the help of Cytoscape 3.01.

Results

General expression profiles of the differentially expressed lncRNAs and mRNAs

We found that 1,476 lncRNA transcripts were differentially expressed between the HCC827-8-1 and paired HCC827 cells, with 703 being upregulated and 773 being downregulated. Among the dysregulated lncRNA transcripts, ENST00000464359 (probe CUST_89973_PI429545380) was the most upregulated, with an FC of 136.11, whereas ENST00000602301 (probe CUST_1064_PI429545384) was the most downregulated, FC being 81.46. According to the absolute value of FC, the dysregulated lncRNA transcripts were stratified into 4 groups: 2 transcripts with FCs >100, 38 transcripts between 10 and 100, 117 transcripts between 5 and 10 and 1,319 transcripts between 2 and 5.

Using the same criteria of lncRNAs, we found that 1,026 mRNA transcripts were dysregulated, with 516 being upregulated and 510 being downregulated. The most upregulated and downregulated mRNA transcripts were FOXP2 (A_24_P16559) and MCOLN2 (A_23_P23639), with FCs of 74.25 and 31.03, respectively. Table II lists the top 20 upregulated and downregulated lncRNAs and mRNAs from our microarray.

Table II

Top 20 upregulated and downregulated lncRNAs and mRNAs in the HCC827 cells compared with the gefitinib-resistant HCC827-8-1 cells.

Table II

Top 20 upregulated and downregulated lncRNAs and mRNAs in the HCC827 cells compared with the gefitinib-resistant HCC827-8-1 cells.

Upregulated lncRNAs
Downregulated lncRNAs
Upregulated mRNAs
Downregulated mRNAs
lncRNAsFClncRNAsFCmRNAsFCmRNAsFC
ENST00000464359136.11 ENST0000060230181.46FOXP274.25MCOLN231.03
NONHSAT130274107.33NONHSAG00229468.92PYDC229.44RNF18327.43
NONHSAT12283347.29 ENST0000054298033.87CHI3L126.90THNSL227.18
NONHSAT12282838.98NONHSAT12225723.86SFTPB19.99PGBD517.46
NONHSAG04858730.11NONHSAG03289622.43GJA116.79FAM26F15.87
uc.225+27.11NONHSAT08224122.04CTGF14.14ZNF65515.23
NONHSAT12282922.16NONHSAT06096520.40CALB113.04IGFBP315.12
NONHSAT12282519.67NONHSAT10117619.36ATP8A112.95ZNF56013.81
NONHSAT12764214.70NONHSAT11861915.69WNT9A10.87RPS6KA613.41
NONHSAT09396814.13 ENST0000042373715.53SPTLC310.79FRMD313.17
NONHSAT06092713.25 ENST0000042469015.04SPOCK19.81BACE213.16
NONHSAT12763812.60NONHSAT12047614.96SPTLC39.48GPR13312.70
FR17395512.50NONHSAG04689914.92SFTPB9.16CYP4F1212.26
NONHSAG05071112.28NONHSAT11862114.67BMP48.77PTGER212.23
XR_254376.111.98NONHSAT10117713.85AQP38.62CKMT1A11.58
NONHSAT10790011.85NONHSAT06078512.80EGR18.27FAM129A11.11
NONHSAT10397211.56NONHSAT08224412.59TOX28.20C17orf10410.12
uc.222+10.73NONHSAT05308612.33GMPR7.92ICAM29.63
TCONS_0001226010.37NONHSAT12225611.48ADORA17.90MUC229.61
FR1652459.81NONHSAT00423211.24CEACAM67.81LOC6439889.53

[i] FC, absolute fold-change (HCC827 cells set as 1).

In the GO pathway analysis, the most common pathways that the dysregulated mRNAs were involved in included chondroitin sulfate metabolic process (GO, biological processes), extracellular space (GO, cellular components) and receptor binding (GO, molecular functions). The most common KEGG pathways involved were graft-versus-host disease, type I diabetes mellitus and viral myocarditis.

In unsupervised hierarchical clustering analysis, the differentially expressed lncRNAs were used to generate a heat map; and they clearly self-segregated into HCC827 and HCC827-8-1 clusters, as shown in Fig. 1.

RT-qPCR validation

To validate the results of the microarray, we chose a total of 7 differentially expressed lncRNA transcripts for RT-qPCR. They could be divided into 2 groups: the first group, including FR165245, which was randomly chosen. In the other group, NONHSAT107900 and NONHSAT082241 were chosen since they were predicted to have cis-regulating potential. ENST00000434951 and ONHSAG031748 were chosen since they were predicted to have trans-regulating potential and involved in constituting the top 100 lncRNA-TF pairs with the most credentiality. ENST00000464359 and ENST00000602301 were chosen for being the most upregulated and downregulated lncRNAs (Table II).

The RT-qPCR results were consistent with that of the microarray, in that all 7 lncRNA transcripts were differentially expressed with the same trend (upregulated or downregulated) and reached statistical significance (P<0.05 for each lncRNA; Student's t-test) as shown in Fig. 2.

lncRNA and mRNA co-expression profiles and lncRNA function prediction

Hundreds of lncRNAs were co-expressed with thousands of mRNAs. For example, ENST00000412387 (CUST_86592_PI429545380) was co-expressed with 3,487 mRNA transcripts and NONHSAT060786 (CUST_61584_PI429545406) with 4,642 mRNA transcripts.

The functions of the differentially expressed lncRNAs were predicted by the GO and KEGG pathway annotations of their co-expressed mRNAs. The lncRNAs were clustered into hundreds of GO and KEGG pathway annotations. Certain pathways are known to be involved in the mechanism of gefitinib resistance, including focal adhesion, cell cycle, cell proliferation and apoptosis (1618). For example, 186 lncRNAs were clustered into apoptosis and 79 in the focal adhesion. As expected, one lncRNA can participate in more than one GO/KEGG pathway involved in the mechanism of gefitinib resistance. For example, ENST00000412387 (CUST_86592_PI429545380) was predicted to have functions in apoptosis, cell cycle, focal adhesion and pathways in cancer.

From the matrix of the lncRNAs and their corresponding KEGG pathway annotations of co-expressed protein-coding genes, we counted and summarized the top 200 annotations with the most credentiality (the lowest P-values). The most frequently predicted functions of the differentially expressed lncRNAs were metabolic pathways, glyoxylate and dicarboxylate metabolism and N-glycan biosynthesis as shown in Fig. 3.

cis-regulation of lncRNAs

A total of 149 lncRNA transcripts with their predicted cis-regulated protein-coding genes were found through accurate genomic mapping, using the above mentioned criteria. Table III lists all the lncRNA transcripts and their potentially cis-regulated mRNA transcripts. For instance, lncRNA NONHSAG010348 was predicted to cis-regulate five mRNAs (corresponding to five mRNA transcripts): COPS7A, TPI1, EMG1, USP5 and LRRC23.

Table III

Dysregulated lncRNA transcripts and their potentially cis-regulated mRNA transcripts.

Table III

Dysregulated lncRNA transcripts and their potentially cis-regulated mRNA transcripts.

lncRNAmRNAPCCP-valueaChrom
ENST00000423737LOC2849300.9907639620.00012822
ENST00000425104SLC35F3−0.8718757290.0235721
ENST00000431813DAPK10.9442543010.0045759
ENST00000438810LOC2849300.9886897390.00019122
ENST00000499202CD27−0.9278158630.00762812
ENST00000499202COPS7A−0.9157417520.01035012
ENST00000499202GAPDH−0.8204314890.04547212
ENST00000499202LTBR−0.9587511210.00251712
ENST00000499202NCAPD2−0.9263499550.00793712
ENST00000499202PLEKHG60.8864770260.01860012
ENST00000508616OR5H15−0.9148574730.0105653
ENST00000510682COMMD100.9311760550.0069425
ENST00000558536C15orf480.9485813590.00389815
ENST00000563635LMO70.9621129370.00212613
ENST00000567305RHBDD1−0.8608507370.0276972
ENST00000605692ATP6V1F0.8698390770.0243107
ENST00000605692CCDC1360.8293613830.0411927
ENST00000605692FAM71F10.8893623320.0176847
ENST00000605692FAM71F20.9199892470.0093467
ENST00000605692HILPDA−0.9831306710.0004247
ENST00000605692IRF50.959230130.0024597
ENST00000605692LOC1001307050.9542480760.0030927
ENST00000606008ZNF350.8904967930.0173303
ENST00000606008ZNF5010.8477410840.0330093
FR009982LOC100131490−0.861442130.02746712
FR070335ADAMTS150.9151578250.01049211
FR070335ST14−0.9111237610.01149711
FR078172SENP30.8918544640.01691117
FR078172WRAP530.8513421230.03150617
FR085258PIM30.897070960.01534622
FR148984COQ10A−0.9189285340.00959212
FR148984ESYT1−0.8591577920.02835812
FR148984GDF11−0.8612430670.02754412
FR148984MYL6−0.9395088530.00537812
FR148984MYL6B−0.8186641020.04634312
FR148984OBFC2B0.9491213750.00381712
FR148984ORMDL2−0.9386719080.00552612
FR148984RAB5B−0.927133640.00777112
FR171954SEMA4B−0.9012778270.01413815
FR230846SERF1B0.9491449880.0038145
FR314507KCNK6−0.9140524410.01076319
FR314507SIPA1L30.9007712560.01428119
FR314507SPRED30.8858028390.01881719
FR331033EFNB10.9676652130.001551X
NONHSAG002294NTNG10.8358370770.0382121
NONHSAG004747LGALS80.9417448360.0049921
NONHSAG006325PPIF−0.8324714470.03974810
NONHSAG006325ZMIZ10.9446931820.00450410
NONHSAG010348COPS7A0.9845172420.00035812
NONHSAG010348EMG10.9313712880.00690312
NONHSAG010348LRRC230.9007119980.01429812
NONHSAG010348TPI10.9535787340.00318212
NONHSAG010348USP50.9189726590.00958212
NONHSAG010749ARNTL20.9173323690.00996812
NONHSAG025545FXYD50.9244982070.00833619
NONHSAG025545FXYD7−0.9270926420.00777919
NONHSAG025545ZNF1810.8938620230.01630019
NONHSAG025545ZNF300.931364710.00690519
NONHSAG025545ZNF3020.8306846540.04057519
NONHSAG028996IL1RN0.9566441730.0027792
NONHSAG034063SERHL20.9056135040.01294322
NONHSAG035134HEMK1−0.9643517030.0018843
NONHSAG035134MAPKAPK30.8174904810.0469253
NONHSAG037825DCAF4L1−0.8868899570.0184674
NONHSAG039533CYP4V20.9674265740.0015744
NONHSAG039533FAM149A0.9863432790.0002784
NONHSAG041216MAN2A1−0.8631093350.0268265
NONHSAG042317C5orf250.8989352460.0148055
NONHSAG042611FOXQ1−0.9080754680.0122876
NONHSAG047857C7orf420.9509708660.0035477
NONHSAG048587MDFIC0.9421373960.0049257
NONHSAG049128GIMAP4−0.8592952310.0283047
NONHSAG049128GIMAP8−0.9003853810.0143907
NONHSAG049128REPIN1−0.9123147960.0111967
NONHSAG052187DCAF10−0.8865590310.0185739
NONHSAG052187RG9MTD30.9748195060.0009439
NONHSAG052737DAPK10.9916510510.0001049
NONHSAG053520COQ40.8162565510.0475419
NONHSAT002986C1orf1900.9441073170.0045991
NONHSAT002986LOC100133124−0.9235787950.0085371
NONHSAT002986NSUN4−0.8568463020.0292731
NONHSAT002986RAD54L−0.8647947910.0261851
NONHSAT004254CYR610.9905489450.0001341
NONHSAT004989GPSM20.8879363220.0181341
NONHSAT005081CSF10.8709134120.0239201
NONHSAT005942ANKRD34A0.911391660.0114291
NONHSAT005942LIX1L0.9696524040.0013671
NONHSAT005942PDZK10.9825570480.0004541
NONHSAT005942RBM8A0.864276350.0263811
NONHSAT005942RNF115−0.8497412060.0321701
NONHSAT006288C1orf510.8886642510.0179031
NONHSAT006288CA140.8951836690.0159041
NONHSAT006288PRPF30.8128733650.0492481
NONHSAT006799LAMTOR20.8775790690.0215631
NONHSAT006799LMNA−0.9400893030.0052761
NONHSAT006799PMF1-BGLAP−0.8654846740.0259251
NONHSAT006799RXFP4−0.8164285910.0474551
NONHSAT006799SEMA4A−0.8655189450.0259121
NONHSAT006799SLC25A44−0.925298960.0081621
NONHSAT006799TMEM790.8643317260.0263601
NONHSAT007139NDUFS2−0.8115609820.0499181
NONHSAT007139NIT10.8505286270.0318431
NONHSAT007139TOMM40L0.9818497590.0004911
NONHSAT007139UFC10.911611810.0113731
NONHSAT007139USP21−0.9085630940.0121591
NONHSAT010406LGALS80.9201352930.0093131
NONHSAT010413LGALS80.9342685150.0063391
NONHSAT010414LGALS80.946257910.0042551
NONHSAT011117AKR1E20.8874634220.01828410
NONHSAT016347DUSP50.9236831220.00851410
NONHSAT016347MXI1−0.9211674330.00907710
NONHSAT016347RBM200.8160152580.04766210
NONHSAT018220NUCB20.9681149910.00150911
NONHSAT018821LDLRAD30.8305635960.04063111
NONHSAT021952C11orf20−0.9393459320.00540711
NONHSAT021952COX8A−0.8492293230.03238411
NONHSAT021952FERMT30.9084908580.01217811
NONHSAT021952FKBP20.8806146320.02052911
NONHSAT021952NAA400.9504928330.00361611
NONHSAT021952OTUB10.968481330.00147411
NONHSAT021952PLCB3−0.9275224930.00768911
NONHSAT021952PRDX50.8736458550.02293911
NONHSAT021952STIP10.977591670.00074811
NONHSAT021952VEGFB−0.988288610.00020511
NONHSAT023878BIRC20.8790255680.02106711
NONHSAT023878C11orf700.9323412020.00671211
NONHSAT026185COPS7A0.9884909360.00019812
NONHSAT026185EMG10.9421367190.00492512
NONHSAT026185LRRC230.9289252180.00739812
NONHSAT026185PTPN60.8312607810.04030712
NONHSAT026185TPI10.9643452320.00188412
NONHSAT026185USP50.9352013290.00616212
NONHSAT028274SCN8A−0.9261301130.00798412
NONHSAT028274SLC4A80.881166990.02034312
NONHSAT028356C12orf440.9897588660.00015712
NONHSAT030840LOC1001294470.8247541250.04337612
NONHSAT031072HSPB80.9555793280.00291612
NONHSAT034980ITGBL10.9607556470.00228013
NONHSAT035015ERCC50.8143105440.04852013
NONHSAT037446ARG20.9289997590.00738314
NONHSAT037446PLEKHH10.9835622960.00040314
NONHSAT037452ARG20.9374911090.00573914
NONHSAT037452PLEKHH10.9749201480.00093614
NONHSAT051795NME40.8853301030.01897016
NONHSAT051795RAB40C0.8892301120.01772516
NONHSAT051795WFIKKN10.8373413880.03753516
NONHSAT054629LOC1001305800.8758323470.02216917
NONHSAT054629PDK20.8953285420.01586117
NONHSAT054629RSAD1−0.875208440.02238817
NONHSAT054629SPATA200.8687692160.02470217
NONHSAT054629TMEM920.9885060980.00019717
NONHSAT056331CCDC400.9387921210.00550517
NONHSAT056331SLC26A110.8684051340.02483617
NONHSAT060927ANGPTL40.9461176540.00427719
NONHSAT060927HNRNPM0.9132056520.01097319
NONHSAT0609272-Mar0.88683030.01848619
NONHSAT060927OR2Z10.9173178570.00997219
NONHSAT060927RAB11B0.953925370.00313519
NONHSAT066483CEACAM30.9905708710.00013319
NONHSAT066483CEACAM60.9748433120.00094119
NONHSAT074144WBP11−0.9246830380.0082952
NONHSAT075349PHOSPHO20.9189786960.0095812
NONHSAT077950FARP20.8429959750.0350402
NONHSAT077950ING50.9197066940.0094122
NONHSAT077950LOC1001296750.9454350060.0043852
NONHSAT082241BACE20.971768780.00118421
NONHSAT083016COL6A20.842420780.03529021
NONHSAT086623LOC1006527690.9475789270.00405022
NONHSAT086623NOL120.8318853480.04001822
NONHSAT086623SH3BP10.9640242530.00191822
NONHSAT086623TRIOBP0.9565193140.00279522
NONHSAT086830FAM83F0.9350607280.00618922
NONHSAT089680KLHDC8B0.9685504570.0014683
NONHSAT089680TCTA0.8204578860.0454593
NONHSAT090112PXK0.983065330.0004283
NONHSAT092654TM4SF4−0.9007487940.0142873
NONHSAT092847GMPS0.8903005990.0173913
NONHSAT093875HRG−0.9371225940.0058063
NONHSAT093875ST6GAL10.9874027810.0002373
NONHSAT093968CCDC500.8845752490.0192153
NONHSAT093968PYDC20.9633236860.0019933
NONHSAT094657FGFRL1−0.9579806040.0026114
NONHSAT094657SPON20.9250433870.0082174
NONHSAT095682FLJ396530.9248899870.0082504
NONHSAT096163LIMCH10.9725497460.0011204
NONHSAT096163UCHL10.9185398860.0096834
NONHSAT096168DCAF4L1−0.8681945180.0249144
NONHSAT096168LIMCH10.9911592280.0001174
NONHSAT099638CYP4V20.9428165220.0048114
NONHSAT099638FAM149A0.9865474380.0002704
NONHSAT099643CYP4V20.9699966190.0013374
NONHSAT099643FAM149A0.9200094580.0093424
NONHSAT101145PTGER40.8767958640.0218345
NONHSAT105337C5orf250.9469659140.0041445
NONHSAT107900GMPR0.9245771540.0083186
NONHSAT107917CAP20.9690620420.0014216
NONHSAT107973KDM1B0.9461712340.0042686
NONHSAT113449C6orf570.8492483290.0323766
NONHSAT114226PRDM10.9383426310.0055856
NONHSAT119003C7orf260.9877821870.0002237
NONHSAT119003RAC10.9358236090.0060467
NONHSAT119451DNAH110.9852353140.0003257
NONHSAT119452DNAH110.9807860370.0005507
NONHSAT119495KLHL7−0.9170449090.0100377
NONHSAT120737CCT6A0.8986100520.0148997
NONHSAT120737GBAS0.9486866830.0038827
NONHSAT120737MRPS170.8118247690.0497837
NONHSAT120737SUMF20.9482637140.0039467
NONHSAT121617ZP30.8615237020.0274367
NONHSAT122253CPSF40.9084759460.0121827
NONHSAT122253ZNF6550.9491037790.0038207
NONHSAT122253ZNF789−0.944337530.0045617
NONHSAT122254CPSF40.9373484720.0057657
NONHSAT122254ZNF6550.9495405970.0037557
NONHSAT122254ZNF789−0.8410965530.0358697
NONHSAT122256CPSF40.948693280.0038817
NONHSAT122256ZNF6550.9823797850.0004637
NONHSAT122256ZNF789−0.8916131020.0169857
NONHSAT122257CPSF40.9792958850.0006397
NONHSAT122257ZNF789−0.904623090.0132117
NONHSAT122826FOXP20.9845249910.0003577
NONHSAT122828FOXP20.9633072590.0019957
NONHSAT122828MDFIC0.9576291770.0026557
NONHSAT122829FOXP20.9722311290.0011467
NONHSAT122829MDFIC0.8903817530.0173667
NONHSAT122833MDFIC0.9253976620.0081417
NONHSAT122928CAPZA20.9594626870.0024327
NONHSAT122928ST70.988162680.0002097
NONHSAT122929CAPZA20.9530178790.0032597
NONHSAT123242ATP6V1F0.9366903820.0058857
NONHSAT123242CCDC1360.8696480720.0243807
NONHSAT123242HILPDA−0.8433723650.0348777
NONHSAT123242LOC1001307050.8435942450.0347817
NONHSAT125539LOC2548960.8739187420.0228438
NONHSAT125539TNFRSF10C0.9761713320.0008458
NONHSAT127008ADHFE10.8325891520.0396948
NONHSAT127008RRS10.8271255230.0422458
NONHSAT129523LY6E0.8684440370.0248228
NONHSAT129523ZFP410.8758384620.0221678
NONHSAT129523ZNF696−0.8851121690.0190418
NONHSAT134920LCN20.9608496580.0022699
NONHSAT134920ODF20.864503310.0262959
NONHSAT136475SMS−0.9760713450.000852X
NONHSAT136568IL1RAPL10.9767887030.000802X
NONHSAT140499ATF7IP20.8817350510.02015316
NONHSAT140500ATF7IP20.9135328780.01089216
NONHSAT140506ATF7IP20.8393424220.03664316
NONHSAT140507ATF7IP20.8358990630.03818416
NONHSAT142843CCDC113−0.9529641780.00326716
NONHSAT142843NDRG40.9836978040.00039616
NONHSAT142849CCDC113−0.9592336930.00245916
NONHSAT145291TXNDC170.9140402220.01076617
NONHSAT145291XAF10.9903861450.00013817
NR_002834.1OBSCN0.8551425670.0299561
NR_002834.1RNF187−0.8656814160.0258511
NR_027459.2SYT140.8714228940.0237351
NR_027621.1GPR112−0.8186158180.046367X
NR_040072.1CCDC113−0.9751616870.00091816
NR_040072.1NDRG40.9896895730.00015916
NR_046396.1LOC100652883−0.8160839920.04762717
NR_046396.1SLC16A13−0.8525238460.03102017
NR_046396.1TXNDC170.930992770.00697917
NR_046396.1XAF10.979163610.00064717
NR_073032.1KCNK6−0.9454788090.00437819
NR_073032.1SIPA1L30.9427748470.00481819
NR_073032.1SPRED30.9550223550.00298919
NR_102308.1MYADM0.8562058520.02952919
NR_102308.1PRKCG0.9501320070.00366819
NR_102308.1TSEN340.8633723160.02672519
NR_103533.1RUNX2−0.8236103520.0439266
TCONS_00001307SRP90.8816152640.0201931
TCONS_00009873CCT5−0.8741911910.0227465
TCONS_000098736-Mar0.8867012630.0185285
TCONS_00012995ANLN−0.8598837050.0280737
TCONS_00013626GIMAP4−0.8157402820.0478007
TCONS_00013626GIMAP8−0.8571743680.0291427
TCONS_00013626REPIN1−0.9315320210.0068717
TCONS_00026596CCDC1650.8774570330.02160518
TCONS_00026596LOC2842190.8702494390.02416118
TCONS_00027489CYP4F120.9685331690.00147019
TCONS_00027489CYP4F30.9791548060.00064719
TCONS_00027489CYP4F80.8315447040.04017619
TCONS_00027489TPM40.8309646320.04044519
TCONS_l2_00015621PSMD10.8697472460.0243442
TCONS_l2_00021700C4orf290.8732078050.0230954

a P-value of PCC (Pearson correlation coefficient); <0.05 as statistically significant; Chrom, chromosome number.

Trans-regulation of the lncRNAs

Using the threshold of P<0.05 and FDR <0.01, 3,429 lncRNA-TF pairs were found, corresponding to 48 TFs. Then, we generated a core network using the top 100 lncRNA-TF pairs with the most credentiality (lowest P-values and FDRs), as shown by Fig. 4 and are listed in Table IV. We can see that most of these potential trans-regulatory lncRNAs participate in pathways regulated by four TFs: E2F transcription factor 1 (E2F1), E2F transcription factor 4 (E2F4), upstream transcription factor 1 (USF1) and transcription factor AP-2γ (TFAP2C). In the core network of lncRNA-TF pairs, E2F1 participates in 31 of the 100 pairs, E2F4 in 10 pairs, USF1 in 17 pairs and TFAP2C in 30 pairs.

Table IV

The top 100 lncRNA-TF pairs with the most credentiality.

Table IV

The top 100 lncRNA-TF pairs with the most credentiality.

lncRNATFP-valueFDR
NONHSAG029883E2F10.0001001990.011522898
NONHSAT098423E2F10.000109950.012644217
FR165245E2F10.0001149540.013219731
NONHSAT001522E2F40.000119770.005932082
TCONS_l2_00015621YY10.0001202820.013591822
NONHSAG035406E2F10.0001205590.013623208
NONHSAT005942E2F10.0001335230.015355119
NONHSAT008438USF10.0001348040.007751239
NONHSAT030840E2F40.0001362620.010899038
NONHSAG050806TFAP2C0.0001408530.008099028
NONHSAT122253TFAP2C0.0001411070.008113627
NONHSAT011117USF10.000141860.008156953
NONHSAT093783E2F10.0001488330.016966919
ENST00000508616E2F10.0001508610.017047245
TCONS_00001679ATF30.0001534480.005882192
NONHSAT001522E2F10.000154750.005932082
NONHSAT081186USF10.0001631550.009218242
NR_027459.2E2F10.0001650250.018977858
ENST00000420044TFAP2C0.0001652340.009500979
NONHSAT092847E2F10.0001702330.00970326
FR009982USF10.0001744820.02006547
ENST00000598393E2F10.0001756420.010099417
NONHSAT028356USF10.0001826840.010504343
ENST00000438810TFAP2C0.0001848360.010628064
NONHSAT030840E2F10.0001895480.010899038
NONHSAT095682TFAP2C0.0001918110.01102914
NONHSAT129523ELK40.0001921260.022094508
NONHSAT010620E2F10.0001945950.02198926
ENST00000533220E2F40.0001951460.011123299
NONHSAT017390E2F10.0002015510.022372179
NONHSAG001064E2F10.0002017770.022800789
NONHSAT140499TFAP2C0.000201950.011612109
FR042098USF10.0002023390.011634494
TCONS_l2_00021700PPARGC1A0.0002168160.024283375
NONHSAT120476TFAP2C0.000223880.012761147
TCONS_00016368E2F40.0002438190.027795322
NONHSAT037766ELK40.0002567040.017378335
NONHSAG031748USF10.0002647720.023070091
FR224199E2F10.0002662610.030619965
NONHSAG003987TFAP2C0.0002750490.015815294
TCONS_l2_00020217E2F10.0002765530.031250531
NONHSAT090951E2F40.0002832690.032575881
NONHSAT138382E2F10.0002912240.032617038
NONHSAT031072TFAP2C0.000299540.017223562
NONHSAT118621TFAP2C0.000299850.017091435
NONHSAT037766E2F10.0003022320.017378335
NONHSAT122928TFAP2C0.0003078540.017701585
NONHSAT026185TFAP2C0.0003081280.017717359
ENST00000412387E2F40.0003109840.035763147
ENST00000602721TFAP2C0.0003112380.017896208
ENST00000434951TFAP2C0.0003189880.020166596
NR_002834.1USF10.0003193640.018363452
NONHSAT030171USF10.000321510.036652133
NONHSAT035083E2F10.0003250070.03672582
NONHSAT060965TFAP2C0.000329510.018946816
NONHSAT093783E2F40.0003349030.019089467
TCONS_00013626E2F10.000340060.039106915
NONHSAG039347E2F10.0003436880.039524106
FR314507ZBTB7A0.0003481420.020018175
NONHSAT022221PPARGC1A0.0003489480.019715575
ENST00000434951TFAP2A0.00035380.020166596
NR_073125.1TFAP2C0.0003551860.020068009
FR070335E2F10.0003602670.041430707
NONHSAT026990USF10.0003759350.021616259
NONHSAT136810TFAP2C0.000376710.021660808
NONHSAT010513E2F10.0003969640.045650912
TCONS_00029195TFAP2C0.0003980870.015818714
NONHSAT119003TFAP2C0.0004032530.022985416
FR055729E2F10.0004037170.0464274
ENST00000417120E2F10.0004064920.045527124
NONHSAG031748E2F40.000408320.023070091
NONHSAT130274TFAP2C0.000408650.023497401
TCONS_00029195TFAP2A0.0004126620.015818714
NONHSAT096163E2F40.0004248830.024218343
NONHSAT056714TFAP2C0.0004259790.024493803
NONHSAG037825USF10.0004277440.024595276
NR_104426.1TFAP2C0.0004312610.024581897
NONHSAT122833TFAP2C0.0004366450.025107059
TCONS_00027489TFAP2C0.0004420760.025419365
NONHSAT093875USF10.0004436360.02550908
NONHSAT140507E2F10.000452960.030884797
NONHSAT023895USF10.0004542180.051326598
NONHSAT077173ATF30.0004578420.026325887
NONHSAT093946TFAP2C0.0004590050.026392802
NONHSAT119495E2F10.0004664420.053640826
NONHSAT136475E2F40.0004664810.026589441
NONHSAT096168TFAP2C0.0004669480.026616051
NONHSAG046899USF10.0004713520.026867058
NONHSAG050557E2F10.0004720680.05428779
NONHSAT032403E2F10.0004781040.054025738
ENST00000561486TFAP2C0.0004794310.027567279
TCONS_l2_00015621 USF10.0004874680.02320976
NONHSAT006799TFAP2C0.0004922180.028302546
NONHSAT034980E2F10.0004935720.056760748
NONHSAT004232TFAP2C0.0005051120.028791388
NONHSAT079625USF10.0005169610.058933596
NONHSAT095682TFAP2A0.0005256170.020148659
NONHSAT077174TFAP2C0.0005258790.024031165
FR280917USF10.000526220.029994527
TCONS_00005559TFAP2A0.0005346750.030743841

[i] P-value was calculated by hypergeometric cumulative distribution function with P<0.05 as statistically significant and FDR calculated under the control of the Benjamini and Hochberg procedure, P<0.05 as statistically significant. These 100 pairs were the top ranking ones with the least P-values and FDRs among all the lncRNA-TF pairs, thus with the highest credentiality. TF, transcription factor.

Discussion

We assessed genome-wide lncRNA expression patterns in gefitinib-resistant human NSCLC HCC827-8-1 cells, compared with their parental HCC827 cells by microarray and explored their possible functions by analyzing their co-expressed protein-coding mRNAs. HCC827-8-1 is an individual clone isolated from resistant HCC827 cells, which is 348-fold more resistant to gefitinib than HCC827 cells with high stability (11). We used HCC827-8-1 cells to study the EGFR-TKI resistance mechanisms and to obtain more accurate results.

The microarray showed that 1,476 lncRNA and 1,026 mRNA transcripts were dysregulated. Then we chose and validated 7 of the differentially expressed lncRNAs by RT-qPCR. Hundreds of lncRNAs were co-expressed with thousands of mRNAs, and some of the lncRNAs may be active in the mechanism of gefitinib resistance through cis-regulation and/or in trans-regulation of these mRNAs.

We predicted the roles of lncRNAs through the co-expressed mRNAs. The most widely used technique for function prediction is the 'guilt by association' method, which aims to interrogate the co-expressed protein-coding genes and relevant bio-pathways (13). Among the hundreds of pathways as-predicted, some are key pathways in the mechanisms of gefitinib resistance, such as focal adhesion, cell proliferation, cell cycle and apoptosis (1618). We compared the differentially expressed mRNAs with other studies, and the high consistency strongly supported our predictions. For example, in the mRNA expression analysis using Agilent SurePrint G3 Human Gene Expression 8×60K array of a gefitinib-resistant lung cancer cell line (19), IGFBP3 was found to be expressed differentially, which is in accordance with our results. Additionally, the constitutive activation of EGFR tyrosine kinase receptor caused by mutation and/or amplification is closely correlated with the development, progression and poor prognosis of NSCLC, through the activation of various downstream signaling pathways. These pathways include the RAS/RAF/MAPK pathway that induces signals associated with proliferative activity, and the PI3K/AKT pathway that results in an anti-apoptotic effect (16). In the present study, the expression levels of AKT and RAS were higher in the HCC827-8-1 cells than that in the HCC827 cells.

cis-regulation is the intrinsic property of ncRNAs, since mRNAs are only effective in dissociation, transportation and translation (5). Moreover, the cis-regulation of nearby protein-coding genes is a common mechanism of lncRNAs (20). In the present study, we identified 149 lncRNA transcripts that were able to cis-regulate the nearby protein-coding genes. Although most of the 149 lncRNAs were uncharacterized, several of the lncRNA-mRNAs attracted our interest. lncRNAs that are related to dysregulated mRNAs known to be involved in EGFR-TKI resistance are valuable for further study. For example, one lncRNA is NONHSAT021952 (CUST_7952_PI429545402), which probably cis-regulates VEGF-B mRNA. As reported, VEGF-B plays important roles in promoting cancer metastasis through the remodeling of microvasculature of NSCLC, which is related to EGFR-TKI resistance (21). NONHSAT021952 (chr11: 63988306-63990592+) has no annotated functions and lies near VEGF-B (chr11: 64006111-64006170+) in the same chromosome. Thus, there may be a new way that nearby lncRNAs regulate the expression of VEGF-B.

Alhough some lncRNAs are proven as cis-regulators, the vast majority of lncRNAs with known functions are actually trans-regulators (22). In contrast to cis-regulation, the trans-regulatory lncRNAs, after isolation from the synthesis site, affect the genes at a long genomic distance away, accumulate and take effect in large regulatory networks (5). We predicted the functions of trans-regulatory lncRNAs through the relevant expression-regulating TFs. In the core network of lncRNA-TF pairs, the lncRNAs were divided into four groups of pathways that are controlled by E2F4, E2F1, USF1 and TFAP2C, respectively.

E2F4 and E2F1 which belong to the E2F family are proposed to be primary transcription repressors. The E2F family is pivotal in the regulation of cell cycle and the action of tumor suppressor proteins. The E2F family is also a target of the transforming proteins of small DNA tumor viruses (23).

E2F1 regulates the transcription of DNA synthesis-related genes and thus, the binding of hypophosphorylated pRB to E2F1 will arrest G1 phase cells (24). Gefitinib treatment downregulates the expression of E2F1 mRNA and protein. The antiproliferative effect of gefitinib is attributed, or at least partially, to the inhibition of E2F1 expression and increases the proportion of G1 phase cells in human LACs (A431 and A549) (25). E2F1 is also downregulated by gefitinib treatment in gefitinib-resistant NSCLC cells with MET amplification, but not in resistant NSCLC cells harboring the T790M mutation of EGFR (26). The above studies indicate that E2F1 is closely associated with the resistance to EGFR-TKIs. The lncRNAs predicted as regulatory targets of E2F1 may take part in the acquisition of gifitinib resistance, which should be validated by further functional tests.

E2F4 takes part in diverse functions, such as cell cycle control, DNA damage repair, apoptosis, mRNA processing and ubiquitination (27,28). The upregulation of E2F4 nuclear expression in breast cancer is related to poor prognosis in patients with larger tumors, recurrence, distant metastasis and poorer outcome (29). E2F4 protein expression in human colorectal tumors is also upregulated by affecting the G1/S phase transition and the proliferation of colon cancer cells (30). In comparison, E2F4 protein expression is down-regulated in sporadic Burkitt lymphoma (sBL) cells but not in immortal B-cells. Besides the ability to reduce E2F1 levels, E2F4 expression in BL cells also selectively weakens the tumorigenic properties and BL cell proliferation and increases the proportions of G2/M cells (31). Since E2F4 has diverse effects in the regulation of tumor progression, it may also play a role in lung cancer.

USF1 is a basic helix-loop-helix (bHLH) TF encoded by distinct genes that are heterodimerized and functionally overlapped (32). USF1 also mediates the transcription of many neuropeptides (3336), surfactant protein A in the lung (37), genes in lung tumors (38) which include hTERT as the immortalizing telomerase subunit (39,40) and CDK4 as the cell cycle regulator (41). USF1 can also regulate neuropeptide genes in lung cancer, particularly arginine vasopressin (AVP) (42,43). USF1 mRNA expression is downregulated in LACs compared with that in normal tissues (44). In lung cancer, USF1 and USF2 dimerize to mediate the transcription via E-box motifs in target genes (45). Moreover, with the E-box motif, USF1 activates GATA5 gene expression (46), which is pivotal in cell differentiation and tissue development in lung cancer (47,48). Since USF1 is closely related to lung cancer occurrence, USF1-regulated lncRNAs may be involved in EGFR-TKI resistance.

TFAP2C, also known as AP2-γ, belongs to the AP-2 TF family, which is composed of five closely related 50-kDa isoforms (49). By regulating the expression of many downstream genes, TFAP2C mediates various cellular processes, including induction, differentiation, survival, proliferation and apoptosis under diverse developmental contexts (5054). The internalization of activated EGFR involves the AP-2 complex (55). TFAP2C is pivotal in human cancer development. For example, through transcriptional regulation of EGFR, TFAP2C mediates the tumor development, cell growth and survival of HER2-amplified breast cancer (56). Moreover, after 10 years of diagnosis, higher TFAP2C scores are correlated with poor overall survival in ERα-positive and endocrine therapy-treated patients (57). Since TFAP2C affects tumorigenesis by mediating EGFR expression, TFAP2C may be an important cause of EGFR-TKI-resistant lung cancer. Thus, the lncRNAs that take part in the TFAP2C-regulated pathways, as predicted, are candidate participants in the mechanisms of EGFR-TKI resistance.

The present study has also a few limitations. Firstly, we only used one cell line to study the mechanisms of EGFR-TKI resistance in vitro. Our results are not integrated but typical, and thus further studies with larger sample size in vivo are needed. Second, we only predicted lncRNA functions indirectly through network and pathway analyses of co-expressed protein-coding genes. The reason is that the majority of the lncRNAs identified to date are not functionally characterized (58). The 'guilt by association'' method for hypothesis generation is a key primary step for functional studies in the future, such as loss/acquisition of functions (59).

In conclusion, nearly 1,500 lncRNAs were differently expressed between the gefitinib-sensitive HCC827 cells and the gefitinib-resistant HCC827-8-1 cells. Many of these lncRNAs play important roles in regulating EGFR-TKI resistance through the cis- and/or trans-regulation of target protein-coding genes. The present study underlies future functional studies on lncRNAs related to EGFR resistance and provides a candidate reservoir for their application as diagnostic and therapeutic targets.

Acknowledgments

The present study was supported by grants from the National Natural Science Foundation of China (81372396).

References

1 

Jackman DM, Miller VA, Cioffredi LA, Yeap BY, Jänne PA, Riely GJ, Ruiz MG, Giaccone G, Sequist LV and Johnson BE: Impact of epidermal growth factor receptor and KRAS mutations on clinical outcomes in previously untreated non-small cell lung cancer patients: Results of an online tumor registry of clinical trials. Clin Cancer Res. 15:5267–5273. 2009. View Article : Google Scholar : PubMed/NCBI

2 

Godin-Heymann N, Bryant I, Rivera MN, Ulkus L, Bell DW, Riese DJ II, Settleman J and Haber DA: Oncogenic activity of epidermal growth factor receptor kinase mutant alleles is enhanced by the T790M drug resistance mutation. Cancer Res. 67:7319–7326. 2007. View Article : Google Scholar : PubMed/NCBI

3 

McDermott U, Pusapati RV, Christensen JG, Gray NS and Settleman J: Acquired resistance of non-small cell lung cancer cells to MET kinase inhibition is mediated by a switch to epidermal growth factor receptor dependency. Cancer Res. 70:1625–1634. 2010. View Article : Google Scholar : PubMed/NCBI

4 

Yamamoto C, Basaki Y, Kawahara A, Nakashima K, Kage M, Izumi H, Kohno K, Uramoto H, Yasumoto K, Kuwano M, et al: Loss of PTEN expression by blocking nuclear translocation of EGR1 in gefitinib-resistant lung cancer cells harboring epidermal growth factor receptor-activating mutations. Cancer Res. 70:8715–8725. 2010. View Article : Google Scholar : PubMed/NCBI

5 

Lee JT: Epigenetic regulation by long noncoding RNAs. Science. 338:1435–1439. 2012. View Article : Google Scholar : PubMed/NCBI

6 

Qiu M, Xu Y, Yang X, Wang J, Hu J, Xu L and Yin R: CCAT2 is a lung adenocarcinoma-specific long non-coding RNA and promotes invasion of non-small cell lung cancer. Tumour Biol. 35:5375–5380. 2014. View Article : Google Scholar : PubMed/NCBI

7 

Gutschner T, Hämmerle M, Eissmann M, Hsu J, Kim Y, Hung G, Revenko A, Arun G, Stentrup M, Gross M, et al: The noncoding RNA MALAT1 is a critical regulator of the metastasis phenotype of lung cancer cells. Cancer Res. 73:1180–1189. 2013. View Article : Google Scholar :

8 

Liu XH, Liu ZL, Sun M, Liu J, Wang ZX and De W: The long non-coding RNA HOTAIR indicates a poor prognosis and promotes metastasis in non-small cell lung cancer. BMC Cancer. 13:4642013. View Article : Google Scholar : PubMed/NCBI

9 

Liu Z1, Sun M, Lu K, Liu J, Zhang M, Wu W, De W, Wang Z and Wang R: The long noncoding RNA HOTAIR contributes to cisplatin resistance of human lung adenocarcinoma cells via downregualtion of p21WAF1/CIP1 expression. PLoS One. 8:e772932013. View Article : Google Scholar

10 

Dong S, Qu X, Li W, Zhong X, Li P, Yang S, Chen X, Shao M and Zhang L: The long non-coding RNA, GAS5, enhances gefitinib-induced cell death in innate EGFR tyrosine kinase inhibitor-resistant lung adenocarcinoma cells with wide-type EGFR via downregulation of the IGF-1R expression. J Hematol Oncol. 8:432015. View Article : Google Scholar : PubMed/NCBI

11 

Wu Y, Yu DD, Hu Y, Cao HX, Yu SR, Liu SW and Feng JF: LXR ligands sensitize EGFR-TKI-resistant human lung cancer cells in vitro by inhibiting Akt activation. Biochem Biophys Res Commun. 467:900–905. 2015. View Article : Google Scholar : PubMed/NCBI

12 

Livak KJ and Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2−ΔΔCT method. Methods. 25:402–408. 2001. View Article : Google Scholar

13 

Guttman M, Amit I, Garber M, French C, Lin MF, Feldser D, Huarte M, Zuk O, Carey BW, Cassady JP, et al: Chromatin signature reveals over a thousand highly conserved large non-coding RNAs in mammals. Nature. 458:223–227. 2009. View Article : Google Scholar : PubMed/NCBI

14 

Efron B and Tibshirani R: Empirical bayes methods and false discovery rates for microarrays. Genet Epidemiol. 23:70–86. 2002. View Article : Google Scholar : PubMed/NCBI

15 

Guttman M, Donaghey J, Carey BW, Garber M, Grenier JK, Munson G, Young G, Lucas AB, Ach R, Bruhn L, et al: lincRNAs act in the circuitry controlling pluripotency and differentiation. Nature. 477:295–300. 2011. View Article : Google Scholar : PubMed/NCBI

16 

Remon J, Morán T, Majem M, Reguart N, Dalmau E, Márquez-Medina D and Lianes P: Acquired resistance to epidermal growth factor receptor tyrosine kinase inhibitors in EGFR-mutant non-small cell lung cancer: A new era begins. Cancer Treat Rev. 40:93–101. 2014. View Article : Google Scholar

17 

Li LH, Wu P, Lee JY, Li PR, Hsieh WY, Ho CC, Ho CL, Chen WJ, Wang CC, Yen MY, et al: Hinokitiol induces DNA damage and autophagy followed by cell cycle arrest and senescence in gefitinib-resistant lung adenocarcinoma cells. PLoS One. 9:e1042032014. View Article : Google Scholar : PubMed/NCBI

18 

Ju L and Zhou C: Association of integrin beta1 and c-MET in mediating EGFR TKI gefitinib resistance in non-small cell lung cancer. Cancer Cell Int. 13:152013. View Article : Google Scholar : PubMed/NCBI

19 

Terai H, Soejima K, Yasuda H, Sato T, Naoki K, Ikemura S, Arai D, Ohgino K, Ishioka K, Hamamoto J, et al: Long-term exposure to gefitinib induces acquired resistance through DNA methylation changes in the EGFR-mutant PC9 lung cancer cell line. Int J Oncol. 46:430–436. 2015.

20 

Ørom UA, Derrien T, Beringer M, Gumireddy K, Gardini A, Bussotti G, Lai F, Zytnicki M, Notredame C, Huang Q, et al: Long noncoding RNAs with enhancer-like function in human cells. Cell. 143:46–58. 2010. View Article : Google Scholar : PubMed/NCBI

21 

Yang X, Zhang Y, Hosaka K, Andersson P, Wang J, Tholander F, Cao Z, Morikawa H, Tegnér J, Yang Y, et al: VEGF-B promotes cancer metastasis through a VEGF-A-independent mechanism and serves as a marker of poor prognosis for cancer patients. Proc Natl Acad Sci USA. 112:E2900–E2909. 2015. View Article : Google Scholar : PubMed/NCBI

22 

Guttman M and Rinn JL: Modular regulatory principles of large non-coding RNAs. Nature. 482:339–346. 2012. View Article : Google Scholar : PubMed/NCBI

23 

Vaishnav YNV, Vaishnav MY and Pant V: The molecular and functional characterization of E2F-5 transcription factor. Biochem Biophys Res Commun. 242:586–592. 1998. View Article : Google Scholar : PubMed/NCBI

24 

Stewart ZA, Westfall MD and Pietenpol JA: Cell-cycle dysregulation and anticancer therapy. Trends Pharmacol Sci. 24:139–145. 2003. View Article : Google Scholar : PubMed/NCBI

25 

Suenaga M, Yamaguchi A, Soda H, Orihara K, Tokito Y, Sakaki Y, Umehara M, Terashi K, Kawamata N, Oka M, et al: Antiproliferative effects of gefitinib are associated with suppression of E2F-1 expression and telomerase activity. Anticancer Res. 26:3387–3391. 2006.PubMed/NCBI

26 

Okabe T, Okamoto I, Tsukioka S, Uchida J, Hatashita E, Yamada Y, Yoshida T, Nishio K, Fukuoka M, Jänne PA, et al: Addition of S-1 to the epidermal growth factor receptor inhibitor gefitinib overcomes gefitinib resistance in non-small cell lung cancer cell lines with MET amplification. Clin Cancer Res. 15:907–913. 2009. View Article : Google Scholar : PubMed/NCBI

27 

Lindeman GJ, Gaubatz S, Livingston DM and Ginsberg D: The subcellular localization of E2F-4 is cell-cycle dependent. Proc Natl Acad Sci USA. 94:5095–5100. 1997. View Article : Google Scholar : PubMed/NCBI

28 

Mikkelsen TS, Ku M, Jaffe DB, Issac B, Lieberman E, Giannoukos G, Alvarez P, Brockman W, Kim TK, Koche RP, et al: Genome-wide maps of chromatin state in pluripotent and lineage-committed cells. Nature. 448:553–560. 2007. View Article : Google Scholar : PubMed/NCBI

29 

Rakha EA, Pinder SE, Paish EC, Robertson JF and Ellis IO: Expression of E2F-4 in invasive breast carcinomas is associated with poor prognosis. J Pathol. 203:754–761. 2004. View Article : Google Scholar : PubMed/NCBI

30 

Garneau H, Paquin MC, Carrier JC and Rivard N: E2F4 expression is required for cell cycle progression of normal intestinal crypt cells and colorectal cancer cells. J Cell Physiol. 221:350–358. 2009. View Article : Google Scholar : PubMed/NCBI

31 

Molina-Privado I, Jiménez-P R, Montes-Moreno S, Chiodo Y, Rodríguez-Martínez M, Sánchez-Verde L, Iglesias T, Piris MA and Campanero MR: E2F4 plays a key role in Burkitt lymphoma tumorigenesis. Leukemia. 26:2277–2285. 2012. View Article : Google Scholar : PubMed/NCBI

32 

Sirito M, Lin Q, Deng JM, Behringer RR and Sawadogo M: Overlapping roles and asymmetrical cross-regulation of the USF proteins in mice. Proc Natl Acad Sci USA. 95:3758–3763. 1998. View Article : Google Scholar : PubMed/NCBI

33 

Nielsen FC, Pedersen K, Hansen TV, Rourke IJ and Rehfeld JF: Transcriptional regulation of the human cholecystokinin gene: Composite action of upstream stimulatory factor, Sp1, and members of the CREB/ATF-AP-1 family of transcription factors. DNA Cell Biol. 15:53–63. 1996. View Article : Google Scholar : PubMed/NCBI

34 

Viney TJ, Schmidt TW, Gierasch W, Sattar AW, Yaggie RE, Kuburas A, Quinn JP, Coulson JM and Russo AF: Regulation of the cell-specific calcitonin/calcitonin gene-related peptide enhancer by USF and the Foxa2 forkhead protein. J Biol Chem. 279:49948–49955. 2004. View Article : Google Scholar : PubMed/NCBI

35 

Paterson JM, Morrison CF, Mendelson SC, McAllister J and Quinn JP: An upstream stimulatory factor (USF) binding motif is critical for rat preprotachykinin-A promoter activity in PC12 cells. Biochem J. 310:401–406. 1995. View Article : Google Scholar : PubMed/NCBI

36 

Hadsell DL, Bonnette S, George J, Torres D, Klimentidis Y, Gao S, Haney PM, Summy-Long J, Soloff MS, Parlow AF, et al: Diminished milk synthesis in upstream stimulatory factor 2 null mice is associated with decreased circulating oxytocin and decreased mammary gland expression of eukaryotic initiation factors 4E and 4G. Mol Endocrinol. 17:2251–2267. 2003. View Article : Google Scholar : PubMed/NCBI

37 

Gao E, Wang Y, Alcorn JL and Mendelson CR: Transcription factor USF2 is developmentally regulated in fetal lung and acts together with USF1 to induce SP-A gene expression. Am J Physiol Lung Cell Mol Physiol. 284:L1027–L1036. 2003. View Article : Google Scholar : PubMed/NCBI

38 

Coulson JM, Edgson JL, Marshall-Jones ZV, Mulgrew R, Quinn JP and Woll PJ: Upstream stimulatory factor activates the vasopressin promoter via multiple motifs, including a non-canonical E-box. Biochem J. 369:549–561. 2003. View Article : Google Scholar

39 

McMurray HR and McCance DJ: Human papillomavirus type 16 E6 activates TERT gene transcription through induction of c-Myc and release of USF-mediated repression. J Virol. 77:9852–9861. 2003. View Article : Google Scholar : PubMed/NCBI

40 

Goueli BS and Janknecht R: Regulation of telomerase reverse transcriptase gene activity by upstream stimulatory factor. Oncogene. 22:8042–8047. 2003. View Article : Google Scholar : PubMed/NCBI

41 

Pawar SA, Szentirmay MN, Hermeking H and Sawadogo M: Evidence for a cancer-specific switch at the CDK4 promoter with loss of control by both USF and c-Myc. Oncogene. 23:6125–6135. 2004. View Article : Google Scholar : PubMed/NCBI

42 

Coulson JM, Fiskerstrand CE, Woll PJ and Quinn JP: E-box motifs within the human vasopressin gene promoter contribute to a major enhancer in small-cell lung cancer. Biochem J. 344:961–970. 1999. View Article : Google Scholar : PubMed/NCBI

43 

Grace CO, Fink G and Quinn JP: Characterization of potential regulatory elements within the rat arginine vasopressin proximal promoter. Neuropeptides. 33:81–90. 1999. View Article : Google Scholar

44 

Khattar NH, Lele SM and Kaetzel CS: Down-regulation of the polymeric immunoglobulin receptor in non-small cell lung carcinoma: Correlation with dysregulated expression of the transcription factors USF and AP2. J Biomed Sci. 12:65–77. 2005. View Article : Google Scholar : PubMed/NCBI

45 

Ocejo-Garcia M, Baokbah TA, Ashurst HL, Cowlishaw D, Soomro I, Coulson JM and Woll PJ: Roles for USF-2 in lung cancer proliferation and bronchial carcinogenesis. J Pathol. 206:151–159. 2005. View Article : Google Scholar : PubMed/NCBI

46 

Chen B, Hsu R, Li Z, Kogut PC, Du Q, Rouser K, Camoretti-Mercado B and Solway J: Upstream stimulatory factor 1 activates GATA5 expression through an E-box motif. Biochem J. 446:89–98. 2012. View Article : Google Scholar : PubMed/NCBI

47 

Kakita T, Hasegawa K, Morimoto T, Kaburagi S, Wada H and Sasayama S: p300 protein as a coactivator of GATA-5 in the transcription of cardiac-restricted atrial natriuretic factor gene. J Biol Chem. 274:34096–34102. 1999. View Article : Google Scholar : PubMed/NCBI

48 

Singh MK, Li Y, Li S, Cobb RM, Zhou D, Lu MM, Epstein JA, Morrisey EE and Gruber PJ: Gata4 and Gata5 cooperatively regulate cardiac myocyte proliferation in mice. J Biol Chem. 285:1765–1772. 2010. View Article : Google Scholar :

49 

Tummala R, Romano RA, Fuchs E and Sinha S: Molecular cloning and characterization of AP-2 epsilon, a fifth member of the AP-2 family. Gene. 321:93–102. 2003. View Article : Google Scholar : PubMed/NCBI

50 

Hoffman TL, Javier AL, Campeau SA, Knight RD and Schilling TF: Tfap2 transcription factors in zebrafish neural crest development and ectodermal evolution. J Exp Zoolog B Mol Dev Evol. 308:679–691. 2007. View Article : Google Scholar

51 

Kuckenberg P, Kubaczka C and Schorle H: The role of transcription factor Tcfap2c/TFAP2C in trophectoderm development. Reprod Biomed Online. 25:12–20. 2012. View Article : Google Scholar : PubMed/NCBI

52 

Li W and Cornell RA: Redundant activities of Tfap2a and Tfap2c are required for neural crest induction and development of other non-neural ectoderm derivatives in zebrafish embryos. Dev Biol. 304:338–354. 2007. View Article : Google Scholar : PubMed/NCBI

53 

Li X, Glubrecht DD and Godbout R: AP2 transcription factor induces apoptosis in retinoblastoma cells. Genes Chromosomes Cancer. 49:819–830. 2010.PubMed/NCBI

54 

Van Otterloo E, Li W, Garnett A, Cattell M, Medeiros DM and Cornell RA: Novel Tfap2-mediated control of soxE expression facilitated the evolutionary emergence of the neural crest. Development. 139:720–730. 2012. View Article : Google Scholar : PubMed/NCBI

55 

Rappoport JZ and Simon SM: Endocytic trafficking of activated EGFR is AP-2 dependent and occurs through preformed clathrin spots. J Cell Sci. 122:1301–1305. 2009. View Article : Google Scholar : PubMed/NCBI

56 

Park JM, Wu T, Cyr AR, Woodfield GW, De Andrade JP, Spanheimer PM, Li T, Sugg SL, Lal G, Domann FE, et al: The role of Tcfap2c in tumorigenesis and cancer growth in an activated Neu model of mammary carcinogenesis. Oncogene. 34:6105–6104. 2015. View Article : Google Scholar : PubMed/NCBI

57 

Perkins SM, Bales C, Vladislav T, Althouse S, Miller KD, Sandusky G, Badve S and Nakshatri H: TFAP2C expression in breast cancer: Correlation with overall survival beyond 10 years of initial diagnosis. Breast Cancer Res Treat. 152:519–531. 2015. View Article : Google Scholar : PubMed/NCBI

58 

Derrien T, Johnson R, Bussotti G, Tanzer A, Djebali S, Tilgner H, Guernec G, Martin D, Merkel A, Knowles DG, et al: The GENCODE v7 catalog of human long noncoding RNAs: Analysis of their gene structure, evolution, and expression. Genome Res. 22:1775–1789. 2012. View Article : Google Scholar : PubMed/NCBI

59 

Rinn JL and Chang HY: Genome regulation by long noncoding RNAs. Annu Rev Biochem. 81:145–166. 2012. View Article : Google Scholar : PubMed/NCBI

Related Articles

Journal Cover

June-2016
Volume 35 Issue 6

Print ISSN: 1021-335X
Online ISSN:1791-2431

Sign up for eToc alerts

Recommend to Library

Copy and paste a formatted citation
x
Spandidos Publications style
Wu Y, Yu D, Hu Y, Yan D, Chen X, Cao HX, Yu S, Wang Z and Feng J: Genome-wide profiling of long non-coding RNA expression patterns in the EGFR-TKI resistance of lung adenocarcinoma by microarray. Oncol Rep 35: 3371-3386, 2016
APA
Wu, Y., Yu, D., Hu, Y., Yan, D., Chen, X., Cao, H. ... Feng, J. (2016). Genome-wide profiling of long non-coding RNA expression patterns in the EGFR-TKI resistance of lung adenocarcinoma by microarray. Oncology Reports, 35, 3371-3386. https://doi.org/10.3892/or.2016.4758
MLA
Wu, Y., Yu, D., Hu, Y., Yan, D., Chen, X., Cao, H., Yu, S., Wang, Z., Feng, J."Genome-wide profiling of long non-coding RNA expression patterns in the EGFR-TKI resistance of lung adenocarcinoma by microarray". Oncology Reports 35.6 (2016): 3371-3386.
Chicago
Wu, Y., Yu, D., Hu, Y., Yan, D., Chen, X., Cao, H., Yu, S., Wang, Z., Feng, J."Genome-wide profiling of long non-coding RNA expression patterns in the EGFR-TKI resistance of lung adenocarcinoma by microarray". Oncology Reports 35, no. 6 (2016): 3371-3386. https://doi.org/10.3892/or.2016.4758