Identification of m6A-modified gene signatures in lung adenocarcinoma tumorigenesis and their potential role in drug resistance
摘要
Lung cancer is one of the most commonly diagnosed cancers. N6-methyladenosine (m6A) modification has a profound impact on RNA translation, splicing, transportation, and stability.
AimsThis research aimed to identify and verify m6A-modified signatures for Lung adenocarcinoma (LUAD) tumorigenesis.
ObjectiveOur previous mRNA-seq and m6A-seq data from 26 pairs of LUAD samples and tumor-adjacent normal tissues are used.
MethodsUnivariate Cox regression analysis and the least absolute shrinkage and selection operator (LASSO) analysis were used to estimate the significance of 37 collected m6A regulators. WGCNA was constructed to identify the genes correlated with LUAD tumorigenesis. Pearson correlation analysis between mRNA-seq and m6A-seq data was used to identify the m6A-correlation genes.
ResultsLASSO-Cox analysis identified 18 m6A 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, m6A-correlation genes, with blue or brown module genes. As a result, 56 m6A-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 m6A modification. In addition, these genes may be involved in regulating drug resistance.
Conclusions56 m6A-modified gene signatures for LUAD tumorigenesis were obtained from Pearson correlation analysis between mRNA-seq and m6A-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 m6A modification-dependent manner.