Acessibilidade / Reportar erro

Evaluation of the influence of genetic variants in Cereblon gene on the response to the treatment of erythema nodosum leprosum with thalidomide

Abstract

BACKGROUND

Erythema nodosum leprosum (ENL) is an acute and systemic inflammatory reaction of leprosy characterised by painful nodules and involvement of various organs. Thalidomide is an immunomodulatory and anti-inflammatory drug currently used to treat this condition. Cereblon (CRBN) protein is the primary target of thalidomide, and it has been pointed out as necessary for the efficacy of this drug in others therapeutics settings.

OBJECTIVES

In this study, we aimed to evaluate the influence of CRBN gene variants on the dose of thalidomide as well as its adverse effects during treatment of ENL.

METHODS

A total of 103 ENL patients in treatment with thalidomide were included in this study. DNA samples were obtained from saliva and molecular analysis of CRBN gene were performed to investigate the variants rs1620675, rs1672770 and rs4183. Different genotypes of CRBN variants were evaluated in relation to their influence on the dose of thalidomide and on the occurrence of adverse effects.

FINDINGS

No association was found between CRBN variants and thalidomide dose variation. However, the genotypes of rs1672770 showed association with gastrointestinal effects (p = 0.040). Moreover, the haplotype DEL/C/T (rs4183/rs1672770/rs1620675) was also associated with gastrointestinal adverse effects (p = 0.015).

MAIN CONCLUSIONS

Our results show that CRBN variants affect the treatment of ENH with thalidomide, especially on the adverse effects related to the drug.

Key words:
CRBN; leprosy; pharmacogenomics; pharmacogenetics; personalised medicine


ENL is a difficult to manage leprosy complication and a potentially disabling condition.11. Chandler DJ, Hansen KS, Mahato B, Darlong J, John A, Lockwood DN. Household costs of leprosy reactions (ENL) in rural India. PLoS Negl Trop Dis. 2015; 9(1): e0003431.,22. 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(5): 868-79. It is an inflammatory reaction characterised by painful nodules on the skin that can ulcerate and by systemic involvement with fever and general malaise and effects on various organs.22. 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(5): 868-79.,33. 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 Tropica. 2018; 183: 134-41. The reaction affects patients with borderline lepromatous (BL) leprosy and lepromatous leprosy (LL) that are associated with a higher bacillary load.11. Chandler DJ, Hansen KS, Mahato B, Darlong J, John A, Lockwood DN. Household costs of leprosy reactions (ENL) in rural India. PLoS Negl Trop Dis. 2015; 9(1): e0003431.,22. 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(5): 868-79.

Thalidomide is effective in treating ENL, rapidly reducing symptoms, such as fever and night sweats, and improving skin lesions.44. Melchert M, List A. The thalidomide saga. Int J Biochem Cell Biol. 2007; 39(7-8): 1489-99.,55. Sheskin J. Thalidomide in the treatment of lepra reactions. Clin Pharmacol Ther. 1965; 6(3): 303-6. Thalidomide is a glutamic acid derivative produced initially in the 1950s as an anticonvulsant and subsequently used as a sedative and antiemetic in early pregnancy.44. Melchert M, List A. The thalidomide saga. Int J Biochem Cell Biol. 2007; 39(7-8): 1489-99.,66. Mercurio A, Adriani G, Catalano A, Carocci A, Rao L, Lentini G, et al. A Mini-review on thalidomide: chemistry, mechanisms of action, therapeutic potential and anti-angiogenic properties in multiple myeloma. Curr Med Chem. 2017; 24: 2736-44. However, its teratogenic effect became known when reports of birth defects in children born to mothers who used thalidomide during pregnancy appeared in 1961.77. Mcbride WG. Thalidomide and congenital abnormalities. The Lancet. 1961; 278(7216): 1358.,88. Lenz W, Pfeiffer RA, Kosenow W, Hayman DJ. Thalidomide and congenital abnormalities. The Lancet. 1962; 279(7219): 45-6. Thalidomide was soon withdrawn from the market; however its efficacy in the treatment of erythema nodosum leprosum (ENL) was accidently discovered few years later, bringing it back for the clinical use even within strict regulations.55. Sheskin J. Thalidomide in the treatment of lepra reactions. Clin Pharmacol Ther. 1965; 6(3): 303-6. By 1999, thalidomide was also proved to be effective in the treatment of multiple myeloma (MM)99. Ito T, Handa H. Cereblon and its downstream substrates as molecular targets of immunomodulatory drugs. Int J Hematol. 2016; 104(3): 293-9.,1010. Singhal S, Mehta J, Desikan R, Ayers D, Roberson P, Eddlemon P, et al. Antitumor activity of thalidomide in refractory multiple myeloma. N Engl J Med. 1999; 341(21): 1565-71. and later for other conditions.

The effectiveness of thalidomide in ENL is initially due to its action on tumour necrosis factor (TNF)-α, but other mechanisms may contribute to its anti-inflammatory effect.1111. Moreira AL, Sampaio EP, Zmuidzinas A, Frindt P, Smith KA, Kaplan G. Thalidomide exerts its inhibitory action on tumor necrosis factor alpha by enhancing mRNA degradation. J Exp Med. 1993; 177(6): 1675-80.,1212. Knobloch J, Jungck D, Koch A. The molecular mechanisms of thalidomide teratogenicity and implications for modern medicine. Curr Mol Med. 2017; 17(2): 108-17. Cereblon protein (CRBN) was described as the primary target of thalidomide, being involved also with the teratogenic potential of the drug.1313. Ito T, Ando H, Suzuki T, et al. Identification of a primary target of thalidomide teratogenicity. Science. 2010; 327(5971): 1345-50. CRBN is part of an E3-ubiquitin ligase complex (CRL4CRBN), acting as a substrate receptor that recognises specific targets for ubiquitination, leading to further degradation by the ubiquitin-proteasome system.1414. Yang J, Huang M, Zhou L, He X, Jiang X, Zhang Y, et al. Cereblon suppresses the lipopolysaccharide-induced inflammatory response by promoting the ubiquitination and degradation of c-Jun. J Biol Chem. 2018; 293(26): 10141-57. More recently, it was shown that CRBN is necessary for the efficacy of thalidomide and its analogs lenalidomide and pomalidomide [named immunomodulatory drugs (IMiDs)] in the treatment of MM.1515. Zhu YX, Braggio E, Shi C-X, Bruins LA, Schmidt JE, van Wier S, et al. Cereblon expression is required for the antimyeloma activity of lenalidomide and pomalidomide. Blood 2011; 118(18): 4771-9.

16. Lopez-Girona A, Mendy D, Ito T, Miller K, Gandhi AK, Kang J, et al. Cereblon is a direct protein target for immunomodulatory and antiproliferative activities of lenalidomide and pomalidomide. Leukemia. 2012; 26(11): 2326-35.

17. Huang SY, Lin CW, Lin HH, et al. Expression of cereblon protein assessed by immunohistochemical staining in myeloma cells is associated with superior response of thalidomide- and lenalidomide-based treatment, but not bortezomib-based treatment, in patients with multiple myeloma. Ann Hematol. 2014; 93(8): 1371-80.

18. Gao S, Wang S, Fan R, Hu J. Recent advances in the molecular mechanism of thalidomide teratogenicity. Biomed Pharmacother. 2020; 127: 110114.
-1919. Huang PA, Beedie SL, Chau CH, Venzon DJ, Gere S, Kazandjian D, et al. Cereblon gene variants and clinical outcome in multiple myeloma patients treated with lenalidomide. Sci Rep. 2019; 9(1): 1-7. Several studies have investigated the thalidomide-CRBN interaction with regard to the teratogenic, immunomodulatory and therapeutic effects of the drug. However, there are no studies on the effect of the thalidomide-CRBN interaction on ENL.

CRBN gene encodes a protein 442 amino acids long, has 11 exons (NM_016302) and is highly conserved; hence, polymorphisms in its coding regions are rare. Therefore, some studies have analysed non-coding regions of the gene that may be associated with the control of gene expression.2020. Szudy-Szczyrek A, Mlak R, Szczyrek M, Chocholska S, Sompor J, Nogalski A, et al. Polymorphisms in the promoter region of the CRBN gene as a predictive factor for the first-line CTD therapy in multiple myeloma patients. Oncotarget. 2018; 9(35): 24054-68.

21. Butrym A, Lacina P, Rybka J, Chaszczewska-Markowska M, Mazur G, Bogunia-Kubik K. Cereblon and IRF4 variants affect risk and response to treatment in multiple myeloma. Arch Immunol Ther Exp (Warsz). 2016; 64(Suppl 1): 151-6.
-2222. Mlak R, Szudy-Szczyrek A, Mazurek M, Szczyrek M, Homa-Mlak I, Mielnik M, et al. Polymorphisms in the promotor region of the CRBN gene as a predictive factor for peripheral neuropathy in the course of thalidomide-based chemotherapy in multiple myeloma patients. Br J Haematol. 2019; 186(5): 695-705. Three variants that flank the exons encoding the thalidomide-binding region of CRBN have been identified by bioinformatics tools as possible modulators of splicing sites with the potential to affect either the expression or activity of the protein.2323. Vianna FSL, Kowalski TW, Tovo-Rodrigues L, Tagliani-Ribeiro A, Godoy BA, Fraga LR, et al. Genomic and in silico analyses of CRBN gene and thalidomide embryopathy in humans. Reprod Toxicol. 2016; 66: 99-106. They are located in regions adjacent to the region encoding the portion of the protein that binds thalidomide, an A>C substitution (rs1620675) and a G>A substitution (rs1672770), both in intron 10, and an insertion/deletion (INS/DEL) of four nucleotides (-/GTTA) in the 3’UTR region adjacent to exon 11 (rs4183). Considering the role of CRBN on thalidomide actions and that these variants might affect the function, activity or expression of the protein, we evaluated the influence of these three variants in CRBN on the variation of the dose of thalidomide and the occurrence of adverse effects on the treatment of ENL.

SUBJECTS AND METHODS

Sample - The sample consisted of 103 ENL patients who were selected from National Reference Center of Sanitary Dermatology Dona Libânia in Fortaleza (state of Ceará, Brazil), Humanized Reference Center of Sanitary Dermatology in Imperatriz and Aquiles Lisboa Hospital in São Luís (Maranhão State) in northeast Brazil, and from the Dermatology Ambulatory of the University of São Paulo in Monte Negro (state of Rondônia, Brazil) in north Brazil.

The patients used thalidomide at different doses and had a follow-up of up to six visits (average of 3.6 months). Data from up to six consultations annotated in the patient’s medical record were analysed with the collection of clinical and demographic information, including sex, age and region of origin, history of leprosy (moment of diagnosis and treatment used) and history of ENL (diagnosis, treatment, adverse effects, history of relapse and dose of medications used).

Genetic analyses - DNA was extracted from saliva samples using the Oragene DNA Extraction Kit (DNA Genotek, Ottawa, Canada), according to the manufacturer’s instructions. A pair of primer was designed to amplify a fragment of 682 base pairs containing the region encompassing the three studied CRBN variants: forward 5’-TGTGGTCTTGGCAACCAGCAATTT-3’ and reverse 5’-ACTGCCGTTCATGCTTGTTTCCT-3’. This region was amplified by polymerase chain reaction (PCR). The fragment obtained was visualised on a 2% agarose gel, purified and sequenced using the same primers.

Sequences were visualised and analysed using CodonCodeAligner®, version 3.0.1 (CodonCode Corporation, Dedham, USA). The hg19 sequence deposited in GenBank was used as the reference sequence. When there was doubt about the variant, sequencing was repeated for confirmation.

Statistical analyses - Chi-square test was used to evaluate Hardy-Weinberg equilibrium for all polymorphisms. Generalised estimating equations method (GEE) was used to evaluate the influence of CRBN variants on thalidomide dose. This method is a repeated measures analysis focused on average changes in response over time and on the impact of covariates on these changes. GEE can model the average response of variables as a linear function of covariates of interest through a transformation or link function and can be used in studies where the data is asymmetric or the data distribution is difficult to verify due to the small-size sample.2424. Liang KY, Zeger SL. Longitudinal data analysis using generalized linear models. Biometrika. 1986; 73(1): 13-22.,2525. Sortica VA, Lindenau JD, Cunha MG, Ohnishi M do, Ventura AMR, Ribeiro-Dos-Santos ÂK, et al. The effect of SNPs in CYP450 in chloroquine/primaquine Plasmodium vivax malaria treatment. Pharmacogenomics. 2016; 17(17): 1903-11. The covariates inserted in the model were place of origin of the patient, concomitant use of multidrug therapy (MDT) for leprosy and use of other medications and other treatments for ENL.

The evaluation of the association of CRBN variants and haplotypes in the occurrence of adverse effects due to thalidomide treatment was based on clinical data using logistic regression. Models with and without correction by gender were used. The INS/C/T (rs4183/rs1672770/rs1620675) haplotype was removed from the association analyses because it presented few events and disturbed the analyses. Peripheral polyneuropathy, although is an adverse effect common to the use of thalidomide, has not been evaluated because it is difficult to distinguish it from polyneuropathy caused by ENL and leprosy. All statistical analyses were performed with SPSS version 20 (www.spss.com).

MLocus tool was used to calculate linkage disequilibrium (LD) for the variants,2626. Long JC, Williams RC, Urbanek M. An E-M algorithm and testing strategy for multiple-locus haplotypes. Am J Hum Genet. 1995; 56(3): 799-810. and haplotypes were inferred using the Bayesian algorithm of the phase 2.1.1 program.2727. Stephens M, Donnelly P. A comparison of Bayesian methods for haplotype reconstruction from population genotype data. Am J Hum Genet. 2003; 73(5): 1162-9.,2828. Stephens M, Smith NJ, Donnelly P. A new statistical method for haplotype reconstruction from population data. Am J Hum Genet. 2001; 68(4): 978-89.

Data availability - The data analysed during the current study are not publicly available to maintain patient confidentiality. Moreover, this type of request has not been previously approved by participants nor the human research committee. This data could, however, be available (anonymously) from the corresponding author on reasonable request.

Ethics - All participants were informed about the research objectives and signed an informed consent form. This study was approved by the Ethics Committee of the Hospital de Clínicas of Porto Alegre under number CAAE 21184413.0.0000.5327. All research was performed in accordance with Brazilian regulations and informed consent was obtained from all participants. This research has been performed in accordance with the Declaration of Helsinki.

RESULTS

In total, 103 ENL patients were included, being 82 (79.6%) male and 79 (76.7%) presenting LL (Table I). Forty-six (44.7%) patients were using MDT and two of these used an alternative MDT regimen. The mean time of treatment with thalidomide was 226 days.

TABLE I
Clinical and demographic characteristics of ENL patients (n = 103)

The genotypic distributions were in Hardy-Weinberg equilibrium and the allelic and genotypic frequencies of the polymorphisms are shown in Table II. A LD [rs4183/rs1672770 (r2 = 0.513), rs4183/rs1620675 (r2 = 0.943), rs1672770/rs1620675 (r2 = 0.567)] between the variants studied was identified (see Supplementary Table) and four haplotypes were identified in the sample (Table III). For all single nucleotide polymorphisms (SNPs), there was no association of the thalidomide dose with time of treatment (Table IV). In addition, no association was found between haplotypes and thalidomide dose.

TABLE II
Genotypic and allelic frequencies of CRBN polymorphisms in ENL patients on treatment with thalidomide
TABLE III
Haplotype frequencies of CRBN
TABLE IV
Analysis of interaction between CRBN genotype and time related to estimated thalidomide dose using the GEE model in ENL treatment

In regard to the relationship of CRBN variants and adverse effects of the use of thalidomide, it was found associations between the genotypes of rs1672770 (p = 0.040) and gastrointestinal effects, which include diarrhoea, vomiting, nausea, constipation and inappetence (Tables V and VI). Haplotype analysis was performed with all SNPs because only one combination of SNPs showed r2>0.8 (rs4183/rs1620675 with r2 = 0.943). The three polymorphisms were maintained in the analysis because when combining rs4183/rs1672770 or rs1672770/rs1620675, it showed r2<0.8. The analysis of the association between haplotypes and adverse effects also showed an association of gastrointestinal adverse effects (p = 0.015) with haplotype DEL/C/T (rs4183/rs1672770/rs1620675) (Tables VII and VIII). These effects were present in eight patients (9.6%) who carried this haplotype.

TABLE V
Analysis of the association of CRBN variants on the occurrence of gastrointestinal adverse effectsa
TABLE VI
Frequency of gastrointestinal adverse effectsa according to CRBN genotypes
TABLE VII
Analysis of the association of CRBN haplotypes on the occurrence of gastrointestinal adverse effectsa
TABLE VIII
Frequency of gastrointestinal adverse effectsa according to CRBN haplotypes

DISCUSSION

This study aimed to identify genetic variants in CRBN gene that might influence in response to the treatment of ENL with thalidomide. We identified that the SNP rs1672770 and the haplotype DEL/C/T (rs4183/rs1672770/rs1620675) were associated with the manifestation of adverse gastrointestinal effects.

CRBN acts as a substrate receptor as part of the E3 ubiquitin ligase complex (CRL4CRBN), which controls the expression of target proteins by their ubiquitination and degradation. It is necessary for the teratogenic effect of thalidomide and also important for the antiproliferative effect of thalidomide and other IMiDs in MM.2929. Gandhi AK, Kang J, Havens CG, Conklin T, Ning Y, Wu L, et al. Immunomodulatory agents lenalidomide and pomalidomide co-stimulate T cells by inducing degradation of T cell repressors Ikaros and Aiolos via modulation of the E3 ubiquitin ligase complex CRL4CRBN. Br J Haematol. 2014; 164(6): 811-21. It is postulated that when thalidomide binds to CRBN, it modifies its function causing teratogenic effects by preventing the degradation proteins and/or by creating neosubtrates for ubiquitination and proteasomal degradation, which play a crucial role in embryonic development.3030. Ito T, Ando H, Handa H. Teratogenic effects of thalidomide: Molecular mechanisms. Cell Mol Life Sci. 2011; 68(9): 1569-79.,3131. Schuster SR, Kortuem KM, Zhu YX, Braggio E, Shi CX, Bruins LA, et al. The clinical significance of cereblon expression in multiple myeloma. Leuk Res. 2014; 38(1): 23-8. In the case of MM, thalidomide binding to CRBN promotes recruitment of the neosubstrates Ikaros (IKZF1) and Aiolos (IKZF3) to the ubiquitin-ligase complex, resulting in increased ubiquitination and degradation of these transcription factors in T cells and MM cells.3131. Schuster SR, Kortuem KM, Zhu YX, Braggio E, Shi CX, Bruins LA, et al. The clinical significance of cereblon expression in multiple myeloma. Leuk Res. 2014; 38(1): 23-8.

32. Dimopoulos K, Fibiger Munch-Petersen H, Winther Eskelund C, Dissing Sjö L, Ralfkiaer E, Gimsing P, et al. Expression of CRBN, IKZF1, and IKZF3 does not predict lenalidomide sensitivity and mutations in the cereblon pathway are infrequent in multiple myeloma. Leuk Lymphoma. 2019; 60(1): 180-8.
-3333. Chamberlain PP, Lopez-Girona A, Miller K, Carmel G, Pagarigan B, Chie-Leon B, et al. Structure of the human cereblon-DDB1-lenalidomide complex reveals basis for responsiveness to thalidomide analogs. Nat Struct Mol Biol. 2014; 21(9): 803-9. In addition, some studies have associated low CRBN mRNA expression with poorer clinical response to IMiDs, suggesting a potential role of CRBN as a predictive biomarker for treatment response.1717. Huang SY, Lin CW, Lin HH, et al. Expression of cereblon protein assessed by immunohistochemical staining in myeloma cells is associated with superior response of thalidomide- and lenalidomide-based treatment, but not bortezomib-based treatment, in patients with multiple myeloma. Ann Hematol. 2014; 93(8): 1371-80.,3131. Schuster SR, Kortuem KM, Zhu YX, Braggio E, Shi CX, Bruins LA, et al. The clinical significance of cereblon expression in multiple myeloma. Leuk Res. 2014; 38(1): 23-8.,3232. Dimopoulos K, Fibiger Munch-Petersen H, Winther Eskelund C, Dissing Sjö L, Ralfkiaer E, Gimsing P, et al. Expression of CRBN, IKZF1, and IKZF3 does not predict lenalidomide sensitivity and mutations in the cereblon pathway are infrequent in multiple myeloma. Leuk Lymphoma. 2019; 60(1): 180-8.,3434. Heintel D, Rocci A, Ludwig H, Bolomsky A, Caltagirone S, Schreder M, et al. High expression of cereblon (CRBN) is associated with improved clinical response in patients with multiple myeloma treated with lenalidomide and dexamethasone. Br J Haematol. 2013; 161(5): 695-700. Since the 1990s, it is known that one of the main effects of thalidomide is to decrease the TNF mRNA half-life, which explains some of its therapeutic effects. Using CRBN knockdown, it has been shown that the inhibitory effect of IMiDs on TNF-α production was also impaired in the silencing of CRBN.1616. Lopez-Girona A, Mendy D, Ito T, Miller K, Gandhi AK, Kang J, et al. Cereblon is a direct protein target for immunomodulatory and antiproliferative activities of lenalidomide and pomalidomide. Leukemia. 2012; 26(11): 2326-35.,3535. Zhu YX, Kortuem KM, Stewart AK. Molecular mechanism of action of immune-modulatory drugs thalidomide, lenalidomide and pomalidomide in multiple myeloma. Leuk Lymphoma. 2013; 54(4): 683-7.

CRBN is composed of 11 exons extending over 30 Kb. Thalidomide binds to CRBN in a region of 104 amino acids (339-442) located in the C-terminal portion, encoded by exons 9, 10 and 11.1313. Ito T, Ando H, Suzuki T, et al. Identification of a primary target of thalidomide teratogenicity. Science. 2010; 327(5971): 1345-50. This gene is extremely conserved and a few polymorphisms was found in the coding region.2323. Vianna FSL, Kowalski TW, Tovo-Rodrigues L, Tagliani-Ribeiro A, Godoy BA, Fraga LR, et al. Genomic and in silico analyses of CRBN gene and thalidomide embryopathy in humans. Reprod Toxicol. 2016; 66: 99-106.,3636. Butrym A, Rybka J, Lacina P, Gebura K, Frontkiewicz D, Bogunia-Kubik K, et al. Polymorphisms within beta-catenin encoding gene affect multiple myeloma development and treatment. Leuk Res. 2015; 39(12): 1462-6. Studies performed with MM showed that variants in non-coding regions of CRBN were associated with response to thalidomide and others IMIDs therapy.2020. Szudy-Szczyrek A, Mlak R, Szczyrek M, Chocholska S, Sompor J, Nogalski A, et al. Polymorphisms in the promoter region of the CRBN gene as a predictive factor for the first-line CTD therapy in multiple myeloma patients. Oncotarget. 2018; 9(35): 24054-68.,2121. Butrym A, Lacina P, Rybka J, Chaszczewska-Markowska M, Mazur G, Bogunia-Kubik K. Cereblon and IRF4 variants affect risk and response to treatment in multiple myeloma. Arch Immunol Ther Exp (Warsz). 2016; 64(Suppl 1): 151-6.,3737. Iskierka-Jazdzewska E, Canzian F, Stepien A, Martino A, Campa D, Stein A, et al. Cereblon (CRBN) gene polymorphisms predict clinical response and progression-free survival in relapsed/refractory multiple myeloma patients treated with lenalidomide: a pharmacogenetic study from the IMMEnSE consortium. Leuk Lymphoma. 2020; 61(3): 699-706.

The genetic variants evaluated in the present study were identified in a previous study as possible splicing sites.2323. Vianna FSL, Kowalski TW, Tovo-Rodrigues L, Tagliani-Ribeiro A, Godoy BA, Fraga LR, et al. Genomic and in silico analyses of CRBN gene and thalidomide embryopathy in humans. Reprod Toxicol. 2016; 66: 99-106. The association of these SNPs and the dose of thalidomide could indicate that polymorphisms in these regions could interfere in the expression of the CRBN gene or in the activity of CRBN, being able to modulate the response to treatment with thalidomide. However, in this study, we were unable to identify an association between these variants and the dose variation of thalidomide over time.

We found an association between the genetic variant rs1672770 (p = 0.040) and the manifestation of gastrointestinal adverse effects. There was also an association of gastrointestinal adverse effects (p = 0.015) with haplotype DEL/C/T (rs4183/rs1672770/rs1620675). These adverse gastrointestinal effects consisted of diarrhoea, vomiting, nausea and constipation. The most common symptom was constipation, a commonly reported side effect of thalidomide.2222. Mlak R, Szudy-Szczyrek A, Mazurek M, Szczyrek M, Homa-Mlak I, Mielnik M, et al. Polymorphisms in the promotor region of the CRBN gene as a predictive factor for peripheral neuropathy in the course of thalidomide-based chemotherapy in multiple myeloma patients. Br J Haematol. 2019; 186(5): 695-705.,3838. Walker SL, Waters MFR, Lockwood DNJ. The role of thalidomide in the management of erythema nodosum leprosum. Lepr Rev. 2007; 78(3): 197-215.,3939. Wines NY, Cooper AJ, Wines MP. Thalidomide in dermatology. Aust J Dermatol. 2002; 43(4): 229-40. The association of adverse effects with the polymorphisms studied may also be related to differences in CRBN expression. Mlak et al.2222. Mlak R, Szudy-Szczyrek A, Mazurek M, Szczyrek M, Homa-Mlak I, Mielnik M, et al. Polymorphisms in the promotor region of the CRBN gene as a predictive factor for peripheral neuropathy in the course of thalidomide-based chemotherapy in multiple myeloma patients. Br J Haematol. 2019; 186(5): 695-705. found an association between CRBN variants (rs6768972, rs16727) and the risk of polyneuropathy and gastrointestinal disorders in patients with MM treated with thalidomide-based regimens. These variants are located in the CRBN promoter region. Thus, they could influence the expression or activity of CRBN.2323. Vianna FSL, Kowalski TW, Tovo-Rodrigues L, Tagliani-Ribeiro A, Godoy BA, Fraga LR, et al. Genomic and in silico analyses of CRBN gene and thalidomide embryopathy in humans. Reprod Toxicol. 2016; 66: 99-106. In the haplotype analysis, only one combination of SNPs showed r2>0.8 (rs4183/rs1620675 with r2 = 0.943), as can be observed in the Supplementary Table. Taking into account this result, our small sample size, and that this is the first study to evaluate CRBN haplotypes, we decided to maintain the logistic regression analysis considering all SNPs. Further studies should be performed in order to elucidate and confirm the haplotypes and associations found here.

In this study, 46 patients were using MDT during the treatment for ENL and only two patients received an alternative MDT regimen. For this reason, they were not divided into a subgroup during the analyses. In addition, the concomitant medications that patients used during treatment for ENL can be diverse. That is why they were not listed in a table. MDT and other drugs were used as covariates in the GEE analysis to avoid bias in the use of these drugs. However, it was not possible to use this correction in the other analyses. The epidemiological data of this study were in agreement with the data found in other studies.22. 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(5): 868-79.,4040. Feuth M, Brandsma JW, Faber WR, Bhattarai B, Feuth T, Anderson AM. Erythema nodosum leprosum in Nepal: a retrospective study of clinical features and response to treatment with prednisolone or thalidomide. Lepr Rev. 2008; 79(3). Most of the individuals with ENL had LL and were male. These frequencies corroborate that men may be more affected by LL and that the bacillary index is a risk factor for the development of the reaction as demonstrated in other studies.22. 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(5): 868-79.,4141. 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(9): e0004065.,4242. 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(5): 384-90. In addition, many patients were on MDT during treatment for ENL, confirming that the reaction manifests mainly in the first year of illness, during MDT.4343. Voorend CG, Post EB. A systematic review on the epidemiological data of erythema nodosum leprosum, a type 2 leprosy reaction. PLoS Negl Trop Dis. 2013; 7(10): e2440.,4444. 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(5): 1905-13.

Some limitations should be considered during the interpretation of our study. It was not possible to obtain samples of skin lesions from these patients. Thus, the DNA samples analysed in this work were obtained from saliva because it presents a methodology that is easier to obtain. However, this does not influence the results presented because it is about the analysis of genetic variants and not the gene expression. Clearly, further studies could focus on the impact of the CRBN expression in the skin lesions before and after the treatment. In addition, the use of clofazimine and in MDT for leprosy, because of its anti-inflammatory effect, and the other drugs may interfere with both the dose reduction of drugs and in the manifestation of adverse effects.22. 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(5): 868-79.,33. 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 Tropica. 2018; 183: 134-41. This is also a retrospective study carried out with clinical data of patients, thalidomide dose, and manifestation of adverse effects obtained through the analysis of medical records of up to six consultations. The lack of standardization in the description of this information or even the absence of a description of the adverse effects cannot be ruled out. Furthermore, as ENL is a chronic condition and its treatment is very long, in this study, the mean treatment time during the collection period was 226 days, but some patients had been in treatment for years and remained on thalidomide treatment after the data collection period. Thus, in some cases, it was not possible to completely monitor the treatment of patients. Peripheral neuropathy, one of the main adverse effects of thalidomide,33. 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 Tropica. 2018; 183: 134-41.,44. Melchert M, List A. The thalidomide saga. Int J Biochem Cell Biol. 2007; 39(7-8): 1489-99. was not evaluated in this study due to its retrospective nature, making it difficult to differentiate neuropathy due to the disease itself or due to the use of thalidomide. Accordingly, the evaluation of the association of variants on the onset of peripheral neuropathy should be performed by a prospective study to reduce potential biases. Another limitation was that this study was carried out with individuals from different regions of the country and the heterogeneous genetic background of the Brazilian population might have been underestimated in the analysis and interpretation of the results. In addition, the limited sample size used in this study can become difficult the identification of minor effects from the genetic variants in different outcomes.

On basis of the results of this study, we identified some genetic variants of CRBN were associated with adverse gastrointestinal effects. These results indicate that such variants could impact the protein and influence the outcome of treatment with thalidomide. This also indicates that CRBN may also be necessary for the action of thalidomide in ENL, as already described in MM.2929. Gandhi AK, Kang J, Havens CG, Conklin T, Ning Y, Wu L, et al. Immunomodulatory agents lenalidomide and pomalidomide co-stimulate T cells by inducing degradation of T cell repressors Ikaros and Aiolos via modulation of the E3 ubiquitin ligase complex CRL4CRBN. Br J Haematol. 2014; 164(6): 811-21. Clearly, more studies evaluating the impact the variants in CRBN, as well as associations with ENL treatment, must be performed in order to confirm this hypothesis. ENL is a chronic and difficult-to-control condition. Thalidomide is an effective drug, but has restrictions on its use due to peripheral neuropathy and its teratogenicity.77. Mcbride WG. Thalidomide and congenital abnormalities. The Lancet. 1961; 278(7216): 1358.,88. Lenz W, Pfeiffer RA, Kosenow W, Hayman DJ. Thalidomide and congenital abnormalities. The Lancet. 1962; 279(7219): 45-6. Therefore, it is important to identify useful biomarkers in predicting treatment response to limit its use to patients who will benefit most from treatment. To our knowledge, this is the first study to evaluate the association of genetic variants of CRBN with thalidomide treatment in ENL patients. There are still many gaps to be filled in our knowledge of the mechanism of action of thalidomide and on how CRBN participates in this process. Thus, this study shows that evaluation of CRBN and its expression may help to understand the action of thalidomide in ENL and perhaps, in the future, be a useful biomarker in ENL.

ACKNOWLEDGEMENTS

To the collaborators who helped in obtaining the specimens at all the participating centres, and the statistician Luciano Santos Pinto Guimarães, for his assistance with statistical analyses. FSLV is the recipient of a CNPq scholarship grant (no. 312960/2021-2).

REFERENCES

  • 1
    Chandler DJ, Hansen KS, Mahato B, Darlong J, John A, Lockwood DN. Household costs of leprosy reactions (ENL) in rural India. PLoS Negl Trop Dis. 2015; 9(1): e0003431.
  • 2
    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(5): 868-79.
  • 3
    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 Tropica. 2018; 183: 134-41.
  • 4
    Melchert M, List A. The thalidomide saga. Int J Biochem Cell Biol. 2007; 39(7-8): 1489-99.
  • 5
    Sheskin J. Thalidomide in the treatment of lepra reactions. Clin Pharmacol Ther. 1965; 6(3): 303-6.
  • 6
    Mercurio A, Adriani G, Catalano A, Carocci A, Rao L, Lentini G, et al. A Mini-review on thalidomide: chemistry, mechanisms of action, therapeutic potential and anti-angiogenic properties in multiple myeloma. Curr Med Chem. 2017; 24: 2736-44.
  • 7
    Mcbride WG. Thalidomide and congenital abnormalities. The Lancet. 1961; 278(7216): 1358.
  • 8
    Lenz W, Pfeiffer RA, Kosenow W, Hayman DJ. Thalidomide and congenital abnormalities. The Lancet. 1962; 279(7219): 45-6.
  • 9
    Ito T, Handa H. Cereblon and its downstream substrates as molecular targets of immunomodulatory drugs. Int J Hematol. 2016; 104(3): 293-9.
  • 10
    Singhal S, Mehta J, Desikan R, Ayers D, Roberson P, Eddlemon P, et al. Antitumor activity of thalidomide in refractory multiple myeloma. N Engl J Med. 1999; 341(21): 1565-71.
  • 11
    Moreira AL, Sampaio EP, Zmuidzinas A, Frindt P, Smith KA, Kaplan G. Thalidomide exerts its inhibitory action on tumor necrosis factor alpha by enhancing mRNA degradation. J Exp Med. 1993; 177(6): 1675-80.
  • 12
    Knobloch J, Jungck D, Koch A. The molecular mechanisms of thalidomide teratogenicity and implications for modern medicine. Curr Mol Med. 2017; 17(2): 108-17.
  • 13
    Ito T, Ando H, Suzuki T, et al. Identification of a primary target of thalidomide teratogenicity. Science. 2010; 327(5971): 1345-50.
  • 14
    Yang J, Huang M, Zhou L, He X, Jiang X, Zhang Y, et al. Cereblon suppresses the lipopolysaccharide-induced inflammatory response by promoting the ubiquitination and degradation of c-Jun. J Biol Chem. 2018; 293(26): 10141-57.
  • 15
    Zhu YX, Braggio E, Shi C-X, Bruins LA, Schmidt JE, van Wier S, et al. Cereblon expression is required for the antimyeloma activity of lenalidomide and pomalidomide. Blood 2011; 118(18): 4771-9.
  • 16
    Lopez-Girona A, Mendy D, Ito T, Miller K, Gandhi AK, Kang J, et al. Cereblon is a direct protein target for immunomodulatory and antiproliferative activities of lenalidomide and pomalidomide. Leukemia. 2012; 26(11): 2326-35.
  • 17
    Huang SY, Lin CW, Lin HH, et al. Expression of cereblon protein assessed by immunohistochemical staining in myeloma cells is associated with superior response of thalidomide- and lenalidomide-based treatment, but not bortezomib-based treatment, in patients with multiple myeloma. Ann Hematol. 2014; 93(8): 1371-80.
  • 18
    Gao S, Wang S, Fan R, Hu J. Recent advances in the molecular mechanism of thalidomide teratogenicity. Biomed Pharmacother. 2020; 127: 110114.
  • 19
    Huang PA, Beedie SL, Chau CH, Venzon DJ, Gere S, Kazandjian D, et al. Cereblon gene variants and clinical outcome in multiple myeloma patients treated with lenalidomide. Sci Rep. 2019; 9(1): 1-7.
  • 20
    Szudy-Szczyrek A, Mlak R, Szczyrek M, Chocholska S, Sompor J, Nogalski A, et al. Polymorphisms in the promoter region of the CRBN gene as a predictive factor for the first-line CTD therapy in multiple myeloma patients. Oncotarget. 2018; 9(35): 24054-68.
  • 21
    Butrym A, Lacina P, Rybka J, Chaszczewska-Markowska M, Mazur G, Bogunia-Kubik K. Cereblon and IRF4 variants affect risk and response to treatment in multiple myeloma. Arch Immunol Ther Exp (Warsz). 2016; 64(Suppl 1): 151-6.
  • 22
    Mlak R, Szudy-Szczyrek A, Mazurek M, Szczyrek M, Homa-Mlak I, Mielnik M, et al. Polymorphisms in the promotor region of the CRBN gene as a predictive factor for peripheral neuropathy in the course of thalidomide-based chemotherapy in multiple myeloma patients. Br J Haematol. 2019; 186(5): 695-705.
  • 23
    Vianna FSL, Kowalski TW, Tovo-Rodrigues L, Tagliani-Ribeiro A, Godoy BA, Fraga LR, et al. Genomic and in silico analyses of CRBN gene and thalidomide embryopathy in humans. Reprod Toxicol. 2016; 66: 99-106.
  • 24
    Liang KY, Zeger SL. Longitudinal data analysis using generalized linear models. Biometrika. 1986; 73(1): 13-22.
  • 25
    Sortica VA, Lindenau JD, Cunha MG, Ohnishi M do, Ventura AMR, Ribeiro-Dos-Santos ÂK, et al. The effect of SNPs in CYP450 in chloroquine/primaquine Plasmodium vivax malaria treatment. Pharmacogenomics. 2016; 17(17): 1903-11.
  • 26
    Long JC, Williams RC, Urbanek M. An E-M algorithm and testing strategy for multiple-locus haplotypes. Am J Hum Genet. 1995; 56(3): 799-810.
  • 27
    Stephens M, Donnelly P. A comparison of Bayesian methods for haplotype reconstruction from population genotype data. Am J Hum Genet. 2003; 73(5): 1162-9.
  • 28
    Stephens M, Smith NJ, Donnelly P. A new statistical method for haplotype reconstruction from population data. Am J Hum Genet. 2001; 68(4): 978-89.
  • 29
    Gandhi AK, Kang J, Havens CG, Conklin T, Ning Y, Wu L, et al. Immunomodulatory agents lenalidomide and pomalidomide co-stimulate T cells by inducing degradation of T cell repressors Ikaros and Aiolos via modulation of the E3 ubiquitin ligase complex CRL4CRBN. Br J Haematol. 2014; 164(6): 811-21.
  • 30
    Ito T, Ando H, Handa H. Teratogenic effects of thalidomide: Molecular mechanisms. Cell Mol Life Sci. 2011; 68(9): 1569-79.
  • 31
    Schuster SR, Kortuem KM, Zhu YX, Braggio E, Shi CX, Bruins LA, et al. The clinical significance of cereblon expression in multiple myeloma. Leuk Res. 2014; 38(1): 23-8.
  • 32
    Dimopoulos K, Fibiger Munch-Petersen H, Winther Eskelund C, Dissing Sjö L, Ralfkiaer E, Gimsing P, et al. Expression of CRBN, IKZF1, and IKZF3 does not predict lenalidomide sensitivity and mutations in the cereblon pathway are infrequent in multiple myeloma. Leuk Lymphoma. 2019; 60(1): 180-8.
  • 33
    Chamberlain PP, Lopez-Girona A, Miller K, Carmel G, Pagarigan B, Chie-Leon B, et al. Structure of the human cereblon-DDB1-lenalidomide complex reveals basis for responsiveness to thalidomide analogs. Nat Struct Mol Biol. 2014; 21(9): 803-9.
  • 34
    Heintel D, Rocci A, Ludwig H, Bolomsky A, Caltagirone S, Schreder M, et al. High expression of cereblon (CRBN) is associated with improved clinical response in patients with multiple myeloma treated with lenalidomide and dexamethasone. Br J Haematol. 2013; 161(5): 695-700.
  • 35
    Zhu YX, Kortuem KM, Stewart AK. Molecular mechanism of action of immune-modulatory drugs thalidomide, lenalidomide and pomalidomide in multiple myeloma. Leuk Lymphoma. 2013; 54(4): 683-7.
  • 36
    Butrym A, Rybka J, Lacina P, Gebura K, Frontkiewicz D, Bogunia-Kubik K, et al. Polymorphisms within beta-catenin encoding gene affect multiple myeloma development and treatment. Leuk Res. 2015; 39(12): 1462-6.
  • 37
    Iskierka-Jazdzewska E, Canzian F, Stepien A, Martino A, Campa D, Stein A, et al. Cereblon (CRBN) gene polymorphisms predict clinical response and progression-free survival in relapsed/refractory multiple myeloma patients treated with lenalidomide: a pharmacogenetic study from the IMMEnSE consortium. Leuk Lymphoma. 2020; 61(3): 699-706.
  • 38
    Walker SL, Waters MFR, Lockwood DNJ. The role of thalidomide in the management of erythema nodosum leprosum. Lepr Rev. 2007; 78(3): 197-215.
  • 39
    Wines NY, Cooper AJ, Wines MP. Thalidomide in dermatology. Aust J Dermatol. 2002; 43(4): 229-40.
  • 40
    Feuth M, Brandsma JW, Faber WR, Bhattarai B, Feuth T, Anderson AM. Erythema nodosum leprosum in Nepal: a retrospective study of clinical features and response to treatment with prednisolone or thalidomide. Lepr Rev. 2008; 79(3).
  • 41
    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(9): e0004065.
  • 42
    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(5): 384-90.
  • 43
    Voorend CG, Post EB. A systematic review on the epidemiological data of erythema nodosum leprosum, a type 2 leprosy reaction. PLoS Negl Trop Dis. 2013; 7(10): e2440.
  • 44
    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(5): 1905-13.
  • Financial support: INAGEMP (CNPq 573993/2008-4, 465549/2014-4, FAPERGS 17/2551.0000521-0 and CAPES), FIPE/HCPA (GPPG 10-0410) and CAPES. FSLV is recipient of a CNPq scholarship grant (312993/2017-0) and FSLV is recipient of a FAPERGS grant (19/2551-0001787-1).

Publication Dates

  • Publication in this collection
    14 Nov 2022
  • Date of issue
    2022

History

  • Received
    10 Feb 2022
  • Accepted
    04 Oct 2022
Instituto Oswaldo Cruz, Ministério da Saúde Av. Brasil, 4365 - Pavilhão Mourisco, Manguinhos, 21040-900 Rio de Janeiro RJ Brazil, Tel.: (55 21) 2562-1222, Fax: (55 21) 2562 1220 - Rio de Janeiro - RJ - Brazil
E-mail: memorias@fiocruz.br