1
|
Hayeck TJ, Li Y, Mosbruger TL, Bradfield JP, Gleason AG, Damianos G, Shaw GTW, Duke JL, Conlin LK, Turner TN, Fernández-Viña MA, Sarmady M, Monos DS. The Impact of Patterns in Linkage Disequilibrium and Sequencing Quality on the Imprint of Balancing Selection. Genome Biol Evol 2024; 16:evae009. [PMID: 38302106 PMCID: PMC10853003 DOI: 10.1093/gbe/evae009] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/23/2023] [Revised: 01/08/2024] [Accepted: 01/12/2024] [Indexed: 02/03/2024] Open
Abstract
Regions under balancing selection are characterized by dense polymorphisms and multiple persistent haplotypes, along with other sequence complexities. Successful identification of these patterns depends on both the statistical approach and the quality of sequencing. To address this challenge, at first, a new statistical method called LD-ABF was developed, employing efficient Bayesian techniques to effectively test for balancing selection. LD-ABF demonstrated the most robust detection of selection in a variety of simulation scenarios, compared against a range of existing tests/tools (Tajima's D, HKA, Dng, BetaScan, and BalLerMix). Furthermore, the impact of the quality of sequencing on detection of balancing selection was explored, as well, using: (i) SNP genotyping and exome data, (ii) targeted high-resolution HLA genotyping (IHIW), and (iii) whole-genome long-read sequencing data (Pangenome). In the analysis of SNP genotyping and exome data, we identified known targets and 38 new selection signatures in genes not previously linked to balancing selection. To further investigate the impact of sequencing quality on detection of balancing selection, a detailed investigation of the MHC was performed with high-resolution HLA typing data. Higher quality sequencing revealed the HLA-DQ genes consistently demonstrated strong selection signatures otherwise not observed from the sparser SNP array and exome data. The HLA-DQ selection signature was also replicated in the Pangenome samples using considerably less samples but, with high-quality long-read sequence data. The improved statistical method, coupled with higher quality sequencing, leads to more consistent identification of selection and enhanced localization of variants under selection, particularly in complex regions.
Collapse
Affiliation(s)
- Tristan J Hayeck
- Division of Genomic Diagnostics, Department of Pathology and Laboratory Medicine, Children's Hospital of Philadelphia, Philadelphia, PA, USA
- Department of Pathology and Laboratory Medicine, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
| | - Yang Li
- Division of Genomic Diagnostics, Department of Pathology and Laboratory Medicine, Children's Hospital of Philadelphia, Philadelphia, PA, USA
- Department of Pathology and Laboratory Medicine, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
| | - Timothy L Mosbruger
- Division of Genomic Diagnostics, Department of Pathology and Laboratory Medicine, Children's Hospital of Philadelphia, Philadelphia, PA, USA
| | | | - Adam G Gleason
- Division of Genomic Diagnostics, Department of Pathology and Laboratory Medicine, Children's Hospital of Philadelphia, Philadelphia, PA, USA
| | - George Damianos
- Division of Genomic Diagnostics, Department of Pathology and Laboratory Medicine, Children's Hospital of Philadelphia, Philadelphia, PA, USA
| | - Grace Tzun-Wen Shaw
- Division of Genomic Diagnostics, Department of Pathology and Laboratory Medicine, Children's Hospital of Philadelphia, Philadelphia, PA, USA
| | - Jamie L Duke
- Division of Genomic Diagnostics, Department of Pathology and Laboratory Medicine, Children's Hospital of Philadelphia, Philadelphia, PA, USA
| | - Laura K Conlin
- Division of Genomic Diagnostics, Department of Pathology and Laboratory Medicine, Children's Hospital of Philadelphia, Philadelphia, PA, USA
- Department of Pathology and Laboratory Medicine, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
| | - Tychele N Turner
- Department of Genetics, Washington University School of Medicine, St. Louis, MO 63110, USA
| | - Marcelo A Fernández-Viña
- Department of Pathology, Stanford University School of Medicine, Palo Alto, CA, USA
- Histocompatibility and Immunogenetics Laboratory, Stanford Blood Center, Palo Alto, CA, USA
| | - Mahdi Sarmady
- Division of Genomic Diagnostics, Department of Pathology and Laboratory Medicine, Children's Hospital of Philadelphia, Philadelphia, PA, USA
- Department of Pathology and Laboratory Medicine, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
| | - Dimitri S Monos
- Division of Genomic Diagnostics, Department of Pathology and Laboratory Medicine, Children's Hospital of Philadelphia, Philadelphia, PA, USA
- Department of Pathology and Laboratory Medicine, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
| |
Collapse
|
2
|
Douillard V, Dos Santos Brito Silva N, Bourguiba-Hachemi S, Naslavsky MS, Scliar MO, Duarte YAO, Zatz M, Passos-Bueno MR, Limou S, Gourraud PA, Launay É, Castelli EC, Vince N. Optimal population-specific HLA imputation with dimension reduction. HLA 2024; 103:e15282. [PMID: 37950640 DOI: 10.1111/tan.15282] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/19/2023] [Revised: 08/29/2023] [Accepted: 10/14/2023] [Indexed: 11/13/2023]
Abstract
Human genomics has quickly evolved, powering genome-wide association studies (GWASs). SNP-based GWASs cannot capture the intense polymorphism of HLA genes, highly associated with disease susceptibility. There are methods to statistically impute HLA genotypes from SNP-genotypes data, but lack of diversity in reference panels hinders their performance. We evaluated the accuracy of the 1000 Genomes data as a reference panel for imputing HLA from admixed individuals of African and European ancestries, focusing on (a) the full dataset, (b) 10 replications from 6 populations, and (c) 19 conditions for the custom reference panels. The full dataset outperformed smaller models, with a good F1-score of 0.66 for HLA-B. However, custom models outperformed the multiethnic or population models of similar size (F1-scores up to 0.53, against up to 0.42). We demonstrated the importance of using genetically specific models for imputing populations, which are currently underrepresented in public datasets, opening the door to HLA imputation for every genetic population.
Collapse
Affiliation(s)
- Venceslas Douillard
- Nantes Université, INSERM, Ecole Centrale Nantes, Center for Research in Transplantation and Translational Immunology, Nantes, France
| | - Nayane Dos Santos Brito Silva
- Nantes Université, INSERM, Ecole Centrale Nantes, Center for Research in Transplantation and Translational Immunology, Nantes, France
- São Paulo State University, Molecular Genetics and Bioinformatics Laboratory, School of Medicine, Botucatu, Brazil
| | - Sonia Bourguiba-Hachemi
- Nantes Université, INSERM, Ecole Centrale Nantes, Center for Research in Transplantation and Translational Immunology, Nantes, France
| | - Michel S Naslavsky
- Human Genome and Stem Cell Research Center, University of São Paulo, São Paulo, Brazil
- Department of Genetics and Evolutionary Biology, Biosciences Institute, University of São Paulo, São Paulo, Brazil
- Hospital Israelita Albert Einstein, São Paulo, Brazil
| | - Marilia O Scliar
- Human Genome and Stem Cell Research Center, University of São Paulo, São Paulo, Brazil
| | - Yeda A O Duarte
- Medical-Surgical Nursing Department, School of Nursing, University of São Paulo, São Paulo, Brazil
- Epidemiology Department, Public Health School, University of São Paulo, São Paulo, Brazil
| | - Mayana Zatz
- Human Genome and Stem Cell Research Center, University of São Paulo, São Paulo, Brazil
- Department of Genetics and Evolutionary Biology, Biosciences Institute, University of São Paulo, São Paulo, Brazil
| | - Maria Rita Passos-Bueno
- Human Genome and Stem Cell Research Center, University of São Paulo, São Paulo, Brazil
- Department of Genetics and Evolutionary Biology, Biosciences Institute, University of São Paulo, São Paulo, Brazil
| | - Sophie Limou
- Nantes Université, INSERM, Ecole Centrale Nantes, Center for Research in Transplantation and Translational Immunology, Nantes, France
| | - Pierre-Antoine Gourraud
- Nantes Université, INSERM, Ecole Centrale Nantes, Center for Research in Transplantation and Translational Immunology, Nantes, France
| | - Élise Launay
- Nantes Université, INSERM, Ecole Centrale Nantes, Center for Research in Transplantation and Translational Immunology, Nantes, France
- Department of Pediatrics and Pediatric Emergency, Hôpital Femme Enfant Adolescent, CHU de Nantes, Nantes, France
| | - Erick C Castelli
- São Paulo State University, Molecular Genetics and Bioinformatics Laboratory, School of Medicine, Botucatu, Brazil
| | - Nicolas Vince
- Nantes Université, INSERM, Ecole Centrale Nantes, Center for Research in Transplantation and Translational Immunology, Nantes, France
| |
Collapse
|
3
|
Kalomoiri M, Prakash CR, Lagström S, Hauschulz K, Ewing E, Shchetynsky K, Kular L, Needhamsen M, Jagodic M. Simultaneous detection of DNA variation and methylation at HLA class II locus and immune gene promoters using targeted SureSelect Methyl-Sequencing. Front Immunol 2023; 14:1251772. [PMID: 37691926 PMCID: PMC10484099 DOI: 10.3389/fimmu.2023.1251772] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/02/2023] [Accepted: 08/02/2023] [Indexed: 09/12/2023] Open
Abstract
The Human Leukocyte Antigen (HLA) locus associates with a variety of complex diseases, particularly autoimmune and inflammatory conditions. The HLA-DR15 haplotype, for example, confers the major risk for developing Multiple Sclerosis in Caucasians, pinpointing an important role in the etiology of this chronic inflammatory disease of the central nervous system. In addition to the protein-coding variants that shape the functional HLA-antigen-T cell interaction, recent studies suggest that the levels of HLA molecule expression, that are epigenetically controlled, also play a role in disease development. However, deciphering the exact molecular mechanisms of the HLA association has been hampered by the tremendous genetic complexity of the locus and a lack of robust approaches to investigate it. Here, we developed a method to specifically enrich the genomic DNA from the HLA class II locus (chr6:32,426,802-34,167,129) and proximal promoters of 2,157 immune-relevant genes, utilizing the Agilent RNA-based SureSelect Methyl-Seq Capture related method, followed by sequencing to detect genetic and epigenetic variation. We demonstrated successful simultaneous detection of the genetic variation and quantification of DNA methylation levels in HLA locus. Moreover, by the detection of differentially methylated positions in promoters of immune-related genes, we identified relevant pathways following stimulation of cells. Taken together, we present a method that can be utilized to study the interplay between genetic variance and epigenetic regulation in the HLA class II region, potentially, in a wide disease context.
Collapse
Affiliation(s)
- Maria Kalomoiri
- Department of Clinical Neuroscience, Center for Molecular Medicine, Karolinska Institutet, Stockholm, Sweden
| | - Chandana Rao Prakash
- Department of Clinical Neuroscience, Center for Molecular Medicine, Karolinska Institutet, Stockholm, Sweden
| | - Sonja Lagström
- Diagnostics and Genomics Group, Agilent Technologies Sweden AB, Sundbyberg, Sweden
| | - Kai Hauschulz
- Diagnostics and Genomics Group, Agilent Technologies Deutschland GmbH, Waldbronn, Germany
| | - Ewoud Ewing
- Department of Clinical Neuroscience, Center for Molecular Medicine, Karolinska Institutet, Stockholm, Sweden
| | - Klementy Shchetynsky
- Department of Clinical Neuroscience, Center for Molecular Medicine, Karolinska Institutet, Stockholm, Sweden
| | - Lara Kular
- Department of Clinical Neuroscience, Center for Molecular Medicine, Karolinska Institutet, Stockholm, Sweden
| | - Maria Needhamsen
- Department of Clinical Neuroscience, Center for Molecular Medicine, Karolinska Institutet, Stockholm, Sweden
| | - Maja Jagodic
- Department of Clinical Neuroscience, Center for Molecular Medicine, Karolinska Institutet, Stockholm, Sweden
| |
Collapse
|
4
|
Khan T, Ledoux IM, Aziz F, Al Ali F, Chin-Smith E, Ata M, Karim MY, Marr N. Associations between HLA class II alleles and IgE sensitization to allergens in the Qatar Biobank cohort. THE JOURNAL OF ALLERGY AND CLINICAL IMMUNOLOGY. GLOBAL 2023; 2:100117. [PMID: 37779520 PMCID: PMC10509938 DOI: 10.1016/j.jacig.2023.100117] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 07/06/2022] [Revised: 01/16/2023] [Accepted: 02/05/2023] [Indexed: 10/03/2023]
Abstract
Background Allergic disorders are the consequence of IgE sensitization to allergens. Population studies have shown that certain human leukocyte antigen (HLA) alleles are associated with increased or decreased risk of developing allergy. Objective We aimed to characterize the relationship between HLA class II allelic diversity and IgE sensitization in an understudied Arab population. Methods We explored associations between IgE sensitization to 7 allergen mixes and mesquite (comprising 41 food or aeroallergens) and 45 common classical HLA class II alleles in a well-defined cohort of 797 individuals representing the general adult population of Qatari nationals and long-term residents. To do so, we performed HLA calling from whole genome sequencing data at 2-field resolution using 2 independent algorithms. We then applied 3 different regression models to assess either each allergen mix independently, in the context of IgE sensitization to other allergens tested, or polysensitization. Results More than half (n = 447) of the study participants showed IgE sensitization to at least 1 allergen, most of them (n = 400) to aeroallergens (Phadiatop). We identified statistically significant negative and positive associations with 24 HLA class II alleles. These have been reported to confer risk or protection from variety of diseases; however, only a few have previously been associated with allergy in other populations. Conclusions Our study reveals several new risk and protective genetic markers for allergen-specific IgE sensitization. This is a first and essential step toward a better understanding of the origins of allergic diseases in this understudied population.
Collapse
Affiliation(s)
- Taushif Khan
- Department of Human Immunology, Sidra Medicine, Doha, Qatar
- Department of Computational Science, The Jackson Laboratory, Farmington, Conn
| | | | - Ferdousey Aziz
- the Department of Pathology, Sidra Medicine, Doha, Qatar
| | - Fatima Al Ali
- Department of Human Immunology, Sidra Medicine, Doha, Qatar
| | | | - Manar Ata
- Department of Human Immunology, Sidra Medicine, Doha, Qatar
| | - Mohammed Yousuf Karim
- the Department of Pathology, Sidra Medicine, Doha, Qatar
- the College of Medicine, Qatar University, Doha, Qatar
| | - Nico Marr
- Department of Human Immunology, Sidra Medicine, Doha, Qatar
- the College of Health and Life Sciences, Hamad Bin Khalifa University, Doha, Qatar
- Institute of Translational Immunology, Brandenburg Medical School, Brandenburg an der Havel, Germany
| |
Collapse
|
5
|
Khan T, Rahman M, Ahmed I, Al Ali F, Jithesh PV, Marr N. Human leukocyte antigen class II gene diversity tunes antibody repertoires to common pathogens. Front Immunol 2022; 13:856497. [PMID: 36003377 PMCID: PMC9393332 DOI: 10.3389/fimmu.2022.856497] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2022] [Accepted: 07/14/2022] [Indexed: 11/13/2022] Open
Abstract
Allelic diversity of human leukocyte antigen (HLA) class II genes may help maintain humoral immunity against infectious diseases. In this study, we investigated germline genetic variation in classical HLA class II genes and employed a systematic, unbiased approach to explore the relative contribution of this genetic variation in the antibody repertoire to various common pathogens. We leveraged a well-defined cohort of 800 adults representing the general Arab population in which genetic material is shared because of the high frequency of consanguineous unions. By applying a high-throughput method for large-scale antibody profiling to this well-defined cohort, we were able to dissect the overall effect of zygosity for classical HLA class II genes, as well as the effects associated with specific HLA class II alleles, haplotypes and genotypes, on the antimicrobial antibody repertoire breadth and antibody specificity with unprecedented resolution. Our population genetic studies revealed that zygosity of the classical HLA class II genes is a strong predictor of antibody responses to common human pathogens, suggesting that classical HLA class II gene heterozygosity confers a selective advantage. Moreover, we demonstrated that multiple HLA class II alleles can have additive effects on the antibody repertoire to common pathogens. We also identified associations of HLA-DRB1 genotypes with specific antigens. Our findings suggest that HLA class II gene polymorphisms confer specific humoral immunity against common pathogens, which may have contributed to the genetic diversity of HLA class II loci during hominine evolution.
Collapse
Affiliation(s)
| | | | | | | | - Puthen Veettil Jithesh
- Research Branch, Sidra Medicine, Doha, Qatar
- College of Health and Life Sciences, Hamad Bin Khalifa University, Doha, Qatar
| | - Nico Marr
- Research Branch, Sidra Medicine, Doha, Qatar
- College of Health and Life Sciences, Hamad Bin Khalifa University, Doha, Qatar
- *Correspondence: Nico Marr,
| |
Collapse
|
6
|
Ås J, Bertulyte I, Eriksson N, Magnusson PK, Wadelius M, Hallberg P. HLA-variants associated with azathioprine-induced pancreatitis in patients with Crohn's disease. Clin Transl Sci 2022; 15:1249-1256. [PMID: 35120281 PMCID: PMC9099136 DOI: 10.1111/cts.13244] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/14/2021] [Revised: 01/14/2022] [Accepted: 01/20/2022] [Indexed: 11/30/2022] Open
Abstract
The immunosuppressant drug azathioprine is associated with a 4% risk of acute pancreatitis in patients with inflammatory bowel disease (IBD). Studies have demonstrated an increased risk in carriers of HLA‐DQA1*02:01 and HLA‐DRB1*07:01. We investigated whether these human leukocyte antigen (HLA) types were associated with azathioprine‐induced pancreatitis also in Swedish patients with IBD, and whether the type of disease affected the association. Nineteen individuals with IBD who developed acute pancreatitis after initiation of azathioprine were genotyped and compared with a population control cohort (n = 4891) and a control group matched for disease (n = 81). HLA‐DQA1*02:01 and HLA‐DRB1*07:01 were in full linkage disequilibrium, and were significantly associated with acute pancreatitis both when cases were compared with population controls (OR 3.97 [95% CI 1.57–9.97], p = 0.0035) and matched controls (OR 3.55 [95% CI 1.23–10.98], p = 0.0275). In a disease‐specific analysis, the correlation was positive in patients with Crohn's disease versus matched controls (OR 9.27 [95% CI 1.86–46.19], p = 0.0066), but not in those with ulcerative colitis versus matched controls (OR 0.69 [95% CI 0.07–6.74], p = 0.749). In patients with Crohn's disease, we estimated the conditional risk of carriers of HLA‐DQA1*02:01‐HLA‐DRB1*07:01 to 7.3%, and the conditional risk of a non‐carrier to 2.2%. We conclude that HLA‐DQA1*02:01‐HLA‐DRB1*07:01 is a marker for increased risk of acute pancreatitis in individuals of Swedish genetic origin, treated with azathioprine for Crohn's disease.
Collapse
Affiliation(s)
- Joel Ås
- Department of Medical Sciences, Uppsala University, Uppsala, Sweden
| | - Ilma Bertulyte
- Department of Medical Sciences, Uppsala University, Uppsala, Sweden
| | | | - Patrik Ke Magnusson
- Swedish Twin Registry, Department of Medical Epidemiology and Biostatistics, Karolinska Institutet, Stockholm, Sweden
| | - Mia Wadelius
- Department of Medical Sciences, Uppsala University, Uppsala, Sweden
| | - Pär Hallberg
- Department of Medical Sciences, Uppsala University, Uppsala, Sweden
| |
Collapse
|
7
|
Anzar I, Sverchkova A, Samarakoon P, Ellingsen EB, Gaudernack G, Stratford R, Clancy T. Personalized
HLA
typing leads to the discovery of novel
HLA
alleles and tumor‐specific
HLA
variants. HLA 2022; 99:313-327. [PMID: 35073457 PMCID: PMC9546058 DOI: 10.1111/tan.14562] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/21/2021] [Revised: 01/08/2022] [Accepted: 01/21/2022] [Indexed: 11/29/2022]
Abstract
Accurate and full‐length typing of the HLA region is important in many clinical and research settings. With the advent of next generation sequencing (NGS), several HLA typing algorithms have been developed, including many that are applicable to whole exome sequencing (WES). However, most of these solutions operate by providing the closest‐matched HLA allele among the known alleles in IPD‐IMGT/HLA Database. These database‐matching approaches have demonstrated very high performance when typing well characterized HLA alleles. However, as they rely on the completeness of the HLA database, they are not optimal for detecting novel or less well characterized alleles. Furthermore, the database‐matching approaches are also not adequate in the context of cancer, where a comprehensive characterization of somatic HLA variation and expression patterns of a tumor's HLA locus may guide therapy and clinical outcome, because of the pivotal role HLA alleles play in tumor antigen recognition and immune escape. Here, we describe a personalized HLA typing approach applied to WES data that leverages the strengths of database‐matching approaches while simultaneously allowing for the discovery of novel HLA alleles and tumor‐specific HLA variants, through the systematic integration of germline and somatic variant calling. We applied this approach on WES from 10 metastatic melanoma patients and validated the HLA typing results using HLA targeted NGS sequencing from patients where at least one HLA germline candidate was detected on Class I HLA. Targeted NGS sequencing confirmed 100% performance for the 1st and 2nd fields. In total, five out of the six detected HLA germline variants were because of Class I ambiguities at the third or fourth fields, and their detection recovered the correct HLA allele genotype. The sixth germline variant let to the formal discovery of a novel Class I allele. Finally, we demonstrated a substantially improved somatic variant detection accuracy in HLA alleles with a 91% of success rate in simulated experiments. The approach described here may allow the field to genotype more accurately using WES data, leading to the discovery of novel HLA alleles and help characterize the relationship between somatic variation in the HLA region and immunosurveillance.
Collapse
Affiliation(s)
- Irantzu Anzar
- NEC OncoImmunity AS, Oslo Cancer Cluster, Ullernchausseen 64/66, 0379 Oslo Norway
| | - Angelina Sverchkova
- NEC OncoImmunity AS, Oslo Cancer Cluster, Ullernchausseen 64/66, 0379 Oslo Norway
| | - Pubudu Samarakoon
- NEC OncoImmunity AS, Oslo Cancer Cluster, Ullernchausseen 64/66, 0379 Oslo Norway
| | | | - Gustav Gaudernack
- Ultimovacs ASA, Oslo Cancer Cluster, Ullernchausseen 64/66 Oslo Norway
| | - Richard Stratford
- NEC OncoImmunity AS, Oslo Cancer Cluster, Ullernchausseen 64/66, 0379 Oslo Norway
| | - Trevor Clancy
- NEC OncoImmunity AS, Oslo Cancer Cluster, Ullernchausseen 64/66, 0379 Oslo Norway
| |
Collapse
|
8
|
Nordin J, Pettersson M, Rosenberg LH, Mathioudaki A, Karlsson Å, Murén E, Tandre K, Rönnblom L, Kastbom A, Cedergren J, Eriksson P, Söderkvist P, Lindblad-Toh K, Meadows JRS. Association of Protective HLA-A With HLA-B∗27 Positive Ankylosing Spondylitis. Front Genet 2021; 12:659042. [PMID: 34335681 PMCID: PMC8320510 DOI: 10.3389/fgene.2021.659042] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/26/2021] [Accepted: 06/09/2021] [Indexed: 11/21/2022] Open
Abstract
Objectives To further elucidate the role of the MHC in ankylosing spondylitis by typing 17 genes, searching for HLA-B∗27 independent associations and assessing the impact of sex on this male biased disease. Methods High-confidence two-field resolution genotyping was performed on 310 cases and 2196 controls using an n-1 concordance method. Protein-coding variants were called from next-generation sequencing reads using up to four software programs and the consensus result recorded. Logistic regression tests were applied to the dataset as a whole, and also in stratified sets based on sex or HLA-B∗27 status. The amino acids driving association were also examined. Results Twenty-five HLA protein-coding variants were significantly associated to disease in the population. Three novel protective associations were found in a HLA-B∗27 positive population, HLA-A∗24:02 (OR = 0.4, CI = 0.2–0.7), and HLA-A amino acids Leu95 and Gln156. We identified a key set of seven loci that were common to both sexes, and robust to change in sample size. Stratifying by sex uncovered three novel risk variants restricted to the female population (HLA-DQA1∗04.01, -DQB1∗04:02, -DRB1∗08:01; OR = 2.4–3.1). We also uncovered a set of neutral variants in the female population, which in turn conferred strong effects in the male set, highlighting how population composition can lead to the masking of true associations. Conclusion Population stratification allowed for a nuanced investigation into the tightly linked MHC region, revealing novel HLA-B∗27 signals as well as replicating previous HLA-B∗27 dependent results. This dissection of signals may help to elucidate sex biased disease predisposition and clinical progression.
Collapse
Affiliation(s)
- Jessika Nordin
- Science for Life Laboratory, Department of Medical Biochemistry and Microbiology, Uppsala University, Uppsala, Sweden.,Science for Life Laboratory, Department of Immunology, Genetics and Pathology, Uppsala University, Uppsala, Sweden
| | - Mats Pettersson
- Science for Life Laboratory, Department of Medical Biochemistry and Microbiology, Uppsala University, Uppsala, Sweden
| | - Lina Hultin Rosenberg
- Science for Life Laboratory, Department of Medical Biochemistry and Microbiology, Uppsala University, Uppsala, Sweden
| | - Argyri Mathioudaki
- Science for Life Laboratory, Department of Medical Sciences, Uppsala University, Uppsala, Sweden
| | - Åsa Karlsson
- Science for Life Laboratory, Department of Medical Biochemistry and Microbiology, Uppsala University, Uppsala, Sweden
| | - Eva Murén
- Science for Life Laboratory, Department of Medical Biochemistry and Microbiology, Uppsala University, Uppsala, Sweden
| | - Karolina Tandre
- Science for Life Laboratory, Department of Medical Sciences, Uppsala University, Uppsala, Sweden
| | - Lars Rönnblom
- Science for Life Laboratory, Department of Medical Sciences, Uppsala University, Uppsala, Sweden
| | - Alf Kastbom
- Department of Rheumatology, University Hospital Linköping, Linköping, Sweden.,Department of Biomedical and Clinical Sciences, Linköping University, Linköping, Sweden
| | - Jan Cedergren
- Department of Rheumatology, University Hospital Linköping, Linköping, Sweden.,Department of Biomedical and Clinical Sciences, Linköping University, Linköping, Sweden
| | - Per Eriksson
- Department of Rheumatology, University Hospital Linköping, Linköping, Sweden.,Department of Biomedical and Clinical Sciences, Linköping University, Linköping, Sweden
| | - Peter Söderkvist
- Department of Biomedical and Clinical Sciences, Linköping University, Linköping, Sweden
| | - Kerstin Lindblad-Toh
- Science for Life Laboratory, Department of Medical Biochemistry and Microbiology, Uppsala University, Uppsala, Sweden.,Broad Institute of MIT and Harvard, Cambridge, MA, United States
| | - Jennifer R S Meadows
- Science for Life Laboratory, Department of Medical Biochemistry and Microbiology, Uppsala University, Uppsala, Sweden
| |
Collapse
|
9
|
Andlauer TFM, Link J, Martin D, Ryner M, Hermanrud C, Grummel V, Auer M, Hegen H, Aly L, Gasperi C, Knier B, Müller-Myhsok B, Jensen PEH, Sellebjerg F, Kockum I, Olsson T, Pallardy M, Spindeldreher S, Deisenhammer F, Fogdell-Hahn A, Hemmer B. Treatment- and population-specific genetic risk factors for anti-drug antibodies against interferon-beta: a GWAS. BMC Med 2020; 18:298. [PMID: 33143745 PMCID: PMC7641861 DOI: 10.1186/s12916-020-01769-6] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 05/14/2020] [Accepted: 08/28/2020] [Indexed: 12/22/2022] Open
Abstract
BACKGROUND Upon treatment with biopharmaceuticals, the immune system may produce anti-drug antibodies (ADA) that inhibit the therapy. Up to 40% of multiple sclerosis patients treated with interferon β (IFNβ) develop ADA, for which a genetic predisposition exists. Here, we present a genome-wide association study on ADA and predict the occurrence of antibodies in multiple sclerosis patients treated with different interferon β preparations. METHODS We analyzed a large sample of 2757 genotyped and imputed patients from two cohorts (Sweden and Germany), split between a discovery and a replication dataset. Binding ADA (bADA) levels were measured by capture-ELISA, neutralizing ADA (nADA) titers using a bioassay. Genome-wide association analyses were conducted stratified by cohort and treatment preparation, followed by fixed-effects meta-analysis. RESULTS Binding ADA levels and nADA titers were correlated and showed a significant heritability (47% and 50%, respectively). The risk factors differed strongly by treatment preparation: The top-associated and replicated variants for nADA presence were the HLA-associated variants rs77278603 in IFNβ-1a s.c.- (odds ratio (OR) = 3.55 (95% confidence interval = 2.81-4.48), p = 2.1 × 10-26) and rs28366299 in IFNβ-1b s.c.-treated patients (OR = 3.56 (2.69-4.72), p = 6.6 × 10-19). The rs77278603-correlated HLA haplotype DR15-DQ6 conferred risk specifically for IFNβ-1a s.c. (OR = 2.88 (2.29-3.61), p = 7.4 × 10-20) while DR3-DQ2 was protective (OR = 0.37 (0.27-0.52), p = 3.7 × 10-09). The haplotype DR4-DQ3 was the major risk haplotype for IFNβ-1b s.c. (OR = 7.35 (4.33-12.47), p = 1.5 × 10-13). These haplotypes exhibit large population-specific frequency differences. The best prediction models were achieved for ADA in IFNβ-1a s.c.-treated patients. Here, the prediction in the Swedish cohort showed AUC = 0.91 (0.85-0.95), sensitivity = 0.78, and specificity = 0.90; patients with the top 30% of genetic risk had, compared to patients in the bottom 30%, an OR = 73.9 (11.8-463.6, p = 4.4 × 10-6) of developing nADA. In the German cohort, the AUC of the same model was 0.83 (0.71-0.92), sensitivity = 0.80, specificity = 0.76, with an OR = 13.8 (3.0-63.3, p = 7.5 × 10-4). CONCLUSIONS We identified several HLA-associated genetic risk factors for ADA against interferon β, which were specific for treatment preparations and population backgrounds. Genetic prediction models could robustly identify patients at risk for developing ADA and might be used for personalized therapy recommendations and stratified ADA screening in clinical practice. These analyses serve as a roadmap for genetic characterizations of ADA against other biopharmaceutical compounds.
Collapse
Affiliation(s)
- Till F M Andlauer
- Department of Neurology, Klinikum rechts der Isar, School of Medicine, Technical University of Munich, Ismaninger Str 22, 81675, Munich, Germany.
- Max Planck Institute of Psychiatry, Kraepelinstr 2-10, 80804, Munich, Germany.
| | - Jenny Link
- Department of Clinical Neuroscience, Karolinska Institutet, Visionsgatan 18, 17176, Stockholm, Sweden
| | - Dorothea Martin
- Department of Neurology, Klinikum rechts der Isar, School of Medicine, Technical University of Munich, Ismaninger Str 22, 81675, Munich, Germany
| | - Malin Ryner
- Department of Clinical Neuroscience, Karolinska Institutet, Visionsgatan 18, 17176, Stockholm, Sweden
| | - Christina Hermanrud
- Department of Clinical Neuroscience, Karolinska Institutet, Visionsgatan 18, 17176, Stockholm, Sweden
| | - Verena Grummel
- Department of Neurology, Klinikum rechts der Isar, School of Medicine, Technical University of Munich, Ismaninger Str 22, 81675, Munich, Germany
| | - Michael Auer
- Department of Neurology, Medical University of Innsbruck, Anichstr 35, 6020, Innsbruck, Austria
| | - Harald Hegen
- Department of Neurology, Medical University of Innsbruck, Anichstr 35, 6020, Innsbruck, Austria
| | - Lilian Aly
- Department of Neurology, Klinikum rechts der Isar, School of Medicine, Technical University of Munich, Ismaninger Str 22, 81675, Munich, Germany
- Institute of Experimental Neuroimmunology, Technical University of Munich, Trogerstr 9, 81675, Munich, Germany
| | - Christiane Gasperi
- Department of Neurology, Klinikum rechts der Isar, School of Medicine, Technical University of Munich, Ismaninger Str 22, 81675, Munich, Germany
| | - Benjamin Knier
- Department of Neurology, Klinikum rechts der Isar, School of Medicine, Technical University of Munich, Ismaninger Str 22, 81675, Munich, Germany
- Institute of Experimental Neuroimmunology, Technical University of Munich, Trogerstr 9, 81675, Munich, Germany
| | - Bertram Müller-Myhsok
- Max Planck Institute of Psychiatry, Kraepelinstr 2-10, 80804, Munich, Germany
- Institute of Translational Medicine, University of Liverpool, Crown Street, Liverpool, L69 3BX, UK
- Munich Cluster for Systems Neurology (SyNergy), Feodor-Lynen-Str. 17, 81377, Munich, Germany
| | | | - Finn Sellebjerg
- DMSC, Department of Neurology, Rigshospitalet, University of Copenhagen, 2100, Copenhagen, Denmark
| | - Ingrid Kockum
- Department of Clinical Neuroscience, Karolinska Institutet, Visionsgatan 18, 17176, Stockholm, Sweden
| | - Tomas Olsson
- Department of Clinical Neuroscience, Karolinska Institutet, Visionsgatan 18, 17176, Stockholm, Sweden
| | - Marc Pallardy
- Inflammation, Microbiome and Immunosurveillance, Université Paris-Saclay, INSERM, Faculté de Pharmacie, rue JB Clément, 92290, Châtenay-Malabry, France
| | - Sebastian Spindeldreher
- Novartis Institutes for Biomedical Research, Novartis Pharma AG, 4056, Basel, Switzerland
- Integrated Biologix GmbH, Steinenvorstadt 33, 4051, Basel, Switzerland
| | - Florian Deisenhammer
- Department of Neurology, Medical University of Innsbruck, Anichstr 35, 6020, Innsbruck, Austria
| | - Anna Fogdell-Hahn
- Department of Clinical Neuroscience, Karolinska Institutet, Visionsgatan 18, 17176, Stockholm, Sweden
| | - Bernhard Hemmer
- Department of Neurology, Klinikum rechts der Isar, School of Medicine, Technical University of Munich, Ismaninger Str 22, 81675, Munich, Germany.
- Munich Cluster for Systems Neurology (SyNergy), Feodor-Lynen-Str. 17, 81377, Munich, Germany.
| |
Collapse
|