Karen A. Hunt

Shared and Distinct Genetic Variants in Type 1 Diabetes and Celiac Disease

BACKGROUND Two inflammatory disorders, type 1 diabetes and celiac disease, cosegregate in populations, suggesting a common genetic origin. Since both diseases are associated with the HLA class II genes on chromosome 6p21, we tested whether non-HLA loci are shared. METHODS We evaluated the association between type 1 diabetes and eight loci related to the risk of celiac disease by genotyping and statistical analyses of DNA samples from 8064 patients with type 1 diabetes, 9339 control subjects, and 2828 families providing 3064 parent-child trios (consisting of an affected child and both biologic parents). We also investigated 18 loci associated with type 1 diabetes in 2560 patients with celiac disease and 9339 control subjects. RESULTS Three celiac disease loci--RGS1 on chromosome 1q31, IL18RAP on chromosome 2q12, and TAGAP on chromosome 6q25--were associated with type 1 diabetes (P<1.00x10(-4)). The 32-bp insertion-deletion variant on chromosome 3p21 was newly identified as a type 1 diabetes locus (P=1.81x10(-8)) and was also associated with celiac disease, along with PTPN2 on chromosome 18p11 and CTLA4 on chromosome 2q33, bringing the total number of loci with evidence of a shared association to seven, including SH2B3 on chromosome 12q24. The effects of the IL18RAP and TAGAP alleles confer protection in type 1 diabetes and susceptibility in celiac disease. Loci with distinct effects in the two diseases included INS on chromosome 11p15, IL2RA on chromosome 10p15, and PTPN22 on chromosome 1p13 in type 1 diabetes and IL12A on 3q25 and LPP on 3q28 in celiac disease. CONCLUSIONS A genetic susceptibility to both type 1 diabetes and celiac disease shares common alleles. These data suggest that common biologic mechanisms, such as autoimmunity-related tissue damage and intolerance to dietary antigens, may be etiologic features of both diseases.

Newly identified genetic risk variants for celiac disease related to the immune response

Our genome-wide association study of celiac disease previously identified risk variants in the IL2–IL21 region. To identify additional risk variants, we genotyped 1,020 of the most strongly associated non-HLA markers in an additional 1,643 cases and 3,406 controls. Through joint analysis including the genome-wide association study data (767 cases, 1,422 controls), we identified seven previously unknown risk regions (P < 5 × 10−7). Six regions harbor genes controlling immune responses, including CCR3, IL12A, IL18RAP, RGS1, SH2B3 (nsSNP rs3184504) and TAGAP. Whole-blood IL18RAP mRNA expression correlated with IL18RAP genotype. Type 1 diabetes and celiac disease share HLA-DQ, IL2–IL21, CCR3 and SH2B3 risk regions. Thus, this extensive genome-wide association follow-up study has identified additional celiac disease risk variants in relevant biological pathways.

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

Using variants from the 1000 Genomes Project pilot European CEU dataset and data from additional resequencing studies, we densely genotyped 183 non-HLA risk loci previously associated with immune-mediated diseases in 12,041 individuals with celiac disease (cases) and 12,228 controls. We identified 13 new celiac disease risk loci reaching genome-wide significance, bringing the number of known loci (including the HLA locus) to 40. We found multiple independent association signals at over one-third of these loci, a finding that is attributable to a combination of common, low-frequency and rare genetic variants. Compared to previously available data such as those from HapMap3, our dense genotyping in a large sample collection provided a higher resolution of the pattern of linkage disequilibrium and suggested localization of many signals to finer scale regions. In particular, 29 of the 54 fine-mapped signals seemed to be localized to single genes and, in some instances, to gene regulatory elements. Altogether, we define the complex genetic architecture of the risk regions of and refine the risk signals for celiac disease, providing the next step toward uncovering the causal mechanisms of the disease.

Muramyl dipeptide and toll-like receptor sensitivity in NOD2-associated Crohn's disease

Both NOD2 (CARD15) alleles are mutated in roughly 15% of patients with Crohns disease, but functional effects are unclear. We analysed the cytokine response of peripheral blood mononuclear cells to muramyl dipeptide (MDP), the ligand for NOD2. MDP induced little TNFalpha or interleukin 1beta, but strong interleukin-8 secretion. MDP also substantially upregulated secretion of TNFalpha and interleukin 1beta induced by toll-like receptor ligands. These effects were abolished by the most common Crohns NOD2 double mutant genotypes at low nanomolar MDP concentrations, and provide the basis to develop a test of NOD2 functional deficiency. In Crohns disease, there are defects in neutrophil recruitment driven by NOD2 and interleukin 8 and in cross talk between the NOD2 and toll-like receptor pathways, which suggests that the immune system fails to receive an early priming signal.

Genome-wide analysis of allelic expression imbalance in human primary cells by high-throughput transcriptome resequencing

Many disease-associated variants identified by genome-wide association (GWA) studies are expected to regulate gene expression. Allele-specific expression (ASE) quantifies transcription from both haplotypes using individuals heterozygous at tested SNPs. We performed deep human transcriptome-wide resequencing (RNA-seq) for ASE analysis and expression quantitative trait locus discovery. We resequenced double poly(A)-selected RNA from primary CD4+ T cells (n = 4 individuals, both activated and untreated conditions) and developed tools for paired-end RNA-seq alignment and ASE analysis. We generated an average of 20 million uniquely mapping 45 base reads per sample. We obtained sufficient read depth to test 1371 unique transcripts for ASE. Multiple biases inflate the false discovery rate which we estimate to be ∼50% for random SNPs. However, after controlling for these biases and considering the subset of SNPs that pass HapMap QC, 4.6% of heterozygous SNP-sample pairs show evidence of imbalance (P < 0.001). We validated four findings by both bacterial cloning and Sanger sequencing assays. We also found convincing evidence for allelic imbalance at multiple reporter exonic SNPs in CD6 for two samples heterozygous at the multiple sclerosis-associated variant rs17824933, linking GWA findings with variation in gene expression. Finally, we show in CD4+ T cells from a further individual that high-throughput sequencing of genomic DNA and RNA-seq following enrichment for targeted gene sequences by sequence capture methods offers an unbiased means to increase the read depth for transcripts of interest, and therefore a method to investigate the regulatory role of many disease-associated genetic variants.

Coeliac disease-associated risk variants in TNFAIP3 and REL implicate altered NF-κB signalling

Objective: Our previous coeliac disease genome-wide association study (GWAS) implicated risk variants in the human leucocyte antigen (HLA) region and eight novel risk regions. To identify more coeliac disease loci, we selected 458 single nucleotide polymorphisms (SNPs) that showed more modest association in the GWAS for genotyping and analysis in four independent cohorts. Design: 458 SNPs were assayed in 1682 cases and 3258 controls from three populations (UK, Irish and Dutch). We combined the results with the original GWAS cohort (767 UK cases and 1422 controls); six SNPs showed association with p<1×10−04 and were then genotyped in an independent Italian coeliac cohort (538 cases and 593 controls). Results: We identified two novel coeliac disease risk regions: 6q23.3 (OLIG3-TNFAIP3) and 2p16.1 (REL), both of which reached genome-wide significance in the combined analysis of all 2987 cases and 5273 controls (rs2327832 p = 1.3×10−08, and rs842647 p = 5.2×10−07). We investigated the expression of these genes in the RNA isolated from biopsies and from whole blood RNA. We did not observe any changes in gene expression, nor in the correlation of genotype with gene expression. Conclusions: Both TNFAIP3 (A20, at the protein level) and REL are key mediators in the nuclear factor kappa B (NF-κB) inflammatory signalling pathway. For the first time, a role for primary heritable variation in this important biological pathway predisposing to coeliac disease has been identified. Currently, the HLA risk factors and the 10 established non-HLA risk factors explain ∼40% of the heritability of coeliac disease.

Negligible impact of rare autoimmune-locus coding-region variants on missing heritability

Genome-wide association studies (GWAS) have identified common variants of modest-effect size at hundreds of loci for common autoimmune diseases; however, a substantial fraction of heritability remains unexplained, to which rare variants may contribute. To discover rare variants and test them for association with a phenotype, most studies re-sequence a small initial sample size and then genotype the discovered variants in a larger sample set. This approach fails to analyse a large fraction of the rare variants present in the entire sample set. Here we perform simultaneous amplicon-sequencing-based variant discovery and genotyping for coding exons of 25 GWAS risk genes in 41,911 UK residents of white European origin, comprising 24,892 subjects with six autoimmune disease phenotypes and 17,019 controls, and show that rare coding-region variants at known loci have a negligible role in common autoimmune disease susceptibility. These results do not support the rare-variant synthetic genome-wide-association hypothesis (in which unobserved rare causal variants lead to association detected at common tag variants). Many known autoimmune disease risk loci contain multiple, independently associated, common and low-frequency variants, and so genes at these loci are a priori stronger candidates for harbouring rare coding-region variants than other genes. Our data indicate that the missing heritability for common autoimmune diseases may not be attributable to the rare coding-region variant portion of the allelic spectrum, but perhaps, as others have proposed, may be a result of many common-variant loci of weak effect.

Synergistic enhancement of Toll-like receptor responses by NOD1 activation

NOD1 is an intracellular pattern‐recognition receptor specific for Gram‐negative peptidoglycan that is important in host response to infections (e.g. Helicobacter pylori and Shigella flexneri). Genetic variation in NOD1 predisposes to asthma and inflammatory bowel disease. Functional responses have not previously been studied in primary human cells. NOD1 activation by low nanomolar concentrations of the specific muropeptide ligand M‐TriDAP induced minimal human peripheral blood mononuclear cell TNF‐α, IL‐1β or IL‐10 secretion, but synergistically increased Toll‐like receptor (TLR)‐induced responses. Synergistic responses were seen across multiple ligands (to TLR1/2, 2/6, 4, 5, 7/8) and a broad range of cytokine secretion (TNF‐α, IL‐1α, IL‐1β, IL‐4, IL‐6, IL‐10, GM‐CSF). Synergy was also observed in the allogeneic mixed lymphocyte reaction. These responses were similar in cells homozygous for Crohns disease‐associated NOD2 mutations. In contrast to cell lines, primary human peripheral blood mononuclear cells respond to NOD1 muropeptides at ∼ 100‐fold lower concentrations. Cross‐talk between cytosolic NOD1 and membrane‐bound TLR enhances responses to the multiple antigens simultaneously presented by a microbe.

Synergy between TLR9 and NOD2 innate immune responses is lost in genetic Crohn's disease

Background: Nucleotide binding oligomerisation domain 2 (NOD2; also known as CARD15) mutations are associated with Crohn’s disease but how mutations cause disease is poorly understood. Innate immune responses are reportedly enhanced by combined NOD2 ligand (muramyl dipeptide, MDP) and Toll-like receptor 4 ligand (TLR4, lipopolysaccharide) stimulation. Intestinal TLR signalling has a dual role—maintaining intestinal homeostasis and protection from injury as well as initiating inflammatory responses. TLR9 is functional in the intestinal epithelium where it is most strongly expressed in Paneth cells. Aims: To study possible interactions between CpG DNA (TLR9 ligand) and MDP using primary human cells of differing NOD2 genotypes. Subjects: NOD2 wild-type healthy controls (n = 7) and NOD2 homozygous Crohn’s disease patients (n = 19), age and sex matched. Methods: Peripheral blood mononuclear cells were stimulated with CpG DNA and MDP. Cytokines were measured by enzyme linked immunosorbent assay. Results: Tumour necrosis factor α (TNF-α) and interleukin 8 (IL-8) responses to CpG DNA were similar in NOD2 wild-type and homozygous mutant cells. Concomitant NOD2 stimulation had a marked synergistic effect on CpG DNA induced TNF-α responses at 10–100 ng/ml MDP. A mean 2.1-fold increase in CpG DNA induced TNF-α responses and a mean 3.7-fold increase in IL-8 responses were observed in NOD2 wild-type cells with 10 ng/ml MDP. This effect was abolished in NOD2 homozygous cells. Conclusions: NOD2 stimulation normally enhances innate immune responses to CpG DNA. This marked synergistic effect is lost in Crohn’s disease patients homozygous for NOD2 mutations, with implications for TLR mediated intestinal homeostasis and inflammation.

Complex nature of SNP genotype effects on gene expression in primary human leucocytes

BackgroundGenome wide association studies have been hugely successful in identifying disease risk variants, yet most variants do not lead to coding changes and how variants influence biological function is usually unknown.MethodsWe correlated gene expression and genetic variation in untouched primary leucocytes (n = 110) from individuals with celiac disease – a common condition with multiple risk variants identified. We compared our observations with an EBV-transformed HapMap B cell line dataset (n = 90), and performed a meta-analysis to increase power to detect non-tissue specific effects.ResultsIn celiac peripheral blood, 2,315 SNP variants influenced gene expression at 765 different transcripts (< 250 kb from SNP, at FDR = 0.05, cis expression quantitative trait loci, eQTLs). 135 of the detected SNP-probe effects (reflecting 51 unique probes) were also detected in a HapMap B cell line published dataset, all with effects in the same allelic direction. Overall gene expression differences within the two datasets predominantly explain the limited overlap in observed cis-eQTLs. Celiac associated risk variants from two regions, containing genes IL18RAP and CCR3, showed significant cis genotype-expression correlations in the peripheral blood but not in the B cell line datasets. We identified 14 genes where a SNP affected the expression of different probes within the same gene, but in opposite allelic directions. By incorporating genetic variation in co-expression analyses, functional relationships between genes can be more significantly detected.ConclusionIn conclusion, the complex nature of genotypic effects in human populations makes the use of a relevant tissue, large datasets, and analysis of different exons essential to enable the identification of the function for many genetic risk variants in common diseases.