help button home button Endocrine Society JCEM
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS

Journal of Clinical Endocrinology & Metabolism , doi:10.1210/jc.2008-0821
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Supplemental Data
Right arrow Submit a related Letter to the Editor
Right arrow Alert me when this article is cited
Right arrow Alert me when eLetters are posted
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow Request Copyright Permission
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Skinningsrud, B.
Right arrow Articles by Undlien, D. E.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Skinningsrud, B.
Right arrow Articles by Undlien, D. E.
Right arrowPubmed/NCBI databases
*Gene*GEO Profiles
*HomoloGene*Protein
*UniGene
Medline Plus Health Information
*Addison's Disease
Related Collections
Right arrow Adrenal and Hypertension
Right arrow Autoimmunity
The Journal of Clinical Endocrinology & Metabolism Vol. 93, No. 9 3310-3317
Copyright © 2008 by The Endocrine Society

Polymorphisms in CLEC16A and CIITA at 16p13 Are Associated with Primary Adrenal Insufficiency

Beate Skinningsrud, Eystein S. Husebye, Simon H. Pearce, David O. McDonald, Kristin Brandal, Anette B. Wolff, Kristian Løvås, Thore Egeland and Dag E. Undlien

Department of Medical Genetics (B.S., K.B., T.E., D.E.U.), Ullevål University Hospital, N-0407 Oslo, Norway; Institute of Medical Genetics (B.S., D.E.U.), University of Oslo, N-0315 Oslo, Norway; Section of Endocrinology (E.S.H., A.B.W., K.L.), Institute of Medicine, University of Bergen, N-5021 Bergen, Norway; Department of Medicine (E.S.H., K.L.), Haukeland University Hospital, N-5021 Bergen, Norway; and Institute of Human Genetics (S.H.P., D.O.M.), Newcastle University, Newcastle upon Tyne NE1 3BZ, United Kingdom

Address all correspondence and requests for reprints to: Beate Skinningsrud, Department of Medical Genetics, Ullevål University Hospital, University of Oslo, Kirkeveien 166, N-0407 Oslo, Norway. E-mail: beate.skinningsrud{at}medisin.uio.no.


    Abstract
 Top
 Abstract
 Introduction
 Patients and Methods
 Results
 Discussion
 References
 
Context/Objectives: It is known that different autoimmune diseases often share the same susceptibility genes. In this study we aimed to investigate if loci found associated with common autoimmune diseases in recent genome-wide association studies also could be susceptibility loci for autoimmune Addison’s disease (primary adrenal insufficiency).

Design/Patients: A total of 139 tagging single nucleotide polymorphisms (SNPs) in 11 candidate genes (IL2, IL21, IL2RA, CLEC2D, CD69, ERBB3, PTPN11, SH2B3, CLEC16A, CIITA, and PTPN2) were genotyped in a case/control study design consisting of Norwegian Addison’s disease patients (n = 332) and Norwegian healthy control individuals (n = 1029). Five SNPs were subsequently selected for analysis in a United Kingdom sample set consisting of Addison’s disease patients (n = 210) and controls (n = 191).

Results: Polymorphisms in CLEC16A and CIITA remained significantly associated with Addison’s disease in the Norwegian sample set at the 0.05 level, even after correction for multiple testing. CLEC16A and CIITA are both located at 16p13, but linkage disequilibrium patterns and logistical regression analyses suggest that SNPs in these two genes are independently associated with Addison’s disease. We were not able to confirm these associations in the United Kingdom material, however, this may well be due to the limited sample size and lack of statistical power.

Conclusion: Two alleles at 16p13 are independently associated with the risk of Addison’s disease in the Norwegian population, suggesting this chromosomal region to harbor common autoimmunity gene(s), CLEC16A and CIITA being possible independent candidates.


    Introduction
 Top
 Abstract
 Introduction
 Patients and Methods
 Results
 Discussion
 References
 
Several studies have shown that many susceptibility genes are common for several different autoimmune disorders. For example, human leukocyte antigen (HLA) class II genes, cytotoxic T lymphocyte-associated protein 4 (CTLA4), and protein tyrosine phosphatase, nonreceptor type 22 (PTPN22) alleles are associated with numerous such diseases (1). Autoimmune Addison’s disease is caused by autoimmune destruction of the adrenal cortex. A distinct feature of autoimmune Addison’s disease is the high frequency (45–70%) of coexisting nonadrenal organ-specific autoimmunity, most often endocrine diseases (2). Addison’s disease with concomitant thyroid autoimmunity and/or type 1 diabetes is the most common syndrome referred to as autoimmune polyendocrine syndrome type II (3). Addison’s disease patients’ strong proneness to autoimmunity makes the disease an interesting model disease in the search for common autoimmunity susceptibility genes. Addison’s disease is a rare autoimmune disease with a prevalence of 14 per 100,000 in Norway (4). Low prevalence has precluded the collection of large Addison’s disease cohorts, and so far no material has been presented in the literature that would have sufficient statistical power to perform genome-wide association (GWA) analysis.

We hypothesized that some recently identified susceptibility loci in other autoimmune diseases could also be involved in autoimmune Addison’s disease. To test this hypothesis, we performed a candidate gene study of selected autoimmunity associated loci. Selection of candidate genes was based on the recently reported genetic associations with type 1 diabetes, rheumatoid arthritis, and Crohn’s disease detected in a GWA scan by The Wellcome Trust Case Control Consortium (WTCCC) (5), and the follow-up study of Todd et al. (6) in type 1 diabetes and Graves’ disease. We selected 11 candidate genes (Table 1Go) for analysis in a total of 332 Norwegian Addison’s disease patients and 1029 healthy controls. Three of these regions were associated with type 1 diabetes, rheumatoid arthritis, and Crohn’s disease (10p15, 12q24, and 18q11) in the WTCCC study, indicating the presence of common autoimmunity susceptibility loci. The other selected regions represented the strongest associations with type 1 diabetes (4q27, 12p13, 12q13, and 16p13), on which 12q13 and 16p13 were also discovered in another GWA study on type 1 diabetes (7, 8). Todd et al. (6) confirmed unequivocally associations of 12q24, 12q13, 16p13, and 18p11 with type 1 diabetes, whereas the other regions all showed smaller effects. Of these, 4q27, 10p15, and 18p11 were also associated with Graves’ disease. Subsequently to the fine mapping in the Norwegian samples, five single nucleotide polymorphisms (SNPs) were selected for analysis in a smaller, statistically less well-powered United Kingdom cohort.


View this table:
[in this window]
[in a new window]

 
TABLE 1. Genes fine mapped with tagSNPs in Norwegian cases and controls

 
Our analyses demonstrated significant associations between Addison’s disease and variants in the two genes: C-type lectin domain family 16, member A, earlier named KIAA0350 (CLEC16A) and class II, major histocompatibility complex, transactivator (CIITA), both located at 16p13. Our findings add to the list of identified susceptibility genes in Addison’s disease that are shared with other autoimmune diseases, which to date includes HLA (9), PTPN22 (10, 11), CTLA4 (12, 13), and cytochrome P450, family 27, subfamily B, polypeptide 1 (14).


    Patients and Methods
 Top
 Abstract
 Introduction
 Patients and Methods
 Results
 Discussion
 References
 
Patients and controls

From Norway 332 patients with Addison’s disease were recruited from a Norwegian registry of organ-specific autoimmune diseases. A total of 1029 anonymous healthy Norwegian blood donors was used as controls. Further details on cases and controls can be found elsewhere (11). The 210 British patients were obtained at the Newcastle upon Tyne Hospitals Trusts, and surrounding district endocrine services. Details of diagnostic criteria have been reported previously (15). The 191 United Kingdom controls were recruited from the local population and had no clinical features or family history of autoimmune disease. Informed consent was given by all participants.

SNP selection

There were 11 genes in seven genomic regions selected for SNP tagging. These were: interleukin-2 (IL2); interleukin-21 (IL21); interleukin-2 receptor {alpha} chain, also known as CD25 (IL2RA); C-type lectin domain family 2 member D (CLEC2D); CD69 antigen (p60, early T-cell activation antigen) (CD69); receptor tyrosine-protein kinase erbB-3 precursor (ERBB3); protein tyrosine phosphatase, non-receptor type 11 (PTPN11); SH2B adaptor protein 3 (SH2B3); CLEC16A; CIITA; and protein tyrosine phosphatase, non-receptor type 2 (PTPN2) (Table 1Go). Tagging SNP (tagSNP) selection was based on the genotyping data from the HapMap phase II CEPH (Utah residents with ancestry from northern and western Europe) (CEU) population (16). We aimed to capture 100% of HapMap SNPs with minor allele frequency (MAF) more than 5% in the selected candidate genes, given a linkage disequilibrium (LD) cutoff at r2 > 0.8. Altogether, 150 SNPs were selected for genotyping in the first sample set of Norwegian cases and controls. Only the most associated SNPs in the five genes that showed some evidence of association in the Norwegian sample set were genotyped in the British cases and controls. These are listed in Table 2Go.


View this table:
[in this window]
[in a new window]

 
TABLE 2. SNPs analyzed in both the Norwegian and United Kingdom materials

 
When we tagged the candidate genes, we aimed to include the SNPs with the strongest associations reported in type 1 diabetes by the WTCCC (5): rs17388568 (Tenr, alias FLJ32741); rs6836189 (IL21); rs3136534 (3' of IL2); rs2104286 (IL2RA); rs3764021 (CLEC2D); rs11052552 (3' of CLEC2D); rs2292239 (ERBB3); rs11171739 (5' of ERBB3); rs3184504 (SH2B3); rs17696736 (between SH2B3 and PTPN11); rs12927773 (5' of SOCS1); rs12708716 (CLEC16A); and rs7200940 (tagging rs9746695 in CLEC16A). Some other SNPs of particular interest were also included; two SNPs in the IL2RA region [rs41295061 (in the NCBI variation database dbSNP replacing ss52580101, as it was initially named in the reference) and rs11594656] recently reported to be correlated with functionality (17); three SNPs in CLEC16A reported associated with type 1 diabetes in another GWA study (rs2903692, rs725613, and rs17673553) (7); and one SNP in the promoter of CIITA (rs3087456) previously found associated with Addison’s disease (18).

Genotyping

The genotyping of the Norwegian sample set was performed with the SNPlex Genotyping System (Applied Biosystems, Foster City, CA). Four SNPs that did not meet the requirement of SNPlex assay design were analyzed by TaqMan SNP Genotyping Assays (Applied Biosystems). The TaqMan technology was also used for genotyping the five SNPs analyzed in the British sample set (Table 2Go), and to reanalyze all samples for the rs8048002 SNP to evaluate any occurrence of genotyping errors in this low-frequency SNP, even though the plots generated by SNPlex were robust. All plots were inspected manually after automatic calling. All TaqMan SNP Genotyping Assays were predeveloped by the manufacturer (http://www.appliedbiosystems.com).

Statistics

Tests for allelic associations were based on 2 x 2 contingency tables of allele frequencies and {chi}2 tests. To correct for multiple comparisons, 10,000 permutations of the data were performed, including all 139 SNPs that were successfully genotyped. These analyses were performed using PLINK version 1.01 (19). The significance level was 0.05 after multiple testing corrections. Overall P values were calculated by treating the allele variable and population variable as covariates in logistical regression analyses. TagSNP selection, Hardy-Weinberg tests, and LD calculations were performed by Haploview versions 3.32 and 4.0 (20). Power calculations (Fig. 1Go), logistical regression analyses, and graphics were done in R version 2.6.1 (21). Forward stepwise logistical regression was performed to evaluate the independency of the associations at 16p13, whereby a SNP was added to the model provided it contributed significantly conditionally on the SNP already present. Genotypes were entered into the logistical regression model as ordered variables coded zero, one, or two, representing the number of occurrences of the risk allele per marker. The equality of odds ratios (ORs) was checked using the Breslow-Day test.


Figure 1
View larger version (7K):
[in this window]
[in a new window]

 
FIG. 1. Power calculations for allelic association of the Norwegian and United Kingdom sample sizes. The solid line represents the Norwegian sample size, whereas the dashed line represents the United Kingdom sample size. The calculations assume a SNP with risk allele frequency of 0.1, which is in complete LD with the causative variant and a significance level of 0.05.

 

    Results
 Top
 Abstract
 Introduction
 Patients and Methods
 Results
 Discussion
 References
 
Genotyping details

Four SNPs were excluded because they did not fulfill criteria for SNPlex assay design, and were untaggable by any other SNP. Seven SNPs were removed due to poor genotyping quality, which included those that deviated from Hardy-Weinberg equilibrium in controls (P < 0.01). In total we were able to analyze 139 of the 150 selected SNPs in the Norwegian sample set, which gave a 93% average coverage (r2 > 0.8) of HapMap SNPs in the 11 tagged genes (Table 1Go). The SNPlex Genotyping technology requires higher DNA quality than the TaqMan SNP Genotyping assays. A total of 63 Norwegian controls had too-low DNA quality for SNPlex analysis. After eliminating these samples from the data analyzed by the SNPlex Genotyping technology, we had an overall 94% genotyping call rate in the Norwegian and British sample sets. Two different control samples were genotyped twice on every 384-well plate. Neither of them showed any genotype discrepancies, nor were genotyping discrepancies detected when rs8048002 was reanalyzed with TaqMan.

SNPs in CLEC16A and CIITA at 16p13 are associated with Addison’s disease

SNPs in five genes showed some evidence of association with Addison’s disease in the Norwegian population, these were: CIITA and CLEC16A at 16p13, IL2RA at 10p13, CLEC2D at 12p13, and ERBB3 at 12q13. All, but IL2RA, showed a deviation in the same direction as previously reported in type 1 diabetes (5). Two of these SNPs, rs8048002 and rs12917716, located in CIITA and CLEC16A, respectively, remained significant at the 0.05 level after correcting for multiple testing of the 139 good-quality SNPs (Table 2Go).

LD at 16p13

Because CLEC16A and CIITA are located only 20 kilobases (kbs) apart, it is conceivable that these two associations reflect a common underlying association and that the association of one of the genes is secondary to LD with one "true" susceptibility gene. To address this question, we first studied the LD pattern in the region (supplemental Fig. S1, which is published as supplemental data on The Endocrine Society’s Journals Online web site at http://jcem.endojournals.org). CLEC16A resides in a 233-kb block of LD that also contains DEXI (dexamethasone-induced transcript), whereas CIITA spans at least two LD blocks. CLEC16A and CIITA display little LD in the Norwegian population, which is illustrated by the low LD between rs8048002 and rs12917716 (D' = 0.08 and r2 = 0.001). Logistical regression analysis also confirmed these two SNPs to contribute with independent effect (P = 0.0003).

In addition to rs8048002, four SNPs in CIITA showed some evidence of association (supplemental Table S1, which is published as supplemental data on The Endocrine Society’s Journals Online web site at http://jcem.endojournals.org). One of them, the promoter polymorphism of CIITA (rs3087456), previously reported associated with Addison’s disease (18) and rheumatoid arthritis (22), was in strong, but incomplete LD with rs8048002 (D' = 0.99 and r2 = 0.21). In Norwegian Addison’s disease patients, the G-risk allele of rs3087456 gave an OR of 1.23 [95% confidence interval (CI) 1.01–1.50], a risk that just reflected a LD hitchhiking effect of the C-risk allele of rs8048002, which had an OR estimate of 1.79 (95% CI 1.31–2.45). Logistical regression analysis also indicated their dependency. When the SNPs initially were chosen for tagging the CIITA gene, rs8048002 was tagging only one additional HapMap SNP, rs6498119, both located in intron 3. However, according to the HapMap data, the LD block harboring these two SNPs also spans the 12-kb upstream promoter region (23). Interestingly, three linked SNPs (rs6498126, rs11074938, and rs4781020) in CIITA with individual P values less than 0.05 were in low LD with rs8048002 (D' ≤ 0.3 and r2 ≤ 0.002) (supplemental Table S1 and supplemental Fig. S1). However, logistical regression analysis showed that they were not significantly independent of the strongest association in CLEC16A, rs12917716.

In CLEC16A, rs12917716 was tagging 31 HapMap SNPs (16) when the initial selection of tagSNPs was performed. There were 10 tagSNPs in CLEC16A, in addition to rs12917716, significantly associated with Addison’s disease before multiple testing corrections (supplemental Table S1). Our genotyping data revealed a high LD (D' > 0.9) between rs12917716 and all these SNPs, except for rs9806963 (D' = 0.37 and r2 = 0.04). These correlated SNPs also embraced the SNPs with the strongest associations with type 1 diabetes reported by the WTCCC (rs12708716) (5) and Hakonarson et al. (7) (rs2903692). However, logistical regression analysis indicated that rs9806963 might contribute with an independent effect of both rs8048002 and rs12917716 (P = 0.05).

Association analysis of the six other genotyped regions

However, none of the other loci demonstrated significant associations with Addison’s disease after multiple testing corrections. SNPs in IL2RA, CLEC2D, and ERBB3 did harbor SNPs with uncorrected significances/borderline significances (Fig. 2Go and Table 2Go). Therefore, it could very well be that SNPs in these genes are Addison’s disease susceptibility variants but that the present study is not powered to detect it. All associations, but IL2RA, went in the same direction in the Norwegian Addison’s disease patients as reported in type 1 diabetes (5). There was no evidence of association between Addison’s disease and SNPs in the other genes analyzed: IL2; IL21; CD69; SH2B3; PTPN11; and PTPN2; or the three individual SNPs each located in Tenr, 5' of SOCS1 and intergenic of SH2B3 and PTPN11 (Fig. 2Go and supplemental Table S1).


Figure 2
View larger version (15K):
[in this window]
[in a new window]

 
FIG. 2. ORs for the SNPs located in seven different loci that were analyzed in 332 Norwegian cases and 1029 controls. The vertical lines indicate 95% CIs.

 
Secondary sample set with United Kingdom origin

The estimated powers to confirm the rs8048002 and rs12917716 allelic associations in the United Kingdom sample set were 59 and 62%, respectively, given the same effect sizes and MAFs as in the Norwegian sample set. Despite this limited statistical power, we decided to analyze five SNPs in the United Kingdom sample set. None of the associations could be confirmed in the United Kingdom sample set. rs7093069 in IL2RA stood out with the most surprising allelic distribution, by showing an opposite direction of the associated allele compared with the Norwegian sample set, i.e. in the same direction as observed in type 1 diabetes (5). The other four SNPs had either an OR close to 1.00, or showed a tendency for deviation in the same direction as observed in the Norwegian material. However, after combining the two cohorts, both rs8048002 and rs12917716 showed evidence of association (Table 2Go).


    Discussion
 Top
 Abstract
 Introduction
 Patients and Methods
 Results
 Discussion
 References
 
Based on findings in GWA studies (5, 7, 8, 24, 25) in related autoimmune diseases, we used a candidate-gene approach and identified association between Addison’s disease and polymorphisms in two genes at 16p13: CIITA and CLEC16A. LD patterns and logistical regression analyses suggest that these significant associations (rs8048002 and rs12917716) are independent.

Two GWA studies (5, 7) and a follow-up study (6) have pointed out SNPs in CLEC16A to be strongly associated with type 1 diabetes. However, there was no evidence for association at this locus either with rheumatoid arthritis or Crohn’s disease in the WTCCC GWA study (5), or with Graves’ disease in the follow-up study by Todd et al. (6). CLEC16A SNPs showed risk estimates with ORs of 0.65 (7) to 0.83 (6) in type 1 diabetes, which is comparable to the effect demonstrated in our Norwegian Addison’s disease patients at rs12917716 with an OR of 0.72. Interestingly, this is the first report of an autoimmune disease besides type 1 diabetes that is associated with a SNP in CLEC16A, and is thereby suggesting the region to harbor a common autoimmunity susceptibility gene.

CLEC16A is a hot candidate of being the real susceptibility gene. The gene is classified as a C-type lectin and encodes a protein of unknown function, which is almost exclusively expressed in immune cells, such as dendritic cells, B lymphocytes, and natural killer cells (http://symatlas.gnf.org/SymAtlas/). C-type lectins are known to have several carbohydrate recognition sites, and thereby be able to discriminate between self and nonself in their involvement in the immune system. Their involvements range from extracellular pathogen recognition to pathogen neutralization inside the cells (26). However, CLEC16A seems to lack crucial domains in carbohydrate recognition, but Todd et al. (6) predicted a conserved immunoreceptor tyrosine-based activation motif encoded by exon 12. Several receptors in the immune system (also some c-type lectins) have immunoreceptor tyrosine-based activation motifs responsible for the signal transduction through binding of SH2 domain-containing proteins (27). There is an obvious necessity to obtain more functional data on the protein encoded by CLEC16A to elucidate its potential effect in immune and autoimmune pathways.

CIITA is the only known gene in the LD block harboring rs8048002. CIITA was included in this study based on its close chromosomal localization to CLEC16A, its protein function as a HLA class II transactivator (28), and the fact that a promoter polymorphism in the gene was previously reported associated with Addison’s disease (18). This promoter polymorphism (rs3087456) has also been studied in other autoimmune diseases, most carefully in rheumatoid arthritis, but with diverging conclusions. A recently published metaanalysis rejected its association with rheumatoid arthritis (29), despite that this SNP or a SNP in tight LD with it seemed to influence the expression of HLA class II genes (22). By selecting tagSNPs spanning CIITA, we identified another SNP located in intron 3, which revealed a stronger association with Addison’s disease in our study than rs3087456 (Breslow-Day test P = 0.05).

CIITA encodes the master control factor for HLA class II expression, as a coactivator interacting within a complex of three other proteins: RFXANK, RFX5, and RFXAP. The absence of HLA class II expression can be due to mutations in CIITA or any of the genes coding for these proteins, and is known as bare lymphocyte syndrome. This syndrome leads to severe immunodeficiency and is usually fatal in early childhood (30). We hypothesize that polymorphism(s) in CIITA could influence minor changes in the level or tissue selectivity of HLA class II expression, and thereby contribute as a predisposing risk factor for Addison’s disease, and possibly other related autoimmune diseases. The fact that the strongest associated SNP resides in a common LD block with SNPs covering the entire promoter region (HapMap data) may also suggest that a causative variant can affect only one of the alternative splice-variants controlled by four different promoters (23). The diverging results on associations with rheumatoid arthritis and other autoimmune disease could be due to the fact that the few studied SNPs were only in partial LD with the real causative variant. Our finding of a stronger association with the intron 3 SNP rs8048002 warrants new studies in rheumatoid arthritis and other autoimmune diseases.

For the other loci analyzed in this study, no convincing associations with Addison’s disease could be detected. However, nine of the 33 SNPs in IL2RA at 10p15 showed some evidence of association in the Norwegian sample set. However, the genotyping results from the United Kingdom sample set showed an opposite direction of the allelic association compared with Norwegian cases and controls, and thereby weakened the hypothesis of an association between IL2RA and Addison’s disease. Two other regions with interesting allelic frequency deviations in the Norwegian cases and controls were 12p13 and 12q13. LD blocks in these regions contain multiple genes. We chose CLEC2D and CD69 (at 12p13), and ERBB3 (at 12q13) as our candidate genes for fine-mapping analysis because these genes were closest to the strongest associated SNP within type 1 diabetes (5), and also because of their known function in immunity (26, 31). The 12q13 locus has been confirmed to be associated with type 1 diabetes (6, 8), whereas a follow-up on 12p13 showed weaker association (6). Anyhow, our skewed allelic distributions between cases and controls in 10p15, 12p13, and 12q13 are inconclusive, and the absence of significance may just be due to the lack of power. This highlights the importance of collecting larger materials of Addison’s disease patients.

Limitations of this study

Our study is underpowered to exclude minor effect susceptibility genes with an OR less than 1.5 (Fig. 1Go). If the studied genes have equivalent effects on Addison’s disease as reported in type 1 diabetes with ORs from 0.8–1.3 (5), we cannot exclude that some of our negative associations are due to inadequate sample sizes, rather than being truly negative. Therefore, this study is predominantly suitable to demonstrate genes with stronger effect in Addison’s disease than type 1 diabetes, as would not be unlikely with the knowledge that Addison’s disease patients are particularly prone to organ-specific autoimmunity. In association studies there is also a question if the genotyped SNPs have been able to capture the causative variant. Our high coverage rates of HapMap SNPs in the candidate genes, supported by the knowledge of a general agreement of LD patterns between the HapMap CEU cohort and European populations (32, 33), would suggest that the vast majority of common SNPs (MAF > 5%) were captured in this study. However, if rarer alleles represent the causative variant(s), deep sequencing should be an approach for future studies before any of the candidate genes are rejected as disease susceptibility factors.

The initial SNP screen was performed on the largest sample size to avoid too many false-negative results when analyzing many SNPs for a disease with low prevalence. With power estimates of 60% to reproduce the significances from the Norwegian sample set in the United Kingdom sample set, overinterpreting the lack of associations in this population might be misleading. Therefore, the results provided from the United Kingdom samples should not be seen as real replications but rather as useful contributions to any subsequent metaanalysis.

In conclusion, we have identified SNPs in CIITA and CLEC16A to be associated with autoimmune Addison’s disease. Further studies with extensive resequencing on large samples are needed to identify the causative polymorphisms at 16p13. Although the data strongly suggest the two associations to be independent, such an approach is also required to exclude the possibility that both SNPs are in LD with a single, as yet unidentified causative polymorphism. This together with functional studies of genetic variants could contribute to an increased understanding of the pathogenesis of Addison’s disease.


    Acknowledgments
 
We thank Øivind Skare, Ph.D., for valuable help.


    Footnotes
 
The study was supported by grants form South-Eastern and Western Regional Health Authorities, the Ullevål University Hospital Scientific Advisory Council (VIRUUS), and The Functional Genomics Programe (FUGE) administrated by The Research Council of Norway.

Disclosure Statement: The authors have nothing to disclose.

First Published Online July 1, 2008

Abbreviations: CEU, CEPH (Utah residents with ancestry from northern and western Europe; CI, confidence interval; CIITA, class II, major histocompatibility complex, transactivator; CLEC2D, C-type lectin domain family 2, member D; CLEC16A, C-type lectin domain family 16, member A; CTLA4, cytotoxic T lymphocyte-associated protein 4; GWA, genome-wide association; HLA, human leukocyte antigen; kb, kilobase; LD, linkage disequilibrium; MAF, minor allele frequency; OR, odds ratio; PTPN2, protein tyrosine phosphatase, nonreceptor type 2; PTPN11, protein tyrosine phosphatase, nonreceptor type 11; PTPN22, protein tyrosine phosphatase, nonreceptor type 22; SH2B3, SH2B adaptor protein 3; SNP, single nucleotide polymorphism; tagSNP, tagging single nucleotide polymorphism; WTCCC, Wellcome Trust Case Control Consortium.

Received April 16, 2008.

Accepted June 23, 2008.


    References
 Top
 Abstract
 Introduction
 Patients and Methods
 Results
 Discussion
 References
 

  1. Brand O, Gough S, Heward J 2005 HLA, CTLA-4 and PTPN22: the shared genetic master-key to autoimmunity? Expert Rev Mol Med 7:1–15[Medline]
  2. Gottlieb PA, Fain PR 2002 Genetic basis of autoimmune adrenal deficiency. Curr Opin Endocrinol Diabet 9:237–243[CrossRef]
  3. Eisenbarth GS, Gottlieb PA 2004 Autoimmune polyendocrine syndromes. N Engl J Med 350:2068–2079[Free Full Text]
  4. Lovas K, Husebye ES 2002 High prevalence and increasing incidence of Addison’s disease in western Norway. Clin Endocrinol (Oxf) 56:787–791[CrossRef][Medline]
  5. Wellcome Trust Case Control Consortium 2007 Genome-wide association study of 14,000 cases of seven common diseases and 3,000 shared controls. Nature 447:661–678[CrossRef][Medline]
  6. Todd JA, Walker NM, Cooper JD, Smyth DJ, Downes K, Plagnol V, Bailey R, Nejentsev S, Field SF, Payne F, Lowe CE, Szeszko JS, Hafler JP, Zeitels L, Yang JH, Vella A, Nutland S, Stevens HE, Schuilenburg H, Coleman G, Maisuria M, Meadows W, Smink LJ, Healy B, Burren OS, Lam AA, Ovington NR, Allen J, Adlem E, Leung HT, Wallace C, Howson JM, Guja C, Ionescu-Tîrgoviçste C, Genetics of Type 1 Diabetes in Finland, Simmonds MJ, Heward JM, Gough SC, Wellcome Trust Case Control Consortium, Dunger DB, Wicker LS, Clayton DG 2007 Robust associations of four new chromosome regions from genome-wide analyses of type 1 diabetes. Nat Genet 39:857–864[CrossRef][Medline]
  7. Hakonarson H, Grant SF, Bradfield JP, Marchand L, Kim CE, Glessner JT, Grabs R, Casalunovo T, Taback SP, Frackelton EC, Lawson ML, Robinson LJ, Skraban R, Lu Y, Chiavacci RM, Stanley CA, Kirsch SE, Rappaport EF, Orange JS, Monos DS, Devoto M, Qu HQ, Polychronakos C 2007 A genome-wide association study identifies KIAA0350 as a type 1 diabetes gene. Nature 448:591–594[CrossRef][Medline]
  8. Hakonarson H, Qu HQ, Bradfield JP, Marchand L, Kim CE, Glessner JT, Grabs R, Casalunovo T, Taback SP, Frackelton EC, Eckert AW, Annaiah K, Lawson ML, Otieno FG, Santa E, Shaner JL, Smith RM, Onyiah CC, Skraban R, Chiavacci RM, Robinson LJ, Stanley CA, Kirsch SE, Devoto M, Monos DS, Grant SF, Polychronakos C 2008 A novel susceptibility locus for type 1 diabetes on Chr12q13 identified by a genome-wide association study. Diabetes 57:1143–1146[Abstract/Free Full Text]
  9. Myhre AG, Undlien DE, Lovas K, Uhlving S, Nedrebo BG, Fougner KJ, Trovik T, Sorheim JI, Husebye ES 2002 Autoimmune adrenocortical failure in Norway autoantibodies and human leukocyte antigen class II associations related to clinical features. J Clin Endocrinol Metab 87:618–623[Abstract/Free Full Text]
  10. Velaga MR, Wilson V, Jennings CE, Owen CJ, Herington S, Donaldson PT, Ball SG, James RA, Quinton R, Perros P, Pearce SH 2004 The codon 620 tryptophan allele of the lymphoid tyrosine phosphatase (LYP) gene is a major determinant of Graves’ disease. J Clin Endocrinol Metab 89:5862–5865[Abstract/Free Full Text]
  11. Skinningsrud B, Husebye ES, Gervin K, Lovas K, Blomhoff A, Wolff AB, Kemp EH, Egeland T, Undlien DE 2008 Mutation screening of PTPN22: association of the 1858T-allele with Addison’s disease. Eur J Hum Genet 16:977–982[CrossRef][Medline]
  12. Kemp EH, Ajjan RA, Husebye ES, Peterson P, Uibo R, Imrie H, Pearce SH, Watson PF, Weetman AP 1998 A cytotoxic T lymphocyte antigen-4 (CTLA-4) gene polymorphism is associated with autoimmune Addison’s disease in English patients. Clin Endocrinol (Oxf) 49:609–613[CrossRef][Medline]
  13. Blomhoff A, Lie BA, Myhre AG, Kemp EH, Weetman AP, Akselsen HE, Huseby ES, Undlien DE 2004 Polymorphisms in the cytotoxic T lymphocyte antigen-4 gene region confer susceptibility to Addison’s disease. J Clin Endocrinol Metab 89:3474–3476[Abstract/Free Full Text]
  14. Jennings CE, Owen CJ, Wilson V, Pearce SH 2005 A haplotype of the CYP27B1 promoter is associated with autoimmune Addison’s disease but not with Graves’ disease in a UK population. J Mol Endocrinol 34:859–863[Abstract/Free Full Text]
  15. Vaidya B, Imrie H, Geatch DR, Perros P, Ball SG, Baylis PH, Carr D, Hurel SJ, James RA, Kelly WF, Kemp EH, Young ET, Weetman AP, Kendall-Taylor P, Pearce SH 2000 Association analysis of the cytotoxic T lymphocyte antigen-4 (CTLA-4) and autoimmune regulator-1 (AIRE-1) genes in sporadic autoimmune Addison’s disease. J Clin Endocrinol Metab 85:688–691[Abstract/Free Full Text]
  16. International HapMap Consortium, et al. 2007 A second generation human haplotype map of over 3.1 million SNPs. Nature 449:851–861[CrossRef][Medline]
  17. Lowe CE, Cooper JD, Brusko T, Walker NM, Smyth DJ, Bailey R, Bourget K, Plagnol V, Field S, Atkinson M, Clayton DG, Wicker LS, Todd JA 2007 Large-scale genetic fine mapping and genotype-phenotype associations implicate polymorphism in the IL2RA region in type 1 diabetes. Nat Genet 39:1074–1082[CrossRef][Medline]
  18. Ghaderi M, Gambelunghe G, Tortoioli C, Brozzetti A, Jatta K, Gharizadeh B, De BA, Pecori GF, Terzolo M, Betterle C, Falorni A 2006 MHC2TA single nucleotide polymorphism and genetic risk for autoimmune adrenal insufficiency. J Clin Endocrinol Metab 91:4107–4111[Abstract/Free Full Text]
  19. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, Maller J, Sklar P, de Bakker PI, Daly MJ, Sham PC 2007 PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet 81:559–575[CrossRef][Medline]
  20. Barrett JC, Fry B, Maller J, Daly MJ 2005 Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics 21:263–265[Abstract/Free Full Text]
  21. R Development Core Team R: 2007 A language and environment for statistical computing. Vienna: R Foundation for Statistical Computing. ISBN 3–900051-07–0. URL http://www.R-project.org
  22. Swanberg M, Lidman O, Padyukov L, Eriksson P, Akesson E, Jagodic M, Lobell A, Khademi M, Borjesson O, Lindgren CM, Lundman P, Brookes AJ, Kere J, Luthman H, Alfredsson L, Hillert J, Klareskog L, Hamsten A, Piehl F, Olsson T 2005 MHC2TA is associated with differential MHC molecule expression and susceptibility to rheumatoid arthritis, multiple sclerosis and myocardial infarction. Nat Genet 37:486–494[CrossRef][Medline]
  23. Reith W, Mach B 2001 The bare lymphocyte syndrome and the regulation of MHC expression. Annu Rev Immunol 19:331–373[CrossRef][Medline]
  24. Hafler DA, Compston A, Sawcer S, Lander ES, Daly MJ, De Jager PL, de Bakker PI, Gabriel SB, Mirel DB, Ivinson AJ, Pericak-Vance MA, Gregory SG, Rioux JD, McCauley JL, Haines JL, Barcellos LF, Cree B, Oksenberg JR, Hauser SL 2007 Risk alleles for multiple sclerosis identified by a genomewide study. N Engl J Med 357:851–862[Abstract/Free Full Text]
  25. van Heel DA, Franke L, Hunt KA, Gwilliam R, Zhernakova A, Inouye M, Wapenaar MC, Barnardo MC, Bethel G, Holmes GK, Feighery C, Jewell D, Kelleher D, Kumar P, Travis S, Walters JR, Sanders DS, Howdle P, Swift J, Playford RJ, McLaren WM, Mearin ML, Mulder CJ, McManus R, McGinnis R, Cardon LR, Deloukas P, Wijmenga C 2007 A genome-wide association study for celiac disease identifies risk variants in the region harboring IL2 and IL21. Nat Genet 39:827–829[CrossRef][Medline]
  26. Zelensky AN, Gready JE 2005 The C-type lectin-like domain superfamily. FEBS J 272:6179–6217[CrossRef][Medline]
  27. Underhill DM, Goodridge HS 2007 The many faces of ITAMs. Trends Immunol 28:66–73[CrossRef][Medline]
  28. Reith W, LeibundGut-Landmann S, Waldburger JM 2005 Regulation of MHC class II gene expression by the class II transactivator. Nat Rev Immunol 5:793–806[CrossRef][Medline]
  29. Bronson PG, Criswell LA, Barcellos LF 2008 The MHC2TA –168A/G polymorphism and risk for rheumatoid arthritis: a meta-analysis of 6861 patients and 9270 controls reveals no evidence for association. Ann Rheum Dis 67:933–936[Abstract/Free Full Text]
  30. Waldburger JM, Masternak K, Muhlethaler-Mottet A, Villard J, Peretti M, Landmann S, Reith W 2000 Lessons from the bare lymphocyte syndrome: molecular mechanisms regulating MHC class II expression. Immunol Rev 178:148–165[CrossRef][Medline]
  31. Jones RB, Gordus A, Krall JA, MacBeath G 2006 A quantitative protein interaction network for the ErbB receptors using protein microarrays. Nature 439:168–174[CrossRef][Medline]
  32. Mueller JC, Lohmussaar E, Magi R, Remm M, Bettecken T, Lichtner P, Biskup S, Illig T, Pfeufer A, Luedemann J, Schreiber S, Pramstaller P, Pichler I, Romeo G, Gaddi A, Testa A, Wichmann HE, Metspalu A, Meitinger T 2005 Linkage disequilibrium patterns and tagSNP transferability among European populations. Am J Hum Genet 76:387–398[CrossRef][Medline]
  33. Smith EM, Wang X, Littrell J, Eckert J, Cole R, Kissebah AH, Olivier M 2006 Comparison of linkage disequilibrium patterns between the HapMap CEPH samples and a family-based cohort of Northern European descent. Genomics 88:407–414[CrossRef][Medline]



This article has been cited by other articles:


Home page
J. Clin. Endocrinol. Metab.Home page
W. Arlt
The Approach to the Adult with Newly Diagnosed Adrenal Insufficiency
J. Clin. Endocrinol. Metab., April 1, 2009; 94(4): 1059 - 1067.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Supplemental Data
Right arrow Submit a related Letter to the Editor
Right arrow Alert me when this article is cited
Right arrow Alert me when eLetters are posted
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow Request Copyright Permission
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Skinningsrud, B.
Right arrow Articles by Undlien, D. E.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Skinningsrud, B.
Right arrow Articles by Undlien, D. E.
Right arrowPubmed/NCBI databases
*Gene*GEO Profiles
*HomoloGene*Protein
*UniGene
Medline Plus Health Information
*Addison's Disease
Related Collections
Right arrow Adrenal and Hypertension
Right arrow Autoimmunity


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Endocrinology Endocrine Reviews J. Clin. End. & Metab.
Molecular Endocrinology Recent Prog. Horm. Res. All Endocrine Journals