Dense genotyping identifies and localizes multiple common and rare variant association signals in celiac disease

Identification of multiple risk variants for ankylosing spondylitis through high-density genotyping of immune-related loci

Ankylosing spondylitis is a common, highly heritable inflammatory arthritis affecting primarily the spine and pelvis. In addition to HLA-B*27 alleles, 12 loci have previously been identified that are associated with ankylosing spondylitis in populations of European ancestry, and 2 associated loci have been identified in Asians. In this study, we used the Illumina Immunochip microarray to perform a case-control association study involving 10,619 individuals with ankylosing spondylitis (cases) and 15,145 controls. We identified 13 new risk loci and 12 additional ankylosing spondylitis-associated haplotypes at 11 loci. Two ankylosing spondylitis-associated regions have now been identified encoding four aminopeptidases that are involved in peptide processing before major histocompatibility complex (MHC) class I presentation. Protective variants at two of these loci are associated both with reduced aminopeptidase function and with MHC class I cell surface expression.

High-density genotyping study identifies four new susceptibility loci for atopic dermatitis

Atopic dermatitis is a common inflammatory skin disease with a strong heritable component. Pathogenetic models consider keratinocyte differentiation defects and immune alterations as scaffolds, and recent data indicate a role for autoreactivity in at least a subgroup of patients. FLG (encoding filaggrin) has been identified as a major locus causing skin barrier deficiency. To better define risk variants and identify additional susceptibility loci, we densely genotyped 2,425 German individuals with atopic dermatitis (cases) and 5,449 controls using the Immunochip array followed by replication in 7,196 cases and 15,480 controls from Germany, Ireland, Japan and China. We identified four new susceptibility loci for atopic dermatitis and replicated previous associations. This brings the number of atopic dermatitis risk loci reported in individuals of European ancestry to 11. We estimate that these susceptibility loci together account for 14.4% of the heritability for atopic dermatitis.

Dense genotyping of immune-related disease regions identifies nine new risk loci for primary sclerosing cholangitis

Primary sclerosing cholangitis (PSC) is a severe liver disease of unknown etiology leading to fibrotic destruction of the bile ducts and ultimately to the need for liver transplantation. We compared 3,789 PSC cases of European ancestry to 25,079 population controls across 130,422 SNPs genotyped using the Immunochip. We identified 12 genome-wide significant associations outside the human leukocyte antigen (HLA) complex, 9 of which were new, increasing the number of known PSC risk loci to 16. Despite comorbidity with inflammatory bowel disease (IBD) in 72% of the cases, 6 of the 12 loci showed significantly stronger association with PSC than with IBD, suggesting overlapping yet distinct genetic architectures for these two diseases. We incorporated association statistics from 7 diseases clinically occurring with PSC in the analysis and found suggestive evidence for 33 additional pleiotropic PSC risk loci. Together with network analyses, these findings add to the genetic risk map of PSC and expand on the relationship between PSC and other immune-mediated diseases.

HLA class II gene associations in African American type 1 diabetes reveal a protective HLA-DRB1*03 haplotype

AIMS

Owing to strong linkage disequilibrium between markers, pinpointing disease associations within genetic regions is difficult in European ancestral populations, most notably the very strong association of the HLA-DRB1*03-DQA1*05:01-DQB1*02:01 haplotype with Type 1 diabetes risk, which is assumed to be because of a combination of HLA-DRB1 and HLA-DQB1. In contrast, populations of African ancestry have greater haplotype diversity, offering the possibility of narrowing down regions and strengthening support for a particular gene in a region being causal. We aimed to study the human leukocyte antigen (HLA) region in African American Type 1 diabetes.

METHODS

Two hundred and twenty-seven African American patients with Type 1 diabetes and 471 African American control subjects were tested for association at the HLA class II genes, HLA-DRB1, HLA-DQA1, HLA-DQB1 and 5147 single nucleotide polymorphisms across the major histocompatibility complex region using logistic regression models. Population admixture was accounted for with principal components analysis.

RESULTS

Single nucleotide polymorphism marker associations were explained by the HLA associations, with the major peak over the class II loci. The HLA association overall was extremely strong, as expected for Type 1 diabetes, even in African Americans in whom diabetes diagnosis is heterogeneous. In addition, there were unique features: the HLA-DRB1*03 haplotype was split into HLA-DRB1*03:01, which confers greatest susceptibility in these samples (odds ratio 3.17, 95% CI 1.72-5.83) and HLA-DRB1*03:02, an allele rarely observed in Europeans, which confers the greatest protection in these African American samples (odds ratio 0.22, 95% CI 0.09-0.55).

CONCLUSIONS

The unique diversity of the African HLA region we have uncovered supports a specific and major role for HLA-DRB1 in HLA-DRB1*03 haplotype-associated Type 1 diabetes risk.

Dense genotyping of immune-related disease regions identifies 14 new susceptibility loci for juvenile idiopathic arthritis

We used the Immunochip array to analyze 2,816 individuals with juvenile idiopathic arthritis (JIA), comprising the most common subtypes (oligoarticular and rheumatoid factor-negative polyarticular JIA), and 13,056 controls. We confirmed association of 3 known JIA risk loci (the human leukocyte antigen (HLA) region, PTPN22 and PTPN2) and identified 14 loci reaching genome-wide significance (P < 5 × 10(-8)) for the first time. Eleven additional new regions showed suggestive evidence of association with JIA (P < 1 × 10(-6)). Dense mapping of loci along with bioinformatics analysis refined the associations to one gene in each of eight regions, highlighting crucial pathways, including the interleukin (IL)-2 pathway, in JIA disease pathogenesis. The entire Immunochip content, the HLA region and the top 27 loci (P < 1 × 10(-6)) explain an estimated 18, 13 and 6% of the risk of JIA, respectively. In summary, this is the largest collection of JIA cases investigated so far and provides new insight into the genetic basis of this childhood autoimmune disease.

Single nucleotide polymorphisms in pathogen recognition receptor genes are associated with susceptibility to meningococcal meningitis in a pediatric cohort

Bacterial meningitis (BM) is a serious infection of the central nervous system, frequently occurring in childhood and often resulting in hearing loss, learning disabilities, and encephalopathy. Previous studies showed that genetic variation in innate immune response genes affects susceptibility, severity, and outcome of BM. The aim of this study is to describe whether single nucleotide polymorphisms (SNPs) in pathogen recognition gene products are associated with susceptibility to develop BM in single genes analysis as well as SNP combinations. Genotype frequencies of seven SNPs, in five immune response genes encoding for Toll-like receptors (TLRs), nucleotide oligomerization domain (NOD) proteins and caspase-1 (CASP1), in 391 children with meningococcal meningitis (MM) and 82 children with pneumococcal meningitis were compared with a large cohort of 1141 ethnically matched healthy controls. Carriage of TLR4 +896 GG mutant predisposed to susceptibility to develop MM (p = 1.2*10⁻⁵, OR  = 9.4, 95% CI  = 3.0–29.2). The NOD2 SNP8 mutant was significantly more frequent in MM patients compared to controls (p = 0.0004, OR  = 12.2, 95% CI  = 2.6–57.8). Combined carriage of TLR2 +2477 and TLR4 +896 mutants was strongly associated with MM (p = 4.2*10⁻⁵, OR  = 8.6, 95% CI  = 2.7–27.3). A carrier trait of TLR4 +896 and NOD2 SNP8 mutants was also strongly associated with susceptibility to develop MM (p = 4.2*10⁻⁵, OR  = 10.6, 95% CI  = 2.9–38.6). This study associates SNPs in TLR4 and NOD2 with susceptibility to develop MM.

Update on the diagnosis and management of refractory coeliac disease

A small subset of coeliac disease (CD) patients experiences persisting or recurring symptoms despite strict adherence to a gluten-free diet (GFD). When other causes of villous atrophy have been excluded, these patients are referred to as refractory celiac disease (RCD) patients. RCD can be divided in two types based on the absence (type I) or presence (type II) of an, usually clonal, intraepithelial lymphocyte population with aberrant phenotype. RCDI usually runs a benign course and may be difficult to be differentiated from uncomplicated, slow responding CD. In contrast, RCDII can be defined as low-grade intraepithelial lymphoma and frequently transforms into an aggressive enteropathy associated T-cell lymphoma with dismal prognosis. This paper describes the clinical characteristics of RCDI and RCDII, diagnostic approach, and the latest insights in treatment options.

GWAS in autoimmune thyroid disease: redefining our understanding of pathogenesis

The ability of the immune system to protect the body from attack by foreign antigens is essential for human survival. The immune system can, however, start to attack the body's own organs. An autoimmune response against components of the thyroid gland affects 2-5% of the general population. Considerable familial clustering is also observed in autoimmune thyroid disease (AITD). Teasing out the genetic contribution to AITD over the past 40 years has helped unravel how immune disruption leads to disease onset. Breakthroughs in genome-wide association studies (GWAS) in the past decade have facilitated screening of a greater proportion of the genome, leading to the identification of a before unimaginable number of AITD susceptibility loci. This Review will focus on the new susceptibility loci identified by GWAS, what insights these loci provide about the pathogenesis of AITD and how genetic susceptibility loci shared between different autoimmune diseases could help explain disease co-clustering within individuals and families. This Review also discusses where future efforts should be focused to translate this step forward in our understanding of the genetic contribution to AITD into a better understanding of disease presentation and progression, and improved therapeutic options.

Genome-wide association studies of allergic diseases

Allergic diseases are complex diseases caused by a combination of genetic and environmental factors. To determine the genetic components of these diseases and to discover the genes and cellular pathways underlying them, a large number of genetic studies have been conducted. Progress in genetics enables us to conduct genome-wide association studies (GWASs), which is a comprehensive and unbiased approach to identify susceptibility loci for multifactorial diseases. Recent GWASs have convincingly detected a large number of loci associated with allergic diseases. Candidate genes in the susceptibility loci suggest roles for epithelial barrier functions, innate-adaptive immunity, IL-1 family signaling, regulatory T cells and the vitamin D pathway in the pathogenesis of allergic diseases. Interestingly, the IL1RL1, HLA, IL13 and C11orf30 regions are overlapping susceptibility loci among atopic dermatitis and asthma or allergic rhinitis. Although a more complete collection of associated genes and pathways is needed, biologic insights revealed by GWASs improve our understanding of the pathophysiology of human allergic diseases and contribute to the development of better treatment and preventive strategies.

Trans-ethnic fine-mapping of lipid loci identifies population-specific signals and allelic heterogeneity that increases the trait variance explained

Genome-wide association studies (GWAS) have identified ∼100 loci associated with blood lipid levels, but much of the trait heritability remains unexplained, and at most loci the identities of the trait-influencing variants remain unknown. We conducted a trans-ethnic fine-mapping study at 18, 22, and 18 GWAS loci on the Metabochip for their association with triglycerides (TG), high-density lipoprotein cholesterol (HDL-C), and low-density lipoprotein cholesterol (LDL-C), respectively, in individuals of African American (n = 6,832), East Asian (n = 9,449), and European (n = 10,829) ancestry. We aimed to identify the variants with strongest association at each locus, identify additional and population-specific signals, refine association signals, and assess the relative significance of previously described functional variants. Among the 58 loci, 33 exhibited evidence of association at P<1×10⁻⁴ in at least one ancestry group. Sequential conditional analyses revealed that ten, nine, and four loci in African Americans, Europeans, and East Asians, respectively, exhibited two or more signals. At these loci, accounting for all signals led to a 1.3- to 1.8-fold increase in the explained phenotypic variance compared to the strongest signals. Distinct signals across ancestry groups were identified at PCSK9 and APOA5. Trans-ethnic analyses narrowed the signals to smaller sets of variants at GCKR, PPP1R3B, ABO, LCAT, and ABCA1. Of 27 variants reported previously to have functional effects, 74% exhibited the strongest association at the respective signal. In conclusion, trans-ethnic high-density genotyping and analysis confirm the presence of allelic heterogeneity, allow the identification of population-specific variants, and limit the number of candidate SNPs for functional studies.

Lipid traits are heritable, but many of the DNA variants that influence lipid levels remain unknown. In a genomic region, more than one variant may affect gene expression or function, and the frequencies of these variants can differ across populations. Genotyping densely spaced variants in individuals with different ancestries may increase the chance of identifying variants that affect gene expression or function. We analyzed high-density genotyped variants for association with TG, HDL-C, and LDL-C in African Americans, East Asians, and Europeans. At several genomic regions, we provide evidence that two or more variants can influence lipid traits; across loci, these additional signals increase the proportion of trait variation that can be explained by genes. At some association signals shared across populations, combining data from individuals of different ancestries narrowed the set of likely functional variants. At PCSK9 and APOA5, the data suggest that different variants influence trait levels in different populations. Variants previously reported to alter gene expression or function frequently exhibited the strongest association at those signals. The multiple signals and population-specific characteristics of the loci described here may be shared by genetic loci for other complex traits.

A population-specific uncommon variant in GRIN3A associated with schizophrenia

BACKGROUND

Genome-wide association studies have successfully identified several common variants showing robust association with schizophrenia. However, individually, these variants only produce a weak effect. To identify genetic variants with larger effect sizes, increasing attention is now being paid to uncommon and rare variants.

METHODS

From the 1000 Genomes Project data, we selected 47 candidate single nucleotide variants (SNVs), which were: 1) uncommon (minor allele frequency < 5%); 2) Asian-specific; 3) missense, nonsense, or splice site variants predicted to be damaging; and 4) located in candidate genes for schizophrenia and bipolar disorder. We examined their association with schizophrenia, using a Japanese case-control cohort (2012 cases and 2781 control subjects). Additional meta-analysis was performed using genotyping data from independent Han-Chinese case-control (333 cases and 369 control subjects) and family samples (9 trios and 284 quads).

RESULTS

We identified disease association of a missense variant in GRIN3A (p.R480G, rs149729514, p = .00042, odds ratio [OR] = 1.58), encoding a subunit of the N-methyl-D-aspartate type glutamate receptor, with study-wide significance (threshold p = .0012). This association was supported by meta-analysis (combined p = 3.3 × 10(-5), OR = 1.61). Nominally significant association was observed in missense variants from FAAH, DNMT1, MYO18B, and CFB, with ORs of risk alleles ranging from 1.41 to 2.35.

CONCLUSIONS

The identified SNVs, particularly the GRIN3A R480G variant, are good candidates for further replication studies and functional evaluation. The results of this study indicate that association analyses focusing on uncommon and rare SNVs are a promising way to discover risk variants with larger effects.

Triggers and drivers of autoimmunity: lessons from coeliac disease

Coeliac disease, an inflammatory disease of the small intestine, shares key features with autoimmune disorders, such as susceptibility genes, presence of autoantibodies and T cell-mediated destruction of specific cells. Strikingly, however, continuous exposure to the exogenous dietary antigen gluten and gluten-specific adaptive immunity are required to maintain immunopathology. These observations challenge the notion that autoimmunity requires adaptive immune activation towards self antigens. Using coeliac disease as an example, we propose that other exogenous factors might be identified as drivers of autoimmune processes, in particular when evidence for T cells with specificity for self antigens driving the disease is lacking.

Identification of multiple independent susceptibility loci in the HLA region in Behçet's disease

Behçet's disease is an inflammatory disease characterized by recurrent oral and genital ulcers and significant organ involvement. Localizing the genetic association between HLA-B*51 and Behçet's disease and exploring additional susceptibility loci in the human leukocyte antigen (HLA) region are complicated by the strong linkage disequilibrium in this region. We genotyped 8,572 variants in the extended HLA locus and carried out imputation and meta-analysis of 24,834 variants in 2 independent Behçet's disease cohorts from 2 ancestry groups. Genotyped SNPs were used to infer classical HLA alleles in the HLA-A, HLA-B, HLA-C, HLA-DQA1, HLA-DQB1 and HLA-DRB1 loci. Our data suggest that the robust HLA-B*51 association in Behçet's disease is explained by a variant located between the HLA-B and MICA genes (rs116799036: odds ratio (OR) = 3.88, P = 9.42 × 10(-50)). Three additional independent genetic associations within PSORS1C1 (rs12525170: OR = 3.01, P = 3.01 × 10(-26)), upstream of HLA-F-AS1 (rs114854070: OR = 1.95, P = 7.84 × 10(-14)) and with HLA-Cw*1602 (OR = 5.38, P = 6.07 × 10(-18)) were also identified and replicated.

Allelic variants in TLR10 gene may influence bilateral affectation and clinical course of Meniere's disease

Toll-like receptors trigger the innate immune response by activating various cell types such us macrophages and lymphocytes. We genotyped SNV of TLR3, TRL7, TLR8 and TLR10 in 863 Spanish and 150 Italian patients with Meniere's disease (MD) and 1,013 controls by using Taqman assays. Real-Time qPCR was used to measure the expression level of TLR10 in peripheral blood leukocytes. The overall dataset showed that the C allele and the CC genotype of rs11096955 in TLR10 gene were more commonly observed in controls than patients (corrected p = 1 × 10(-3), OR = 0.68 [95 % confidence interval, 0.54-0.84] for CC genotype; corrected p = 1.5 × 10(-5), OR = 0.75 [0.66-0.85] for allele C). Moreover, the CC genotype was more frequent in patients with uni- (19 %) than bilateral sensorineural hearing loss (SNHL) (13 %). Logistic regression demonstrated that the time since the onset of MD, Tumarkin crises, hearing stage and rs11096955 were independent factors influencing the risk of bilateral SNHL. In addition, rs11096955 influenced hearing loss progression in patients with bilateral MD. No change in expression of TLR10 was observed according to CC, CA or AA genotypes. Our data suggest that allelic variants of TLR10 gene may influence the susceptibility and time-course of hearing loss of MD in the European population.

Role of NOD1 polymorphism in susceptibility and clinical progression of rheumatoid arthritis

OBJECTIVE

One of the disease hallmarks of RA is progressive cartilage and bone destruction in the joints. The exact mechanism underlying this disease process is largely unknown. Nod1, an intracellular pattern recognition receptor expressed by the innate immune system, has been previously shown to display anti-inflammatory effects in experimental arthritis. Furthermore, an insertion/deletion polymorphism in NOD1 has been demonstrated to modulate cytokine responses of immune cells. In this study, the effect of the insertion/deletion polymorphism in NOD1 on RA susceptibility and severity was assessed.

METHODS

Ex vivo stimulation of primary immune cells and osteoclasts with microbial triggers was performed to measure cytokine responses and osteoclast-specific gene expression in relation to the NOD1 genotype. In total, 1047 RA patients from two centres were genotyped for the NOD1 polymorphism and compared with 431 healthy controls. Clinical scores of joint inflammation and destruction were correlated with the NOD1 genotype.

RESULTS

Functional analysis revealed increased production of pro-inflammatory cytokines in cells from individuals bearing the NOD1 +32656 insertion allele. Furthermore, osteoclast bone resorption activity was elevated, as reflected by increased expression of the lysosomal protease cathepsin K. However, the insertion allele of the NOD1 +32656 polymorphism was not associated with either susceptibility to, or clinical parameters of, inflammation or bone destruction in RA patients.

CONCLUSION

These findings demonstrate that the NOD1 polymorphism modulates pro-inflammatory cytokine responses induced through Toll-like receptor or Nod-like receptor ligands. Nevertheless, these effects of genetic variation in NOD1 appear to be redundant in RA susceptibility and severity.

Review of Ets1 structure, function, and roles in immunity

The Ets1 transcription factor is a member of the Ets gene family and is highly conserved throughout evolution. Ets1 is known to regulate a number of important biological processes in normal cells and in tumors. In particular, Ets1 has been associated with regulation of immune cell function and with an aggressive behavior in tumors that express it at high levels. Here we review and summarize the general features of Ets1 and describe its roles in immunity and autoimmunity, with a focus on its roles in B lymphocytes. We also review evidence that suggests that Ets1 may play a role in malignant transformation of hematopoietic malignancies including B cell malignancies.

The emerging spectrum of allelic variation in schizophrenia: current evidence and strategies for the identification and functional characterization of common and rare variants

After decades of halting progress, recent large genome-wide association studies (GWAS) are finally shining light on the genetic architecture of schizophrenia. The picture emerging is one of sobering complexity, involving large numbers of risk alleles across the entire allelic spectrum. The aims of this article are to summarize the key genetic findings to date and to compare and contrast methods for identifying additional risk alleles, including GWAS, targeted genotyping and sequencing. A further aim is to consider the challenges and opportunities involved in determining the functional basis of genetic associations, for instance using functional genomics, cellular models, animal models and imaging genetics. We conclude that diverse approaches will be required to identify and functionally characterize the full spectrum of risk variants for schizophrenia. These efforts should adhere to the stringent standards of statistical association developed for GWAS and are likely to entail very large sample sizes. Nonetheless, now more than any previous time, there are reasons for optimism and the ultimate goal of personalized interventions and therapeutics, although still distant, no longer seems unattainable.

Detecting rare variants for psychiatric disorders using next generation sequencing: a methods primer

Recent advances in massively parallel sequencing (MPS) have had an extensive impact on research in medical genomics. In particular, the analysis of rare variants using MPS promises to lead to a better understanding of complex disorders. Nevertheless, for meaningful studies that address the genetic basis for neuropsychiatric disorders, at least hundreds of patient samples have to be analyzed. This undertaking is still not feasible for single research groups on a whole-genome scale and in individual samples. Thus, researchers increasingly employ strategies for reducing the amount of sequencing efforts, such as target enrichment and non-barcoded sample pooling. This review provides an overview of current technologies, discusses options for reduced experimental designs, and illustrates the successful application of the presented methodologies in a recent study of panic disorder patients. Thereby, it aims to introduce the emerging field of MPS into neuropsychiatric research and might serve as a guide for further studies.

Rare variants in TMEM132D in a case-control sample for panic disorder

Genome-wide association studies have identified common variants associated with common diseases. Most variants, however, explain only a small proportion of the estimated heritability, suggesting that rare variants might contribute to a larger extent to common diseases than assumed to date. Here, we use next-generation sequencing to test whether such variants contribute to the risk for anxiety disorders by re-sequencing 40 kb including all exons of the TMEM132D locus which we have previously shown to be associated with panic disorder and anxiety severity measures. DNA from 300 patients suffering from anxiety disorders, mostly panic disorder (84.7%), and 300 healthy controls was screened for the presence of genetic variants using next-generation re-sequencing in a pooled approach. Results were verified by individual re-genotyping. We identified 371 variants of which 247 had not been reported before, including 15 novel non-synonymous variants. The majority, 76% of these variants had a minor allele frequency less than 5%. While we did not identify additional common variants in TMEM132D associated with panic disorders, we observed an overrepresentation of presumably functional coding variants in healthy controls as compared to cases as well as a higher rate of private coding variants in cases, with one non-synonymous coding variant present in four patients but not in any of the matched controls nor in over 5,500 individuals of different ethnic origins from publicly available re-sequencing datasets. Our data suggest that not only common but also putatively functional and/or rare variants within TMEM132D might contribute to the risk to develop anxiety disorders.

Identification of 15 new psoriasis susceptibility loci highlights the role of innate immunity

To gain further insight into the genetic architecture of psoriasis, we conducted a meta-analysis of 3 genome-wide association studies (GWAS) and 2 independent data sets genotyped on the Immunochip, including 10,588 cases and 22,806 controls. We identified 15 new susceptibility loci, increasing to 36 the number associated with psoriasis in European individuals. We also identified, using conditional analyses, five independent signals within previously known loci. The newly identified loci shared with other autoimmune diseases include candidate genes with roles in regulating T-cell function (such as RUNX3, TAGAP and STAT3). Notably, they included candidate genes whose products are involved in innate host defense, including interferon-mediated antiviral responses (DDX58), macrophage activation (ZC3H12C) and nuclear factor (NF)-κB signaling (CARD14 and CARM1). These results portend a better understanding of shared and distinctive genetic determinants of immune-mediated inflammatory disorders and emphasize the importance of the skin in innate and acquired host defense.

High-density genetic mapping identifies new susceptibility loci for rheumatoid arthritis

Using the Immunochip custom SNP array, which was designed for dense genotyping of 186 loci identified through genome-wide association studies (GWAS), we analyzed 11,475 individuals with rheumatoid arthritis (cases) of European ancestry and 15,870 controls for 129,464 markers. We combined these data in a meta-analysis with GWAS data from additional independent cases (n = 2,363) and controls (n = 17,872). We identified 14 new susceptibility loci, 9 of which were associated with rheumatoid arthritis overall and five of which were specifically associated with disease that was positive for anticitrullinated peptide antibodies, bringing the number of confirmed rheumatoid arthritis risk loci in individuals of European ancestry to 46. We refined the peak of association to a single gene for 19 loci, identified secondary independent effects at 6 loci and identified association to low-frequency variants at 4 loci. Bioinformatic analyses generated strong hypotheses for the causal SNP at seven loci. This study illustrates the advantages of dense SNP mapping analysis to inform subsequent functional investigations.

Seven newly identified loci for autoimmune thyroid disease

Autoimmune thyroid disease (AITD), including Graves' disease (GD) and Hashimoto's thyroiditis (HT), is one of the most common of the immune-mediated diseases. To further investigate the genetic determinants of AITD, we conducted an association study using a custom-made single-nucleotide polymorphism (SNP) array, the ImmunoChip. The SNP array contains all known and genotype-able SNPs across 186 distinct susceptibility loci associated with one or more immune-mediated diseases. After stringent quality control, we analysed 103 875 common SNPs (minor allele frequency >0.05) in 2285 GD and 462 HT patients and 9364 controls. We found evidence for seven new AITD risk loci (P < 1.12 × 10(-6); a permutation test derived significance threshold), five at locations previously associated and two at locations awaiting confirmation, with other immune-mediated diseases.

Immunochip analyses identify a novel risk locus for primary biliary cirrhosis at 13q14, multiple independent associations at four established risk loci and epistasis between 1p31 and 7q32 risk variants

To further characterize the genetic basis of primary biliary cirrhosis (PBC), we genotyped 2426 PBC patients and 5731 unaffected controls from three independent cohorts using a single nucleotide polymorphism (SNP) array (Immunochip) enriched for autoimmune disease risk loci. Meta-analysis of the genotype data sets identified a novel disease-associated locus near the TNFSF11 gene at 13q14, provided evidence for association at six additional immune-related loci not previously implicated in PBC and confirmed associations at 19 of 22 established risk loci. Results of conditional analyses also provided evidence for multiple independent association signals at four risk loci, with haplotype analyses suggesting independent SNP effects at the 2q32 and 16p13 loci, but complex haplotype driven effects at the 3q25 and 6p21 loci. By imputing classical HLA alleles from this data set, four class II alleles independently contributing to the association signal from this region were identified. Imputation of genotypes at the non-HLA loci also provided additional associations, but none with stronger effects than the genotyped variants. An epistatic interaction between the IL12RB2 risk locus at 1p31and the IRF5 risk locus at 7q32 was also identified and suggests a complementary effect of these loci in predisposing to disease. These data expand the repertoire of genes with potential roles in PBC pathogenesis that need to be explored by follow-up biological studies.

Performance and robustness of penalized and unpenalized methods for genetic prediction of complex human disease

A central goal of medical genetics is to accurately predict complex disease from genotypes. Here, we present a comprehensive analysis of simulated and real data using lasso and elastic-net penalized support-vector machine models, a mixed-effects linear model, a polygenic score, and unpenalized logistic regression. In simulation, the sparse penalized models achieved lower false-positive rates and higher precision than the other methods for detecting causal SNPs. The common practice of prefiltering SNP lists for subsequent penalized modeling was examined and shown to substantially reduce the ability to recover the causal SNPs. Using genome-wide SNP profiles across eight complex diseases within cross-validation, lasso and elastic-net models achieved substantially better predictive ability in celiac disease, type 1 diabetes, and Crohn's disease, and had equivalent predictive ability in the rest, with the results in celiac disease strongly replicating between independent datasets. We investigated the effect of linkage disequilibrium on the predictive models, showing that the penalized methods leverage this information to their advantage, compared with methods that assume SNP independence. Our findings show that sparse penalized approaches are robust across different disease architectures, producing as good as or better phenotype predictions and variance explained. This has fundamental ramifications for the selection and future development of methods to genetically predict human disease.

Population-based case-control association studies

This unit provides an overview of the design and analysis of population-based case-control studies of genetic risk factors for complex disease. Considerations specific to genetic studies are emphasized. The unit reviews basic study designs differentiating case-control studies from others, presents different genetic association strategies (candidate gene, genome-wide association, and high-throughput sequencing), introduces basic methods of statistical analysis for case-control data and approaches to combining case-control studies, and discusses measures of association and impact. Admixed populations, controlling for confounding (including population stratification), consideration of multiple loci and environmental risk factors, and complementary analyses of haplotypes, genes, and pathways are briefly discussed. Readers are referred to basic texts on epidemiology for more details on general conduct of case-control studies.

HLA-DQA1 and HLA-DQB1 in Celiac disease predisposition: practical implications of the HLA molecular typing

Celiac disease (CD) is a multifactorial disorder with an estimated prevalence in Europe and USA of 1:100 and a female:male ratio of approximately 2:1. The disorder has a multifactorial etiology in which the triggering environmental factor, the gluten, and the main genetic factors, Human Leukocyte Antigen (HLA)-DQA1 and HLA-DQB1 loci, are well known. About 90-95% of CD patients carry DQ2.5 heterodimers, encoded by DQA1*05 and DQB1*02 alleles both in cis or in trans configuration, and DQ8 molecules, encoded by DQB1*03:02 generally in combination with DQA1*03 variant. Less frequently, CD occurs in individuals positive for the DQ2.x heterodimers (DQA1≠*05 and DQB1*02) and very rarely in patients negative for these DQ predisposing markers. HLA molecular typing for Celiac disease is, therefore, a genetic test with a negative predictive value. Nevertheless, it is an important tool able to discriminate individuals genetically susceptible to CD, especially in at-risk groups such as first-degree relatives (parents, siblings and offspring) of patients and in presence of autoimmune conditions (type 1 diabetes, thyroiditis, multiple sclerosis) or specific genetic disorders (Down, Turner or Williams syndromes).

Age-dependent variation of genotypes in MHC II transactivator gene (CIITA) in controls and association to type 1 diabetes

The major histocompatibility complex class II transactivator (CIITA) gene (16p13) has been reported to associate with susceptibility to multiple sclerosis, rheumatoid arthritis and myocardial infarction, recently also to celiac disease at genome-wide level. However, attempts to replicate association have been inconclusive. Previously, we have observed linkage to the CIITA region in Scandinavian type 1 diabetes (T1D) families. Here we analyze five Swedish T1D cohorts and a combined control material from previous studies of CIITA. We investigate how the genotype distribution within the CIITA gene varies depending on age, and the association to T1D. Unexpectedly, we find a significant difference in the genotype distribution for markers in CIITA (rs11074932, P=4 × 10(-5) and rs3087456, P=0.05) with respect to age, in the collected control material. This observation is replicated in an independent cohort material of about 2000 individuals (P=0.006, P=0.007). We also detect association to T1D for both markers, rs11074932 (P=0.004) and rs3087456 (P=0.001), after adjusting for age at sampling. The association remains independent of the adjacent T1D risk gene CLEC16A. Our results indicate an age-dependent variation in CIITA allele frequencies, a finding of relevance for the contrasting outcomes of previously published association studies.

CSK regulatory polymorphism is associated with systemic lupus erythematosus and influences B-cell signaling and activation

C-src tyrosine kinase, Csk, physically interacts with the intracellular phosphatase Lyp (PTPN22) and can modify the activation state of downstream Src kinases, such as Lyn, in lymphocytes. We identified an association of Csk with systemic lupus erythematosus (SLE) and refined its location to an intronic polymorphism rs34933034 (OR 1.32, p = 1.04 × 10⁻⁹). The risk allele is associated with increased CSK expression and augments inhibitory phosphorylation of Lyn. In carriers of the risk allele, B cell receptor (BCR)-mediated activation of mature B cells, as well as plasma IgM, are increased. Moreover, the fraction of transitional B cells is doubled in the cord blood of carriers of the risk allele compared to non-risk haplotypes due to an expansion of the late transitional cells, a stage targeted by selection mechanisms. This suggests that the Lyp-Csk complex increases susceptibility to lupus at multiple maturation and activation points of B cells.

Human genetics offers an emerging picture of common pathways and mechanisms in autoimmunity

In genetic studies of autoimmune and inflammatory diseases, one clear finding that has emerged from genome-wide association studies is that a substantial fraction of variation modifying risk in one disease also contributes mediate risk to multiple, additional autoimmune and inflammatory diseases. The unexpected magnitude of this overlap presents the unique opportunity to dissect the pathogenic mechanisms underlying multiple disease states in the expectation that this may lead to both more sensitive diagnostics and novel therapies. Here, we review the current evidence for this shared genetic architecture and, based on these data, outline models for shared pathways, the underlying hypotheses for them, how these models can be tested and validated.

Dense fine-mapping study identifies new susceptibility loci for primary biliary cirrhosis

We genotyped 2,861 cases of primary biliary cirrhosis (PBC) from the UK PBC Consortium and 8,514 UK population controls across 196,524 variants within 186 known autoimmune risk loci. We identified 3 loci newly associated with PBC (at P<5×10(-8)), increasing the number of known susceptibility loci to 25. The most associated variant at 19p12 is a low-frequency nonsynonymous SNP in TYK2, further implicating JAK-STAT and cytokine signaling in disease pathogenesis. An additional five loci contained nonsynonymous variants in high linkage disequilibrium (LD; r2>0.8) with the most associated variant at the locus. We found multiple independent common, low-frequency and rare variant association signals at five loci. Of the 26 independent non-human leukocyte antigen (HLA) signals tagged on the Immunochip, 15 have SNPs in B-lymphoblastoid open chromatin regions in high LD (r2>0.8) with the most associated variant. This study shows how data from dense fine-mapping arrays coupled with functional genomic data can be used to identify candidate causal variants for functional follow-up.

Common variants at the MHC locus and at chromosome 16q24.1 predispose to Barrett's esophagus

Barrett's esophagus is an increasingly common disease that is strongly associated with reflux of stomach acid and usually a hiatus hernia, and it strongly predisposes to esophageal adenocarcinoma (EAC), a tumor with a very poor prognosis. We report the first genome-wide association study on Barrett's esophagus, comprising 1,852 UK cases and 5,172 UK controls in the discovery stage and 5,986 cases and 12,825 controls in the replication stage. Variants at two loci were associated with disease risk: chromosome 6p21, rs9257809 (Pcombined=4.09×10(-9); odds ratio (OR)=1.21, 95% confidence interval (CI)=1.13-1.28), within the major histocompatibility complex locus, and chromosome 16q24, rs9936833 (Pcombined=2.74×10(-10); OR=1.14, 95% CI=1.10-1.19), for which the closest protein-coding gene is FOXF1, which is implicated in esophageal development and structure. We found evidence that many common variants of small effect contribute to genetic susceptibility to Barrett's esophagus and that SNP alleles predisposing to obesity also increase risk for Barrett's esophagus.

Pathophysiology of celiac disease

Celiac disease results from the interplay of genetic, environmental, and immunologic factors. An understanding of the pathophysiology of celiac disease, in which the trigger (wheat, rye, and barley) is known, will undoubtedly reveal basic mechanisms that underlie other autoimmune diseases (eg, type 1 diabetes) that share many common pathogenic perturbations. This review describes seminal findings in each of the 3 domains of the pathogenesis of celiac disease, namely genetics, environmental triggers, and immune dysregulation, with a focus on newer areas of investigation such as non-HLA genetic variants, the intestinal microbiome, and the role of the innate immune system.

GWAS implicates a role for quantitative immune traits and threshold effects in risk for human autoimmune disorders

Genome wide association studies in human autoimmune disorders have provided a long list of alleles with rather modest degrees of risk. A large fraction of these associations are probably owing to either quantitative differences in gene expression or amino acid changes that regulate quantitative aspects of the immune response. While functional studies are still lacking for most of these associations, we present examples of autoimmune disease risk alleles that influence quantitative changes in lymphocyte activation, cytokine signaling and dendritic cell function. The analysis of immune quantitative traits associated with autoimmune loci is clearly going to be an important component of understanding the pathogenesis of autoimmunity. This will require both new and more efficient ways of characterizing the normal immune system, as well as large population resources in which genotype-phenotype correlations can be convincingly demonstrated. Future development of new therapies will depend on understanding the mechanistic underpinnings of immune regulation by these new risk loci.

Limited evidence for parent-of-origin effects in inflammatory bowel disease associated loci

Background

Genome-wide association studies of two main forms of inflammatory bowel diseases (IBD), Crohn’s disease (CD) and ulcerative colitis (UC), have identified 99 susceptibility loci, but these explain only ∼23% of the genetic risk. Part of the ‘hidden heritability’ could be in transmissible genetic effects in which mRNA expression in the offspring depends on the parental origin of the allele (genomic imprinting), since children whose mothers have CD are more often affected than children with affected fathers. We analyzed parent-of-origin (POO) effects in Dutch and Indian cohorts of IBD patients.

Methods

We selected 28 genetic loci associated with both CD and UC, and tested them for POO effects in 181 Dutch IBD case-parent trios. Three susceptibility variants in NOD2 were tested in 111 CD trios and a significant finding was re-evaluated in 598 German trios. The UC-associated gene, BTNL2, reportedly imprinted, was tested in 70 Dutch UC trios. Finally, we used 62 independent Indian UC trios to test POO effects of five established Indian UC risk loci.

Results

We identified POO effects for NOD2 (L1007fs; OR = 21.0, P-value = 0.013) for CD; these results could not be replicated in an independent cohort (OR = 0.97, P-value = 0.95). A POO effect in IBD was observed for IL12B (OR = 3.2, P-value = 0.019) and PRDM1 (OR = 5.6, P-value = 0.04). In the Indian trios the IL10 locus showed a POO effect (OR = 0.2, P-value = 0.03).

Conclusions

Little is known about the effect of genomic imprinting in complex diseases such as IBD. We present limited evidence for POO effects for the tested IBD loci. POO effects explain part of the hidden heritability for complex genetic diseases but need to be investigated further.

What have we learned from six years of GWAS in autoimmune diseases, and what is next?

Genome-wide association studies (GWAS) have discovered hundreds of common genetic variants that predispose humans to autoimmune diseases, opening up unprecedented potential for elucidating the pathways and processes of disease. To understand the role of these variants in susceptibility, we need to derive mechanistic insight by integration of genetic results with other biological data types and also with careful functional studies. In many cases, such studies have highlighted coherent biological processes at a high level and elucidated specific mechanisms that contribute to autoimmunity and inflammation. The understanding of the genetic component of autoimmune etiology will become more complete as fine-mapping and sequencing data become readily available. A comprehensive catalog of human immune phenotypes could provide a functional basis for assessing genetic influence on immune function and variation in response to therapeutic interventions, as well as for rationally designing new targeted therapeutics.

Functional variants in NFKBIE and RTKN2 involved in activation of the NF-κB pathway are associated with rheumatoid arthritis in Japanese

Rheumatoid arthritis is an autoimmune disease with a complex etiology, leading to inflammation of synovial tissue and joint destruction. Through a genome-wide association study (GWAS) and two replication studies in the Japanese population (7,907 cases and 35,362 controls), we identified two gene loci associated with rheumatoid arthritis susceptibility (NFKBIE at 6p21.1, rs2233434, odds ratio (OR) = 1.20, P = 1.3×10⁻¹⁵; RTKN2 at 10q21.2, rs3125734, OR = 1.20, P = 4.6×10⁻⁹). In addition to two functional non-synonymous SNPs in NFKBIE, we identified candidate causal SNPs with regulatory potential in NFKBIE and RTKN2 gene regions by integrating in silico analysis using public genome databases and subsequent in vitro analysis. Both of these genes are known to regulate the NF-κB pathway, and the risk alleles of the genes were implicated in the enhancement of NF-κB activity in our analyses. These results suggest that the NF-κB pathway plays a role in pathogenesis and would be a rational target for treatment of rheumatoid arthritis.

Rheumatoid arthritis (RA) is a chronic autoimmune disease affecting approximately 1% of the general adult population. More than 30 susceptibility loci for RA have been identified through genome-wide association studies (GWAS), but the disease-causal variants at most loci remain unknown. Here, we performed replication studies of the candidate loci of our previous GWAS using Japanese cohorts and identified variants in NFKBIE and RTKN2 gene loci that were associated with RA. To search for causal variants in both gene regions, we first examined non-synonymous (ns)SNPs that alter amino-acid sequences. As NFKBIE and RTKN2 are known to be involved in the NF-κB pathway, we evaluated the effects of nsSNPs on NF-κB activity. Next, we screened in silico variants that may regulate gene transcription using publicly available epigenetic databases and subsequently evaluated their regulatory potential using in vitro assays. As a result, we identified multiple candidate causal variants in NFKBIE (2 nsSNPs and 1 regulatory SNP) and RTKN2 (2 regulatory SNPs), indicating that our integrated in silico and in vitro approach is useful for the identification of causal variants in the post–GWAS era.

Celiac disease: an immunological jigsaw

Celiac disease (CD) is a chronic enteropathy induced by dietary gluten in genetically predisposed people. The keystone of CD pathogenesis is an adaptive immune response orchestrated by the interplay between gluten and MHC class II HLA-DQ2 and DQ8 molecules. Yet, other factors that impair immunoregulatory mechanisms and/or activate the large population of intestinal intraepithelial lymphocytes (IEL) are indispensable for driving tissue damage. Herein, we summarize our current understanding of the mechanisms and consequences of the undesirable immune response initiated by gluten peptides. We show that CD is a model disease to decipher the role of MHC class II molecules in human immunopathology, to analyze the mechanisms that link tolerance to food proteins and autoimmunity, and to investigate how chronic activation of IEL can lead to T cell lymphomagenesis.

Targeted modification of wheat grain protein to reduce the content of celiac causing epitopes

The prolamin peptides in wheat gluten and in the homologous storage proteins of barley and rye cause painful chronic erasure of microvilli of the small intestine epithelium in celiac patients. If untreated, it can lead to chronic diarrhea, abdominal distension, osteoporosis, weight-loss due to malabsorption of nutrients, and anemia. In addition to congenital cases, life-long exposure to gluten proteins in bread and pasta can also induce development of celiac sprue in adults. To date, the only effective treatment is life-long strict abstinence from the staple food grains. Complete exclusion of dietary gluten is, however, difficult due to use of wheat in many foods, incomplete labeling and social constraints. Thus, finding alternative therapies for this most common foodborne disease remained an active area of research, which has led to many suggestions in last few years. The pros and cons associated with these therapies were reviewed in the present communication. As different celiac patients are immunogenic to different members of the undigestible proline/glutamine rich peptides of ~149 gliadins and low molecular weight glutenin subunits as well as the six high molecular weight glutenin subunits, an exhaustive digestion of the immunogenic peptides in the stomach, duodenum, jejunum, and ileum of celiacs is required. In view of the above, we evaluated the capacity of cereal grains to synthesize and store the enzymes prolyl endopeptidase from Flavobacterium meningosepticum and the barley cysteine endoprotease B2, which in combination are capable of detoxifying immunogenic gluten peptides in a novel treatment of celiac disease.

Refractory celiac disease: from bench to bedside

Refractory celiac disease is defined by the persistence of symptoms of malnutrition and intestinal villous atrophy for more than 6-12 months despite strict gluten-free diet in celiac patients. Diagnosis of this rare condition is made after excluding other causes of chronic small intestinal inflammation and villous atrophy and inadvertent intake of gluten. Over the past 15 years, multidisciplinary approaches have been developed to assess the mechanism of resistance to the diet, and two distinct entities have been delineated. Type II refractory celiac disease (RCD) can be defined as a low-grade intraepithelial lymphoma. RCD II is characterised by a massive accumulation of abnormal IEL that display an aberrant hybrid NK/T cell phenotype, contain clonal T cell rearrangement(s) and can mediate a cytolytic attack of the gut epithelium. This condition has a severe prognosis, largely due to the frequent transformation of RCDII IEL into overt aggressive enteropathy-type-associated T cell lymphoma. In contrast, in type I RCD, intestinal lymphocytes have a normal phenotype, and this generally milder condition remains often difficult to differentiate from uncomplicated CD except for the resistance to gluten-free diet (GFD). Several mechanisms may underlie resistance to gluten. Herein, we review the distinctive characteristics of RCD I and RCD II, the mechanisms underlying the onset of resistance to GFD, the risk of developing high grade lymphoma and possible clues to improve their treatment.

Genetic variations in interleukin-12 related genes in immune-mediated diseases

The interleukin-12 (IL-12) family comprises a group of heterodimeric cytokines and their respective receptors that play key roles in immune responses. A growing number of autoimmune diseases has been found to be associated with genetic variation in these genes. Based on their respective associations with the IL-12 genes, autoimmune diseases appear to cluster in two groups that either show strong associations with the Th1/Th17 pathway (as indicated by genetic association with IL12B and IL23R) or the Th1/IL-35 pathway as the consequence of their association with polymorphisms in the IL12A gene region. The genetic associations are described in relation to what is known of the functionality of these genes in the various diseases. Comparing association data for gene families in different diseases may lead to better insight in the function of the genes in the onset and course of the disease.

The chbG gene of the chitobiose (chb) operon of Escherichia coli encodes a chitooligosaccharide deacetylase

The chb operon of Escherichia coli is involved in the utilization of the β-glucosides chitobiose and cellobiose. The function of chbG (ydjC), the sixth open reading frame of the operon that codes for an evolutionarily conserved protein is unknown. We show that chbG encodes a monodeacetylase that is essential for growth on the acetylated chitooligosaccharides chitobiose and chitotriose but is dispensable for growth on cellobiose and chitosan dimer, the deacetylated form of chitobiose. The predicted active site of the enzyme was validated by demonstrating loss of function upon substitution of its putative metal-binding residues that are conserved across the YdjC family of proteins. We show that activation of the chb promoter by the regulatory protein ChbR is dependent on ChbG, suggesting that deacetylation of chitobiose-6-P and chitotriose-6-P is necessary for their recognition by ChbR as inducers. Strains carrying mutations in chbR conferring the ability to grow on both cellobiose and chitobiose are independent of chbG function for induction, suggesting that gain of function mutations in ChbR allow it to recognize the acetylated form of the oligosaccharides. ChbR-independent expression of the permease and phospho-β-glucosidase from a heterologous promoter did not support growth on both chitobiose and chitotriose in the absence of chbG, suggesting an additional role of chbG in the hydrolysis of chitooligosaccharides. The homologs of chbG in metazoans have been implicated in development and inflammatory diseases of the intestine, indicating that understanding the function of E. coli chbG has a broader significance.

Clinical practice : coeliac disease

Coeliac disease (CD) is an immune-mediated systemic condition elicited by gluten and related prolamines in genetically predisposed individuals and characterised by gluten-induced symptoms and signs, specific antibodies, a specific human leukocyte antigen (HLA) type and enteropathy. The risk of coeliac disease is increased in first-degree relatives, certain syndromes including Down syndrome and autoimmune disorders. It is thought to occur in 1 in 100-200 individuals, but still only one in four cases is diagnosed. Small-bowel biopsy is no longer deemed necessary in a subgroup of patients, i.e. when all of the following are present: typical symptoms or signs, high titres of and transglutaminase antibodies, endomysial antibodies, and HLA-type DQ2 or DQ8. In all other cases, small-bowel biopsy remains mandatory for a correct diagnosis. Therapy consists of a strictly gluten-free diet. This should result in complete disappearance of symptoms and of serological markers. Adequate follow-up is considered essential.

CONCLUSION

Although small-bowel biopsy may be omitted in a minority of patients, small-bowel biopsy is essential for a correct diagnosis of CD in all other cases. Diagnostic work-up should be completed before treatment with gluten-free diet instituted.

The immunopathogenesis of celiac disease reveals possible therapies beyond the gluten-free diet

Celiac disease is a T cell-mediated autoimmune inflammatory disease of the small intestine that is activated by gluten. The diagnosis of celiac disease is challenging as patients display a wide range of symptoms and some are asymptomatic. A lifelong gluten-free diet is the only currently approved treatment of celiac disease. Although the diet is safe and effective, the compliance rates and patient acceptance vary. Furthermore, many patients treated with a gluten-free diet continue to be mildly to severely symptomatic with persistent histological abnormalities, and a small number of patients develop refractory celiac disease. New therapeutic adjuncts and potential alternatives to the gluten-free diet could improve the treatment options for these patients. Advances in understanding the immunopathogenesis of celiac disease have suggested several types of therapeutic strategies that may augment or supplant the gluten-free diet. Some of these strategies attempt to decrease the immunogenicity of gluten-containing grains by manipulating the grain itself or by using oral enzymes to break down immunogenic peptides that normally remain intact during digestion. Other strategies focus on preventing the absorption of these peptides, preventing tissue transglutaminase from rendering gluten peptides more immunogenic, or inhibiting their binding to celiac disease-specific antigen-presenting molecules. Strategies that limit T cell migration to the small intestine or that reestablish mucosal homeostasis and tolerance to gluten antigens are also being explored. Additionally, it is vital to develop new therapeutic options for refractory celiac disease patients. This review highlights therapeutic strategies that may ultimately improve the health and well-being of individuals with celiac disease.