Development of a SNP set for human identification: A set with high powers of discrimination which yields high genetic information from naturally degraded DNA samples in the Thai population.

This study describes the development of a SNP typing system for human identification in the Thai population, in particular for extremely degraded DNA samples. A highly informative SNP marker set for forensic identification was identified, and a multiplex PCR-based Invader assay was developed. Fifty-one highly informative autosomal SNP markers and three sex determination SNP markers were amplified in two multiplex PCR reactions and then detected using Invader assay reactions. The average PCR product size was 71 base pairs. The match probability of the 54-SNP marker set in 124 Thai individuals was 1.48×10(-21), higher than that of STR typing, suggesting that this 54-SNP marker set is beneficial for forensic identification in the Thai population. The selected SNP marker set was also evaluated in 90 artificially degraded samples, and in 128 naturally degraded DNA samples from real forensic casework which had shown no profiles or incomplete profiles when examined using a commercial STR typing system. A total of 56 degraded samples (44%) achieved the matching probability (PM) equivalent to STR gold standard analysis (successful genotyping of 44 SNP markers) for human identification. These data indicated that our novel 54-SNP marker set provides a very useful and valuable approach for forensic identification in the Thai population, especially in the case of highly to extremely degraded DNA. In summary, we have developed a set of 54 Thai-specific SNPs for human identification which have higher discrimination power than STR genotyping. The PCRs for these 54 SNP markers were successfully combined into two multiplex reactions and detected with an Invader assay. This novel SNP genotyping system also yields high levels of genetic information from naturally degraded samples, even though there are much more difficult to recover than artificially degraded samples.

Identification of a susceptibility locus for severe adolescent idiopathic scoliosis on chromosome 17q24.3.

Adolescent idiopathic scoliosis (AIS) is the most common spinal deformity, affecting around 2% of adolescents worldwide. Genetic factors play an important role in its etiology. Using a genome-wide association study (GWAS), we recently identified novel AIS susceptibility loci on chromosomes 10q24.31 and 6q24.1. To identify more AIS susceptibility loci relating to its severity and progression, we performed GWAS by limiting the case subjects to those with severe AIS. Through a two-stage association study using a total of ∼12,000 Japanese subjects, we identified a common variant, rs12946942 that showed a significant association with severe AIS in the recessive model (P=4.00 × 10(-8), odds ratio [OR]=2.05). Its association was replicated in a Chinese population (combined P=6.43 × 10(-12), OR = 2.21). rs12946942 is on chromosome 17q24.3 near the genes SOX9 and KCNJ2, which when mutated cause scoliosis phenotypes. Our findings will offer new insight into the etiology and progression of AIS.

Impact of allele copy number of polymorphisms in FCGR3A and FCGR3B genes on susceptibility to ulcerative colitis.

BACKGROUND: Polymorphisms in the Fcγ receptor genes have been implicated in several autoimmune diseases, including ulcerative colitis (UC). However, most of these reports had not taken into account the effect of copy number variation at this region. METHODS: We investigated the combined effect of allele and gene copy number of FCGR3A-158F/V and FCGR3B-NA1/NA2 on susceptibility to UC. Study subjects were composed of a total of 752 Japanese patients with UC and 2062 Japanese control subjects. To estimate allele copy number of the 2 polymorphisms, we integrated the results of PCR-based real-time Invader assay (PCR-RETINA) that measures the allelic ratio and Taqman assay that detects the total copy number. We analyzed the associations of allele copy number with UC using logistic regression model. RESULTS: Gene and allele copy numbers of FCGR3A and FCGR3B were successfully determined in more than 99.5% of the study subjects. Allele copy number of FCGR3A-158F/V demonstrated significant association with susceptibility to UC (P = 0.02), although each single-nucleotide polymorphism and copy number variation alone did not show significant association. Although allele copy number of FCGR3B-NA1/NA2 (P = 0.002) also showed significant association with UC susceptibility, this association seemed to reflect the effect of FCGR3B gene copy number. Subsequent haplotype analyses revealed a strong association of a haplotype FCGR2A-131H/R and copy number of FCGR3B gene (P = 6.5 × 10). CONCLUSIONS: Allele copy number of FCGR3A-158F/V and FCGR3B gene copy number were associated with UC susceptibility. Our results suggest that organizing handling of immune complex by FCGR3A, FCGR3B, and FCGR2A may play a crucial role in the pathogenesis of UC.

Genetic variants in GPR126 are associated with adolescent idiopathic scoliosis.

Adolescent idiopathic scoliosis (AIS) is the most common pediatric skeletal disease. We previously reported a locus on chromosome 10q24.31 associated with AIS susceptibility in Japanese using a genome-wide association study (GWAS) consisting of 1,033 cases and 1,473 controls. To identify additional AIS-associated loci, we expanded the study by adding X-chromosome SNPs in the GWAS and increasing the size of the replication cohorts. Through a stepwise association study including 1,819 cases and 25,939 controls, we identified a new susceptibility locus on chromosome 6q24.1 in Japanese (P = 2.25 × 10(-10); odds ratio (OR) = 1.28). The most significantly associated SNP, rs6570507, was in GPR126 (encoding G protein-coupled receptor 126). Its association was replicated in Han Chinese and European-ancestry populations (combined P = 1.27 × 10(-14); OR = 1.27). GPR126 was highly expressed in cartilage, and the knockdown of gpr126 in zebrafish caused delayed ossification of the developing spine. Our results should provide insights into the etiology and pathogenesis of AIS.

A genome-wide association study of HCV-induced liver cirrhosis in the Japanese population identifies novel susceptibility loci at the MHC region.

BACKGROUND & AIMS: We performed a genome-wide association study (GWAS) of hepatitis C virus (HCV)-induced liver cirrhosis (LC) to identify predictive biomarkers for the risk of LC in patients with chronic hepatitis C (CHC). METHODS: A total of 682 HCV-induced LC cases and 1045 CHC patients of Japanese origin were genotyped by Illumina Human Hap 610-Quad bead Chip. RESULTS: Eight SNPs which showed possible associations (p<1.0 × 10(-5)) at the GWAS stage were further genotyped using 936 LC cases and 3809 CHC patients. We found that two SNPs within the major histocompatibility complex (MHC) region on chromosome 6p21, rs910049 and rs3135363, were significantly associated with the progression from CHC to LC (pcombined=9.15 × 10(-11) and 1.45 × 10(-10), odds ratio (OR)=1.46 and 1.37, respectively). We also found that HLA-DQA1(*)0601 and HLA-DRB1(*)0405 were associated with the progression from CHC to LC (p=4.53 × 10(-4) and 1.54 × 10(-4) with OR=2.80 and 1.45, respectively). Multiple logistic regression analysis revealed that rs3135363, rs910049, and HLA-DQA1(*)0601 were independently associated with the risk of HCV-induced LC. In addition, individuals with four or more risk alleles for these three loci have a 2.83-fold higher risk for LC than those with no risk allele, indicating the cumulative effects of these variations. CONCLUSIONS: Our findings elucidated the crucial roles of multiple genetic variations within the MHC region as prognostic/predictive biomarkers for CHC patients.

A genome-wide association study identifies 2 susceptibility Loci for Crohn's disease in a Japanese population.

BACKGROUND & AIMS: Crohn's disease is an inflammatory bowel disease induced by multiple genetic and environmental factors. Genome-wide association studies have identified genetic factors that affect the risk for Crohn's disease in European populations, but information from other ethnic groups is scarce. We therefore investigated genetic factors associated with Crohn's disease in the Japanese population. METHODS: We performed a genome-wide association study with 372 individuals with Crohn's disease (cases) and 3389 controls, all from the Japanese population. To confirm identified associations, we performed a replication study with an independent panel of 1151 Crohn's disease cases and 15,800 controls. We also performed an association analysis using genome-wide genotype imputation in the discovery cohort. RESULTS: We confirmed associations of Crohn's disease with variants in MHC (rs7765379, P = 2.11 × 10(-59)), TNFSF15 (rs6478106, P = 3.87 × 10(-45)), and STAT3 (rs9891119, P = 2.24 × 10(-14)). We identified 2 new susceptibility loci: on chromosome 4p14 (rs1487630, P = 2.40 × 10(-11); odds ratio, 1.33), and in the SLC25A15-ELF1-WBP4 region on 13q14 (rs7329174 in ELF1, P = 5.12 × 10(-9); odds ratio, 1.27). CONCLUSIONS: In a genome-wide association study, we identified 2 new susceptibility loci for Crohn's disease in a Japanese population. These findings could increase our understanding of the pathogenesis of Crohn's disease.

Whole-genome sequencing of liver cancers identifies etiological influences on mutation patterns and recurrent mutations in chromatin regulators.

Hepatocellular carcinoma (HCC) is the third leading cause of cancer-related death worldwide. We sequenced and analyzed the whole genomes of 27 HCCs, 25 of which were associated with hepatitis B or C virus infections, including two sets of multicentric tumors. Although no common somatic mutations were identified in the multicentric tumor pairs, their whole-genome substitution patterns were similar, suggesting that these tumors developed from independent mutations, although their shared etiological backgrounds may have strongly influenced their somatic mutation patterns. Statistical and functional analyses yielded a list of recurrently mutated genes. Multiple chromatin regulators, including ARID1A, ARID1B, ARID2, MLL and MLL3, were mutated in ∼50% of the tumors. Hepatitis B virus genome integration in the TERT locus was frequently observed in a high clonal proportion. Our whole-genome sequencing analysis of HCCs identified the influence of etiological background on somatic mutation patterns and subsequent carcinogenesis, as well as recurrent mutations in chromatin regulators in HCCs.

Novel deletion mutations of OPTN in amyotrophic lateral sclerosis in Japanese.

Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disease characterized by selective motor neuron death in the brain and spinal cord. Many disease genes for ALS have been identified; however, each disease gene is responsible for very small fractions of ALS. Recently, mutations of the gene encoding optineurin (OPTN) are reported in familial and sporadic ALS. OPTN is also responsible for a small number of ALS, 3.8% of familial and 0.29% of sporadic ALS in Japanese. The low prevalence may be an underestimation due to incomplete screening of the mutation. To examine OPTN mutations more extensively, we screened the OPTN deletions using a quantitative PCR system. We examined 710 Japanese ALS subjects who had previously been found to have no OPTN mutations by a screening using a PCR-direct sequence strategy. We identified 3 kinds of deletions in 5 patients; one was homozygous, and the remaining were heterozygous. All deletions occurred due to the Alu-mediated recombination and are expected to result in null alleles. Our results suggest that the OPTN deletion mutation in ALS is not infrequent and the prevalence of the OPTN mutation in Japanese sporadic ALS is considerably high.

New pharmacogenetic test for detecting an HLA-A*31: 01 allele using the InvaderPlus assay.

BACKGROUND & AIMS: Carbamazepine (CBZ) is widely used for the treatment of epilepsy and other neurological disorders. However, 3-5% of CBZ-treated individuals suffer from cutaneous adverse drug reactions (cADRs). Recently, in a genome-wide association study, HLA-A*31:01 has been reported to be a strong genetic marker for CBZ-induced cADRs in both Japanese and European populations. As most of the available methods for HLA genotyping are laborious, the development of a simple and rapid genotyping method for HLA-A*31:01 is desirable from the viewpoint of a clinical pharmacogenetic test. METHODS: More than 1700 sequences for HLA-A alleles were obtained from the MHC database of the National Center for Biotechnology Information (dbMHC). Several HLA-A*31:01-discriminating single-nucleotide polymorphisms were selected. These SNPs were used for sequence-specific primer PCR (SSP-PCR) and for the target site of the Invader reaction. By combining SSP-PCR with a target-specific Invader reaction, we designed two sets of primers/probes for HLA-A*31:01 allele detection. The performance of both sets was evaluated using 90 Asian HapMap samples. Further evaluation was carried out using another 376 Japanese samples and 90 CEU (European) and 90 YRI (African) HapMap samples. RESULTS: Our assay specifically detected an HLA-A*31:01 allele in a total of 466 individuals of the Asian population. Furthermore, the assay correctly identified HLA-A*31:01-positive carriers from the CEU and the YRI population, respectively, implying that the assay has potential for application to other ethnic groups. CONCLUSION: We developed a new HLA-A*31:01-detecting method by a combination of SSP-PCR with target-specific InvaderPlus technology. As our assay is rapid and accurate, it is hoped that this method will be used in a pharmacogenetic test in a clinical setting to avoid CBZ-induced cADRs.

Genetic differences in the two main groups of the Japanese population based on autosomal SNPs and haplotypes.

Although the Japanese population has a rather low genetic diversity, we recently confirmed the presence of two main clusters (the Hondo and Ryukyu clusters) through principal component analysis of genome-wide single-nucleotide polymorphism (SNP) genotypes. Understanding the genetic differences between the two main clusters requires further genome-wide analyses based on a dense SNP set and comparison of haplotype frequencies. In the present study, we determined haplotypes for the Hondo cluster of the Japanese population by detecting SNP homozygotes with 388,591 autosomal SNPs from 18,379 individuals and estimated the haplotype frequencies. Haplotypes for the Ryukyu cluster were inferred by a statistical approach using the genotype data from 504 individuals. We then compared the haplotype frequencies between the Hondo and Ryukyu clusters. In most genomic regions, the haplotype frequencies in the Hondo and Ryukyu clusters were very similar. However, in addition to the human leukocyte antigen region on chromosome 6, other genomic regions (chromosomes 3, 4, 5, 7, 10 and 12) showed dissimilarities in haplotype frequency. These regions were enriched for genes involved in the immune system, cell-cell adhesion and the intracellular signaling cascade. These differentiated genomic regions between the Hondo and Ryukyu clusters are of interest because they (1) should be examined carefully in association studies and (2) likely contain genes responsible for morphological or physiological differences between the two groups.

A genome-wide association study identified AFF1 as a susceptibility locus for systemic lupus eyrthematosus in Japanese.

Systemic lupus erythematosus (SLE) is an autoimmune disease that causes multiple organ damage. Although recent genome-wide association studies (GWAS) have contributed to discovery of SLE susceptibility genes, few studies has been performed in Asian populations. Here, we report a GWAS for SLE examining 891 SLE cases and 3,384 controls and multi-stage replication studies examining 1,387 SLE cases and 28,564 controls in Japanese subjects. Considering that expression quantitative trait loci (eQTLs) have been implicated in genetic risks for autoimmune diseases, we integrated an eQTL study into the results of the GWAS. We observed enrichments of cis-eQTL positive loci among the known SLE susceptibility loci (30.8%) compared to the genome-wide SNPs (6.9%). In addition, we identified a novel association of a variant in the AF4/FMR2 family, member 1 (AFF1) gene at 4q21 with SLE susceptibility (rs340630; P = 8.3×10(-9), odds ratio = 1.21). The risk A allele of rs340630 demonstrated a cis-eQTL effect on the AFF1 transcript with enhanced expression levels (P<0.05). As AFF1 transcripts were prominently expressed in CD4(+) and CD19(+) peripheral blood lymphocytes, up-regulation of AFF1 may cause the abnormality in these lymphocytes, leading to disease onset.

Common variants at CDKAL1 and KLF9 are associated with body mass index in east Asian populations.

Obesity is a disorder with a complex genetic etiology, and its epidemic is a worldwide problem. Although multiple genetic loci associated with body mass index, the most common measure of obesity, have been identified in European populations, few studies have focused on Asian populations. Here we report a genome-wide association study and replication studies with 62,245 east Asian subjects, which identified two new body mass index-associated loci in the CDKAL1 locus at 6p22 (rs2206734, P = 1.4 × 10(-11)) and the KLF9 locus at 9q21 (rs11142387, P = 1.3 × 10(-9)), as well as several previously reported loci (the SEC16B, BDNF, FTO, MC4R and GIPR loci, P < 5.0 × 10(-8)). We subsequently performed gene-gene interaction analyses and identified an interaction (P = 2.0 × 10(-8)) between a SNP in the KLF9 locus (rs11142387) and one in the MSTN (also known as GDF8) locus at 2q32 (rs13034723). These findings should provide useful insights into the etiology of obesity.

Dose-adjustment study of tamoxifen based on CYP2D6 genotypes in Japanese breast cancer patients.

CYP2D6 is a key enzyme responsible for the metabolism of tamoxifen to active metabolites, endoxifen, and 4-hydroxytamoxifen. The breast cancer patients who are heterozygous and homozygous for decreased-function and null alleles of CYP2D6 showed lower plasma concentrations of endoxifen and 4-hydroxytamoxifen compared to patients with homozygous-wild-type allele, resulting in worse clinical outcome in tamoxifen therapy. We recruited 98 Japanese breast cancer patients, who had been taking 20 mg of tamoxifen daily as adjuvant setting. For the patients who have one or no normal allele of CYP2D6, dosages of tamoxifen were increased to 30 and 40 mg/day, respectively. The plasma concentrations of tamoxifen and its metabolites were measured at 8 weeks after dose-adjustment using liquid chromatography-tandem mass spectrometry. Association between tamoxifen dose and the incidence of adverse events during the tamoxifen treatment was investigated. In the patients with CYP2D6*1/*10 and CYP2D6*10/*10, the mean plasma endoxifen levels after dose increase were 1.4- and 1.7-fold higher, respectively, than those before the increase (P < 0.001). These plasma concentrations of endoxifen achieved similar level of those in the CYP2D6*1/*1 patients receiving 20 mg/day of tamoxifen. Plasma 4-hydroxytamoxifen concentrations in the patients with CYP2D6*1/*10 and CYP2D6*10/*10 were also significantly increased to the similar levels of the CYP2D6*1/*1 patients according to the increasing tamoxifen dosages (P < 0.001). The incidence of adverse events was not significantly different between before and after dose adjustment. This study provides the evidence that dose adjustment is useful for the patients carrying CYP2D6*10 allele to maintain the effective endoxifen level.

A genome-wide association study identifies locus at 10q22 associated with clinical outcomes of adjuvant tamoxifen therapy for breast cancer patients in Japanese.

Although many association studies of polymorphisms in candidate genes with the clinical outcomes of breast cancer patients receiving adjuvant tamoxifen therapy have been reported, genetic factors determining individual response to tamoxifen are not fully understood. To identify genetic polymorphisms associated with clinical outcomes of patients with tamoxifen treatment, we conducted a genome-wide association study (GWAS). We studied 462 Japanese patients with hormone receptor-positive, invasive breast cancer receiving adjuvant tamoxifen therapy. Of them, 240 patients were analyzed by genome-wide genotyping using the Illumina Human610-Quad BeadChips, and two independent sets of 105 and 117 cases were used for replication studies. In the GWAS, we detected significant associations with recurrence-free survival at 15 single-nucleotide polymorphisms (SNPs) on nine chromosomal loci (1p31, 1q41, 5q33, 7p11, 10q22, 12q13, 13q22, 18q12 and 19p13) that satisfied a genome-wide significant threshold (log-rank P= 2.87 × 10(-9)-9.41 × 10(-8)). Among them, rs10509373 in C10orf11 gene on 10q22 was significantly associated with recurrence-free survival in the replication study (log-rank P= 2.02 × 10(-4)) and a combined analysis indicated a strong association of this SNP with recurrence-free survival in breast cancer patients treated with tamoxifen (log-rank P= 1.26 × 10(-10)). Hazard ratio per C allele of rs10509373 was 4.51 [95% confidence interval (CI), 2.72-7.51; P= 6.29 × 10(-9)]. In a combined analysis of rs10509373 genotype with previously identified genetic makers, CYP2D6 and ABCC2, the number of risk alleles of these three genes had cumulative effects on recurrence-free survival among 345 patients receiving tamoxifen monotherapy (log-rank P= 2.28 × 10(-12)). In conclusion, we identified a novel locus associated with recurrence-free survival in Japanese breast cancer patients receiving adjuvant tamoxifen therapy.

PlatinumCNV: a Bayesian Gaussian mixture model for genotyping copy number polymorphisms using SNP array signal intensity data.

We present a statistical model for allele-specific patterns of copy number polymorphisms (CNPs) in commercial single nucleotide polymorphism (SNP) array data. This model is based on the observation that fluorescent signal intensities tend to cluster into clouds of similar allele-specific copy number (ASCN) genotypes at each SNP locus. To capture the tendency of this clustering to be made vague by instrumental errors, our model allows for cluster memberships to overlap each other, according to a Bayesian Gaussian mixture model (GMM). This approach is flexible, allowing for both absolute scale differences and X/Y scale imbalances of fluorescent signal intensities. The resulting model is also robust toward unobserved ASCN genotypes, which can be problematic for ordinary GMMs. We illustrated the utility of the model by applying it to commercial SNP array intensity data obtained from the Illumina HumanHap 610K platform. We retrieved more than 4,000 allele-specific CNPs, though 99% of them showed rather simple allele-specific CNP patterns with only a single aneuploid haplotype among the normal haplotypes. The genotyping accuracy was assessed by two approaches, quantitative PCR and replicated subjects. The results of both of these approaches demonstrated mean genotyping error rates of 1%. We demonstrated a preliminary genome-wide association study of three hematological traits. The result exhibited that it could form the foundation for new, more effective statistical methods for the mapping of both disease genes and quantitative trait loci with genome-wide CNPs. The methods described in this work are implemented in a software package, PlatinumCNV, available on the Internet.

A genome-wide association study identifies common variants near LBX1 associated with adolescent idiopathic scoliosis.

Adolescent idiopathic scoliosis is a pediatric spinal deformity affecting 2-3% of school-age children worldwide(1). Genetic factors have been implicated in its etiology(2). Through a genome-wide association study (GWAS) and replication study involving a total of 1,376 Japanese females with adolescent idiopathic scoliosis and 11,297 female controls, we identified a locus at chromosome 10q24.31 associated with adolescent idiopathic scoliosis susceptibility. The most significant SNP (rs11190870; combined P = 1.24 × 10(-19); odds ratio (OR) = 1.56) is located near LBX1 (encoding ladybird homeobox 1). The identification of this susceptibility locus provides new insights into the pathogenesis of adolescent idiopathic scoliosis.

Genome-wide association study identifies two susceptibility loci for exudative age-related macular degeneration in the Japanese population.

Age-related macular degeneration (AMD), the leading cause of irreversible blindness in the world, is a complex disease caused by multiple environmental and genetic risk factors. To identify genetic factors that modify the risk of exudative AMD in the Japanese population, we conducted a genome-wide association study and a replication study using a total of 1,536 individuals with exudative AMD and 18,894 controls. In addition to CFH (rs800292, P = 4.23 × 10(-15)) and ARMS2 (rs3750847, P = 8.67 × 10(-29)) loci, we identified two new susceptibility loci for exudative AMD: TNFRSF10A-LOC389641 on chromosome 8p21 (rs13278062, combined P = 1.03 × 10(-12), odds ratio = 0.73) and REST-C4orf14-POLR2B-IGFBP7 on chromosome 4q12 (rs1713985, combined P = 2.34 × 10(-8), odds ratio = 1.30). Fine mapping revealed that rs13278062, which is known to alter TNFRSF10A transcriptional activity, had the most significant association in 8p21 region. Our results provide new insights into the pathophysiology of exudative AMD.

Identification of nine novel loci associated with white blood cell subtypes in a Japanese population.

White blood cells (WBCs) mediate immune systems and consist of various subtypes with distinct roles. Elucidation of the mechanism that regulates the counts of the WBC subtypes would provide useful insights into both the etiology of the immune system and disease pathogenesis. In this study, we report results of genome-wide association studies (GWAS) and a replication study for the counts of the 5 main WBC subtypes (neutrophils, lymphocytes, monocytes, basophils, and eosinophils) using 14,792 Japanese subjects enrolled in the BioBank Japan Project. We identified 12 significantly associated loci that satisfied the genome-wide significance threshold of P<5.0×10(-8), of which 9 loci were novel (the CDK6 locus for the neutrophil count; the ITGA4, MLZE, STXBP6 loci, and the MHC region for the monocyte count; the SLC45A3-NUCKS1, GATA2, NAALAD2, ERG loci for the basophil count). We further evaluated associations in the identified loci using 15,600 subjects from Caucasian populations. These WBC subtype-related loci demonstrated a variety of patterns of pleiotropic associations within the WBC subtypes, or with total WBC count, platelet count, or red blood cell-related traits (n = 30,454), which suggests unique and common functional roles of these loci in the processes of hematopoiesis. This study should contribute to the understanding of the genetic backgrounds of the WBC subtypes and hematological traits.

A genome-wide association study of chronic hepatitis B identified novel risk locus in a Japanese population.

Hepatitis B virus (HBV) infection is a major health issue worldwide which may lead to hepatic dysfunction, liver cirrhosis and hepatocellular carcinoma. To identify host genetic factors that are associated with chronic hepatitis B (CHB) susceptibility, we previously conducted a two-stage genome-wide association study (GWAS) and identified the association of HLA-DP variants with CHB in Asians; however, only 179 cases and 934 controls were genotyped using genome-wide single nucleotide polymorphism (SNP) arrays. Here, we performed a second GWAS of 519 747 SNPs in 458 Japanese CHB cases and 2056 controls. After adjustment with the previously identified variants in the HLA-DP locus (rs9277535), we detected strong associations at 16 loci with P-value of <5 × 10(-5). We analyzed these loci in three independent Japanese cohorts (2209 CHB cases and 4440 controls) and found significant association of two SNPs (rs2856718 and rs7453920) within the HLA-DQ locus (overall P-value of 5.98 × 10(-28) and 3.99 × 10(-37)). Association of CHB with SNPs rs2856718 and rs7453920 remains significant even after stratification with rs3077 and rs9277535, indicating independent effect of HLA-DQ variants on CHB susceptibility (P-value of 1.52 × 10(-21)- 2.38 × 10(-30)). Subsequent analyses revealed DQA1*0102-DQB1*0604 and DQA1*0101-DQB1*0501 [odds ratios (OR) =0.16, and 0.39, respectively] as protective haplotypes and DQA1*0102-DQB1*0303 and DQA1*0301-DQB1*0601 (OR = 19.03 and 5.02, respectively) as risk haplotypes. These findings indicated that variants in antigen-binding regions of HLA-DP and HLA-DQ contribute to the risk of persistent HBV infection.