Epigallocatechin‑3‑Gallate Modulates Long Non‑Coding RNA and Mrna Expression Profiles in Lung Cancer Cells
Total Page:16
File Type:pdf, Size:1020Kb
MOLECULAR MEDICINE REPORTS 19: 1509-1520, 2019 Epigallocatechin‑3‑gallate modulates long non‑coding RNA and mRNA expression profiles in lung cancer cells DONG-LI HU1,2, GUO WANG1,2, JING YU1,2, LI-HUA ZHANG1,2, YUAN-FEI HUANG1,2, DAN WANG1,2 and HONG-HAO ZHOU1,2 1Department of Clinical Pharmacology, Xiangya Hospital, Central South University, Changsha, Hunan 410008; 2Hunan Key Laboratory of Pharmacogenetics, Institute of Clinical Pharmacology, Central South University, Changsha, Hunan 410078, P.R. China Received March 12, 2018; Accepted October 26, 2018 DOI: 10.3892/mmr.2019.9816 Abstract. (-)-Epigallocatechin-3-gallate (EGCG), a Introduction major constituent of green tea, is a potential anticancer agent, but the molecular mechanisms of its effects are Lung cancer is one of the most common solid tumors and has not well-understood. The present study was conducted to very high global incidence and mortality rates (1). Although examine the mechanism of EGCG in lung cancer cells. therapeutic treatment of lung cancer has made great progress, Alterations in long non-coding RNAs (lncRNAs) and the prognosis of lung cancer patients remains unsatisfactory and mRNAs were investigated in lung cancer cells treated treatment causes a number of side effects (2). Natural dietary with EGCG by lncRNA microarray analysis. Furthermore, molecules are promising candidates for novel therapeutic agents the functions and signaling pathways regulated by EGCG in cancer. Green tea is one of the most widely consumed bever- were predicted by bioinformatics analysis. A total of ages, worldwide (3). Furthermore, (-)-epigallocatechin-3-gallate 960 lncRNAs and 1,434 mRNAs were significantly (EGCG), the main component of green tea polyphenols, was altered following EGCG treatment. These lncRNAs were demonstrated to have various biological activities and demon- distributed across nearly all human chromosomes and the strates potential as a chemical and therapeutic agent in several mRNAs were involved in the cell cycle and the mitotic cell diseases, including various cancers (3-8) and atherosclerosis (9). cycle process. Through a combination of microarray and Long non-coding RNAs (lncRNAs) are defined as a bioinformatics analysis, 20 mRNAs predicted to serve a key class of transcripts longer than 200 nucleotides (nt) in the role in the EGCG regulation were identified, and certain cytoplasm and nucleus (10-13). LncRNAs, characterized regulatory networks involving EGCG-regulated lncRNAs by the complexity and diversity of their sequences, have were predicted. In conclusion, EGCG affects the expression been implicated in a wide spectrum of cellular activities of various lncRNAs and mRNAs in the cells, therefore and diseases (14,15). Emerging studies have demonstrated affecting cell functions. The results of the present study that certain lncRNAs are frequently abnormally regulated provide an insight into the mechanism of EGCG, which may in several types of cancer and serve important roles in the be useful for therapeutic development. occurrence and development of tumors (16-19). Accumulating evidence has confirmed the activity of EGCG in lung cancer therapy and many lncRNAs and mRNAs are involved in regu- lating tumorigenesis (10,20-22). However, microarray analysis Correspondence to: Dr Guo Wang, Department of Clinical of the differential expression profiles of lncRNAs and mRNAs Pharmacology, Xiangya Hospital, Central South University, in lung cancer cells treated with EGCG has not been reported. 110 Xiangya Road, Changsha, Hunan 410008, P.R. China In the present study, the EGCG-regulated lncRNAs and E-mail: [email protected] mRNAs were investigated in lung cancer cells by bioinfor- matics analysis. Based on the microarray analysis, a set of Abbreviations: EGCG, (-)-epigallocatechin-3-gallate; lncRNAs and mRNAs was identified whose expression levels lncRNAs, long non-coding RNAs; nt, nucleotides; MTT, were significantly modulated by EGCG in lung cancer cells, 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide; some of which are novel. However, the extent to which these GO, gene ontology; KEGG, Kyoto Encyclopedia of Genes and Genomes; FDR, false discovery rate; RT-qPCR, reverse genes are associated with lung cancer cells through EGCG transcription-quantitative polymerase chain reaction; TF, remains unclear and needs to be addressed in the near future. transcription factor Materials and methods Key words: (-)-epigallocatechin-3-gallate, long non-coding RNAs, microarray, lung cancer, bioinformatics analysis Cell culture and cell proliferation assay. The human lung cancer cell lines A549 (lung adenocarcinoma) and NCI-H460 (large cell lung carcinoma) were cultured in Dulbecco's 1510 HU et al: MICROARRAY ANALYSIS OF EGCG IN LUNG CANCER CELLS modified Eagle's medium containing 10% fetal bovine serum (www geneontology.org). The GO database provides a network (Gibco; Thermo Fisher Scientific, Inc., Waltham, MA, USA). of three structured definition terms that describe the properties Cells were plated in 96-well plates (5,000 cells/well) and were of a gene product (25). Pathway analysis of the differentially incubated overnight (37˚C, 5% CO2), followed by treatment expressed genes was performed according to the Kyoto with EGCG (20, 40, 80, 160, 320 µM) and 0 µM EGCG (control Encyclopedia of Genes and Genomes (KEGG) database group) for 24 or 48 h respectively. The MTT-based assay (www.genome.jp/kegg/). Specifically, Fisher's exact test and (Sangon Biotech Co., Ltd., Shanghai, China) was performed to χ2 tests were performed to classify the GO category and select determine viable cell numbers at 37˚C and DMSO was used to the significant pathway. The false discovery rate (FDR) was dissolve the formazan. Absorbance at 570 nm was measured calculated and the threshold of significance was defined as and the 50% maximal inhibitory concentration (IC50) of EGCG P<0.05. FunDO database (http://fundo.nubic.northwestern. for A549 and NCI-H460 cells were calculated by Graphpad edu/) was used for disease analysis. Prism version 6.0 (GraphPad Software, Inc., CA, USA). Network analysis of EGCG target genes. Search Tool for RNA extraction and quality control. Total RNA was extracted Interactions of Chemicals (STITCH; http://stitch.embl. from 1x106 A549 or NCI-H460 cells treated with different de/) is a database of interactions and correlations between concentration of EGCG and control group respectively, compounds and genes. The EGCG direct target genes were using TRIzol reagent (Invitrogen; Thermo Fisher Scientific, searched and analyzed by using STITCH according to the Inc.), purified with a mirVana miRNA Isolation kit (Ambion; reference (26). Search Tool for the Retrieval of Interacting Thermo Fisher Scientific, Inc.) according to the manufac- Genes (STRING, https://string-db.org/) is a database of turer's protocol, and was quantified using a spectrophotometer protein-protein interactions (27). The interaction of EGCG (NanoDrop ND-1000; NanoDrop, Wilmington, DE, USA). direct target gene and differential expression oncogene The RNA integrity of each sample was assessed by capillary affected by EGCG was determined in the STRING database, electrophoresis using the RNA 6000 Nano Lab-on-a-Chip kit and only the interactions with a combined score > 0.4 were and Bioanalyzer 2100 (Agilent Technologies, Santa Clara, CA, considered as significant. Finally, the network of compound USA). Only RNA extracts with RNA integrity values >6 were and gene interactions was mapped. subjected to further analysis. Analysis of lncRNA binding transcription factors. We RNA amplification, labeling, and hybridization. comple- downloaded the transcriptional factors for each of the mentary (c)DNA labeled with a fluorescent dye (Cy5 or 20 differential expression oncogenes most affected by EGCG Cy3-deoxycytidine triphosphate; CapitalBio Technology Co., treatment from the UCSC Genome Browser (http://genome. Ltd., Beijing, China) was produced by Eberwine's linear RNA ucsc.edu/cgi-bin/hgTracks?hgsid=698152999_7PLacukNAM amplification method as previously described (23), and subse- sjS2GSMIFtIxAic9jW). Then, the lncRNA-transcriptional quently RNase H enzymatic reaction (37˚C for 45 min, and factor-oncogene regulatory network was constructed by using followed by 95˚C for 5 min). The labeled cDNAs were purified the LncPro database (http://bioinfo.bjmu.edu.cn/lncpro/) as using a Capital Bioc RNA Amplification and Labeling kit previously described (28). (CapitalBio Corporation, Beijing, China) and then hybridized with specific probes (CapitalBio Technology Co., Ltd.) in a Reverse transcription‑quantitative polymerase chain reaction hybridization oven (Xinghua Analytical Instrument Factory, (RT‑qPCR). The Total RNA of 1x106 A549 cells treated with Jiangsu, China) overnight at 45˚C. 80 µM EGCG for 48 h and the control group cells was extracted as described above respectively. The PrimeScript™ RT Microarray analysis. The expression levels of lncRNAs and Reagent kit (Takara Bio, Inc., Otsu, Japan) was used for cDNA mRNAs were determined using Gene Spring software V13.0 synthesis and genomic DNA removal. qPCR was performed (Agilent Technologies). The differentially expressed lncRNAs using SYBR® Premix Ex Taq™ (Takara Bio, Inc., Otsu, Japan) and mRNAs between the EGCG-treated and control groups and in an ABI 7000 instrument (Applied Biosystems; Thermo were identified based on the threshold values of ≥2 and ≤‑2 Fisher Scientific, Inc.) for 40 cycles (95˚C for 15 sec, 60˚C for fold-change and paired Student's t-test P<0.05.