Abstract
Host genetics play a vital role in determining clinical outcomes of hepatitis B virus (HBV) infection. To identify novel susceptibility loci to HBV progression, we carried out a genome-wide association study with DNA pooling. This study assessed the relationship between 8 highly-ranked SNPs selected from our DNA pool and disease progression of HBV infection in two independent case-control studies. The first population included 628 asymptomatic HBV carriers (AsC) and 1729 progressed HBV carriers recruited from Hubei Province in south China. The second population was composed of 226 AsC and 215 progressed HBV carriers recruited from Shandong Province in north China. Of the 8 SNPs, variant rs11866328 (G/T), located in the glutamate receptor ionotropic N-methyl D-aspartate 2A (GRIN2A) gene, was replicated and had significant associations with disease progression of HBV infection in the DNA pooling stage both in the Hubei (OR 1.65; 95% CI 1.34,2.02; p=1.96×10−6; additive model), and in the Shandong (OR 1.73; 95% CI 1.14,2.65; p=1.00×10−2; additive model) population. Polymorphism rs11866328 in the GRIN2A gene might be a genetic variant underlying the susceptibility of HBV carriers to disease progression.
Introduction
Materials and Methods
Study participants
Two different populations were enrolled in this study. The first population was unrelated southern Han Chinese recruited from Hubei Province, and the second was unrelated northern Chinese Han recruited from Shandong Province. They were collected from September 2007 to March 2010. Each population was composed of three subgroups: naïve control, AsC, and progressed HBV carriers (acute liver failure [ALF], liver cirrhosis [LC], and HCC). The Hubei population included 210 naïve controls, 628 AsC, and 1729 progressed HBV carriers (261 with ALF, 631 with LC, and 837 with HCC), recruited from Wuhan Tongji Hospital, Wuhan Union Hospital, and Traditional Chinese Medicine Hospital of Hubei Province. The Shandong population was composed of 210 naïve controls, 226 AsC, and 215 progressed HBV carriers (2 with ALF, 155 with LC, and 58 with HCC), which was recruited from Qingdao Hospital for Infectious Disease and the Affiliated Hospital of Binzhou Medical College. All HBV carriers were positive for both HBsAg and anti-HBcAb for at least 6 mo. The criteria for AsC, ALF, LC, and HCC complied with a previous description (20,21), with minor modifications (Supplementary Table 1; see online supplementary material at
Informed written consent was obtained from all participants, and the study was performed with the approval of the ethical committee of Huazhong University of Science and Technology.
DNA isolation, DNA pooling, and genotyping
Genomic DNA was isolated from EDTA-anticoagulated blood using a QuickGene DNA whole blood kit S (DB-S) with QuickGene-Mini80 equipment (Fujifilm, Tokyo, Japan). Affymetrix Genome-Wide Human Mapping SNP6.0 Arrays were performed for DNA pooling, which were constructed by pooling 120 ng of DNA from each participant. Four independent pools were created: case A was the ALF group (n=86), case B was the LC group (n=88), case C was the HCC group (n=90), and case D was the AsC group (n=66), that was considered as control. Twelve chips (each pool was replicated in triplicate) were finished according to the manufacturer's instructions. The detailed process of DNA pooling has been submitted to GEO (
Statistical analysis
The Hardy-Weinberg equilibrium of genotype frequencies was evaluated online using the website of the Technical University Munich (
Results
Characteristics of HBV carriers are shown in Table 1. Although an effort was made to obtain a good match for gender and age between AsC and progressed HBV carriers, there were more men in the progressed HBV carrier group than the AsC group (p<0.01) in our hospital-based case-control study. The subjects in the progressed HBV group were younger in Hubei Province (p=0.01), and older in Shandong Province (p<0.01). More individuals had first-degree relatives with HBV infection in the progressed HBV carrier group than the AsC group (p<0.01). A large difference in alcohol consumption between AsC and progressed HBV carriers was not evident. Genotype distributions of the participants are shown in Table 2. All genotypic distributions of variants in naive controls, AsC, and progressed HBV carriers were consistent with Hardy-Weinberg equilibrium (p≥0.05).
Compared with AsC (p=0.01).
Compared with AsC p<0.01).
HBV family history: individuals who have known parents or siblings who are HBV carriers.
Defined as alcohol consumption of >40 g/wk for men and >20 g/wk for women.
The differences of clinical characteristics were calculated using chi-square test.
M, male; F, female; SD, standard deviation; no., number; AsC, asymptomatic HBV carriers; HBeAg, hepatitis B early antigen.
Number of minor homozygote/number of heterozygote/number of major homozygote.
Dominant model.
Additive model.
Total numbers for each polymorphism vary because not all samples were successfully genotyped.
SNP, single nucleotide polymorphism; Control, asymptomatic HBV carriers; Case, progressed HBV carriers.
p Values, odds ratios (OR) and 95% confidence intervals (CI) were calculated by logistic regression compared with asymptomatic HBV carriers, adjusting for gender, age and alcohol consumption.
Of the 8 SNPs, 3 SNPs: rs11866328 (OR 1.65; 95% CI 1.34,2.02; p=1.96×10−6; additive model), rs7861010 (OR 1.31; 95% CI 1.06,1.61; p=1.20×10−2; additive model), and rs12206945 (OR 1.74; 95% CI 1.21,2.51; p=3.00×10−3; dominant model) were replicated and had strong correlations with disease progression of HBV infection, as in the DNA pooling stage in the Hubei population, based on multivariable logistic regression analysis with adjustment for age, sex, and alcohol use (Table 2). The three strongly correlated SNPs were further genotyped in an additional independent case-control population recruited from Shandong Province in northern China using the same TaqMan assay. On the basis of logistic regression adjusting for sex, age, and alcohol consumption, only SNP rs11866328 was found to affect susceptibility to HBV progression (OR 1.73; 95% CI 1.14,2.65; p=1.00×10−2; additive model) (Table 2). Combining the two case-control studies, the association of rs11866328 with HBV progression reached a p value of 1.60×10−8 (OR 1.68; 95% CI 1.40,2.02; additive model) (Table 3). In stratification analysis of SNP rs11866328 in all studied HBV carriers, the results remained significant in males (OR 1.61; 95% CI 1.30,1.99; p=1.10×10−5; additive model), females (OR 1.78; 95% CI 1.25,2.53; p=1.00×10−3; additive model), and patients 35 y or older (OR 1.63; 95% CI 1.35,1.97; p=3.10×10−7; additive model) (Table 3). As shown in Table 4, GG genotype was a risk factor for ALF (p=1.70×10−2), LC (p=6.60×10−7), and HCC (p=4.90×10−4), in the Hubei population. This variant might perform a crucial role in determining the outcomes of HBV infection.
The number of progressed HBV carriers for SNP rs11866328 was 1723 because not all samples were successfully genotyped.
p Values, odds ratios (OR), and 95% confidence intervals (CI) were calculated by logistic regression compared with asymptomatic HBV carriers, adjusting for gender, age, and alcohol consumption.
Recessive model (GG/TG+TT); the additive model was p=0.11.
The number of progressed HBV carriers for SNP rs11866328 was 1723 because not all samples were successfully genotyped.
p Values, odds ratios (OR), and 95% confidence intervals (CI) were calculated by logistic regression compared with asymptomatic HBV carriers, adjusting for gender, age, and alcohol consumption.
AsC, asymptomatic HBV carrier; ALF, acute liver failure; LC, liver cirrhosis; HCC, hepatocellular carcinoma; HBV, hepatitis B virus.
Discussion
We performed this study using DNA pooling in the screening stage. DNA pooling has been confirmed to be an effective and efficient method to select candidate susceptibility loci for follow-up by individual genotyping (23,25 –27). The results of DNA pooling were analyzed using silhouettes, which were introduced in 1987 as a general graphical aid for interpretation and validation of cluster analysis (24,28).
There was a strong difference in the distribution of sex and age between our AsC and progressed HBV carriers (Table 1), hence we controlled for the effect of sex and age in all analyses by logistic regression. To decrease the bias of gender on the effect estimates, we conducted a stratification analysis for sex. Interestingly, the association remained significant for male and female patients (Table 3). Long-term cohort studies show that advanced age is a risk factor for HBV persistence and disease progression (29 –31), so we limited our AsC subjects to those 35 y or older to minimize the deviation of age. In stratification analysis for age (≥35 y), subjects bearing the GG genotype showed a strong correlation with adverse outcomes of HBV infection (Table 3). This is in accord with our observation that the variant rs11866328 is associated with disease progression in HBV carriers. There were more individuals with parents or siblings who were HBV carriers in the progressed HBV carrier group, which is consistent with the fact that genetic factors play an important role in HBV progression.
The genotype frequencies of rs11866328 in our naive subjects (TT: 3.1%, TG: 35.0%, GG: 61.9%) were similar to the frequency distributions seen in chronic hepatitis B patients reported in dbSNP (TT: 2.2%, TG: 31.1%, GG: 66.7%;
The GRIN2A gene is located on chromosome 16p13.2, which encodes the N-methyl-D-aspartate (NMDA) glutamate receptor NR2A subunit. NMDA receptors (NMDARs) are a major class of excitatory neurotransmitter receptors in the central nervous system. They form glutamate-gated ion channels that are highly permeable to calcium, and mediate activity-dependent synaptic plasticity (34). SNP rs11866328 is located in the tenth intron of the GRIN2A gene, and does not alter the amino acid sequence of the protein. Further studies are needed to more precisely map and identify the causative polymorphism, and to clarify which genes drive the genetic association. Further functional studies are required to establish the role of this gene in the pathogenesis of disease progression of HBV infection. We realize that the ideal model for studying disease progression in HBV carriers is a community-based, long-term cohort study. From this initial population an HBV carrier cohort can be established, and time-to-progression analyses may then be conducted.
Conclusion
We revealed the important role of SNP rs11866328 on the gene GRIN2A in disease progression of HBV infection in two independent case-control studies in Han Chinese HBV carriers. To our knowledge, this is the first report showing that rs11866328 may affect an individual's susceptibility to HBV progression. Our research indicates that SNP rs11866328 may be a genetic variant underlying chronic HBV carriers' susceptibility to disease progression.
Footnotes
Acknowledgments
This study was supported by the National Basic Research Program of China (973 Program) grant 2007CB512900, and the Chinese National Natural Science Foundation grant 30872237.
Author Disclosure Statement
No competing financial interests exist.
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
