Juliane Russlies1
Anke Fähnrich1,2
Mareike Witte1,3
Junping Yin4
Sandrine Benoit5,6Regine Gläser6,7†
Claudia Günter6,8
Rüdiger Eming6,9Jeanette Erdmann2
Damian Gola10
Yask Gupta1Maike Marleen Holtsche3Johannes S. Kern6,11,12
Inke R. König10
Dimitra Kiritsi6,11
Wolfgang Lieb13,14
Christian D. Sadik3Miklós Sárdy6,15,16
Franziska Schauer6,11Nina van Beek3Anke Weidinger6,7Margitta Worm6,17
Detlef Zillikens3,6
Enno Schmidt1,3,6
Hauke Busch1,2
Saleh M. Ibrahim1,6
Misa Hirose1*- 1Luebeck Institute of Experimental Dermatology, University of Luebeck, Luebeck, Germany
- 2Institute for Cardiogenetics, University of Luebeck, Luebeck, Germany
- 3Department of Dermatology, University of Luebeck, Luebeck, Germany
- 4Research Center Borstel, Leibniz-Center for Medicine and Bioscience, Borstel, Germany
- 5Department of Dermatology, Venereology and Allergology, University Hospital Würzburg, Würzburg, Germany
- 6The German Autoimmune Bullous Disease Genetic Study Group, Germany
- 7Department of Dermatology, Venereology and Allergology, University Hospital Schleswig-Holstein, Kiel, Germany
- 8Department of Dermatology, University Hospital of Dresden, Dresden, Germany
- 9Department of Dermatology and Allergology, Phillips-Universität Marburg, Marburg, Germany
- 10Institute of Medical Biometry and Statistics, University of Luebeck, Luebeck, Germany
- 11Department of Dermatology, Faculty of Medicine, Medical Center-University of Freiburg, University of Freiburg, Freiburg, Germany
- 12Dermatology Department, Faculty of Medicine, Dentistry and Health Sciences, The Royal Melbourne Hospital, University of Melbourne, Melbourne, VIC, Australia
- 13Institute of Epidemiology, Christian-Albrecht University of Kiel, Kiel, Germany
- 14Popgen Biobank, Christian-Albrecht University of Kiel, Kiel, Germany
- 15Department of Dermatology, Venereology and Dermatooncology, Semmelweis University, Budapest, Hungary
- 16Department of Dermatology and Allergy, University Hospital, LMU Munich, Munich, Germany
- 17Department of Dermatology, Venereology and Allergology, Allergy Center Charité, Charité Medical University Berlin, Berlin, Germany
Bullous pemphigoid (BP) is the most prevalent autoimmune skin blistering disease and is characterized by the generation of autoantibodies against the hemidesmosomal proteins BP180 (type XVII collagen) and BP230. Most intriguingly, BP is distinct from other autoimmune diseases because it predominantly affects elderly individuals above the age of 75 years, raising the question why autoantibodies and the clinical lesions of BP emerges mostly in this later stage of life, even in individuals harboring known putative BP-associated germline gene variants. The mitochondrial genome (mtDNA) is a potential candidate to provide additional insights into the BP etiology; however, the mtDNA has not been extensively explored to date. Therefore, we sequenced the whole mtDNA of German BP patients (n = 180) and age- and sex-matched healthy controls (n = 188) using next generation sequencing (NGS) technology, followed by the replication study using Sanger sequencing of an additional independent BP (n = 89) and control cohort (n = 104). While the BP and control groups showed comparable mitochondrial haplogroup distributions, the haplogroup T exhibited a tendency of higher frequency in BP patients suffering from neurodegenerative diseases (ND) compared to BP patients without ND (50%; 3 in 6 BP with haplogroup T). A total of four single nucleotide polymorphisms (SNPs) in the mtDNA, namely, m.16263T>C, m.16051A>G, and m.16162A>G in the D-loop region of the mtDNA, and m.11914G>A in the mitochondrially encoded NADH:ubiquinone oxidoreductase core subunit 4 gene (MT-ND4), were found to be significantly associated with BP based on the meta-analysis of our NGS data and the Sanger sequencing data (p = 0.0017, p = 0.0129, p = 0.0076, and p = 0.0132, respectively, Peto's test). More specifically, the three SNPs in the D-loop region were negatively, and the SNP in the MT-ND4 gene was positively associated with BP. Our study is the first to interrogate the whole mtDNA in BP patients and controls and to implicate multiple novel mtDNA variants in disease susceptibility. Studies using larger cohorts and more diverse populations are warranted to explore the functional consequences of the mtDNA variants identified in this study on immune and skin cells to understand their contributions to BP pathology.
Introduction
A number of studies to identify candidate genes in autoimmune blistering skin diseases, particularly bullous pemphigoid (BP) have been conducted to date. These studies predominantly identified associated gene polymorphisms in immune system-related genes, e.g., HLA region (1–6), Fc gamma receptor genes (7), and cytokine genes (8). In addition to these nuclear-encoded genes, our group has recently shown that polymorphisms in a gene encoded in the mitochondrial genome, MT-ATP8, are associated with BP (9). Considering the nature of autoimmune diseases, the causal factors of BP are not only limited to genetics, but also involve environmental factors. For example, our recent findings showed that the composition of skin microbiota was altered in BP patients compared to healthy controls (10).
The mitochondrial genome (mtDNA) is a circular DNA molecule with a length of ~16 kilobase pairs. The mtDNA encodes 13 protein-coding genes, 22 transfer RNA genes, and 2 ribosomal RNA genes (11, 12). Multiple copies of the mtDNA exist in a single mitochondrion. All of the 13 mtDNA-encoded proteins consist of subunits of the oxidative phosphorylation (OXPHOS) complexes, which are responsible for cellular energy production in the form of ATP, as well as for the production of reactive oxygen species (ROS) as a by-product of the OXPHOS reaction. The mtDNA is polymorphic, and variations in the mtDNA are known to be associated with alterations in mitochondrial functions (13). The mtDNA variations in humans are categorized into the following three groups: (1) recent maternally inherited deleterious mutations; (2) ancient adaptive polymorphisms; and (3) somatic mutations that accumulate during development and in tissues with age (14). Recent maternally inherited deleterious mtDNA mutations have been well-described in rare mitochondrial disorders, including familial mitochondrial encephalomyopathy (15) and Leber's hereditary optic neuritis (16). Ancient adaptive polymorphisms are commonly used to establish haplogroup ancestry, as these variations are believed to have occurred during the migration of human ancestors in order to adapt different environments (e.g., nutritional availability and climates) (17). Furthermore, several ancient polymorphisms in the mtDNA have been reported to be associated with common diseases, including chronic inflammation and autoimmune diseases (18, 19). These are not surprising because the fate and the function of immune cells are largely determined by cellular metabolism (immunometabolism), which is to a large extent controlled by mitochondrial functions (20–22). One example of such immune cell types is regulatory T cells, which have been reported to be involved in the pathology of BP (23–25), and their differentiation is determined by the levels of fatty acid synthesis, one of the mitochondrial functions (26). Longevity and aging have also been associated with certain mtDNA polymorphisms (27–29). Another unique characteristic of the mtDNA is the higher frequency of somatic mutations in aging compared to the nuclear genome (30), indicating the presence of variations identified only in elderly people. In fact, BP is the most prevalent autoimmune blistering skin disease and predominantly affects the elderly population, i.e., usually in late 70s (31, 32). As mentioned above, our group has recently demonstrated changes in the skin microbiota composition having been observed in BP patients compared to healthy controls (10). Recently, certain mtDNA haplogroups have been reported to be associated with the abundance of certain bacterial taxa (33). Consistent with these findings, our group recently identified that variations in the mtDNA are associated with the composition of microbiota in the gut (34) and the skin of mice (unpublished), suggesting that the variations in the mtDNA observed in BP patients contribute to a shift of the skin microbiota composition, which in turn enhances susceptibility to the disease. All of the abovementioned characteristics of the mtDNA support its potential involvement in BP.
Therefore, we have explored here the whole mitochondrial genome by next generation sequencing technology in German BP patients and their age- and sex-matched controls.
Materials and Methods
Study Cohorts
DNA samples for the NGS discovery study and the Sanger sequencing replication study were obtained from the German AIBD Genetics Study Group and PopGen Biobank. BP patients were diagnosed by clinicians at the participating centers. All included patients satisfied all of the following criteria; (i) a compatible clinical presentation, (ii) the detection of linear deposits of IgG and/or C3 at the dermal-epidermal junction by direct immunofluorescence (IF) microscopy of a perilesional skin biopsy, and (iii) the detection of serum autoantibodies against BP180 NC16A and/or BP230 by ELISA, according to the guideline of the German Dermatological Society for the diagnosis of BP (35). DNA samples from a total of 270 BP patients (180 for the NGS and 90 for a replication study using Sanger sequencing) and 294 controls (188 for the NGS and 106 for Sanger sequencing) were used in this study. The description of the cohort is summarized in Supplementary Table 1.
Of the samples tested by the NGS, 180 BP patients and 144 controls were evaluated for their genetic ancestry using their genome-wide SNP data, which were previously obtained using Affymetrix Biobank Axiom ArrayTM (Thermo Fisher Scientific, MA, USA). Principal components analysis (PCA) (Supplementary Figure 1) plots of the 180 BP and 144 control samples almost overlapped, suggesting that the BP patients and controls in this study belonged to the same population.
Of the 180 BP patients whose mtDNA was sequenced, clinical history of neurodegenerative disease (dementia and Parkinson's disease) was available in 58 patients. Of the 58 BP patients, data on their anti-BP180 NC16A IgG titer were available in 49 patients.
The study was approved by the ethical committees of the University of Lübeck (10-026 and 15-051) and the individual study centers.
Next Generation Sequencing of the Whole Mitochondrial Genome
Genomic DNA samples were processed for library preparation, as previously described in the Human mtDNA Genome protocol for Illumina Sequencing Platform (36). In brief, two primer sets [MTL-F1 (AAAGCACATACCAAGGCCAC) and MTL-R1 (TTGGCTCTCCTTGCAAAGTT); MTL-F2 (TATCCGCCATCCCATACATT) and MTL-R2 (AATGTTGAGCCGTAGATGCC)] were used to amplify the mtDNA by long-range PCR. Library preparation was performed using a Nextera XT DNA Library Preparation Kit (Illumina Inc., CA, USA), and the 10-pM library was sequenced on the Illumina MiSeq sequencing platform (2 × 150 bp paired-end reads) (Illumina Inc.).
NGS Data Analysis
Our previously described data analysis method (37) was modified and adapted for human mtDNA analysis. After quality control, the reads were mapped to the revised Cambridge Reference Sequence (rCRS; NC_012920.1) using Burrows-Wheeler Aligner bwa version 0.705 (38), and bam files were generated. Duplicated reads generated during PCR were removed using Markduplicates (Picard tools version 1.119) (39), and indels were realigned using IndelRealigner (Genome analysis tool kit version 3.3) (40). The processed bam files were assessed for frequency and base quality (≥30) for each reference and alternate base in the mtDNA using pysamstats (version 0.24.3) (41). When the frequency of the alternate allele compared with the reference allele was >90%, it was considered as homoplasmic mutation, whereas the 10–90% range was considered as heteroplasmy. Additionally, bam files were manually inspected for the presence of mutations and indels using IGV software (42). mtDNA variants were annotated using MSeqDR mvTool (43), a DNA Web resource for comprehensive variant annotation.
Mitochondrial Haplogroup Analysis
Mitochondrial haplogroup assignment was conducted using HaploGrep 2 (44). In brief, HaploGrep weighs each polymorphism present in PhyloTree17 (45) based on its informativeness to define haplogroups. The set of SNPs in the input file were classified as informative or remaining (not informative). A score is given based on the weights of the informative SNPs, and the offset was determined based on the number of remaining SNPs.
Replication Study (Sanger Sequencing)
DNA samples for the replication study were prepared using standard DNA extraction kits (Qiagen, Hilden, Germany). SNP regions were amplified by standard PCR. The primers used for the PCR reaction are listed in Supplementary Table 2. The PCR products were sent to Genewiz (Essex, UK) for Sanger sequencing, and the obtained data were analyzed using the freely available software Unipro UGENE (46).
Statistical Analysis
Data from the mtSNP association study and the mitochondrial haplogroup association study were analyzed using R package “exact2 × 2,” which provides a non-central confidence interval matching the two-sided Fisher's exact test based on the principle of likelihood estimation (47).
The meta-analysis was conducted using Peto's method (48) from the R package “metaphor” (49). This method provides a weighted estimate of the log odds ratio under a fixed-effects model. We used Hommel's method (50) to give strong control of the family-wise error rate, i.e., the probability of at least one type I error, by adjusting each p-value obtained from Peto's method.
Statistical analyses for other studies were performed using GraphPad Prism 6 (GraphPad Software, San Diego, CA, USA). Statistical tests used for the analysis are indicated in the figure legends.
Results
BP Patients in the Mitochondrial Haplogroup T May Exhibit a Higher Co-incidence With Neurodegenerative Conditions
Mitochondrial haplogroup analysis in this German population revealed that 47.28% of the sequenced individuals belonged to the haplogroup H, which is the major haplogroup in Europeans, followed by haplogroup U (18.21%), haplogroup J (10.05%), and haplogroup T (8.42%). When the data were analyzed for disease association, there was no association between BP status and the mtDNA haplogroups (p = 0.7963, Fisher's exact test, Figure 1A, Supplementary Table 3).
Figure 1. Mitochondrial haplogroup distribution in BP patients and controls. (A) BP patients and controls showed similar mitochondrial haplogroup distributions. (B) Mitochondrial haplogroup distribution when BP patients were stratified for clinical history of neurodegenerative diseases (ND). BP patients with the haplogroup T background tended to have higher incidence of ND.
Among the 180 BP patients sequenced for the whole mtDNA, 58 cases had clinical history of neurodegenerative diseases (ND; e.g., Parkinson's disease and dementia), consistent with recent studies that reported an association between neurodegenerative diseases and BP in different populations (51). We analyzed the mitochondrial haplogroup associations in these 58 BP patients and found that BP patients in the mitochondrial haplogroup T tend to have higher risk of ND (Figure 1B, p = 0.1448, Fisher's exact test). We additionally evaluated 49 BP patients with available records of autoantibody (anti-BP180 NC16A IgG) levels. While no associations between the autoantibody titers and mitochondrial haplogroups were observed (Figure 2A), BP patients with concurrent ND showed significantly higher variation in autoantibody titers compared to BP patients without ND (Figure 2B, p = 0.027, Mann–Whitney test).
Figure 2. (A) The levels of anti-BP180 antibodies in individuals with each haplogroup. The mean values of anti-BP180 IgG were 223.8 U/L in haplogroup H, 101 U/L in haplogroup U, 326.7 U/L in haplogroup T, 106.4 U/L in haplogroup J, 144 U/L in haplogroup I, 113 U/L in haplogroup V, 27 U/L in haplogroup L, 114 U/L in haplogroup W. (B) The levels of the autoantibodies BP180 NC16A were highly varied in BP patients who also suffer from neurodegenerative diseases. P = 0.0272, Mann–Whitney test.
Novel Candidate SNPs in the mtDNA-Associated With BP in Germans
A total of 1,010 SNPs in the mtDNA (mtSNPs) were identified in this study. The mtSNPs associated with BP (exploratory p < 0.1) are listed in Table 1. Of the five listed variants, the three mtSNPs; m.16263T>C, m.11914G>A, and m.15904C>T; were selected for the replication study using Sanger sequencing in the replication cohort of an additional independent 90 BP and 106 control samples. In this replication study, two relevant mtDNA sequences covering all three mtSNPs were PCR-amplified (Supplementary Table 2) and processed for Sanger sequencing. The Sanger sequencing data confirmed the validity of the NGS results obtained from the individuals carrying a variant from each of the three mtSNPs. In addition to the analysis of the three targeted SNPs, the design of the primers used for the Sanger sequencing replication study enabled us to evaluate the sequencing data of the locus ranging from m.15800 to m.16290 (covering a part of the displacement-loop region; D-loop region), as well as the locus ranging from m.11870 to m.12137 (covering a part of the mitochondrially encoded NADH:ubiquinone oxidoreductase core subunit 4 gene; MT-ND4) in the replication cohorts.
Table 1. Five candidate BP-associated SNPs in the mtDNA identified by the next-generation sequencing were selected for the replication study.
The meta-analysis of the NGS and the Sanger sequencing data revealed that m.16263T>C, m.11914G>A, m.16051A>G, and m.16162A>G were significantly associated with BP in Germans. Variants in the three mtSNPs in the D-loop region were more frequent in controls, while the variant in the m.11914G>A in the MT-ND4 gene was more frequent in BP (Table 2). Meta-analysis revealed no significant association between m.15904C>T and BP (p = 0.2629, Peto's test).
Table 2. Meta-analysis of NGS data and Sanger sequencing data of BP patients and controls.
Discussion
In this study, we performed NGS of the whole mitochondrial genomes of 180 BP patients and 188 controls, followed by a replication study using Sanger sequencing of independent samples of up to 90 BP patients and 104 controls from Germany. With cohorts tested in this study, we estimated the power of our study to detect mitochondrial variants. For the NGS study to detect SNPs to forward to replication, we used threshold of alpha = 0.1. Based on this, we estimated that for variants with frequencies of about 5% in controls, odds ratios of about 2.6 and higher would be detectable at a power of 80%; if the variants are more frequent with about 10% in controls, odds ratios of about 2.2 and higher are detectable reliability. For the replication study, Sanger sequencing on 90 independent cases and 104 independent controls was performed; applying a replication significance level of 0.05, this study was well-powered with a power of at least 80% to detect odds ratios of at least 4.2. This might explain why we were not able to replicate the variant m.15904C>T because the power for this SNP might have been too low. The study cohort is the largest for this disease worldwide, however, it is still relatively small in size compared to those used in genetic studies of other chronic inflammatory skin diseases, e.g., psoriasis or atopic dermatitis. This is primarily due to the rarity of the disease. Nevertheless, this study is highly valuable to the field as it is the first to interrogate the whole mitochondrial genome by sequencing in autoimmune blistering skin diseases.
The whole mtDNA sequencing data were analyzed for mitochondrial haplogroup associations with BP. As mentioned above, mitochondrial haplogroups are defined by specific combinations of ancient adaptive polymorphisms in the mtDNA and often reflect on mitochondrial functionality to adapt to the specific environments according to the geographical locations where our ancestors migrated. Therefore, mitochondrial haplogroups have been used to define ethnic origins on mostly prehistoric time scales. In fact, as previously shown (52), combinations of the adaptive polymorphisms in the mtDNA altered the mitochondrial functions, which are causal for common complex diseases in certain populations.
Our findings showed that the haplogroup distribution between BP and controls in this study was comparable. However, interestingly, when we correlated the mitochondrial haplogroup and the presence of neurodegenerative diseases (ND; i.e., dementia and Parkinson's disease) in BP patients, the analysis revealed that the haplogroup T appeared more frequent in BP patients suffering from ND comorbidity. Haplogroup T is known to be associated with ND (53, 54). In parallel, the levels of anti-BP180 antibodies in individuals with haplogroup T exhibited a tendency of higher levels though the result is exploratory (average values of 326.7 U/L in individuals with haplogroup T, while the average values in all BP patients were 175.6 U/L). Together with the positive association between the levels of anti-BP180 antibodies and the presence of ND, the link between mitochondrial haplogroup, neurodegenerative disease and the levels of anti-BP180 autoantibodies among BP patients may be plausible. The haplogroup J, in contrast, did not share the same tendency, even though haplogroup J and T belong to the same sub-cradle of the haplogroup JT. These results are in agreement with those of previous studies demonstrating that the haplogroup J is protective against Parkinson's disease (18, 55), which suggested that the mtDNA variants defining the haplogroup J may be protective from concurrent ND in BP patients. Furthermore, the average levels of autoantibodies in individuals with haplogroup J were 106.4 U/L, which were lower than the average values of all BP patients in this analysis. Nevertheless, to confirm this exploratory observation, further studies should analyze a larger patient cohort with available clinical histories.
Next, whole mtDNA sequencing data were analyzed for associations between the single nucleotide polymorphisms in the mtDNA (mtSNPs) and BP. The analysis identified five top candidate mtSNPs. The meta-analysis of the discovery study (whole mtDNA NGS data) and the replication study (Sanger sequencing data of partial mtDNA region) revealed four mtSNPs that were significantly associated with BP. The three mtSNPs located in the non-coding D-loop region, namely, m.16051A>G, m.16162A>G, and m.16263T>C, were all enriched in the controls, while the mtSNP m11914 in the MT-ND4 gene was more frequent in BP patients. The D-loop is a non-coding region of the mtDNA, including the replication origin of the H (heavy) strand (OriH), the promoters for transcription of the H and L (light) strand (HSP and LSP) and two hyper variable segments (HVS1: m.16,024-m.16,383, and HVS2: m.57-m.372), constituting the most variable regions in the mtDNA (56). All three mtSNPs in the D-loop region are located in HVS1. Functional consequences of the variants in this non-coding D-loop region, particularly the HVS, remain unknown to date. Interestingly, most carriers of a variant from any of these three mtSNPs in the D-loop region belong to haplogroup H; two of the 81 BP patients belonged to haplogroup H, and 14 of 93 controls belonged to haplogroup H. Given that the three mtSNPs serve as defining SNPs for haplogroup H-subgroups, i.e., m.16051A>G for H1a3 and H2a2a1c (45), the functional relevance of these mtSNPs may be linked with other haplogroup H-subgroups-defining SNPs. Considering the limited sample size of patients with a rare disease, a more finely-branched subgroup haplogroup analysis cannot be conducted in this study. This also applies to the mtSNP m.11914G>A. The A variant enriched in the BP group is a synonymous mutation in the MT-ND4 gene, suggesting that its functional relevance is unlikely. The variant in m.11914G>A is also a defining SNP for several haplogroups-subgroups U, K, T, and H (45). Three of 12 BP patients and one of 15 controls belonging to haplogroup K carried the A variant in m.11914G>A.
Previously, our group has reported an association between a rare variant in the MT-ATP8 gene and BP in Germans (9). The current study is an extension of the previous study to explore the whole mtDNA in a German cohort, including other individuals than those enrolled in the previous study. The BP-associated mtSNP, m.8519G>A, which was identified in our previous study, was found in one BP patient among newly sequenced 137 BP patients and 20 controls in the current study (p = 1.000, Fisher's exact test). The meta-analysis of the NGS data and the previously published Sanger sequencing data still showed a positive association between m.8519G>A and BP (p = 0.0151, Peto's test, odds ratio 7.3454, 95% confidence interval 1.4718–36.6594).
We utilized peripheral blood DNA samples in this study. Somatic mutations accumulate over time in tissues with age (17), and the levels of the mutant mtDNA, i.e., levels of heteroplasmy, can be different between various tissues and organs (30, 57). Therefore, an evaluation of the whole mtDNA sequence in other tissues of importance, i.e., skin samples obtained from BP patients and controls, is warranted. Within BP skin samples, site different samples, i.e., peri-lesional and unaffected, are also of great interest to evaluate differential levels of the mtDNA mutations within the same individual. Such age-dependent and tissue-specific changes in the mtDNA may elucidate the pathways in late-onset diseases such as BP. Thus, the interpretation of the results in this study needs to be cautious as the results obtained from this study using peripheral blood DNA might not be the same as those using tissue DNA. As aforementioned, BP is a multifactorial disease, and the involvement of immune cells and skin microbiota in the disease has been proposed (10, 58). Both the immune system and the composition of microbiota are altered in aging (59–62), interact with each other (63), and are associated with mitochondrial functions (21, 64). Complex interactions between these age-related alterations and age- and tissue-specific mtDNA variants could provide key insights to the pathways in diseases with complex traits. To date, only a few studies using mice proposed an association between mtDNA variants, microbiota and clinical phenotypes have been reported (34, 65), but none in humans.
In summary, we investigated genetic variants in the whole mtDNA genome in German BP patients and their age- and sex-matched controls, which is currently the largest available study cohort for BP worldwide. Our findings showed that the maternally inherited natural variants in the mtDNA are associated with BP, which predominantly affects the elderly population. Therefore, more complex interactions between nuclear genome variants and mtDNA variants, as well as aging, are likely to be involved in the pathogenesis of BP. To identify the functionally relevant mtDNA variants in BP, studies with the larger sample sizes and analysis of the mtDNA genome in the skin should be conducted.
Data Availability
Sequencing data used for this study is available under controlled access from the European Genome-phenome Archive with the Study ID EGAS00001003932 (https://ega-archive.org/studies/EGAS00001003932).
Ethics Statement
This study was conducted in accordance with the recommendations of institutional review board of the University of Lübeck, and the review boards of the collaborating centers. The protocol was approved by the institutional review board of the University of Lübeck (File No. 10-026 and 15-051), and the review boards of the collaborating centers. All subjects provided written informed consent in accordance with the Declaration of Helsinki. Samples and demographic data of patients and controls were collected in adherence to ethics and German privacy protection regulations.
Author Contributions
JR, MWi, and JY performed the experiment. SB, RG, CG, RE, MMH, JK, DK, CS, MS, FS, NB, AW, MWo, and DZ recruited the patients and control cohorts. ES initiated and coordinated the German AIBD Genetics Study group. WL provided POPGEN control cohort DNA samples. AF, HB, JE, and YG performed the NGS data process. JR, AF, HB, DG, and IK conducted statistical analysis. JR, SI, and MH interpreted the data and wrote the manuscript with support of all co-authors. SI and MH designed the study.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
The authors thank everyone involved in the collection of blood samples, as well as all volunteers and patients, who donated blood for the study. We also thank Miriam Freitag, Paul Schilf, Christian Schwarm, and Stefanie Wichmann for advice and/or excellent technical assistance. AF and HB acknowledge support through the high-performance computer cluster of the University of Lübeck OmicsCluster. This work was supported by the Deutsche Forschungsgemeinschaft (EXC306, KFO303, MP-II, and EXC 22167-390884018).
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2019.02200/full#supplementary-material
References
1. Esmaili N, Mortazavi H, Chams-Davatchi C, Daneshpazhooh M, Damavandi MR, Aryanian Z, et al. Association between HLA-DQB1*03:01 and Bullous pemphigoid in Iranian patients. Iran J Immunol. (2013) 10:1–9.
2. Zakka LR, Reche P, Ahmed AR. Role of MHC Class II genes in the pathogenesis of pemphigoid. Autoimmun Rev. (2011) 11:40–7. doi: 10.1016/j.autrev.2011.07.002
3. Gao XH, Winsey S, Li G, Barnardo M, Zhu XL, Chen HD, et al. HLA-DR and DQ polymorphisms in bullosu pemphigoid from northern China. Clin Exp Dermatol. (2002) 27:319–21. doi: 10.1046/j.1365-2230.2002.01037.x
4. Okazaki A, Miyagawa S, Yamashina Y, Kitamura W, Shirai T. Polymorphisms of HLA-DR and -DQ genes in Japanese patients with bullous pemphigoid. J Dermatol. (2000) 27:149–56. doi: 10.1111/j.1346-8138.2000.tb02141.x
5. Delgado JC, Turbay D, Yunis EJ, Yunis JJ, Morton ED, Bhol K, et al. A common major histocompatibility complex class II allele HLA-DQB1* 0301 is present in clinical variants of pemphigoid. Proc Natl Acad Sci USA. (1996) 93:8569–71. doi: 10.1073/pnas.93.16.8569
6. Seignalet J, Guillot B, Guilhou JJ, Meynadier J. Probable association between HLA-DR5 and bullous pemphigoid. Tissue Antigens. (1987) 30:190–1. doi: 10.1111/j.1399-0039.1987.tb01620.x
7. Weisenseel P, Martin S, Partscht K, Messer G, Prinz JC. Relevance of the low-affinity type of the Fcgamma-receptor IIIa-polymorphism in bullous pemphigoid. Arch Dermatol Res. (2007) 299:163–4. doi: 10.1007/s00403-007-0755-8
8. Chang YT, Liu HN, Yu CW, Lin MW, Huang CH, Chen CC, et al. Cytokine gene polymorphisms in bullous pemphigoid in a Chinese population. Br J Dermatol. (2006) 154:79–84. doi: 10.1111/j.1365-2133.2005.06938.x
9. Hirose M, Schilf P, Benoit S, Eming R, Gläser R, Homey B, et al. Polymorphisms in the mitochondrially encoded ATP synthase 8 gene are associated with susceptibility to bullous pemphigoid in the German population. Exp Dermatol. (2015) 24:715–7. doi: 10.1111/exd.12732
10. Miodovnik M, Künstner A, Langan EA, Zillikens D, Gläser R, Sprecher E, et al. A distinct cutaneous microbiota profile in autoimmune bullous disease patients. Exp Dermatol. (2017) 26:1221–7. doi: 10.1111/exd.13357
11. Anderson S, Bankier AT, Barrell BG, de Bruijn MH, Coulson AR, Drouin J, et al. Sequence and organization of the human mitochondrial genome. Nature. (1981) 290:457–65. doi: 10.1038/290457a0
12. Bibb MJ, Van Etten RA, Wright CT, Walberg MW, Clayton DA. Sequence and gene organization of mouse mitochondrial DNA. Cell. (1981) 26:167–80. doi: 10.1016/0092-8674(81)90300-7
13. Wallace DC. Mitochondrial DNA variation in human radiation and disease. Cell. (2015) 163:33–8. doi: 10.1016/j.cell.2015.08.067
14. Wallace DC. Mitochondrial genetic medicine. Nat Genet. (2018) 50:1642–9. doi: 10.1038/s41588-018-0264-z
15. Wallace DC, Zheng XX, Lott MT, Shoffner JM, Hodge JA, Kelley RI, et al. Familial mitochondrial encephalomyopathy (MERRF): genetic, pathophysiological, and biochemical characterization of a mitochondrial DNA disease. Cell. (1988) 55:601–10. doi: 10.1016/0092-8674(88)90218-8
16. Wallace DC, Singh G, Lott MT, Hodge JA, Schurr TG, Lezza AM, et al. Mitochondrial DNA mutation associated with Leber's hereditary optic neuropathy. Science. (1988) 242:1427–30. doi: 10.1126/science.3201231
17. Wallace DC. A mitochondrial paradigm of metabolic and degenerative diseases, aging, and cancer: a dawn for evolutionary medicine. Annu Rev Genet. (2005) 39:359–407. doi: 10.1146/annurev.genet.39.110304.095751
18. Hudson G, Gomez-Duran A, Wilson IJ, Chinnery PF. Recent mitochondrial DNA mutations increase the risk of developing common late-onset human diseases. PLoS Genet. (2014) 10:e1004369. doi: 10.1371/journal.pgen.1004369
19. Taylor RW, Turnbull DM. Mitochondrial DNA mutations in human disease. Nat Rev Genet. (2005) 6:389–402. doi: 10.1038/nrg1606
20. Buck MD, Sowell RT, Kaech SM, Pearce EL. Metabolic instruction of immunity. Cell. (2017) 169:570–86. doi: 10.1016/j.cell.2017.04.004
21. Mills EL, Kelly B, O'Neill LAJ. Mitochondria are the powerhouses of immunity. Nat Immunol. (2017) 18:488–98. doi: 10.1038/ni.3704
22. Freitag J, Berod L, Kamradt T, Sparwasser T. Immunometabolism and autoimmunity. Immunol Cell Biol. (2016) 94:925–34. doi: 10.1038/icb.2016.77
23. Arakawa M, Dainichi T, Ishii N, Hamada T, Karashima T, Nakama T, et al. Lesional Th17 cells and regulatory T cells in bullous pemphigoid. Exp Dermatol. (2011) 20:1022–4. doi: 10.1111/j.1600-0625.2011.01378.x
24. Antiga E, Quaglino P, Volpi W, Pierini I, Del Bianco E, Bianchi B, et al. Regulatory T cells in skin lesions and blood of patients with bullous pemphigoid. J Eur Acad Dermatol Venereol. (2014) 28:222–30. doi: 10.1111/jdv.12091
25. Haeberle S, Wei X, Bieber K, Goletz S, Ludwig RJ, Schmidt E, et al. Regulatory T-cell deficiency leads to pathogenic bullous pemphigoid antigen 230 autoantibody and autoimmune bullous disease. J Allergy Clin Immunol. (2018) 142:1831–42.e7. doi: 10.1016/j.jaci.2018.04.006
26. Berod L, Friedrich C, Nandan A, Freitag J, Hagemann S, Harmrolfs K, et al. De novo fatty acid synthesis controls the fate between regulatory T and T helper 17 cells. Nat Med. (2014) 20:1327–33. doi: 10.1038/nm.3704
27. Tanaka M, Gong JS, Zhang J, Yoneda M, Yagi K. Mitochondrial genotype associated with longevity. Lancet. (1998) 351:185–6. doi: 10.1016/S0140-6736(05)78211-8
28. Ivanova R, Lepage V, Charron D, Schächter F. Mitochondrial genotype associated with French Caucasian centenarians. Gerontology. (1998) 44:349. doi: 10.1159/000022041
29. De Benedictis G, Rose G, Carrieri G, De Luca M, Falcone E, Passarino G, et al. Mitochondrial DNA inherited variants are associated with successful aging and longevity in humans. FASEB J. (1999) 13:1532–6. doi: 10.1096/fasebj.13.12.1532
30. Wallace DC, Chalkia D. Mitochondrial DNA genetics and the heteroplasmy conundrum in evolution and disease. Cold Spring Harb Perspect Biol. (2013) 5:a021220. doi: 10.1101/cshperspect.a021220
31. Joly P, Baricault S, Sparsa A, Bernard P, Bédane C, Duvert-Lehembre S, et al. Incidence and mortality of bullous pemphigoid in France. J Invest Dermatol. (2012) 132:1998–2004. doi: 10.1038/jid.2012.35
32. Hübner F, Recke A, Zillikens D, Linder R, Schmidt E. Prevalence and age distribution of pemphigus and pemphigoid diseases in Germany. J Invest Dermatol. (2016) 136:2495–8. doi: 10.1016/j.jid.2016.07.013
33. Ma J, Coarfa C, Qin X, Bonnen PE, Milosavljevic A, Versalovic J, et al. mtDNA haplogroup and single nucleotide polymorphisms structure human microbiome communities. BMC Genomics. (2014) 15:257. doi: 10.1186/1471-2164-15-257
34. Hirose M, Künstner A, Schilf P, Sünderhauf A, Rupp J, Jöhren O, et al. Mitochondrial gene polymorphism is associated with gut microbial communities in mice. Sci Rep. (2017) 7:15293. doi: 10.1038/s41598-017-15377-7
35. Schmidt E, Goebeler M, Hertl M, Sárdy M, Sitaru C, Eming R, et al. S2k guideline for the diagnosis of pemphigus vulgaris/foliaceus and bullous pemphigoid. J Dtsch Dermatol Ges. (2015) 13:713–27. doi: 10.1111/ddg.12612
36. Human mtDNA Genome for the Illumina Sequencing Platform. Available online at: http://emea.support.illumina.com/content/dam/illumina-support/documents/documentation/chemistry_documentation/samplepreps_legacy/human-mtdna-genome-guide-15037958-01.pdf (accessed January 30, 2019).
37. Hirose M, Schilf P, Gupta Y, Zarse K, Künstner A, Fähnrich A, et al. Low-level mitochondrial heteroplasmy modulates DNA replication, glucose metabolism and lifespan in mice. Sci Rep. (2018) 8:5872. doi: 10.1038/s41598-018-24290-6
38. Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. (2009) 25:1754–60. doi: 10.1093/bioinformatics/btp324
39. Picard Tools - By Broad Institute. Available online at: http://broadinstitute.github.io/picard/ (accessed July 18, 2016).
40. McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, et al. The genome analysis toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. (2010) 20:1297–303. doi: 10.1101/gr.107524.110
41. alimanfoo/pysamstats. GitHub. Available online at: https://github.com/alimanfoo/pysamstats (accessed July 18, 2016).
42. Robinson JT, Thorvaldsdóttir H, Winckler W, Guttman M, Lander ES, Getz G, et al. Integrative genomics viewer. Nat Biotech. (2011) 29:24–6. doi: 10.1038/nbt.1754
43. Shen L, Attimonelli M, Bai R, Lott MT, Wallace DC, Falk MJ, et al. MSeqDR mvTool: a mitochondrial DNA Web and API resource for comprehensive variant annotation, universal nomenclature collation, and reference genome conversion. Hum Mutat. (2018) 39:806–10. doi: 10.1002/humu.23422
44. Weissensteiner H, Pacher D, Kloss-Brandstätter A, Forer L, Specht G, Bandelt H-J, et al. HaploGrep 2: mitochondrial haplogroup classification in the era of high-throughput sequencing. Nucleic Acids Res. (2016) 44:W58–63. doi: 10.1093/nar/gkw233
45. van Oven M, Kayser M. Updated comprehensive phylogenetic tree of global human mitochondrial DNA variation. Hum Mutat. (2009) 30:E386–94. doi: 10.1002/humu.20921
46. Okonechnikov K, Golosova O, Fursov M. Unipro UGENE: a unified bioinformatics toolkit. Bioinformatics. (2012) 28:1166–7. doi: 10.1093/bioinformatics/bts091
47. Fay MP. Confidence intervals that match Fisher's exact or Blaker's exact tests. Biostatistics. (2010) 11:373–4. doi: 10.1093/biostatistics/kxp050
48. Yusuf S, Peto R, Lewis J, Collins R, Sleight P. Beta blockade during and after myocardial infarction: an overview of the randomized trials. Prog Cardiovasc Dis. (1985) 27:335–71. doi: 10.1016/S0033-0620(85)80003-7
49. Viechtbauer W. Conducting meta-analyses in R with the metafor Package. J Stat Softw. (2010) 36:1–48. doi: 10.18637/jss.v036.i03
50. Hommel G. A stagewise rejective multiple test procedure based on a modified Bonferroni test. Biometrika. (1988) 75:383–6. doi: 10.1093/biomet/75.2.383
51. Kibsgaard L, Rasmussen M, Lamberg A, Deleuran M, Olesen AB, Vestergaard C. Increased frequency of multiple sclerosis among patients with bullous pemphigoid: a population-based cohort study on comorbidities anchored around the diagnosis of bullous pemphigoid. Br J Dermatol. (2017) 176:1486–91. doi: 10.1111/bjd.15405
52. Ji F, Sharpley MS, Derbeneva O, Alves LS, Qian P, Wang Y, et al. Mitochondrial DNA variant associated with Leber hereditary optic neuropathy and high-altitude Tibetans. Proc Natl Acad Sci USA. (2012) 109:7391–6. doi: 10.1073/pnas.1202484109
53. Tranah GJ, Nalls MA, Katzman SM, Yokoyama JS, Lam ET, Zhao Y, et al. Mitochondrial DNA sequence variation associated with dementia and cognitive function in the elderly. J Alzheimers Dis. (2012) 32:357–2. doi: 10.3233/JAD-2012-120466
54. Ross OA, McCormack R, Maxwell LD, Duguid RA, Quinn DJ, Barnett YA, et al. mt4216C variant in linkage with the mtDNA TJ cluster may confer a susceptibility to mitochondrial dysfunction resulting in an increased risk of Parkinson's disease in the Irish. Exp Gerontol. (2003) 38:397–405. doi: 10.1016/S0531-5565(02)00266-8
55. Huerta C, Castro MG, Coto E, Blázquez M, Ribacoba R, Guisasola LM, et al. Mitochondrial DNA polymorphisms and risk of Parkinson's disease in Spanish population. J Neurol Sci. (2005) 236:49–54. doi: 10.1016/j.jns.2005.04.016
56. Nicholls TJ, Minczuk M. In D-loop: 40 years of mitochondrial 7S DNA. Exp Gerontol. (2014) 56:175–81. doi: 10.1016/j.exger.2014.03.027
57. Jenuth JP, Peterson AC, Shoubridge EA. Tissue-specific selection for different mtDNA genotypes in heteroplasmic mice. Nat Genet. (1997) 16:93–5. doi: 10.1038/ng0597-93
58. Schmidt E, Zillikens D. Pemphigoid diseases. Lancet. (2013) 381:320–32. doi: 10.1016/S0140-6736(12)61140-4
59. Claesson MJ, Cusack S, O'Sullivan O, Greene-Diniz R, de Weerd H, Flannery E, et al. Composition, variability, and temporal stability of the intestinal microbiota of the elderly. Proc Natl Acad Sci USA. (2011) 108(Suppl. 1):4586–91. doi: 10.1073/pnas.1000097107
60. O'Toole PW, Jeffery IB. Gut microbiota and aging. Science. (2015) 350:1214–5. doi: 10.1126/science.aac8469
61. Franceschi C, Bonafè M, Valensin S, Olivieri F, De Luca M, Ottaviani E, et al. Inflamm-aging: an evolutionary perspective on immunosenescence. Ann NY Acad Sci. (2000) 908:244–54. doi: 10.1111/j.1749-6632.2000.tb06651.x
62. Franceschi C, Campisi J. Chronic inflammation (inflammaging) and its potential contribution to age-associated diseases. J Gerontol A Biol Sci Med Sci. (2014) 69(Suppl. 1):S4–9. doi: 10.1093/gerona/glu057
63. Belkaid Y, Segre JA. Dialogue between skin microbiota and immunity. Science. (2014) 346:954–9. doi: 10.1126/science.1260144
64. Saint-Georges-Chaumet Y, Edeas M. Microbiota-mitochondria inter-talk: consequence for microbiota-host interaction. Pathog Dis. (2016) 74:ftv096. doi: 10.1093/femspd/ftv096
Keywords: mitochondrial DNA, mitochondrial haplogroup, polymorphisms, autoimmune skin diseases, bullous pemphigoid, mitochondrial function, next generation sequencing
Citation: Russlies J, Fähnrich A, Witte M, Yin J, Benoit S, Gläser R, Günter C, Eming R, Erdmann J, Gola D, Gupta Y, Holtsche MM, Kern JS, König IR, Kiritsi D, Lieb W, Sadik CD, Sárdy M, Schauer F, van Beek N, Weidinger A, Worm M, Zillikens D, Schmidt E, Busch H, Ibrahim SM and Hirose M (2019) Polymorphisms in the Mitochondrial Genome Are Associated With Bullous Pemphigoid in Germans. Front. Immunol. 10:2200. doi: 10.3389/fimmu.2019.02200
Received: 31 May 2019; Accepted: 30 August 2019;
Published: 22 November 2019.
Edited by:
Attila Mócsai, Semmelweis University, HungaryReviewed by:
Takashi Hashimoto, Osaka City University, JapanWataru Nishie, Hokkaido University, Japan
Teruki Dainichi, Kyoto University, Japan
Copyright © 2019 Russlies, Fähnrich, Witte, Yin, Benoit, Gläser, Günter, Eming, Erdmann, Gola, Gupta, Holtsche, Kern, König, Kiritsi, Lieb, Sadik, Sárdy, Schauer, van Beek, Weidinger, Worm, Zillikens, Schmidt, Busch, Ibrahim and Hirose. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Misa Hirose, bWlzYS5oaXJvc2UmI3gwMDA0MDt1a3NoLmRl
†ORCID: Regine Gläser orcid.org/0000-0002-5819-777X