Variations in the MHC Region Confer Risk to Esophageal Squamous Cell Carcinoma on the Subjects from High-Incidence Area in Northern China

Abstract

Background

The human major histocompatibility complex (MHC) is the most important region in vertebrate genome, and is crucial in innate immunity. Recent studies have demonstrated the possible role of polymorphisms in the MHC region to high risk for esophageal squamous cell carcinoma (ESCC). Our previous genome-wide association study (GWAS) has indicated that the MHC region may confer important risk loci for ESCC, but without further fine mapping. The aim of this study is to further identify the risk loci in the MHC region for ESCC in Chinese population.

Methods

Conditional logistic regression analysis (CLRA) was performed on 24 single nucleotide polymorphisms (SNPs) within the MHC region, which were obtained from the genetically matched 937 cases and 692 controls of Chinese Han population. The identified promising SNPs were further correlated with clinical and clinicopathology characteristics. Immunohistochemistry was performed to explore the protein expression pattern of the related genes in ESCC and neighboring normal tissues.

Results

Of the 24 promising SNPs analyzed, we identified three independent SNPs in the MHC region associated with ESCC: rs35399661 (P = 6.07E-06, OR = 1.71, 95%CI = 1.36–2.17), rs3763338 (P = 1.62E-05, OR = 0.63, 95%CI = 0.50–0.78) and rs2844695 (P = 7.60E-05, OR = 0.74, 95%CI = 0.64–0.86). These three SNPs were located at the genes of HLA-DQA1, TRIM27, and DPCR1, respectively. Further analyses showed that rs2844695 was preferentially associated with younger ESCC cases (P = 0.009). The positive immunostaining rates both for HLA-DQA1 and TRIM27 were much higher in ESCC tissues than in neighboring normal tissues (69.4% vs. 26.8% for HLA-DQA1 and 77.6% vs. 47.8% for TRIM27, P<0.001). Furthermore, the overexpression of HLA-DQA1 is correlated significantly with age (P = 0.001) and family history (P<0.001).

Conclusion

This study for the first time provides evidence that multiple genetic factors within the MHC region confer risk to ESCC on the subjects from high-risk area in northern China.

About the Authors:

Fang-Fang Shen

Contributed equally to this work with: Fang-Fang Shen, Wen-Bin Yue, Fu-You Zhou

Affiliations Henan Key Laboratory for Esophageal Cancer Research of the First Affiliated Hospital, Zhengzhou University, Zhengzhou, Henan, China, Cancer Research Center and the Third Affiliated Hospital, Xinxiang Medical University, Xinxiang, Henan, China

Wen-Bin Yue

Contributed equally to this work with: Fang-Fang Shen, Wen-Bin Yue, Fu-You Zhou

Affiliation: Henan Key Laboratory for Esophageal Cancer Research of the First Affiliated Hospital, Zhengzhou University, Zhengzhou, Henan, China

Fu-You Zhou

Contributed equally to this work with: Fang-Fang Shen, Wen-Bin Yue, Fu-You Zhou

Affiliation: Department of Thoracic Surgery, Anyang Tumor Hospital, Anyang, Henan, China

Ying Pan

Affiliation: Cancer Research Center and the Third Affiliated Hospital, Xinxiang Medical University, Xinxiang, Henan, China

Xue-Ke Zhao

Affiliation: Henan Key Laboratory for Esophageal Cancer Research of the First Affiliated Hospital, Zhengzhou University, Zhengzhou, Henan, China

Yan Jin

Affiliation: Henan Key Laboratory for Esophageal Cancer Research of the First Affiliated Hospital, Zhengzhou University, Zhengzhou, Henan, China

Xin Song

Affiliation: Henan Key Laboratory for Esophageal Cancer Research of the First Affiliated Hospital, Zhengzhou University, Zhengzhou, Henan, China

Bei Li

Affiliation: Henan Key Laboratory for Esophageal Cancer Research of the First Affiliated Hospital, Zhengzhou University, Zhengzhou, Henan, China

Xue-Na Han

Affiliation: Henan Key Laboratory for Esophageal Cancer Research of the First Affiliated Hospital, Zhengzhou University, Zhengzhou, Henan, China

Sa Tang

Affiliation: Henan Key Laboratory for Esophageal Cancer Research of the First Affiliated Hospital, Zhengzhou University, Zhengzhou, Henan, China

Yan Li

Affiliation: Henan Key Laboratory for Esophageal Cancer Research of the First Affiliated Hospital, Zhengzhou University, Zhengzhou, Henan, China

Guo Yuan

Affiliations Henan Key Laboratory for Esophageal Cancer Research of the First Affiliated Hospital, Zhengzhou University, Zhengzhou, Henan, China, Cancer Research Center and the Third Affiliated Hospital, Xinxiang Medical University, Xinxiang, Henan, China

Li-Sha Chen

Affiliations Henan Key Laboratory for Esophageal Cancer Research of the First Affiliated Hospital, Zhengzhou University, Zhengzhou, Henan, China, Cancer Research Center and the Third Affiliated Hospital, Xinxiang Medical University, Xinxiang, Henan, China

Ya-Li Liu

Affiliations Henan Key Laboratory for Esophageal Cancer Research of the First Affiliated Hospital, Zhengzhou University, Zhengzhou, Henan, China, Cancer Research Center and the Third Affiliated Hospital, Xinxiang Medical University, Xinxiang, Henan, China

Yan-Long Hu

Affiliations Henan Key Laboratory for Esophageal Cancer Research of the First Affiliated Hospital, Zhengzhou University, Zhengzhou, Henan, China, Cancer Research Center and the Third Affiliated Hospital, Xinxiang Medical University, Xinxiang, Henan, China

Xiu-Min Li

Affiliation: Cancer Research Center and the Third Affiliated Hospital, Xinxiang Medical University, Xinxiang, Henan, China

Jing-Li Ren

Affiliation: Department of Pathology, the Second Affiliated Hospital of Zhengzhou University, Zhengzhou, Henan, China

Li-Dong Wang

* E-mail: [email protected]

Affiliation: Henan Key Laboratory for Esophageal Cancer Research of the First Affiliated Hospital, Zhengzhou University, Zhengzhou, Henan, China

Introduction

Esophageal cancer (EC) is one of the most aggressive gastrointestinal malignancies and prevalent in the developing world [1]. The rural regions around the Taihang Mountain in the junction of Henan, Hebei, and Shanxi provinces in northern China have been well documented as the highest incidence areas for EC [2]-[4] _ENREF_1_ENREF_1. There are two main histological types for EC, squamous cell carcinoma (ESCC) and adenocarcinoma, each with distinct etiological and pathological characteristics [5]_ENREF_4. ESCC is the predominant histological type worldwide,especially in northern China, comprising more than 90% of all EC cases [6]_ENREF_5. The pathogenesis of ESCC is complex, involving both environmental and genetic risk factors. Recent genome-wide association studies (GWAS) have identified couple of susceptibility loci for ESCC in Han Chinese [7]–[9]. Assessing the individual genetic susceptibility can help identify high risk patients with a potential better benefit from the surveillance programs [10]. Our previous GWAS [9] has indicated that major histocompatibility complex (MHC) regions may confer important risk loci for ESCC, but without further fine mapping.

The MHC, human leukocyte antigen (HLA) system, is the most important region in vertebrate genome with respect to autoimmunity, and is crucial in innate immunity [11] _ENREF_9. MHC genes are located in a gene-dense region of the human genome on Chromosome 6p21.31 and the complex is organized from centromere to telomere as classes I, II, and III. Lack of HLA presentation has been proposed to contribute to the immune evasion of cancer cells in some cancers including EC [12]_ENREF_10. Although several studies have demonstrated that the genetic variants in MHC region is associated with an increased risk to ESCC [12]–[14], the association of GWAS with MHC in ESCC has not been well characterized.

The present study, thus, was undertaken to fine mapping the association signal within MHC region with susceptibility to ESCC with conditional logistic regression analysis (CLRA) on our previous GWAS of ESCC in Chinese population [9]. The identified promising SNPs were further correlated with gender, age, family history, alcohol consumption, smoking status and clinicopathology characteristics. Immunohistochemistry was performed to explore the protein expression pattern of the related genes in ESCC and neighboring normal tissues.

Materials and Methods

Ethics Statement

The study was approved by the ethical review committee of Zhengzhou University and conducted according to Declaration of Helsinki principles. Written informed consent was obtained from all the participants.

Recruitment of patients and controls

All the 1,089 patients with ESCC and the 1,763 normal control subjects were from the high incidence area in northern China and obtained from Endoscopic Screening Centers within multiple hospitals for early detection of ESCC. The enrollment criteria for the cases and controls in this study were the same as described in our previous GWAS study [9]_ENREF_11_ENREF_11. In brief, there were no restrictions on age, gender, ethnicity, risk factor exposures or cancer-stages for cases or controls; local residents in the same geographic areas; all the cases were histopathologically confirmed ESCC and no prior history of treatment for ESCC other than surgery; the tumor location, gross type, degree of differentiation, the regional lymph node metastasis and clinicopathological stage were recorded based on the International Statistical Classification of Diseases and Related Health Problems (10th Revision, WHO, 2006, online); all the controls subjects in this study, without any family history of ESCC, were undergone esophageal endoscopic examination to exclude the early ESCC.

SNP genotyping

Genomic DNAs were extracted from peripheral blood leukocytes using a standard procedure using Flexi Gene DNA kits (QIAGEN, Germany) as previously described [9]. Briefly, DNA concentration was normalized to 50 ng/μl (diluted in 10 mM Tris/L mM EDTA) with a Nanodrop Spectrophotometer (ND-1000). Approximately 200 ng of genomic DNA was used for genotyping analysis. Each sample was whole-genome amplified, fragmented, precipitated and resuspended in appropriate hybridization buffer. Denatured samples were hybridized on prepared Illumina Human 610-Quad BeadChips. After hybridization, the BeadChips oligonucleotides were extended by a single labeled base, which was detected by fluorescence imaging with an Illumina Bead Array Reader. Normalized bead intensity data obtained for each sample were loaded into the Illumina BeadStudio 3. The genome-wide genotyping analysis was conducted using Illumina Human 610-Quad BeadChips at the Key Laboratory of Dermatology (Anhui Medical University), Ministry of Education, China, Hefei, Anhui, China. All genotyping passed through our quality control procedure and there was no significant deviation from Hardy-Weinberg proportions in control population.

As previous described [9], the GWAS was carried out in 1,089 ESCC cases and 1,763 controls of the subjects in northern China using Illumina Human 610-Quad BeadChips. After SNPs and samples based quality control filtering, 506,666 SNPs as well as 1,077 cases and 1,733 controls were left for further analyses.

Principal components analysis

As described in our previous GWAS [9], principal components analysis (PCA) was used twice in our data. In brief, firstly, PCA of our 2,810 GWAS samples (1,077 cases and 1,733 controls) alone or in combination with 206 reference HapMap samples was performed. The HapMap samples were drawn from: Yoruba in Ibadan, Nigeria (YRI) (n = 57), Japanese in Tokyo, Japan (JPT) (n = 44), Han Chinese in Beijing, China (CHB) (n = 45) and CEPH (Utah residents with ancestry from northern and western Europe) (CEU) (n = 60). This analysis identified 1,181 genetically unmatched subjects. After removal these samples, we got 1629 matched subset samples. Second, 1,629 GWAS samples (937 cases and 692 controls) alone or in combination with 206 reference HapMap samples was performed. The genome-wide χ2 inflation factor λ was reduced from 1.505 to 1.075 after removal of the 1181 unmatched samples, suggesting that population structure was not a major confounder in the matched-sample analysis.

Selection criteria of SNPs

All the 6,252 SNPs in MHC region (Chr. 6: 25–34 Mb) were obtained for further study based on the 1,629 samples (937 cases and 692 controls). Among these SNPs, 1049 SNPs failed in missingness test >0.1; 2,607 SNPs failed in frequency test (MAF<0.05), and 163 SNPs failed in Hardy-Weinberg equilibrium (HWE) test (p<0.001). Finally, 3,515 SNPs in 937 ESCC cases and 692 controls were left for downstream analysis. Simultaneously, we selected the first-choice set of promising SNPs for further analysis based on the considerations of each SNP in terms of the P value, frequency of minor allele and biological characteristics of the related genes in carcinogenesis. Finally, 24 SNPs from 17 loci were selected for downstream analysis.

Clinicopathological characteristics measurement

Smoking status was defined as never, former and current smokers. Alcohol consumption was defined as nondrinkers who reported never drinking and drinkers who drank regularly in the 1 year prior to age at diagnosis for cases. Subjects were asked to report whether any first-, second-, and third-degree relatives had been diagnosed with other upper gastrointestinal tract cancers, including EC, gastric cancer or gastric cardia adenocarcinoma. Positive family history of EC was defined as occurrence at least 2 upper gastrointestinal tract cancer cases within three degree relatives, which was the same as we previously reported [9]. The tumor location was recorded based on the International Statistical Classification of Diseases and Related Health Problems (10th Revision, WHO, 2006, online) and was classified as upper (20 to 25 cm from incisors), middle (>25 to 30 cm from incisors) and lower thoracic (>30 to 40 cm from incisors) (7th Edition of the AJCC).

Immunohistochemical staining in ESCC and neighboring normal tissues

To understand the protein expression pattern in the tissues, serial paraffin block tissue sections were obtained from of the surgically resected ESCC specimen and neighboring normal tissues. Sixty-seven ESCC cases with 60 matched normal tissues (7 pieces of tissue were off-chip during the histological dehydration process) for TRIM27 immunostaining and 49 ESCC cases with 41 matched normal tissues (8 pieces of tissue were off-chip during the histological dehydration process) for HLA-DQA1 immunostaining (ProteinTech Group, Inc., Chicago, USA). All the patients were from Linzhou, Henan Province, the high incidence area for EC in northern China and had not received treatment other than surgery. The avidin-biotin-peroxidase complex (ABC) method was applied as previously described [15]_ENREF_12. In brief, after dewaxing, inactivating, endogenous peroxidase activity and blocking cross-reactivity with pre-immune serum, the sections were incubated over night at 4°C with the primary antibodies (antibody for TRIM27 was diluted at 1∶150 and the one for HLA-DQA1 was 1∶100). Localization of the primary antibodies was achieved by subsequent incubation of a biotinylated anti-primary antibody, an avidin-biotin complex conjugated to horseradish peroxidase, and diaminobenzidine (Vectastain Elite Kit). The slides were washed 3 times with phosphate-buffered saline (PBS) after incubation. As negative controls, some slides were subjected to normal serum blocking and omission of the primary antibody. Intense nuclear or cytoplasm staining was the criterion for a “positive” reaction. We applied the criteria established by our laboratory previously [15] to describe the patterns of positive result as follows: “scattered”, in which only some isolated positive cells were identified; “papillary”, where immunostain-positive cells were identified only in the esophageal epithelial papillary area; “focal”, where wide clusters of positive cells were seen in some areas of the epithelia; and “diffuse”, in which the sheets of positive cells were found throughout most areas of the lesions. Two co-authors, the experienced pathologists (L-D W and J-L R) were specified to read the immunostaining slides.

Statistical analysis

The conditional logistic analysis (CLRA) with backwards selection was applied in this study to identify the strongest associated SNP by the SAS statistical software package (SAS Institute Inc., Cary, NC). All the 24 SNPs identified from the GWAS study were taken as independent parameters in this CLRA model. ORs and 95% CIs were calculated as estimates of the relative risk for variables. Data were analyzed with PLINK 1.06, STATISTICAL ANALYSIS SYSTEM (SAS) version 8.1 and SPSS 17.0 (SPSS Inc, Chicago, IL). The Bonferroni corrected P value based on multiple testing of 3,515 SNPs and was set at 0.05. In all stratification analysis, comparisons of allele frequencies were made using the Pearson's χ2 test and the statistical significance was set at P<0.05.

Results

SNP genotyping and quality control and CLRA

After strict quality control to the SNPs, the 24 promising SNPs selected for association testing are listed in Table 1 along with the Bonferroni corrected P values based on multiple testing of 3,515 SNPs. Since SNP rs17533090 was highly correlated with SNP rs35399661 and could be completely represented by the latter (D’ = 1/r2 = 1 in HapMap CHB data), only the rs35399661 locus was kept for further statistical analysis. Of the 24 promising SNPs analyzed, we identified three independent SNPs in the MHC region that associated with ESCC through CLRA: rs35399661 (P = 6.07E-06, OR = 1.71, 95%CI = 1.35–2.17), rs3763338 (P = 1.62E-05, OR = 0.63, 95%CI = 0.50–0.78) and rs2844695 (P = 7.60E-05, OR = 0.74, 95%CI = 0.64–0.86) (Table 1). The Bonferroni corrected P values for rs35399661, rs3763338 and rs2844695 were 0.021, 0.057 and 0.267, respectively (Table 1). These three SNPs in fact represent three independent loci containing multiple candidate genes which lie within each linkage disequilibrium (LD) block on chromosome 6p21.31 (Figures 1-3), which were located at the genes of HLA-DQA1, TRIM27, and DPCR1 (Table 1). Signal I surrounding rs35399661 covers two (HLA-DQA1 and HLA-DRB1) genes (Figure 1); Signal II surrounding rs3763338 covers four (TRIM27, ZNF311, OR2W1 and OR2B3) genes (Figure 2); Signal III surrounding rs2844695 covers two (DPCR1 and C6 or f205) genes (Figure 3).

[Figure omitted. See PDF.]

Figure 1. The LD pattern of signal I at 6p21.31 in Asian population (CHB+JPT).

https://doi.org/10.1371/journal.pone.0090438.g001

[Figure omitted. See PDF.]

Figure 2. The LD pattern of signal II at 6p21.31 in Asian population (CHB+JPT).

https://doi.org/10.1371/journal.pone.0090438.g002

[Figure omitted. See PDF.]

Figure 3. The LD pattern of signal III at 6p21.31 in Asian population (CHB+JPT).

https://doi.org/10.1371/journal.pone.0090438.g003

[Figure omitted. See PDF.]

Table 1. Summary of association results and conditional logistic regression analysis of the 24 SNPs.

https://doi.org/10.1371/journal.pone.0090438.t001

The effects of SNPs on clinical characteristics

To further correlate the three promising SNPs identified above with clinicopathological characteristics of ESCC, we performed stratified analyses of rs35399661, rs3763338 and rs2844695 by gender, age, family history, alcohol consumption, smoking status, tumor location, gross type, cancer cell differentiation, regional lymph node metastasis and pathological staging (Tables 2–3). The χ2 test of each subgroup showed the rs2844695 was preferentially associated with younger ESCC cases (P = 0.009, OR = 1.32, 95%CI = 1.07–1.62) (Table 2). However, the other clinicopathological characteristics did not show significant association with these three promising SNPs (Tables 2–3), indicating the potential values of these SNPs in clinic.

[Figure omitted. See PDF.]

Table 2. Stratified analyses of rs35399661, rs3763338 and rs2844695 by the clinical characteristics in ESCC.

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

[Figure omitted. See PDF.]

Table 3. Stratified analyses of rs35399661, rs3763338 and rs2844695 by the clinicopathology characteristics in ESCC.

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

Immunohistochemical analysis for TRIM27 and HLA-DQA1 protein

To further demonstrated the clinical relevance of the two SNPs showing the association with ESCC, expression of TRIM27 (for rs3763338) and HLA-DQA1 (for rs35399661) proteins in ESCC and the neighboring normal tissues was determined using immunohistochemical analysis (Figure 4). Marked expression of TRIM27 was observed in 77.6% of ESCC sections as compared with 47.8% of neighboring normal tissues (P<0.001). HLA-DQA1 staining was seen in 69.4% of ESCC samples versus 26.8% of neighboring normal tissues (P<0.001). Expression of TRIM27 and HLA-DQA1 was analyzed with the stratification of different sex, age, family history, high/low-incidence areas, the regional lymph node metastasis and tumor location (Tables 4 and 5). The overexpression of HLA-DQA1 was correlated significantly with age (P = 0.001) and family history (P<0.001) (Table 5). There was no correlation between TRIM27 expression and clinical characteristics (Table 4).

[Figure omitted. See PDF.]

Figure 4. The hematoxylin-eosin staining (A, 10X) and immunohistochemical staining for TRIM27 (B, 10X and E, 20X) and HLA-DQA1 (D, 20X and F, 20X) on ESCC tissue.

C (10X) was the negative control of A (omission of the first antibody application). A, B and C were the serial tissue sections from the same ESCC patient. E and F were also from the same patient.

https://doi.org/10.1371/journal.pone.0090438.g004

[Figure omitted. See PDF.]

Table 4. Relationship of TRIM27 protein immunostaining and clinical characteristics of ESCC.

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

[Figure omitted. See PDF.]

Table 5. Relationship of HLA-DQA1 protein immunostaining and clinical characteristics of ESCC.

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

Discussion

Three GWAS of ESCC conducted among Chinese Han populations have been published [7]–[9]. Though previous studies have investigated the role of individual genetic variation in the etiology of esophageal cancers, the effect of high density SNP genotypes of the MHC region in the etiology of ESCC is still unclear, and little genome-wide data are available. Strengths of this study include that it is conducted in high-incidence areas with highly homogeneous population and that the controls are well-matched with ESCC cases with similar environmental exposures.

In the present study, we have performed the first GWAS of ESCC in the MHC region on the subjects from high risk area in northern China and found three important independent susceptibility loci containing three biologically interesting candidate genes, i.e., HLA-DQA1, TRIM27 and DPCR1. TRIM27, also known as ret finger protein (RFP), which is characterized by a conserved RING finger, a B-box, and a coiled-coil domain (together called RBCC) [17], [18], encodes a member of the tripartite motif (TRIM) family [19]_ENREF_18. Currently, the biological function of RFP has not been characterized; however, many of the TRIM family members sharing the RBCC moiety participate in the control of cell survival [19]_ENREF_18. It has been shown that the RBCC moiety is required for the transforming capacities of these TRIM oncogenes_ENREF_19. The role of RFP influences apoptotic pathways [19]_ENREF_18 and acquire oncogenic activity when fused to kinases by chromosomal rearrangements [20]_ENREF_20. Compared to many normal tissues, cancer cells are highly sensitized to apoptotic signals, and survive only because they have acquired lesions [21]_ENREF_21. Indeed, immunohistochemical analysis for TRIM27 in our study demonstrated the higher positive rate of the TRIM27 protein in ESCC samples than in neighboring normal epithelia (P<0.001). Therefore, RFP is likely to play an important role in esophageal carcinogenesis through modulating apoptotic pathways.

HLA-DQA1 belongs to the HLA class II alpha chain paralogues. The class II molecule is a heterodimer consisting of an alpha (DQA) and a beta chain (DQB), both anchored in the membrane. HLA-DQA1 plays a central role in the immune system by presenting peptides derived from extracellular proteins [22]_ENREF_22. Class II molecules are expressed in antigen presenting cells (B lymphocytes, dendritic cells and macrophages). Within the DQ molecule both the alpha chain and the beta chain contain the polymorphisms specifying the peptide binding specificities, resulting in up to four different molecules. MHC class II molecules bind peptides derived from proteins that have entered the endocytic pathway and presented them at the cell surface for interaction with CD4+ T cells [23]_ENREF_23. The presence of HLA-DQ expression in premalignant lesions and on some tumor cells appears to confer an advantage to the host in terms of restricted tumor growth [24]_ENREF_24 and survival [24] _ENREF_24through their role as initiators of CD4+ T helper cell responses against the tumor. Simultaneously, our preliminary immunohistochemical analysis showed that the positive rate of HLA-DQA1 staining in ESCC samples was apparently higher than the neighboring normal epithelia. The results further supported the potential role of HLA-DQA1 in esophageal carcinogenesis.

DPCR1, also known as diffuse panbronchiolitis critical region 1, located between HLA-B and HLA-A on chromosome 6p21.33, is classified as one of the MHC class I molecules. The DPCR1 gene may contain markers for diagnosis of diffuse pan-bronchiolitis, a bronchiolar disease that affects human airways [25]. However, to our knowledge, the association between DPCR1 variations and risk of ESCC has not yet been investigated. Our findings may be meaningful for future studies validating the function of DPCR1. Interestingly, the LD pattern of Signal III (Figure 3) showed that rs2844695 was an independent locus.

Previous studies have indicated that smoking, alcohol consumption, aging and gender are risk factors for ESCC [26], [27]_ENREF_27_ENREF_25. Our results demonstrated that rs2844695 was preferentially associated with younger ESCC cases. The other factors (gender, age, family history, alcohol consumption, smoking, tumor location, gross type, degree of differentiation, the regional lymph node metastasis and pathological stage) did not significantly alter the effects ofrs35399661 and rs3763338 on the risk to ESCC, indicating that the three SNPs identified as most promising from our study could provide orthologous information to existing clinicophthological covariates.

We also examined the correlation of TRIM27 and HLA-DQA1 expression in ESCC with clinical characteristics. Our results showed that the expression of TRIM27 and HLA-DQA1 was higher in ESCC tissues than in neighboring normal tissues. The results further supporting the potential role of TRIM27 and HLA-DQA1 in esophageal carcinogenesis. The expression of HLA-DQA1 was correlated with age and family history of ESCC with the older ESCC cases having higher HLA-DQA1 expression than the younger ones. Furthermore, the percentage of patients expressing HLA-DQA1 is higher for negative family history than positive family history. The result suggested that HLA-DQA1 may be potential in inhibiting the malignant properties in ESCC.

The limitation of the present study is the absence of replication in another independent samples for rs35399661 (HLA-DQA1 genes), rs3763338 (TRIM27 genes) and rs2844695 (DPCR1 genes). Interestingly, recent studies have demonstrated that polymorphisms occurred in HLA-DQA1 increase the risk and prognosis to lung squamous cell carcinoma [28] and gastric cancer [29]. Further replication studies would be desirable to elucidate the role of these SNPs on the susceptibility, and possibly in the prognosis for ESCC.

Conclusion

This study for the first time provides evidence that multiple genetic factors within the MHC region confer risk to ESCC on the subjects from high-risk area in northern China. Further dissection of the roles of these loci will likely to lead to insights into the etiology of this rapidly evolving and fatal cancer.

Acknowledgments

We thank all the patients and their family members whose contributions made this work possible. We acknowledge personnel and the working arrangements with the Genotyping at Key Laboratory of Dermatology (Anhui Medical University), Ministry of Education, China, Hefei, Anhui for genotyping and statistical and bioinformatics analyses. We also thank Dr. Liang Dan Sun at Anhui Medical University for help with the manuscript preparation and helpful discussion. We are also grateful to Ph.D Le Ming Shi in Fudan University, Shanghai, China for his carefully edit of this manuscript.

Author Contributions

Conceived and designed the experiments: LDW. Performed the experiments: FFS WBY FYZ YP XKZ YJ XS BL XNH ST YL GY LSC YLL YLH XML JLR. Analyzed the data: LDW FFS XS. Contributed reagents/materials/analysis tools: LDW FFS WBY FYZ XKZ XS XNH ST YL GY LSC YLL YLH XML. Wrote the paper: LDW FFS.

Citation: Shen F-F, Yue W-B, Zhou F-Y, Pan Y, Zhao X-K, Jin Y, et al. (2014) Variations in the MHC Region Confer Risk to Esophageal Squamous Cell Carcinoma on the Subjects from High-Incidence Area in Northern China. PLoS ONE 9(3): e90438. https://doi.org/10.1371/journal.pone.0090438

References

1. Jemal A, Bray F (2011) Center MM, Ferlay J, Ward E, et al (2011) Global cancer statistics. CA Cancer J Clin 61: 69–90.

2. Haocai L, Yu SL (1983) Esophageal cancer in migrants from high- or low-risk areas in China. Ecol Dis 2: 249–253.

3. Yang CS (1980) Research on esophageal cancer in China: a review. Cancer Res 40: 2633–2644.

4. Tran GD, Sun XD, Abnet CC, Fan JH, Dawsey SM, et al. (2005) Prospective study of risk factors for esophageal and gastric cancers in the Linxian general population trial cohort in China. Int J Cancer 113: 456–463.

5. Guo W, Jiang YG (2009) Current gene expression studies in esophageal carcinoma. Curr Genomics 10: 534–539.

6. Daly JM, Fry WA, Little AG, Winchester DP, McKee RF, et al. (2000) Esophageal cancer: results of an American College of Surgeons Patient Care Evaluation Study. J Am Coll Surg 190: 562–572 discussion 572–563.

7. Wu C, Kraft P, Zhai K, Chang J, Wang Z, et al. (2012) Genome-wide association analyses of esophageal squamous cell carcinoma in Chinese identify multiple susceptibility loci and gene-environment interactions. Nat Genet 44: 1090–1097.

8. Abnet CC, Freedman ND, Hu N, Wang Z, Yu K, et al. (2010) A shared susceptibility locus in PLCE1 at 10q23 for gastric adenocarcinoma and esophageal squamous cell carcinoma. Nat Genet 42: 764–767.

9. Wang LD, Zhou FY, Li XM, Sun LD, Song X, et al. (2010) Genome-wide association study of esophageal squamous cell carcinoma in Chinese subjects identifies susceptibility loci at PLCE1 and C20orf54. Nat Genet 42: 759–763.

10. Dura P, van Veen EM, Salomon J, te Morsche RH, Roelofs HM, et al. (2013) Barrett associated MHC and FOXF1 variants also increase esophageal carcinoma risk. Int J Cancer 133: 1751–1755.

11. Horton R, Wilming L, Rand V, Lovering RC, Bruford EA, et al. (2004) Gene map of the extended human MHC. Nat Rev Genet 5: 889–899.

12. Liu Q, Hao C, Su P, Shi J (2009) Down-regulation of HLA class I antigen-processing machinery components in esophageal squamous cell carcinomas: association with disease progression. Scand J Gastroenterol 44: 960–969.

13. Nie Y, Yang G, Song Y, Zhao X, So C, et al. (2001) DNA hypermethylation is a mechanism for loss of expression of the HLA class I genes in human esophageal squamous cell carcinomas. Carcinogenesis 22: 1615–1623.

14. Rajendra S, Ackroyd R, Karim N, Mohan C, Ho JJ, et al. (2006) Loss of human leucocyte antigen class I and gain of class II expression are early events in carcinogenesis: clues from a study of Barrett's oesophagus. J Clin Pathol 59: 952–957.

15. Wang LD, Hong JY, Qiu SL, Gao H, Yang CS (1993) Accumulation of p53 protein in human esophageal precancerous lesions: a possible early biomarker for carcinogenesis. Cancer Res 53: 1783–1787.

16. Price AL, Patterson NJ, Plenge RM, Weinblatt ME, Shadick NA, et al. (2006) Principal components analysis corrects for stratification in genome-wide association studies. Nat Genet 38: 904–909.

17. Takahashi M, Buma Y, Iwamoto T, Inaguma Y, Ikeda H, et al. (1988) Cloning and expression of the ret proto-oncogene encoding a tyrosine kinase with two potential transmembrane domains. Oncogene 3: 571–578.

18. Saurin AJ, Borden KL, Boddy MN, Freemont PS (1996) Does this have a familiar RING? Trends Biochem Sci 21: 208–214.

19. Dho SH, Kwon KS (2003) The Ret finger protein induces apoptosis via its RING finger-B box-coiled-coil motif. J Biol Chem 278: 31902–31908.

20. Takahashi M, Cooper GM (1987) ret transforming gene encodes a fusion protein homologous to tyrosine kinases. Mol Cell Biol 7: 1378–1385.

21. Vousden KH, Lu X (2002) Live or let die: the cell's response to p53. Nat Rev Cancer 2: 594–604.

22. Jordanova ES, Philippo K, Giphart MJ, Schuuring E, Kluin PM (2003) Mutations in the HLA class II genes leading to loss of expression of HLA-DR and HLA-DQ in diffuse large B-cell lymphoma. Immunogenetics 55: 203–209.

23. Germain RN (1994) MHC-dependent antigen processing and peptide presentation: providing ligands for T lymphocyte activation. Cell 76: 287–299.

24. Lazaris AC, Theodoropoulos GE, Davaris PS, Panoussopoulos D, Nakopoulou L, et al. (1995) Heat shock protein 70 and HLA-DR molecules tissue expression. Prognostic implications in colorectal cancer. Dis Colon Rectum 38: 739–745.

25. Matsuzaka Y, Tounai K, Denda A, Tomizawa M, Makino S, et al. (2002) Identification of novel candidate genes in the diffuse panbronchiolitis critical region of the class I human MHC. Immunogenetics 54: 301–309.

26. Hiyama T, Yoshihara M, Tanaka S, Chayama K (2007) Genetic polymorphisms and esophageal cancer risk. Int J Cancer 121: 1643–1658.

27. Morita M, Kumashiro R, Kubo N, Nakashima Y, Yoshida R, et al. (2010) Alcohol drinking, cigarette smoking, and the development of squamous cell carcinoma of the esophagus: epidemiology, clinical findings, and prevention. Int J Clin Oncol 15: 126–134.

28. Kohno T, Kunitoh H, Mimaki S, Shiraishi K, Kuchiba A, et al. (2011) Contribution of the TP53, OGG1, CHRNA3, and HLA-DQA1 genes to the risk for lung squamous cell carcinoma. J Thorac Oncol 6: 813–817.

29. Huang LM, Cheng Y, Yu DK, Zhai K, Tan W, et al. (2012) [Association between HLA-DQA1 gene copy number polymorphisms and susceptibility to gastric cancer]. Zhonghua Zhong Liu Za Zhi 34: 269–271.