CYP21A2 polymorphisms in patients with autoimmune Addison s disease, and linkage disequilibrium to HLA risk alleles

Similar documents
HLA Antigen Expression in Autoimmune Endocrinopathies

BDC Keystone Genetics Type 1 Diabetes. Immunology of diabetes book with Teaching Slides

HLA-DR, HLA-DQB1 and PTPN22 gene polymorphism: association with age at onset for autoimmune diabetes

ASSESSMENT OF THE RISK FOR TYPE 1 DIABETES MELLITUS CONFERRED BY HLA CLASS II GENES. Irina Durbală

Lack of association of IL-2RA and IL-2RB polymorphisms with rheumatoid arthritis in a Han Chinese population

Handling Immunogenetic Data Managing and Validating HLA Data

HLA Complex Genetics & Biology

CYP21A2 Mutations Found in Congenital Adrenal Hyperplasia Patients in the California Population

ADDISON S disease is the consequence of the destruction,

ZHANG Qin, YANG Yun-mei, LV Xue-ying

SALSA MLPA KIT P050-B2 CAH

Normal overall mortality rate in Addison s disease, but young patients are at risk of premature death

Advances in Genetics Endocrine Research. Site: The study was conducted at a university referral center.

An association analysis of the HLA gene region in latent autoimmune diabetes in adults

Anti 21 - Hydroxylase Assay RIA

Completing the CIBMTR Confirmation of HLA Typing Form (Form 2005)

Early onset of polyglandular failure is associated with HLA-DRB1*03

Award Number: W81XWH TITLE: CYP1B1 Polymorphism as a Risk Factor for Race-Related Prostate Cancer

Whole-genome detection of disease-associated deletions or excess homozygosity in a case control study of rheumatoid arthritis

Investigation of Programmed Cell Death-1 (PD-1) Gene Variations at Positions PD1.3 and PD1.5 in Iranian Patients with Non-small Cell Lung Cancer

Antigen Presentation to T lymphocytes

Part XI Type 1 Diabetes

Diversity and Frequencies of HLA Class I and Class II Genes of an East African Population

Association between atopic dermatitis-related single nucleotide polymorphisms rs and psoriasis vulgaris in a southern Chinese cohort

ANALYSIS OF IL17 AND IL17RA POLYMORPHISMS IN SPANISH PSORIASIS PATIENTS: ASSOCIATION WITH RISK FOR DISEASE.

Drug Metabolism Disposition

Department of Pediatrics, Kawasaki Medical School, Kurashiki , Japan 3)

Myoglobin A79G polymorphism association with exercise-induced skeletal muscle damage

Relationship Between HLA-DMA, DMB Alleles and Type 1 Diabetes in Chinese

The Human Major Histocompatibility Complex

Results. Introduction

Endocrine autoimmune disease: genetics become complex

THE ASSOCIATION OF SINGLE NUCLEOTIDE POLYMORPHISMS IN INTRONIC REGIONS OF ISLET CELL AUTOANTIGEN 1 AND TYPE 1 DIABETES MELLITUS.

Topic (Final-03): Immunologic Tolerance and Autoimmunity-Part II

Functional studies of novel mutations in CYP21A2 detected. in Norwegian patients with congenital adrenal hyperplasia

Endocrine Surgery. Characteristics of the Germline MEN1 Mutations in Korea: A Literature Review ORIGINAL ARTICLE. The Korean Journal of INTRODUCTION

Influence of interleukin-18 gene polymorphisms on acute pancreatitis susceptibility in a Chinese population

Profiling HLA motifs by large scale peptide sequencing Agilent Innovators Tour David K. Crockett ARUP Laboratories February 10, 2009

the HLA complex Hanna Mustaniemi,

SALSA MLPA probemix P241-D2 MODY mix 1 Lot D As compared to version D1 (lot D1-0911), one reference probe has been replaced.

CYP19 gene polymorphisms and the susceptibility to breast cancer in Xinjiang Uigur women

Effects of age-at-diagnosis and duration of diabetes on GADA and IA-2A positivity

Association between interleukin-17a polymorphism and coronary artery disease susceptibility in the Chinese Han population

Research: Genetics HLA class II gene associations in African American Type 1 diabetes reveal a protective HLA-DRB1*03 haplotype

DEFINITION OF HIGH RISK TYPE 1 DIABETES HLA-DR AND HLA- DQ TYPES USING ONLY THREE SINGLE NUCLEOTIDE POLYMORPHISMS

Association between the -77T>C polymorphism in the DNA repair gene XRCC1 and lung cancer risk

HLA-A*26 and Susceptibility of Iranian Patients with Non-Hodgkin Lymphoma

SHORT COMMUNICATION. K. Lukacs & N. Hosszufalusi & E. Dinya & M. Bakacs & L. Madacsy & P. Panczel

Development of a Multiplex CYP21A2 Genotyping Assay for Congenital Adrenal Hyperplasia Screening

Association between the CYP11B2 gene 344T>C polymorphism and coronary artery disease: a meta-analysis

To test the possible source of the HBV infection outside the study family, we searched the Genbank

BST227 Introduction to Statistical Genetics. Lecture 4: Introduction to linkage and association analysis

Breast Cancer Susceptibility and HLA Complex: Source of Missing Heritability?

TITLE: CYP1B1 Polymorphism as a Risk Factor for Race-Related Prostate Cancer

Lack of Association between Endoplasmic Reticulum Stress Response Genes and Suicidal Victims

Genetics and Genomics in Medicine Chapter 8 Questions

Sex-specific association of PTPN T with type 1 diabetes but not with Hashimoto s thyroiditis or Addison s disease in the German population

From DEPARTMENT OF MOLECULAR MEDICINE AND SURGERY Karolinska Institutet, Stockholm, Sweden GENETIC ANALYSIS OF TYPE 1 DIABETES.

A rapid screening for steroid 21-hydroxylase mutations in patients with congenital adrenal hyperplasia

Lack of association between rare mutations of the SIAE gene and rheumatoid arthritis in a Han Chinese population

SALSA MLPA probemix P241-D2 MODY mix 1 Lot D2-0716, D As compared to version D1 (lot D1-0911), one reference probe has been replaced.

Relationship between vitamin D (1,25-dihydroxyvitamin D3) receptor gene polymorphisms and primary biliary cirrhosis risk: a meta-analysis

Association of MHC class III gene polymorphisms with ER-positive breast cancer in a Chinese Han population

Allelic and Haplotype Frequencies of the p53 Polymorphisms in Brain Tumor Patients

Original Article. C18orf1 located on chromosome 18p11.2 may confer susceptibility to schizophrenia

Genetics and Genomics in Endocrinology

Cover Page. The handle holds various files of this Leiden University dissertation.

Genome-wide Association Analysis Applied to Asthma-Susceptibility Gene. McCaw, Z., Wu, W., Hsiao, S., McKhann, A., Tracy, S.

Polymorphism of the PAI-1gene (4G/5G) may be linked with Polycystic Ovary Syndrome and associated pregnancy disorders in South Indian Women

Original Article The programmed death-1 gene polymorphism (PD-1.5 C/T) is associated with non-small cell lung cancer risk in a Chinese Han population

HLA and antigen presentation. Department of Immunology Charles University, 2nd Medical School University Hospital Motol

Proposal form for the evaluation of a genetic test for NHS Service Gene Dossier

Insulin Resistance. Biol 405 Molecular Medicine

Cover Page. The handle holds various files of this Leiden University dissertation.

Retrospective Genetic Analysis of Efficacy and Adverse Events in a Rheumatoid Arthritis Population Treated with Methotrexate and Anti-TNF-α

Detection of adrenocortical autoantibodies in Addison s disease with a peroxidaselabelled protein A technique

HLA Class II Variants in Chinese Breast Cancer Patients

HLA and antigen presentation. Department of Immunology Charles University, 2nd Medical School University Hospital Motol

IL10 rs polymorphism is associated with liver cirrhosis and chronic hepatitis B

Graves' disease and HLA association

Autoimmune diseases. Autoimmune diseases. Autoantibodies. Autoimmune diseases relatively common

The major histocompatibility complex (MHC) is a group of genes that governs tumor and tissue transplantation between individuals of a species.

Cover Page. The handle holds various files of this Leiden University dissertation.

Association between ERCC1 and ERCC2 polymorphisms and breast cancer risk in a Chinese population

Immunology - Lecture 2 Adaptive Immune System 1

2) Cases and controls were genotyped on different platforms. The comparability of the platforms should be discussed.

SHORT COMMUNICATION Insulin gene VNTR genotype associates with frequency and phenotype of the autoimmune response to proinsulin.

California Association for Medical Laboratory Technology. Distance Learning Program

Supplementary information

Patients: Adult KPD patients (n 384) were followed longitudinally in a research clinic.

During the hyperinsulinemic-euglycemic clamp [1], a priming dose of human insulin (Novolin,

Allele and Haplotype Frequencies of Human Leukocyte Antigen-A, -B, -C, -DRB1, and -DQB1 From Sequence- Based DNA Typing Data in Koreans

Protein tyrosine phosphatase 1B is not a major susceptibility gene for type 2 diabetes mellitus or obesity among Pima Indians

Deliverable 2.1 List of relevant genetic variants for pre-emptive PGx testing

GENOME-WIDE ASSOCIATION STUDIES

Immunological Tolerance

Table 1 Functional polymorphisms identified by XGEN group, Center for Pharmacogenomics in OSU College of Medicine.

Belgian Thyroid Club

HLA Mismatches. Professor Steven GE Marsh. Anthony Nolan Research Institute EBMT Anthony Nolan Research Institute

Influence of the c.1517g>c genetic variant in the XRCC1 gene on pancreatic cancer susceptibility in a Chinese population

Transcription:

Page of Accepted Preprint first posted on September 0 as Manuscript EJE--0 CYPA polymorphisms in patients with autoimmune Addison s disease, and linkage disequilibrium to HLA risk alleles Ingeborg Brønstad, Beate Skinningsrud, Eirik Bratland, Kristian Løvås,, Dag Undlien,, Eystein Sverre Husebye, & Anette Susanne Bøe Wolff 0 Department of Clinical Science, University of Bergen, 0 Bergen, Norway Department of Medical Genetics, Oslo University Hospital, 00 Oslo, Norway Department of Medicine, Haukeland University Hospital, 0 Bergen, Norway Institute of Medical Genetics, University of Oslo, 0 Oslo, Norway Corresponding author: Ingeborg Brønstad Department of Clinical Science, University of Bergen Laboratory building, th floor 0 Bergen, Norway e-mail: Ingeborg.Bronstad@k.uib.no Telephone: + 0 0 Short title: CYPA in Autoimmune Addison s disease Keywords: CYPA; Autoimmune Addison s disease; Human Leukocyte Antigen; -Hydroxylase Word count: 0 Copyright 0 European Society of Endocrinology.

Page of 0 0 Abstract Objective: Steroid -hydroxylase (OH), encoded by CYPA, is the major autoantigen in autoimmune Addison s disease (AAD). CYPA is located in the region of the human leukocyte antigen (HLA) complex on chromosome p., which harbours several risk alleles for AAD. The objective was to investigate whether CYPA gene variants confer risk of AAD independently of other risk alleles in the HLA loci. Design: DNA samples from Norwegian patients with AAD and 0 healthy controls (HC) previously genotyped for the HLA-A, -B, -DRB, and -DQB and MICA loci were used for genotyping of CYPA. Methods: Genotyping of CYPA was performed by direct sequencing. Linkage of CYPA to the HLA loci was assessed using UNPHASED version.0.0 and PHASE version.. Results: Heterozygotes of the single nucleotide polymorphisms (SNPs) rs, rs, rs, rs and rs were detected significantly more frequently in AAD patients compared to HC (P < 0.00), but all SNPs were in linkage disequilibrium (LD) with high-risk HLA-DRB haplotypes. rsc protected against AAD (Odds ratio = 0., % confidence interval [0.0 0.0], P =. e-0). This SNP was not in LD with HLA loci (P = 0.0), but did not increase protection when considering the effect of HLA-DRB alleles. Mutations causing congenital adrenal hyperplasia were found in heterozygosity in <. % of the cases in both groups. Conclusion: Genotypes of CYPA associated with AAD are in LD with the main autoimmune AAD risk loci HLA-DRB and do not constitute an independent genetic susceptibility locus.

Page of 0 Introduction Autoimmune Addison s disease (AAD) is caused by the destruction of hormone producing cells in the adrenal cortex by autoreactive immunological mechanisms. AAD often occurs together with other organ specific autoimmune disorders, as part of autoimmune polyendocrine syndromes types and (APS- and-, respectively),. APS- is a rare monogenic autosomal recessive disease caused by mutations in the autoimmune regulator (AIRE) gene. Isolated AAD and APS-, however, are thought to be caused by susceptibility variants of multiple genes interacting with environmental triggers. To date, the human leukocyte antigen (HLA) class II haplotypes DRB*0:0 and DRB*0:0 have been shown to be the strongest predisposing genetic factors for AAD with odds ratio (OR) of. and., respectively. Risk is particularly high when these haplotypes are combined (OR = ). The HLA class I genes HLA-A and B (OR =. of HLA-B*0) and MHC-class I related chain A (MICA) also harbour AAD risk alleles -. In addition, several other genes related to the immune system have been associated with AAD, i.e. cytotoxic T lymphocyte antigen- (CTLA-),, protein tyrosine phosphatase non-receptor type (PTPN) 0,, MHC class II transactivator (CIITA), C-type lectin domain family, member A (CLECA), cytochrome P0, family, subfamily B, polypeptide (CYPB ), and programmed death ligand (PD-L). 0 The major autoantigen in AAD is the enzyme steroid -hydroxylase (OH), which is expressed mainly in the adrenal cortex, and circulating autoantibodies (Ab) against OH are present in more than % of patients with AAD in various European cohorts, -. OH is encoded by CYPA, and mutations in this gene are the main cause of congenital adrenal hyperplasia (CAH) 0. Polymorphisms and mutations in CYPA are common,

Page of mainly due to the presence of the structurally related pseudogene CYPAP. Variants in CYPA may potentially have an impact on the immunogenicity of OH and thereby influence the chance of an autoreactive response being directed towards it. Examples of such effects are seen in type diabetes where autoantibody frequencies towards the zinc transporter eight (SLC0A) is dependent on which isoform is expressed. Moreover, different isoforms of the thyroid stimulating hormone receptor (TSHR) vary in their immunogenicity,. The expression of CYPA in thymic medullary cells could potentially be dependent on polymorphisms thereby influencing negative selection of OH-specific T cells. 0 We therefore asked if CYPA is a susceptibility gene for AAD. Previously, two studies have addressed CYPA polymorphisms in relation to risk of developing AAD without finding any disease-specific variants,. However, both studies were underpowered; the largest cohort included patients with AAD. In the present study we analysed the much larger well-characterised Norwegian cohort (n = ). Since CYPA is located on chromosome p. in the HLA-III locus, which is in strong linkage disequilibrium (LD) with the HLA-I and II loci, such a large sample size is required to determine whether potential CYPA variants add to the risk of AAD independently of the HLA loci. 0 Subjects and methods Subjects A total of patients with AAD (% females, % males; mean age years) recruited from the National Norwegian Registry of patients with Addison s disease and 0 anonymous Norwegian blood donors controls (HC) were available for genetic analyses. OH

Page of autoantibody index values as determined by radioimmunoassay (RIA) were available for all patients.aad was diagnosed in patients with adrenal insufficiency with either a positive OH-Ab test, or adrenal insufficiency combined with another endocrine autoimmune disease. Of the AAD patients, % were positive for OH-Ab, and % had APS- (Addison s disease plus autoimmune thyroid disease and/or type diabetes), while % had isolated AAD. 0 Ethics All included patients and blood donors signed a written consent form. The study was approved by the Regional Committee for Medical Ethics of Western Norway, and performed according to the Helsinki Declaration. 0 Genetic analyses DNA was isolated from whole blood using standard commercial kits. The sequence of the whole CYPA gene, including 0 exons and introns, were analysed in patients with AAD. Since the antigenicity of OH most likely would be located in the coding regions, coding polymorphisms and rare variants detected in the patients were sought for in a larger material of AAD patients (n = ) and HC (n = 0). The CYPA gene was amplified in two different fragments using primers discriminating the pseudogene from the functional CYPA gene. The primers specific for CYPA were targeted to the base pair (bp) deletion region in exon and the exon cluster region, which are the major markers for distinguishing the pseudogene from CYPA. The first primer pair amplified a fragment starting with the promoter region and ending with the point mutation region in exon (fragment ); the second primer pair amplified the gene from exon spanning the bp

Page of 0 deletion region to 0 bp of the end from exon 0 (fragment ). Expand High Fidelity PCR system (Roche; Basel, Switzerland) was used for amplifying PCR fragment ; the thermal cycling protocol consisted of an initial denaturation step at ºC for minutes, followed by 0 cycles of denaturation at ºC for seconds, annealing at ºC for 0 seconds and elongation at ºC for. minutes. Then the cycle was repeated 0 times with an increment of seconds for each cycle. The PCR was terminated by an elongation step at ºC for minutes. AmpliTaq Gold with Gene Amp PCR system (Life Technologies, Waltham, Massachusetts, USA ) was used for amplifying PCR fragment ; the thermal cycling conditions were initial denaturation at ºC for 0 min, followed by cycles of ºC for 0 sec, ºC for 0 sec and ºC for. min. The final elongation step was run at ºC for min. Variants and single nucleotide polymorphisms (SNP) of CYPA were identified by direct sequencing using standard conditions on a ABI 0 DNA-Analyzer. Primers (purchased from Eurogentec; Seraing, Liège, Belgium) of PCR and sequencing reactions have been described previously -. Sequencing data were aligned to the CYPA gene sequence (http://www.ncbi.nlm.nih.gov/nuccore/m.). As a quality control to check that the pseudogene was not amplified in the PCR product, we checked that no pseudogene variants were present in the sequences. Copy number of the CYPA gene was determined in AAD patients and HC by duplex real time PCR method 0. The genotyping of HLA-A, -B, -DRB and DQB, and MICA was described previously by. 0 Statistical analyses Chi-square and Fisher s exact tests for determination of significant differences between genotype and allele frequencies of patients and HC, grouped according to age at diagnosis, patients with or without type diabetes or thyroid disease and gender; and calculations of OR with % confidence intervals (CI) were performed by IBM SPSS Statistics 0 using

Page of 0 significance level of P < 0.0. Calculation of Hardy-Weinberg equilibrium (p > 0.00 cutoff) and LD map of the CYPA gene SNPs were conducted by Haploview version.. Global association tests with conditional analysis for allele main effects were performed by UNPHASED version.0.0 to determine if associations with CYPA were independent (P < 0.0) of the HLA loci or not. Haplotype frequencies were generated by PHASE version.. D Agostino and Pearson omnibus normality test was used to test whether the values of OH index and age at diagnosis were distributed in a Gaussian manner; and Kruskal-Wallis test with Dunn s multiple comparison was used for testing differences between CYPA/HLA genotypes or haplotypes and OH antibody index values and mean age at diagnosis using the Graph Pad Prism version.0. Results 0 SNPs and haplotypes of CYPA in AAD patients and healthy controls The common exon variants p._0insl/g.0_0insctg (rs), p.k0r/g.a>g (rs), p.ns/g.a>g (rs), the intron variant g.c>a (rs), and the intron variant g. C>T (rs), were frequently detected in both AAD and HC (Table ). Conversely, p.st/g. G>C (rs) appeared rarely (Table ). Several other silent and intron variants were also detected in both AAD patients and HC (listed in Supplementary Table ), but were not further investigated in the large cohort. For the rs, rs and rs SNPs, the heterozygote genotype occurred significantly (P 0.00) more often in patients with AAD than in HC, compared to the

Page of homozygote genotype of the major alleles (Table ). Haplotype reconstruction of CYPA revealed that AA and +GG (in the order rs/rs/rs, from now denoted as the AA and +GG haplotypes) mainly occurred together (Supplementary Table ), with a strong LD (r > 0., D > 0., Fig. a and b). Frequencies of these haplotypes were significantly different between AAD patients and HC (P = 0.00), where. % AAD patients and. % HC had the haplotype AA, and. % AAD patients and. % HC had the haplotype +GG. Other combinations of these SNPs were seen in. % of AAD patients and.% of HC. 0 Regarding the intron variants rs and rs, the heterozygote genotype and the minor alleles were also more common in AAD than in HC (Table ). The frequency of rsc was.% in the AAD group, compared to.% in HC (OR =., [. -.], P = 0.00). The allele frequencies of rst were.% in AAD and.% in HC (OR =. [.0 -.], P = 0.00). The majority of the rst alleles were seen together with the +GG haplotype, but never with the AA haplotype (Supplementary Table ), and was almost in complete LD with all these three coding SNPs (D > 0., Fig. b). 0 The rsc seemed to be a protective variant to AAD, where the allele frequency of the C allele was.% in HC compared to.% in AAD patients (OR = 0., [0.0 0.0], P =. e-0, Table ). The rsc allele was only seen together with the +GG variant of the rs/rs/rs haplotype (Supplementary Table ), which was in complete LD (D =, Fig. b). The rsc allele in the AAD patients was mainly seen ( of 0 alleles) in the haplotype +CGTCG (in order rs/rs/rs/rs/rs/rs). In the control group, the haplotypes +AGTCG and +AGCCG were more frequent ( and of rsc alleles, respectively, Supplementary Table ).

Page of 0 Are the CYPA associations independent of HLA I and II loci? To this end a global regression analysis was performed to determine if AAD associations of CYPA were independent (P < 0.0) of high risk HLA loci. Table shows that all the CYPA SNPs associations, except rs (P = 0.0), were lost when HLA-DRB was set as conditional marker, suggesting that the HLA-type determines the AAD-risk. The SNPs rs, rs and rs showed independence from the HLA-DQB locus (P = 0.0, P = 0.000 and P = 0.0, respectively), while the rs, rs and rs SNPs were in LD with HLA-DQB (Table ). All the exon variants of CYPA further showed independency to the HLA-A, HLA-B and MICA while the two intron variants were both in LD with HLA-A and MICA. The rs was in LD with HLA-B, whereas rs was not (P = 0.0) (Table ). 0 We next wanted to investigate how CYPA alleles are linked to the HLA-DRB risk alleles. Supplementary Table shows the major CYPA haplotypes (the SNPs rs/rs/rs/rs with haplotype AGA and +GGG) in combination with the most common HLA-DRB haplotypes generated by PHASE. The HLA-DRB*00 mainly occurred together with the +GGG haplotype (i.e.. % and. % of AAD and HC, respectively), while the HLA-DRB*0:0 haplotype was most commonly seen with the AGA variant (. % and. % of AAD and HC cases, respectively (Supplementary Table ). The HLA-DRB* and *0 haplotypes also occurred most frequently together with the +GGG CYPA haplotype, while the HLA-DRB*0:0 was most frequently seen together with AGA (Supplementary Table ). Furthermore, the individuals with heterozygous HLA- DRB*0:0/*0:0 high-risk genotype, were also heterozygote for CYPA AA/+GG in. % AAD patients (n = ) and. % HC (n = ), respectively.

Page 0 of For rs, the HLA-DRB*0:0 haplotype was most frequently reconstructed together with the minor C allele, while the majority of HLA-DRB*0:0, -* and -*0 haplotypes were seen together with the A allele. For HLA-DRB*0:0 both rs alleles were more equally distributed (Supplementary Table ). For rs SNP, the minor T allele was mainly seen together with the HLA-DRB*0:0 haplotype (Supplementary Table ). 0 The majority of the rsc alleles in AAD ( of 0) were reconstructed by PHASE together with the high-risk HLA-DRB*0:0, while the rest with HLA-DRB*0:0 ( of 0) and HLA-DRB*:0 ( of 0). For the HC group, the rsc allele was mainly seen together with the haplotypes HLA-DRB*:0 ( of ),-*0 ( of ) and -*0:0 ( of ). The protective effect of rsc was lost (p > 0.0) when considering the effect of the HLA-DRB haplotypes (Table ). 0 CYPA variants related to biochemical and clinical parameters Although all of the patients in the present study were diagnosed with AAD, % lacked OH-Ab at the time of testing. Our analyses revealed that the frequencies of heterozygous variants of CYPA were higher in the OH-Ab positive group (Supplementary Table ). The mean OH-Ab levels at diagnosis were also higher in patients with the AAD associated genotypes of CYPA, although not statistically significant (Supplementary Table ). Furthermore, the HLA-DRB*0:0 allele was found in. % and. % of OH-Ab positives and negatives, respectively. Likewise, for HLA-DRB*0:0, the allele frequency was. % and. % for OH-Ab positives and negatives. HLA-DRB*0:0 and -* were more common in the OH-Ab negative group at % and 0. %, respectively, compared to. % and. % in the OH-Ab positive group. Finally, the HLA- 0

Page of DRB*0:0/*0:0 genotype was found only in OH-Ab positives. The mean OH-Ab levels at diagnosis were significantly higher in patients with the HLA-DRB*0:0 and - *0:0 haplotypes compared to other haplotypes, and similarly for the HLA- DRB*0:0/*0:0 genotype (p < 0.000). Patients with the HLA-B*0 had significantly higher OH Ab indices at diagnoses compared to patients with other HLA-B haplotypes (p = 0.0), while the OH Ab indices were significantly higher in patients with the HLA-A*0 compared to HLA-A*0 (p = 0.00). The results are summarized in Supplementary Table. 0 The AAD associated CYPA variants were more frequent in patients with onset of the disease before 0 years age than in those with later onset. In addition, the protective variant rsc was not detected in patients with disease onset before 0 years (Supplementary Table ). There was no association between autoimmune thyroid disease of the APS - patients and the CYPA variants. The rst allele was, however, more frequent in APS- patients with type diabetes than in AAD patients without type diabetes (OR =., [.0 -.0], P = 0.0), but there was no such association to the other CYPA variants. There was no association between gender and CYPA variants. 0 Rare variants of CYPA in AAD patients and healthy controls A rare heterozygous coding variant, p. AT/g.G>A, was found in /0 patients with AAD and in /0 HC. All four were also heterozygous for the rare haplotype HLA-B*0 which was absent in other patients and controls. Known disease causing mutations frequently reported in patients with CAH, were also found as heterozygous in both AAD patients and HC: The I splice variant (rsg) was detected in /0 AAD patients and / HC. The

Page of base pair deletion in exon (rs00) was found in one AAD patient and in one HC. The non-classical CAH variant p.vl (rs) was found in /0 AAD patients and / HC, while p.ps (rs) was found in / AAD patients and / HC. No copy number variations were detected in AAD patients or HC; they all had two copies of CYPA. 0 Discussion It is not known why certain proteins, such as OH become autoantigens in organ-specific autoimmune diseases. However, common features are that autoantigens are intracellular enzymes with restricted tissue expression. One possibility is that only altered or modified selfproteins can induce autoreactive immune responses, and that genetic polymorphisms may represent such alterations. By direct sequencing of CYPA in AAD patients we did not find any disease specific risk alleles, but heterozygous genotypes of common CYPA variants were significantly more frequent in AAD patients than in HC. However, we found that all the CYPA SNPs associations, except for rs, were dependent on the HLA-DRB locus, which is the strongest single genetic risk factor of AAD,,,. Therefore, the higher heterozygous genotype frequency of the two main haplotypes of CYPA alleles, AA and +GG, merely reflects the high risk genotype of HLA-DRB*0:0/*0:0. The protective rsc association could also be explained by the linkage to HLA protective variants. 0 The non-coding AAD-associated CYPA variants rs and rs were also dependent of MICA and HLA-A or B. The rst allele was further associated with type diabetes of the APS- patients included in the study. Interestingly, the rsc and

Page of rst alleles have both shown to be in LD with components of the. ancestral haplotype. 0 Our data showed that all of our patients with the DRB*0:0/*0:0 genotype were positive for OH Ab and had the highest OH Ab indices compared to other genotypes. This is consistent of HLA-DRB*0:0/*0:0 genotype association with expression of OH Ab. The corresponding CYPA genotype did not add to the high OH Ab association of HLA- DRB*0:0/*0:0, which underlines the strength of HLA-DRB risk factors. Similarly, we here found a trend for association betweencypa genotypes and age of disease onset. This is in agreement with previous findings of a significant lower mean age of onset in patients who carries the HLA-DRB*0:0/*0:0 genotype. 0 A hypothesis proposes that low protein expression due to genetic polymorphisms might increase the probability of autoreactive T cells to escape from the thymus during negative selection. This has been suggested for the variable number of tandem repeats polymorphism in the gene encoding insulin and its association to production of Ab and autoreactive T cells to insulin in type diabetes, 0. Although the HLA-DRB locus proves to be the highest predisposing factor of AAD, the significance of CYPA heterozygote genotypes is not clear since it is not known if the gene is inherited co-dominantly like the HLA genes. In conclusion, sequencing of the exons and introns of CYPA did not reveal any disease specific risk variants associated with AAD. Instead, CYPA polymorphisms are in LD with the high risk haplotypes of the HLA-DRB loci, and do not independently add to increased risk or protection, but merely reflects the effect of HLA-DRB. Declaration of interest: The authors declare no conflict of interest.

Page of Funding The study was supported by grants from the Regional Health Authorities of Western Norway, Bergen Medical Research Foundation, The Norwegian Research council (project no 0) and the EU FP project Euradrenal (grant number 0). I.B is a PhD student supported by University of Bergen. 0 Acknowledgements We thank Hajirah Muneer Elin Theodorsen and Elisabeth Tombra Halvorsen for excellent organizing of DNA sampling and antibody analyses. We also thank all the patients and their medical doctors who have contribute to the National Norwegian Registry of patients with Addison s disease (ROAS).

Page of References 0 0. Arlt W & Allolio B. Adrenal insufficiency. Lancet 00 -.. Erichsen MM, Lovas K, Skinningsrud B, Wolff AB, Undlien DE, Svartberg J, Fougner KJ, Berg TJ, Bollerslev J, Mella B, Carlson JA, Erlich H & Husebye ES. Clinical, immunological, and genetic features of autoimmune primary adrenal insufficiency: observations from a Norwegian registry. J Clin Endocrinol Metab 00-0.. Husebye ES, Perheentupa J, Rautemaa R & Kampe O. Clinical manifestations and management of patients with autoimmune polyendocrine syndrome type I. J Intern Med 00 -.. Skinningsrud B, Lie BA, Lavant E, Carlson JA, Erlich H, Akselsen HE, Gervin K, Wolff AB, Erichsen MM, Lovas K, Husebye ES & Undlien DE. Multiple loci in the HLA complex are associated with Addison's disease. J Clin Endocrinol Metab 0 E0-0.. Gambelunghe G, Falorni A, Ghaderi M, Laureti S, Tortoioli C, Santeusanio F, Brunetti P & Sanjeevi CB. Microsatellite polymorphism of the MHC class I chainrelated (MIC-A and MIC-B) genes marks the risk for autoimmune Addison's disease. J Clin Endocrinol Metab 0-0.. Park YS, Sanjeevi CB, Robles D, Yu L, Rewers M, Gottlieb PA, Fain P & Eisenbarth GS. Additional association of intra-mhc genes, MICA and DS, with Addison's disease. Tissue Antigens 00 0 -.. Triolo TM, Baschal EE, Armstrong TK, Toews CS, Fain PR, Rewers MJ, Yu L, Miao D, Eisenbarth GS, Gottlieb PA & Barker JM. Homozygosity of the polymorphism MICA. identifies extreme risk of progression to overt adrenal insufficiency among

Page of 0 0 -hydroxylase antibody-positive patients with type diabetes. J Clin Endocrinol Metab 00 -.. Blomhoff A, Lie BA, Myhre AG, Kemp EH, Weetman AP, Akselsen HE, Huseby ES & Undlien DE. Polymorphisms in the cytotoxic T lymphocyte antigen- gene region confer susceptibility to Addison's disease. J Clin Endocrinol Metab 00 -.. Brozzetti A, Marzotti S, Tortoioli C, Bini V, Giordano R, Dotta F, Betterle C, De Bellis A, Arnaldi G, Toscano V, Arvat E, Bellastella A, Mantero F, Falorni A & Italian Addison N. Cytotoxic T lymphocyte antigen- Ala polymorphism is a genetic marker of autoimmune adrenal insufficiency: Italian association study and meta-analysis of European studies. Eur J Endocrinol 00 -. 0. Roycroft M, Fichna M, McDonald D, Owen K, Zurawek M, Gryczynska M, Januszkiewicz-Lewandowska D, Fichna P, Cordell H, Donaldson P, Nowak J & Pearce S. The tryptophan 0 allele of the lymphoid tyrosine phosphatase (PTPN) gene predisposes to autoimmune Addison's disease. Clin Endocrinol (Oxf) 00 0 -.. Skinningsrud B, Husebye ES, Gervin K, Lovas K, Blomhoff A, Wolff AB, Kemp EH, Egeland T & Undlien DE. Mutation screening of PTPN: association of the Tallele with Addison's disease. Eur J Hum Genet 00 -.. Skinningsrud B, Husebye ES, Pearce SH, McDonald DO, Brandal K, Wolff AB, Lovas K, Egeland T & Undlien DE. Polymorphisms in CLECA and CIITA at p are associated with primary adrenal insufficiency. J Clin Endocrinol Metab 00 0-.

Page of 0 0. Jennings CE, Owen CJ, Wilson V & Pearce SH. A haplotype of the CYPB promoter is associated with autoimmune Addison's disease but not with Graves' disease in a UK population. J Mol Endocrinol 00 -.. Lopez ER, Zwermann O, Segni M, Meyer G, Reincke M, Seissler J, Herwig J, Usadel KH & Badenhoop K. A promoter polymorphism of the CYPB gene is associated with Addison's disease, Hashimoto's thyroiditis, Graves' disease and type diabetes mellitus in Germans. Eur J Endocrinol 00 -.. Mitchell AL, Cordell HJ, Soemedi R, Owen K, Skinningsrud B, Wolff AB, Ericksen M, Undlien D, Husebye E & Pearce SH. Programmed death ligand (PD-L) gene variants contribute to autoimmune Addison's disease and Graves' disease susceptibility. J Clin Endocrinol Metab 00 -.. Winqvist O, Karlsson FA & Kampe O. -Hydroxylase, a major autoantigen in idiopathic Addison's disease. Lancet -.. Betterle C, Scarpa R, Garelli S, Morlin L, Lazzarotto F, Presotto F, Coco G, Masiero S, Parolo A, Albergoni MP, Favero R, Barollo S, Salva M, Basso D, Chen S, Rees Smith B, Furmaniak J & Mantero F. Addison's disease: a survey on patients in Padova. Eur J Endocrinol 0 -.. Falorni A, Laureti S, Nikoshkov A, Picchio ML, Hallengren B, Vandewalle CL, Gorus FK, Tortoioli C, Luthman H, Brunetti P & Santeusanio F. -hydroxylase autoantibodies in adult patients with endocrine autoimmune diseases are highly specific for Addison's disease. Belgian Diabetes Registry. Clin Exp Immunol 0 -.. Seissler J, Schott M, Steinbrenner H, Peterson P & Scherbaum WA. Autoantibodies to adrenal cytochrome P0 antigens in isolated Addison's disease and autoimmune polyendocrine syndrome type II. Exp Clin Endocrinol Diabetes 0 0-.

Page of 0 0 0. Speiser PW, Azziz R, Baskin LS, Ghizzoni L, Hensle TW, Merke DP, Meyer- Bahlburg HF, Miller WL, Montori VM, Oberfield SE, Ritzen M, White PC & Endocrine S. Congenital adrenal hyperplasia due to steroid -hydroxylase deficiency: an Endocrine Society clinical practice guideline. J Clin Endocrinol Metab 00-0.. White PC, New MI & Dupont B. Structure of human steroid -hydroxylase genes. Proc Natl Acad Sci U S A -.. Wenzlau JM, Liu Y, Yu L, Moua O, Fowler KT, Rangasamy S, Walters J, Eisenbarth GS, Davidson HW & Hutton JC. A common nonsynonymous single nucleotide polymorphism in the SLC0A gene determines ZnT autoantibody specificity in type diabetes. Diabetes 00 -.. Dechairo BM, Zabaneh D, Collins J, Brand O, Dawson GJ, Green AP, Mackay I, Franklyn JA, Connell JM, Wass JA, Wiersinga WM, Hegedus L, Brix T, Robinson BG, Hunt PJ, Weetman AP, Carey AH & Gough SC. Association of the TSHR gene with Graves' disease: the first disease specific locus. Eur J Hum Genet 00-0.. Brand OJ, Barrett JC, Simmonds MJ, Newby PR, McCabe CJ, Bruce CK, Kysela B, Carr-Smith JD, Brix T, Hunt PJ, Wiersinga WM, Hegedus L, Connell J, Wass JA, Franklyn JA, Weetman AP, Heward JM & Gough SC. Association of the thyroid stimulating hormone receptor gene (TSHR) with Graves' disease. Hum Mol Genet 00 0-.. Nikoshkov A, Falorni A, Lajic S, Laureti S, Wedell A, Lernmark K & Luthman H. A conformation-dependent epitope in Addison's disease and other endocrinological autoimmune diseases maps to a carboxyl-terminal functional domain of human steroid -hydroxylase. J Immunol -.

Page of 0 0. Peterson P, Partanen J, Aavik E, Salmi H, Pelkonen R & Krohn KJ. Steroid - hydroxylase gene polymorphism in Addison's disease patients. Tissue Antigens -.. Wedell A & Luthman H. Steroid -hydroxylase deficiency: two additional mutations in salt-wasting disease and rapid screening of disease-causing mutations. Hum Mol Genet -0.. Wilson RC, Wei JQ, Cheng KC, Mercado AB & New MI. Rapid deoxyribonucleic acid analysis by allele-specific polymerase chain reaction for detection of mutations in the steroid -hydroxylase gene. J Clin Endocrinol Metab 0-0.. Lau IF, Soardi FC, Lemos-Marini SH, Guerra Jr G, Jr., Baptista MT & De Mello MP. H+C insertion in the CYP gene: a novel frameshift mutation in a Brazilian patient with the classical form of -hydroxylase deficiency. J Clin Endocrinol Metab 00-0. 0. Parajes S, Quinterio C, Dominguez F & Loidi L. A simple and robust quantitative PCR assay to determine CYPA gene dose in the diagnosis of -hydroxylase deficiency. Clin Chem 00 -.. Barrett JC, Fry B, Maller J & Daly MJ. Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics 00 -.. Dudbridge F. Likelihood-based association analysis for nuclear families and unrelated subjects with missing genotype data. Hum Hered 00 -.. Stephens M, Smith NJ & Donnelly P. A new statistical method for haplotype reconstruction from population data. Am J Hum Genet 00 -.. Wedell A. Congenital adrenal hyperplasia. Clin Biochem 0 0-0.. Baker PR, Baschal EE, Fain PR, Triolo TM, Nanduri P, Siebert JC, Armstrong TK, Babu SR, Rewers MJ, Gottlieb PA, Barker JM & Eisenbarth GS. Haplotype analysis

Page 0 of 0 0 discriminates genetic risk for DR-associated endocrine autoimmunity and helps define extreme risk for Addison's disease. J Clin Endocrinol Metab 00 E- 0.. Yu L, Brewer KW, Gates S, Wu A, Wang T, Babu SR, Gottlieb PA, Freed BM, Noble J, Erlich HA, Rewers MJ & Eisenbarth GS. DRB*0 and DQ alleles: expression of -hydroxylase autoantibodies and risk of progression to Addison's disease. J Clin Endocrinol Metab -.. Blasko B, Banlaki Z, Gyapay G, Pozsonyi E, Sasvari-Szekely M, Rajczy K, Fust G & Szilagyi A. Linkage analysis of the CA/CB copy number variation and polymorphisms of the adjacent steroid -hydroxylase gene in a healthy population. Mol Immunol 00 -.. Barker JM, Ide A, Hostetler C, Yu L, Miao D, Fain PR, Eisenbarth GS & Gottlieb PA. Endocrine and immunogenetic testing in individuals with type diabetes and - hydroxylase autoantibodies: Addison's disease in a high-risk population. J Clin Endocrinol Metab 00 0 -.. Durinovic-Bello I, Wu RP, Gersuk VH, Sanda S, Shilling HG & Nepom GT. Insulin gene VNTR genotype associates with frequency and phenotype of the autoimmune response to proinsulin. Genes Immun 00 -. 0. Pugliese A, Zeller M, Fernandez A, Jr., Zalcberg LJ, Bartlett RJ, Ricordi C, Pietropaolo M, Eisenbarth GS, Bennett ST & Patel DD. The insulin gene is transcribed in the human thymus and transcription levels correlated with allelic variation at the INS VNTR-IDDM susceptibility locus for type diabetes. Nat Genet -. 0

Page of

Page of Legend to figures: Figure a) r linkage plot of common variants of CYPA. b) D linkage plot of common variants of CYPA.

Page of Table Differences in genotype frequencies (%) of the most common variants in CYPA detected in patients with autoimmune Addison s disease (AAD) and healthy controls (HC). P values were calculated by Chi square test for x tables and Fisher s exact test on x tables. CYPA SNP Genotype P - value OR % CI P - value rs +/+ +/- -/- +/+ vs +/- AAD () 0 (.) (.) 0.00... 0.000 HC (.) (.) 0 (.0) rs AA AC CC AA vs AC AAD (.) () (.) 0.000.. -. 0.000 HC 0 (.) (.) (.0) rs GG GA AA GG vs GA AAD 0 (.) 0 (.) (.) 0.00... 0.00 HC (.0) (.) (.) rs CC CT TT CC vs CT AAD (.0) (.) 0 (0.) 0.0000.0... e- HC (.) (0.) (.) rs GG GC CC G vs C AAD 0 (.) 0 (.) 0 (0). e- 0. 0.0 0.0. e-0 HC (.0) (.) (.) rs GG GA AA GG vs GA AAD (.) (.) (.0) 0.00... 0.000 HC (.) (.0) (.)

Fisher s exact test Page of

Page of Table Global regression analysis of the HLA and CYPA loci for determination of independent associations in autoimmune Addison s disease calculated by UNPHASED. Test locus Analysis for conditional locus (P - values) A B MICA DRB DQB rs rs rs rs rs rs HLA-A 0.0 0.00 0. 0.0. e- 0.00. e- 0.00. e-.0 e- HLA-B. e-. e- 0.0. e-. e-.0 e-0. e-. e-.0 e-0. e- MICA.0 e- 0.000 0.00. e-. e-. e-. e-.0. e-. e-. e- HLA-DRB. e-. e-. e-. e-. e-. e-. e-0. e-. e-. e- HLA-DQB. e-0. e-. e- 0.0. e-. e-. e-. e-.0 e-0. e- rs. e-. e-. e- 0.0 0.0 0.00 0.0. e- 0.0 0. rs 0.0 0. 0. 0. 0.0 0.000 0.000 0.0 0.00 0.000 rs 0.000. e- 0.000 0. 0.0 0. 0.00. e- 0. 0. rs 0. 0.0 0. 0. 0.. e- 0.. e- 0.000. e- rs.0 e-. e-. e- 0.0 0.000.0 e-. e-. e-. e-. e- rs.0 e-.0 e-. e- 0. 0.0 0. 0.00 0.. e- 0.0

Page of Table Haplotype reconstruction of protective and high risk alleles of HLA-DRB and rs. P values were calculated by Fisher s exact test. HLA-DRB rs AAD HC Odds Ratio % Confidence Interval P-Value *0 C *0 G *00 C *00 G *0 C 0 *0 G *00 G 0 *00 C 0 0.0 0.0. 0. 0.0 0...0 0. 0.00.0 0.0 0. 0.0. 0.

Page of xmm ( x DPI)