A Comprehensive Study of the Association between LEPR Gene rs1137101 Variant and Risk of Digestive System Cancers

Digestive system (DS) cancers include tumors of the mouth, esophagus, liver, stomach, bile duct, colon, rectum, pancreas, and gallbladder. GLOBOCAN estimates that DS cancers accounted for approximately 19.3 million newly diagnosed cancers and 10.0 million deaths worldwide in 2020. Colorectal, gastric, liver, and esophageal cancers are the four most common DS cancers, ranked 3rd, 5th, 6th, and 7th in global incidence and 2nd, 4th, 3rd, and 6th in overall mortality, respectively^[1]. Although novel therapeutic strategies and precise diagnostics have been developed, the outcome for patients with DS cancer remains unsatisfactory^[2]. More importantly, due to the largely asymptomatic nature of most types of DS cancers, patients are often diagnosed at a late stage, resulting in limited treatment options and poor prognosis.

Leptin is a 16-kD cytokine-like peptide hormone secreted primarily by the white adipose tissue^[3]. It exerts its biological functions through a specific homodimeric receptor called leptin receptor (LEPR), which is expressed in multiple organs^[4]. The long form of the leptin receptor, Ob-Rb, is the predominant isoform of LEPR. Interestingly, activation of Ob-Rb by leptin binding has been shown to stimulate the janus kinase 2/signal transducer and activator of transcription 3 (JAK2/STAT3), phosphoinositide 3-kinase/mammalian target of rapamycin (PI3K/mTOR), and mitogen-activated protein kinase/extracellular signal-regulated kinase (MAPK/ERK) signaling pathways^[5-8]. Numerous lines of evidence suggest that these key signaling cascades play crucial roles in multiple aspects of tumorigenesis; thus, the LEPR gene is regarded as a plausible candidate gene for cancer susceptibility. Indeed, functional studies have revealed that LEPR is involved in promoting angiogenesis, facilitating cell proliferation and inhibiting epithelial cell apoptosis^[9-11].

LEPR is located on chromosome 1p31 (position_38: 65420652-65637493). More than 3,117 genetic variants have been identified in the LEPR gene region. Current knowledge suggests that genetic variants with regulatory function are more likely to be the causative loci of complex human diseases. Matsuoka et al. reported seven common missense variants (Lys109Arg, Gln223Arg, Ser343Ser, Ser492Thr, Lys656Asn, Ala976Asp, and Pro1019Pro) that cause amino acid changes in the Ob-Rb coding region^[12]. Of the seven genetic variants, Gln223Arg (rs1137101) has been extensively examined for its association with cancer risk. The interesting question is why rs1137101 has attracted the most attention. Here, we used the SNPinfo software (https://snpinfo.niehs.nih.gov/) to predict that rs1137101 may act as an exonic splicing enhancer or silencer of LEPR mRNA, which modifies the efficiency of RNA splicing and thereby alters the function of LEPR.

The genetic effect of rs1137101 has been investigated mainly in DS cancers (oral^[13-16], esophageal^[17,18], gastric^[19-21], colorectal^[22-27], and liver^[28-31]) and breast cancer^[32], as well as in a few other types (endometrial^[33] and lung cancer^[34]). Nevertheless, no conclusive results have been published, which calls for further investigation. To replicate and further explore the association between rs1137101 and risk of DS cancer, we conducted a case–control study to investigate the genetic effect of rs1137101 on the risk of gastric, colorectal, and liver cancers in a Hubei Chinese population.

On a different note, although the association between rs1137101 and overall cancer risk has been previously determined in several meta-analyses^[35-37], a meta-analysis specific to the association between rs1137101 and DS cancer risk is lacking. Given that articles related to the relationship between rs1137101 and DS cancers are continually published and new experimental data were obtained in the present study, we performed a comprehensive meta-analysis to determine the contribution of rs1137101 to DS cancer risk.

This hospital-based case–control study included 1,727 patients with cancer and 800 healthy controls. Of these patients, 460 were diagnosed with gastric cancer, 787 with colorectal cancer, and 480 with liver cancer; the diagnoses were confirmed histopathologically. The patients were recruited from the Hubei Cancer Hospital and Wuhan Xinzhou District People’s Hospital and the healthy controls were selected from cancer-free individuals who visited the Wuhan Xinzhou District People’s Hospital for medical checkups. All enrolled participants were unrelated Chinese Han individuals living in Hubei Province. Information on age, sex, smoking status, and drinking status was collected from all participants. Regular smoking was defined as smoking at least one cigarette per day on average over one year or having quit smoking for less than one year. Regular drinking was defined as drinking at least 100 mL of alcohol per day on average over one year or having quit drinking for less than one year. This study was conducted after obtaining approval from the ethical committee of Wuhan University of Technology, and signed informed consent forms were obtained from all participants.

Peripheral blood samples from all participants 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. The DNA samples were stored at −20 °C before use. Genotyping of the LEPR gene rs1137101 variant was performed using a polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) assay. The DNA fragments containing rs1137101 were amplified using the primers 5'-TGTGCCAACAGCCAAACTCA-3' (forward) and 5'-CCCCAGGACTACATCTACCATCA-3' (reverse). The PCR products were digested with the restriction enzyme Msp I (Thermo Scientific, USA) for 4 h at 37 °C. The digested fragments were then separated on a 2.0% agarose gel and the RFLP bands were visualized under ultraviolet light using Gel-Red staining. The rs1137101 AA genotype was identified by the presence of a 243-bp band; the AG genotype by the presence of three bands of 243, 171, and 72 bp; and the GG genotype by the presence of two bands of 171 and 72 bp. To ensure the accuracy of the PCR-RFLP genotyping results, we performed Sanger sequencing on 20% of randomly selected DNA samples, and the results were 100% concordant. Supplementary Figure S1 (available in www.besjournal.com) shows the genotyping diagrams of rs1137101 obtained using PCR-RFLP and Sanger sequencing.

Figure S1.  Genotyping diagrams of the LEPR gene rs1137101 variant using Sanger sequencing (A) and PCR-RFLP assay (B). PCR-RFLP, polymerase chain reaction-restriction fragment length polymorphism.

All statistical analyses were conducted using SPSS 15.0 (Chicago, IL, USA). The chi-square (χ²) test was used to compare the distributions of age, sex, smoking status, and drinking status between the patients and healthy controls. The genotypic frequency of rs1137101 in healthy controls was examined for deviation from Hardy-Weinberg equilibrium (HWE). Unconditional logistic regression analysis under allele, carrier, homozygote, heterozygote, recessive, and dominant models was used to assess the correlation between rs1137101 and cancer risk. Statistical significance was set at P < 0.05, and an adjusted threshold of 0.0084 (0.05/6) by Bonferroni correction was used for multiple comparisons.

A thorough literature search of the PubMed, EMBASE, ISI, Web of Science, CNKI, and Wanfang databases up to August 31, 2023 was conducted to retrieve eligible studies related to rs1137101 and DS cancer risk. The following terms and keywords were searched: “leptin receptor” OR “LEPR” or “obese receptor” OR “OBR” AND “rs1137101” OR “Q223R” OR “p. Gln223Arg” or “c. 668A>G” and “cancer” or “tumor” OR “carcinoma” OR “malignancy” OR “neoplasms”. To identify additional eligible studies, we screened the references listed in the retrieved articles. The inclusion criteria were as follows: (1) case–control study, (2) rs1137101 and risk of DS cancer (oral, esophageal, gastric, liver, bile duct, colorectal, pancreatic, and gallbladder cancer), (3) human study, (4) extractable genotyping data, and (5) no departure from HWE in the controls. Supplementary Figure S2 (available in www.besjournal.com) shows a flowchart of the process for literature retrieval for this meta-analysis. Statistical analyses, including heterogeneity tests, pooled odds ratios (ORs), publication bias, and sensitivity analysis, were performed using STATA 14.0 (Stata Corp, College Station, Texas). P < 0.05 indicated the significance level and the Bonferroni correction was adopted for multiple comparisons. In addition, trial sequential analysis (TSA) was conducted to examine whether the statistical results were reliable and conclusive. The elucidation of TSA results was described in our previous study^[38].

Figure S2.  Flow diagram of the literature review process for this meta-analysis.

The principal characteristics of the participants are presented in Table 1. The distributions of age, sex, smoking status, and drinking status did not differ significantly between patients with cancer and healthy controls (P > 0.05). The average age of healthy controls and patients with gastric, colorectal, and liver cancer were 56.16 ± 10.85, 57.36 ± 10.46, 57.87 ± 10.76, and 57.57 ± 10.38 years, respectively. No statistical differences were observed in the average age between patients (with gastric, colorectal, or liver cancer) and healthy controls. These results indicate that the data from the patients with cancer and healthy controls can be compared for study purposes.

The association between rs1137101 and risk of gastric, colorectal, or liver cancer was examined in a Hubei Chinese population, and the statistical results are presented in Table 2. The observed genotypic distributions of rs1137101 were consistent with HWE in the healthy controls (P = 0.778), suggesting that this enrolled cohort was representative of the population. rs1137101 was nominally associated with gastric cancer risk at a significance level of 0.05 but the association did not remain after Bonferroni correction. No association was observed between the rs1137101 polymorphism and colorectal cancer risk. Of note, the presence of rs1137101 was significantly associated with liver cancer risk, even after Bonferroni correction. Specifically, compared to the G allele, the rs1137101 A allele conferred a lower risk of liver cancer (A vs. G, P = 0.002, OR = 0.66, 95% confidence intervals [CI]: 0.50–0.85). Individuals carrying the AG genotype and those carrying the AA + AG genotypes had reduced risk of liver cancer compared to individuals carrying the GG genotype (AG vs. GG, P = 0.006, OR = 0.67, 95% CI: 0.50–0.89; AA + AG vs. AG, P = 0.003, OR = 0.65, 95% CI: 0.49–0.86).

A total of 24 eligible studies were included in this meta-analysis. Two separate authors (HU Wei Qiong and ZHOU Wei Guang) carefully collected essential information from the eligible studies, which is presented in Supplementary Table S1 (available in www.besjournal.com). As shown in Table 3, no significant association was identified between rs1137101 and DS cancer risk. Interestingly, the stratified analysis by ethnicity also revealed lack of any significant association between rs1137101 and DS cancer risk in the Asian or Caucasian populations. In contrast, the stratified analysis by country suggested that rs1137101 was significantly associated with DS cancer risk in the Chinese population. Chinese individuals with the AA genotype were less likely to present with DS cancers than those with the AG, GG, and AG+GG genotypes (AA vs. AG, P < 0.001, OR = 0.62, 95% CI: 0.50–0.78; AA vs. GG, P < 0.001, OR = 0.53, 95% CI: 0.42–0.67; AA vs. AG + GG, P < 0.001, OR = 0.57, 95% CI: 0.46–0.70). Lastly, no significant association between rs1137101 and DS cancer risk was identified in the stratified analysis based on PCR-RFLP genotyping.

A stratified analysis according to cancer type was conducted, and the pooled results are presented in Table 4. We found that rs1137101 had no significant effect on the risk of oral, esophageal, or colorectal cancer. However, a marked effect of rs1137101 on the risk of gastric and liver cancers was observed in the total population. The AA genotype carriers were less likely to have gastric cancer than the GG (AA vs. GG, P = 0.006, OR = 0.44, 95% CI: 0.24–0.79) or AG+GG genotype carriers (AA vs. AG+GG, P = 0.008, OR = 0.45, 95% CI: 0.25–0.82). Moreover, the A allele showed a protective role compared to the G allele in liver carcinogenesis (A vs. G, P = 0.003, OR = 0.74, 95% CI: 0.61–0.91). The individuals with the AA genotype had a lower risk of presenting with liver cancer than those with the AG (AA vs. AG, P < 0.001, OR = 0.56, 95% CI: 0.42–0.75), GG (AA vs. GG, P < 0.001, OR = 0.41, 95% CI: 0.30–0.56), or AG+GG genotypes (AA vs. AG+GG, P < 0.001, OR = 0.50, 95% CI: 0.38–0.65). The individuals with the AG genotype were also less likely to have liver cancer than those with the GG genotype (AG vs. GG, P = 0.005, OR = 0.80, 95% CI: 0.69–0.94). Interestingly, the association between rs1137101 and liver cancer risk was confirmed in the Chinese population (AA vs. AG, AA vs. GG, AG vs. GG, and AA vs. AG+GG).

Under the genetic model of AA vs. GG, forest plots, funnel plots, sensitivity analysis, and TSA were separately performed for the meta-analysis on the association between rs1137101 and DS cancer risk in the total population (Figure 1A, 2A, 3A, and 4A) and Chinese population (Figures 1B, 2B, 3B, and 4B); and for the meta-analysis of the association of rs1137101 with gastric cancer risk in the total population (Figures 1C, 2C, 3C, and 4C) and with liver cancer risk in the total population (Figures 1D, 2D, 3D, and 4D). The Begg (Tables 3 and 4) and Egger tests (Tables 3 and 4) demonstrated no evidence of publication bias in any of the genetic models. The funnel plots of AA vs. GG also excluded the occurrence of publication bias. The sensitivity analysis showed that the pooled ORs did not significantly change after systematically removing individual studies. As shown in Figure 4A, the cumulative Z-curve did not cross Z = 1.96 or the trial sequential monitoring boundary (TSMB) before reaching the required information size line, providing firm evidence for the lack of association between rs1137101 and DS cancer risk in the total population. The TSA results in Figures 4B, 4C, and 4D show that the cumulative Z-curves crossed both Z = 1.96 and TSMB before reaching the required information size line, confirming the real contribution of rs1137101 to the risk of DS cancers in the Chinese population and risk of gastric and liver cancers in the total population.

Figure 1.  Forest plots of AA vs. GG for the meta-analysis of the association of rs1137101 with digestive-system cancer risk in the total population (A) and the Chinese population (B); and for the meta-analysis of the association of rs1137101 with gastric cancer risk (C) and liver cancer risk (D) in the total population. OR, odd ratio; CI, confidence interval.

Figure 4.  Trial sequential analysis of AA vs. GG for the meta-analysis on the relationship between rs1137101 and digestive-system cancer risk in the total population (A) and the Chinese population (B); and for the meta-analysis of the association of rs1137101 with gastric cancer risk (C) and liver cancer risk (D) in the total population.

DS cancers are a global health problem and a serious threat to human life. Detection at an early stage may be beneficial to patient outcome by improving the 5-year survival rate and may even lead to a complete cure. LEPR mediates key signaling pathways implicated in carcinogenesis. Specifically, the review by Rodrigues et al.^[7] shows that leptin binding to Ob-Rb in the cell membrane induces homodimerization of the receptor followed by JAK-2 phosphorylation and subsequent phosphorylation of Tyr1138 in Ob-Rb. These events lead to the activation of STAT3 via phosphorylation, thereby regulating the expression of several cancer-associated genes^[39]. Through STAT3, leptin can also trigger important cancer signaling pathways such as PI3K/mTOR. Additionally, through the phosphorylation of Tyr985 in Ob-Rb, leptin can activate Src homology 2 domain-containing phosphatase 2 (SHP2), triggering the MAPK/ERK pathway^[40]. This evidence strongly suggests that LEPR is a plausible candidate gene for cancer susceptibility^[10]. Consequently, the identification of certain causative variants in the LEPR gene would be of great help in the early diagnosis and treatment of DS cancers.

The LEPR gene rs1137101, a potential functional variant, has been extensively investigated for its contribution to the risk of DS cancers (21 studies up to August 31, 2023)^[13-31], breast cancer (24 studies up to February 10, 2022), and endometrial^[33] and lung cancers^[34]. The findings of the studies on rs1137101 and DS cancer risk are summarized in Supplementary Table S1. rs1137101 was significantly associated with oral cancer risk in the Caucasian but not in the Asian population, and with esophageal cancer risk in the Asian but not in the Caucasian population. A significant association between rs1137101 and gastric cancer risk has been previously identified in the Chinese population^[19-21] but the present study failed to confirm this finding. Lin et al. previously reported that the presence of rs1137101 may decrease colorectal cancer risk^[27], although no such association was noted in other studies^[22-26] or in the present study. A significant association between rs1137101 and liver cancer risk was observed in the Egyptian population^[30] and in Chinese populations from Henan and Hubei provinces^[28], but not in the Chinese population from Fujian or Guangxi provinces^[29,31]. These discrepancies may be attributed to the following reasons. First, a tissue‐specific pattern may exist for the regulatory network of rs1137101 in LEPR function in different types of DS cancers. Second, the interaction of rs1137101 with different genetic backgrounds may lead to different susceptibilities to cancer among different ethnic populations. Third, the recruited Chinese participants are from different geographical areas of China and differ in their living environments and lifestyles. Fourth, the insufficient sample size may have resulted in inadequate statistical strength.

Several meta-analyses have attempted to estimate the association between rs1137101 and cancer risk more accurately. A recent meta-analysis of 44 case–control studies by Rong et al. revealed no significant association between rs1137101 and overall cancer risk^[37]. Another recent meta-analysis by Zhu et al., which included 20 case-control studies, revealed a significant correlation between rs1137101 and breast cancer risk in African populations^[41]. Since the studies related to rs1137101 and breast cancer make up nearly half of all related studies, they likely affect the estimation of the association between rs1137101 and overall cancer risk. Similarly, the studies related to rs1137101 and DS cancers account for a high proportion of all related studies. Thus, meta-analyses are necessary to determine the effect of rs1137101 on DS cancer risk. After pooling the data from previous related studies and the present study, a significant association between rs1137101 and DS cancer risk was only identified in the Chinese population. In the stratified analysis by cancer type, rs1137101 was significantly associated with the risk of gastric and liver cancer but not oral, esophageal, or colorectal cancer.

Several limitations should be considered when interpreting these results. First, the recruitment of participants was hospital-based, which may potentially have contributed to a selection bias. Second, the association of rs1137101 with gastric cancer risk involved only a Chinese population; confirmatory research in other ethnic groups is needed. Third, considering the small number of included studies for each specific DS cancer in this meta-analysis, replication studies and updated meta-analyses are required for a more accurate estimation of the effect of rs1137101 on DS cancer risk. Lastly, no evidence is currently available to show that rs1137101 is functionally significant and contributes to cancer development in vitro and in vivo. Future functional studies are thus needed to gain a better understanding of the role that LEPR gene rs1137101 plays in cancer pathogenesis.

In conclusion, our case–control study suggests that the LEPR gene rs1137101 variant is associated with liver cancer risk in the Hubei Chinese population. The subsequent meta-analysis revealed that rs1137101 may contribute to susceptibility to DS cancers (especially gastric and liver cancer) in the Chinese population. Additional well-designed studies with larger sample sizes and involving different ethnic groups are required to confirm the present findings.