Evaluation of Polymorphisms in Toll-Like Receptor Genes as Biomarkers of the Response to Treatment of Erythema Nodosum Leprosum

Introduction

Leprosy is a chronic infectious disease caused by Mycobacterium leprae. It especially injures the skin and peripheral nervous system, resulting in disabilities and characteristic deformities (1). Since leprosy has different clinical features, there are several systems of classification of the disease. The most used classification systems are Ridley-Jopling (based on clinical, pathological, bacilloscopic, and immunological criteria) and the World Health Organization (WHO) (a simplified operational classification to facilitate leprosy treatment) criteria. In the Ridley-Jopling system, patients are classified as indeterminate, tuberculoid, borderline, and lepromatous. In WHO system, patients are classified as paucibacillary (PB) and multibacillary (MB), according to clinical manifestation, being the latter more severe (2, 3). Around 50% of lepromatous leprosy (LL) and 5–10% of borderline-lepromatous (BL) patients develop inflammation conditions related do leprosy. One of them, called erythema nodosum leprosum (ENL), is characterized by painful subcutaneous erythematous nodules that can ulcerate and be associated to fever, general malaise and systemic inflammation (2, 4). ENL affects patients with poor cellular immune responses but with preserved humoral responses, and its immunopathological mechanism is not yet fully understood, although high levels of the pro-inflammatory cytokine tumoral necrosis factor-alpha (TNF-α) have been associated with the pathogenesis of the reaction (2, 5, 6). In order to manage the ENL symptoms, thalidomide and prednisone are drugs frequently used. Thalidomide has a rapid action in the symptoms control within 24–48 h, mainly due to the its anti-inflammatory action by inhibiting gene expression of TNF-α (2, 7, 8). In addition, the use of systemic corticosteroids may be necessary for patients with moderate, severe, and recurrent ENL (9, 10). They act in inflammation by suppressing or inhibiting the activation of the transcription nuclear factor kappa beta (NF-kβ), which regulates genes that encode various pro-inflammatory cytokines, including TNF-α (2, 4).

Interestingly, it has been postulated that the high levels of the pro-inflammatory cytokine observed in the ENL could be related to continuous recognition of M. leprae by Toll-like receptors (TLRs) and activation of an immune response (2, 11). After infection, M. leprae bacilli are initially recognized by several innate immune receptors, including TLRs. TLRs are membrane and endosomal receptors capable of interacting with several types of molecules, such as proteins and lipopolysaccharides (4, 12, 13). To date, ten members of the TLR family have been identified in humans, which are expressed in a variety of tissues and cell types throughout the body (14, 15). Each receptor recognizes different patterns and specific ligands of microbial components, having a pivotal role on the initiation of innate immune response (16, 17). Several studies indicate that recognition of mycobacteria by TLRs represents an essential step to produce an immune response, qualified to protect against infection (9, 18–21). TLRs 1, 2, 4, and 6 have been identified as capable of recognizing M. leprae, and, after interacting with the ligand, they trigger a wide antimicrobial and inflammatory response. The innate immune response influences the type of adaptive immune response induced by the expression of genes, which are mostly regulated by interaction of the TLR (22). Transcription factors activated by the TLR signaling pathways stimulate the expression of cytokines, such TNF-α and Interleukin 1 (IL-1), chemokines, and endothelial adhesion molecules (15, 23). Studies have revealed that thalidomide and its derivatives are effective inhibitors of TLR4-induced TNF-α production, and that the inhibitory effect of thalidomide on TLR-induced cytokines and type 1 IFN production may potentially contribute to the effectiveness of this drug in treatment of ENH (24–27). These findings indicate that the treatment alters the toll-like induced immune response. Thus, the recognition of M. leprae by the different TLRs determines the immune response profile of individuals and may be associated with the response to ENL treatment.

Polymorphisms in TLRs genes have been described and some of them might modify the levels of expression of the receptors, as well as their interaction with ligands or others receptors. Moreover, genetic variants in TLR1, TLR2, TLR4, and TLR6 genes have been associated with autoimmune, infectious and inflammatory diseases, such as leprosy, pulmonary tuberculosis, atopic dermatitis, and systemic lupus erythematosus (6, 15, 22, 23, 28, 29). The single nucleotide polymorphism (SNP) rs4833095 of the TLR1 gene results in an amino acid alteration from asparagine to serine. This alteration results in decreased expression of the TLR1 receptor in the immune system cells, altering the recognition of M. leprae (30, 31). Despite TLR4 SNP, rs1927914, being located on the 5′ untranslated region (5′ UTR) of the gene, it can influence in regulatory processes, such as transcription factors' binding, altering the gene expression. It is believed that changes in the binding motifs of the transcription factors influence cellular responses related to TLRs (32, 33). SNPs in TLR2 and TLR6 were previously associated with leprosy or tuberculosis (14, 28, 34), but their specific roles are less understood. Due to the TLR relevance in inducing an immune response and previous association with infectious disease, in this study we investigated the association between thalidomide and prednisone doses with polymorphisms in TLR genes over time in the response to ENL treatment, in order to understand whether they could be useful as therapeutic molecular biomarkers.

Materials and Methods

Ethical Issues

All participants were informed about the research objectives and signed an informed consent form. This study was approved by the Ethics Committee of Hospital de Clínicas de Porto Alegre under CAAE 8989519.5.0000.5327.

Sample

The sample consisted of 162 patients with ENL that were selected in seven Brazilian centers: at Hospital de Clínicas of Porto Alegre (n = 4), Dermatology Service of Santa Casa Hospital of Porto Alegre (n = 4) and at Sanitary Dermatology Clinic (n = 47), in Porto Alegre, South Brazil; at National Reference Center of Sanitary Dermatology Dona Libania in Fortaleza (n = 15), at Humanized Reference Center of Sanitary Dermatology (n = 13) in Imperatriz and Aquiles Lisboa Hospital in São Luís (n = 65) in Northeast Brazil, and at Dermatology Ambulatory of the University of São Paulo in Monte Negro (n = 14) in the north of Brazil. The inclusion criteria were ENL patients between 18 and 85 years-old, treated with thalidomide and/or prednisone. Patients using non-steroidal anti-inflammatory drugs were excluded.

Clinical and Demographic Data Analysis

Up to ten visits previously registered in the patient's medical record were analyzed, with the collection of demographic data, including: sex, age, and place of origin, according to Brazilian political region; history of leprosy (moment of diagnosis and treatment); and history of ENL [diagnosis, types of treatment, adverse effects, history of relapse, and dose (mg) of medications used]. ENL was classified as acute, when a single episode lasting < 24 weeks; recurrent, if a second episode or following episode 28 days or more after the end of treatment; and chronic, when constant episodes were observed in more than 24 weeks (35).

Genetic Analyses

Saliva and peripheral blood were used as biological material to DNA extraction. Genomic DNA was extracted from saliva samples with Oragene^® kits (DNA Genotek), whilst peripheral blood DNA was obtained with the Flexigene^® Blood Kit (QiagenTM) extraction kit, according to manufacturers' instructions.

Single nucleotide polymorphisms (SNPs) were selected based on minor allele frequency in each gene for Brazilian population and from the literature data regarding to the genetic association studies. The rs4833095 of TLR1 (C_44103606_10), rs3804099 of TLR2 (C_22274563_10), rs1927914 of TLR4 (C_2704048_10), and rs5743810 of TLR6 (C_1180648_20) SNPs were genotyped using the Real-Time PCR technique, with TaqMan assays (Applied Biosystems USA).

Statistical Analyses

Clinical categorical variables were described in total number and percentage frequency. The medication doses were described using median and 25 and 75 percentiles. The Chi-Square Test was used to evaluate Hardy-Weinberg Equilibrium for all genetic polymorphisms and to compare the allelic frequencies of the polymorphisms, between the studied sample and the public database genetic variants in the Brazilian population (ABraOM).

To assess the association of variation in average doses of thalidomide and prednisone over time as a function of SNP, time and interaction between SNP and time, models of Generalized Estimating Equations (GEE) stratified by region (South and Northeast) were fitted. A total of 16 models were generated and the concomitant use of multidrug therapy (MDT) for treatment leprosy, use of other medications and other treatments for ENL were included as adjustment variables. The GEE method is a repeated measures analysis, focused on average changes in response over time and on the impact of covariates on these changes. The models were used with unstructured work correlation matrix, robust estimator of the covariance matrix and gamma distribution with logarithmic link function. The time variable is treated as quantitative in the models, so whenever that an interaction of the SNP with time is significant, the differences in the mean responses predicted by the model for the different categories of SNPs will be presented over time through a figure. In our analysis, if it is detected an interaction significant statistically between time and genetic variant analyzed, this result can be interpreted that average doses of the drug used in the ENL treatment are different depending on the moment of the treatment and the genotype of the individual. Therefore, every time this interaction is detected, we will describe it as a time-dependent association. SAS Studio was used to perform the Clopper-Pearson exact test. All other statistical analyzes were performed with SPSS version 18 (SPSS, www.spss.com, IBM, USA).

Results

In this study, 162 patients with ENL were included, of which 123 (75.9%) were male and 144 (88.9%) had chronic ENL. About 75% (n = 121) of the patients presented LL, and 53.1% (n = 86) were using multidrug therapy during the treatment for ENL. Patients in the South (n = 45) and Northeast (n = 84) regions used other medications during treatment for ENL. The most common adverse effects were registered in the central nervous (14.8%), gastrointestinal (16.7%) and locomotor (24.2%) systems. The demographic and clinical characteristics of the sample are presented in Table 1. It was possible identify a large difference among the Brazilian geographical regions regarding to the clinical features of the patients and the treatment used.

TABLE 1

Table 1. Clinical and demographic characteristics of ENL patients.

Genotype distributions were tested for Hardy-Weinberg Equilibrium and the allelic and genotypic frequencies of the polymorphisms are shown in Table 2. TLR4 (rs1927914) polymorphism was the only that was not in Hardy-Weinberg Equilibrium, representing an alteration in the allelic frequencies in the study population.

TABLE 2

Table 2. Genotype and allele frequency of toll-like receptors.

Table 3 shows the average doses of thalidomide and prednisone, stratified by the patient's region of origin. To represent the evolution of the treatment, the averages of medication doses were described in three time periods, in the 1st consultation, between 75 and 105 days after the 1st consultation (about 3 months after the start of treatment), and between 165 and 195 days after the 1st consultation, ~6 months after starting treatment. In the North region, patients started medication with higher mean doses of thalidomide and prednisone than patients in the Northeast and South regions. For all regions, it was possible to identify a small reduction in the average dose of medications during the period evaluated. Despite this reduction, 41 patients were still using thalidomide, and 32 patients were still receiving prednisone 6 months after starting treatment.

TABLE 3

Table 3. Average doses of medications by region and selected treatment periods.

GEE analyses were performed to assess the association of polymorphisms with the variation of the dose of thalidomide and prednisone over the time. Taking into account that triple interactions between the patient's region of origin, SNPs and time were observed, data were analyzed separately by geographical Brazilian region. Regarding treatment with thalidomide, we identified that, in the South region, the association of the genotypes of TLR1/rs4833095, TLR2/rs3804099, and TLR4/rs1927914 with the dose variation of the drug was time-dependent (interaction effect: p < 0.05; Table 4). For example, individuals with the CC genotype of TLR1/rs4833095 experienced a dose reduction throughout treatment (Figure 1A). In contrast, individuals with the AA genotype of TLR4/rs1927914 had increases in the medication dose during treatment (Figure 1D). In addition, we identified that, in the Northeast region, the association of the TLR2/rs3804099 (Figure 1C) and TLR6/rs5783810 (Figure 1E) genotypes with the variation of thalidomide dose also were time-dependent. For TLR2/rs3804099 we identified that individuals with CC and TT genotypes started treatment with a higher mean dose of thalidomide, but showed a greater reduction in mean dose over the treatment period when compared to CT heterozygotes. As for TLR6/rs5783810, individuals with AA genotype showed an increase in the average dose of thalidomide along the time of treatment evaluated, while individuals with the AG or GG genotype showed a reduction in the average dose. Regarding treatment with prednisone, we were able to identify that, in the South region; the association of the genotypes of TLR1/rs4833095, TLR2/rs3804099, and TLR6/rs5783810 with the dose variation of the drug was time-dependent. As for the Northeast region, only the association of the TLR6/rs5783810 genotypes with the dose variation was time-dependent (Table 4). Despite these interactions, we observed that, the difference in the average dose of both drugs, between the genetic variants evaluated, did not present impact on the ENL treatment. For example, the biggest difference in the average dose of thalidomide, between the genetic variants, was 52 mg among patients with genotype AA and GG of TLR6/rs5783810 in the South, representing only half of a tablet of the drug (100 mg) (Supplementary Table 1). We did not perform GEE analyses for the North region data due to the small sample size.

TABLE 4

Table 4. P-values for the association of genotypes of TLR polymorphism, time and genotype^*time interaction on dose of thalidomide or prednisone, as obtained in eight separate analyses for each drug and SNP using the Generalized Estimates Equation model (GEE).

FIGURE 1

Figure 1. Lines charts of the thalidomide dose averages with Toll-like receptor genotypes discriminated by the times evaluated. (A) TLR1/ rs4833095 in the South region: blue line, CC genotype; red line, CT genotype; green line, TT genotype. (B) TLR2/rs3804099 in the South region: blue line, CC genotype; red line, CT genotype; green line, TT genotype. (C) TLR2/rs3804099 in the Northeast region: blue line, CC genotype; red line, CT genotype; green line, TT genotype. (D) TLR4/rs1927914 in the South region: blue line, AA genotype; red line, AG genotype; green line, GG genotype. (E) TLR6/ rs5783810 in the Northeast region: blue line, AA genotype; red line, AG genotype; green line, GG genotype.

Discussion

This study aimed to evaluate the association of genetic variants of TLRs genes with the response to the treatment of ENL with thalidomide and prednisone. We found that the association of TLRs genes polymorphisms with thalidomide and prednisone dose variation was time-dependent. These interactions between time and polymorphism mean that the association of the genetic variant with the dose of drug varies throughout the treatment, i.e., the association with the genetic variant is different depending on the moment of the treatment. Another interesting finding was the huge difference among Brazilian geographical regions regarding the treatment and clinical features of patients. This could be related to the different sample sizes, the heterogeneous genetic background of the Brazilian population and/or the treatment scheme used. Furthermore, during the statistical analyses, triple interactions between the patient's region of origin, SNPs, and time were observed. Thus, we decided to split our analysis by region. In spite of these regional differences, patients enrolled in this study were mainly male, diagnosed with lepromatous leprosy, chronic ENL, and using multidrug therapy and other medications during the treatment of ENL. The predominance of men is in accordance with the literature, which reports that proportionally more men are diagnosed with leprosy and ENL than women, being men more frequently affected by severe disabilities (36, 37). The greater number of patients with lepromatous leprosy, which is the form mainly associated with ENL, could indicate that the bacillary index is a risk factor for the development of the reaction as reported in other studies (36, 38, 39). In addition, the concomitant use of MDT in the present sample corroborate that the reaction manifests itself mainly in the first year of treatment of the disease (24, 40). Together, the findings support the hypothesis that the amount of soluble antigen during treatment, due to bacterial death, causes excessive activation of the immune system, inducing an inflammatory state, especially in those patients with high bacillary load at the beginning of therapy (2, 11).

Different constituents of M. leprae have been characterized as potent TLR ligands and stimulators of the immune system, such as the triacylated lipoproteins that have been indicated as activators of TLR2/TLR1 heterodimers and TLR2 homodimer (41). The TLR1 gene encodes the receptor that recognizes M. leprae glycolipids and lipopolysaccharides and activates NF-kβ, which induces the expression of pro-inflammatory genes. The rs4833095 SNP involves a change from asparagine to serine at the position 248 of the protein, altering the expression of the TLR1 receptor in cells of the immune system (30, 31). Individuals with the T allele were shown to have decreased expression of TLR1 on the leukocyte surface. A defect in cell traffic is believed to affect the transport of the receptor to the cell membrane, altering its expression and recognition of M. leprae. Thus, the immune response generated becomes insufficient (31, 42). These data may indicate that individuals homozygous for the T allele of this SNP may show less signaling and inflammatory activity in response to M. leprae, resulting in the need for lower doses of anti-inflammatory drugs. Corroborating this hypothesis, in the present study we identified that, in the South region, despite having the dose reduced over the course of treatment, individuals with the CC genotype started treatment with an average higher dose of thalidomide (Figure 1A). This could indicate that individuals homozygous for the C allele of this SNP have greater inflammatory activity. However, we were unable to identify dose variation over time for individuals with CT and TT genotype. In addition, the dose reduction for individuals of the CC genotype was not significant for clinical management.

The SNP rs3804099 of TLR2 is a synonymous polymorphism (43, 44). Synonymous variants can be functional by affecting the expression and activity of the receptor, by altering the transcription factor binding and/or the stability of the mRNA (43, 45). However, the exact mechanism related to the variant and TLR2 activity remains unclear (45). Kang and Chae reported that the TT genotype of the rs3804099/TLR2 polymorphism was associated with susceptibility to leprosy. The authors suggested that this variant (T) could affect the signaling function of the receptor, resulting in a deficient immune response (46, 47). In agreement, Santana et al. reported that the T allele of this SNP was positively associated with an increased risk of developing leprosy per se in a population in Northeast Brazil. In addition, in the population studied by Santana et al., carriers of the T allele were associated with higher serum levels of IL-17 and greater production of IL-6 (14). In this study, we identified that, in the South region, individuals with TT genotype started treatment with an average dose of thalidomide higher than those with CC genotype (Figure 1B). It may indicate that individuals homozygous for the T allele have an immunological profile with greater inflammatory activity, requiring higher doses of the drug. Moreover, they showed a dose reduction over the represented treatment period, while for individuals with CC genotype, we were unable to identify dose variation over the treatment period. Despite this, the reduction in the average dose, about 0.3% in the average dose/day, does not represent a clinical difference for medical management and the hypothesis about the role of T allele remains to be clarified.

The rs1927914 polymorphism of TLR4 was not in Hardy-Weinberg Equilibrium. The Hardy-Weinberg principle is based on population genetics (48). Hardy-Weinberg deviation has been proposed as a measure of disease association when analyzing a disease group, since the individuals are a non-random selection of the population based on a phenotype of interest (49–52). In fact, leprosy has been associated with rs1927914 polymorphism (14), and that association could be explaining the Hardy-Weinberg imbalance found in our results. This finding could represent a bias; however, we did not remove it from of our analysis because the study sample was selected with leprosy patients, and therefore, it is possible that the risk allele is enriched in the entire sample. Moreover, we evaluated response to the treatment of ENL, hence the association between SNP and leprosy does not interfere in our analysis. The SNP rs1927914 of TLR4 is located in the 5′-UTR. Thus, although the functional role of the SNP remains unknown, it is possible that it can influence gene expression, interfering with transcription factor binding and regulating the activity of the promoter (32, 53). Despite the classic TLR4 ligand being lipopolysaccharide (LPS), this receptor demonstrated the ability to recognize M. leprae and M. tuberculosis (9, 54). Santana et al. (14) observed significant differences in the levels of cytokines IL-17 and IL-1β in leprosy patients carrying the genotype AA, which produced more of both cytokines compared to individuals carrying the genotypes AG or GG together (14). In the present study, the results were different, because individuals with the GG genotype had an increase in the average dose of thalidomide over the course of treatment, when compared with individuals with the AA genotype. Clearly, other functional and molecular experiments with a genetic background different sample should be performed in order to understand the association between SNP rs1927914 and time in response to ENL treatment.

SNP rs5783810 is located in a coding region of TLR6 and it promotes a change from serine to proline at the position 249 in the extracellular domain of the protein. This amino acid exchange results in the functional deterioration of TLR6 and predisposes individuals to dysregulation of the innate immune system (55, 56). Mattos et al. (20) suggested that TLR6-dependent signaling in interactions between M. leprae and Schwann cells plays a critical role favoring phagocytosis and signaling, to induce lipid droplet biogenesis in infected cells, thus providing mycobacterial survival and persistence of the microorganism in the nerve (20). The SNP rs5743810 of TLR6 has been associated with susceptibility to tuberculosis among ethnic African, European or Hispanic groups (30). It has been proposed that this polymorphism influences the recognition of ligands and reduces cell signaling, with the A allele and the AG and AA genotypes, demonstrating to have a protective effect against the development of tuberculosis (30). In this study, individuals with AA genotype showed an increase in the average dose of thalidomide along the time of treatment, while individuals with the AG or GG genotype showed a reduction in the average dose over the course of treatment in the Northeast region. Lower cell signaling induced by the G allele of this SNP could imply less inflammatory response, resulting in the need to administer lower mean doses of the medication in individuals with the GG genotype. However, the group of individuals with the AA genotype is composed by a lower number of individuals, making the results interpretation difficult. Thus, further studies should be carried out to clarify the association between the SNP rs5743810 and the time in response to the treatment of ENL.

For prednisone treatment, we identified associations between polymorphisms in TLR genes and treatment time (Table 4). In the South region, we observed a negative association between the dose of the drug and the treatment time for the genotypes of TLR1/rs4833095 (CT and TT) and TLR2/rs3804099 (CC and CT). That is, there was a reduction in the average dose of the drug throughout the treatment. In addition, this association has also been identified in individuals with genotypes AA and AG of TLR6/rs5783810. In the Northeast, in all genotypes of TLR6/rs5783810, there was a negative association between dose and treatment time. However, in none of the interactions the mean dose variation was clinically significant, being with little or no impact on the clinical conduct of ENL.

Taking into account all results, we observed a reduction in the dose averages over the evaluated period. However, after around 6 months of treatment, 41 of the 134 patients who used thalidomide were still using the medication, whilst 32 of the 103 initial patients using prednisone remained receiving the drug. These data confirm that many individuals have chronic ENL, being treated with thalidomide and corticosteroids for prolonged periods of months. In addition, we showed that, in the three periods of time evaluated, patients of the North region received an average dose of both drugs higher than patients of the South region (Table 3). This finding could be related to the difference in sample size between the regions and/or the genetic heterogeneity of the Brazilian population (57, 58). However, it can also be related to the treatment regimen with thalidomide, since currently there is no standard dosage or indication for thalidomide dose reduction in the ENL in Brazil, and each center can adopt different and local protocols related to this medication (2, 39). Nevertheless, prolonged use of thalidomide and prednisone and/or high doses can cause serious adverse effects. Here, we grouped the adverse effects registered throughout of the treatment according to the symptoms. The most common were central nervous, gastrointestinal, and locomotor systems related. All of these effects could be associated with the use of thalidomide, since they are known causes of discontinuing the use of the drug (2, 59). However, we were unable to differentiate between adverse effects due to the use of thalidomide or prednisone, as well as those caused by leprosy. Intriguingly, we found no reports of sedation, the most common adverse effect associated with the use of thalidomide (60). This can be due to the lack of description of this effect in the patients' medical records, rather than the lack occurrence of the effect. Indeed, we noticed that the records of adverse effects were different when we stratifying the sample by region. Thus, the number of individuals with reports of specific adverse effects was significantly reduced in all the regions, making any causal inference about association between those effects and genetic factors likely spurious. For these reasons, we decided not to assess the association of TLRs with the occurrence of adverse effects related to the treatment of ENL with thalidomide and prednisone.

In this sense, there are other limitations in our study that must be taking into account in regard to the results interpretation here presented. Although we have had used a statistic model that controls confounding variables that could influence on the doses of the drugs, it is possible that other confounding factors may not have been considered in these models. One example is related to patient's origin. They were collected in different regions of the country and the heterogeneous genetic background and ancestry of the Brazilian population (57) may have been underestimated in the analysis and interpretation of the results. In addition, the limited sample size of this study can make it difficult to identify small effects of genetic variants in different outcomes for specific ancestries. Furthermore, this may also have influenced the large difference found between regions regarding to the average doses of thalidomide. Currently there is a wide therapeutic dose range of thalidomide for ENL in Brazil (100–400 mg/day) (61), and each center can adopt different protocols related to this medication and severity of disease. The analysis stratified by regions minimizes the potential biases of the different clinical settings. Finally, we conducted a retrospective study in which the patient's clinical data, therapeutic conduct and related adverse effects were obtained from medical records. Therefore, some important information could be missing in our analyses, because of the study design. The lack of standardization in the description of such information or the absence of reports cannot be excluded as well.

From the results of this study, we identified genetic variants in the TLR genes that may be associated with the response to treatment with thalidomide and prednisone. The association with the genotypes was dependent on the moment of treatment and patient's region of origin. However, the association of the SNPs on dose variation was not enough to suggest a change in the therapeutic conduct according to patient's genotype (Supplementary Table 1). ENL is a complex, chronic, and difficult to control condition, and all drugs used in its treatment require careful consideration (2). Thalidomide is an effective drug, but has use restrictions related to its teratogenicity and induced peripheral neuropathy (2, 62). Similarly, the prolonged use of corticoids causes several life-threatening adverse effects (63). Therefore, it is important to identify biomarkers that indicate response to treatment, in order to restrict the use of medications only to patients who will benefit the most of them. Here, we did not find biomarkers with clinical relevance to be used in the monitoring the treatment of ENL with thalidomide or prednisone. However, there are still many gaps to be filled in our knowledge of the interaction mechanism between the different TLRs and M. leprae and how these genes can interfere in the course and treatment of the disease. Nevertheless, despite TLR relevance in bacillus recognition, inducing an immune response and the association between genetic variants and ENL treatment showed here, genetic polymorphisms in TLRs may not be useful biomarkers in response to the treatment of ENL. Further studies analyzing the association of SNPs with the dose variation over time should be carried out to confirm this hypothesis, in larger samples with more homogeneous genetic features and assessment of ancestry.

Data Availability Statement

Data cannot be shared publicly because of specific authorization was not asked to the patients. Requests to access the datasets should be directed to Hospital de Clinicas de Porto Alegre Institutional Ethics Committee, cep@hcpa.edu.br.

Ethics Statement

The studies involving human participants were reviewed and approved by Ethics Committee of Hospital de Clínicas de Porto Alegre-HCPA. The patients/participants provided their written informed consent to participate in this study.

Author Contributions

MM-F contributed to obtaining the specimens, designing and conducting the experiments, performing the statistical analyses, and writing the manuscript. PC contributed to devising the concept and designing the experiments. TK, LS-F, RB, RV, and LE contributed to designing the experiment and correcting the manuscript. PM, MS, and LC contributed to obtaining the specimens and correcting the manuscript. SC-J and SC contributed to designing the statistical analyses and writing the manuscript. FV contributed to devising the concept, designing the experiments, supervising the analyses, writing and correcting the manuscript. All authors discussed the results and contributed to scientifically to the manuscript.

Funding

This work was supported by Instituto Nacional de Genética Médica Populacional (INAGEMP) (Grant Nos. CNPq 573993/2008-4 and FAPERGS 17/2551.0000521-0), Fundo de Incentivo à Pesquisa e Eventos (FIPE) of the Hospital de Clínicas de Porto Alegre (HCPA) (Grant No. 2019-0155), Coordenação Brasileira de Aperfeiçoamento de Pessoal de Nível Superior (CAPES), and Fundação de Amparo à Pesquisa do Rio Grande do Sul- FAPERGS (Grant No. 19/2551-0001787-1).

Conflict of Interest

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.

Publisher's Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

We would like to thank the collaborators of the participating centers who helped in obtaining the specimens.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmed.2021.713143/full#supplementary-material

References

1. Britton WJ, Lockwood DNJ. Leprosy. Lancet. (2004) 363:1209–19. doi: 10.1016/S0140-6736(04)15952-7

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Costa P do SS, Fraga LR, Kowalski TW, Daxbacher ELR, Schuler-Faccini L, Vianna FSL. Erythema nodosum leprosum: update and challenges on the treatment of a neglected condition. Acta Trop. (2018) 183:134–41. doi: 10.1016/j.actatropica.2018.02.026

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Gurung P, Gomes CM, Vernal S, Leeflang MMG. Diagnostic accuracy of tests for leprosy: a systematic review and meta-analysis. Clin Microbiol Infect. (2019) 25:1315–27. doi: 10.1016/j.cmi.2019.05.020

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Tartey S, Takeuchi O. Pathogen recognition and Toll-like receptor targeted therapeutics in innate immune cells. Int Rev Immunol. (2017) 36:57–73. doi: 10.1080/08830185.2016.1261318

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Mendonça VA, Costa RD, De Melo GEBA, Antunes CM, Teixeira AL. Imunologia da hanseníase. An Bras Dermatol. (2008) 83:343–50. doi: 10.1590/S0365-05962008000400010

CrossRef Full Text | Google Scholar

6. Fonseca AB de L, Simon M do V, Cazzaniga RA, de Moura TR, de Almeida RP, Duthie MS, et al. The influence of innate and adaptative immune responses on the differential clinical outcomes of leprosy. Infect Dis Poverty. (2017) 6:1–8. doi: 10.1186/s40249-016-0229-3

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Fischer M. Leprosy – an overview of clinical features, diagnosis, and treatment. J Ger Soc Dermatology. (2017) 15:801–27. doi: 10.1111/ddg.13301

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Kim EJ, Lee JG, Kim JY, Song SH, Joo DJ, Huh KH, et al. Enhanced immune-modulatory effects of thalidomide and dexamethasone co-treatment on T cell subsets. Immunology. (2017) 152:628–37. doi: 10.1111/imm.12804

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Polycarpou A, Holland MJ, Karageorgiou I, Eddaoudi A, Walker SL, Willcocks S, et al. Mycobacterium leprae activates toll-like receptor-4 signaling and expression on macrophages depending on previous bacillus Calmette-Guerin vaccination. Front Cell Infect Microbiol. (2016) 6:72. doi: 10.3389/fcimb.2016.00072

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Rupasree Y, Naushad SM, Rajasekhar L, Uma A, Kutala VK. Association of TLR4 (D299G, T399I), TLR9−1486T>C, TIRAP S180L and TNF-α promoter (-1031,−863,−857) polymorphisms with risk for systemic lupus erythematosus among South Indians. Lupus. (2015) 24:50–7. doi: 10.1177/0961203314549792

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Schreuder PAM, Noto S, Richardus JH. Epidemiologic trends of leprosy for the 21st century. Clin Dermatol. (2016) 34:24–31. doi: 10.1016/j.clindermatol.2015.11.001

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Anthoney N, Foldi I, Hidalgo A. Toll and toll-like receptor signalling in development. Dev. (2018) 145:1–6. doi: 10.1242/dev.156018

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Antunes DE, Goulart IMB, Lima MIS, Alves PT, Tavares PCB, Goulart LR. Differential expression of IFN- γ, IL-10, TLR1, and TLR2 and their potential effects on downgrading leprosy reaction and erythema nodosum leprosum. J Immunol Res. (2019) 2019:3405103. doi: 10.1155/2019/3405103

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Santana N, de L, Rêgo JL, Oliveira JM, de Almeida LF, Braz M, Machado LMM, et al. Polymorphisms in genes TLR1, 2 and 4 are associated with differential cytokine and chemokine serum production in patients with leprosy. Mem Inst Oswaldo Cruz. (2017) 112:260–8. doi: 10.1590/0074-02760160366

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Abbas AK, Lichtman AH, Pobe JS. Imunologia Celular e Molecular. Rio de Janeiro: Elsevier (2015). p. 823 a 860.

16. Polycarpou A, Walker SL, Lockwood DNJ. A systematic review of immunological studies of erythema nodosum leprosum. Front Immunol. (2017) 8:233. doi: 10.3389/fimmu.2017.00233

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Takeda K, Akira S. Toll-Like receptors. Curr Protoc Immunol. (2015) 2015:14.12.1-14.12.10. doi: 10.1002/0471142735.im1412s109

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Krutzik SR, Ochoa MT, Sieling PA, Uematsu S, Ng YW, Legaspi A, et al. Activation and regulation of Toll-like receptors 2 and 1 in human leprosy. Nat Med. (2003) 9:525–32. doi: 10.1038/nm864

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Oliveira RB, Ochoa MT, Sieling PA, Rea TH, Rambukkana A, Sarno EN, et al. Expression of toll-like receptor 2 on human schwann cells: a mechanism of nerve damage in leprosy. Infect Immun. (2003) 71:1427–33. doi: 10.1128/IAI.71.3.1427-1433.2003

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Mattos KA, Oliveira VGC, D'Avila H, Rodrigues LS, Pinheiro RO, Sarno EN, et al. TLR6-driven lipid droplets in mycobacterium leprae- infected schwann cells: immunoinflammatory platforms associated with bacterial persistence. J Immunol. (2011) 187:2548–58. doi: 10.4049/jimmunol.1101344

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Mattos KA, D'Avila H, Rodrigues LS, Oliveira VGC, Sarno EN, Atella GC, et al. Lipid droplet formation in leprosy: toll-like receptor-regulated organelles involved in eicosanoid formation and Mycobacterium leprae pathogenesis. J Leukoc Biol. (2010) 87:371–84. doi: 10.1189/jlb.0609433

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Brightbill HD, Modlin RL. Toll-like receptors: molecular mechanisms of the mammalian immune response. Immunology. (2000) 101:1–10. doi: 10.1046/j.1365-2567.2000.00093.x

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Medzhitov R, Janeway CA. Innate immunity: impact on the adaptive immune response. Curr Opin Immunol. (1997) 9:4–9. doi: 10.1016/S0952-7915(97)80152-5

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Dias AA, Silva CO, Santos JPS, Batista-Silva LR, Acosta CCD, Fontes ANB, et al. DNA sensing via TLR-9 constitutes a major innate immunity pathway activated during erythema nodosum leprosum. J Immunol. (2016) 197:1905–13. doi: 10.4049/jimmunol.1600042

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Millrine D, Miyata H, Tei M, Dubey P, Nyati K, Nakahama T, et al. Immunomodulatory drugs inhibit TLR4-induced type-1 interferon production independently of Cereblon via suppression of the TRIF/IRF3 pathway. Int Immunol. (2016) 28:307–15. doi: 10.1093/intimm/dxw005

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Millrine D, Kishimoto T. Review A brighter side to thalidomide : its potential use in immunological disorders. Trends Mol Med. (2017) xx:1–14. doi: 10.1016/j.molmed.2017.02.006

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Negera E. The effects of prednisolone treatment on cytokine expression in patients with erythema nodosum leprosum reactions. Front Immunol. (2018) 9:189. doi: 10.3389/fimmu.2018.00189

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Bochud PY, Hawn TR, Siddiqui MR, Saunderson P, Britton S, Abraham I, et al. Toll-like receptor 2 (TLR2) polymorphisms are associated with reversal reaction in leprosy. J Infect Dis. (2008) 197:253–61. doi: 10.1086/524688

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Suryadevara NC, Neela VSK, Kovvali S, Pydi SS, Jain S, Siva Sai KSR, et al. Genetic association of G896A polymorphism of TLR4 gene in leprosy through family-based and case-control study designs. Trans R Soc Trop Med Hyg. (2013) 107:777–82. doi: 10.1093/trstmh/trt084

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Ma X, Liu Y, Gowen BB, Graviss EA, Clark AG, Musser JM. Full-exon resequencing reveals toll-like receptor variants contribute to human susceptibility to tuberculosis disease. PLoS ONE. (2007) 2:e1318. doi: 10.1371/journal.pone.0001318

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Niitsuma ENA, Da Rocha Fernandes G, Lana FCF. The TLR1 gene is associated with higher protection from leprosy in women. PLoS ONE. (2018) 13:1–10. doi: 10.1371/journal.pone.0205234

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Li T, Zhang X, Sang L, Li XT, Sun HY, Yang J, et al. The interaction effects between TLR4 and MMP9 gene polymorphisms contribute to aortic aneurysm risk in a Chinese Han population. BMC Cardiovasc Disord. (2019) 19:1–11. doi: 10.1186/s12872-019-1049-8

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Kim JS, Kim SK, Park JY, Kim YG, Moon JY, Lee SH, et al. Significant association between toll-like receptor gene polymorphisms and posttransplantation diabetes mellitus. Nephron. (2016) 133:279–86. doi: 10.1159/000446570

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Shey MS, Randhawa AK, Bowmaker M, Smith E, Scriba TJ, De Kock M, et al. Single nucleotide polymorphisms in toll-like receptor 6 are associated with altered lipopeptide- and mycobacteria-induced interleukin-6 secretion. Genes Immun. (2010) 11:561–72. doi: 10.1038/gene.2010.14

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Walker SL, Lebas E, Doni SN, Lockwood DNJ, Lambert SM. The mortality associated with erythema nodosum leprosum in ethiopia: a retrospective hospital-based study. PLoS Negl Trop Dis. (2014) 8:e2690. doi: 10.1371/journal.pntd.0002690

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Guerra JG, Penna GO, Castro LCM De, Martelli CMT, Stefani MMA, Costa MB. Erythema nodosum leprosum case series report: clinical profile, immunological basis and treatment implemented in health services. Rev Soc Bras Med Trop. (2004) 37:384–90. doi: 10.1590/S0037-86822004000500003

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Ye Q, Chen B, Tong Z, Nakamura S, Sarria R, Costabel U, et al. Thalidomide reduces IL-18, IL-8 and TNF-α release from alveolar macrophages in interstitial lung disease. Eur Respir J. (2006) 28:824–31. doi: 10.1183/09031936.06.00131505

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Pocaterra L, Jain S, Reddy R, Muzaffarullah S, Torres O, Suneetha S, et al. Clinical course of erythema nodosum leprosum: an 11-year cohort study in Hyderabad, India. Am J Trop Med Hyg. (2006) 74:868–79. doi: 10.4269/ajtmh.2006.74.868

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Walker SL, Balagon M, Darlong J, Doni SN, Hagge DA, Halwai V, et al. ENLIST 1: an international multi-centre cross-sectional study of the clinical features of erythema nodosum leprosum. PLoS Negl Trop Dis. (2015) 9:e0004065. doi: 10.1371/journal.pntd.0004065

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Voorend CGN, Post EB. A systematic review on the epidemiological data of erythema nodosum leprosum, a type 2 leprosy reaction. PLoS Neglected Tropical Diseases. (2013) 7:e2440. doi: 10.1371/journal.pntd.0002440

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Geluk A. Correlates of immune exacerbations in leprosy. Semin Immunol. (2018) 39:111–8. doi: 10.1016/j.smim.2018.06.003

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Uciechowski P, Imhoff H, Lange C, Meyer CG, Browne EN, Kirsten DK, et al. Susceptibility to tuberculosis is associated with TLR1 polymorphisms resulting in a lack of TLR1 cell surface expression. J Leukoc Biol. (2011) 90:377–88. doi: 10.1189/jlb.0409233

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Mirkamandar E, Nemati M, Hayatbakhsh MM, Bassagh A, Khosravimashizi A, Jafarzadeh A. Association of a single nucleotide polymorphism in the TLR2 gene (rs3804099), but not in the TLR4 gene (rs4986790), with Helicobacter pylori infection and peptic ulcer. Turkish J Gastroenterol. (2018) 29:283–91. doi: 10.5152/tjg.2018.17484

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Zhang J, Zhao Z, Zhong H, Wu L, Zhou W, Peng W, et al. Importance of common TLR2 genetic variants on clinical phenotypes and risk in tuberculosis disease in a Western Chinese population. Infect Genet Evol. (2018) 60:173–80. doi: 10.1016/j.meegid.2018.02.031

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Press D. No genetic relationship between TLR2 rs4696480, rs3804100, and rs3804099 gene polymorphisms. Onco Targets Ther. (2017) 10:2325–33. doi: 10.2147/OTT.S121618

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Noman ASM, Koide N, Hassan FI.-E-Khuda I, Dagvadorj J, Tumurkhuu G, et al. Thalidomide inhibits lipopolysaccharide-induced tumor necrosis factor-α production via down-regulation of MyD88 expression. Innate Immun. (2009) 15:33–41. doi: 10.1177/1753425908099317

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Kang TJ, Chae GT. Detection of Toll-like receptor 2 (TLR2) mutation in the lepromatous leprosy patients. FEMS Immunol Med Microbiol. (2001) 31:53–8. doi: 10.1111/j.1574-695X.2001.tb01586.x

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Namipashaki A, Razaghi-moghadam Z, Ansari-pour N. The essentiality of reporting hardy-weinberg equilibrium calculations in population-based genetic association studies (debate). Cell J. (2015) 17:187–92. doi: 10.22074/cellj.2016.3711

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Domingo R, Ellis MC, Fullan A, Hinton LM, Jones NL, Kimmel BE, et al. A novel MHC class I-like gene is mutated in patients with hereditary haemochromatosis. Nat Genet. (1996) 13:399–408. doi: 10.1038/ng0896-399

PubMed Abstract | CrossRef Full Text | Google Scholar

50. Lee W. Searching for disease-susceptibility loci by testing for Hardy-Weinberg disequilibrium in a gene bank of affected individuals. Am J Epidemiol. (2003) 158:397–400. doi: 10.1093/aje/kwg150

PubMed Abstract | CrossRef Full Text | Google Scholar

51. Nielsen DM, Ehm MG, Weir BS. Detecting marker-disease association by testing for Hardy-Weinberg disequilibrium at a marker locus. Am J Hum Genet. (1999) 63:1531–40. doi: 10.1086/302114

PubMed Abstract | CrossRef Full Text | Google Scholar

52. Song K, Elston RC. A powerful method of combining measures of association and Hardy – Weinberg disequilibrium for ýne-mapping in case-control studies. Stat Med. (2006) 25:105–26. doi: 10.1002/sim.2350

PubMed Abstract | CrossRef Full Text | Google Scholar

53. Gu L, Huang J, Liang B, Chen Q, Xie J, Yang J, et al. TLR4 polymorphisms affect stroke risk and inflammatory response in Chinese ischemic stroke patients. Neurol Sci. (2018) 39:127–33. doi: 10.1007/s10072-017-3151-y

PubMed Abstract | CrossRef Full Text | Google Scholar

54. Hart BE, Tapping RI. Genetic diversity of toll-like receptors and immunity to M. leprae infection. J Trop Med. (2012) 2021:415057. doi: 10.1155/2012/415057

PubMed Abstract | CrossRef Full Text | Google Scholar

55. Hoffjan S, Stemmler S, Parwez Q, Petrasch-Parwez E, Arinir U, Rhode G, et al. Evaluation of the toll-like receptor 6 Ser249Pro polymorphism in patients with asthma, atopic dermatitis and chronic obstructive pulmonary disease. BMC Med Genet. (2005) 6:1–6. doi: 10.1186/1471-2350-6-34

PubMed Abstract | CrossRef Full Text | Google Scholar

56. Semlali A, Almutairi M, Rouabhia M, Parine NR, Amri A Al, Al-Numair NS, et al. Novel sequence variants in the TLR6 gene associated with advanced breast cancer risk in the Saudi Arabian population. PLoS ONE. (2018) 13:1–13. doi: 10.1371/journal.pone.0203376

PubMed Abstract | CrossRef Full Text | Google Scholar

57. Fernanda Lima-Costa M, Rodrigues LC, Barreto ML, Gouveia M, Horta BL, Mambrini J, et al. Genomic ancestry and ethnoracial self-classification based on 5,871 community-dwelling Brazilians (The Epigen Initiative). Sci Rep. (2015) 5:9812. doi: 10.1038/srep09812

PubMed Abstract | CrossRef Full Text | Google Scholar

58. Kehdy FSG, Gouveia MH, Machado M, Magalhães WCS, Horimoto AR, Horta BL, et al. Origin and dynamics of admixture in Brazilians and its effect on the pattern of deleterious mutations. Proc Natl Acad Sci USA. (2015) 112:8696–701. doi: 10.1073/pnas.1504447112

PubMed Abstract | CrossRef Full Text | Google Scholar

59. Walker SL, Waters MFR, Lockwood DNJ. The role of thalidomide in the management of erythema nodosum leprosum. Lepr Rev. (2007) 78:197–215. doi: 10.47276/lr.78.3.197

PubMed Abstract | CrossRef Full Text | Google Scholar

60. Franks ME, Macpherson GR, Figg WD. Thalidomide. Encycl Toxicol. (2014) 363:523–6. doi: 10.1016/S0140-6736(04)16308-3

PubMed Abstract | CrossRef Full Text | Google Scholar

61. Ministério da Saúde. Guia prático sobre a hanseníase. Secretaria de Vigilância em Saúde Departamento de Vigilância e Doenças Transmissíveis. (2017). 70 p. Available online at: https://bvsms.saude.gov.br/bvs/publicacoes/guia_pratico_hanseniase.pdf (accessed November 25, 2021).

62. Mcbride WG. Thalidomide and congenital abnormalities. Lancet. (1961) 278:1358. doi: 10.1016/S0140-6736(61)90927-8

CrossRef Full Text | Google Scholar

63. Sugumaran DST. Leprosy reactions - complications of steroid therapy. Int J Lepr Other Mycobact Dis. (1998) 66:10–5.

PubMed Abstract | Google Scholar