An Investigation of the Effects of B7-H4 Gene rs10754339 and miR-125a Gene rs12976445 on Cancer Susceptibility

JIN Yu Chen ,DONG Li Juan ,YANG Qin Yue ,XIONG Wei Ning ,WANG Wei Yi ,FENG Xian Hong ,YU Wei,HUANG Wei,#,and CHEN Bi Feng,#

1.Department of Biological Science and Technology,School of Chemistry,Chemical Engineering and Life Sciences,Wuhan University of Technology,Wuhan 430081,Hubei,China;2.Department of Biotechnology,Institute of WUTAMU,Wuhan University of Technology,Wuhan 430081,Hubei,China;3.Department of Clinical Laboratory,Wuhan Xinzhou District People’s Hospital,Wuhan 430400,Hubei,China;4.Department of Air Crew Service,School of Civil Aviation,Wuhan Business and Trade Vocational College,Wuhan 430062,Hubei,China

Abstract Objective To investigate the effects of the B7-H4 gene rs10754339 and miR-125a gene rs12976445 on cancer susceptibility through a case-control study and meta-analysis.Methods A total of 1,490 cancer patients (lung/gastric/liver/: 550/460/480) and 800 controls were recruited in this case-control study.The meta-analysis was performed by pooling the data from previous related studies and the present study.Results The results of this study showed that in the Hubei Han Chinese population,the rs10754339 gene was significantly associated with the risk of lung and gastric cancer but not liver cancer,and the rs12976445 gene was significantly associated with the risk of lung cancer but not liver or gastric cancer.The meta-analysis results indicated that rs10754339 and rs12976445 contributed to cancer susceptibility in the Chinese population and also revealed a significant association between rs10754339 and breast cancer risk,as well as between rs12976445 and lung cancer risk.Conclusion The B7-H4 gene rs10754339 and miR-125a gene rs12976445 may be the potential genetic markers for cancer susceptibility in the Chinese population,which should be validated in future studies with larger sample sizes in other ethnic populations.

Key words: B7-H4 gene;miR-125a gene;rs10 754339;rs12 976445;Cancer susceptibility

INTRODUCTION

Cancer has become a serious threat to public health and the economy worldwide[1].Lung cancer has been the highest cause of cancer-related morbidity and mortality,closely followed by gastric cancer and liver cancer[2].Cancer is a complex disease involving multiple factors,including environmental factors(e.g.,alcohol intake,tobacco use) and genetic factors,as well as their interactions[3].Although clinical treatment of cancer has improved greatly,the overall outcome is still unsatisfactory,and the five-year survival rate is still low[4].Thus,the identification of reliable genetic biomarkers is of great importance for risk prediction and early detection of cancer.

Immunity plays an important role in carcinogenesis,and chronic inflammation is recognized as a condition favoring the development of cancer[5].The B7-H4 is an immune checkpoint that plays an essential role in carcinogenesis,and it is able to affect T-cell functions (cytotoxicity,cytokine secretion,and cell proliferation),epithelial-tomesenchymal transition,and cancer stemness[6-9].The miR-125a is an inflamma-miRNA involved in the regulation of innate and adaptive immune responses and regulates immune and inflammation processes in human cancers[10].Aberrant expressedB7-H4gene andmiR-125agene occur in many malignant tumors,whereas the underlying mechanism remains largely unknown.

Single nucleotide polymorphism (SNP),a welldefined genetic variation,has been extensively applied in the medical field.Regulatory SNPs (rSNPs)are a group of SNPs that participate in the regulation of gene expression and protein structural/functional behavior[11].Among the enormous SNPs in theB7-H4gene andmiR-125agene,two rSNPs,namely,rs10754339[12]and rs12976445[13],have been highlighted in the literature[14-23].However,the results were contradictory rather than conclusive.To further explore their contribution to carcinogenesis,we conducted a case-control study to investigate the effects of rs12976445 and rs10754339 on susceptibility to lung,liver,and gastric cancer in the Chinese population of Hubei Province.Additionally,we estimated the association of rs12976445/rs10754339 and cancer susceptibilityviaa metaanalysis by combining the data of previous related studies and those of the present study.

MATERIALS AND METHODS

Studied Subjects

A total of 1,490 cancer patients (lung/gastric/liver/: 550/460/480) and 800 controls were enrolled in this study.All cancer patients were confirmed histopathologically and recruited from the Hubei Cancer Hospital and Wuhan Xinzhou District People’s Hospital.Individuals to be included in the control were selected from cancer-free individuals who visited the Wuhan Xinzhou District People’s Hospital for regular physical examinations.All participants were biologically unrelated Han Chinese living in the Hubei Province of China.Each participant signed an informed consent form.This study was approved by the ethical committee of the Wuhan University of Technology.

Preparation of DNA Samples

The peripheral blood samples of all studied subjects were collected into blood vacuum tubes containing ethylenediaminetetraacetic acid (EDTA)and stored at 4 °C.Genomic DNA was extracted from the blood samples using the TIANamp Blood DNA Kit(DP348;TianGen Biotech,Beijing,China) within one week,according to the manufacturer’s instructions.Finally,the DNA samples were stored at -20 °C before use.

Genotyping of rs10754339 and rs12976445

The genotyping of rs10754339 and rs12976445 was performed by the polymerase chain reactionrestriction fragment length polymorphism (PCR-RFLP)method.Table 1 shows the PCR primer sequences and the restriction enzymes used.The accuracy of the PCRRFLP was confirmed by sequencing the PCR products from 10% randomly selected samples,and there were no differences in results between the two methods.Supplementary Figure S1 (available in www.besjournal.com) displays the genotyping diagrams of rs10754339 and rs12976445 by PCR-RFLP and Sanger sequencing.

Table 1.Oligonucleotides and restriction enzymes of rs10754339 and rs12976445 for genotyping

Statistical Analysis

SPSS 15.0 was used to perform all the statistical analyses.χ2test was applied to determine the statistical differences in age,gender,smoking status,and drinking status between cancer patients and the control group.The genotypic frequencies of rs10754339 and rs12976445 were tested for deviation from the Hardy-Weinberg equilibrium(HWE) in the control group.Furthermore,the logistic regression analysis under six genetic comparisons was used to estimate the effects of rs10754339 and rs12976445 on cancer susceptibility.The level of significance was set atP＜ 0.05,and the Bonferroni correction was applied for multiple comparisons (P＜0.0084,0.05/6).

Meta-Analysis

Related literature up to March of 2022 was searched in PubMed,ISI Web of Science,CNKI,and Wanfang databases with the terms “rs10754339 and cancer/tumor”or “rs12976445 and cancer/tumor.”Supplementary Figure S2 (available in www.besjournal.com) shows the flow diagram of the literature retrieving process.References listed in retrieved articles were checked for missing information.For a study to be included in the metaanalysis,it had to meet the following criteria: (1) it is a human study,(2) it is an investigation of rs10754339/rs12976445 and cancer risk,(3) it is a case-control study,(4) frequency data are extractable,and (5) the control group conforms to the HWE.The STATA 14.0 was used to conduct all the statistical analyses,including the heterogeneity test,pooledORs,publication bias,and sensitivity analysis.P＜ 0.05 was set as the significance level,and Bonferroni correction was applied for multiple comparisons (P＜ 0.0084,0.05/6).Trial sequential analysis (TSA) was used to assess whether the evidence of this meta-analysis was enough and conclusive[24],and the interpretation of TSA results has been described in our previous article[25].

RESULTS

Principal Characteristics of the Studied Subjects

The principal characteristics of the age,gender,smoking status,and alcohol status of studied subjects are presented in Table 2.The distributions of these four variables between cancer patients and the control group did not differ significantly,suggesting the cases and controls were wellmatched in this study.

Table 2.Characteristics of participants in this study

Effects of rs10754339 and rs12976445 on Cancer Susceptibility

Table 3 shows the allelic and genotypic distributions of rs10754339 and rs12976445,as well as their contribution to cancer susceptibility.No deviation from HWE was observed for rs10754339 or rs12976445 in the control group (P＞ 0.05).In this study,rs10754339 was significantly associated with the risk of lung and gastric cancer but not with liver cancer.The A allele of rs10754339 was associated with a lower risk of lung cancer (Avs.G,P=0.002,OR=0.68,95%CI: 0.54–0.87) and gastric cancer (Avs.G,P=0.006,OR=0.71,95%CI: 0.55–0.91) than the G allele.Concordantly,carriers of the AA genotype were less inclined to suffer from lung cancer than carriers of the AG/AG+GG genotype(AAvs.AG,P=0.006,OR=0.69,95%CI: 0.53–0.90;AAvs.AG+GG,P=0.003,OR=0.67,95%CI:0.52–0.87),and were less inclined to suffer from gastric cancer than AG+GG genotype carriers (AAvs.AG+GG,P=0.007,OR=0.68,95%CI:0.52–0.90).Additionally,the rs12976445 was significantly associated with the risk of lung cancer but not liver or gastric cancer.Specifically,the T allele of rs12976445 conferred a higher risk for lung cancer than the C allele (Tvs.C,P=0.001,OR=1.48,95%CI: 1.17–1.87),and individuals carrying the genotype with at least one T allele had a higher lung cancer risk than that CC carriers (TT+TCvs.CC,P=0.004,OR=1.46,95%CI: 1.13–1.88).

Meta-Analysis of the Association between rs10754339 and Cancer Susceptibility

Supplementary Table S1 (available in www. besjournal.com) shows the main features of the included studies in the meta-analysis of rs10754339 and cancer susceptibility.The statistical results are presented in Table 4 and Supplementary Figure S3(available in www.besjournal.com).It was observed that rs10754339 was significantly associated with total cancer risk under Avs.G (P=0.004,OR=0.78,95%CI: 0.66–0.93),AAvs.AG (P=0.004,OR=0.76,95%CI: 0.62–0.91),and AAvs.AG+GG (P=0.004,OR=0.75,95%CI: 0.62–0.91).Interestingly,the ethnicity-stratified analysis showed that the association between rs10754339 and total cancer risk especially existed in the Chinese population.Furthermore,there was a significant association between rs10754339 and breast cancer risk.

Table 4.Meta-analysis of the association between rs10754339 and cancer risk

Meta-Analysis of the Association between rs12976445 and Cancer Susceptibility

The main features of the included studies in the meta-analysis of rs12976445 and cancer susceptibility are presented in Supplementary Table S1.The statistical results are shown in Table 5 and Supplementary Figure S4 (available in www.besjournal.com).No significant association between rs12976445 and total cancer risk was observed.However,the stratified analysis by ethnicity revealed that rs12976445 was significantly associated with total cancer risk in the Chinese population under Tvs.C,TTvs.TC,TTvs.CC,TTvs.TC+CC and TT+TCvs.CC.Moreover,the stratified analysis by cancer type identified a significant association between rs12976445 and lung cancer risk under TTvs.TC,TTvs.CC and TTvs.TC+CC.

Table 5.Meta-analysis of the association between rs12976445 and cancer risk

Sensitivity Analysis,Publication Bias,and TSA

Removal of any single case-control study from the analysis of rs10754339 and rs12976445 in the allele model did not significantly affect the consistency of pooledORs (Figures 1–2).Moreover,potential publication bias was assessed for rs10754339 and rs12976445,and the results did not show any evidence of publication bias in any of the genetic comparisons (Tables 4–5),suggesting the statistical results were credible.As shown in Figures 3–4,the cumulative Z-curves crossed both the traditional boundary (Z=1.96) and trial sequential monitoring boundary in the TSA of rs10754339 and cancer risk,as well as in the TSA of rs12976445 and cancer risk,indicating that true positive results could be obtained.

DISCUSSION

B7-H4,an immune checkpoint,and miR-125a,an immunoregulator,have attracted much attention in cancer immunology research.The rs10754339,located in the 3'-UTR of theB7-H4gene,has been reported to impair the regulation of miR-506-3p on B7-H4[12].The rs12976445 within the promoter region ofmiR-125aaffects miR-125a expression,and bioinformatic analysis predicted a possible GATA-1 binding site by the presence of rs12976445[13].Therefore,we evaluated the effects of the rs10754339 and rs12976445 genotypes on the expression levels of theB7-H4gene and themiR-125agene,respectively.Expression quantitative trait loci (eQTL) analysis on the GTEx portal website(http://www.gtexportal.org/) revealed thatB7-H4mRNA differentially expressed in pancreas tissues depending on the three genotypes (P＜ 0.001) with AA tissues showing the highest expression level.ThemiR-125amRNA is differentially expressed in whole blood depending on the three genotypes (P＜ 0.001),with TT samples showing the highest expression level (Supplementary Figure S5,available in www.besjournal.com).These findings indicated that theB7-H4gene rs10754339 and themiR-125agene rs12976445 are two functional variants that may affect individual susceptibility to cancer.

The association of theB7-H4gene rs10754339 polymorphism with breast cancer has been repeatedly investigated[14-17],and studies have recorded a significant positive association except for ?zg?z et al.[16].This discrepancy might be mainly attributed to the inadequate power in some studies caused by insufficient sample size.Indeed,in the study of ?zg?z et al.(31 breast cancer patients/30 healthy women)[16],rs10754339 showed an association in the same direction (ORs ＞ 1 for G allele and AG genotype) as reported in the other three studies[14,15,17].Additionally,?zg?z et al.revealed that rs10754339 was a genetic marker for bladder cancer[18].This study has added to the list of cancers associated with rs10754339 by showing that rs10754339 was significantly associated with the risk of lung and gastric cancer.The link between rs10754339 and lung cancer susceptibility might be attributed to the G ＞ A transition inhibiting miR-506-3p duringB7-H4expression[11],leading to the aberrant expression of B7-H4,which finally leads to the occurrence of lung cancerviaAMPK/mTOR signaling[9,26].Although statistically not significant,we observed a similar tendency that the A allele and AA genotype were more frequent in liver cancer patients than in normal individuals,suggesting that rs10754339 has the potential to be a predictive biomarker for liver cancer,but a confirmatory study with a greater number of participants is required.

Figure 1.Sensitivity analysis of rs10754339 and overall cancer risk under A vs. G in the total population(A),Chinese population (B) and Caucasian population (C).Sensitivity analysis of rs10754339 and breast cancer risk under A vs.G (D).

Previous studies have explored the contribution ofmiR -125agene rs12976445 polymorphism to colorectal cancer (CRC)[19],prostate cancer(PCa)[20,21],breast cancer[22],and lung cancer[23].Consistent with the findings of Sun et al.[23],we also found a connection between rs12976445 and lung cancer risk.However,we found no significant association of rs12976445 with susceptibility to CRC,breast cancer,liver cancer,or gastric cancer.Interestingly,Hakimian et al.did not find an association between rs12976445 and PCa susceptibility[21],whereas Damodaran et al.came to the opposite conclusion[20].Two possible reasons may explain the discrepancy.First,the sample sizes of these two studies were too small (300 subjects and 200 subjects,respectively),which might have reduced the statistical power and increased the margin of error.Therefore,their findings may not be completely convincing.Secondly,the different living environments,lifestyles,and genetic backgrounds of different ethnic descents (Iranians and Indians)might also have led to the discrepancy.

To improve the statistical power and resolve the discordant results of previous studies,a metaanalysis was conducted to further evaluate the effects of rs10754339 and rs12976445 on cancer susceptibility.We found that rs10754339 was significantly associated with total cancer risk,especially in the Chinese population,and a strong connection between rs10754339 and breast cancer risk was also observed.Additionally,rs12976445 was significantly associated with total cancer risk only in the Chinese population.The cancer type-stratified analysis showed that rs12976445 was closely related to lung cancer risk,which confirmed the findings of our case-control study.However,no association was observed between rs12976445 and prostate cancer risk and between rs12976445 and total cancer risk based on the PCR-RFLP method.

There are several limitations of this study.First,the studied subjects were all enrolled in the hospital,which could not fully rule out selection bias.Secondly,our findings on the association between rs10754339/rs12976445 and the risk of liver,lung,and gastric cancer apply only to the Hubei Han Chinese;thus,further confirmatory studies in other ethnic populations are required.Thirdly,although rs12976445 was predicted to regulate the expression ofmiR -125aby affecting the binding affinity of GATA-1,the direct relationship between rs12976445 and miR-125a expression still needs to be validated by a functional study.Finally,owing to the limited number of studies for specific cancer types,the cancer type-stratified analysis was only conducted for rs10754339 and breast cancer,rs12976445 and prostate cancer,and rs12976445 and lung cancer.Therefore,replication studies and updated metaanalysis are needed to gain a better understanding of the correlation between rs10754339/rs12976445 and cancer susceptibility.

Figure 2.Sensitivity analysis of rs12976445 and overall cancer risk under T vs. C in the total population(A),in the Asian population (B),based on PCR-RFLP (C),in the Chinese population (D) and in the Iranian population (E).Sensitivity analysis of rs12976445 and specific cancer risk (F) lung cancer;(G) prostate cancer.

In conclusion,our case-control study indicated that theB7-H4gene rs10754339 is associated with the susceptibility to lung and gastric cancer in the Chinese population and that themiR -125agene rs12976445 may be a genetic factor of lung cancer in the Chinese population.Meta-analysis indicated that both rs10754339 and rs12976445 may be contributing to cancer susceptibility in the Chinese population.However,additional studies with larger sample sizes in different ethnic populations are required to confirm our present findings.

Figure 3.Trial sequential analysis of rs10754339 and overall cancer risk in the total population [(A) A vs.G,(B) AA vs.AG,and (C) AA vs.AG+GG] and the Chinese population [(D) A vs.G,(E) AA vs.AG,and (F) AA vs.AG+GG],trial sequential analysis of rs10754339 and breast cancer risk [(G) A vs.G,(H) AA vs.AG,and(I) AA vs.AG+GG].

Figure 4.Trial sequential analysis of rs12976445 and overall cancer risk in the Chinese population under T vs.C (A),TT vs.CC (B),TT+TC vs.CC (C),and TT vs.TC+CC (D).

CONFLICT OF INTEREST

The authors declare that they have no conflict of interest in this work.

ACKNOWLEDGEMENTS

The authors thank all the participants and investigators enrolled in this study.This work was supported by the Fundamental Research Funds for the Central Universities (WUT: 2020IB029).

AUTHOR CONTRIBUTIONS

CHEN Bi Feng and HUANG Wei designed and supervised the study.FENG Xian Hong collected the blood samples.JIN Yu Chen,DONG Li Juan,and YANG Qin Yue performed the research study and collected the data.XIONG Wei Ning,WANG Wei Yi,and Yu Wei analyzed the data.JIN Yu Chen and DONG Li Juan wrote the manuscript.CHEN Bi Feng and HUANG Wei revised the manuscript critically and made final approval of the manuscript.