Associations Between CAMKK1 Polymorphism rs7214723 and the Prognosis of Patients With Lung Cancer

Background: The 5-year survival rate of patients with lung cancer in China is less than 20% and predicting their prognosis is challenging. We investigated the association between a common non-synonymous single nucleotide polymorphism (SNP), rs7214723, in the Ca²⁺/calmodulin-dependent protein kinase kinase 1 (CAMKK1) gene and the prognosis of patients with lung cancer.

Methods: Genomic DNA was extracted from the blood samples of 839 patients with lung cancer, recruited from Changhai Hospital (n = 536) and Taizhou Institute of Health Sciences (n = 352), and genotyped using the SNPscan technique. The association between patient prognosis and the genotypic data for CAMKK1 was analyzed using a multivariate Cox proportional hazards model adjusted for multiple potential confounders. The CRISPR/Cas9 gene-editing system was used to introduce point mutations in the CAMKK1 rs7214723 of A549 and NCI-H358 cells. Subsequently, Cell proliferation and migration ability were assessed with the Cell Counting Kit-8 and scratch assay. The Annexin V-FITC apoptosis detection kit was used to detect cell apoptosis.

Results: The CAMKK1 rs7214723 recessive CC genotype conferred significantly better overall survival (CC vs. TT + TC: adjusted hazard ratio = 0.78, 95% confidence interval [CI], 0.61-1.00, P = 0.049) than the TT + TC genotypes. Stratified analysis showed that the CAMKK1 rs7214723 CC genotype and recessive CC genotype conferred a significantly decreased risk of death in patients who were male, had a smoking history, or had stage III + IV cancer, compared with the TT and TT + TC genotypes. Relative to the TT + TC genotypes, the rs7214723 recessive CC genotype was also associated with a decreased risk of death in patients aged < 60 years (CC vs. TT + TC: adjusted hazard ratio = 0.59, 95% CI, 0.37-0.93, P = 0.024) and patients with squamous cell carcinoma (CC vs. TT + TC: adjusted hazard ratio = 0.65, 95% CI, 0.44-0.98, P = 0.038). Remarkably, CRISPR/Cas9-guided single nucleotide editing demonstrated that CAMKK1 rs7214723 T > C mutation significantly inhibits cell proliferation and migration and promotes cell apoptosis.

Conclusions: CAMKK1 SNP rs7214723 may be a significant prognostic factor for the risk of death among patients with lung cancer.

Introduction

Despite significant research efforts in the diagnosis and treatment of lung cancer in recent decades, reliable prediction of prognosis remains challenging. The 5-year survival rate of patients with lung cancer remains lower than 20% (1, 2). Although the most important prognostic determinant of lung cancer is the tumor, lymph node, and metastasis (TNM) staging system, the wide range of survival rates among patients with non-small cell lung cancer with the same stage indicates considerable heterogeneity of tumors that possibly differ in biology, which consequently confers different prognoses (3, 4). Smoking is the most significant factor associated with the risk and prognosis of lung cancer. Approximately 80%-90% of lung cancer cases are caused by smoking. However, only 15% of smokers develop lung cancer (5, 6). The other known risk factors for lung cancer include environmental factors such as air pollution and occupational exposure to carcinogens (7, 8). Moreover, the main reasons for lung cancer development are genetic and epigenetic damage caused by environmental carcinogens (7). Increasing evidence reveals that both the risk and prognosis of lung cancer are influenced by the interaction between environmental risk factors and individual genetic factors (6, 9). Therefore, establishing accurate prognostic indicators is crucial for improving survival among patients with lung cancer, whose prognosis depends on several factors.

The Ca²⁺/calmodulin-dependent protein kinases (CaMKs) family member CaMK kinase (CaMKK) has two subtypes: CaMKK1 and CaMKK2 (10, 11). Increased Ca²⁺ concentration enhances the activity of CaMKK, which phosphorylates and activates specific downstream protein kinases, including Ca²⁺/calmodulin-dependent protein kinases I (CaMK I), CaMK IV, and 5’-AMP-activated protein kinase (AMPK), which subsequently mediate multiple Ca²⁺ signaling cascades (12–14). These cascade reactions are not only involved in the development of neurons (15–17) and the regulation of mitochondrial morphology (18); they also regulate transcriptional activation through phosphorylation of cAMP-response element-binding proteins and serum response factors (19–21). It also affects transcriptional activation of glycolytic genes (22). Inhibiting the activity of CaMK/CaMKK can also promote cell apoptosis (12, 23, 24). Williams et al. (25) have demonstrated that the activation of CaMK II and CaMK IV inhibits cell cycle progression in small cell lung cancer (SCLC) cells. In addition, Hülsmann et al. (26) determined that AMPK activation contributes to the inhibition and apoptosis of growing human lung cancer cells. Therefore, we hypothesize that CAMKK1 may also play an essential role in the development of lung cancer.

Single nucleotide polymorphisms (SNPs) is a common form of genetic variation in humans. A large number of studies have explored the effects of multiple lung cancer-related gene polymorphisms on the treatment efficacy and prognosis of patients with lung cancer (27–31). The findings of Rudd’s (32) and Chen’s (33) studies on 1529 British and 320 Chinese patients, respectively, with lung cancer showed that the SNP, rs7214723, in the CAMKK1 gene can increase the risk of lung cancer. Furthermore, activation of the CaMKK/CaMK IV pathway in hepatocellular carcinoma (HCC) leads to increased secretion of high-mobility group box 1, which significantly promotes the progression of HCC and leads to poor prognosis; in contrast, CaMKK/CaMK IV inhibitors can significantly reduce the invasion and metastasis of HCC cells in vitro (34). To further investigate whether CAMKK1 SNP rs7214723 is associated with lung cancer progression, we collected peripheral blood samples from 839 patients with a confirmed lung cancer diagnosis before treatment initiation and performed genotyping and patient follow-up. The association between CAMKK1 SNP rs7214723 as well as the prognosis and survival of patients with lung cancer was evaluated. Moreover, to further demonstrate the effects of CAMKK1 SNP rs7214723 on Cell viability and apoptosis, the CRISPR/Cas9 gene editing system was used to introduce point mutations in CAMKK1 rs7214723 in A549 and NCI-H358 cells. Cell proliferation and migration were assessed with the Cell Counting Kit-8 and scratch assay, respectively. The annexin V-FITC apoptosis detection kit was used to detect cell apoptosis.

Materials and Methods

Patients

A total of 888 patients with primary lung cancer were enrolled from January to November 2009; subsequently, 49 patients were excluded because of incomplete data. Some members of our cohort were recruited from Changhai Hospital, which is affiliated to the Naval Military Medical University (Second Military Medical University) (n = 536); the others were recruited from the Taizhou Institute of Health Sciences, Fudan University (n = 352). The inclusion criterion was primary lung cancer diagnosis by histopathological examination, with no history of malignancy in other organs. No age or sex limitations were set in this study. Medical records were used to obtain patients’ clinical data, and follow-up data were collected through telephone interviews. The ethics committee of the School of Life Sciences, Fudan University approved this study. Informed consent was obtained from all study subjects.

SNP Selection and Genotyping

Prior to treatment initiation, 5-mL blood was collected from each patient. The QIAamp DNA Blood Mini Kit (Qiagen, 51106) was used to extract genomic DNA, while a 2 × 48-plex SNPscan TM kit (Genesky Biotechnologies, G0104) was used for genotyping, as previously described (35, 36). SNPscan is a proprietary multiplex SNP genotyping system that allows for the simultaneous genotyping of 48, 96, 144, or 192 SNPs per sample—all using a single tube/sample. Moreover, SNPscan discriminates alleles through a highly specific ligation reaction. We determined the genotyping quality using a detailed procedure comprising a successful call rate set at > 95%, duplicate genotype calling, having samples of internal positive controls, and performing the Hardy–Weinberg equilibrium (HWE) test. In addition, the genotyping assays were performed by personnel who were blinded to the patients’ clinical information.

Cell Lines, Media, and Cell Culture Conditions

The cancer cell lines, A549 (ATCC, CRM-CCL-185) and NCI-H358 (ATCC, CRL-5807), were employed for the in vitro assays and for the gene editing study. The cells were cultured in RPMI-1640 medium (Gibco, C22400500BT) supplemented with 10% fetal bovine serum (Hyclone, SH30084.03HI), penicillin (100 IU/mL), and streptomycin (50 μg/mL) at 37°C with 5% CO₂. The medium renewal interval was obtained from the ATCC animal cell culture guide, downloaded from the ATCC website.

Genome Editing Using CRISPR-Cas9

The sgRNA was designed using the CRISPR/Cas9 Design Tool (http://crispr.mit.edu), targeting the exon 12 of the CAMKK1 gene. The designed primers were annealed and joined to form double-stranded DNA (dsDNA). The dsDNA was inserted via BbsI (NEB R0539S) and T4 DNA Ligase into the pSpCas9 (BB)-2A-Puro (PX459) V2.0 plasmid (Addgene 62988)—named PX459.CAMKK1. The A549 and NCI-H358 cells were lipofected (Lipo8000TM, Beyotime, C0533) with the PX459.CAMKK1 plasmid and the designed ssDNA; the resulting cells were named A549 CAMKK1-W and NCI-H358 CAMKK1-W, respectively. Transfected cells were positively selected with puromycin dihydrochloride (Beyotime, ST551) 24 hours after lipofection. Positive clones were selected for polymerase chain reaction (PCR) and, subsequently, the PCR products were sequenced.

Cell Viability Assays

Cell proliferation was assessed using the Cell Counting Kit-8 (SAB biotech, CP002-1) according to the manufacturer’s instructions. Cells were seeded in 96-well plates at a density of 1 × 10³ cells/well in 100 μL cell medium and cultured at 37°C with 5% CO₂, 10 μL of CCK-8 solution was added to the cell culture after 24 hours. After 3–4 hours of incubation, cell proliferation was evaluated by spectrophotometry, using an Epoch microplate spectrophotometer (BioTek, USA), at an absorbance of 450 nm. The experiment was repeated three times. Cell migration was assessed with the scratch assay. About 5 × 10³ cells were seeded into a 6-well plate and cultured for 24 h. The cells were then scratched vertically with a pipette tip and then washed with phosphate-buffered saline (PBS) three times to remove the detached cells. Serum-free medium was added, and the cells were cultured at 37°C with 5% CO₂. After 24 hours, samples were taken, photographs were acquired, and control wells were set for each group of experiments; the assay was repeated three times.

Apoptosis Detection

The annexin V-FITC apoptosis detection kit (Beyotime, C1062S) was used to detect cell apoptosis. Cells were plated in six-well plates at a density of 5 × 10⁵ cells/well. The apoptosis rate was analyzed by flow cytometry after 24 h of incubation; control wells were set for each group of experiments, and the assay was repeated three times.

Statistical Analysis

Pearson’s chi-square test was used to check for HWE. The overall survival was calculated from the sample collection date to the date of either death from any cause or the last follow-up visit. The Kaplan–Meier method was used to estimate the median survival time, and differences between groups were tested using the log-rank test. The hazard ratio (HR) and the 95% confidence interval (CI) with age, sex and hospital adjustments were estimated using multivariate Cox regression analysis. Four SNP genetic models (allele, genotype, dominant, and recessive) were analyzed. In addition, stratified analyses were performed according to age, sex, smoking status, malignant cancer family history, TNM stage, and lung cancer histologic type. An independent sample t-test was used for comparison between two groups of cells. The statistical significance was considered at P < 0.05. All respective tests were two-sided. Statistical analyses were performed using R version 3.6.2 (Vienna, Austria).

Results

Patient Demographics and Clinical Characteristics

The follow-up period was set from the start of enrollment until November 15, 2019. After excluding 49 patients because of incomplete clinical information, the data from 839 patients were analyzed. The study sample was an ethnically homogeneous group of Han Chinese individuals. A total of 668 (79.6%) patients died, 103 (12.3%) survived for longer than 5 years, and 68 (8.1%) were lost to follow-up. A total of 610 (72.7%) patients were male, 524 (62.5%) were aged ≥ 60 years, 582 (69.4%) had a history of smoking, and 302 (36%) had a family history of malignant cancer. As for the cancer subtypes, 367 (43.7%) patients were diagnosed with adenocarcinoma, 282 (33.6%) with squamous cell carcinoma (SCC), 72 (8.6%) with SCLC, and 118 (14.1%) with other cancer types, including adenosquamous carcinoma (ASC), large cell carcinoma, carcinosarcoma (CS) and mucoepidermoid carcinoma (MEC). There were 154 (18.4%) patients who were diagnosed with stage I and stage II disease, and 625 (74.5%) patients had stage III and stage IV disease (Table 1 and Figure 1).

TABLE 1

Table 1 Distribution of characteristics in Chinese patients with lung cancer and prognosis analysis.

FIGURE 1

Figure 1 Patient demographics and clinical characteristics. We evaluated patients with lung cancer in terms of their clinical features and performed strategic analysis. For the 839 patients with primary lung cancer (excluding 49 patients because of insufficient clinical information), we evaluated their clinical features including smoking status, tumor-node-metastasis stage, and histology. Statistical analysis of demographic and clinical characteristics, single nucleotide polymorphisms (SNPs), and prognosis were performed. Stratified analysis of SNPs and prognosis were also performed. ^*Other carcinomas include ASC, LCC, CS and MST.

Association Between Patient Characteristics and Lung Cancer Outcomes

As shown in Table 1, the median survival time was significantly lower among males than females (34.27 vs. 40.17 months; P = 0.01), among patients aged ≥ 60 years than those aged < 60 years (33.2 vs. 40.87 months; P = 0.003), and among smokers than non-smokers (33.9 vs. 41.03 months; P < 0.001). Moreover, patients with advanced tumor stage had a significantly shorter median survival time than patients with early-stage tumors (29.4 vs. 113.93 months; P < 0.001) (Table 1). Association between patient characteristics and lung cancer outcomes in different hospitals showed no statistically significant difference (Supplementary Table S1).

Association Between CAMKK1 Polymorphism and Lung Cancer Prognosis

There were 341 TT, 400 TC, and 91 CC genotypes at CAMKK1 rs7214723, with a genotype detection rate of 99.17% (Table 2). The genotype frequency at rs7214723 T/C of CAMKK1 was consistent with HWE (P = 0.102), indicating that the investigated population was in a genetic balance (i.e., the population survey data were credible). Moreover, 1082 T alleles and 582 C alleles were detected at the rs7214723 locus of CAMKK1. The frequency of alleles T and C were 65.16% (864/1326) and 34.84% (462/1326), respectively, among the group of patients who died. Compared with the TT+TC genotypes, the risk of death was lower among individuals with the recessive CC genotype (adjusted HR for CC =0.78; 95% CI: 0.61-1.00, P = 0.049) (Table 2).

TABLE 2

Table 2 Association between CAMKK1 polymorphisms and prognosis in Chinese patients with lung cancer.

Association Between CAMKK1 Polymorphism and Lung Cancer Prognosis Stratified by Patient Characteristics

We observed that among male patients, compared with genotype TT, genotype CC conferred a lower risk of death than genotype TT (P = 0.026) (Table 3), as did the recessive genotype CC relative to genotypes TT+TC (P = 0.024) (Table 4). In patients less than 60 years of age, the recessive genotype CC had a lower prognostic risk of death than genotypes TT+TC (P = 0.024). Among patients who were smokers, compared with genotype TT, genotype CC conferred a lower risk of death than genotype TT (P = 0.026) (Table 3), as did the recessive genotype CC compared with genotypes TT+TC (P = 0.028) (Table 4). For patients with SCC, the rs7214723 recessive genotype CC decreased the prognostic risk of death relative to TT+TC (P = 0.038) (Table 4). In patients with stage III + IV disease, the rs7214723 genotype CC had a lower prognostic risk of death than genotype TT (P = 0.049) (Table 3). The recessive genotype CC conferred a lower prognostic risk of death than genotypes TT+TC (P = 0.020) (Table 4).

TABLE 3

Table 3 Association between CAMKK1 SNP rs7214723 in genotype models and prognosis in Chinese patients with lung cancer.

TABLE 4

Table 4 Association between CAMKK1 polymorphisms in gene model and prognosis in Chinese patients with lung cancer.

The impact of the CAMKK1 gene polymorphism rs7214723 on the survival and prognosis of patients with lung cancer are shown in Figure 2. Among all 839 patients with available survival data, there was a trend toward longer survival in patients with the rs7214723 CC genotype who were male (P = 0.0193, by log-rank test), had a smoking history (P = 0.0199), had stage III + IV disease (P = 0.0224), or were young (aged < 60 years, P = 0.0494) when compared with those with the TT+TC genotypes (Figure 2). Additionally, the dominant genotypes TC+CC showed no significant association with the prognostic risk of death relative to genotype TT (Supplementary Table S2).

FIGURE 2

Figure 2 Impacts of CAMKK1 gene polymorphism rs7214723 on the survival and prognosis of patients with lung cancer who are (A) male, (B) have smoking history, (C) with stage III + IV disease, and (D) young (aged < 60 years). The P-value for each analysis is indicated.

Effects of CAMKK1 rs7214723 on Cell Viability in A549 and NCI-H358 Cells

The CCK-8 assay was used to assess the effects of CAMKK1 rs7214723 on cell proliferation in A549 and NCI-H358 cells. As shown in Figure 3, cell growth was suppressed by CAMKK1 rs7214723 compared with the control (both P < 0.05). A scratch assay was performed to measure the effect of CAMKK1 rs7214723 on the migration of A549 and NCI-H358 cells. The result showed that liposome treatment markedly decreased the motility of A549 and NCI-H358 cells, as determined by the migration area (Figure 3).

FIGURE 3

Figure 3 CAMKK1 rs7214723 influences proliferation, migration and apoptosis in A549 and NCI-H358 cells. (A) Effects of CAMKK1 rs7214723 on proliferation in A549 and NCI-H358 cells. All data presented the means ± SD from three independent experiments. **P < 0.01 compared with the wild-type groups. (a) Effects of CAMKK1 rs7214723 on proliferation in A549 cells. (b) Effects of CAMKK1 rs7214723 on proliferation in NCI-H358 cells. (B) Effects of CAMKK1 rs7214723 on migration in A549 and NCI-H358 cells. (C) Effects of CAMKK1 rs7214723 on apoptosis in A549 and NCI-H358 cells. Both P < 0.05 indicted a significant difference, compared with the wild-type groups.

Effects of CAMKK1 rs7214723 on Apoptosis in A549 and NCI-H358 Cells

Annexin V-FITC staining assay was performed to determine whether CAMKK1 rs7214723 is associated with lung cancer progression, mediated by affecting the apoptosis process. A549 CAMKK1-w and NCI-H358 CAMKK1-w groups had remarkably increased apoptosis (both P < 0.05; Figure 3) when compared with the control. In summary, rs7214723 mutation of CAMKK1 can affect the proliferation, migration, and apoptosis of A549 and NCI-H358 cells.

Discussion

In this study, we evaluated the association between the SNP rs7214723 in CAMKK1 and prognosis using blood samples from 839 Chinese patients with lung cancer. We found that CAMKK1 rs7214723 T > C decreased the prognostic risk of death. Further stratified analysis showed that the risk of death significantly decreased among male patients, those aged < 60 years, those with a smoking history, those with SCC, and those with stage III + IV disease. Furthermore, CRISPR/Cas9-guided single nucleotide editing and cell functional experiment demonstrated that CAMKK1 rs7214723 T > C significantly inhibits A549 and NCI-H358 cell proliferation and migration and promotes cell apoptosis. To the best of our knowledge, this is the first study to investigate the association between CAMKK1 gene polymorphism rs7214723 and the prognosis of patients with lung cancer.

Previous studies have mainly focused on analyzing the association between the CAMKK1 gene SNP rs7214723 and the risk of lung cancer. Rudd (32) conducted a large-scale genome-wide association study among British patients, analyzing 1529 cases and 2707 controls, and found that the rs7214723 (E375G) polymorphism in CAMKK1 was a susceptibility locus of lung cancer in the UK population. Similarly, Chen (33) verified this conclusion in a study of 320 Chinese patients with lung cancer and 320 individuals in the control group, with the results showing that the rs7214723 T > C increased lung cancer risk among Chinese patients. However, no studies have reported on the association between the CAMKK1 SNP rs7214723 and the prognosis of lung cancer. This study stratified the aforementioned analyses by epidemiological factors such as pathological type, smoking status, and disease stage and explored the impact of the loci on the prognosis of patients with lung cancer. Our results show that the SNP rs7214723 was positively correlated with the benign prognosis of lung cancer patients in the allele and recessive models.

CAMKK1 consists of distinct N-terminal and C-terminal domains and a central kinase domain, followed by a regulatory domain consisting of overlapping autoinhibitory and calmodulin-binding domains (37). CAMKK1 is also a serine/threonine kinase. Therefore, the expression of downstream protein kinases that it activates is significantly related to the progression of human malignancies (38). Takai et al. (12) revealed that high CaMK I and CaMK II expression in endometrial cancer is related to its increased malignant potential and that increased CaMK IV expression is significantly associated with advanced endometrial cancer (39). Similarly, CaMK I and CaMKK are involved in the control of cell cycle progression of MCF-7 human breast cancer (40). In addition, by inhibiting the CaMK IV signaling pathway, the differentiation of human neuroblastoma cells induced by retinoic acid can be accelerated (23). Studies have shown that CaMK IV can be expressed in SCLC cells, and that the inactivation of this kinase can inhibit the proliferation of SCLC (25). The activation of AMPK, which is also induced by CaMKK1, can inhibit the growth and promote the apoptosis of lung cancer cells. Accordingly, it is inferred that the activation of AMPK signaling in some tumors increases the sensitivity of certain targeted therapies (26). Our study reached a similar conclusion through CRISPR/CAS9-guided single nucleotide mutation and cell functional experiment. The results showed that CAMKK1 rs7214723 T > C significantly inhibited the proliferation and migration of A549 and NCI-H358 cells, and inhibited cell apoptosis. This may be related to the excessive activation of CAMKK1 and its downstream AMPK signaling pathway caused by CAMKK1 rs7214723 T > C. These findings all prove the important role of CAMKK1 in cancer development.

The mutation SNP rs7214723 is a T to C transition resulting in the substitution of glutamate (E) to glycine (G) at the amino acid position 375. Our results indicate that rs7214723 T > C is associated with a better prognosis among patients with lung cancer, which may be due to the amino acid change in the kinase domain of CaMKK1. E375G reduces the activation ability of CaMKK1 on the downstream kinase and weakens the effect of its downstream pathway in terms of promoting the growth of lung cancer cells and inhibiting their apoptosis, which ultimately leads to a better prognosis.

In the stratified analysis, we observed that the significant association between SNP rs7214723 and the prognosis of patients with lung cancer was evident among male patients and those with history of smoking. More than 200 compounds in tobacco smoke are classified as lung carcinogens, which have been confirmed to be closely associated with increased lung cancer incidence and mortality (41, 42). Many components in tobacco smoke exert biological effects by binding the aryl hydrocarbon receptor, which is a ubiquitously expressed transcription factor (43, 44). In addition, the CaMKK/CaMK V cascade also regulates transcription by phosphorylating transcription factors such as the cAMP-response element-binding protein and the serum response factors that affect transcriptional activation (19–21). In our study, the results of the stratified analysis suggest that the CAMKK1 gene SNP rs7214723 reduced the prognostic risk of death among patients with smoking history. We hypothesize that although the SNP rs7214723 has the effect of regulating the transcriptional activity of CAMKK1, its mutation has little effect. We also hypothesized that the smoking factor increases the effect of SNP rs7214723 on the function of CAMKK1 to a certain extent, thereby delaying the progression of the lung cancer tumor, and subsequently leading to a better prognosis. Nevertheless, the specific mechanism for this remains to be clarified further. Several studies have shown that the risk of death after lung cancer diagnosis is significantly higher in men than in women (45–47). In this study, the SNP rs7214723 decreased the risk of prognostic death from lung cancer among male patients than among female patients. Smoking is a critical risk factor because men smoke more frequently than women (8, 48).

Previous reports have conducted genome-wide association studies, which can detect millions of SNPs throughout the entire genome and identify associations between SNPs and complex diseases. However, the size and complexity of the data from this method are difficult to manage, and some crucial genes could be overlooked (49). The candidate-gene associations explored in this study focused on the important locus of the CAMKK1 gene rs7214723 to further investigate the effect on lung cancer prognosis of the SNP at this locus. We conducted a large sample size study, with a total of 839 patients with lung cancer in the survival cohort. Our results confirm that the SNP rs7214723 is associated with the prognosis and survival of patients with lung cancer, reducing the probability of false positives, improving the accuracy of the analysis, and providing a reliable estimate.

This study had some limitations. First, this result should be validated with independent studies or meta-analysis, and more extensive prospective studies are needed in the future to further evaluate the association between the CAMKK1 SNP rs7214723 and lung cancer prognosis. Second, the subjects were of Han Chinese ethnicity, and inherent selection bias could not be completely excluded. This limits the generalization of our results to other populations. Finally, the number of SNPs analyzed in this study was limited. However, the correlation between the SNP rs7214723 and lung cancer prognosis was clear.

In conclusion, our results provided the first evidence that CAMKK1 rs7214723 T > C is associated with poor prognosis among patients with lung cancer, playing a more significant role depending on sex, age, smoking status, histology and TNM stage. These findings have potential clinical significance for predicting the prognosis of patients with lung cancer and in the formulation of new therapeutic strategies.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding authors.

Ethics Statement

The studies involving human participants were reviewed and approved by the ethics committee of the School of Life Sciences, Fudan University. The patients/participants provided their written informed consent to participate in this study.

Author Contributions

Conception and design: all authors. Administrative support: LJ, MX, FL, and JJW. Provision of study materials or patients: JF, YL, MX, FL, and JJW. Collection and assembly of data: HZ, ZZ, YL, XH, and CX. Cell functional experiments: BC, ZZ, and MX. Data analysis and interpretation: HZ and YL. Manuscript writing: HZ, ZZ, and JJW. Final approval of manuscript: all authors. All authors agree to be accountable for the content of the work. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by Shanghai Committee of Science and Technology (grant number 20Z11901002), National Natural Science Foundation of China (grant number 81372236, 81360359, and 82160847), and Key Project of the National Science and Technology Pillar Program (grant number 2011BAI09B00).

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.

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

The authors acknowledge all the volunteer donors involved in this study and the Changhai Hospital and Taizhou Institute of Health Sciences staff for supporting during sample collection. We thank Dr. Hao Hu, Wei Feng, and Zhunyi Gao for their assistance in recruiting the subjects. Dr. Zhenhong Zhao and Yang Bao for their laboratory assistance.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2021.757484/full#supplementary-material

Abbreviations

AMPK, 5’-AMP-activated protein kinase; ASC, adenosquamous carcinoma; CaMK, Ca²⁺/calmodulin-dependent protein kinase; CaMKK, Ca2+/calmodulin-dependent protein kinase kinase; CI, confidence interval; CS, carcinosarcoma; HCC, hepatocellular carcinoma; HR, adjusted hazard ratio; HWE, Hardy-Weinberg equilibrium; LCC, large cell carcinoma; MEC, mucoepidermoid carcinoma; MST, median survival time; SCC, squamous cell carcinoma; SCLC, small cell lung cancer; SNP, single nucleotide polymorphism.

References

1. Jemal A, Bray F, Center MM, Ferlay J, Ward E, Forman D. Global Cancer Statistics. CA Cancer J Clin (2011) 61(2):69–90. doi: 10.3322/caac.20107

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Siegel R, Naishadham D, Jemal A. Cancer Statistics, 2012. CA Cancer J Clin (2012) 62(1):10–29. doi: 10.3322/caac.20138

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Yang Y, Fei F, Song Y, Li X, Zhang Z, Fei Z, et al. Polymorphisms of Epcam Gene and Prognosis for non-Small-Cell Lung Cancer in Han Chinese. Cancer Sci (2014) 105(1):89–96. doi: 10.1111/cas.12318

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Chen HY, Yu SL, Chen CH, Chang GC, Chen CY, Yuan A, et al. A Five-Gene Signature and Clinical Outcome in Non-Small-Cell Lung Cancer. N Engl J Med (2007) 356(1):11–20. doi: 10.1056/NEJMoa060096

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Jin Y, Xu H, Zhang C, Kong Y, Hou Y, Xu Y, et al. Combined Effects of Cigarette Smoking, Gene Polymorphisms and Methylations of Tumor Suppressor Genes on non Small Cell Lung Cancer: A Hospital-Based Case-Control Study in China. BMC Cancer (2010) 10:422. doi: 10.1186/1471-2407-10-422

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Hecht SS. Cigarette Smoking and Lung Cancer: Chemical Mechanisms and Approaches to Prevention. Lancet Oncol (2002) 3(8):461–9. doi: 10.1016/s1470-2045(02)00815-x

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Sekido Y, Fong KM, Minna JD. Molecular Genetics of Lung Cancer. Annu Rev Med (2003) 54:73–87. doi: 10.1146/annurev.med.54.101601.152202

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Dalton SO, Frederiksen BL, Jacobsen E, Steding-Jessen M, Østerlind K, Schüz J, et al. Socioeconomic Position, Stage of Lung Cancer and Time Between Referral and Diagnosis in Denmark, 2001-2008. Br J Cancer (2011) 105(7):1042–8. doi: 10.1038/bjc.2011.342

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Network CGAR. Comprehensive Molecular Profiling of Lung Adenocarcinoma. Nature (2014) 511(7511):543–50. doi: 10.1038/nature13385

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Wayman GA, Lee YS, Tokumitsu H, Silva AJ, Soderling TR. Calmodulin-Kinases: Modulators of Neuronal Development and Plasticity. Neuron (2008) 59(6):914–31. doi: 10.1016/j.neuron.2008.08.021

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Soderling TR, Stull JT. Structure and Regulation of Calcium/Calmodulin-Dependent Protein Kinases. Chem Rev (2001) 101(8):2341–52. doi: 10.1021/cr0002386

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Takai N, Ueda T, Nasu K, Yamashita S, Toyofuku M, Narahara H. Targeting Calcium/Calmodulin-Dependence Kinase I and II as a Potential Anti-Proliferation Remedy for Endometrial Carcinomas. Cancer Lett (2009) 277(2):235–43. doi: 10.1016/j.canlet.2008.12.018

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Racioppi L, Means AR. Calcium/Calmodulin-Dependent Kinase IV in Immune and Inflammatory Responses: Novel Routes for an Ancient Traveller. Trends Immunol (2008) 29(12):600–7. doi: 10.1016/j.it.2008.08.005

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Matsushita M, Nairn AC. Characterization of the Mechanism of Regulation of Ca2+/Calmodulin-Dependent Protein Kinase I by Calmodulin and by Ca2+/Calmodulin-Dependent Protein Kinase Kinase. J Biol Chem (1998) 273(34):21473–81. doi: 10.1074/jbc.273.34.21473

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Wayman GA, Impey S, Marks D, Saneyoshi T, Grant WF, Derkach V, et al. Activity-Dependent Dendritic Arborization Mediated by Cam-Kinase I Activation and Enhanced CREB-Dependent Transcription of Wnt-2. Neuron (2006) 50(6):897–909. doi: 10.1016/j.neuron.2006.05.008

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Takemoto-Kimura S, Ageta-Ishihara N, Nonaka M, Adachi-Morishima A, Mano T, Okamura M, et al. Regulation of Dendritogenesis via a Lipid-Raft-Associated Ca2+/Calmodulin-Dependent Protein Kinase CLICK-III/Camkigamma. Neuron (2007) 54(5):755–70. doi: 10.1016/j.neuron.2007.05.021

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Horigane S, Ageta-Ishihara N, Kamijo S, Fujii H, Okamura M, Kinoshita M, et al. Facilitation of Axon Outgrowth via a Wnt5a-Camkk-Camkiα Pathway During Neuronal Polarization. Mol Brain (2016) 98(8):1–11. doi: 10.1186/s13041-016-0189-3

CrossRef Full Text | Google Scholar

18. Han XJ, Lu YF, Li SA, Kaitsuka T, Sato Y, Tomizawa K, et al. Cam Kinase I Alpha-Induced Phosphorylation of Drp1 Regulates Mitochondrial Morphology. J Cell Biol (2008) 182(3):573–85. doi: 10.1083/jcb.200802164

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Nakanishi A, Hatano N, Fujiwara Y, Sha’ri A, Takabatake S, Akano H, et al. AMP-Activated Protein Kinase-Mediated Feedback Phosphorylation Controls the Ca(2+)/Calmodulin (Cam) Dependence of Ca(2+)/Cam-Dependent Protein Kinase Kinase β. J Biol Chem (2017) 292(48):19804–13. doi: 10.1074/jbc.M117.805085

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Miranti CK, Ginty DD, Huang G, Chatila T, Greenberg ME. Calcium Activates Serum Response Factor-Dependent Transcription by a Ras- and Elk-1-Independent Mechanism That Involves a Ca2+/Calmodulin-Dependent Kinase. Mol Cell Biol (1995) 15(7):3672–84. doi: 10.1128/mcb.15.7.3672

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Bito H, Deisseroth K, Tsien RW. CREB Phosphorylation and Dephosphorylation: A Ca(2+)- and Stimulus Duration-Dependent Switch for Hippocampal Gene Expression. Cell (1996) 87(7):1203–14. doi: 10.1016/s0092-8674(00)81816-4

PubMed Abstract | CrossRef Full Text | Google Scholar

22. McArdle J, Schafer XL, Munger J. Inhibition of Calmodulin-Dependent Kinase Kinase Blocks Human Cytomegalovirus-Induced Glycolytic Activation and Severely Attenuates Production of Viral Progeny. J Virol (2011) 85(2):705–14. doi: 10.1128/jvi.01557-10

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Feliciano DM, Edelman AM. Repression of Ca2+/Calmodulin-Dependent Protein Kinase IV Signaling Accelerates Retinoic Acid-Induced Differentiation of Human Neuroblastoma Cells. J Biol Chem (2009) 284(39):26466–81. doi: 10.1074/jbc.M109.027680

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Kawamura K, Grabowski D, Krivacic K, Hidaka H, Ganapathi R. Cellular Events Involved in the Sensitization of Etoposide-Resistant Cells by Inhibitors of Calcium-Calmodulin-Dependent Processes. Role for Effects on Apoptosis, DNA Cleavable Complex, and Phosphorylation. Biochem Pharmacol (1996) 52(12):1903–9. doi: 10.1016/s0006-2952(96)00692-2

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Williams CL, Phelps SH, Porter RA. Expression of Ca2+/Calmodulin-Dependent Protein Kinase Types II and IV, and Reduced DNA Synthesis Due to the Ca2+/Calmodulin-Dependent Protein Kinase Inhibitor KN-62 (1-[N,O-Bis(5-Isoquinolinesulfonyl)-N-Methyl-L-Tyrosyl]-4-Phenyl Piperazine) in Small Cell Lung Carcinoma. Biochem Pharmacol (1996) 51(5):707–15. doi: 10.1016/s0006-2952(95)02393-3

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Hülsmann HJ, Rolff J, Bender C, Jarahian M, Korf U, Herwig R, et al. Activation of AMP-Activated Protein Kinase Sensitizes Lung Cancer Cells and H1299 Xenografts to Erlotinib. Lung Cancer (2014) 86(2):151–7. doi: 10.1016/j.lungcan.2014.09.001

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Kim MJ, Kang HG, Lee SY, Jeon HS, Lee WK, Park JY, et al. AKT1 Polymorphisms and Survival of Early Stage Non-Small Cell Lung Cancer. J Surg Oncol (2012) 105(2):167–74. doi: 10.1002/jso.22071

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Zhang H, Li Y, Guo S, Wang Y, Wang H, Lu D, et al. Effect of ERCC2 Rs13181 and Rs1799793 Polymorphisms and Environmental Factors on the Prognosis of Patients With Lung Cancer. Am J Transl Res (2020) 12(10):6941–53.

PubMed Abstract | Google Scholar

29. Qian J, Liu H, Gu S, Wu Q, Zhao X, Wu W, et al. Genetic Variants of the MDM2 Gene Are Predictive of Treatment-Related Toxicities and Overall Survival in Patients With Advanced NSCLC. Clin Lung Cancer (2015) 16(5):e37–53. doi: 10.1016/j.cllc.2015.02.001

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Shao L, Zuo X, Yang Y, Zhang Y, Yang N, Shen B, et al. The Inherited Variations of a P53-Responsive Enhancer in 13q12.12 Confer Lung Cancer Risk by Attenuating TNFRSF19 Expression. Genome Biol (2019) 20(1):103. doi: 10.1186/s13059-019-1696-1

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Zhang H, Fu L. The Role of ALDH2 in Tumorigenesis and Tumor Progression: Targeting ALDH2 as a Potential Cancer Treatment. Acta Pharm Sin B (2021) 11(6):1400–11. doi: 10.1016/j.apsb.2021.02.008

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Rudd MF, Webb EL, Matakidou A, Sellick GS, Williams RD, Bridle H, et al. Variants in the GH-IGF Axis Confer Susceptibility to Lung Cancer. Genome Res (2006) 16(6):693–701. doi: 10.1101/gr.5120106

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Chen D, Zhong F, Chen Y. Association of Calcium/Calmodulin-Dependent Protein Kinase Kinase1 Rs7214723 Polymorphism With Lung Cancer Risk in a Chinese Population. Biosci Rep (2017) 37(4):1–5. doi: 10.1042/bsr20170762

CrossRef Full Text | Google Scholar

34. Chen S, Dong Z, Yang P, Wang X, Jin G, Yu H, et al. Hepatitis B Virus X Protein Stimulates High Mobility Group Box 1 Secretion and Enhances Hepatocellular Carcinoma Metastasis. Cancer Lett (2017) 394:22–32. doi: 10.1016/j.canlet.2017.02.011

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Yin J, Wang L, Shi Y, Shao A, Tang W, Wang X, et al. Interleukin 17A Rs4711998 a>G Polymorphism was Associated With a Decreased Risk of Esophageal Cancer in a Chinese Population. Dis Esophagus: Off J Int Soc Dis Esophagus (2014) 27(1):87–92. doi: 10.1111/dote.12045

CrossRef Full Text | Google Scholar

36. Spinola M, Meyer P, Kammerer S, Falvella FS, Boettger MB, Hoyal CR, et al. Association of the PDCD5 Locus With Lung Cancer Risk and Prognosis in Smokers. J Clin Oncol: Off J Am Soc Clin Oncol (2006) 24(11):1672–8. doi: 10.1200/jco.2005.04.4339

CrossRef Full Text | Google Scholar

37. Racioppi L, Means AR. Calcium/Calmodulin-Dependent Protein Kinase Kinase 2: Roles in Signaling and Pathophysiology. J Biol Chem (2012) 287(38):31658–65. doi: 10.1074/jbc.R112.356485

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Tokumitsu H, Enslen H, Soderling TR. Characterization of a Ca2+/Calmodulin-Dependent Protein Kinase Cascade. Molecular Cloning and Expression of Calcium/Calmodulin-Dependent Protein Kinase Kinase. J Biol Chem (1995) 270(33):19320–4. doi: 10.1074/jbc.270.33.19320

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Shang S, Takai N, Nishida M, Miyazaki T, Nasu K, Miyakawa I. Camkiv Expression Is Associated With Clinical Stage and PCNA-Labeling Index in Endometrial Carcinoma. Int J Mol Med (2003) 11(2):181–6.

PubMed Abstract | Google Scholar

40. Rodriguez-Mora OG, LaHair MM, McCubrey JA, Franklin RA. Calcium/Calmodulin-Dependent Kinase I and Calcium/Calmodulin-Dependent Kinase Kinase Participate in the Control of Cell Cycle Progression in MCF-7 Human Breast Cancer Cells. Cancer Res (2005) 65(12):5408–16. doi: 10.1158/0008-5472.Can-05-0271

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Xue J, Yang S, Seng S. Mechanisms of Cancer Induction by Tobacco-Specific NNK and NNN. Cancers (Basel) (2014) 6(2):1138–56. doi: 10.3390/cancers6021138

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Yoshino I, Maehara Y. Impact of Smoking Status on the Biological Behavior of Lung Cancer. Surg Today (2007) 37(9):725–34. doi: 10.1007/s00595-007-3516-6

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Iqbal J, Sun L, Cao J, Yuen T, Lu P, Bab I, et al. Smoke Carcinogens Cause Bone Loss Through the Aryl Hydrocarbon Receptor and Induction of Cyp1 Enzymes. Proc Natl Acad Sci USA (2013) 110(27):11115–20. doi: 10.1073/pnas.1220919110

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Kung MH, Yukata K, O’Keefe RJ, Zuscik MJ. Aryl Hydrocarbon Receptor-Mediated Impairment of Chondrogenesis and Fracture Healing by Cigarette Smoke and Benzo(a)Pyrene. J Cell Physiol (2012) 227(3):1062–70. doi: 10.1002/jcp.22819

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Henschke CI, Yip R, Miettinen OS. Women’s Susceptibility to Tobacco Carcinogens and Survival After Diagnosis of Lung Cancer. Jama (2006) 296(2):180–4. doi: 10.1001/jama.296.2.180

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Sachs E, Sartipy U, Jackson V. Sex and Survival After Surgery for Lung Cancer: A Swedish Nationwide Cohort. Chest (2020) 159(5):2029–39. doi: 10.1016/j.chest.2020.11.010

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Salmerón D, Chirlaque MD, Isabel Izarzugaza M, Sánchez MJ, Marcos-Gragera R, Ardanaz E, et al. Lung Cancer Prognosis in Spain: The Role of Histology, Age and Sex. Respir Med (2012) 106(9):1301–8. doi: 10.1016/j.rmed.2012.06.006

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Hsu LH, Chu NM, Liu CC, Tsai SY, You DL, Ko JS, et al. Sex-Associated Differences in Non-Small Cell Lung Cancer in the New Era: Is Gender an Independent Prognostic Factor? Lung Cancer (2009) 66(2):262–7. doi: 10.1016/j.lungcan.2009.01.020

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Hardy J, Singleton A. Genomewide Association Studies and Human Disease. N Engl J Med (2009) 360(17):1759–68. doi: 10.1056/NEJMra0808700

PubMed Abstract | CrossRef Full Text | Google Scholar