Identification of Loci at 1q21 and 16q23 That Affect Susceptibility to Inflammatory Bowel Disease in Koreans.

Recent genome-wide association studies have identified more than 200 regions that affect susceptibility to inflammatory bowel disease (IBD). However, identified common variants account for only a fraction of IBD heritability and largely have been identified in populations of European ancestry. We performed a genome-wide association study of susceptibility loci in Korean individuals, comprising a total of 1505 IBD patients and 4041 controls. We identified 2 new susceptibility loci for IBD at genome-wide significance: rs3766920 near PYGO2-SHC1 at 1q21 and rs16953946 in CDYL2 at 16q23. In addition, we confirmed associations, in Koreans, with 28 established IBD loci (P < 2.16 × 10-4). Our findings support the complementary value of genetic studies in different populations.

An association study of TOLL and CARD with leprosy susceptibility in Chinese population.

Previous genome-wide association studies (GWASs) identified multiple susceptibility loci that have highlighted the important role of TLR (Toll-like receptor) and CARD (caspase recruitment domain) genes in leprosy. A large three-stage candidate gene-based association study of 30 TLR and 47 CARD genes was performed in the leprosy samples of Chinese Han. Of 4363 SNPs investigated, eight SNPs showed suggestive association (P < 0.01) in our previously published GWAS datasets (Stage 1). Of the eight SNPs, rs2735591 and rs4889841 showed significant association (P < 0.001) in an independent series of 1504 cases and 1502 controls (Stage 2), but only rs2735591 (next to BCL10) showed significant association in the second independent series of 938 cases and 5827 controls (Stage 3). Rs2735591 showed consistent association across the three stages (P > 0.05 for heterogeneity test), significant association in the combined validation samples (Pcorrected = 5.54 × 10(-4) after correction for 4363 SNPs tested) and genome-wide significance in the whole GWAS and validation samples (P = 1.03 × 10(-9), OR = 1.24). In addition, we demonstrated the lower expression of BCL10 in leprosy lesions than normal skins and a significant gene connection between BCL10 and the eight previously identified leprosy loci that are associated with NFκB, a major regulator of downstream inflammatory responses, which provides further biological evidence for the association. We have discovered a novel susceptibility locus on 1p22, which implicates BCL10 as a new susceptibility gene for leprosy. Our finding highlights the important role of both innate and adaptive immune responses in leprosy.

Identification of IL18RAP/IL18R1 and IL12B as leprosy risk genes demonstrates shared pathogenesis between inflammation and infectious diseases.

Of eight leprosy susceptibility loci identified by genome-wide association studies, five have been implicated in Crohn disease, suggesting a common genetic fingerprint between leprosy and inflammatory bowel disease (IBD). Here, we conducted a multiple-stage genetic association study of 133 IBD susceptibility loci in multiple leprosy samples (totaling 4,971 leprosy cases and 5,503 controls) from a Chinese population and discovered two associations at rs2058660 on 2q12.1 (p = 4.57 × 10(-19); odds ratio [OR] = 1.30) and rs6871626 on 5q33.3 (p = 3.95 × 10(-18); OR = 0.75), implicating IL18RAP/IL18R1 and IL12B as susceptibility genes for leprosy. Our study reveals the important role of IL12/IL18-mediated transcriptional regulation of IFN-γ production in leprosy, and together with previous findings, it demonstrates the shared genetic susceptibility between infectious and inflammatory diseases.

Genomewide association study of leprosy.

BACKGROUND: The narrow host range of Mycobacterium leprae and the fact that it is refractory to growth in culture has limited research on and the biologic understanding of leprosy. Host genetic factors are thought to influence susceptibility to infection as well as disease progression. METHODS: We performed a two-stage genomewide association study by genotyping 706 patients and 1225 controls using the Human610-Quad BeadChip (Illumina). We then tested three independent replication sets for an association between the presence of leprosy and 93 single-nucleotide polymorphisms (SNPs) that were most strongly associated with the disease in the genomewide association study. Together, these replication sets comprised 3254 patients and 5955 controls. We also carried out tests of heterogeneity of the associations (or lack thereof) between these 93 SNPs and disease, stratified according to clinical subtype (multibacillary vs. paucibacillary). RESULTS: We observed a significant association (P<1.00x10(-10)) between SNPs in the genes CCDC122, C13orf31, NOD2, TNFSF15, HLA-DR, and RIPK2 and a trend toward an association (P=5.10x10(-5)) with a SNP in LRRK2. The associations between the SNPs in C13orf31, LRRK2, NOD2, and RIPK2 and multibacillary leprosy were stronger than the associations between these SNPs and paucibacillary leprosy. CONCLUSIONS: Variants of genes in the NOD2-mediated signaling pathway (which regulates the innate immune response) are associated with susceptibility to infection with M. leprae.

Identification of Ten Additional Susceptibility Loci for Ulcerative Colitis Through Immunochip Analysis in Koreans.

BACKGROUND: Recent genetic association studies identified more than 160 susceptibility loci for inflammatory bowel disease in Caucasian populations, but studies in Asian populations are limited. We have previously reported 3 loci associated with Korean ulcerative colitis (UC). METHODS: Using the Immunochip custom single nucleotide polymorphisms (SNP) array designed for dense genotyping of 186 known disease loci from 12 immune-mediated diseases, we analyzed 705 patients with UC and 1178 controls for 536,821 SNPs (89,057 genotyped and 447,764 imputed) in the discovery stage followed by replication in additional 980 affected individuals and 2694 controls in a Korean population. RESULTS: We confirmed the associations of 10 known UC risk loci in Koreans: rs76418789 in IL23R (combined P = 1.25 × 10), rs4728142 in IRF5 (combined P = 3.17 × 10), rs1830610 near JAK2 (combined P = 2.28 × 10), rs1555791 near TNFRSF14 (combined P = 1.62 × 10), rs880790 between IL10-IL19 (combined P = 3.73 × 10), rs10185424 between IL1R2-IL1R1 (combined P = 1.54 × 10), rs6478108 in TNFSF15 (combined P = 9.28 × 10), rs861857 between UBE2L3-YDJC (combined P = 3.05 × 10), rs1801274 in FCGR2A (discovery P = 1.54 × 10), and rs17085007 between GPR12-USP12 (discovery P = 3.64 × 10). The percentage of phenotype variance explained by the 13 risk loci (including 3 previously reported loci) was 5.61% in Koreans (on the liability scale, population prevalence = 0.0308%). CONCLUSIONS: Our study increased the number of UC susceptibility loci in Koreans to 13 and highlighted the extensive sharing of genetic risk across populations of UC.

Discovery of six new susceptibility loci and analysis of pleiotropic effects in leprosy.

Genome-wide association studies (GWAS) have led to the discovery of several susceptibility loci for leprosy with robust evidence, providing biological insight into the role of host genetic factors in mycobacterial infection. However, the identified loci only partially explain disease heritability, and additional genetic risk factors remain to be discovered. We performed a 3-stage GWAS of leprosy in the Chinese population using 8,313 cases and 16,017 controls. Besides confirming all previously published loci, we discovered six new susceptibility loci, and further gene prioritization analysis of these loci implicated BATF3, CCDC88B and CIITA-SOCS1 as new susceptibility genes for leprosy. A systematic evaluation of pleiotropic effects demonstrated a high tendency for leprosy susceptibility loci to show association with autoimmunity and inflammatory diseases. Further analysis suggests that molecular sensing of infection might have a similar pathogenic role across these diseases, whereas immune responses have discordant roles in infectious and inflammatory diseases.

Identification of new susceptibility loci for IgA nephropathy in Han Chinese.

IgA nephropathy (IgAN) is one of the most common primary glomerulonephritis. Previously identified genome-wide association study (GWAS) loci explain only a fraction of disease risk. To identify novel susceptibility loci in Han Chinese, we conduct a four-stage GWAS comprising 8,313 cases and 19,680 controls. Here, we show novel associations at ST6GAL1 on 3q27.3 (rs7634389, odds ratio (OR)=1.13, P=7.27 × 10(-10)), ACCS on 11p11.2 (rs2074038, OR=1.14, P=3.93 × 10(-9)) and ODF1-KLF10 on 8q22.3 (rs2033562, OR=1.13, P=1.41 × 10(-9)), validate a recently reported association at ITGAX-ITGAM on 16p11.2 (rs7190997, OR=1.22, P=2.26 × 10(-19)), and identify three independent signals within the DEFA locus (rs2738058, P=1.15 × 10(-19); rs12716641, P=9.53 × 10(-9); rs9314614, P=4.25 × 10(-9), multivariate association). The risk variants on 3q27.3 and 11p11.2 show strong association with mRNA expression levels in blood cells while allele frequencies of the risk variants within ST6GAL1, ACCS and DEFA correlate with geographical variation in IgAN prevalence. Our findings expand our understanding on IgAN genetic susceptibility and provide novel biological insights into molecular mechanisms underlying IgAN.

Combined linkage and family-based association analysis improves candidate gene detection in Genetic Analysis Workshop 18 simulation data.

Because the genotype-phenotype correlation information is investigated differently by linkage and association analyses, various efforts have been made to model linkage and association jointly. However, joint modeling methods are usually computationally intensive; hence they cannot currently accommodate large pedigrees with dense markers. This article proposes a simple method to combine the linkage and association evidence with the aim of improving the detection power of disease susceptibility genes. Our detection power comparisons show that the combined linkage-association p values can improve remarkably the causal gene detection power in Genetic Analysis Workshop 18 simulation data.

Genome-wide linkage, exome sequencing and functional analyses identify ABCB6 as the pathogenic gene of dyschromatosis universalis hereditaria.

BACKGROUND: As a genetic disorder of abnormal pigmentation, the molecular basis of dyschromatosis universalis hereditaria (DUH) had remained unclear until recently when ABCB6 was reported as a causative gene of DUH. METHODOLOGY: We performed genome-wide linkage scan using Illumina Human 660W-Quad BeadChip and exome sequencing analyses using Agilent SureSelect Human All Exon Kits in a multiplex Chinese DUH family to identify the pathogenic mutations and verified the candidate mutations using Sanger sequencing. Quantitative RT-PCR and Immunohistochemistry was performed to verify the expression of the pathogenic gene, Zebrafish was also used to confirm the functional role of ABCB6 in melanocytes and pigmentation. RESULTS: Genome-wide linkage (assuming autosomal dominant inheritance mode) and exome sequencing analyses identified ABCB6 as the disease candidate gene by discovering a coding mutation (c.1358C>T; p.Ala453Val) that co-segregates with the disease phenotype. Further mutation analysis of ABCB6 in four other DUH families and two sporadic cases by Sanger sequencing confirmed the mutation (c.1358C>T; p.Ala453Val) and discovered a second, co-segregating coding mutation (c.964A>C; p.Ser322Lys) in one of the four families. Both mutations were heterozygous in DUH patients and not present in the 1000 Genome Project and dbSNP database as well as 1,516 unrelated Chinese healthy controls. Expression analysis in human skin and mutagenesis interrogation in zebrafish confirmed the functional role of ABCB6 in melanocytes and pigmentation. Given the involvement of ABCB6 mutations in coloboma, we performed ophthalmological examination of the DUH carriers of ABCB6 mutations and found ocular abnormalities in them. CONCLUSION: Our study has advanced our understanding of DUH pathogenesis and revealed the shared pathological mechanism between pigmentary DUH and ocular coloboma.

Chromosome 2p14 is linked to susceptibility to leprosy.

BACKGROUND: A genetic component to the etiology of leprosy is well recognized but the mechanism of inheritance and the genes involved are yet to be fully established. METHODOLOGY: A genome-wide single nucleotide polymorphism (SNP) based linkage analysis was carried out using 23 pedigrees, each with 3 to 7 family members affected by leprosy. Multipoint parametric and non-parametric linkage analyses were performed using MERLIN 1.1.1. PRINCIPAL FINDINGS: Genome-wide significant evidence for linkage was identified on chromosome 2p14 with a heterogeneity logarithm of odds (HLOD) score of 3.51 (rs1106577) under a recessive model of inheritance, while suggestive evidence was identified on chr.4q22 (HLOD 2.92, rs1349350, dominant model), chr. 8q24 (HLOD 2.74, rs1618523, recessive model) and chr.16q24 (HLOD 1.93, rs276990 dominant model). Our study also provided moderate evidence for a linkage locus on chromosome 6q24-26 by non-parametric linkage analysis (rs6570858, LOD 1.54, p = 0.004), overlapping a previously reported linkage region on chromosome 6q25-26. CONCLUSION: A genome-wide linkage analysis has identified a new linkage locus on chromosome 2p14 for leprosy in Pedigrees from China.

ARVCF genetic influences on neurocognitive and neuroanatomical intermediate phenotypes in Chinese patients with schizophrenia.

OBJECTIVE: There are notable similarities between velocardiofacial syndrome and schizophrenia in terms of neurocognitive deficits and brain structural abnormalities. These similarities have supported the role of the armadillo repeat gene deleted in velocardiofacial syndrome (ARVCF) as a susceptibility gene in schizophrenia. This study investigated the relationships between haplotypes of the ARVCF gene and specific intermediate phenotypes in schizophrenia. We hypothesized that ARVCF gene haplotypes influence caudate nucleus volume, fractional anisotropy, and neurocognitive functioning in schizophrenia. METHOD: Between May 2006 and November 2009, 200 Chinese participants (125 patients with DSM-IV diagnosis of schizophrenia and 75 controls) were genotyped using blood samples, and a subset of 166 participants (99 patients with DSM-IV diagnosis of schizophrenia and 67 controls) underwent structural magnetic resonance imaging, diffusion tensor imaging, and completed neuropsychological testing. RESULTS: The haplotype T-G-A-T-T-G-G-C-T-G-T (ARVCF-Hap1) was significantly associated with fractional anisotropy of the caudate nucleus and executive functioning in patients. Specifically, patients with more copies of ARVCF-Hap1 have lower white matter integrity in caudate nucleus (P = .0008) and greater perseverative errors (P = .00003) on the Wisconsin Card Sorting Test. A trend of lower caudate volume (P = .015) in patients with more copies of ARVCF-Hap1 was also observed. CONCLUSIONS: These findings are consistent with known ARVCF gene effects on neurodevelopment in terms of cellular arrangement, migration, and intracellular signaling involving the striatum and may involve interactions with other brain networks such as prefrontal cortex, and they underscore the importance of imaging-genetic studies to elucidate the genetic influences underlying intermediate phenotypes in complex neurobehavioral disorders.

Genetic variants in ER cofactor genes and endometrial cancer risk.

Given that the transcriptional regulatory activity of estrogen receptor (ER) is modulated by its biochemical cofactors, genetic variation within the ER cofactor genes may alter cellular response to estrogen exposure and consequently modify the risk for endometrial cancer. We genotyped 685 tagging SNPs within 60 ER cofactor genes in 564 endometrial cancer cases and 1,510 controls from Sweden, and tested their associations with the risk of endometrial cancer. We investigated the associations of individual SNPs by using a trend test as well as multiple SNPs within a gene or gene complex by using multi-variant association analysis. No significant association was observed for any individual SNPs or genes, but a marginal association of the cumulative genetic variation of the NCOA2 complex as a whole (NCOA2, CARM1, CREBBP, PRMT1 and EP300) with endometrial cancer risk was observed (P(adjusted) = 0.033). However, the association failed to be replicated in an independent European dataset of 1265 cases and 5190 controls (P = 0.71). The results indicate that common genetic variants within ER cofactor genes are unlikely to play a significant role in endometrial cancer risk in European population.

A comprehensive association analysis of homocysteine metabolic pathway genes in Singaporean Chinese with ischemic stroke.

BACKGROUND: The effect of genetic factors, apart from 5,10-methylenetetrahydrofolate reductase (MTHFR) polymorphisms, on elevated plasma homocysteine levels and increasing ischemic stroke risk have not been fully elucidated. We conducted a comprehensive analysis of 25 genes involved in homocysteine metabolism to investigate association of common variants within these genes with ischemic stroke risk. METHODOLOGY/PRINCIPAL FINDINGS: The study was done in two stages. In the initial study, SNP and haplotype-based association analyses were performed using 147 tagging Single Nucleotide Polymorphisms (SNPs) in 360 stroke patients and 354 non-stroke controls of Singaporean Chinese ethnicity. Joint association analysis of significant SNPs was then performed to assess the cumulative effect of these variants on ischemic stroke risk. In the replication study, 8 SNPs were selected for validation in an independent set of 420 matched case-control pairs of Singaporean Chinese ethnicity. SNP analysis from the initial study suggested 3 risk variants in the MTRR, SHMT1 and TCN2 genes which were moderately associated with ischemic stroke risk, independent of known stroke risk factors. Although the replication study failed to support single-SNP associations observed in the initial study, joint association analysis of the 3 variants in combined initial and replication samples revealed a trend of elevated risk with an increased number of risk alleles (Joint P(trend) = 1.2×10(-6)). CONCLUSIONS: Our study did not find direct evidence of associations between any single polymorphisms of homocysteine metabolic pathway genes and ischemic stroke, but suggests that the cumulative effect of several small to moderate risk variants from genes involved in homocysteine metabolism may jointly confer a significant impact on ischemic stroke risk.

Identification of two new loci at IL23R and RAB32 that influence susceptibility to leprosy.

We performed a genome-wide association study with 706 individuals with leprosy and 5,581 control individuals and replicated the top 24 SNPs in three independent replication samples, including a total of 3,301 individuals with leprosy and 5,299 control individuals from China. Two loci not previously associated with the disease were identified with genome-wide significance: rs2275606 (combined P = 3.94 × 10(-14), OR = 1.30) on 6q24.3 and rs3762318 (combined P = 3.27 × 10(-11), OR = 0.69) on 1p31.3. These associations implicate IL23R and RAB32 as new susceptibility genes for leprosy. Furthermore, we identified evidence of interaction between the NOD2 and RIPK2 loci, which is consistent with the biological association of the proteins encoded by these genes (NOD2-RIPK2 complex) in activating the NF-κB pathway as a part of the host defense response to infection. Our findings have expanded the biological functions of IL23R by uncovering its involvement in infectious disease susceptibility and suggest a potential involvement of autophagocytosis in leprosy pathogenesis. The IL23R association supports previous observations of the marked overlap of susceptibility genes for leprosy and Crohn's disease, implying common pathogenesis mechanisms.

A genome-wide association study in Han Chinese identifies multiple susceptibility loci for IgA nephropathy.

We performed a two-stage genome-wide association study of IgA nephropathy (IgAN) in Han Chinese, with 1,434 affected individuals (cases) and 4,270 controls in the discovery phase and follow-up of the top 61 SNPs in an additional 2,703 cases and 3,464 controls. We identified associations at 17p13 (rs3803800, P = 9.40 × 10(-11), OR = 1.21; rs4227, P = 4.31 × 10(-10), OR = 1.23) and 8p23 (rs2738048, P = 3.18 × 10(-14), OR = 0.79) that implicated the genes encoding tumor necrosis factor (TNFSF13) and α-defensin (DEFA) as susceptibility genes. In addition, we found multiple associations in the major histocompatibility complex (MHC) region (rs660895, P = 4.13 × 10(-20), OR = 1.34; rs1794275, P = 3.43 × 10(-13), OR = 1.30; rs2523946, P = 1.74 × 10(-11), OR = 1.21) and confirmed a previously reported association at 22q12 (rs12537, P = 1.17 × 10(-11), OR = 0.78). We also found that rs660895 was associated with clinical subtypes of IgAN (P = 0.003), proteinuria (P = 0.025) and IgA levels (P = 0.047). Our findings show that IgAN is associated with variants near genes involved in innate immunity and inflammation.

A genome-wide association study in Han Chinese identifies new susceptibility loci for ankylosing spondylitis.

To identify susceptibility loci for ankylosing spondylitis, we performed a two-stage genome-wide association study in Han Chinese. In the discovery stage, we analyzed 1,356,350 autosomal SNPs in 1,837 individuals with ankylosing spondylitis and 4,231 controls; in the validation stage, we analyzed 30 suggestive SNPs in an additional 2,100 affected individuals and 3,496 controls. We identified two new susceptibility loci between EDIL3 and HAPLN1 at 5q14.3 (rs4552569; P = 8.77 × 10(-10)) and within ANO6 at 12q12 (rs17095830; P = 1.63 × 10(-8)). We also confirmed previously reported associations in Europeans within the major histocompatibility complex (MHC) region (top SNP, rs13202464; P < 5 × 10(-324)) and at 2p15 (rs10865331; P = 1.98 × 10(-8)). We show that rs13202464 within the MHC region mainly represents the risk effect of HLA-B*27 variants (including HLA-B*2704, HLA-B*2705 and HLA-B*2715) in Chinese. The two newly discovered loci implicate genes related to bone formation and cartilage development, suggesting their potential involvement in the etiology of ankylosing spondylitis.

A genome-wide association study of nasopharyngeal carcinoma identifies three new susceptibility loci.

To identify genetic susceptibility loci for nasopharyngeal carcinoma (NPC), a genome-wide association study was performed using 464,328 autosomal SNPs in 1,583 NPC affected individuals (cases) and 1,894 controls of southern Chinese descent. The top 49 SNPs from the genome-wide association study were genotyped in 3,507 cases and 3,063 controls of southern Chinese descent from Guangdong and Guangxi. The seven supportive SNPs were further confirmed by transmission disequilibrium test analysis in 279 trios from Guangdong. We identified three new susceptibility loci, TNFRSF19 on 13q12 (rs9510787, Pcombined=1.53x10(-9), odds ratio (OR)=1.20), MDS1-EVI1 on 3q26 (rs6774494, Pcombined=1.34x10(-8), OR=0.84) and the CDKN2A-CDKN2B gene cluster on 9p21 (rs1412829, Pcombined=4.84x10(-7), OR=0.78). Furthermore, we confirmed the role of HLA by revealing independent associations at rs2860580 (Pcombined=4.88x10(-67), OR=0.58), rs2894207 (Pcombined=3.42x10(-33), OR=0.61) and rs28421666 (Pcombined=2.49x10(-18), OR=0.67). Our findings provide new insights into the pathogenesis of NPC by highlighting the involvement of pathways related to TNFRSF19 and MDS1-EVI1 in addition to HLA molecules.

Genome-wide association study of follicular lymphoma identifies a risk locus at 6p21.32.

To identify susceptibility loci for non-Hodgkin lymphoma subtypes, we conducted a three-stage genome-wide association study. We identified two variants associated with follicular lymphoma at 6p21.32 (rs10484561, combined P = 1.12 x 10(-29) and rs7755224, combined P = 2.00 x 10(-19); r(2) = 1.0), supporting the idea that major histocompatibility complex genetic variation influences follicular lymphoma susceptibility. We also found confirmatory evidence of a previously reported association between chronic lymphocytic leukemia/small lymphocytic lymphoma and rs735665 (combined P = 4.24 x 10(-9)).