Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Genetic variability in LMP2 and LMP7 is associated with the risk of esophageal squamous cell carcinoma in the Kazakh population but is not associated with HPV infection

  • Lan Yang ,

    Contributed equally to this work with: Lan Yang, Yu Ji, Ling Chen

    Roles Conceptualization, Methodology, Writing – review & editing

    Affiliation Department of Pathology and Key Laboratory for Xinjiang Endemic and Ethnic Diseases, Shihezi University School of Medicine, Shihezi, Xinjiang, China

  • Yu Ji ,

    Contributed equally to this work with: Lan Yang, Yu Ji, Ling Chen

    Roles Conceptualization, Formal analysis, Investigation, Methodology, Visualization, Writing – original draft, Writing – review & editing

    Affiliation Department of Pathology and Key Laboratory for Xinjiang Endemic and Ethnic Diseases, Shihezi University School of Medicine, Shihezi, Xinjiang, China

  • Ling Chen ,

    Contributed equally to this work with: Lan Yang, Yu Ji, Ling Chen

    Roles Data curation

    Affiliation Department of Pathology and Key Laboratory for Xinjiang Endemic and Ethnic Diseases, Shihezi University School of Medicine, Shihezi, Xinjiang, China

  • Mei Li,

    Roles Investigation

    Affiliation Department of Pathology and Key Laboratory for Xinjiang Endemic and Ethnic Diseases, Shihezi University School of Medicine, Shihezi, Xinjiang, China

  • Fei Wu,

    Roles Investigation

    Affiliation Department of Pathology and Key Laboratory for Xinjiang Endemic and Ethnic Diseases, Shihezi University School of Medicine, Shihezi, Xinjiang, China

  • Jianming Hu,

    Roles Data curation

    Affiliation Department of Pathology and Key Laboratory for Xinjiang Endemic and Ethnic Diseases, Shihezi University School of Medicine, Shihezi, Xinjiang, China

  • Jinfang Jiang,

    Roles Data curation, Resources

    Affiliation Department of Pathology and Key Laboratory for Xinjiang Endemic and Ethnic Diseases, Shihezi University School of Medicine, Shihezi, Xinjiang, China

  • Xiaobin Cui,

    Roles Formal analysis, Software

    Affiliation Department of Pathology and Key Laboratory for Xinjiang Endemic and Ethnic Diseases, Shihezi University School of Medicine, Shihezi, Xinjiang, China

  • Yunzhao Chen,

    Roles Formal analysis

    Affiliation Department of Pathology and Key Laboratory for Xinjiang Endemic and Ethnic Diseases, Shihezi University School of Medicine, Shihezi, Xinjiang, China

  • Lijuan Pang,

    Roles Formal analysis

    Affiliation Department of Pathology and Key Laboratory for Xinjiang Endemic and Ethnic Diseases, Shihezi University School of Medicine, Shihezi, Xinjiang, China

  • Yutao Wei,

    Roles Validation

    Affiliation The First Affiliated Hospital, Shihezi University School of Medicine, Shihezi, Xinjiang, China

  • Feng Li

    Roles Conceptualization, Funding acquisition, Project administration, Supervision

    lifeng7855@126.com

    Affiliations Department of Pathology and Key Laboratory for Xinjiang Endemic and Ethnic Diseases, Shihezi University School of Medicine, Shihezi, Xinjiang, China, Department of Pathology, Beijing ChaoYang Hospital, Capital Medical University, Beijing, China

Abstract

The Kazakh population in Xinjiang Province in northwestern China exhibits a high incidence of esophageal squamous cell carcinoma (ESCC). Although the etiology of esophageal carcinoma (EC) has not been elucidated, there are reports of the involvement of an immunologic mechanism. In the current study, 268 Kazakh ESCC patients and 500 age- and sex-matched control subjects were recruited. DNA was extracted from paraffin-embedded tumor specimens from the patients and peripheral blood lymphocytes from the controls and used for LMP2/LMP7 genotyping. Polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) analysis was performed to detect LMP2/LMP7 gene single-nucleotide polymorphisms (SNPs). We found a clear increased risk of ESCC in the Kazakh population for the heterozygous LMP2 R/C genotype and the homozygous C/C genotype (OR = 1.470, 95%CI = 1.076–2.008, p = 0.015 forLMP2R/C; OR = 2.048, 95% CI = 1.168–3.591, p = 0.011 for LMP2 C/C). Conversely, the heterozygous LMP7 Q/K polymorphism was found to decrease the risk of ESCC in this population (OR = 0.421, 95% CI = 0.286–0.621, p = 8.83×10−6). Moreover, LMP2 R/C+C/C genotype was associated with increased tumor invasion depth (p = 0.041). Haplotype analysis showed that haplotype A, which includes wild-type homozygous LMP2/TAP1 and mutant LMP7, decreases susceptibility to ESCC in the Kazakh population; in contrast, haplotype E, which includes wild-type homozygous LMP2/LMP7/TAP1, acts as a risk factor for increased susceptibility to ESCC. This is the first study to report that the heterozygous LMP2 R/C and homozygous C/C genotypes increase susceptibility to ESCC in the Kazakh population and that the heterozygous LMP7 Q/K genotype decreases susceptibility to ESCC in this population. Nevertheless, neither LMP2 nor LMP7 was associated with human papillomavirus (HPV) infection. Understanding LMP2/LMP7 genetic variability will provide a new therapeutic perspective for Kazakh patients with ESCC.

Introduction

Esophageal carcinoma (EC) is considered one of the most malignant cancers worldwide [15]. More than 90% of ECs are either esophageal squamous cell carcinomas (ESCCs) or adenocarcinomas (EACs). It is reported that EAC is markedly more prevalent in western countries; in contrast, EC cases are almost entirely ESCC in eastern and African countries [2, 6, 7]. Some etiologic and risk factors contribute to this phenomenon. In western countries and areas with a low incidence of ESCC, excess alcohol intake in combination with smoking greatly increases the risk of EC [8]. In developing countries, diets deficient in vitamins as well as intake of carcinogens have been speculated as risk factors for ESCC [9]. The Kazakh population, a minority group in the Xinjiang region of northwest China that still maintains its traditional way of life (such as living in a dry, hot climate and consuming a vitamin-deficient diet), exhibits the highest ESCC mortality rate among different minorities in China. Indeed, the death rate in this population can be up to 68.95/100,000 [10]. Therefore, it is urgent to identify ways to reduce the mortality rate in this population. Recently, some researchers showed that endogenous factors, such as the immune system and genetic factors, may be related to ESCC. Moreover, many genetic variations, especially polymorphisms, are prevalent in the Kazakh population, some of which are associated with an increased risk of cancer [11].

The low-molecular-weight polypeptide (LMP) system recognizes intracellular infection and plays a crucial role in immunological surveillance in humans through the MHC-I molecule and cytotoxic T-lymphocyte (CTL) pathway [12, 13]. LMP genes are located in the MHC class II region of chromosome 6p21.3 [14]. The genes LMP2 and LMP7 encode two proteasome subunits that are considered to function in the degradation of cytosolic proteins and the production of antigenic peptides [15]. After antigenic peptides are degraded by LMP2 and LMP7, they are transferred into the endoplasmic reticulum (ER) lumen by the transporter associated with antigen processing (TAP) protein, loaded onto MHC class I molecules, and presented to cytotoxic T lymphocytes [16]. Thus, because LMPs participates in antigen presentation, they play vital roles in autoimmune diseases and tumorigenesis [1721]. A series of previous studies found that polymorphisms in LMP2 codon 60 (LMP2–60) (Arg–Cys) and LMP7 codon 145 (LMP7–145) (Gln–Lys) are responsible for functional evolution and thus reduce the ability of these molecules to function in antigen processing [22, 23]. Accordingly, these genetic variants in LMP2/LMP7 have been linked to the occurrence, development, and prognosis of many diseases, including viral infection, autoimmune disease, and malignant tumors [1720].

Many studies to date have found that LMP2/LMP7 polymorphism affects susceptibility to numerous cancers. For instance, the LMP7–145 gene polymorphism contributes to increased susceptibility to gastric cancer and ovarian cancer [11, 24]. LMP7–145 lysine carriers in the Han population in Anyang, Henan, China, were also found to be more susceptible to EC [22]. In contrast, no statistical correlation has yet been found between LMP2–60 polymorphic genotypes and the risk of ovarian or lung cancers [11, 25]. However, the relationship between LMP2/LMP7 polymorphism and the risk of ESCC in the Kazakh population remains unclear. Therefore, studying LMP2/LMP7 genetic variability in Kazakh patients with ESCC is important for providing a new therapeutic perspective.

Materials and methods

Specimens

Two hundred sixty-eight formalin-fixed and paraffin-embedded ESCC tissues were obtained from 268 Kazakh patients who had been treated with esophagectomy without any additional therapy, such as radiation and chemotherapy, during the period from January 1990 to March 2014 at People’s Hospital of Xinjiang Uygur Autonomous Region, First Affiliated Hospital of Shihezi University Medical College in Xinjiang and some hospitals in the Yili Kazakh Autonomous Region. All 268 of the cases originated from the ESCC-high-risk Kazakh population in Xinjiang. The age of the study group ranged from 23 years to 79 years, with a mean of 53.14±8.20 years. A standardized questionnaire was used to collect information on demographic characteristics (gender, age, area of birth and residence) and life-style variables, including consumption of hot tea, alcohol and tobacco(S3 Table),we defined drinkers as drinking more than 20g per day,greater than 5 times a week,drinking time is greater than 10 years. Smokers are defined as smoking more than 1 cigarettes one day, smoking time is more than 6 months. In order to ensure the reliability of the results of the survey and the data collected by the authenticity, we carried out the analysis of reliability and validity of the questionnaire before the formal investigation, and according to the analysis results, then screened second issues, adjusted the overall structure of the questionnaire.Pathological diagnoses of tumor-node-metastasis (TNM) stages were estimated for all cases using guidelines provided in Cancer Stage Manual 7th Edition (2009) issued by the American Joint Committee on Cancer/Union for International Cancer Control (AJCC/UICC). For controls, we randomly selected 500 Kazakh volunteers from the same geographic area who had no cancer history and were gender- and age-matched with the cases. Blood samples from these 500 controls were collected, and DNA was extracted for LMP gene genotyping. All participants in the study signed an informed consent form. First Affiliated Hospital of Shihezi University School of Medicine Institutional Ethics Review Board approved this research. We conducted the experiment in 2015.

Genotyping of LMP2/ LMP7 polymorphisms and haplotype construction

DNA was extracted from paraffin sections and blood samples using standard proteinase K digestion and a tissue DNA extraction kit following the manufacturer’s instructions (QIAGEN, Hilden, Germany). The two previously documented polymorphisms in the LMP2/LMP7 coding regions are listed in S1 Table. The polymerase chain reaction (PCR) primers used to assess the polymorphisms mentioned above are listed in S1 Table. PCR was performed in a total volume of 25 μl containing 1 μl DNA, 2.5 μl 10× PCR buffer (100 mM Tris, pH 7.8, 100 mM NaCl, 10 mM ethylenediaminetetraacetic acid (EDTA) and 0.5% sodium dodecyl sulfate (SDS)), 1.8 μl 25 mM MgCl2, 0.5 ml each dNTP, 10 pmol primers and 1.5 U Taq Polymerase (Promega, Madison, WI, USA). For PCR amplification, the reaction was initially maintained at 94°C for 2 min, followed by 30 cycles of the following conditions: denaturation at 94°C for 40 s, annealing at 57.5°C for 40 s, elongation at 72°C for 40 s. A final extension at 72°C was performed for 10 min. Gene-Amp PCR System 9700 was used. The PCR products were purified using a MultiScreen-PCR purifying plate (Millipore Company, USA).

For restriction fragment length polymorphism (RFLP) analysis, 330-bp and 351-bp PCR products for the LMP2–60 and LMP7–145 polymorphisms were digested using the specific restriction endonucleases HhaI and BsmI (New England Biolabs, Ipswich, MA, USA). The PCR products were digested under the conditions recommended by the manufacturer and then electrophoresed through a 3% agarose gel.

The wild-type genotype of the LMP2–60 polymorphism yielded bands of 79 and 251 bp, and the heterozygous genotype yielded bands of 79, 251 and 330 bp (S1 Fig). Regarding LMP7–145, bands of 146 and 205 bp for the Gln/Gln genotype, a band of 351 bp for the Lys/Lys genotype, and bands of 146, 205 and 351 bp for the Gln/Lys genotype were obtained (S2 Fig). Fragment sizes were evaluated by comparison with known markers (Fermentas, Inc, Glen Burnie, MD). DNA direct sequencing was used to confirm the genotypes of LMP2–60 and LMP7–145 (S1 Fig and S2 Fig), and the two results were 100% concordant. LMP2/LMP7 and TAP1 gene polymorphism haplotypes were constructed using PHASE [26] software.

Detection of HPV DNA

DNA was extracted from 268 paraffin-embedded tissue samples from Kazakh patients using standard proteinase K digestion and a tissue DNA extraction kit (QIAGEN, Hilden, Germany). We used a human β-actin primer set as an internal control for PCR: forward, 5’-CAGACACCATGGTGCACCTGAC-3’, and reverse, 5’-CCAATAGGCAGAGAGAGTCAGTG-3’. We used non-degenerate primers to determine HPV infection (GP 5+/6+, forward primer 5’-TTGGATCCTTTGTACTGTGGTAGATACTAC-3’, and reverse primer 5’-TTGGATCCGAAAAATAAACTGTAAATCATATTC-3’), which amplify a 150-bp fragment of the L1 gene in a large range of HPV types. For HPV16 detection, E7 was amplified with forward primer 5’-GATGAAATAGATGGTCCAGC-3’ and reverse primer 5’-GCTTTGTACGCACAACCGAGC-3’. For each PCR reaction, the final volume was 25 μl and included 5 μl extracted concentrated DNA. The reaction was first carried out at 95°C for 10 min, followed by 40 cycles at 94°C for 30 s, 42°C for 90 s, 72°C for 30 s, and a final extension at 72°C for 5 min. The amplicons were then denatured and subjected to hybridization. The samples were tested in triplicate, and β-actin was used as an internal control.

Statistical analysis

All statistical analyses were performed using SPSS version 17.0 (U.S.). Hardy-Weinberg equilibrium (HWE) for the cases and controls was evaluated using chi-square tests to compare the observed and expected genotype frequencies (S2 Table). Allele and genotype frequencies were compared by contingency table analysis using Fisher’s exact test if the number of observations in any cell was less than five. The p values were corrected based on the number of specimens tested and the number of comparisons performed. Student’s t-test and the χ2 test were utilized to calculate the significance of differences in continuous variables between LMP genotypes. The relationship between single-nucleotide polymorphisms (SNPs) in LMP2/LMP7 and the risk of ESCC was estimated by determining the odds ratios (ORs) at 95% confidence intervals (CIs) using logistic regression analyses for crude ORs and adjusted ORs when adjusting for age and sex. Logistic regression was also applied to analyze the HPV status and risk of ESCC.

Results

Characteristics of the study population

The characteristics of the 268 cases and 500 controls included in this study are summarized in Table 1. The age was 53.14 ± 8.20 (mean ± SD) years for the ESCC cases and 54.59±10.75 (mean ± SD) years for the controls. No significant difference was observed for age and gender (p = 0.095 and p = 0.639, respectively). Statistical analysis showed that smoking is closely linked to the risk of ESCC in the Kazakh population (p = 0.027), though drinking was not significantly associated (p = 0.393). We also evaluated data for the case group with regard to differentiation status, tumor location, depth of invasion and lymph node metastasis.

Genotype frequencies of LMP2/LMP7 gene polymorphisms in cases and controls

The genotype frequencies of LMP2/LMP7 gene polymorphisms are presented in Table 2. In single-locus analyses of genotype frequencies, the frequencies for LMP2 genotypes were 41.4% (R/R), 48.5% (R/C), and 10.1% (C/C) among the cases and 52.2% (R/R), 41.6% (R/C), and 6.2% (C/C) among the controls. A significant difference was observed in this study. Relative to the R/R allele, R/C and C/C genotypes were associated with increased risk of ESCC in the Kazakh population (OR = 1.470, 95%CI = 1.076–2.008, p = 0.015 for R/C, OR = 2.048, 95%CI = 1.168–3.591, p = 0.011 for C/C, Table 2). The C allele also showed an increased risk of ESCC in the Kazakh population (OR = 1.413, 95% CI = 1.127–1.772, p = 0.003). In the dominant model, the RC+CC genotype was associated with a statistically increased risk of ESCC in the Kazakh population when the RR genotype was used as a reference group (RR versus RC+CC OR = 1.545, 95% CI = 1.145–2.085, p = 0.004, Table 3). However, in the recessive model, in which the RR+RC genotype was used as the reference group, the CC genotype was not associated with risk of ESCC in the Kazakh population (RR+RC versus CC OR = 1.695, 95% CI = 0.989–2.905, p = 0.053, Table 3).

thumbnail
Table 2. Genotype frequencies of LMP2/LMP7 gene polymorphisms in cases and controls.

https://doi.org/10.1371/journal.pone.0186319.t002

thumbnail
Table 3. Stratification analyses between LMP2 polymorphism and Kazakh ESCC patient clinicopathological parameters.

https://doi.org/10.1371/journal.pone.0186319.t003

A significant difference was also observed for LMP7, with genotype frequencies of 81.3% (Q/Q), 14.9% (Q/K), and 3.73% (K/K) among the cases and 68.4% (Q/Q), 29.8% (Q/K), and 1.8% (K/K) among the controls. Compared with the Q/Q genotype, the Q/K genotype decreased the risk of ESCC in the Kazakh population (OR = 0.421, 95%CI = 0.286–0.621, p<0.001,Table 2). The K allele also decreased the risk of ESCC (OR = 0.633, 95% CI = 0.462–0.869, p = 0.004). The dominant model, for which the QQ genotype was used as a reference group, showed that the QK+KK genotype was associated with a statistically increased risk of ESCC in the Kazakh population (QQ versus QK+KK (OR = 0.496, 95%CI = 0.346–0.172, p<0.001), Table 3). In contrast, the KK genotype was not associated with the risk of ESCC in the Kazakh population in the recessive model (QQ+QK versus KK (OR = 1.684, 95% CI = 0.675–4.201, p = 0.259), Table 2). We used the expression of quantitative trait loci (eQTL) data which have been reported (https://gtexportal.org/home/) to study the effects of SNP locus of LMP2/LMP7 gene on gene expression, the results showed that SNP locus of LMP2 in esophageal tissues affected the expression of LMP2 (p = 1.30×10–7, S3 Fig). The SNP locus of LMP7 was not related to the expression of LMP7.This phenomenon may be due to a same sense mutation in the SNP locus of LMP7. Although the base substitution have been replaced, the protein levels have not been altered, and amino acids have not been replaced.

Correlations between clinicopathological parameters and LMP2/LMP7 polymorphism in Kazakh patients with ESCC

To evaluate the contribution of confounding factors, including gender /age/ smoking and drinking to the risk of ESCC in the Kazakh population, stratification analyses were performed to investigate the potential effect of genetic variants of LMP2/LMP7 with ESCC risk in population subgroups (Table 3 and Table 4). This analysis showed that the LMP2–60 R/C and R/C+C/C genotypes were associated with an increased risk of ESCC in Kazakh patients aged ≥57 years but not in those aged <57 years. A significantly elevated risk was also identified in male subjects. In contrast, the LMP7–145 Q/K and QK+KK genotypes were associated with a decreased risk of ESCC in male patients aged ≥57 years but not in male patients aged <57 years. This analysis explained that the LMP2–60 R/C and R/C+C/C genotypes were associated with an increased risk of ESCC in Kazakh patients with smoking but not in those without smoking(p<0.05). There was no statistically significant with drinker and nondrinker(p>0.05). Equally contrary, the LMP7–145 Q/K and QK+KK genotypes were associated with a decreased risk of ESCC with smoking people but not in those without smoking(p<0.05). Nevertheless, the result indicated that LMP7–145 Q/K and QK+KK genotypes were associated with a decreased risk of ESCC who don't drink but not in those with drinking(p<0.05). When Kazakh ESCC patients were divided into two subgroups, T1/T2 and T3/T4, according to the AJCC TNM classification of EC, the RC+CC genotype was significantly associated with ESCC tumor depth (OR = 1.679, 95% CI = 1.023–2.756, p = 0.041). However, there was no significant association between LMP2/LMP7 polymorphism and ESCC in the Kazakh population with respect to clinical pathological parameters (p>0.05).

thumbnail
Table 4. Stratification analyses between LMP7 polymorphism and clinicopathological parameters of Kazakh ESCC patients.

https://doi.org/10.1371/journal.pone.0186319.t004

Linkage disequilibrium and haplotype analysis of LMP2 /LMP7 and TAP1

The association between TAP1 polymorphisms and ESCC in the Kazakh population was shown in a previous study (40). Thus, haplotype analysis was performed to further evaluate the combined effects of LMP2/LMP7 and TAP1 polymorphisms on the risk of ESCC in this population. Among the eight haplotypes derived from the observed genotypes, A and E were the most common in both the cases and controls (A: 33.77% of patients and 45.5% of controls, E: 28.73% of patients and 17.6% of controls) (Table 5). The results showed that haplotype A, which includes wild-type homozygous LMP2/TAP1 and mutant LMP7, decreases susceptibility to ESCC in the Kazakh population (OR = 0.37, 95%CI = 0.23–0.59, p<0.001), whereas haplotype E, which includes wild-type homozygous LMP2/LMP7/TAP1, acts as a risk factor to increase susceptibility to ESCC in the Kazakh population (OR = 3.67, 95% CI = 1.85–7.28, p<0.001). The results are summarized in Table 6.

thumbnail
Table 5. Distribution of estimated haplotype frequencies for LMP/TAP genes in cases and controls.

https://doi.org/10.1371/journal.pone.0186319.t005

thumbnail
Table 6. Association of haplotypes with polymorphisms of LMP2, LMP7 and TAP1-2 with risk of ESCC in the Kazakh population.

https://doi.org/10.1371/journal.pone.0186319.t006

Relationship between HPV infection in Kazakh patients with ESCC and LMP2/LMP7 genotypes

The results of genotype frequency analysis for LMP2/LMP7 gene polymorphisms in HPV-positive and HPV-negative patients are summarized in Table 3 and Table 4. There was no association between LMP2/LMP7 and HPV infection in Kazakh patients with ESCC (p>0.05).

Discussion

ESCC is a polygenic hereditary disease that can be affected by genetic factors [2730]. In our study, the association between LMP2 and LMP7 polymorphisms and the risk of ESCC in the Kazakh population was investigated.

According to our results, heterozygous LMP2 R/C and homozygous LMP2 C/C genotypes are risk factors for ESCC in the Kazakh population. This finding differed from that of Xiang Ma[24] and Bangwei Cao’s [22], which revealed no statistical correlation between LMP2–60 polymorphism and ovarian and gastric cancers.

However, one study did find that LMP2 plays a role in the growth of multiple myeloma and acute myeloid leukemia but that LMP7 is not associated with these hematological malignancies [31]. Nevertheless, in our study, the heterozygous LMP7 Q/K genotype was associated with significantly decreased risk of ESCC in the Kazakh population. This finding was also not consistent with previous results showing that LMP7–145 increased susceptibility to gastric cancer, ovarian cancer and ESCC in a Chinese Han population [11, 22, 24]. The discrepancy between these results and our findings may be attributed to the fact that the same polymorphism may have diverse genetic effects on different types of cancer, ethnicities and races [32, 33]. When considered together, Kazakh patients possessing both the LMP2 R/C+C/C genotype and the wild-type homozygous LMP7 Q/Q genotype exhibited a higher risk of developing ESCC, further illustrating the existence of linkage disequilibrium in LMP polymorphisms.

In subgroup analysis, the LMP2 R/C genotype was associated with increased tumor invasion depth in ESCC, suggesting that genetic variation in immune response genes may play a role in ESCC progression. Additionally, no association between LMP7–145 polymorphism genotypes and ESCC clinicopathological characteristics of the Kazakh population was observed in our study. Song [11] revealed that LMP7 polymorphism increases the risk of lymph node and tumor distant metastasis in ovarian cancer.

The LMP2–60 C/C and R/C+C/C genotypes were associated with an increased risk of ESCC in Kazakh patients aged ≥57 years but not in those aged <57 years. This indicates that carcinogenesis is an accumulation of genetic events during aging [34]. Additionally, the increased risk of ESCC in the elderly also demonstrates that gene–environment interactions may be involved in carcinogenesis. A significantly elevated risk was identified in male subjects, consistent with the results of some studies suggesting that female sex hormone exposure may play a protective role in the development of this cancer. Nonetheless, the LMP7–145 Q/K and QK+KK genotypes were associated with a decreased risk of ESCC in male patients aged ≥57 years but not in male patients aged <57 years.

Our previous studies found that TAP1 polymorphisms could increase susceptibility to ESCC in the Kazakh population [35]. Moreover, TAP1, LMP2 and LMP7 are all located within a narrow region of the class II cluster of the major histocompatibility complex on chromosome 6 [14]. In our study, eight haplotypes were constructed to analyze linkage disequilibrium for LMP2, LMP7 and TAP1, and our results indicated that haplotypes A and E were more common in both the patients and control individuals. Haplotype A, including wild-type homozygous LMP2/TAP1 and mutant LMP7, decreased the risk of ESCC in the Kazakh population. Moreover, haplotype E, including wild-type homozygous LMP2/LMP7/TAP1, increased the risk of ESCC in this population. These results indicate that the LMP7–145 (Gln–Lys) gene variation may weaken the susceptibility of the haplotype to ESCC in the Kazakh population. However, further studies are needed to clarify this association.

To date, many studies have confirmed that HPV is present in ESCC samples, with infection rates ranging from 17.1 [36] to 78.11% [22]. One study reported an infection rate of 100%, which was confirmed by detection of the HPV16 E6 and E7 genes in a high-risk area in Anyang, Henan [37]. Some reports even found that HPV infection is an important factor contributing to the high incidence of ESCC in the Kazakh population in Xinjiang [35]. Indeed, a study on cervical cancer showed that the genetic makeup of host immune response genes may be responsible for biological variability to HPV infection [38]. Our study sought to assess a connection between LMP2/LMP7 polymorphism and HPV infection, though we found no association between these variables. Our results indicate that HPV infection and LMP2/LMP7 polymorphisms may function independently to influence the risk of ESCC in the Kazakh population, a conclusion that supports the research of Cao [22] in Anyang, China.

In conclusion, to the best of our knowledge, the present study is the first to demonstrate that LMP2/LMP7 polymorphisms are associated with the risk of ESCC in the Kazakh population, a Chinese ethnic minority. LMP2 can be considered a risk factor for ESCC and LMP7 as a protective factor against ESCC in the Kazakh population. Additional studies with larger sample sizes and the inclusion of various populations, as well as functional studies, are required to verify the present preliminary findings.

Supporting information

S1 Fig. LMP2 gene enzyme digestion and sequencing map.

https://doi.org/10.1371/journal.pone.0186319.s001

(PDF)

S2 Fig. LMP7 gene enzyme digestion and sequencing map.

https://doi.org/10.1371/journal.pone.0186319.s002

(PDF)

S3 Fig. Bioinformatics tools to analyze LMP2 gene expression.

https://doi.org/10.1371/journal.pone.0186319.s003

(PDF)

S1 Table. Primer sets used for amplification and sequencing of LMP/TAP genes.

https://doi.org/10.1371/journal.pone.0186319.s004

(PDF)

S2 Table. Hardy-Weinberg test for cases and controls.

https://doi.org/10.1371/journal.pone.0186319.s005

(PDF)

S3 Table. Epidemiological questionnaire in Chinese and in English.

https://doi.org/10.1371/journal.pone.0186319.s006

(PDF)

References

  1. 1. Parkin DM, Bray F, Ferlay J, Pisani P. Global cancer statistics, 2002. CA: a cancer journal for clinicians. 2005;55(2):74–108. Epub 2005/03/12. pmid:15761078.
  2. 2. Lambert R, Hainaut P. The multidisciplinary management of gastrointestinal cancer. Epidemiology of oesophagogastric cancer. Best practice & research Clinical gastroenterology. 2007;21(6):921–45. Epub 2007/12/12. doi: pmid:18070696.
  3. 3. Kamangar F, Dores GM, Anderson WF. Patterns of cancer incidence, mortality, and prevalence across five continents: defining priorities to reduce cancer disparities in different geographic regions of the world. Journal of clinical oncology: official journal of the American Society of Clinical Oncology. 2006;24(14):2137–50. Epub 2006/05/10. doi: pmid:16682732.
  4. 4. Zheng S, Vuitton L, Sheyhidin I, Vuitton DA, Zhang Y, Lu X. Northwestern China: a place to learn more on oesophageal cancer. Part one: behavioural and environmental risk factors. European journal of gastroenterology & hepatology. 2010;22(8):917–25. Epub 2010/06/04. pmid:20520561.
  5. 5. Morita M, Kumashiro R, Kubo N, Nakashima Y, Yoshida R, Yoshinaga K, et al. Alcohol drinking, cigarette smoking, and the development of squamous cell carcinoma of the esophagus: epidemiology, clinical findings, and prevention. International journal of clinical oncology. 2010;15(2):126–34. Epub 2010/03/13. pmid:20224884.
  6. 6. Messmann H. Squamous cell cancer of the oesophagus. Best practice & research Clinical gastroenterology. 2001;15(2):249–65. Epub 2001/05/18. pmid:11355914.
  7. 7. Kamangar F, Chow WH, Abnet CC, Dawsey SM. Environmental causes of esophageal cancer. Gastroenterology clinics of North America. 2009;38(1):27–57, vii. Epub 2009/03/31. pmid:19327566; PubMed Central PMCID: PMCPmc2685172.
  8. 8. Enzinger PC, Mayer RJ. Esophageal cancer. The New England journal of medicine. 2003;349(23):2241–52. Epub 2003/12/06. pmid:14657432.
  9. 9. Sur M, Cooper K. The role of the human papilloma virus in esophageal cancer. Pathology. 1998;30(4):348–54. Epub 1998/12/05. pmid:9839308.
  10. 10. Chang F, Syrjanen S, Wang L, Syrjanen K. Infectious agents in the etiology of esophageal cancer. Gastroenterology. 1992;103(4):1336–48. Epub 1992/10/01. pmid:1327935.
  11. 11. Song L, Ma N, Han L, Yan H, Yan B, Yuan Z, et al. Association between LMP2/LMP7 genetic variability and the metastasis risk of ovarian cancer in Chinese women in Beijing. Human immunology. 2014;75(3):239–44. Epub 2014/01/01. pmid:24374040.
  12. 12. Grotzke JE, Lewinsohn DM. Role of CD8+ T lymphocytes in control of Mycobacterium tuberculosis infection. Microbes and infection / Institut Pasteur. 2005;7(4):776–88. Epub 2005/04/13. pmid:15823514.
  13. 13. Rafiei A, Kiutade Y. Identification of Mycobacterium tuberculosis CTL epitopes restricted by HLA-A*0201 in HHD mice. Iranian biomedical journal. 2007;11(1):23–31. Epub 2007/12/07. pmid:18051701.
  14. 14. Monaco JJ. A molecular model of MHC class-I-restricted antigen processing. Immunology today. 1992;13(5):173–9. Epub 1992/05/01. pmid:1386516.
  15. 15. Schutt P, Schutt B, Switala M, Bauer S, Stamatis G, Opalka B, et al. Prognostic relevance of soluble human leukocyte antigen-G and total human leukocyte antigen class I molecules in lung cancer patients. Human immunology. 2010;71(5):489–95. Epub 2010/02/17. pmid:20156510.
  16. 16. Ceppellini R, Frumento G, Ferrara GB, Tosi R, Chersi A, Pernis B. Binding of labelled influenza matrix peptide to HLA DR in living B lymphoid cells. Nature. 1989;339(6223):392–4. Epub 1989/06/01. pmid:2786148.
  17. 17. Camarena A, Aquino-Galvez A, Falfan-Valencia R, Sanchez G, Montano M, Ramos C, et al. PSMB8 (LMP7) but not PSMB9 (LMP2) gene polymorphisms are associated to pigeon breeder's hypersensitivity pneumonitis. Respiratory medicine. 2010;104(6):889–94. Epub 2010/02/16. pmid:20153157.
  18. 18. Casp CB, She JX, McCormack WT. Genes of the LMP/TAP cluster are associated with the human autoimmune disease vitiligo. Genes and immunity. 2003;4(7):492–9. Epub 2003/10/11. pmid:14551602.
  19. 19. Vargas-Alarcon G, Gamboa R, Zuniga J, Fragoso JM, Hernandez-Pacheco G, Londono J, et al. Association study of LMP gene polymorphisms in Mexican patients with spondyloarthritis. Human immunology. 2004;65(12):1437–42. Epub 2004/12/18. pmid:15603870.
  20. 20. Deshpande A, Wheeler CM, Hunt WC, Peyton CL, White PS, Valdez YE, et al. Variation in HLA class I antigen-processing genes and susceptibility to human papillomavirus type 16-associated cervical cancer. The Journal of infectious diseases. 2008;197(3):371–81. Epub 2008/02/06. pmid:18248301.
  21. 21. Kramer U, Illig T, Grune T, Krutmann J, Esser C. Strong associations of psoriasis with antigen processing LMP and transport genes TAP differ by gender and phenotype. Genes and immunity. 2007;8(6):513–7. Epub 2007/06/22. pmid:17581627.
  22. 22. Cao B, Tian X, Li Y, Jiang P, Ning T, Xing H, et al. LMP7/TAP2 gene polymorphisms and HPV infection in esophageal carcinoma patients from a high incidence area in China. Carcinogenesis. 2005;26(7):1280–4. Epub 2005/03/19. pmid:15774487.
  23. 23. Xu C, Qi S, Gao L, Cui H, Liu M, Yang H, et al. Genetic polymorphisms of LMP/TAP gene and hepatitis B virus infection risk in the Chinese population. Journal of clinical immunology. 2007;27(5):534–41. Epub 2007/05/26. pmid:17525827.
  24. 24. Ma X, Yang C, Tang R, Xu Z, Zhang Z, Wang Y, et al. Association between LMP2 and LMP7 gene polymorphisms and the risk of gastric cancer: A case-control study. Oncology letters. 2015;10(1):509–17. Epub 2015/07/15. pmid:26171060; PubMed Central PMCID: PMCPmc4487101.
  25. 25. Cao B, Dong S, Su Q, Yan B, Zhen H, Yuan Z, et al. Lack of association between the unique LMP2 gene polymorphism and the outcome of lung cancer in a population of Chinese Han nationality. Human immunology. 2012;73(5):580–4. Epub 2012/01/21. pmid:22261389.
  26. 26. Stephens M, Smith NJ, Donnelly P. A new statistical method for haplotype reconstruction from population data. American journal of human genetics. 2001;68(4):978–89. Epub 2001/03/20. pmid:11254454; PubMed Central PMCID: PMCPmc1275651.
  27. 27. Wang Y, Vogelsang M, Schafer G, Matejcic M, Parker MI. MicroRNA polymorphisms and environmental smoke exposure as risk factors for oesophageal squamous cell carcinoma. PloS one. 2013;8(10):e78520. Epub 2013/11/10. pmid:24205249; PubMed Central PMCID: PMCPmc3804537.
  28. 28. Cui XB, Chen YZ, Pang XL, Liu W, Hu JM, Li SG, et al. Multiple polymorphisms within the PLCE1 are associated with esophageal cancer via promoting the gene expression in a Chinese Kazakh population. Gene. 2013;530(2):315–22. Epub 2013/08/29. pmid:23981775.
  29. 29. Cui XB, Pang XL, Li S, Jin J, Hu JM, Yang L, et al. Elevated expression patterns and tight correlation of the PLCE1 and NF-kappaB signaling in Kazakh patients with esophageal carcinoma. Medical oncology. 2014;31(1):791. Epub 2013/12/07. pmid:24307345.
  30. 30. Qin JM, Yang L, Chen B, Wang XM, Li F, Liao PH, et al. Interaction of methylenetetrahydrofolate reductase C677T, cytochrome P4502E1 polymorphism and environment factors in esophageal cancer in Kazakh population. World journal of gastroenterology. 2008;14(45):6986–92. Epub 2008/12/06. pmid:19058336; PubMed Central PMCID: PMCPmc2773864.
  31. 31. Ozbas-Gerceker F, Bozman N, Kok S, Pehlivan M, Yilmaz M, Pehlivan S, et al. Association of an LMP2 polymorphism with acute myeloid leukemia and multiple myeloma. Asian Pacific journal of cancer prevention: APJCP. 2013;14(11):6399–402. Epub 2014/01/01. pmid:24377540.
  32. 32. Yu Z, Liu Q, Huang C, Wu M, Li G. The interleukin 10 -819C/T polymorphism and cancer risk: a HuGE review and meta-analysis of 73 studies including 15,942 cases and 22,336 controls. Omics: a journal of integrative biology. 2013;17(4):200–14. Epub 2013/04/12. pmid:23574339; PubMed Central PMCID: PMCPmc3615694.
  33. 33. Zhang YM, Zhou XC, Xu Z, Tang CJ. Meta-analysis of epidemiological studies of association of two polymorphisms in the interleukin-10 gene promoter and colorectal cancer risk. Genetics and molecular research: GMR. 2012;11(3):3389–97. Epub 2012/10/20. pmid:23079832.
  34. 34. Milne AN, Carvalho R, Morsink FM, Musler AR, de Leng WW, Ristimaki A, et al. Early-onset gastric cancers have a different molecular expression profile than conventional gastric cancers. Modern pathology: an official journal of the United States and Canadian Academy of Pathology, Inc. 2006;19(4):564–72. Epub 2006/02/14. pmid:16474375.
  35. 35. Zou N, Yang L, Chen L, Li T, Jin T, Peng H, et al. Heterozygote of TAP1 Codon637 decreases susceptibility to HPV infection but increases susceptibility to esophageal cancer among the Kazakh populations. Journal of experimental & clinical cancer research: CR. 2015;34:70. Epub 2015/07/25. pmid:26205887; PubMed Central PMCID: PMCPmc4514451.
  36. 36. de Villiers EM, Lavergne D, Chang F, Syrjanen K, Tosi P, Cintorino M, et al. An interlaboratory study to determine the presence of human papillomavirus DNA in esophageal carcinoma from China. International journal of cancer. 1999;81(2):225–8. Epub 1999/04/03. pmid:10188723.
  37. 37. Li T, Lu ZM, Chen KN, Guo M, Xing HP, Mei Q, et al. Human papillomavirus type 16 is an important infectious factor in the high incidence of esophageal cancer in Anyang area of China. Carcinogenesis. 2001;22(6):929–34. Epub 2001/05/29. pmid:11375901.
  38. 38. Gostout BS, Poland GA, Calhoun ES, Sohni YR, Giuntoli RL 2nd, McGovern RM, et al. TAP1, TAP2, and HLA-DR2 alleles are predictors of cervical cancer risk. Gynecologic oncology. 2003;88(3):326–32. Epub 2003/03/22. pmid:12648582.