Wenhan Cai
Miao Jing1
Hua Guo
Zhiqiang Xue- 1Medical School of Chinese PLA, Beijing, China
- 2Department of Thoracic Surgery, The First Medical Centre, Chinese PLA General Hospital, Beijing, China
This study focused on the epigenetic alterations of DNA methylation and miRNAs for lung adenocarcinoma (LUAD) diagnosis and treatment using bioinformatics analyses. DNA methylation data and mRNA and miRNA expression microarray data were obtained from The Cancer Genome Atlas (TCGA) database. The differentially methylated genes (DMGs), differentially expressed genes (DEGs), and differentially expressed miRNAs were analyzed by using the limma package. The DAVID database performed GO and KEGG pathway enrichment analyses. Using STRING and Cytoscape, we constructed the protein–protein interaction (PPI) network and achieved visualization. The online analysis tool CMap was used to identify potential small-molecule drugs for LUAD. In LUAD, 607 high miRNA-targeting downregulated genes and 925 low miRNA-targeting upregulated genes, as well as 284 hypermethylated low-expression genes and 315 hypomethylated high-expression genes, were obtained. They were mainly enriched in terms of pathways in cancer, neuroactive ligand–receptor interaction, cAMP signaling pathway, and cytosolic DNA-sensing pathway. In addition, 40 upregulated and 84 downregulated genes were regulated by both aberrant alternations of DNA methylation and miRNAs. Five small-molecule drugs were identified as a potential treatment for LUAD, and five hub genes (SLC2A1, PAX6, LEP, KLF4, and FGF10) were found in PPI, and two of them (SLC2A1 and KLF4) may be related to the prognosis of LUAD. In summary, our study identified a series of differentially expressed genes associated with epigenetic alterations of DNA methylation and miRNA in LUAD. Five small-molecule drugs and five hub genes may be promising drugs and targets for LUAD treatment.
1 Introduction
Lung cancer is a fatal malignancy featuring the highest incidence (11.6%) and mortality (18.4%) globally. Lung adenocarcinoma (LUAD), increasing yearly, is the most common histological type of lung cancer (Bray et al., 2018). LUAD often arises peripherally and has adenoid tissue differentiation and mucin production in histology (Chen et al., 2014). It has been documented that the pathogenesis of LUAD involves inflammation, oxidative stress, mitochondrial dysfunction, changes in lipid metabolism, and epigenetic changes (Lindskog et al., 2014; Tan et al., 2016; Zhang C. et al., 2021). Unlike other subtypes, LUAD has a high rate of gene mutations. Although targeted therapy has improved the survival rate and quality of life of nearly 60% of LUAD patients with corresponding driver gene mutations in recent years, drug resistance is still inevitable. The long-term survival rate of LUAD is still not satisfactory (Kulasingam and Diamandis, 2008; Hirsch et al., 2017). A large amount of literature studies show that environmental factors and genetic and epigenetic factors will affect the occurrence and development of lung adenocarcinoma (Fabrizio et al., 2020; Gong et al., 2020; Wang et al., 2020; Huang et al., 2021; Shi et al., 2021). Although detailed knowledge about the processes of initiation and progression of LUAD is still unknown and remains a major stumbling block on the road to LUAD treatment, robust and accurate development of biomarkers will greatly facilitate early diagnosis and treatment of biological characteristics of LUAD. Therefore, there is an urgent need to identify new therapeutic targets and some chemicals of LUAD.
Carcinogenesis is a complex process involving genetic and epigenetic changes. Abnormal genetic and epigenetic changes are the hallmarks of cancer. Epigenetic modifications can modify the gene expression without altering the DNA sequence. In cancer, deviant epigenetic regulation includes miRNA gene silencing, DNA methylation, mRNA and non-coding RNA methylation, histone methylation, and histone acetylation (Maruyama et al., 2011). The aforementioned processes are closely related and affect protein synthesis. Interference with each operation may lead to dysfunction.
miRNAs are small non-coding RNA sequences about 19–23 nucleotides in length, which are highly conserved in regulating post-translational modifications (Bartel, 2004). miRNAs can exhibit carcinogenic effects or suppressor tumors by regulating target genes. These two miRNAs are termed oncomiR and tumor suppressor (TS) miRNA, respectively
miRNAs can show carcinogenic effects or suppressor tumors by regulating target genes. These two miRNAs are termed oncomiR and tumor suppressor (TS) miRNA, respectively (Zhang et al., 2014). They have emerged as promising biomarkers for diagnostic, therapeutic, and prognostic applications due to their association with LUAD (Gu et al., 2017; Wang et al., 2019; Yuan et al., 2019). For instance, miR-196b-5p displays high expressions, whereas its target gene RSPO2 (R-Spodin 2) is expressed low in the cancer tissues and normal in para-cancer tissues, promoting proliferation and migration and invasion of LUAD (Xu and Xu, 2020).
DNA methylation is a genetic modification that does not change the DNA sequence (Santos et al., 2005). DNA methylation is associated with the subtypes and prognosis of multiple tumors, including LUAD (Fleischer et al., 2017; Long et al., 2019; Ding et al., 2020; Xu et al., 2020). Shen et al. (2019) discovered that the hypermethylation of HOXA9 and hypomethylation of TULP2, CCND1, and KRTAP8-1 could be used as biomarkers for the early detection of LUAD in the undetermined lung nodules.
As yet, although a large number of studies have demonstrated the abnormal DNA methylations or the global methylation level and miRNA level in LUAD, the comprehensive regulatory network and pathways analyses of DNA methylation levels and miRNA epigenetic alterations have not yet been conducted.
This study systematically analyzed the data on DNA methylation microarrays, miRNA expression microarrays, and mRNA expression profiling microarrays from TCGA database to identify the core genes and pathways that lead to the occurrence and development of LUAD via epigenetic regulation.
2 Materials and Methods
2.1. Microarray Data
In this study, the data on DNA methylation microarrays (including 437 LUAD and 29 adjacent normal tissue samples), miRNA expression microarrays (including 483 LUAD and 45 adjacent normal tissue samples), and mRNA expression profiling microarrays (including 497 LUAD and 54 adjacent normal tissue samples) were obtained from TCGA database (https://portal.gdc.cancer.gov/).
2.2. Data Process
The Perl script (Perl version 5.18.4) was used to process expression data to obtain mRNA and miRNA matrix. R (version 4.0.2) and Bioconductor packages were used to preprocess the raw gene expression profiles, including background correction, normalization, and logarithmic conversion. Differentially methylated probes (DMPs), differentially expressed miRNAs (DEMs), and differentially expressed genes (DEGs) were performed by using the limma package in R. DMPs were screened with P. adjust <0.05 and ∣logFC | >0.2 as the cut-off criteria. DEMs were screened with P. adjust <0.05 and ∣logFC |>2 as the cut-off criteria, and DEGs were screened with P. adjust <0.05 and ∣logFC |>1 as the threshold. Draw Venn Diagram online software (http://bioinformatics.psb.ugent.be/webtools/Venn/) was used to find overlapping genes from DMPs, DEMs, and DEGs. Aberrant methylated and expressed genes were overlapped to obtain hypermethylated low-expression genes and hypomethylated high-expression genes. Subsequently, high miRNA-targeting downregulated genes and low miRNA-targeting upregulated genes were obtained via overlapping potential targets of DEMs and DEGs. The Kaplan–Meier plotter database (https://kmplot.com/analysis/index.php?p=service&cancer=lung) was used for the survival analysis of hub genes.
2.3. Prediction of Potential Targets of miRNAs and Construction of the miRNA–mRNA Network
The targets of DEGs were predicted by microT-CDS online software of DIANA TOOLS (http://diana.imis.athena-innovation.gr/DianaTools/index.php?r=microT_CDS/index) and the miRWalk database (http://mirwalk.umm.uni-heidelberg.de/). In addition, the Cytoscape tool (v3.7.2) was used to construct the entire miRNA–mRNA regulatory network.
2.4. Functional and Pathway Enrichment Analysis
Gene ontology (GO) analyses, including the biological process (BP), cellular component (CC), and molecular function (MF), were conducted for the upregulated genes, downregulated genes, hypermethylation-low-expression genes, and hypomethylation-high-expression genes selected by DAVID (https://david.ncifcrf.gov/). Subsequently, we performed the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses for the high miRNA-targeting downregulated genes, low miRNA-targeting upregulated genes, hypermethylation-low-expression genes, and hypomethylation-high-expression genes. All analyses were performed with p < 0.05 as the screening condition.
2.5. Protein–Protein Interaction Network Construction and Module Analysis
We used the Search Tool of the Retrieval of Interacting Genes (STRING) online tool to perform PPI networks of hypermethylation-low-expression genes and hypomethylation-high-expression genes, respectively. cytoHubba in Cytoscape software was used to obtain hub genes within the PPI (top 10 nodes ranked by degree). The functional and pathway enrichment analysis of the genes in each module was performed by DAVID with p < 0.05 as the threshold.
2.6. Real-Time Quantitative PCR
PC9 and BEAS-2B cell lines were purchased from Zhejiang Meisen Cell Technology Co., Ltd. (MeisenCTCC). The PC9 cell line was cultured using 1640 + 10% FBS+1% anti-anti. The BEAS-2B cell line was cultured using BEGM, with 1% anti-anti added to the culture. Real-time Quantitative PCR was performed using Bio-Rad CFX96. After cell culture, the cells were washed three times with iced PBS. The RNA isolater Total RNA Extraction Reagent (Vazyme) was used to isolate the total RNA from cells. Then, 1 μg of total RNA and HiScript III-RT SuperMix for qPCR (Vazyme) were used for reverse transcription, according to the manufacturer’s instructions. Amplification reactions were set up in 20 μL volume containing ChamQ Universal SYBR qPCR Master Mix (Vazyme) and amplification primers according to the manufacturer’s instructions. The primer sequences used for real-time PCR are listed as follows. An amount of 5ng of cDNA was used in each amplification reaction.
The primer sequences for PCR amplification were as follows: SLC2A1, forward: 5′-TCTGGCATCAACGCTGTCTTC-3′ and reverse: 5′-CGATACCGGAGCCAATGGT-3′; PAX6, forward: 5′-TGGGCAGGTATTACGAGACTG-3′ and reverse: 5′-ACTCCCGCTTATACTGGGCTA-3′; LEP forward: 5′-TGCCTTCCAGAAACGTGATCC-3′ and reverse: 5′-CTCTGTGGAGTAGCCTGAAGC-3′; KLF4 forward: 5′-CGGACATCAACGACGTGAG-3′ and reverse: 5′-GACGCCTTCAGCACGAACT-3′; FGF forward: 5′-CAGTAGAAATCGGAGTTGTTGCC-3′ and reverse: 5′-TGAGCCATAGAGTTTCCCCTTC-3′; and β actin forward: 5′-CATGTACGTTGCTATCCAGGC-3′ and reverse: 5′-CTCCTTAATGTCACGCACGAT-3′.
2.7. Drug Exploration in CMap
The Connectivity Map (CMap) database (https://www.broadinstitute.org/) contains gene expression profiles of human cells treated with small bioactive molecules. Researchers can use CMap to identify connections among small molecules that share a physiological process, chemicals, and actions and then predict the potential drugs (Lamb et al., 2006). We used the CMap database to identify potential small-molecule drugs that reverse or induce DEGs’ modified expression in LUAD cell lines (mean range from −0.5 to 0.5 and p < 0.01).
2.8. Statistical Analyses
All the results were analyzed and processed by GraphPad Prism 8 software. The unpaired t-test was used for statistical analysis, and the data were expressed as mean ± standard deviation. p < 0.05 was considered statistically significant; ∗ meant p < 0.05, ∗∗ meant p < 0.01, ∗∗∗ meant p < 0.001, and ∗∗∗∗ meant p < 0.0001.
3 Results
3.1. Identification of Abnormal Methylated Differentially Expressed Genes in LUAD
The characteristics of mRNA and miRNA transcriptome profiling and DNA methylation profiling based on the TCGA database are shown in Supplementary Table S1. In mRNA expression profiling microarrays of TCGA, a total of 12972 DEGs were screened in cancer tissue samples from LUAD, including 10076 upregulated genes and 2896 downregulated genes. Simultaneously, 13 high-expressed miRNAs and 18 low-expressed miRNAs were identified in miRNA expression microarrays of TCGA database. Supplementary Table S2 records the characteristics of the top five differentially expressed miRNAs and their potential target DEGs. As to DNA methylation microarrays, 2405 hypermethylated genes and 2155 hypomethylated genes were found.
Finally, 607 high miRNA-targeting downregulated genes and 925 low miRNA-targeting upregulated genes were screened via overlapping target genes of DEMs and DEGs (Figures 1A,B). In addition, 284 hypermethylation-low-expression genes and 315 hypomethylation-high-expression genes by overlapping abnormal methylation and regulated genes were identified (Figures 1C,D).
FIGURE 1. Identification of target genes of differentially expressed miRNAs and mRNA, as well as aberrantly methylated differentially expressed genes between cancer and adjacent samples from LUAD patients. (A,B) Target genes of differentially expressed miRNAs and mRNAs; miRNA target genes were predicted by microT-CDS online software and the miRWalk database, respectively. (C,D) Aberrantly methylated differentially expressed genes.
Of all DMGs, 51.76% were hypermethylated, and 48.33% were hypomethylated (Figure 2A). Moreover, Figures 2B,C suggested that the DEGs and DMGs (top 20 upregulated and top 20 downregulated genes, as well as top 20 hypermethylation and top 20 hypomethylation genes) can be differentiated between LUAD and normal samples.
FIGURE 2. DMG and DEG analysis in cancer and adjacent samples from LUAD patients. (A) Volcano plot of DMGs; the percentages of hypermethylated and hypomethylated DMGs are displayed on top. (B) Heatmap of top 20 DEGs (10 upregulated genes and 10 downregulated genes). (C) Heatmap of top 20 DMGs (10 hypermethylated genes and 10 hypomethylated genes); red indicates that the expression of genes is relatively upregulated, or the level of methylation is hypermethylated, and blue indicated that the expression of genes is relatively downregulated, or the level of methylation is hypomethylated.
3.2. DEGs Associated With Altered Targeting miRNAs
3.2.1. Low-Expression miRNAs and Upregulated Genes
For low-expression miRNAs and upregulated genes, 168 GO terms were screened with the thresholds of p < 0.05, which were mainly associated with the regulation of transcription and cell adhesion (Figure 3A). The most enriched KEGG pathways were the neuroactive ligand–receptor interaction, PI3K-Akt signaling pathway, focal adhesion, protein digestion, and absorption, and ECM–receptor interaction. The KEGG enrichment chart of low-expression miRNAs and upregulated genes is shown in Figures 3C,D. It also shows the miRNA-mRNA network of the 925 upregulated genes.
FIGURE 3. Visualized regulatory network and enrichment bubble graph for miRNA-targeting DEGs. (A) Enrichment bubble graph for upregulated genes. (B) Enrichment bubble graph for downregulated genes. (C) KEGG enrichment bubble chart of low-expression miRNAs and upregulated genes and high-expression miRNAs and downregulated genes. (D) Regulatory network graph of 18 low-expression miRNAs. (E) Regulatory network graph of 13 high-expression miRNAs.
3.2.2. High-Expression miRNAs and Downregulated Genes
A total of 607 high-expression miRNAs and downregulated genes were enriched in 169 GO terms with the thresholds of p < 0.05, which were mainly associated with the regulation of transcription, signal transduction, and cell adhesion (Supplementary Table S3, Figure 3B). The most enriched KEGG pathways were pathways in cancer, HTLV-I infection, neuroactive ligand–receptor interaction, focal adhesion, and the Rap1 signaling pathway. The KEGG enrichment chart of high-expression miRNAs and downregulated genes is shown in Figure 3C. Meanwhile, we constructed the miRNA-mRNA network to reveal further significant miRNA/mRNAs regulated in LUAD progression (Figure 3E).
3.3. DEGs Associated With Altered DNA Methylation
3.3.1. Hypermethylation and Low-Expression Genes
Functional enrichment analysis of hypermethylation and low-expression genes suggested that 141 GO terms were recognized with the thresholds of p < 0.05, such as the regulation of transcription, signal transduction, and cell adhesion (Figure 4A). The most enriched KEGG pathways were the neuroactive ligand–receptor interaction, pathways in cancer, cAMP signaling pathway, signaling pathways regulating pluripotency of stem cells, and cell adhesion molecules (CAMs) (Supplementary Table S4, Figure 4C). In total, 191 nodes and 486 edges are shown in the PPI network (Supplementary Figure S1A).
FIGURE 4. Protein-protein interaction (PPI) network and enrichment bar graph for methylation-related DEGs. (A) Enrichment bubble graph for hypermethylated downregulated genes. (B) Enrichment bubble graph for hypomethylated upregulated genes. (C) KEGG enrichment bubble chart of hypermethylated downregulated genes and hypomethylated upregulated genes. (D) PPI network of hypermethylation and low-expression genes between cancer and adjacent samples from LUAD patients (only displays the 10 hub genes and the expanded subnetwork identified by the cytoHubba app in Cytoscape). (E) PPI network of hypomethylation and high-expression genes between cancer and adjacent samples from LUAD patients (only displays the 10 hub genes and the expanded subnetwork identified by the cytoHubba app in Cytoscape).
WNT3A, WNT2, SOX17, CD34, PAX6, GATA2, GNG11, BDNF, GATA6, and KLF4 were identified as hub genes by the degree rank with the cytoHubba app in Cytoscape (Figure 4D, Supplementary Table S5). In these 10 hub genes, WNAT3A is calculated with the highest degree (degree = 24). The Kaplan–Meier survival analysis showed that low-expression WNT3, SOX17, GATA2, GATA6, and KLF4 were correlated significantly with poor OS (Figure 5A).
FIGURE 5. (Continued).
3.3.2. Hypomethylation and High-Expression Genes
As for hypomethylation and high-expression genes, 39 GO terms were identified with the thresholds of p < 0.05 (Figure 4B). KEGG pathway analysis recognized enriched cytosolic DNA-sensing pathway and bile secretion (Supplementary Table S4, Figure 4C). In total, 133 nodes and 238 edges were shown in the PPI network (Supplementary Figure S1B).
GAPDH, IVL, IL17A, SPRR1B, SPRR1A, SPRR3, SPRR2G, S100A7, and SPRR2A were identified as hub genes by degree rank with the cytoHubba app (Figure 4E, Supplementary Table S5). GAPDH is calculated with the highest degree (degree = 23). The Kaplan–Meier survival analysis showed that high-expression GAPDH, SPRR1B, SPRR1A, SPRR3, and S100A7 were all significant with poor OS (Figure 5B).
3.4. DEGs Associated With Both Abnormal miRNA and DNA Methylation
We found that several DEGs were regulated by both abnormal miRNA and DNA methylation. It suggested that these DEGs might be of vital importance in the occurrence and development of LUAD. A total of 84 genes such as FAT4, KLF4, and EPB41L3 were downregulated under the regulation of both increased miRNA and hypermethylation (Figure 6A). Coincidentally, 40 genes such as SUGCT, RNF43, and UGT2B15 were upregulated under the regulation of both decreased miRNA and hypomethylation (Figure 6B). The modulatory miRNA and binding sites, as well as the DNA methylation site, cg ID, and its relation to CpG island, are summarized in Supplementary Table S6. In total, 86 genes, including 29 hypomethylation miRNA-targeting upregulated genes and 57 hypermethylation miRNA-targeting downregulated genes (11 hypomethylated miRNA-targeting upregulated genes and 26 hypermethylated miRNA-targeting downregulated genes without corresponding Affymetrix Probe Set ID on GPL96 cannot be used for CMap), were submitted to the CMap online tool to predict potential drugs in the therapy for LUAD depending on the expression alteration. By ranking the p-value in the ascending order and filtering the mean range from −0.5 to 0.5, five small-molecule chemicals were identified as latent treatment options for LUAD (Table 1). Furthermore, a PPI network for all the abnormal expressed genes, including 40 upregulated genes and 84 downregulated (63 nodes and 61 edges) genes were constructed (Figure 6C). Five hub genes were identified for further analysis, including SLC2A1 with up-regulated expression levels under both low miRNA and hypomethylation regulation, as well as four genes with downregulated expression levels under both high miRNA and hypermethylation regulation of PAX6, LEP, KLF4, and FGF10. The Kaplan–Meier survival analysis showed that high-expression SLC2A1 and low-expression KLF4 were all significant with poor OS (Figure 6D). To verify the difference in the expression of the five hub genes in TCGA database, we used qRT-PCR to evaluate the expression of the five hub genes at the transcription level and found that the expression levels of SLC2A1 mRNA (p < 0.0001, Figure 7), LEP mRNA (p < 0.05, Figure 7), and FGF4 mRNA ((p < 0.05, Figure 7) in the PC9 cell line were significantly higher than those in the BEAS-2B cell line. The expression levels of KLF4 mRNA (p < 0.0001, Figure 7) and PAX6 mRNA (p < 0.001, Figure 7) in the PC9 cell line were significantly lower than those in the BEAS-2B cell line. A CpG island prediction has been proceeded, and the results are shown in Figure 8A. The JASPER database predicted the sequence of four possible transcription factors of SLC2A1 and KLF4, as shown in Figures 8B,C. However, further clinical trials are required to verify these findings.
FIGURE 6. Details for all the overlapped genes. (A,B) Venn graph for all the overlapped genes including 84 downregulated genes and 40 upregulated genes, respectively. (C) PPI network of all the overlapped genes including 84 downregulated genes and 40 upregulated genes. (D) Kaplan–Meier analysis of patients with low (black curve) and high (red curve) expressions of SLC2A1 and KLF4.
TABLE1. Five chemicals were predicted as putative therapeutic agents for LUAD.
FIGURE 7. Comparison of the expressions of the five hub genes between the PC9 cell line and BEAS-2B cell lines.
FIGURE 8. Details of the five screened genes. (A) Epigenetic regulatory patterns of both miRNA and DNA methylation-regulating gene expression, as well as CpG island prediction analysis of five screened genes. (B) The JASPER database predicted the sequence of two possible transcription factors of SLC2A1. (C) The JASPER database predicted the sequence of two possible transcription factors of KLF4.
4 Discussion
CpG island-specific methylation in the promoter region of genes is associated with gene silencing (Morgan et al., 2018), which changes the expression of downstream hub genes and promotes abnormal cell proliferation (Kulis and Esteller, 2010). DNA methylation and miRNA expression can make a real difference in LUAD by up- or downregulating gene expressions (Herbst et al., 2018; He et al., 2021). Aberrant DNA methylation and miRNA expression can be regarded as impactful biomarkers to distinguish LUAD from normal samples (Ren et al., 2019; Sherafatian and Arjmand, 2019), which would be helpful in diagnosis, assessment of treatment, and prediction of prognosis (Ye et al., 2021). In this present study, data on DNA methylation microarrays, miRNA expression microarrays, and mRNA expression profiling microarrays (the aforementioned data are all obtained from TCGA database) were methodically analyzed, which compare the differential profiling between cancer and adjacent samples from LUAD patients. Hub genes and core pathways have been enriched to screen pivotal events in epigenetic alteration regulated by DNA methylation and miRNA.
A total of 607 high miRNA-targeted downregulated genes were identified through overlapping targets of DEMs and DEGs. The GO analysis showed that these 607 genes are primarily enriched in the cellular component in LUAD, reminding us of the potential regulation of membrane-related metabolism in LUAD. Furthermore, for molecular function, these genes were significantly enriched in protein binding, which indicated an interaction of any protein or protein complex in LUAD. As for KEGG pathway analysis, the target genes were most enriched in pathways in cancer, which suggested that these genes may participate in the tumorigenesis in LUAD. Previous research has indicated that hsa-miR-1269a had the most target genes among the 13 high-expression miRNAs, including NEGR1, ITGA8, CLDN18, JAM2, and JAM3, associated with cell adhesion. Cell adhesion molecules are a type of membrane surface glycoprotein molecules involved in regulating inflammatory response and promoting the metastasis of LUAD (Liu et al., 2021). Loss of cell adhesion is one of the characteristics of epithelial-to-mesenchymal transition (EMT), and the low expression of adhesion molecules is associated with distant metastasis in LUAD (Kim et al., 2013).
A total of 925 low miRNA targeted upregulated genes by overlapping targets of DEMs and DEGs were finally exhibited. The GO term analysis indicated that the upregulated genes were primarily enriched in integral components of membrane, positive regulation of transcription from RNA polymerase II promoter, and calcium ion binding, which indicated a regulatory role in RNA translation and transcription. The previous study has shown that the activation of Ca2+ in cells may be related to tumorigenicity and metastasis in LUAD (Li et al., 2018). KEGG analysis revealed pathways including the neuroactive ligand–receptor interaction, PI3K-Akt signaling pathway, focal adhesion, and protein digestion and absorption. GABA receptors are regulated by neuroactive steroids and are considered to control cell proliferation (Watanabe et al., 2006). A PI3K-Akt signaling pathway is a key signal medium that activates EMT-induced transcription factors (Karimi Roshan et al., 2019). According to the research, hsa-let-7c-5p upregulated 236 genes, including COL1A1, COL24A1, LAMA1, ITGA2, and other genes, and mainly enriched in focal adhesion is associated with EMT (Wu et al., 2021).
Until now, 284 hypermethylation and low expression genes were obtained via overlapping strategies of DMGs and DEGs. KEGG pathway analysis showed that hypermethylation-induced disorder of Neuroactive ligand-receptor interaction and Pathways in cancer might cause LUAD. The PPI network of hypermethylation and low-expression genes shows their functional connections; not only the top 10 hub genes among them but also five genes related to prognosis, such as WNT3, SOX17, GATA2, GATA6, and KLF4, were also selected. Interestingly, eight out of the top 10 hub genes were enriched in the biological process of positive regulation of transcription from the RNA polymerase II promoter.
As for 315 low-methylation and high-expression genes, overlapping hypomethylation and upregulation in LUAD, GO, and KEGG pathway analysis showed enrichment in the innate immune response and cytosolic DNA-sensing pathway. A series of studies indicated that the cytosolic DNA-sensing pathway was associated with antitumor immunity (Amouzegar et al., 2021; Verrier and Langevin, 2021). Therefore, hypomethylation-induced aberrance of high-expression genes may affect the antitumor immunity and promote the progression of LUAD. GAPDH, SPRR1B, SPRR1A, SPRR3, and S100A7 associated with the prognosis of LUAD are obtained through the PPI network. GAPDH is indeed the internal reference used in PCR and Western blot analysis, but GAPDH has also been shown to be dysregulated in the lung, kidney, breast, stomach, glioma, liver, colorectal, melanoma, prostate, pancreatic, and bladder cancers, and GAPDH is generally upregulated in many types of cancer. GAPDH could be utilized as a reference gene for normalizing lung cell lines, while it was de-regulated in non–small cell lung cancer specimens (Schmidt et al., 2005; Nguewa et al., 2008). The de-regulation of GAPDH in tumor tissues or cells demonstrated that the utilization of GAPDH as a reference gene/protein should be chosen very carefully (Guo et al., 2013). In normal tissues, small proline-rich proteins (SPRRs) are involved in the structural integrity of the cornified cell envelope (Patel et al., 2003), and the upregulation of SPRRs is also common under various pathophysiological conditions. Compelling evidence shows that SPRR downregulates p53 and promotes EMT (Demetris et al., 2008; Mizuguchi et al., 2012). A large number of studies have confirmed that SPRR is related to the progression of a variety of tumors (Carregaro et al., 2013). SPRR1 B activates the MAPK signaling pathway involved in LUAD proliferation and metastasis (Zhang Z. et al., 2021), but the influence of other genes of the SPRR family on the progression of LUAD is poorly understood, and further experimental elucidation is needed.
Interestingly, the abnormal expression of DEGs may be regulated by combining the epigenetic alterations of DNA methylation and miRNA. Forty genes such as SUGCT, RNF43, and UGT2B15 were upregulated due to the regulation of both decreased DNA methylation and miRNAs, while under the modulation of both increased DNA methylation and miRNAs, 84 genes including FAT4, KLF4, and EPB41L3 were downregulated. GO analysis for 40 low miRNA-targeting high-expression hypomethylation genes identified enrichment in cell adhesion, glucose homeostasis, and cellular response to interleukin-1. The three upregulated genes (CCL7, ADAMTS12, and PCK1) are involved in the cellular response to tumor necrosis factor. Moreover, the aforementioned genes are also involved in the glucagon signaling pathway, insulin resistance, bile secretion, and adipocytokine signaling pathway. For 84 high miRNA-targeting and low-expression hypermethylation genes, GO analysis screened the most significantly enriched CC, BP, and MF are integral components of membrane, positive regulation of transcription from the RNA polymerase II promoter, and sequence-specific DNA binding, respectively. Two significant results, neuroactive ligand–receptor interaction and signaling pathways regulating pluripotency of stem cells, were retrieved from the KEGG pathway analysis. Three of the five hub genes (PAX6, LEP, and KLF4) are involved in the aforementioned two pathways; the low expression of KLF4 is associated with poor prognosis. In addition, qRT-PCR was performed to verify the differential expressions of the five hub genes in LUAD. The mRNA expression level of SLC2A1 in the PC9 cell line was observed to be significantly higher than that in the BEAS-2B cell line, and the expression level of mRNA of KLF4 in the PC9 cell line was significantly lower than that in the BEAS-2B cell line. The results of qRT-PCR are consistent with those of bioinformatics analysis, meaning that SLC2A1 may be an oncogene in LUAD, while KLF4 may be a tumor suppressor gene.
Although chemotherapy, targeted therapy, and immunotherapy have brought hope to LUAD patients, drug resistance is still inevitable. It is urgent to find new therapeutic targets, explore new drugs, or reuse existing drugs; online databases can help us predict drugs. By far, the effectiveness of the CMap database has been confirmed by a large number of studies due to its practical value in drug prediction (Aramadhaka et al., 2013; Wang et al., 2016). From the CMap database, five compounds, including mecamylamine, LM-1685, 5182598, tetracycline, and aminoglutethimide, may have significant therapeutic effects on LUAD. Mecamylamine is a nicotinic acetylcholine receptor (nAChR) antagonist; research by Zhu et al. (2003) showed that mecamylamine could reverse the increase in VEGF and circulating endothelial progenitor cells (EPC) caused by secondhand smoke, thereby inhibiting tumor growth and angiogenesis. LM-1685 is a kind of selective COX-2 inhibitor, which induces cancer cell apoptosis and cell cycle arrest and inhibits tumor angiogenesis (Wu et al., 2004; Grosch et al., 2006; Liggett et al., 2014). It was reported that the selective COX-2 inhibitor might enhance the effect of conventional antitumor treatments by intensifying the sensitivity of lung cancer cells to NK cell-mediated cytotoxicity (Kim et al., 2020). It is observed that 5182598 has been reported to be an effective anti-tumor drug from the group of benzylisoquinoline alkaloids (Cordell et al., 2001).
Our research still has some shortcomings. The drug prediction results from the CMap database require a large number of rigorous clinical trials to corroborate their availability in the treatment of LUAD. In addition, the effects of both abnormal DNA methylation and miRNA expression on gene expression also need to be verified by corresponding experiments.
This study indicated that a cavalcade of abnormal methylated differentially expressed genes is related to the epigenetic changes of DNA methylation and miRNAs in LUAD. In total, 607 high miRNA-targeting downregulated genes and 925 low miRNA-targeting upregulated genes were identified by overlapping targets of DEMs and DEGs, which were enriched in the pathways in cancer and the PI3K-Akt signaling pathway, respectively. Furthermore, 284 hypermethylated downregulated genes and 315 hypomethylated upregulated genes obtained by overlapping DMGs and DEGs were associated with the neuroactive ligand–receptor interaction and cytosolic DNA-sensing pathway. Interestingly, 40 genes were upregulated under the co-regulation of hypomethylation and decreased miRNA, while 84 were downregulated under the co-regulation of hypermethylation and increased miRNA. Five small-molecule drugs were identified as potential therapeutic agents for LUAD. Finally, from these genes, SLC2A1, PAX6, LEP, KLF4, and FGF10 were identified as hub genes, especially SLC2A1 and KLF4, which were related to the prognosis of LUAD, and might be used as biomarkers for the precise diagnosis and treatment of LUAD.
Data Availability Statement
Publicly available datasets were analyzed in this study. This data can be found at: TCGA.
Author Contributions
WC was responsible for the statistical analysis and wrote the manuscript. MJ and HG contributed to the review and revision of the manuscript. JW and ZX were accountable for the design. WH and HG performed the experiments. All authors read and approved the final manuscript.
Funding
This work was supported by the National Natural Science Foundation of China (No. 62076254).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The reviewer YW declared a shared affiliation with the authors to the handling editor at time of review.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We are grateful to all the members for their generous participation.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2022.817552/full#supplementary-material
References
Amouzegar, A., Chelvanambi, M., Filderman, J. N., Storkus, W. J., and Luke, J. J. (2021). STING Agonists as Cancer Therapeutics. Cancers 13, 2695. doi:10.3390/cancers13112695
Aramadhaka, L. R., Prorock, A., Dragulev, B., Bao, Y., and Fox, J. W. (2013). Connectivity Maps for Biosimilar Drug Discovery in Venoms: The Case of Gila Monster Venom and the Anti-diabetes Drug Byetta. Toxicon 69, 160–167. doi:10.1016/j.toxicon.2013.03.018
Bray, F., Ferlay, J., Soerjomataram, I., Siegel, R. L., Torre, L. A., and Jemal, A. (2018). Global Cancer Statistics 2018: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA A Cancer J. Clin. 68, 394–424. doi:10.3322/caac.21492
Carregaro, F., Stefanini, A. C. B., Henrique, T., and Tajara, E. H. (2013). Study of Small Proline-Rich Proteins (SPRRs) in Health and Disease: a Review of the Literature. Arch. Dermatol Res. 305, 857–866. doi:10.1007/s00403-013-1415-9
Chen, Z., Fillmore, C. M., Hammerman, P. S., Kim, C. F., and Wong, K.-K. (2014). Non-small-cell Lung Cancers: a Heterogeneous Set of Diseases. Nat. Rev. Cancer 14, 535–546. doi:10.1038/nrc3775
Cordell, G. A., Quinn-Beattie, M. L., and Farnsworth, N. R. (2001). The Potential of Alkaloids in Drug Discovery. Phytother. Res. 15, 183–205. doi:10.1002/ptr.890
Demetris, A. J., Specht, S., Nozaki, I., Lunz, J. G., Stolz, D. B., Murase, N., et al. (2008). Small Proline-Rich Proteins (SPRR) Function as SH3 Domain Ligands, Increase Resistance to Injury and Are Associated with Epithelial-Mesenchymal Transition (EMT) in Cholangiocytes. J. Hepatology 48, 276–288. doi:10.1016/j.jhep.2007.09.019
Ding, W., Feng, G., Hu, Y., Chen, G., and Shi, T. (2020). Co-occurrence and Mutual Exclusivity Analysis of DNA Methylation Reveals Distinct Subtypes in Multiple Cancers. Front. Cell Dev. Biol. 8, 20. doi:10.3389/fcell.2020.00020
Fabrizio, F. P., Mazza, T., Castellana, S., Sparaneo, A., and Muscarella, L. A. (2020). Epigenetic Scanning of KEAP1 CpG Sites Uncovers New Molecular-Driven Patterns in Lung Adeno and Squamous Cell Carcinomas. Antioxidants 9, 904. doi:10.3390/antiox9090904
Fleischer, T., Tekpli, X., Tekpli, X., Mathelier, A., Wang, S., Nebdal, D., et al. (2017). DNA Methylation at Enhancers Identifies Distinct Breast Cancer Lineages. Nat. Commun. 8, 1379. doi:10.1038/s41467-017-00510-x
Gong, H., Li, Y., Yuan, Y., Li, W., Zhang, H., Zhang, Z., et al. (2020). EZH2 Inhibitors Reverse Resistance to Gefitinib in Primary EGFR Wild-type Lung Cancer Cells. BMC Cancer 20, 1189. doi:10.1186/s12885-020-07667-7
Grösch, S., Maier, T. J., Schiffmann, S., and Geisslinger, G. (2006). Cyclooxygenase-2 (COX-2)-independent Anticarcinogenic Effects of Selective COX-2 Inhibitors. J. Natl. Cancer Inst. 98, 736–747. doi:10.1093/jnci/djj206
Gu, Y., Liu, S., Zhang, X., Chen, G., Liang, H., Yu, M., et al. (2017). Oncogenic miR-19a and miR-19b Co-regulate Tumor Suppressor MTUS1 to Promote Cell Proliferation and Migration in Lung Cancer. Protein Cell 8, 455–466. doi:10.1007/s13238-017-0393-7
Guo, C., Liu, S., and Sun, M.-Z. (2013). Novel Insight into the Role of GAPDH Playing in Tumor. Clin. Transl. Oncol. 15, 167–172. doi:10.1007/s12094-012-0924-x
He, J., Fu, Y., Hu, J., Chen, J., and Lou, G. (2021). Hypomethylation-Mediated AGR2 Overexpression Facilitates Cell Proliferation, Migration, and Invasion of Lung Adenocarcinoma. Cmar Vol. 13, 5177–5185. doi:10.2147/CMAR.S304869
Herbst, R. S., Morgensztern, D., and Boshoff, C. (2018). The Biology and Management of Non-small Cell Lung Cancer. Nature 553, 446–454. doi:10.1038/nature25183
Hirsch, F. R., Scagliotti, G. V., Mulshine, J. L., Kwon, R., Curran, W. J., Wu, Y.-L., et al. (2017). Lung Cancer: Current Therapies and New Targeted Treatments. Lancet 389, 299–311. doi:10.1016/S0140-6736(16)30958-8
Huang, J., Zhang, Q., Shen, J., Chen, X., and Ma, S. (2021). Multi-omics Analysis Identifies Potential Mechanisms of AURKB in Mediating Poor Outcome of Lung Adenocarcinoma. Aging 13, 5946–5966. doi:10.18632/aging.202517
Karimi Roshan, M., Soltani, A., Soleimani, A., Rezaie Kahkhaie, K., Afshari, A. R., and Soukhtanloo, M. (2019). Role of AKT and mTOR Signaling Pathways in the Induction of Epithelial-Mesenchymal Transition (EMT) Process. Biochimie 165, 229–234. doi:10.1016/j.biochi.2019.08.003
Kim, H., Yoo, S. B., Sun, P., Jin, Y., Jheon, S., Lee, C. T., et al. (2013). Alteration of the E-Cadherin/β-Catenin Complex Is an Independent Poor Prognostic Factor in Lung Adenocarcinoma. Korean J. Pathol. 47, 44–51. doi:10.4132/KoreanJPathol.2013.47.1.44
Kim, J., Noh, M., Hur, D., Kim, B., Kim, Y., and Lee, H.-K. (2020). Celecoxib Upregulates ULBP-1 E-xpression in L-ung C-ancer C-ells via the JNK/PI3K S-ignaling P-athway and I-ncreases S-usceptibility to N-atural K-iller C-ell C-ytotoxicity. Oncol. Lett. 20, 1. doi:10.3892/ol.2020.12142
Kulasingam, V., and Diamandis, E. P. (2008). Strategies for Discovering Novel Cancer Biomarkers through Utilization of Emerging Technologies. Nat. Rev. Clin. Oncol. 5, 588–599. doi:10.1038/ncponc1187
Kulis, M., and Esteller, M. (2010). DNA Methylation and Cancer. Adv. Genet. 70, 27–56. doi:10.1016/B978-0-12-380866-0.60002-2
Lamb, J., Crawford, E. D., Peck, D., Modell, J. W., Blat, I. C., Wrobel, M. J., et al. (2006). The Connectivity Map: Using Gene-Expression Signatures to Connect Small Molecules, Genes, and Disease. Science 313, 1929–1935. doi:10.1126/science.1132939
Li, Y., Yu, W.-K., Chen, L., Chan, Y.-s., Liu, D., Fong, C.-C., et al. (2018). Electrotaxis of Tumor-Initiating Cells of H1975 Lung Adenocarcinoma Cells Is Associated with Both Activation of Stretch-Activated Cation Channels (SACCs) and Internal Calcium Release. Bioelectrochemistry 124, 80–92. doi:10.1016/j.bioelechem.2018.03.013
Liggett, J. L., Zhang, X., Eling, T. E., and Baek, S. J. (2014). Anti-tumor Activity of Non-steroidal Anti-inflammatory Drugs: Cyclooxygenase-independent Targets. Cancer Lett. 346, 217–224. doi:10.1016/j.canlet.2014.01.021
Lindskog, C., Fagerberg, L., Hallström, B., Edlund, K., Hellwig, B., Rahnenführer, J., et al. (2014). The Lung‐specific Proteome Defined by Integration of Transcriptomics and Antibody‐based Profiling. FASEB J. 28, 5184–5196. doi:10.1096/fj.14-254862
Liu, T.-T., Li, R., Huo, C., Li, J.-P., Yao, J., Ji, X.-l., et al. (2021). Identification of CDK2-Related Immune Forecast Model and ceRNA in Lung Adenocarcinoma, a Pan-Cancer Analysis. Front. Cell Dev. Biol. 9, 682002. doi:10.3389/fcell.2021.682002
Long, J., Chen, P., Lin, J., Bai, Y., Yang, X., Bian, J., et al. (2019). DNA Methylation-Driven Genes for Constructing Diagnostic, Prognostic, and Recurrence Models for Hepatocellular Carcinoma. Theranostics 9, 7251–7267. doi:10.7150/thno.31155
Maruyama, R., Choudhury, S., Kowalczyk, A., Bessarabova, M., Beresford-Smith, B., Conway, T., et al. (2011). Epigenetic Regulation of Cell Type-specific Expression Patterns in the Human Mammary Epithelium. PLoS Genet. 7, e1001369. doi:10.1371/journal.pgen.1001369
Mizuguchi, Y., Specht, S., Lunz, J. G., Isse, K., Corbitt, N., Takizawa, T., et al. (2012). SPRR2A Enhances P53 Deacetylation through HDAC1 and Down Regulates P21 Promoter Activity. BMC Mol. Biol. 13, 20. doi:10.1186/1471-2199-13-20
Morgan, A. E., Davies, T. J., and Mc Auley, M. T. (2018). The Role of DNA Methylation in Ageing and Cancer. Proc. Nutr. Soc. 77, 412–422. doi:10.1017/S0029665118000150
Nguewa, P. A., Agorreta, J., Blanco, D., Lozano, M., Gomez-Roman, J., Sanchez, B. A., et al. (2008). Identification of Importin 8 (IPO8) as the Most Accurate Reference Gene for the Clinicopathological Analysis of Lung Specimens. BMC Mol. Biol. 9, 103. doi:10.1186/1471-2199-9-103
Patel, S., Kartasova, T., and Segre, J. A. (2003). Mouse Sprr Locus: a Tandem Array of Coordinately Regulated Genes. Mamm. Genome 14, 140–148. doi:10.1007/s00335-002-2205-4
Ren, Z. P., Hou, X. B., Tian, X. D., Guo, J. T., Zhang, L. B., Xue, Z. Q., et al. (2019). Identification of Nine Micro RNA S as Potential Biomarkers for Lung Adenocarcinoma. FEBS Open Bio 9, 315–327. doi:10.1002/2211-5463.12572
Santos, K. F., Mazzola, T. N., and Carvalho, H. F. (2005). The Prima Donna of Epigenetics: the Regulation of Gene Expression by DNA Methylation. Braz J. Med. Biol. Res. 38, 1531–1541. doi:10.1590/s0100-879x2005001000010
Schmidt, B., Engel, E., Carstensen, T., Weickmann, S., John, M., Witt, C., et al. (2005). Quantification of Free RNA in Serum and Bronchial Lavage: a New Diagnostic Tool in Lung Cancer Detection? Lung Cancer 48, 145–147. doi:10.1016/j.lungcan.2004.09.013
Shen, N., Du, J., Zhou, H., Chen, N., Pan, Y., Hoheisel, J. D., et al. (2019). A Diagnostic Panel of DNA Methylation Biomarkers for Lung Adenocarcinoma. Front. Oncol. 9, 1281. doi:10.3389/fonc.2019.01281
Sherafatian, M., and Arjmand, F. (2019). Decision Tree-based C-lassifiers for L-ung C-ancer D-iagnosis and S-ubtyping U-sing TCGA miRNA E-xpression D-ata. Oncol. Lett. 18, 2125–2131. doi:10.3892/ol.2019.10462
Shi, R., Bao, X., Unger, K., Sun, J., Lu, S., Manapov, F., et al. (2021). Identification and Validation of Hypoxia-Derived Gene Signatures to Predict Clinical Outcomes and Therapeutic Responses in Stage I Lung Adenocarcinoma Patients. Theranostics 11, 5061–5076. doi:10.7150/thno.56202
Tan, Q., Cui, J., Huang, J., Ding, Z., Lin, H., Niu, X., et al. (2016). Genomic Alteration during Metastasis of Lung Adenocarcinoma. Cell Physiol. Biochem. 38, 469–486. doi:10.1159/000438644
Verrier, E. R., and Langevin, C. (2021). Cyclic Guanosine Monophosphate-Adenosine Monophosphate Synthase (cGAS), a Multifaceted Platform of Intracellular DNA Sensing. Front. Immunol. 12, 637399. doi:10.3389/fimmu.2021.637399
Wang, J., He, L., Tang, Y., Li, D., Yang, Y., and Zeng, Z. (2020). Development and Validation of a Nomogram with an Epigenetic Signature for Predicting Survival in Patients with Lung Adenocarcinoma. aging 12, 23200–23216. doi:10.18632/aging.104090
Wang, J., Li, M., Wang, Y., and Liu, X. (2016). Integrating Subpathway Analysis to Identify Candidate Agents for Hepatocellular Carcinoma. Ott 9, 1221–1230. doi:10.2147/OTT.S97211
Wang, S. S., Fang, Y. Y., Huang, J. C., Liang, Y. Y., Guo, Y. N., Pan, L. J., et al. (2019). Clinical Value of microRNA-198-5p D-ownregulation in L-ung A-denocarcinoma and its P-otential P-athways. Oncol. Lett. 18, 2939–2954. doi:10.3892/ol.2019.10610
Watanabe, M., Maemura, K., Oki, K., Shiraishi, N., Shibayama, Y., and Katsu, K. (2006). Gamma-aminobutyric Acid (GABA) and Cell Proliferation: Focus on Cancer Cells. Histol. Histopathol. 21, 1135–1141. doi:10.14670/HH-21.1135
Wu, T., Leng, J., Han, C., and Demetris, A. J. (2004). The Cyclooxygenase-2 Inhibitor Celecoxib Blocks Phosphorylation of Akt and Induces Apoptosis in Human Cholangiocarcinoma Cells. Mol. Cancer Ther. 3, 299–307.
Wu, Y., Liu, L., Shen, X., Liu, W., and Ma, R. (2021). Plakophilin-2 Promotes Lung Adenocarcinoma Development via Enhancing Focal Adhesion and Epithelial-Mesenchymal Transition. Cmar Vol. 13, 559–570. doi:10.2147/CMAR.S281663
Xu, F., He, L., Zhan, X., Chen, J., Xu, H., Huang, X., et al. (2020). DNA Methylation-Based Lung Adenocarcinoma Subtypes Can Predict Prognosis, Recurrence, and Immunotherapeutic Implications. Aging 12, 25275–25293. doi:10.18632/aging.104129
Xu, Q., and Xu, Z. (2020). miR-196b-5p Promotes Proliferation, Migration and Invasion of Lung Adenocarcinoma Cells via Targeting RSPO2. Cmar Vol. 12, 13393–13402. doi:10.2147/CMAR.S274171
Ye, G.-C., Liu, Y.-F., Huang, L., Zhang, C.-Y., Sheng, Y.-L., Wu, B., et al. (2021). Key microRNAs and Hub Genes Associated with Poor Prognosis in Lung Adenocarcinoma. Aging 13, 3742–3762. doi:10.18632/aging.202337
Yuan, L., Bing, Z., Yan, P., Li, R., Wang, C., Sun, X., et al. (2019). Integrative Data Mining and Meta-Analysis to Investigate the Prognostic Role of microRNA-200 Family in Various Human Malignant Neoplasms: A Consideration on Heterogeneity. Gene 716, 144025. doi:10.1016/j.gene.2019.144025
Zhang, C., Chen, X., Chen, Y., Cao, M., Tang, J., Zhong, B., et al. (2021a). The PITX Gene Family as Potential Biomarkers and Therapeutic Targets in Lung Adenocarcinoma. Med. Baltim. 100, e23936. doi:10.1097/MD.0000000000023936
Zhang, Y., Yang, Q., and Wang, S. (2014). MicroRNAs: a New Key in Lung Cancer. Cancer Chemother. Pharmacol. 74, 1105–1111. doi:10.1007/s00280-014-2559-9
Zhang, Z., Shi, R., Xu, S., Li, Y., Zhang, H., Liu, M., et al. (2021b). Identification of Small Proline‐rich Protein 1B ( SPRR1B ) as a Prognostically Predictive Biomarker for Lung Adenocarcinoma by Integrative Bioinformatic Analysis. Thorac. Cancer 12, 796–806. doi:10.1111/1759-7714.13836
Keywords: DNA methylation, miRNA, epigenetics, lung adenocarcinoma, mRNA
Citation: Cai W, Jing M, Wen J, Guo H and Xue Z (2022) Epigenetic Alterations of DNA Methylation and miRNA Contribution to Lung Adenocarcinoma. Front. Genet. 13:817552. doi: 10.3389/fgene.2022.817552
Received: 18 November 2021; Accepted: 26 April 2022;
Published: 31 May 2022.
Edited by:
Ritu Kulshreshtha, Indian Institute of Technology Delhi, IndiaReviewed by:
Neeti Sharma, Ajeenkya D Y Patil University, IndiaYuanyong Wang, People ‘s Liberation Army General Hospital, China
Copyright © 2022 Cai, Jing, Wen, Guo and Xue. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Zhiqiang Xue, xuezhiqiang301@126.com