Skip to main content
Download PDF

MTF1 genetic variants are associated with lung cancer risk in the Chinese Han population

BMC Cancer volume 24, Article number: 778 (2024) Cite this article

Abstract

Background

Metal-regulatory transcription factor 1 (MTF1), a conserved metal-binding transcription factor in eukaryotes, regulates the proliferation of cancer cells by activating downstream target genes and then participates in the formation and progression of tumors, including lung cancer (LC). The expression level of MTF1 is down-regulated in LC, and high expression of MTF1 is associated with a good prognosis of LC. However, the association between MTF1 polymorphism and LC risk has not been explored.

Methods

The genotyping of MTF1 Single nucleotide polymorphisms (SNPs) including rs473279, rs28411034, rs28411352, and rs3748682 was identified by the Agena MassARRAY system among 670 healthy controls and 670 patients with LC. The odds ratio (OR) and 95% confidence intervals (CI) were calculated by logistics regression to assess the association of these SNPs with LC risk.

Results

MTF1 rs28411034 (OR 1.22, 95% CI 1.03–1.45, p = 0.024) and rs3748682 (OR 1.24, 95% CI 1.04–1.47, p = 0.014) were associated with higher LC susceptibility overall. Moreover, the effect of rs28411034 and rs3748682 on LC susceptibility was observed in males, subjects with body mass index (BMI) ≥ 24 kg/m2, smokers, drinkers, and patients with lung squamous carcinoma (OR and 95% CI > 1, p < 0.05). Besides, rs28411352 (OR 0.73, 95% CI 0.55–0.97, p = 0.028,) showed protective effect for reduced LC risk in drinkers.

Conclusions

We were first who reported that rs28411034 and rs3748682 tended to be relevant to increased LC susceptibility among the Chinese Han population. These results of this study could help to recognize the pathogenic mechanisms of the MTF1 gene in LC progress.

Peer Review reports

Introduction

Globally, lung cancer (LC) is the second most frequently diagnosed cancer with an estimated 2.2 million (11.4%) new cancer cases and the leading cause of cancer death with an estimated 1.8 million (18.0%) deaths [1]. LC is the leading cause of cancer-related death in many developed countries [2]. In China, LC is also a type of primary cancer with high morbidity and mortality [3]. According to the LC statistics of China National Lung Cancer Research Center in 2015, the occurrence rate of LC in males is more than twice that of females (China's age-standardized mortality rate: 40.11 per 100,000 men and 16.54 per 100,000 women) [4]. Smokers, people exposed to smoke, oil fields, toxic occupational sites, people with family history, radiation exposure, and people with chronic lung diseases are all high-risk groups for LC screening [2]. In addition, genetic variants have been reported to explain about 12% to 21% of the heritability of LC [5, 6]. Previously, some association studies have identified multiple risk loci for LC development [7,8,9]. However, the genetic factors related to LC susceptibility have not been entirely recognized and there are still a large number of potential LC risk loci that need further study, especially in the Chinese population.

Metal-regulatory transcription factor 1 (MTF1) is a conserved metal-binding transcription factor in eukaryotes that binds to conserved DNA sequence motifs. MTF1 regulates the proliferation of normal and cancer cells by activating downstream target genes and then participates in the formation and progression of tumors [10, 11]. MTF1 has tumorigenic effects and promotes epithelial-mesenchymal transition (EMT) in ovarian tumor metastasis [12]. MTF1 may be used as a diagnostic indicator of gastric cancer and is linked to a good prognosis [13]. MTF1 is down-regulated in most cancers, including lung adenocarcinoma and lung squamous cell carcinoma, and higher expression of MTF1 may predict a better prognosis for LC patients [14, 15]. Arsenic exposure is a worldwide health concern associated with an increased risk of lung cancer but arsenic trioxide (AsIII) is also an effective chemotherapeutic agent, Inactivation of ZnT1 or its transcriptional activator MTF1 resulted in considerable AsIII resistance [16]. A recent report showing that a novel Zn2+ chelator, LOR-253, can inhibit lung xenograft growth, proliferation and angiogenesis in association with changes in MTF-1 protein levels [17]. Gene expression profiles of A549 cultures treated with one of these water-soluble zinc ionophores, PCI-5002, reveal the activation of stress response pathways under the control of MTF-1 and heat shock transcription factors [18]. However, the association between MTF1 polymorphism and LC risk has not been studied.

Here, we aim to study the contribution of MTF1 gene variation to LC risk in the Chinese Han population, to investigate the interaction between MTF1 gene polymorphism and demographic characteristics on LC risk, and to explore the relationship between MTF1 polymorphisms and clinical features of LC patients.

Methods

Study population

A cohort of 1340 individuals from the First Affiliated Hospital of Xi’an Jiaotong University included 670 healthy controls and 670 LC patients with LC, all of whom were genetically unrelated individuals from the Han ethnic group. All patients were confirmed as primary LC by medical and histopathological examination according to the World Health Organization (WHO) classifications. Exclusion criteria included receiving radiotherapy or chemotherapy before the study, and a history of malignancy, chronic pulmonary diseases, and other serious organ diseases. All controls were randomly employed from the same hospital physical examination center. The control group met the inclusion criteria: 1) no personal or family history of malignant tumors, 2) no respiratory disease, infectious disease, or immune disease. These healthy controls were matched with LC patients on age, gender, body mass index (BMI), cigarette, and alcohol consumption status to minimize the effect of these confounders. Demographical and clinical data were collected through standard questionnaires and medical records. This study conformed to the Helsinki Declaration and was permitted by the Ethical Committee of the First People's Hospital of Yunnan Province (No. 2016LH036). All individuals signed informed consent forms.

Genotyping

MTF1 single nucleotide polymorphisms (SNPs) including rs473279, rs28411034, rs28411352, and rs3748682 were randomly selected for genotyping based on 1) the physical position and Variant table of the MTF1 gene on the chromosome 1:37,809,574–37859592 through the e!GRCh38.p14 (http://asia.ensembl.org/Homo_sapiens/Info/Index) database, 2) minor allele frequency (MAF) > 0.01 from the dbSNP database for genotyping, 3) principles of MassARRAY primer design, 4) HWE > 0.05, and the detection rate for genotyping > 99.5% in our study, and 5) the genotype of these SNPs associated with the MTF1 expression according to in silico analysis, as Supplementary Fig. 1. The GTEx Portal database (https://gtexportal.org/home/)[19], HaploReg v4.1 (https://pubs.broadinstitute.org/mammals/haploreg/haploreg.php)[20], and RegulomeDB (https://regulome.stanford.edu/regulome-search/)[21] bioinformatics tools were applied to predict the potential functions of these polymorphisms.

Peripheral blood samples (5 mL) of all participants were taken into EDTA vacutainer tubes, and genomic DNA was isolated by GoldMag DNA Extraction Kit (GoldMag Co. Ltd, Xi′an, China), and stored at –20˚C. The MassARRAY platform is based on MALDI-TOF (matrix-assisted laser desorption/ionization-time of flight) mass spectrometry in a high-throughput and cost-effective manner [22]. The genotyping steps for MassARRAY iPLEX were based on manufacturer’s protocol, as following: 1) DNA templates containing SNPs were amplified by PCR, and PCR products were treated with alkaline phosphatase to neutralize unincorporated nucleotides. 2) A single-base extension reaction was then performed to extend the PCR fragments by one base into the SNPs. 3) The resin purification reaction was performed, and the purified resin extension products were transferred to SpectroCHIP Assay using the purpose-built dispenser Agena Bioscience Nanodispenser RS1000. 4) Due to the different bases of polymorphism sites, the different terminal bases of extension products will lead to the difference of molecular weight after extension. The size of the product molecular weight was detected using MALDI-TOFmass spectrometry analysis. The Agena Biosicence Assay Design Suite V2.0 software (http://agenacx.com/online-tools) was used to design the extended primer. The MassARRAY Nanodispenser (Agena Bioscience, San Diego, CA, USA) and the MassARRAY iPLEX platform (Agena Bioscience, San Diego, CA, USA) were used to genotype, and then Agena Bioscience TYPER software (version 4.0) was used to analyze the data. The primers were detailed in Supplementary Table 1. Positive (Human Genome standard sample) and negative (H2O) controls were designed for each batch of genotyping experiments in the same 384-well plate. In addition, randomly selected 5% of samples were re-genotyping, with 100% consistency.

Data analysis

The Kolmogorov–Smirnov test is a non-parametric test that assesses whether a sample comes from a population with a specific distribution, such as the normal distribution. Continuous variables was evaluated for normality using the Kolmogorov–Smirnov test. The Mann–Whitney U test, also known as the Wilcoxon rank-sum test, is a non-parametric test used to determine whether there is a significant difference between the distributions of two independent samples. It is an alternative to the independent samples t-test when the data does not meet the assumptions of normality. Continuous variables with non-normal distribution as median with interquartile range (IQR) were compared using Mann–Whitney U test. Categorical variables were compared between the study groups using χ2. The genotype frequency of these variants for Hardy–Weinberg equilibrium (HWE) was detected by the χ2 test. To evaluate the relationship of these SNPs to LC risk, logistics regression adjusted for variables such as age, sex, BMI, smoking, and drinking was used for determining the odds ratio (OR) and 95% confidence intervals (CI). False-positive report probability (FPRP) analysis was used to evaluate the noteworthy associations of the significant findings to reduce the numbers of false-positive findings [23]. FPRP threshold is 0.2, and the prior probability is 0.1, which is used to evaluate the significant association of significant findings. If the FPRP value is less than the preset threshold (0.2) at the prior probability is 0.1, it means that the false positive rate of the positive result is lower than the expected value, and the positive result is noteworthy. Furthermore, a multifactorial dimensionality reduction (MDR) analysis (https://sourceforge.net/projects/mdr) [24, 25] for the influence of SNP-SNP interaction on LC risk was performed by MDR 3.0.2. Using this method, multilocus genotypes are classified into high-risk and low-risk groups, effectively reducing the genotype predictors from n dimensions to one dimension. The new, one-dimensional multilocus genotype variable is evaluated for its ability to classify and predict disease status through cross-validation (CV). The pairwise linkage disequilibrium (LD) were measured by the Lewontin’s coefficient D’ using the Haploview v4.2 software and the correlation of MTF1 haplotypes with LC risk was evaluated by logistic regression model. Data were studied by SPSS version 20.0 (SPSS Inc., Chicago, IL, USA), and p < 0.05 was significant. The statistical significance of the Bonferroni-corrected p values would be set at p < 0.05/(4 × 4 × 5).

Results

Study population characteristics

Table 1 displayed the frequency distributions of age, sex, smoking, alcohol consumption, and BMI between 670 patients with LC (61 years, male/female: 472/194) and their 670 matched non-cancer healthy controls (61 years, male/female: 464/206). The table also listed the clinical data of LC patients, including pathological type, stage, and lymph node metastasis. No significant differences in age (p = 0.371), sex (p = 0.547), smoking (p = 0.291) drinking (p = 0.229) status, and BMI (p = 0.993) between the two groups were observed. Adenocarcinoma accounted for 46.4% (311/670), followed by squamous cell carcinoma (31.6%, 212/670).

Table 1 Characteristics of patients with lung cancer and healthy controls
Full size table

Association between MTF1 SNPs and susceptibility to LC

The information on MTF1 SNPs (rs473279, rs28411034, rs28411352, and rs3748682) was presented in Table 2, including Chr: position, role, allele, MAFs (> 0.05) of cases and controls, p for HWE (p > 0.05), and the potential function of these SNPs. Moreover, rs28411034 (OR 1.20, 95% CI 1.02–1.41, p = 0.032) and rs3748682 (OR 1.22, 95% CI 1.03–1.43, p = 0.020) tended to have increased LC susceptibility in the allele model. Here, none of the SNPs was significant after the Bonferroni correction.

Table 2 Basic information of candidate SNPs of MTF1 in the study
Full size table

By HaploReg annotation, we found that the possible function of these SNPs might be related to SiPhy cons, promoter/ enhancer histone marks, DNAse, motifs changed, proteins bound, GRASP quantitative trait locus (QTL) hits, and selected expression quantitative trait locus (eQTL) hits. The results of RegulomeDB displayed that these selected SNPs were related to eQTL/caQTL, transcription factor binding, and/or chromatin accessibility peak. Based on the GTEx Portal database, the genotypes of rs473279 (p = 4.55e-6), rs28411034 (p = 5.96e-6), rs28411352 (p = 6.88e-37), and rs3748682 (p = 7.05e-19) in MTF1 were associated with the mRNA expression in lung tissues (Supplementary Fig. 1). These results might suggest that these polymorphisms may be involved in LC carcinogenic by affecting the expression or function of MTF1, which provides a theoretical basis for subsequent mechanistic studies.

The genotype frequency distribution of these SNPs between two groups was shown in Table 3. Genotypes distribution analysis showed that rs28411034 A-allele carriers (GA + AA genotypes) were more frequent in the case group than in controls (55.8% vs. 50.3%). Genetic model results indicated that rs28411034 was linked with increased LC susceptibility under the dominant (OR 1.26, 95% CI 1.01–1.56, p = 0.038) and log-additive (OR 1.22, 95% CI 1.03–1.45, p = 0.024) models. Moreover, rs3748682 C-allele carriers (TC + CC genotypes) were more frequent in cases than in controls (56.0% vs. 49.6%). MTF1 rs3748682 was also related to the higher LC risk under the codominant (OR 1.28, 95% CI 1.02–1.60, p = 0.046), dominant (OR 1.30, 95% CI 1.05–1.62, p = 0.016), and log-additive (OR 1.24, 95% CI 1.04–1.47, p = 0.014) models. Here, none of the SNPs was significant after the Bonferroni correction.

Table 3 The association of MTF1 genetic polymorphisms and lung cancer risk
Full size table

Association of MTF1 polymorphisms with the LC risk in subgroup analysis by demographic characteristics

Stratified analyses were performed based on demographic features (sex, age, BMI, cigarette, and alcohol consumption) to investigate the combined effects of MTF1 polymorphisms and these factors on LC risk, as shown in Fig. 1 and Supplementary Table 2.

Fig. 1
figure 1

Forest map for the association of MTF1 variants with the LC risk in subgroup analysis by demographic characteristics. The red dots denote the SNP/model associated with an increased risk of lung cancer, green dots represent the SNP/model related to a decreased risk, and black dots indicate the SNP/model with no significant correlation to lung cancer risk. Furthermore, the dots signify the OR values, while the bars represent the 95%CI. LC: lung cancer; SNP: single nucleotide polymorphism; OR: odds ratio; 95% CI: 95% confidence interval; BMI: mass body index

Full size image

Gender: MTF1 rs28411034 (dominant and log-additive) and rs3748682 (codominant, dominant, and log-additive) were correlated with LC risk in males.

Age: No significant association of these SNPs with LC risk was observed in the subjects aged ≤ 50 years, subjects aged 61–70 years, and subjects aged > 70 years. We found that rs473279 (codominant and recessive) was associated with the higher LC risk in the subjects aged 51–60 years.

BMI: Rs28411034 (codominant, dominant, recessive, and log-additive) and rs3748682 (codominant, dominant, and log-additive) might be conferred to LC risk in the subjects with BMI ≥ 24 kg/m2.

Smoking: The risk-increasing association of rs28411034 (dominant and log-additive) and rs3748682 (codominant, dominant, and log-additive) with LC risk was also observed in smokers.

Drinking: Among drinkers, rs28411034 (codominant, dominant, and log-additive) and rs3748682 (codominant, dominant: and log-additive) might confer the increased LC risk. Whereas, rs28411352 (log-additive) was a protective factor against the susceptibility to LC in drinkers.

Relationship of MTF1 variants to clinical characteristics of LC patients in subgroup analysis

The association between MTF1 variants and clinical features of LC patients was examined in subgroups (Table 4 and Supplementary Table 3). Stratified by histological type (Table 4), rs28411034 (codominant, dominant, and log-additive) and rs3748682 (codominant, dominant, and log-additive) contributed to the increasing lung squamous carcinoma risk, but not lung adenocarcinoma. No noteworthy association of MTF1 polymorphisms with stage and lymph node metastasis of LC patients was found (Supplementary Table 3).

Table 4 The association of MTF1 genetic polymorphisms and pathology type risk of lung cancer
Full size table

FPRP analysis

FPRP analysis (Table 5) exhibited the positive results of rs28411034 and rs3748682 for LC susceptibility in the overall analysis with 0.1 prior probability level and FPRP < 0.2. The effects of rs28411034 and rs3748682 on LC risk in men, BMI 24 kg/m2, drinkers, and lung cancer patients were consistent with a significant association of FPRP levels less than 0.2 at a prior probability level of 0.1. In addition, rs3748682 was also significantly associated with a positive outcome for LC risk in smokers, with an FPRP value of < 0.2 despite a prior probability level of 0.1.

Table 5 False-positive report probability for the associations of MTF1 variants with the risk of lung cancer
Full size table

SNP-SNP interaction in the risk of LC

The results of MDR analysis for SNP-SNP interactions showed that rs3748682 was the best single-locus model (cross-validation consistency [CVC]: 10/10; testing balanced accuracy: 0.5328). The best multi-locus model was the combined mode of rs473279, rs28411034, and rs3748682 (CVC: 10/10; testing balanced accuracy: 0.5463, Table 6). The dendrogram (Fig. 2A) showed that these loci with strong interactions were close together on the branches, while these loci with weak interactions were far apart from each other. The Fruchterman-Reingold plot (Fig. 2B) displayed that rs3748682 was the most significant single-locus factor for LC susceptibility, with an information gain of 0.54%.

Table 6 MDR analysis of SNP-SNP interactions
Full size table
Fig. 2
figure 2

SNP-SNP interaction dendrogram (a) and Fruchterman-Reingold (b). Green and blue represent redundancy or association. Values in nodes represent the information gains of individual attribute (main effects). Values between nodes are information gains of each pair of attributes (interaction effects). SNP: single nucleotide polymorphism

Full size image

Association between MTF1 haplotypes and the risk of LC

Moreover, LD and haplotype analysis was performed to estimate the association between MTF1 haplotypes and the risk of LC. As shown in Fig. 3, a high LD block was composed of four MTF1 polymorphisms (rs473279, rs28411034, rs28411352, rs3748682) which formed four haplotypes (Crs473279Ars28411034Crs28411352Crs3748682, Trs473279Grs28411034Crs28411352Trs3748682, Crs473279Grs28411034Crs28411352 Trs3748682, and Crs473279Grs28411034Trs28411352Trs3748682). Furthermore, the haplotype frequency distribution was shown in Table 7. We noted that there was a significant association of Crs473279Grs28411034Crs28411352Trs3748682 (OR 0.76, 95% CI 0.61–0.94, p = 0.013) and Crs473279Grs28411034Trs28411352Trs3748682 (OR 0.79, 95% CI 0.63–0.99, p = 0.045) haplotype with the reduced LC risk.

Fig. 3
figure 3

Haplotype block map for eight SNPs in the MTF1 gene. The numbers inside the diamonds indicate the D′ for pairwise analyses

Full size image
Table 7 Haplotype analysis for the effect of MTF1 haplotypes on the risk of LC
Full size table

Discussion

In our study, we first reported that rs28411034 and rs3748682 tended to have a higher LC susceptibility overall among the Chinese Han population. In the subgroup analysis of demographic characteristics, the effect of rs28411034 and rs3748682 on LC susceptibility was found in males, subjects with BMI ≥ 24 kg/m2, smokers, drinkers, and patients with lung squamous carcinoma. Besides, rs28411352 showed protective effect for reduced LC risk in drinkers. These results might contribute to understanding the pathogenesis of the MTF1 gene in LC progression.

MTF1 is an important transcription factor for heavy metal reactions and is associated with the reduction of oxidative and hypoxic stress in cells [26]. In the presence of p53 in breast cancer cells, MTF1 can be activated by zinc and copper [27]. Previously, MTF1 knockout reduced the proliferation, migration, and invasion of two types of ovarian cancer cells [12]. MTF1 is reported to be related to angiogenesis, and cell invasion [14]. Moreover, in vitro experiments have shown that MTF1 knockdown prevents liver cancer cell proliferation, and promotes cell death [14]. The higher expression of MTF1 may predict a better prognosis for LC patients after chemotherapy [14]. Lung adenocarcinoma cells lacking MTF1 are more sensitive to oxidative stress [15]. Little has been reported about the contribution of MTF1 variants to the susceptibility of tumors. This study is the first to show that rs28411034 and rs3748682 are related to the increased LC susceptibility in the Chinese Han population. Bioinformatics analysis suggested that the possible function of these SNPs might be related to SiPhy cons, promoter/ enhancer histone marks, DNAse, proteins bound, GRASP QTL hits, selected eQTL hits, transcription factor binding, and/or chromatin accessibility peak. Based on the GTEx Portal database, the genotypes of rs473279, rs28411034, rs28411352, and rs3748682 in MTF1 were associated with the mRNA expression in lung tissue. All these SNPs were located in the 3'-UTR region of MTF1 gene. We hypothesize that these SNPs may be related to miRNA binding sites, potentially affecting the expression levels of the MTF1 gene by altering the binding efficiency of miRNA. Our findings imply that these polymorphisms could be implicated in the carcinogenesis of LC by potentially modulating the expression or functionality of the MTF1 gene, thereby laying a theoretical foundation for future mechanistic investigations.

With the population aging, LC incidence is on the rise [28]. Internationally, the incidence of LC is higher in men than in women [29, 30]. Recent studies have shown that in most countries, morbidity among women is increasing, while morbidity and mortality among men are declining [31]. These studies suggest that age and sex are confounding factors in the genetic association with LC. MTF1 rs28411034 and rs3748682 were related to LC risk in males. Moreover, rs473279 was associated with an increased LC susceptibility in subjects aged 51–60 years. BMI is a risk factor for LC and is negatively correlated with LC risk [32]. According to the Chinese BMI classification [33, 34], overweight was defined as ≥ 24 kg/m2. Therefore, we used BMI of 24 kg/m2 as a cut-off. Here, we found that rs28411034 and rs3748682 may contribute to LC risk in the subjects with BMI ≥ 24 kg/m2. Smoking is a well-known major risk factor for LC [35]. The role of alcohol consumption in LC occurrence is controversial. Some studies have shown that alcohol consumption increases LC risk [36], and others displayed that moderate alcohol consumption is a protective factor in LC development [37]. Stratified results showed that rs28411034 and rs3748682 were associated with LC risk in smokers and drinkers. Moreover, rs28411352 was related to the reduced LC risk in drinkers. Our study demonstrates the synergistic role of MTF1 polymorphisms and demographic characteristics in the occurrence of LC. Our results propose that the association between MTF1 polymorphism and LC risk are age-, sex-, BMI-, smoking-, and drinking-specific.

Lung cancer, a multifaceted condition influenced by genetic and environmental factors, may have its risk factors uncovered through polygenic or SNP-SNP interaction studies. Notably, the MDR method is recognized for its effectiveness in identifying SNP-SNP/gene–gene interactions in the absence of individual gene effects, which is crucial in case–control studies of complex diseases [38] In this study, we applied MDR to analyze the interactions among four specific SNPs. Our findings indicate that rs3748682 was the best single-locus model and the best multi-locus model was the combined mode of rs473279, rs28411034, and rs3748682.

Haplotypes, as fundamental genetic variations and units of inheritance, can influence phenotypes either directly by altering promoter activity and protein structure or indirectly through linkage with nearby untyped causal variations [39]. Consequently, haplotype association studies are vital for elucidating the etiology of complex phenotypes. In our research, haplotype analysis displayed the significant association of Crs473279Grs28411034Crs28411352Trs3748682 and Crs473279Grs28411034Trs28411352Trs3748682 haplotypes with the reduced LC risk. These hinted us that MTF1 haplotypes could be a potential risk factor for LC occurrence.

There are some limitations to our study. First, all the samples were from one hospital, so there was some selection bias and it was not representative of the entire population. Therefore, further large-scale genetic studies in different populations are required to verify our findings. Second, our study only assesses the association between four SNPs in MTF1 and LC risk, and a large number of loci remain to be studied. Third, this study discovered that rs28411034 and rs3748682 were related to LC risk, but their specific mechanisms in lung carcinogenesis need to be further explored in a complete functional experiment. Fourth, none of the SNPs was significant after the Bonferroni correction (p < 0.05/(4 × 4 × 5)). This may be in that the Bonferroni correction adjusts the value of alpha based on the number of tests performed and is thus conservative; in some cases, truly significant differences may be deemed non-significant as a result of type II errors. Therefore, we did not use the Bonferroni correction in our study.Whatever, this study has exploratory value for the relationship between MTF1 polymorphisms and the risk of LC. Our findings need to be further confirmed in a larger sample.

Conclusion

Taken together, our study suggests that MTF1 rs28411034 and rs3748682 tended to be associated with increased LC susceptibility among the Chinese Han population, especially in males, subjects with BMI ≥ 24 kg/m2, smokers, drinkers, and patients with lung squamous carcinoma. These findings help to shed light on the underlying mechanism of the MTF1 gene in LC progression. However, functional studies and larger patient population studies are needed to confirm our conclusions.

Availability of data and materials

All data generated or analyzed during this study are included in this published article.

The datasets generated and/or analysed during the current study are available in the Zenodo repository, https://zenodo.org/record/8384945.

References

  1. Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, Bray F. Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2021;71(3):209–49.

    Article  PubMed  Google Scholar 

  2. Nooreldeen R, Bach H. Current and future development in lung cancer diagnosis. Int J Mol Sci. 2021;22(16):8661.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Chen P, Liu Y, Wen Y, Zhou C. Non-small cell lung cancer in China. Cancer Commun (Lond). 2022;42(10):937–70.

    Article  PubMed  Google Scholar 

  4. Gao S, Li N, Wang S, Zhang F, Wei W, Li N, Bi N, Wang Z, He J. Lung cancer in People’s Republic of China. J Thorac Oncol. 2020;15(10):1567–76.

    Article  PubMed  Google Scholar 

  5. Dai J, Shen W, Wen W, Chang J, Wang T, Chen H, Jin G, Ma H, Wu C, Li L, et al. Estimation of heritability for nine common cancers using data from genome-wide association studies in Chinese population. Int J Cancer. 2017;140(2):329–36.

    Article  CAS  PubMed  Google Scholar 

  6. Sampson JN, Wheeler WA, Yeager M, Panagiotou O, Wang Z, Berndt SI, Lan Q, Abnet CC, Amundadottir LT, Figueroa JD, et al. Analysis of heritability and shared heritability based on genome-wide association studies for thirteen cancer types. J Natl Cancer Inst. 2015;107(12):djv279.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Dai J, Lv J, Zhu M, Wang Y, Qin N, Ma H, He YQ, Zhang R, Tan W, Fan J, et al. Identification of risk loci and a polygenic risk score for lung cancer: a large-scale prospective cohort study in Chinese populations. Lancet Respir Med. 2019;7(10):881–91.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Zhou W, Liu G, Hung RJ, Haycock PC, Aldrich MC, Andrew AS, Arnold SM, Bickeböller H, Bojesen SE, Brennan P, et al. Causal relationships between body mass index, smoking and lung cancer: univariable and multivariable Mendelian randomization. Int J Cancer. 2021;148(5):1077–86.

    Article  CAS  PubMed  Google Scholar 

  9. Wang C, Dai J, Qin N, Fan J, Ma H, Chen C, An M, Zhang J, Yan C, Gu Y, et al. Analyses of rare predisposing variants of lung cancer in 6,004 whole genomes in Chinese. Cancer Cell. 2022;40(10):1223-1239.e1226.

    Article  CAS  PubMed  Google Scholar 

  10. Si M, Lang J. The roles of metallothioneins in carcinogenesis. J Hematol Oncol. 2018;11(1):107.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Helal NS, Moussa MM, Aboushousha T, Aboul-Ezz MA, Magdy M. Differential expression of metallothionein and p21 in gastric cancer and some precursor lesions. Eur Rev Med Pharmacol Sci. 2022;26(9):3100–8.

    CAS  PubMed  Google Scholar 

  12. Ji L, Zhao G, Zhang P, Huo W, Dong P, Watari H, Jia L, Pfeffer LM, Yue J, Zheng J. Knockout of MTF1 inhibits the epithelial to mesenchymal transition in ovarian cancer cells. J Cancer. 2018;9(24):4578–85.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. He J, Jiang X, Yu M, Wang P, Fu L, Zhang G, Cai H. MTF1 has the potential as a diagnostic and prognostic marker for gastric cancer and is associated with good prognosis. Clin Transl Oncol. 2023;25(11):3241-51.

  14. Song L, Zeng R, Yang K, Liu W, Xu Z, Kang F. The biological significance of cuproptosis-key gene MTF1 in pan-cancer and its inhibitory effects on ROS-mediated cell death of liver hepatocellular carcinoma. Discov Oncol. 2023;14(1):113.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Liu Y, Lin W, Yang Y, Shao J, Zhao H, Wang G, Shen A. Role of cuproptosis-related gene in lung adenocarcinoma. Front Oncol. 2022;12:1080985.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Sobh A, Loguinov A, Yazici GN, Zeidan RS, Tagmount A, Hejazi NS, Hubbard AE, Zhang L, Vulpe CD. Functional profiling identifies determinants of arsenic trioxide cellular toxicity. Toxicol Sci. 2019;169(1):108–21.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Shi Y, Amin K, Sato BG, Samuelsson SJ, Sambucetti L, Haroon ZA, Laderoute K, Murphy BJ. The metal-responsive transcription factor-1 protein is elevated in human tumors. Cancer Biol Ther. 2010;9(6):469–76.

    Article  CAS  PubMed  Google Scholar 

  18. Magda D, Lecane P, Wang Z, Hu W, Thiemann P, Ma X, Dranchak PK, Wang X, Lynch V, Wei W, et al. Synthesis and anticancer properties of water-soluble zinc ionophores. Cancer Res. 2008;68(13):5318–25.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. GTEx Consortium. The Genotype-Tissue Expression (GTEx) project. Nat Genet. 2013;45(6):580–5.

  20. Ward LD, Kellis M. HaploReg v4: systematic mining of putative causal variants, cell types, regulators and target genes for human complex traits and disease. Nucleic Acids Res. 2016;44(D1):D877-881.

    Article  CAS  PubMed  Google Scholar 

  21. Dong S, Boyle AP. Predicting functional variants in enhancer and promoter elements using RegulomeDB. Hum Mutat. 2019;40(9):1292–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Ellis JA, Ong B. The MassARRAY(®) system for targeted SNP genotyping. Methods Mol Biol. 2017;1492:77–94.

    Article  CAS  PubMed  Google Scholar 

  23. Wacholder S, Chanock S, Garcia-Closas M, El Ghormli L, Rothman N. Assessing the probability that a positive report is false: an approach for molecular epidemiology studies. J Natl Cancer Inst. 2004;96(6):434–42.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Motsinger AA, Ritchie MD. The effect of reduction in cross-validation intervals on the performance of multifactor dimensionality reduction. Genet Epidemiol. 2006;30(6):546–55.

    Article  PubMed  Google Scholar 

  25. Zhang H, Zhang Z, Zhang J, Xu L, Ye Z, Hao Y, Cai B, Zhou S, Liu K, Sun L, et al. Fine-mapping of ABO gene identifies two novel SNPs associated with large artery atherosclerotic stroke in a Chinese Han population. Mol Neurobiol. 2017;54(3):2107–13.

    Article  CAS  PubMed  Google Scholar 

  26. Tavera-Montañez C, Hainer SJ, Cangussu D, Gordon SJV, Xiao Y, Reyes-Gutierrez P, Imbalzano AN, Navea JG, Fazzio TG, Padilla-Benavides T. The classic metal-sensing transcription factor MTF1 promotes myogenesis in response to copper. FASEB journal : official publication of the Federation of American Societies for Experimental Biology. 2019;33(12):14556–74.

    Article  PubMed  Google Scholar 

  27. Ostrakhovitch EA, Olsson PE, von Hofsten J, Cherian MG. P53 mediated regulation of metallothionein transcription in breast cancer cells. J Cell Biochem. 2007;102(6):1571–83.

    Article  CAS  PubMed  Google Scholar 

  28. Cao W, Chen HD, Yu YW, Li N, Chen WQ. Changing profiles of cancer burden worldwide and in China: a secondary analysis of the global cancer statistics 2020. Chin Med J. 2021;134(7):783–91.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Barta JA, Powell CA, Wisnivesky JP. Global epidemiology of lung cancer. Ann Glob Health. 2019;85(1):8.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Stapelfeld C, Dammann C, Maser E. Sex-specificity in lung cancer risk. Int J Cancer. 2020;146(9):2376–82.

    Article  CAS  PubMed  Google Scholar 

  31. Huang J, Deng Y, Tin MS, Lok V, Ngai CH, Zhang L, Lucero-Prisno DE 3rd, Xu W, Zheng ZJ, Elcarte E, et al. Distribution, risk factors, and temporal trends for lung cancer incidence and mortality: a global analysis. Chest. 2022;161(4):1101–11.

    Article  PubMed  Google Scholar 

  32. You D, Wang D, Wu Y, Chen X, Shao F, Wei Y, Zhang R, Lange T, Ma H, Xu H, et al. Associations of genetic risk, BMI trajectories, and the risk of non-small cell lung cancer: a population-based cohort study. BMC Med. 2022;20(1):203.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Pan XF, Wang L, Pan A. Epidemiology and determinants of obesity in China. Lancet Diabetes Endocrinol. 2021;9(6):373–92.

    Article  PubMed  Google Scholar 

  34. Zhou BF. Predictive values of body mass index and waist circumference for risk factors of certain related diseases in Chinese adults–study on optimal cut-off points of body mass index and waist circumference in Chinese adults. Biomedical and environmental sciences : BES. 2002;15(1):83–96.

    PubMed  Google Scholar 

  35. Schuller HM. The impact of smoking and the influence of other factors on lung cancer. Expert Rev Respir Med. 2019;13(8):761–9.

    Article  CAS  PubMed  Google Scholar 

  36. Bunjaku J, Lama A, Pesanayi T, Shatri J, Chamberlin M, Hoxha I. Lung cancer and lifestyle factors: umbrella review. Hematol Oncol Clin North Am. 2024;38(1):171–84.

    Article  PubMed  Google Scholar 

  37. Chen C, Hu Q, Wang J, Wen T, Zhu C, Tan W, Chen X, Zhao Q, Wang W, Cao H, et al. Habitual consumption of alcohol with meals and lung cancer: a Mendelian randomization study. Annals of translational medicine. 2021;9(3):263.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Ritchie MD, Hahn LW, Roodi N, Bailey LR, Dupont WD, Parl FF, Moore JH. Multifactor-dimensionality reduction reveals high-order interactions among estrogen-metabolism genes in sporadic breast cancer. Am J Hum Genet. 2001;69(1):138–47.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Servin B, Stephens M. Imputation-based analysis of association studies: candidate regions and quantitative traits. PLoS Genet. 2007;3(7):e114.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

The authors thank the First Affiliated Hospital of Xi’an Jiaotong University for providing samples for our study.

Funding

This work was supported by the Applied Basic Research Foundation of Yunnan Province (CN) [grant number: 2017FE468 (-125) and 202001AY070001-111], the Open Project of the Clinical Medicine Center of the First People's Hospital of Yunnan Province [grant number: 2021LCZXXF-XY12] and CAMS Innovation Fund for Medical Sciences (CIFMS) [grant number: 2016-I2M-3–024], Yunnan Province Clinical Research Center for Hematologic Disease [grant number: 2022LCZXKF-XY04], Yunnan Province Clinical Research Center for Hematologic Disease (2023YJZX-XY02) and the Kunming University of Science and Technology & the First People's Hospital of Yunnan Province Joint Special Project on Medical Research (KUST-KH2022028Y), Yunnan Fundamental Research Projects (202401AY070001-109), Yunnan Provincial Clinical Medical Center for Blood Diseases and Thrombosis Prevention and Treatment.

Author information

Authors and Affiliations

  1. Department of Respiratory Medicine, The First Affiliated Hospital of School of Medicine of Xi’an Jiaotong University, Yanta District, No. 277, Yanta West Road, Xi’an, 710061, Shaanxi, China

    Yujing Cheng & Mingwei Chen

  2. Department of Blood Transfusion, The First People’s Hospital of Yunnan Province, The Afiliated Hospital of Kunming University of Science and Technology, Kunming, 650032, Yunnan, China

    Yujing Cheng, Chan Zhang, Qi Li, Xin Yang & Wanlu Chen

  3. Department of Blood Transfusion, The First People’s Hospital of Qujing City, Qujing, 655099, Yunnan, China

    KunHua He

Authors
  1. Yujing Cheng
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Chan Zhang
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Qi Li
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Xin Yang
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Wanlu Chen
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. KunHua He
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Mingwei Chen
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

YC and MC designed this study protocol and drafted the manuscript; CZ and QL performed the DNA extraction and genotyping; XY and WC performed the data analysis; YC and KH polished the manuscript, MC conceived and supervised the study. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Mingwei Chen.

Ethics declarations

Ethics approval and consent to participate

Approval for this study was taken from the Ethics Committee of the First Affiliated Hospital of Xi’an Jiaotong University and the written informed consent was signed by each participant. All experiments were in line with the Declaration of Helsinki.

Consent for publication

Not applicable (NA).

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

12885_2024_12516_MOESM1_ESM.tif

Supplementary Material 1: Suppl_Figure 1. The violin plot for the association between the genotypes of MTF1 variants and the mRNA expression in the lung tissue. Data were from GTEx Portal database (https://gtexportal.org/home/).

Supplementary Material 2.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Cheng, Y., Zhang, C., Li, Q. et al. MTF1 genetic variants are associated with lung cancer risk in the Chinese Han population. BMC Cancer 24, 778 (2024). https://doi.org/10.1186/s12885-024-12516-y

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12885-024-12516-y

Keywords

BMC Cancer

ISSN: 1471-2407

Contact us