Linkage disequilibrium of a type 1 diabetes susceptibility locus with a regulatory IL12B allele

Type 1 diabetes (T1D; or insulin-dependent diabetes mellitus, IDDM) is an autoimmune disease with both genetic and environmental components. In addition to the human leukocyte antigen (HLA) complex, the single major genetic contributor of susceptibility, an unknown number of other unidentified genes are required to mediate disease. Although many loci conferring susceptibility to T1D have been mapped, their identification has proven problematic due to the complex nature of this disease. Our strategy for finding T1D susceptibility genes has been to test for human homologues of loci implicated in diabetes-prone NOD (non-obese diabetic) mice, together with application of biologically relevant stratification methods. We report here a new susceptibility locus, IDDM18, located near the interleukin-12 (IL-12)p40 gene, IL12B. Significant bias in transmission of IL12B alleles was observed in affected sibpairs and was confirmed in an independent cohort of simplex families. A single base change in the 3' UTR showed strong linkage disequilibrium with the T1D susceptibility locus. The IL12B 3' UTR alleles showed different levels of expression in cell lines. Variation in IL-12p40 production may influence T-cell responses crucial for either mediating or protecting against this and other autoimmune diseases.     

Characterization of single-nucleotide polymorphisms in coding regions of human genes.

A major goal in human genetics is to understand the role of common genetic variants in susceptibility to common diseases. This will require characterizing the nature of gene variation in human populations, assembling an extensive catalogue of single-nucleotide polymorphisms (SNPs) in candidate genes and performing association studies for particular diseases. At present, our knowledge of human gene variation remains rudimentary. Here we describe a systematic survey of SNPs in the coding regions of human genes. We identified SNPs in 106 genes relevant to cardiovascular disease, endocrinology and neuropsychiatry by screening an average of 114 independent alleles using 2 independent screening methods. To ensure high accuracy, all reported SNPs were confirmed by DNA sequencing. We identified 560 SNPs, including 392 coding-region SNPs (cSNPs) divided roughly equally between those causing synonymous and non-synonymous changes. We observed different rates of polymorphism among classes of sites within genes (non-coding, degenerate and non-degenerate) as well as between genes. The cSNPs most likely to influence disease, those that alter the amino acid sequence of the encoded protein, are found at a lower rate and with lower allele frequencies than silent substitutions. This likely reflects selection acting against deleterious alleles during human evolution. The lower allele frequency of missense cSNPs has implications for the compilation of a comprehensive catalogue, as well as for the subsequent application to disease association.     

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.     

Genome-wide association study and mouse model identify interaction between RET and EDNRB pathways in Hirschsprung disease

Genetic studies of Hirschsprung disease, a common congenital malformation, have identified eight genes with mutations that can be associated with this condition. Mutations at individual loci are, however, neither necessary nor sufficient to cause clinical disease. We conducted a genome-wide association study in 43 Mennonite family trios using 2,083 microsatellites and single-nucleotide polymorphisms and a new multipoint linkage disequilibrium method that searches for association arising from common ancestry. We identified susceptibility loci at 10q11, 13q22 and 16q23; the gene at 13q22 is EDNRB, encoding a G protein-coupled receptor (GPCR) and the gene at 10q11 is RET, encoding a receptor tyrosine kinase (RTK). Statistically significant joint transmission of RET and EDNRB alleles in affected individuals and non-complementation of aganglionosis in mouse intercrosses between Ret null and the Ednrb hypomorphic piebald allele are suggestive of epistasis between EDNRB and RET. Thus, genetic interaction between mutations in RET and EDNRB is an underlying mechanism for this complex disorder.     

Two independent alleles at 6q23 associated with risk of rheumatoid arthritis

To identify susceptibility alleles associated with rheumatoid arthritis, we genotyped 397 individuals with rheumatoid arthritis for 116,204 SNPs and carried out an association analysis in comparison to publicly available genotype data for 1,211 related individuals from the Framingham Heart Study. After evaluating and adjusting for technical and population biases, we identified a SNP at 6q23 (rs10499194, approximately 150 kb from TNFAIP3 and OLIG3) that was reproducibly associated with rheumatoid arthritis both in the genome-wide association (GWA) scan and in 5,541 additional case-control samples (P = 10(-3), GWA scan; P < 10(-6), replication; P = 10(-9), combined). In a concurrent study, the Wellcome Trust Case Control Consortium (WTCCC) has reported strong association of rheumatoid arthritis susceptibility to a different SNP located 3.8 kb from rs10499194 (rs6920220; P = 5 x 10(-6) in WTCCC). We show that these two SNP associations are statistically independent, are each reproducible in the comparison of our data and WTCCC data, and define risk and protective haplotypes for rheumatoid arthritis at 6q23.     

Systematic meta-analyses and field synopsis of genetic association studies in schizophrenia: the SzGene database

In an effort to pinpoint potential genetic risk factors for schizophrenia, research groups worldwide have published over 1,000 genetic association studies with largely inconsistent results. To facilitate the interpretation of these findings, we have created a regularly updated online database of all published genetic association studies for schizophrenia ('SzGene'). For all polymorphisms having genotype data available in at least four independent case-control samples, we systematically carried out random-effects meta-analyses using allelic contrasts. Across 118 meta-analyses, a total of 24 genetic variants in 16 different genes (APOE, COMT, DAO, DRD1, DRD2, DRD4, DTNBP1, GABRB2, GRIN2B, HP, IL1B, MTHFR, PLXNA2, SLC6A4, TP53 and TPH1) showed nominally significant effects with average summary odds ratios of approximately 1.23. Seven of these variants had not been previously meta-analyzed. According to recently proposed criteria for the assessment of cumulative evidence in genetic association studies, four of the significant results can be characterized as showing 'strong' epidemiological credibility. Our project represents the first comprehensive online resource for systematically synthesized and graded evidence of genetic association studies in schizophrenia. As such, it could serve as a model for field synopses of genetic associations in other common and genetically complex disorders.     

Resequencing of positional candidates identifies low frequency IL23R coding variants protecting against inflammatory bowel disease

Genome-wide association studies (GWAS) have identified dozens of risk loci for many complex disorders, including Crohn's disease. However, common disease-associated SNPs explain at most ～20% of the genetic variance for Crohn's disease. Several factors may account for this unexplained heritability, including rare risk variants not adequately tagged thus far in GWAS. That rare susceptibility variants indeed contribute to variation in multifactorial phenotypes has been demonstrated for colorectal cancer, plasma high-density lipoprotein cholesterol levels, blood pressure, type 1 diabetes, hypertriglyceridemia and, in the case of Crohn's disease, for NOD2 (refs. 14,15). Here we describe the use of high-throughput resequencing of DNA pools to search for rare coding variants influencing susceptibility to Crohn's disease in 63 GWAS-identified positional candidate genes. We identify low frequency coding variants conferring protection against inflammatory bowel disease in IL23R, but we conclude that rare coding variants in positional candidates do not make a large contribution to inherited predisposition to Crohn's disease.     

Genome-wide association analysis identifies 20 loci that influence adult height

Adult height is a model polygenic trait, but there has been limited success in identifying the genes underlying its normal variation. To identify genetic variants influencing adult human height, we used genome-wide association data from 13,665 individuals and genotyped 39 variants in an additional 16,482 samples. We identified 20 variants associated with adult height (P < 5 x 10(-7), with 10 reaching P < 1 x 10(-10)). Combined, the 20 SNPs explain approximately 3% of height variation, with a approximately 5 cm difference between the 6.2% of people with 17 or fewer 'tall' alleles compared to the 5.5% with 27 or more 'tall' alleles. The loci we identified implicate genes in Hedgehog signaling (IHH, HHIP, PTCH1), extracellular matrix (EFEMP1, ADAMTSL3, ACAN) and cancer (CDK6, HMGA2, DLEU7) pathways, and provide new insights into human growth and developmental processes. Finally, our results provide insights into the genetic architecture of a classic quantitative trait.     

Genetic variation in the 5q31 cytokine gene cluster confers susceptibility to Crohn disease

Linkage disequilibrium (LD) mapping provides a powerful method for fine-structure localization of rare disease genes, but has not yet been widely applied to common disease. We sought to design a systematic approach for LD mapping and apply it to the localization of a gene (IBD5) conferring susceptibility to Crohn disease. The key issues are: (i) to detect a significant LD signal (ii) to rigorously bound the critical region and (iii) to identify the causal genetic variant within this region. We previously mapped the IBD5 locus to a large region spanning 18 cM of chromosome 5q31 (P<10(-4)). Using dense genetic maps of microsatellite markers and single-nucleotide polymorphisms (SNPs) across the entire region, we found strong evidence of LD. We bound the region to a common haplotype spanning 250 kb that shows strong association with the disease (P< 2 x 10(-7)) and contains the cytokine gene cluster. This finding provides overwhelming evidence that a specific common haplotype of the cytokine region in 5q31 confers susceptibility to Crohn disease. However, genetic evidence alone is not sufficient to identify the causal mutation within this region, as strong LD across the region results in multiple SNPs having equivalent genetic evidence-each consistent with the expected properties of the IBD5 locus. These results have important implications for Crohn disease in particular and LD mapping in general.     

A predominant role for the HLA class II region in the association of the MHC region with multiple sclerosis

Genetic susceptibility to multiple sclerosis is associated with genes of the major histocompatibility complex (MHC), particularly HLA-DRB1 and HLA-DQB1 (ref. 1). Both locus and allelic heterogeneity have been reported in this genomic region. To clarify whether HLA-DRB1 itself, nearby genes in the region encoding the MHC or combinations of these loci underlie susceptibility to multiple sclerosis, we genotyped 1,185 Canadian and Finnish families with multiple sclerosis (n = 4,203 individuals) with a high-density SNP panel spanning the genes encoding the MHC and flanking genomic regions. Strong associations in Canadian and Finnish samples were observed with blocks in the HLA class II genomic region (P < 4.9 x 10(-13) and P < 2.0 x 10(-16), respectively), but the strongest association was with HLA-DRB1 (P < 4.4 x 10(-17)). Conditioning on either HLA-DRB1 or the most significant HLA class II haplotype block found no additional block or SNP association independent of the HLA class II genomic region. This study therefore indicates that MHC-associated susceptibility to multiple sclerosis is determined by HLA class II alleles, their interactions and closely neighboring variants.     

Genetics of gene expression in primary immune cells identifies cell type-specific master regulators and roles of HLA alleles

Trans-acting genetic variants have a substantial, albeit poorly characterized, role in the heritable determination of gene expression. Using paired purified primary monocytes and B cells, we identify new predominantly cell type-specific cis and trans expression quantitative trait loci (eQTLs), including multi-locus trans associations to LYZ and KLF4 in monocytes and B cells, respectively. Additionally, we observe a B cell-specific trans association of rs11171739 at 12q13.2, a known autoimmune disease locus, with IP6K2 (P = 5.8 × 10(-15)), PRIC285 (P = 3.0 × 10(-10)) and an upstream region of CDKN1A (P = 2 × 10(-52)), suggesting roles for cell cycle regulation and peroxisome proliferator-activated receptor γ (PPARγ) signaling in autoimmune pathogenesis. We also find that specific human leukocyte antigen (HLA) alleles form trans associations with the expression of AOAH and ARHGAP24 in monocytes but not in B cells. In summary, we show that mapping gene expression in defined primary cell populations identifies new cell type-specific trans-regulated networks and provides insights into the genetic basis of disease susceptibility.     

Chromosome 6q25 is linked to susceptibility to leprosy in a Vietnamese population

Leprosy, a chronic infectious disease caused by Mycobacterium leprae, affects an estimated 700,000 persons each year. Clinically, leprosy can be categorized as paucibacillary or multibacillary disease. These clinical forms develop in persons that are intrinsically susceptible to leprosy per se, that is, leprosy independent of its specific clinical manifestation. We report here on a genome-wide search for loci controlling susceptibility to leprosy per se in a panel of 86 families including 205 siblings affected with leprosy from Southern Vietnam. Using model-free linkage analysis, we found significant evidence for a susceptibility gene on chromosome region 6q25 (maximum likelihood binomial (MLB) lod score 4.31; P = 5 x 10(-6)). We confirmed this by family-based association analysis in an independent panel of 208 Vietnamese leprosy simplex families. Of seven microsatellite markers underlying the linkage peak, alleles of two markers (D6S1035 and D6S305) showed strong evidence for association with leprosy (P = 6.7 x 10(-4) and P = 5.9 x 10(-5), respectively).     

A functional switch from lung cancer resistance to susceptibility at the Pas1 locus in Kras2LA2 mice

Pulmonary adenoma susceptibility 1 (Pas1) is the major mouse lung cancer susceptibility locus on chromosome 6 (ref. 1). Kras2 is a common target of somatic mutation in chemically induced mouse lung tumors and is a candidate Pas1 gene. M. spretus mice (SPRET/Ei) carry a Pas1 resistance haplotype for chemically induced lung tumors. We demonstrate that the SPRET/Ei Pas1 allele is switched from resistance to susceptibility by fixation of the parental origin of the mutant Kras2 allele. This switch correlates with low expression of endogenous Kras2 in SPRET/Ei. We propose that the Pas1 modifier effect is due to Kras2, and that a sensitive balance between the expression levels of wild-type and mutant alleles determines lung tumor susceptibility. These data demonstrate that cancer predisposition should also be considered in the context of somatic events and could have major implications for the design of human association studies to identify cancer susceptibility genes.     

Sequence variants in IL10, ARPC2 and multiple other loci contribute to ulcerative colitis susceptibility

Inflammatory bowel disease (IBD) typically manifests as either ulcerative colitis (UC) or Crohn's disease (CD). Systematic identification of susceptibility genes for IBD has thus far focused mainly on CD, and little is known about the genetic architecture of UC. Here we report a genome-wide association study with 440,794 SNPs genotyped in 1,167 individuals with UC and 777 healthy controls. Twenty of the most significantly associated SNPs were tested for replication in three independent European case-control panels comprising a total of 1,855 individuals with UC and 3,091 controls. Among the four consistently replicated markers, SNP rs3024505 immediately flanking the IL10 (interleukin 10) gene on chromosome 1q32.1 showed the most significant association in the combined verification samples (P = 1.35 x 10(-12); OR = 1.46 (1.31-1.62)). The other markers were located in ARPC2 and in the HLA-BTNL2 region. Association between rs3024505 and CD (1,848 cases, 1,804 controls) was weak (P = 0.013; OR = 1.17 (1.01-1.34)). IL10 is an immunosuppressive cytokine that has long been proposed to influence IBD pathophysiology. Our findings strongly suggest that defective IL10 function is central to the pathogenesis of the UC subtype of IBD.     

Screening a large reference sample to identify very low frequency sequence variants: comparisons between two genes

Most human sequence variation is in the form of single-nucleotide polymorphisms (SNPs). It has been proposed that coding-region SNPs (cSNPs) be used for direct association studies to determine the genetic basis of complex traits. The success of such studies depends on the frequency of disease-associated alleles, and their distribution in different ethnic populations. If disease-associated alleles are frequent in most populations, then direct genotyping of candidate variants could show robust associations in manageable study samples. This approach is less feasible if the genetic risk from a given candidate gene is due to many infrequent alleles. Previous studies of several genes demonstrated that most variants are relatively infrequent (<0.05). These surveys genotyped small samples (n<75) and thus had limited ability to identify rare alleles. Here we evaluate the prevalence and distribution of such rare alleles by genotyping an ethnically diverse reference sample that is more than six times larger than those used in previous studies (n=450). We screened for variants in the complete coding sequence and intron-exon junctions of two candidate genes for neuropsychiatric phenotypes: SLC6A4, encoding the serotonin transporter; and SLC18A2, encoding the vesicular monoamine transporter. Both genes have unique roles in neuronal transmission, and variants in either gene might be associated with neurobehavioral phenotypes.     

Association analyses identify six new psoriasis susceptibility loci in the Chinese population

We extended our previous genome-wide association study for psoriasis with a multistage replication study including 8,312 individuals with psoriasis (cases) and 12,919 controls from China as well as 3,293 cases and 4,188 controls from Germany and the United States and 254 nuclear families from the United States. We identified six new susceptibility loci associated with psoriasis in the Chinese study containing the candidate genes ERAP1, PTTG1, CSMD1, GJB2, SERPINB8 and ZNF816A (combined P < 5 × 10??) and replicated one locus, 5q33.1 (TNIP1-ANXA6), previously reported (combined P = 3.8 × 10?21) in the European studies. Two of these loci showed evidence for association in the German study at ZNF816A and GJB2 with P = 3.6 × 10?3 and P = 7.9 × 10?3, respectively. ERAP1 and ZNF816A were associated with type 1 (early onset) psoriasis in the Chinese Han population (test for heterogeneity P = 6.5 × 10?3 and P = 1.5 × 10?3, respectively). Comparisons with the results of previous GWAS of psoriasis highlight the heterogeneity of disease susceptibility between the Chinese and European populations. Our study identifies new genetic susceptibility factors and suggests new biological pathways in psoriasis.     

Genetic variation in DLG5 is associated with inflammatory bowel disease

Crohn disease and ulcerative colitis are two subphenotypes of inflammatory bowel disease (IBD), a complex disorder resulting from gene-environment interaction. We refined our previously defined linkage region for IBD on chromosome 10q23 and used positional cloning to identify genetic variants in DLG5 associated with IBD. DLG5 encodes a scaffolding protein involved in the maintenance of epithelial integrity. We identified two distinct haplotypes with a replicable distortion in transmission (P = 0.000023 and P = 0.004 for association with IBD, P = 0.00012 and P = 0.04 for association with Crohn disease). One of the risk-associated DLG5 haplotypes is distinguished from the common haplotype by a nonsynonymous single-nucleotide polymorphism 113G-->A, resulting in the amino acid substitution R30Q in the DUF622 domain of DLG5. This mutation probably impedes scaffolding of DLG5. We stratified the study sample according to the presence of risk-associated CARD15 variants to study potential gene-gene interaction. We found a significant difference in association of the 113A DLG5 variant with Crohn disease in affected individuals carrying the risk-associated CARD15 alleles versus those carrying non-risk-associated CARD15 alleles. This is suggestive of a complex pattern of gene-gene interaction between DLG5 and CARD15, reflecting the complex nature of polygenic diseases. Further functional studies will evaluate the biological significance of DLG5 variants.