Extremely high frequency of autoimmune-predisposing alleles in medieval specimens

Extremely high frequency of autoimmune-predisposing alleles in medieval specimens first-author †E-mail:witas@usk4.umed.lodz.pl †E-mail:witas@usk4.umed.lodz.pl witas@usk4.umed.lodz.pl Department of Molecular Biology, Medical University of Lodz, Sporna 36/50, PL-91738 Lodz, Poland Department of Molecular Biology, Medical University of Lodz, Sporna 36/50, PL-91738 Lodz, Poland The precise etiology and reasons for the increase in incidence of autoimmune disorders still remain unclear, and although both genetic and environmental factors have been proven to shape individual predisposition, it is not known which of the factors, if not both, is responsible for the boom observed during the last decades. In order to establish whether a higher frequency of autoimmune-predisposing alleles may explain this increase we took advantage of ancient DNA methodology to establish the genetic predisposition, conferred by cytotoxic T lymphocyte associated antigen-4 ( CTLA4 ) +49A/G and human leukocyte antigens ( HLA ) DQB1 57 , in population inhabiting Poland in the Middle Ages. After successful typing of 42 individuals from a 12th~14th’s century archeological burial site, we found that frequencies of the predisposing alleles in the medieval population were higher than they are at present, suggesting thus that the recently observed incidence increase results most probably from factors of other than genetic nature. Ancient DNA (aDNA) Autoimmunity Cytotoxic T lymphocyte associated antigen-4 ( CTLA4 ) gene HLA DQB1 Type 1 diabetes (T1D) INTRODUCTION 1 It is estimated that up to 3%~5% of the world’s population is affected by at least one disorder that results from autoimmunity (Marrack et al., 2001 2001 ), and the incidence of at least some of the autoimmune diseases has been growing steadily for the last decades (Onkamo et al., 1999 1999 ; Ivarsson et al., 2000 2000 ; Haynes et al., 2004 2004 ). It is estimated that autoimmune-related disorders are the third (after cancer and atherosclerosis) cause of mortality and morbidity in developed countries (Chatenoud, 2006 2006 ). A number of loci have been thus far proven to exert an influence on predisposition towards autoaggressive responses (Aune et al., 2004 2004 ), among which the most influential are HLA (human leukocyte antigens) genes, CTLA4 (cytotoxic T lymphocyte associated antigen-4) and PTPN22 (protein tyrosine phosphatase non-receptor type 22) (Pearce and Merriman, 2006 2006 ). However, studies on monozygotic twins and affected families indicate that autoimmunity is triggered by other than genetic factors. There have been many possible environmental factors found to influence autoimmune response, such as certain drugs, pathogens, pollutants, diet, newly introduced chemicals or stress (Gerstein, 1994 1994 ; Longnecker and Daniels, 2001 2001 ; Roivainen, 2006 2006 ; Dahlquist, 2006 2006 ). Nevertheless, none of the genetic factors, nor an environmental one, was proven to be responsible for the observed boom of the diseases incidence. It is feasible, that since the living conditions improved and more affected individuals reach procreation age due to advances in medical treatment, the ‘autoimmune’ alleles are preserved in the gene pool, making populations more prone to this kind of disorders. To find out whether this might be at least partly the case, we decided to study the balance in the case of two predisposing alleles’ frequencies in past populations, by means of ancient DNA (aDNA) analysis, in order to compare with their contemporary presence. We chose CTLA4 +49A/G and HLA DQB 57 polymorphisms—the only confirmed for the Polish population sites involved in the onset of certain autoimmune mediated diseases and found to influence predisposition towards type 1 diabetes (T1D) and coeliac disease among others. We decided to analyse a population inhabiting Poland in the Middle Ages. Taking advantage of such a relatively old group ensures that the frequency of genotypes was not influenced by any major human movements or chemicals introduced to the environment during the later periods of industrial changes. MATERIALS AND METHODS 1 Materials 2 Analyzed medieval DNA was isolated from teeth of over 100 non-related individuals, collected on an archeological burial site in Stary Brzesc Kujawski (SBK-4) central Poland, dated back to 12th~14th century. The reference groups (Witas et al., 1998 1998 ; Witas and Młynarski, 1999 1999 ) comprised individuals inhabiting contemporarily the same region of Poland. Methods 2 Each tooth was mechanically and chemically cleaned, powdered and after pre-treatment with EDTA (ethylene diamine tetraacetic acid), proteinase K and PTB (N-phenacylthiazolium bromide), subjected to DNA isolation in MagNA Pure® Compact nucleic acid extraction system (Roche). Sequences of interest were subsequently amplified from isolated DNA fragments and identified using RFLP (restriction fragment length polimorphism) analysis ( CTLA4 +49A/G) or ASA-PCR (allele-specific amplification, polymerase chain reaction) ( HLA DQB1 57 ). Alternatively, amplification products were sequenced. Primers used in PCR were designed to amplify as short sequences as possible (116 and 100 bp respectively), which are more likely to be extracted from aDNA. As necessary with all aDNA isolation and analytical procedures, special attention was given not to contaminate the samples while collecting, and to provide controlled purity of environment for analysis. All work was done in a laminar air flow workstation, in rooms specially dedicated to aDNA work, where no contemporary DNA had ever been analyzed. All personnel involved wore disposable clothing and all disposables and reagents used throughout isolation and amplification procedures had been carefully selected from those found to be most reliable (Schmidt et al., 1995 1995 ). Only those results were considered that could be confirmed in at least two independent DNA isolation procedures from different samples of the same specimen. RESULTS AND DISCUSSION 1 Our preliminary results for CTLA4 +49A/G polymorphism (Table 1 1 ) indicate that the frequency of the predisposing G allele was higher in medieval Poland than it is contemporarily observed. Similarly, the frequency value of the protecting Asp codon at HLA DQB1 57 (Table 2 2 ) is closer to contemporary diabetic patients than to contemporary control group. Testing for trend in frequency distribution of predisposing alleles clearly confirmed the direction of changes from contemporary control through medieval close in value to contemporary diabetics: P < 0.001 for linear frequency trend of GG vs remaining genotypes, allele G vs allele A frequency; XX vs Asp and Asp vs non-Asp alleles. This strongly suggests that the increase in autoimmune diseases incidence observed within the last decades, especially T1D, does not result from changes in predisposing alleles’ frequencies of the studied loci, at least in Polish population. Moreover, the studied HLA DQB1 57 and CTLA4 +49A/G alleles’ frequencies seem to confer even higher risk for the medieval population. The results remain in accordance with other authors’ findings. Hermann et al.( 2003 2003 ) found the frequency of HLA DRB1 - DQA1 - DQB1 T1D protecting haplotypes higher in contemporary Finland when compared to mid-20th century population, despite 2.5-fold increase in the disease incidence, which may, similarly to our results, point at environment changes as a reason for the recent boom of autoimmune diseases increase. Further work is being currently done to increase studied group quantity, to raise the number of studied sites, and to determine not only the aforementioned but other autoimmune-predisposing loci frequencies in past populations in order to add to the discussion on potential factors involved in the increase of autoimmune disorders incidence. According to Witas et al.( 1998 1998 ) Group AA (%) AG (%) GG (%) A (%) G (%) Contemporary control * * ( n =223) 40.8 51.6 7.6 66.6 33.4 Medieval population ( n =42) 23.8 52.4 23.8 50.0 50.0 Contemporary diabetics * * ( n =207) 26.5 49.8 23.7 51.4 48.6 According to Witas and Młynarski ( 1999 1999 ) X: Amino acid other than Asp coded at HLA DQB1 57 Group Asp/Asp (%) Asp/X (%) X/X (%) Asp (%) Contemporary control * * ( n =78) 27.0 47.0 26.0 50.64 Medieval population ( n =64) 14.1 31.2 54.7 29.69 Contemporary diabetics * * ( n =75) 0.0 25.0 75.0 12.67 References TM Aune , K Maas , J Parker , JH Moore , NJ Olsen . Profiles of gene expression in human autoimmune disease . Cell Biochem Biophys , 2004 . 40 ( 2 ): 81 - 96 . 15054216 DOI: 10.1385/CBB:40:2:081 http://doi.org/10.1385/CBB:40:2:081 . Profiles of gene expression in human autoimmune disease L Chatenoud . Immune therapies of autoimmune diseases: are we approaching a real cure? . Curr Opin Immunol , 2006 . 18 ( 6 ): 710 - 717 . 17011768 DOI: 10.1016/j.coi.2006.09.004 http://doi.org/10.1016/j.coi.2006.09.004 . Immune therapies of autoimmune diseases: are we approaching a real cure? G Dahlquist . Can we slow the rising incidence of childhood-onset autoimmune diabetes? The overload hypothesis . Diabetologia , 2006 . 49 ( 1 ): 20 - 24 . 16362279 DOI: 10.1007/s00125-005-0076-4 http://doi.org/10.1007/s00125-005-0076-4 . Can we slow the rising incidence of childhood-onset autoimmune diabetes? The overload hypothesis HC Gerstein . Cow’s milk exposure and type I diabetes mellitus . Diabetes Care , 1994 . 17 ( 1 ): 13 - 19 . 8112184 DOI: 10.2337/diacare.17.1.13 http://doi.org/10.2337/diacare.17.1.13 . Cow’s milk exposure and type I diabetes mellitus A Haynes , C Bower , MK Bulsara , TW Jones , EA Davis . Continued increase in the incidence of childhood type 1 diabetes in a population-based Australian sample (1985~2002) . Diabetologia , 2004 . 47 ( 5 ): 866 - 870 . 15095039 DOI: 10.1007/s00125-004-1385-8 http://doi.org/10.1007/s00125-004-1385-8 . Continued increase in the incidence of childhood type 1 diabetes in a population-based Australian sample (1985~2002) R Hermann , M Knip , R Veijola , O Simell , AP Laine , HK Akerblom , PH Groop , C Forsblom , K Pettersson-Fernholm , J Ilonen . FinnDiane Study Group Temporal changes in the frequencies of HLA genotypes in patients with type 1 diabetes-indication of an increased environmental pressure? . Diabetologia , 2003 . 46 ( 3 ): 420 - 425 . 12687342 . Temporal changes in the frequencies of HLA genotypes in patients with type 1 diabetes-indication of an increased environmental pressure? A Ivarsson , LA Persson , L Nyström , H Ascher , B Cavell , L Danielsson , A Dannaeus , T Lindberg , B Lindquist , L Stenhammer , 等 . Epidemic of coeliac disease in Swedish children . Acta Paediatr , 2000 . 89 ( 2 ): 165 - 171 . 10709885 DOI: 10.1080/080352500750028771 http://doi.org/10.1080/080352500750028771 . Epidemic of coeliac disease in Swedish children MP Longnecker , JL Daniels . Environmental contaminants as etiological factors for diabetes . Environ Health Persp , 2001 . 109 871 - 876 . DOI: 10.2307/3454649 http://doi.org/10.2307/3454649 . Environmental contaminants as etiological factors for diabetes P Marrack , J Kappler , B Kotzin . Autoimmune disease: why and where in occurs . Nat Med , 2001 . 7 ( 8 ): 899 - 905 . 11479621 DOI: 10.1038/90935 http://doi.org/10.1038/90935 . Autoimmune disease: why and where in occurs P Onkamo , S Vaananen , M Karvonen , J Tuomilehto . Worldwide increase in incidence of type 1 diabetes—the analysis of the data on published incidence trends . Diabetologia , 1999 . 42 ( 12 ): 1395 - 1403 . 10651256 DOI: 10.1007/s001250051309 http://doi.org/10.1007/s001250051309 . Worldwide increase in incidence of type 1 diabetes—the analysis of the data on published incidence trends SHS Pearce , TR Merriman . Genetic progress towards the molecular basis of autoimmunity . Trends Mol Med , 2006 . 12 ( 2 ): 90 - 98 . 16412690 DOI: 10.1016/j.molmed.2005.12.005 http://doi.org/10.1016/j.molmed.2005.12.005 . Genetic progress towards the molecular basis of autoimmunity M Roivainen . Enteroviruses: new findings on the role of enteroviruses in type 1 diabetes . Int J Biochem Cell Biol , 2006 . 38 ( 5-6 ): 721 - 725 . 16226050 DOI: 10.1016/j.biocel.2005.08.019 http://doi.org/10.1016/j.biocel.2005.08.019 . Enteroviruses: new findings on the role of enteroviruses in type 1 diabetes T Schmidt , S Hummel , B Herrmann . Evidence of contamination in PCR laboratory disposables . Naturwis-senschaften , 1995 . 82 ( 9 ): 423 - 431 . DOI: 10.1007/s001140050212 http://doi.org/10.1007/s001140050212 . Evidence of contamination in PCR laboratory disposables HW Witas , W Młynarski . Type I diabetes and asthma: alteration of HLA-DQβ 57 codon is associated with pathological and physiological status . Biomed Lett , 1999 . 59 15 - 26 . . Type I diabetes and asthma: alteration of HLA-DQβ57 codon is associated with pathological and physiological status HW Witas , W Młynarski , R Sychowski . Polymorphism located in promoter of CTLA-4 gene is not associated with type 1 diabetes in contrast to A-G transition within exon 1 . Biomed Lett , 1998 . 58 163 - 168 . . Polymorphism located in promoter of CTLA-4 gene is not associated with type 1 diabetes in contrast to A-G transition within exon 1 We are deeply indebted to Prof. E. Żądzińska for providing us with dated archaeological specimens and Prof. C. Watała for careful statistical analysis. 10.1631/jzus.2007.B0512 Project (No. 3 P05E 055 25) supported by the Ministry of Education and Science of Poland 2007-03-03 2007-05-09 2022-01-05 Copyright © 2007, Journal of Zhejiang University Science 2007 Journal of Zhejiang University-SCIENCE B (Biomedicine and Biotechnology) 07 08