Abstract
BACKGROUND: Lung cancer is one of the most commonly diagnosed cancers. N(6)-methyladenosine (m(6)A) modification has a profound impact on RNA translation, splicing, transportation, and stability. AIMS: This research aimed to identify and verify m(6)A-modified signatures for Lung adenocarcinoma (LUAD) tumorigenesis. OBJECTIVE: Our previous mRNA-seq and m(6)A-seq data from 26 pairs of LUAD samples and tumor-adjacent normal tissues are used. METHODS: Univariate Cox regression analysis and the least absolute shrinkage and selection operator (LASSO) analysis were used to estimate the significance of 37 collected m(6)A regulators. WGCNA was constructed to identify the genes correlated with LUAD tumorigenesis. Pearson correlation analysis between mRNA-seq and m(6)A-seq data was used to identify the m(6)A-correlation genes. RESULTS: LASSO-Cox analysis identified 18 m(6)A significant regulators. The top 3 regulators, including METTL16, FTO, and SRSF10, and their downstream genes which were reported in the literature were analysed to confirm their role in LUAD tumorigenesis. Blue and brown coexpression modules were chosen as key modules for LUAD tumorigenesis. At last, we intersected Lasso-downstream genes, m(6)A-correlation genes, with blue or brown module genes. As a result, 56 m(6)A-modified gene signatures were obtained. Among them, AKAP9, PLXNB2, BRPF3, HPS4, EXOC7, and KLF6 have an inconsistent expression in protein and mRNA levels, probably due to m(6)A modification. In addition, these genes may be involved in regulating drug resistance. CONCLUSIONS: 56 m(6)A-modified gene signatures for LUAD tumorigenesis were obtained from Pearson correlation analysis between mRNA-seq and m(6)A-seq data, along with LASSO and WGCNA analysis. Among them, AKAP9, PLXNB2, BRPF3, HPS4, EXOC7 and KLF6 play a crucial role in LUAD tumorigenesis in an m(6)A modification-dependent manner.