- Research
- Open access
- Published:
Whole mitogenome sequencing uncovers a relation between mitochondrial heteroplasmy and leprosy severity
Human Genomics volume 17, Article number: 110 (2023)
Abstract
Background
In recent years, the mitochondria/immune system interaction has been proposed, so that variants of mitochondrial genome and levels of heteroplasmy might deregulate important metabolic processes in fighting infections, such as leprosy.
Methods
We sequenced the whole mitochondrial genome to investigate variants and heteroplasmy levels, considering patients with different clinical forms of leprosy and household contacts. After sequencing, a specific pipeline was used for preparation and bioinformatics analysis to select heteroplasmic variants.
Results
We found 116 variants in at least two of the subtypes of the case group (Borderline Tuberculoid, Borderline Lepromatous, Lepromatous), suggesting a possible clinical significance to these variants. Notably, 15 variants were exclusively found in these three clinical forms, of which five variants stand out for being missense (m.3791T > C in MT-ND1, m.5317C > A in MT-ND2, m.8545G > A in MT-ATP8, m.9044T > C in MT-ATP6 and m.15837T > C in MT-CYB). In addition, we found 26 variants shared only by leprosy poles, of which two are characterized as missense (m.4248T > C in MT-ND1 and m.8027G > A in MT-CO2).
Conclusion
We found a significant number of variants and heteroplasmy levels in the leprosy patients from our cohort, as well as six genes that may influence leprosy susceptibility, suggesting for the first time that the mitogenome might be involved with the leprosy process, distinction of clinical forms and severity. Thus, future studies are needed to help understand the genetic consequences of these variants.
Background
Mitochondria are cytoplasmic organelles that participate in several processes in cellular functioning in humans, including different types of cell death, control of calcium levels, regulation of the immune system, metabolic cell signaling and generation of cellular energy in the form of Adenosine Triphosphate (ATP) by tricarboxylic acid (TCA) cycle and oxidative phosphorylation (OXPHOS) [1,2,3].
Given their evolutionary origin, mitochondria have their own genetic material (mtDNA), double-stranded circular molecules located in the mitochondrial matrix and associated with the inner membrane of the organelle [4, 5]. The human mitogenome is 16,569 bp in length, with 37 genes—13 of OXPHOS-associated polypeptides, 22 of transfer RNA (tRNA) and two of ribosomal RNA (rRNA)—in addition to non-coding regions, which include the displacement loop (D-loop) region [6,7,8,9].
Mitochondria play an extremely important part in the immune system, such as proliferation in the energy supply for the synthesis of signaling and effector molecules, as well as acting directly on signaling pathways for the activation of these cells through intermediary molecules [10]. These include mtDNA, which might act in the pathogenesis process by Mycobacterium genus [11], and mitochondrial reactive oxygen species (mtROS), which play a central role in the process of NLRP3 inflammasome regulation and activity. This molecular complex is crucial in the process of releasing proinflammatory oxytocins, such as IL-1b and IL-18 [12, 13].
Importantly, different cells have a variable number of mtDNA copies that can lead to a state called mitochondrial heteroplasmy, which can be characterized as the presence of two or more mtDNA variants in varying proportions within individual organisms [8, 14, 15]. Heteroplasmy is a normal part of healthy human biology, but it is also relevant in disease processes, with the level of heteroplasmy being crucial for the expression of specific pathological phenotypes [2, 8, 15]. In addition, given the importance of mitochondria to the immune system, the accumulation of mutations can lead to mitochondrial dysfunction, which in turn might be responsible for cellular dysregulation, leading to the development and aggravation of multiple infectious diseases [16, 17] such as leprosy.
Leprosy, or Hansen’s disease, is a chronic granulomatous bacterial infection that primarily affects skin and peripheral nerves. The etiological agent is the obligate intracellular bacteria Mycobacterium leprae, which produces a broad spectrum of the illness, while diffuse lepromatous leprosy may also be caused by Mycobacterium lepromatosis, a new species described in 2008 [18]. Regardless, host factors that regulate susceptibility to the diverse clinical forms of the disease are important, but largely unknown [19, 20]. Leprosy remains a serious public health problem in various parts of the world and, in 2021, more than 140,000 new cases were reported globally, a higher number compared to 127,396 cases reported in 2020, but still a lower number compared to 202,488 cases reported in 2019. However, these data should be observed with great caution because this decrease is probably due to less detection during the COVID-19 pandemic [21, 22].
Leprosy classification is complex and may include clinical, histopathological, microbiological and immunological features. The Ridley–Jopling system classifies leprosy as a spectral disease: in one extreme, there is the polar tuberculoid form (TT), with a low bacterial load, mainly cell-mediated immunity and minor production of specific antibodies. In the other extreme, there is the polar lepromatous form (LL), in which patients show high bacterial load and respond to infection with elevated production of antibodies, as well as lower or absent M. leprae-specific cell-mediated immunity. Between these two polar forms, there is the clinically unstable borderline spectrum: borderline-tuberculoid (BT), borderline-borderline (BB) and borderline-lepromatous (BL), with BB being the most unstable form [23].
The complete genome sequence of the M. leprae contains 3,268,210 bp and has an average G + C content of 57.8% [24]. By being able to absorb host cell carbon, only about half of the bacillus genome contains functional protein-coding genes [25]. For this reason, M. leprae has a dependence on the host’s energy production and nutritional products, resulting in parasitic life adaptation, undoubtedly involving the main function of the mitochondrion, due to the cell signaling pathways in which this organelle participates and connects its metabolism to meet their nutrient demands [26, 27].
Notably, mitochondria have important functions in the regulation of novel immune signaling pathways exerting control over redox metabolism, energy flow, apoptosis, xenophagy and activating inflammasomes. For instance, it has already been shown that leprosy patients have a differential expression of non-coding RNAs such as piRNAs compared to clinically healthy people in the clinical form, as well as in the clinical spectrum of the disease. One of the differential expressions is related to the activation of anti-apoptotic pathways, evidence of the pathogen's interference in the host's mitochondria [28]. In addition, these organelles directly influence intracellular pathogens that attempt to invade their space; inhibition of mitochondrial energy metabolism likely emerges as a novel and overlooked mechanism developed by M. leprae to evade xenophagy and the host immune response [25].
In this scenario, mutations that affect mitochondrial functions might influence the host response to this infection, leading to multiple possibilities in leprosy development and outcome. These mutations might be present in the host's mitochondrial genome. Therefore, we sequenced the whole mitochondrial genome to investigate variants and their heteroplasmy levels in the context of leprosy. To the best of our knowledge, this is the first study to perform such genomic approach regarding infection by M. leprae.
Methods
Sampling
Blood samples were obtained from patients affected by leprosy (n = 33, case group) and healthy household contacts with leprosy patients (n = 37, control group), all residents of Pará state, Brazil. The case group was composed of borderline lepromatous (BL) leprosy (n = 12), lepromatous (LL) (n = 11) and borderline tuberculoid (BT) leprosy (n = 10). This study adhered to the Declaration of Helsinki and was approved by the Ethics Committee of Institute of Health Sciences at the Federal University of Pará (CEP-ICS/UFPA n. 197/07), and all participants read and signed an informed consent form.
As inclusion and exclusion criteria, the samples from case group were selected from patients affected by leprosy who had a positive clinical and laboratory diagnosis. The samples from control group were selected from healthy household contacts who had a negative clinical and laboratory diagnosis. All participants were recruited at Dr. Marcello Candia Reference Unit in Sanitary Dermatology of the State of Pará (URE) located in Marituba, Pará, Brazil.
Clinical and laboratory diagnosis
The diagnosis of leprosy was conducted with the well-accepted clinical signs and symptoms based on the Ridley–Jopling classification, including detection of hypopigmented or infiltrated skin lesions with loss of sensation assessed with standard graded Semmes–Weinstein monofilaments and the palpation of peripheral nerves to identify characteristic pain associated with inflammation or swelling, as previously described [29, 30].
To establish laboratory parameters, antibody titers of three antigens (NDO-BSA, LID-1 and NDO-LID) were evaluated with molecular detection of RLEP by qPCR in leprosy patients and contacting patients [29]. The cutoff values of antibody titers were determined using previously described criteria, and the cutoffs for anti-NDO-BSA and anti-LID-1 were 0.295 and, for anti-NDO-LID, the cutoff was 0.475 [29].
To determine qPCR positivity, a standard curve was prepared from purified M. leprae DNA extracted from nude mouse footpads, and then five standard dilution points were included in each plate, considering the samples as positive when the fluorescent signal crossed the automatically calculated threshold line [29]. Amplifications with cycle threshold (Ct) ≤ 45 were considered positive for RLEP [29].
DNA extraction
DNA was extracted by phenol–chloroform method [31]. Quantification of the extracted material was performed with the NanoDrop 1000 spectrophotometer (Thermo Fisher Scientific, Wilmington, DE, USA).
Amplification and sequencing
Amplification of mtDNA from the total DNA was performed by conventional PCR with specific primers, as described by Cavalcante et al. [1], to cover the entire mitochondrial genome. To verify the quality of the amplification, the samples were applied to a 1% agarose gel and, later, measured in a Qubit 2.0 fluorometer for the library preparation (Thermo Fisher Scientific). Sequencing of the complete mitochondrial genome was performed using Nextera XT DNA Library Preparation Kit (Illumina Inc., Chicago, IL, USA) to prepare the libraries and MiSeq Reagent Kit V3 (600-cycles) (Illumina) for sequencing on the MiSeq System (Illumina), according to the manufacturer’s instructions. During the preparation of the libraries, High Sensitivity D1000 ScreenTape was used on the Agilent 2200 TapeStation System (Agilent Technologies, Santa Clara, CA, USA) to assess the quality of the genetic material.
Bioinformatics and statistical analyses
After sequencing, we updated the pipeline for bioinformatics analysis previously described [1]. The paired-end sequencing reads (.fastq files) were trimmed with Trimmomatic v.0.39 [32] to remove leading low quality (LEADING:10), trailing low quality (TRAILING:10) and to scan reads with a 3-base wide sliding window, cutting when the average quality per base drops below 10 (SLIDINGWINDOW:3:10) and those reads with length less than 36nt were discarded. After trimming, fastq files were aligned with the human reference mtDNA sequence—revised Cambridge reference sequence (rCRS)—using Burrows-Wheeler Alignment tool (BWA, v.0.7) [33]. SAMtools (v.1.15.1) [34] were used for mapping and sorting sequences, while Picard was used to mark the duplicated reads (v.2.27.5, available at https://github.com/broadinstitute/picard).
After preprocessing the sequences in the aforementioned steps, paired-end.bam files were submitted to mutserve for SNP calling, SNP annotation and heteroplasmy detection (https://mitoverse.readthedocs.io/mutserve/mutserve/). For SNP calling, we performed mutserve for each sample with the following quality parameters: mapping quality = 20, base quality = 20 and alignment quality = 30. SNP annotation was based on the rCRS genome annotation (available at https://github.com/seppinho/mutserve). Mutserve outputs SNPs in.vcf file format that was used as input for inferences of mitochondrial haplogroups by HaploGrep (v2) [35]. To reinforce the reliability of the results in the analyses of the variants, only those that passed rigorous additional filters were considered (filter = PASS; coverage ≥ 545 per variant; presence of heteroplasmy). Then, three databases were used for identification of the found variants based on their position: dbSNP (https://www.ncbi.nlm.nih.gov/snp/) [36], ClinVar (https://www.ncbi.nlm.nih.gov/clinvar/) [37] and gnomAD Browser (https://gnomad.broadinstitute.org/) [38], in addition to literature search. R language [39] was employed for statistical analysis and graph generation with the following packages: ggplot2 [40] and UpSetR [41]. P-value < 0.05 was considered as statistically significant.
Results
Characterization of the cohort
After processing, the mean sequencing coverage of all samples was 1489×, and two low-quality samples were excluded from the analyses. In the studied cohort, we observed similar average age for both case and control groups, as well as for the case subgroups (Table 1). Interestingly, there were more females than males in the LL subgroup, in comparison with the other two analyzed clinical forms (BL and BT), differing from a previous study by our research group that reported more males with this clinical form [42]. However, in the case group considering all clinical forms, we observed more males than females, corroborating a meta-analysis that investigated leprosy in multiple populations, including Brazil [43].
By analyzing the distribution of mitochondrial macro-haplogroups in our cohort, we found that H2, of European (EUR) ancestry, is the most frequent individual haplogroup in both case (32.35%) and control (21.62%), but the Native American (NAT) haplogroups together account for the largest share in both groups (46.88% in case and 43.24% in control) (Table 2). It is also noteworthy that NAT and African (AFR) ancestries presented more diversity of haplogroups than the others.
Regarding the leprosy subtypes present in our case group (BT, BL and LL), we found three ancestries based on mitochondrial haplogroups, with different distributions among the subtypes (Fig. 1A). Notably, the Native American ancestry was as frequent as the European ancestry in LL and different in BT, while it was four times more frequent than the other ancestries in BL. Then, an analysis unifying the subtypes into poles was carried out, considering tuberculoid pole (T, with BT, n = 10) and lepromatous pole (L, with BL + LL, n = 23), to assess the mitochondrial ancestry of these individuals considering both poles (Fig. 1B). NAT remained four times more frequent in L than in T, and nearly two times more frequent than EUR and AFR within L pole; in T pole, EUR was nearly two times more frequent than the other two ancestries.
Distribution of mitochondrial variants
Furthermore, we investigated the distribution of the general mitochondrial heteroplasmic variants present in the study cohort, to assess the variants that were exclusive or shared among the different groups and subgroups, that is, the intersection between groups, as seen in Fig. 2. We found 116 variants to be present in at least two of the subtypes of the case group, but not in the control group. In addition, 1332 variants were found exclusively in one leprosy subtype.
Particularly, 15 mitochondrial variants were exclusively found in all three case subgroups, suggesting they might play a role in the leprosy process: m.3791T > C (MT-ND1, no dbSNP identification found), m.5317C > A (MT-ND2, no dbSNP identification found), m.8455C > T (MT-ATP8, rs1603221490), m.8503T > C (MT-ATP8, rs1556423476), m.8545G > A (MT-ATP8, rs1603221578), m.9044T > C (MT-ATP6, no dbSNP identification found), m.9103T > C (MT-ATP6, rs1603222077), m.12879T > C (MT-ND5, rs1556424182), m.13512A > G (MT-ND5, rs878930809), m.14721G > A (MT-TE, rs1603224843), m.14860C > T (MT-CYB, no dbSNP identification found), m.14905G > A (MT-CYB, rs193302983), m.14941A > G (MT-CYB, rs1603224969), m.15019T > C (MT-CYB, no dbSNP identification found) and m.15837T > C (MT-CYB, no dbSNP identification found). Thus, four are in Complex I genes, five in Complex III, five in Complex V and one in tRNA. Out of these, some variants in the protein-coding genes stand out for being missense and presenting a relevant pathogenicity prediction: m.3791T > C in MT-ND1 (probably damaging), m.5317C > A in MT-ND2 (possibly damaging), m.8545G > A in MT-ATP8 (benign), m.9044T > C in MT-ATP6 (probably damaging) and m.15837T > C in MT-CYB (benign). It is noteworthy that, although most of the variants found in the cohort were detected at low levels of heteroplasmy, recent single-cell analyses have shown that even low-level heteroplasmy (< 5%) can alter transcription levels of nuclear genes involved in ATP synthesis and important cellular processes [44].
Regardless, to analyze the distribution of heteroplasmic mutations in leprosy poles and the intersection between groups, and to exclude variants with low levels of heteroplasmy, we filtered the same variants for those with heteroplasmy levels > 5% and < 95%, to control possible artifacts and false positives (Fig. 3). After this filter, it is possible to notice that 26 variants stand out as shared by both T and L poles.
Of these 26 variants, it is noteworthy that most of them were concentrated in rRNA genes (MT- RNR2) and genes encoding OXPHOS complexes, particularly in complexes I (MT-ND1, MT-ND5), III (MT-CYB) and IV (MT- CO2 and MT- CO3), illustrated in Table 3. This distribution could suggest that at the poles of leprosy there may be a dysregulation of mitochondrial gene expression, as well as in the oxidative environment and energy production, and that these genes may develop roles in disease processes.
Analysis of heteroplasmy levels
To investigate mitochondrial heteroplasmy in our cohort, we assessed the levels in which this state was presented in each region of the mitogenome, according to the filtering of > 5% and < 95% (Fig. 4). Interestingly, it seems that tRNAs and control regions (CR) display higher heteroplasmy levels in healthy controls than in leprosy patients, while protein-coding genes have a prominent heteroplasmy variation in the BT subgroup with a widespread pattern in comparison with BL and LL in all OXPHOS protein-coding genes (MT-ND1, MT-ND2, MT-ND3, MT-ND4, MT-ND4L, MT-ND5, MT-ND6, MT-CYB, MT-CO1, MT-CO2, MT-CO3, MT-ATP6 and MT-ATP8).
Moreover, Fig. 5 shows that T pole and L pole present a different pattern of heteroplasmy levels, with L pole being more like healthy controls and T pole being more distant than these states. D-loop presents increased levels of heteroplasmy in T pole when compared to L pole and HC, suggesting that the regulation of mitochondria might be altered in this unstable state.
Strikingly, MT-DLOOP1, MT-ND1 and MT-CYB have a much higher mean of heteroplasmy levels in the control group in comparison with the leprosy clinical forms; not as high, but the heteroplasmy mean levels of MT-CO1 and MT-TG are also elevated when compared to the other groups. In T pole, it should be noted that the MT- DLOOP2, MT-TC, MT-ND3, MT-ND4L and MT-TS2 also present a higher mean level than the other analyzed groups. Considering that T pole, represented by the BT subtype, is clinically more unstable than the L pole, it is not surprising to see multiple regions with increased heteroplasmy levels in this subgroup, indicating an active inflammatory process with different immune responses [45, 46].
Furthermore, by analyzing the overall presence of variants and their respective heteroplasmy levels in case and control groups, we found statistical significance in seven regions: one control region, one rRNA and five genes that encode Complexes I, III and IV (Fig. 6). In most of these scenarios, there were more heteroplasmic variants in leprosy patients than in healthy controls, reinforcing the idea that these variants may contribute to the process of susceptibility to leprosy.
When we considered only the category of mitochondrial regions (coding genes, CR, rRNA and tRNA), we found that all categories still presented statistical significance regarding the presence of heteroplasmic variants (Fig. 7). Notably, there are more variants in the case than in the control, and these variants have diverse heteroplasmy levels, particularly in the rRNA genes. This might reflect the relevant presence of heteroplasmic variants in Complexes I, III, IV, as well as in 16S rRNA, as shown in Fig. 6.
Discussion
Although leprosy is an important public health problem, it is still neglected and overlooked in research studies, in particular the mechanisms by which the host mitogenome can influence the profile of susceptibility to the disease, especially in the northern region of Brazil [47].
Regarding the results of the analysis of the ancestry profile of the individuals in the cohort, they were expected given the formation process of the Brazilian population and the fact that our cohort is from the North region of the country, in which the Native American ancestry is particularly frequent [48, 49]. Furthermore, the results in Fig. 1 suggest that Native American ancestry could have an influence on the development of different types of leprosy subtypes, particularly BL, upon M. leprae infection and reinforced the analysis that mitochondrial ancestry might influence the developed pole.
As for the general distribution of variants illustrated in Fig. 2, no previous studies were found in the global literature with most of these. However, m.15837T > C has been identified in breast nipple aspirate fluid in breast cancer [50]. In addition, we only found a few variants with dbSNP identification, suggesting the remaining variants to be unreported in different databases. Considering three of these variants have been predicted to have a damaging potential, we recommend more studies to clarify their impact in OXPHOS (Complexes I and V) that might affect disease processes such as leprosy.
It is important to emphasize that these exclusive mutations are divided into five differential groups: MT-ND (MT-ND1, MT-ND2 and MT-ND5), MT-ATP6, MT-ATP8, MT-CYB and MT-TE genes (Fig. 8). The MT-ND genes codify NADH dehydrogenase proteins, composing Complex I. This protein complex transfers the energy generated from the oxidation of NADH to NAD+ and the transfer of protons from the mitochondrial matrix to the intermembrane space through flavin mononucleotide (FMN) and seven to nine iron-sulfur (Fe-S) clusters for the reduction of ubiquinone, the first electron acceptor [51,52,53,54].
The NAD + generated by Complex I is used in numerous metabolic reactions by NADH-linked dehydrogenases, including components of the fatty acid oxidation cycle. NADPH is a crucial reductant used in lipid anabolism including synthesis of important components of the mycobacterial cell wall [54] and represents a major electron donor feeding the respiratory chain. Therefore, the NADH molecules generated by 3b-HSD activity could supply, at least in part, electrons to the respiratory chain contributing to M. leprae ATP synthesis [54].
The 3b-HSD activity generates the electron donors NADH and NADPH that, respectively, fuel the M. leprae respiratory chain and provide reductive power for the biosynthesis of the dominant bacterial cell wall lipids and phenolic glycolipid (PGL)-I [55], the M. leprae-specific antigen first reported in 1980 [56], initially tested as a tool for leprosy serodiagnosis [57, 58], but also helping the diagnosis and prediction of relapses [59].
Unfortunately, PGL-I shows low ability to detect true positives, so it is now a well-known biomarker of M. leprae exposition and has been used to recognize hidden leprosy cases [60]. Its association with molecular tools or spatial epidemiology is a strong tool to give direction to public policies that lead to an increase in the detection of cases among students and contacts of patients and that, if timely diagnosed, prevent the physical disability of patients, while also breaking the chain of transmission.
Mutations in the MT-ATP8 gene (a component of Complex V) have been described in reactive oxygen species (ROS) generation [61]. The relationship between M. leprae and the increased damage caused by ROS remains unclear, but previous studies have suggested a correlation between bacillary load and oxidative stress [62]. Under chronic inflammatory conditions, ROS reduce activation signals to the T cell and impair the immune response against the pathogen [63]. However, the delicately maintained physiological balance of oxidative stress is changed in favor of ROS from phagocyte and by the treatment (multidrug therapy—MDT), both increasing production of ROS [64, 65] and decreasing antioxidants [63, 66]. Recently, a study verified that elderly patients with leprosy have higher ROS than younger patients [67].
The MT-ATP6 gene encodes the ATP6 subunit of mitochondrial ATP synthase (Complex V), which catalyzes the last stage of the electron transport chain (ETC), a key molecular process for the normal axonal function of the central and peripheral nervous system [68, 69]. Mutations in MT-ATP6 have been linked to episodic muscle weakness [70] and Charcot-Marie-Tooth axonal disease [71]. This disorder damages peripheral nerves and affects both sensory and motor nerves (nerves that trigger an impulse from contraction in muscles) in the arms, hands, legs and feet. Peripheral neuropathies form an integral part of the symptomatology of leprosy and include numbness, painlessness, anesthesia, hypoesthesia and patchy motor deficits, paresthesias (pins and needles), pain (allodynia and dysesthesias), impairment of temperature perception followed by touch and pain; sensory loss, wasting and weakness gradually occur in involved nerve territories and partial involvement [72].
The MT-CYB gene codifies a subunit of cytochrome b oxidase (Complex III), involved in oxidative phosphorylation [73], and it has been described as less expressed in leprosy patients when compared to non-leprosy population [74]. The M. leprae, like M. ulcerans, lost genes encoding cytochrome b oxidase [24] during reductive evolution, making M. leprae extremely sensitive to QcrB inhibitors, such as the drug Telacebec [75]. In this context, it should be noted that the established MDT for leprosy was efficient in reducing most cases worldwide, but after 50 years of the same treatment regimen, the increase in resistant cases is the most critical problem [76]. Hence, QcrB inhibitors could represent a new class of bactericidal drugs for leprosy [77], due to their high potency against M. ulcerans [78].
The MT-TE gene belongs to the set of tRNAs encoded by mtDNA, being important for the biosynthesis of mitochondrial proteins and one of the major causes of disorders in the genome [6, 79, 80]. Pathogenic variants on mitochondrial tRNAs cause a wide range of disease phenotypes, with energy-intensive tissues such as neuromuscular and nervous tissues being particularly vulnerable, with progressive neurological deficits being the most prominent and often the most disabling feature of the disease [79, 81]. Mutations in the MT-TE gene have already been associated with the development of diabetes and myopathies, as well as early-onset cataracts, ataxia and progressive paraparesis [81, 82], but there are no reports in the literature about the variant found in our study.
Furthermore, when analyzing the distribution of variants after filtering, 26 variants stand out to be shared in both leprosy poles (Fig. 3). Overall, 16 variants presented identification in dbSNP (Table 3), but only three variants have already been described in the ClinVar as being associated with diseases, such as Leigh syndrome (m.8027G > A) and familial breast cancer (m.14783T > C, m.14905G > A). Two variants were characterized as missense, belonging to MT-ND1 (m.4248T > C) and MT-CO2 (m.8027G > A) genes, both being OXPHOS regulators. Interestingly, no previous studies were found in the global literature on most of these variants.
The heteroplasmy analyses suggest that heteroplasmy across the mitochondrial genome can occur differently in each region depending on the affected or unaffected state and the clinical form that leprosy presents (Fig. 4). Therefore, heteroplasmy seems to influence the oxidative environment in the development of this disease. This corroborates previous studies that suggested that the heterogeneity of mtDNA copies might increase during pathophysiological processes and might even be a potential target for therapies of different inflammatory diseases [6, 7, 14]. In leprosy, mitochondrial activity impairment and mtDNA content decrease have been reported, particularly in OXPHOS proteins [25]. This could be related to the observed widespread state of heteroplasmy, especially in the unstable BT clinical form.
The different profile of heteroplasmy between mitochondrial genes at the L, T and HC poles, demonstrated in Fig. 5, suggests that heteroplasmy may influence the type of host immune response to M. leprae. The response to leprosy can be classified into two main types: type 1 reaction (T1R) or reverse reaction that occurs mainly in the unstable forms of the disease (BT, BB, BL) and in the TT form and the type 2 reaction (T2R) or erythema leprosy nodosis (ENL) that occurs mainly in BL and LL patients with high bacillary burdens [83, 84].
The different clinical manifestations of leprosy and the histopathology of the two polar forms of leprosy are also determined by the adaptive immune response, the classic paradigm being that the L pole is associated with a Th2 immune response and the T pole with a Th1 immune response [83, 84]. Our findings demonstrate that the heteroplasmic profile of the HC group is more like individuals belonging to the L pole, which suggests that this group tends to be associated with a Th2-type humoral response.
Recent studies have demonstrated that intact mitochondria from human plasma cells harbor immunologically active membrane-associated proteins, including CD270 and programmed cell death ligand 1 (PD-L1) [85]. In addition, intact human plasma mitochondria and immunologically active surface proteins have been associated with upregulation of activated CD4 + T cells and CD8 + T cells and reduced concentrations of pro-inflammatory cytokines [85].
In this perspective, these results suggest that the rate of differential heteroplasmic mutations between genes may be influencing the dysregulation of mitochondrial functions, which may induce a distinct immune response among individuals, dividing them in the clinical forms of leprosy. Importantly, it should be noted that, for many infectious diseases, host genetic factors have long been considered a major contributor to variations in individual susceptibility and immune response [84].
Conclusion
Here, we reported the whole mitochondrial genome sequencing of leprosy patients and healthy unrelated household controls from a Brazilian population. Importantly, this is likely to be the first study to report a clear association between mitochondrial heteroplasmy and leprosy. This unprecedented approach showed that a higher number of mitochondrial variants and diverse heteroplasmy levels are significantly found in the leprosy patients from our cohort, suggesting for the first time that the mitochondrial genome, and particularly heteroplasmy, may be involved in the inflammatory response observed in the peripheral nerve trunks and in the skin of leprosy cases, as well as in defining the immunological response that determines the clinical polarization of leprosy.
Furthermore, we identified 26 heteroplasmic variants shared between the T and L poles that are present in the MT-RNR2, MT-ND1, MT-ND5, MT-CYB, MT-CO2 and MT-CO3 genes, suggesting that these genes may be correlated with the susceptibility and severity of leprosy. The infectious process of leprosy has a direct impact on mitochondrial functions, and their dysregulation is objectively involved in the pathophysiological process of the disease, in the viability of the pathogen in the organism, in the host's homeostasis to the immune response and in the therapeutic response. Thus, further clinical studies with larger cohorts with sequencing of other tissues and/or functional studies regarding the found variants must be carried out to expand knowledge about the pathogenesis of M. leprae and its possible impacts on the mitochondrial genetic profile of affected individuals and, thus, suggest potential biomarkers for leprosy.
Availability of data and materials
All raw sequences are deposited at the European Nucleotide Archive (ENA) under the accession number PRJEB59275.
References
Cavalcante GC, Marinho ANR, Anaissi AK, Vinasco-Sandoval T, Ribeiro-dos-Santos A, Vidal AF, et al. Whole mitochondrial genome sequencing highlights mitochondrial impact in gastric cancer. Sci Rep. 2019;9:15716.
Andrieux P, Chevillard C, Cunha-Neto E, Nunes JPS. Mitochondria as a cellular hub in infection and inflammation. IJMS. 2021;22:11338.
Yan C, Duanmu X, Zeng L, Liu B, Song Z. Mitochondrial DNA: distribution, mutations, and elimination. Cells. 2019;8:379.
Chapman J, Ng YS, Nicholls TJ. The maintenance of mitochondrial DNA integrity and dynamics by mitochondrial membranes. Life. 2020;10:164.
Roger AJ, Muñoz-Gómez SA, Kamikawa R. The origin and diversification of mitochondria. Curr Biol. 2017;27:R1177–92.
Nguyen NNY, Kim SS, Jo YH. Deregulated mitochondrial DNA in diseases. DNA Cell Biol. 2020;39:1385–400.
Nissanka N, Moraes CT. Mitochondrial DNA heteroplasmy in disease and targeted nuclease-based therapeutic approaches. EMBO Rep. 2020;21:e49612.
Pérez-Amado CJ, Bazan-Cordoba A, Hidalgo-Miranda A, Jiménez-Morales S. Mitochondrial heteroplasmy shifting as a potential biomarker of cancer progression. IJMS. 2021;22:7369.
van den Ameele J, Li AYZ, Ma H, Chinnery PF. Mitochondrial heteroplasmy beyond the oocyte bottleneck. Semin Cell Dev Biol. 2020;97:156–66.
Wang Y, McLean AS. The role of mitochondria in the immune response in critical illness. Crit Care. 2022;26:80.
Patrick KL, Watson RO. Mitochondria: powering the innate immune response to mycobacterium tuberculosis infection. Infect Immun. 2021;89:e00687-20.
Chen M, Zhu X, Ran L, Lang H, Yi L, Mi M. Trimethylamine-N-Oxide induces vascular inflammation by activating the NLRP3 inflammasome through the SIRT3-SOD2-mtROS signaling pathway. JAHA. 2017;6:e006347.
Swanson KV, Deng M, Ting JPY. The NLRP3 inflammasome: molecular activation and regulation to therapeutics. Nat Rev Immunol. 2019;19:477–89.
Stefano GB, Kream RM. Mitochondrial DNA heteroplasmy in human health and disease. Biomed Rep. 2016;4:259–62.
Wilton PR, Zaidi A, Makova K, Nielsen R. A population phylogenetic view of mitochondrial heteroplasmy. Genetics. 2018;208:1261–74.
Prates Mori M, de Souza-Pinto NC. Role of mitochondrial dysfunction in the pathophysiology of DNA repair disorders: mitochondrial role in DNA repair disorders. Cell Biol Int. 2018;42:643–50.
Roca-Bayerri C, Robertson F, Pyle A, Hudson G, Payne BAI. Mitochondrial DNA damage and brain aging in human immunodeficiency virus. Clin Infect Dis. 2021;73:e466–73.
Han XY, Seo Y-H, Sizer KC, Schoberle T, May GS, Spencer JS, et al. A new Mycobacterium species causing diffuse lepromatous leprosy. Am J Clin Pathol. 2008;130:856–64.
Eichelmann K, González González SE, Salas-Alanis JC, Ocampo-Candiani J. Leprosy. An update: Definition, pathogenesis, classification, diagnosis, and treatment. Actas Dermo-Sifiliográficas. 2013;104:554–63.
Misch EA, Berrington WR, Vary JC, Hawn TR. Leprosy and the human genome. Microbiol Mol Biol Rev. 2010;74:589–620.
Cáceres-Durán MÁ. Epidemiological behavior of leprosy in several Latin American countries, 2011–2020 Comportamento epidemiológico da hanseníase em vários países da América Latina, 2011–2020. Rev Panam Salud Publica. 2022;46:e14.
Weekly Epidemiological Record, 2021, vol. 96, 36 [full issue] [Internet]. World Health Organization = Organisation mondiale de la Santé; 2021; p. 421–44. https://apps.who.int/iris/handle/10665/345048
Ridley DS, Jopling WH. Classification of leprosy according to immunity. A five-group system. Int J Lepr Other Mycobact Dis. 1966;34:255–73.
Cole ST, Eiglmeier K, Parkhill J, James KD, Thomson NR, Wheeler PR, et al. Massive gene decay in the leprosy bacillus. Nature. 2001;409:1007–11.
Oliveira MF, Medeiros RCA, Mietto BS, Calvo TL, Mendonça APM, Rosa TLSA, et al. Reduction of host cell mitochondrial activity as Mycobacterium leprae’s strategy to evade host innate immunity. Immunol Rev. 2021;301:193–208.
de Souza FG, Cavalcante GC. Mitochondria in Mycobacterium infection: From the immune system to mitochondrial haplogroups. IJMS. 2022;23:9511.
Wang D, Su L-Y, Zhang A-M, Li Y-Y, Li X-A, Chen L-L, et al. Mitochondrial DNA copy number, but not haplogroup, confers a genetic susceptibility to Leprosy in Han Chinese from Southwest China. PLOS ONE. 2012;7:e38848.
Pinto P, da Silva MB, Moreira FC, Bouth RC, Gobbo AR, Sandoval TV, et al. Leprosy piRnome: exploring new possibilities for an old disease. Sci Rep. 2020;10:12648.
Gobbo AR, Bouth RC, Moraes TMP, Pinto P, da Costa PF, Barreto JG, et al. NDO-BSA, LID-1, and NDO-LID antibody responses for infection and RLEP by quantitative PCR as a confirmatory test for early Leprosy diagnosis. Front Trop Dis. 2022;3:850886.
Villarroel MF, Orsini MBP, Lima RC, Antunes CMF. Comparative study of the cutaneous sensation of leprosy-suspected lesions using Semmes–Weinstein monofilaments and quantitative thermal testing. Lepr Rev. 2007;78:102–9.
Sambrook J. Molecular cloning: a laboratory manual (2nd edn). (No Title). 1989
Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30:2114–20.
Li H, Durbin R. Fast and accurate long-read alignment with Burrows–Wheeler transform. Bioinformatics. 2010;26:589–95.
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The sequence alignment/map format and SAMtools. Bioinformatics. 2009;25:2078–9.
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.
Sherry ST, Ward M, Sirotkin K. dbSNP-database for single nucleotide polymorphisms and other classes of minor genetic variation. Genome Res. 1999;9:677–9.
Landrum MJ, Lee JM, Benson M, Brown GR, Chao C, Chitipiralla S, et al. ClinVar: improving access to variant interpretations and supporting evidence. Nucleic Acids Res. 2018;46:D1062–7.
Karczewski KJ, Francioli LC, Tiao G, Cummings BB, Alföldi J, Wang Q, et al. The mutational constraint spectrum quantified from variation in 141,456 humans. Nature. 2020;581:434–43.
Core Team R. R: A language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2014.
Wickham H, Chang W, Wickham MH. Package ‘ggplot2.’ Create elegant data visualisations using the grammar of graphics version. 2016;2:1–189.
Lex A, Gehlenborg N, Strobelt H, Vuillemot R, Pfister H. UpSet: visualization of intersecting sets. IEEE Trans Visual Comput Graphics. 2014;20:1983–92.
Salgado CG, Pinto P, Bouth RC, Gobbo AR, Messias ACC, Sandoval TV, et al. miRNome expression analysis reveals new players on leprosy immune physiopathology. Front Immunol. 2018;9:463.
de Paula HL, de Souza CDF, Silva SR, Martins-Filho PRS, Barreto JG, Gurgel RQ, et al. Risk factors for physical disability in patients with leprosy: a systematic review and meta-analysis. JAMA Dermatol. 2019;155:1120.
Nie Y, Murley A, Golder Z, Rowe JB, Allinson K, Chinnery PF. Heteroplasmic mitochondrial DNA mutations in frontotemporal lobar degeneration. Acta Neuropathol. 2022;143:687–95.
Hudson G, Carelli V, Spruijt L, Gerards M, Mowbray C, Achilli A, et al. Clinical expression of leber hereditary optic neuropathy is affected by the mitochondrial DNA–Haplogroup background. Am J Hum Genet. 2007;81:228–33.
Wang J, Peng L-Y, You C-P, Li Q-L, Wen M, Liu S-J, et al. Minifish mtDNA has abundance of repeat sequences and inefficient replication in vitro. CMM. 2014;14:1299–307.
Soares AMM, Corrêa RDGCF, dos Santos KCB, Figueiredo IA, Paiva MDFL, de Aquino DMC. Leprosy cases diagnosed by contacts examination in a hyperendemic capital city of northeastern Brazil. Anais Brasileiros de Dermatologia. 2021;96:510–3.
Andrade RB, Amador MAT, Cavalcante GC, Leitão LPC, Fernandes MR, Modesto AAC, et al. Estimating Asian contribution to the Brazilian population: a new application of a validated set of 61 ancestry informative markers. G3 Genes Genomes Genet. 2018;8:3577–82.
Schaan AP, Costa L, Santos D, Modesto A, Amador M, Lopes C, et al. mtDNA structure: the women who formed the Brazilian Northeast. BMC Evol Biol. 2017;17:185.
Zhu W. Mitochondrial DNA mutations in breast cancer tissue and in matched nipple aspirate fluid. Carcinogenesis. 2004;26:145–52.
Bergman O, Ben-Shachar D. Mitochondrial oxidative phosphorylation system (OXPHOS) deficits in schizophrenia: possible interactions with cellular processes. Can J Psychiatry. 2016;61:457–69.
Cavalcante GC, Magalhães L, Ribeiro-dos-Santos Â, Vidal AF. Mitochondrial epigenetics: non-coding RNAs as a novel layer of complexity. IJMS. 2020;21:1838.
Mimaki M, Wang X, McKenzie M, Thorburn DR, Ryan MT. Understanding mitochondrial complex I assembly in health and disease. Biochim Biophys Acta BBA Bioenerg. 2012;1817:851–62.
Minnikin DE, Kremer L, Dover LG, Besra GS. The methyl-branched fortifications of mycobacterium tuberculosis. Chem Biol. 2002;9:545–53.
Rosa TLSA, Marques MAM, DeBoard Z, Hutchins K, Silva CAA, Montague CR, et al. Reductive power generated by Mycobacterium leprae through cholesterol oxidation contributes to lipid and ATP synthesis. Front Cell Infect Microbiol. 2021;11:709972.
Brennan PJ, Barrow WW. Evidence for species-specific lipid antigens in Mycobacterium leprae. Int J Lepr Other Mycobact Dis. 1980;48:382–7.
Levis WR, Meeker HC, Schuller-Levis GB, Gillis TP, Marino LJ, Zabriskie J. Serodiagnosis of leprosy: relationships between antibodies to Mycobacterium leprae phenolic glycolipid I and protein antigens. J Clin Microbiol. 1986;24:917–21.
van Hooij A, Tjon Kon Fat EM, da Silva MB, Carvalho Bouth R, Cunha Messias AC, Gobbo AR, et al. Evaluation of immunodiagnostic tests for Leprosy in Brazil, China and Ethiopia. Sci Rep. 2018;8:17920.
Devides AC, Rosa PS, de Faria Fernandes Belone A, Coelho NMB, Ura S, Silva EA. Can anti-PGL-1 and anti-NDO-LID-1 antibody titers be used to predict the risk of reactions in leprosy patients? Diagn Microbiol Infect Dis. 2018;91:260–5.
da Silva MB, Li W, Bouth RC, Gobbo AR, Messias ACC, Moraes TMP, et al. Latent leprosy infection identified by dual RLEP and anti-PGL-I positivity: implications for new control strategies. PLOS ONE. 2021;16:e0251631.
Weiss H, Wester-Rosenloef L, Koch C, Koch F, Baltrusch S, Tiedge M, et al. The mitochondrial Atp8 mutation induces mitochondrial ROS generation, secretory dysfunction, and β-Cell mass adaptation in conplastic B6-mtFVB mice. Endocrinology. 2012;153:4666–76.
Swathi M, Tagore R. Study of oxidative stress in different forms of leprosy. Indian J Dermatol. 2015;60:321.
Vijayaraghavan R, Paneerselvam C. Erythrocyte antioxidant enzymes in multibacillary leprosy patients. 2011 [cited 2023 Mar 28]; https://imsear.searo.who.int/jspui/handle/123456789/161761
Yu BP. Cellular defenses against damage from reactive oxygen species. Physiol Rev. 1994;74:139–62.
Vijayaraghavan R, Suribabu CS, Oommen PK, Panneerselvam C. Vitamin E reduces reactive oxygen species mediated damage to bio-molecules in leprosy during multi-drug therapy. Curr Trends Biotechnol Pharm. 2009;3:428–39.
Sunita G. Research: role of antioxidant vitamins in immune function in leprosy. Pharm Glob Int J Compr Pharm. 2011;2:1–3.
da Silva PHL, de Castro KKG, Mendes MA, Calvo TL, Leal JMP, Hacker MDAVB, et al. Increased oxidative stress in elderly leprosy patients is related to age but not to bacillary load. PLOS Negl Trop Dis. 2021;15:e0009214.
Dirican E, Savrun ŞT, Aydın İE, Gülbay G, Karaman Ü. Analysis of mitochondrial DNA cytochrome-b (CYB) and ATPase-6 gene mutations in COVID-19 patients. J Med Virol. 2022;94:3138–46.
Pezeshkpour G, Krarup C, Buchthal F, DiMauro S, Bresolin N, McBurney J. Peripheral neuropathy in mitochondrial disease. J Neurol Sci. 1987;77:285–304.
Panosyan FB, Tawil R, Herrmann DN. Episodic weakness and Charcot–marie–tooth disease due to a mitochondrial MT-ATP6 mutation. Muscle Nerve. 2017;55:922–7.
Pitceathly RDS, Murphy SM, Cottenie E, Chalasani A, Sweeney MG, Woodward C, et al. Genetic dysfunction of MT-ATP6 causes axonal Charcot-Marie-Tooth disease. Neurology. 2012;79:1145–54.
Khadilkar SV, Patil SB, Shetty VP. Neuropathies of leprosy. J Neurol Sci. 2021;420:117288.
Tió-Coma M, Kiełbasa SM, van den Eeden SJF, Mei H, Roy JC, Wallinga J, et al. Blood RNA signature RISK4LEP predicts leprosy years before clinical onset. EBioMedicine. 2021;68:103379.
Guerreiro LTA, Robottom-Ferreira AB, Ribeiro-Alves M, Toledo-Pinto TG, Brito TR, Rosa PS, et al. Gene expression profiling specifies chemokine, mitochondrial and lipid metabolism signatures in leprosy. PLOS ONE. 2013;8:e64748.
Lahiri R, Adams LB, Thomas SS, Pethe K. Sensitivity of Mycobacterium leprae to Telacebec. Emerg Infect Dis. 2022;28:749–51.
van Hooij A, Tjon Kon Fat EM, van den Eeden SJF, Wilson L, da Silva MB, Salgado CG, et al. Field-friendly serological tests for determination of M. leprae-specific antibodies. Sci Rep. 2017;7:8868.
Pethe K, Bifani P, Jang J, Kang S, Park S, Ahn S, et al. Discovery of Q203, a potent clinical candidate for the treatment of tuberculosis. Nat Med. 2013;19:1157–60.
Scherr N, Bieri R, Thomas SS, Chauffour A, Kalia NP, Schneide P, et al. Targeting the Mycobacterium ulcerans cytochrome bc1:aa3 for the treatment of Buruli ulcer. Nat Commun. 2018;9:5370.
Richter U, McFarland R, Taylor RW, Pickett SJ. The molecular pathology of pathogenic mitochondrial tRNA variants. FEBS Lett. 2021;595:1003–24.
Wong L-JC, Chen T, Wang J, Tang S, Schmitt ES, Landsverk M, et al. Interpretation of mitochondrial tRNA variants. Genet Med. 2020;22:917–26.
Lax NZ, Gnanapavan S, Dowson SJ, Alston CL, He L, Polvikoski TM, et al. Early-onset cataracts, spastic paraparesis, and ataxia caused by a novel mitochondrial tRNAGlu (MT-TE) gene mutation causing severe complex I deficiency: a clinical, molecular, and neuropathologic study. J Neuropathol Exp Neurol. 2013;72:164–75.
Ding Y, Gao B, Huang J. Mitochondrial cardiomyopathy: the roles of mt-tRNA mutations. JCM. 2022;11:6431.
Hooij A, Geluk A. In search of biomarkers for leprosy by unraveling the host immune response to Mycobacterium leprae. Immunol Rev. 2021;301:175–92.
Mi Z, Liu H, Zhang F. Advances in the immunology and genetics of leprosy. Front Immunol. 2020;11:567.
Stefano GB, Kream RM. Mitochondrial DNA heteroplasmy as an informational reservoir dynamically linked to metabolic and immunological processes associated with COVID-19 neurological disorders. Cell Mol Neurobiol. 2022;42:99–107.
Acknowledgements
We thank all participants of the study and the funding agencies of the project. We also thank Eric Baia for technical support.
Funding
This research was funded by Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq), Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES/Brazil), Fundação Amazônia de Amparo a Estudos e Pesquisas (FAPESPA) and Pró-Reitoria de Pesquisa e Pós-Graduação da Universidade Federal do Pará (PROPESP/UFPA). This work is part of Rede de Pesquisa em Genômica Populacional Humana (Biocomputacional—Protocol no. 3381/2013/CAPES). VALE S.A. 27756/2019, Fulbright Scholar to Brazil 2019–2020 (JS) and the Heiser Program of the New York Community Trust for Research in Leprosy (JB, MS, CS and JS) grants P15-000827, P16-000796 and P18-000250. The funders had no role in study design, data collection, analysis, interpretation or writing of the report.
Author information
Authors and Affiliations
Contributions
CGS, ÂRS and GCC designed the study; CGS, ARG, MBS, PFC and PP collected the data and revised the manuscript; MNSS, FGS, CSS and GCC performed the bench work; GSA, FGS and GCC performed the analyses; CGS and ÂRS acquired the resources; FGS, GSA, AHGP, MAC-D and GCC wrote and revised the manuscript; CGS, ÂRS and GCC supervised the project. All authors have read and agreed to the published version of the manuscript.
Corresponding authors
Ethics declarations
Ethics approval and consent to participate
The study was conducted in accordance with the Declaration of Helsinki and approved by the Ethics Committee of Institute of Health Sciences at the Federal University of Pará (ICS/UFPA n. 197/07). Informed consent was obtained from all subjects involved in the study.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
de Souza, F.G., da Silva, M.B., de Araújo, G.S. et al. Whole mitogenome sequencing uncovers a relation between mitochondrial heteroplasmy and leprosy severity. Hum Genomics 17, 110 (2023). https://doi.org/10.1186/s40246-023-00555-8
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s40246-023-00555-8