Skip to main content

Blood gene expression of Toll-like receptors in SLE patients with lupus nephritis or neuropsychiatric systemic lupus erythematosus

Abstract

Background

To determine differences in the blood innate gene expression signatures of systemic lupus erythematosus (SLE) patients across various organ manifestations and disease activity, with a focus on lupus nephritis (LN) and central nervous system (CNS) involvement.

Methods

Toll-like receptor family (TLR 1–10) mRNA expression was investigated in peripheral blood mononuclear cells from patients with SLE (n = 74) and healthy controls (n = 34). We compared patients with histologically confirmed active LN or neuropsychiatric systemic lupus erythematosus (NPSLE) with patients without these symptoms. The expression of TLR mRNA was determined by RT‒qPCR using a high-throughput SmartChip Real-Time-qPCR system (WaferGen). Multivariate analysis and nonparametric statistics were used for data analysis to assess the associations between TLRs and disease activity and severity.

Results

TLR4 (0.044 vs. 0.081, p = 0.012) was upregulated and TLR10 (0.009 vs. 0.006, p = 0.0007) was downregulated in the whole cohort of SLE patients compared to healthy controls. A comparison of the active LN group with participants without kidney involvement revealed increased expression of TLR2 (0.078 vs. 0.03, p = 0.009), and TLR5 (0.035 vs. 0.017, p = 0.03). Moreover, a significant difference was observed in TLR9 expression between inactive LN and the control group (0.014 vs. 0.009, p = 0.01), together with borderline correlation in TLR2 expression (0.04 vs. 0.03, p = 0.06). Receiver operating characteristic (ROC) curve analysis revealed that TLR1 and TLR2 expression were the best potential diagnostic markers for active LN. The NPSLE group showed upregulation of TLR1 (0.088 vs. 0.048, p = 0.01), TLR4 (0.173 vs. 0.066, p = 0.0003) and TLR6 (0.087 vs. 0.036, 0.007). Our correlation analysis supported the close relationships among the expression of individual TLRs in the whole lupus cohort and its subgroups.

Conclusion

Our study revealed differences in TLR expression between a lupus cohort and healthy controls. Additionally, our analysis provides insight into specific TLR expression in cases with severe organ manifestations, such as LN and NPSLE. The multiple mutual relationships of TLRs demonstrate the activation of innate immunity in SLE and suggest promising targets for future therapies or diagnostics.

Background

The pathophysiology of SLE is complex and involves interactions between factors and cells involved in innate and adaptive immunity. Innate immune receptors such as Toll-like receptors (TLRs) have the unique ability to recognise pathogen-associated molecular patterns (PAMPs), damage-associated molecular patterns (DAMPs) and self-molecules in autoimmune diseases, thus inducing activation of the immune system by linking innate and adaptive immune responses [1]. TLRs are classified according to their location. Extracellular TLRs (TLR1,2,4,5,6) respond to bacterial or fungal PAMPs and DAMPs, while intracellular TLRs (TLR3,7,8,9) are located on the endosomal membrane and respond to pathogen nucleic acids as well as self-nucleic acids [2]. They are usually expressed by cells of the immune system but can also be expressed by a specific organ tissue [3].

A number of studies have demonstrated that different members of the TLR family are involved in the pathogenesis of SLE [4, 5]. The expression of TLRs varies in different clinical manifestations of the disease [6]. Some TLRs (TLR7 and TLR8) are potential therapeutic targets for monoclonal antibodies [7].

The role of innate immunity in SLE development is already known. However, there is limited knowledge regarding the differences and similarities in TLR family members based on specific types of organ damage. This study focused on the possible differences between the expression of TLR family members in SLE patients with severe disease manifestations, such as active lupus nephritis (LN), and patients with neuropsychiatric systemic lupus erythematosus (NPSLE). This study analysed the gene expression of TLR family members in a clinically defined group of SLE patients to evaluate the potential of TLRs for the diagnosis of organ involvement.

Methods

Study subjects

The study was designed as an observational cross-sectional single-centre analysis of a cohort of patients from the Rheumatology Department of the University Hospital in Olomouc, Czech Republic. The primary objective of the study was to compare the expression of TLRs in patients with SLE and healthy controls. The secondary objectives were to investigate differences in TLR expression between patients with kidney or central nervous system (CNS) involvement and to compare these results with disease activity scores, major antibody levels and disease-related organ damage.

The depersonalised and anonymous data and samples were collected and stored after signed written informed consent was obtained for data collection from all participants and healthy controls. The study was approved by the local ethical committee (Multicenter Research Ethics Committee of the University Hospital Olomouc and Faculty of Medicine and Dentistry, Palacky University Olomouc, Czech Republic). The study was performed in accordance with the Helsinki Declaration of 1964 and its later amendments.

All SLE patients enrolled in the study met the 2019 EULAR/ACR classification criteria for SLE [8]. Disease activity was quantified using the Systemic Lupus Eryhematosus Disease Activity Index 2000 (SLEDAI-2 K) [9]. The degree of organ damage that has accumulated since disease onset was assessed using the Systemic Lupus International Collaborating Clinics American College of Rheumatology Damage Index (SLICC/ACR DI) [10].

Subgroups were created based on the clinical manifestations of the disease. The diagnosis of LN was confirmed by kidney biopsy. Patients were divided into two groups according to the activity of the LN, as assessed by the Renal Domain Score of the Systemic Lupus Erythematosus Disease Activity Index (SLEDAI-R). The groups were named the Active LN-SLE group, with SLEDAI-R scores ≥ 4, and the Inactive LN-SLE group. The SLEDAI-R was calculated by summing the SLEDAI-2 K item scores for proteinuria > 0.5 g/day, hematuria and pyuria (both > 5 cells/HPF, infection excluded) and cellular casts [11]. Patients without kidney impairment were classified into the No LN-SLE group. Another subgroup, called CNS-SLE, consisted of patients with NPSLE and CNS involvement according to the ACR criteria [12]. Patients with both active LN and NPSLE were excluded from the subgroup analysis. Furthermore, the expression of TLR was evaluated in relation to the administration of corticosteroids or antimalarials in the entire cohort of patients with SLE.

The control group consisted of age- and sex-matched healthy volunteers. In all healthy controls, recent vaccinations, infections, and the use of immunosuppressive drugs as well as the presence of inflammatory autoimmune diseases in first- or second-degree relatives were excluded by questionnaire.

Sample processing

Peripheral blood mononuclear cells (PBMCs) were isolated from peripheral blood collected in K3EDTA tubes by Ficoll density gradient centrifugation (Sigma‒Aldrich, Germany) and stored in Tri reagent (Sigma‒Aldrich, Germany) at − 80 °C until analysis. Total RNA was extracted using a Direct-zol RNA kit (Zymo Research, USA) according to the manufacturer’s recommendations. After reverse transcription with a Transcriptor First Strand cDNA Synthesis Kit (Roche, Switzerland), qPCR was performed in a 100-nl reaction volume containing LightCycler 480 SYBR Green I Master mix (Roche, Switzerland) using a high-throughput SmartChip Real-Time-qPCR system (WaferGen, USA), and the relative mRNA expression was calculated using phosphoglycerate kinase 1 as a reference gene as reported previously [13].

Statistical analysis

IBM SPSS Statistics for Windows, Version 23.0 Armonk, NY: IBM Corp. statistical software was used for statistical processing. Multivariate analysis (Andrews curves, Associate rules, Spearman correlation coefficient and Kruskal‒Wallis test) was used for data analysis. The Mann-Whitney U-test was used to compare TLR expression according to treatment. All tests were performed at the 0.05 level of significance. Differences for which the p value was less than 0.05 were considered to be statistically significant and are indicated by an asterisk. According to the receiver operating characteristic (ROC) curve analysis, areas under the ROC curve (AUCs) greater than 0.7 are considered to be of fair discriminating value for predicting diagnostic tests. Optimal cut-off values were established based on Youden’s J statistic [14, 15].

Results

Participants

A total of 74 white patients fulfilled the inclusion criteria and were enrolled in the study. The group of healthy controls consisted of 34 volunteers. The demographic and clinical characteristics are presented in the Table 1., while the current and previous medications, along with the main comorbidities, is presented in Additional file Table 4. The group of patients with active LN, hereafter referred to as the Active LN-SLE group, included 11 patients (mean age 35 years, min–max 19–57 years, female/male 10/1) and 14 patients (mean age 35 years, min–max 22–54 years, female/male 11/3) with histologically confirmed but inactive LN at the time of sampling were considered the Inactive LN-SLE group. The most prevalent class of LN among patients with SLE and renal impairment, as determined by kidney biopsy, was class IV in 44% of cases and class II in 28%. The distribution of the histological classes of LN in subgroups is presented in Additional file Table 5. A total of 47 patients with SLE (No LN-SLE group) did not have kidney impairment (mean age 43 years, min–max 20–74 years, female/male 43/4). No statistically significant difference has been identified between antibody levels in the aforementioned study subgroups (Additional file Table 6). Six patients were diagnosed with NPSLE (mean age 47 years, min–max 20–60 years, female/male 6/0). A total of nine patients (12% of the SLE cohort) were not currently undergoing antimalarial treatment. All of these patients were included in the No LN-SLE group. Furthermore, 22 patients (30% of the SLE cohort) were not currently receiving corticosteroid treatment. Of these, three were from the Inactive LN-SLE group (21% of the Inactive LN-SLE group). The remaining 19 patients (40% of the No LN-SLE group) were part of the No LN-SLE group.

TLR family mRNA expression in patients with SLE versus healthy controls

The relative mRNA expression of the TLR family members was compared between the SLE patients and healthy controls. The aim was to identify differences between these groups. Compared to healthy controls, the relative mRNA expression of TLR4 was upregulated in the SLE group (0.044 vs. 0.081, p = 0.012; Table 2A, Fig. 1). Conversely, TLR10 expression was significantly lower in the SLE group than in the healthy control group (0.009 vs. 0.006, p = 0.0007; Table 2A, Fig. 1). Furthermore, increased expression of TLR6 and TLR7 and decreased expression of TLR3 were observed, although the differences did not reach statistical significance.

TLR family expression in active LN-SLE, inactive LN-SLE and no LN-SLE patients

To investigate the clinical manifestation-specific innate immune gene expression pattern, we compared SLE patients without renal impairment with patients with active and inactive LN according to the SLEDAI-R. A comparison of the active LN group with participants without kidney involvement revealed increased expression of TLR2 (0.078 vs. 0.03, p = 0.009), and TLR5 (0.035 vs. 0.017, p = 0.03). Moreover, a significant difference was observed in TLR9 expression between inactive LN and the control group (0.014 vs. 0.009, p = 0.01), together with borderline correlation in TLR2 expression (0.04 vs. 0.03, p = 0.06; Table 2B, Fig. 2).

Prediction of active LN based on TLR mRNA expression

ROC analysis identified TLR1 and TLR2 as potential predictive markers for active LN among all of the studied TLRs in the entire cohort of 74 patients with active LN versus those without. TLR1 had an AUC of 0.707 (95% CI 0.556–0.857, p = 0.020), and TLR2 had an AUC of 0.715 (95% CI 0.564–0.866, p = 0.015) (Fig. 3, Additional file Fig. 5). The optimal cut-off value for TLR1 mRNA expression in active LN was 0.0636, with a specificity of 0.787 and sensitivity of 0.615. The optimal cut-off value for TLR2 mRNA expression was 0.0282, with a specificity of 0.607 and sensitivity of 0.769. In the case of a summation of TLR1 and TLR2 values for the purpose of improved LN prediction, the sum is 0.095, exhibiting high sensitivity (0.846) but less specificity (0.689).

In conclusion, if a patient tests positive for TLR1 and/or TLR2 (at least one test), there is a high probability of active LN presence (0.923; 95% CI (0.640–0.998)). If a patient tests negative for both TLR1 and TLR2, there is a high probability that kidney damage is not caused by active LN (0.885; 95% CI (0.778–0.953)).

Correlation between TLR mRNA expression and disease activity score

Multivariate analysis was used to assess the associations between TLR1-10 mRNA expression and the SLEDAI-2 K score, SLICC/ACR DI score, complement C3 and C4 levels and anti-ds DNA antibody level in all patients with SLE. This analysis confirmed a statistically significant association between TLR9 levels and C4. Higher TLR9 mRNA levels were associated with active LN (0.0103, p = 0.035) and low complement C4 levels (-0.0759, p = 0.039) (Additional file Fig. 6). A correlation between TLR4 mRNA expression and the SLEDAI-2 K score was also demonstrated (Additional file Fig. 7). However, when the Spearman correlation was applied, no correlation was found. No other significant correlations were detected between TLR1-10 mRNA expression and the activity score or complement or antibody levels.

Correlation between TLR mRNA expression and treatment

A significant difference was observed between patients who were off and on corticosteroid therapy in TLR7 and TLR8 expressions. Patients who received corticosteroid treatment exhibited statistically significantly higher levels of TLR7 (0.011 vs. 0.023, p = 0.039) and TLR8 (0.031 vs. 0.05, p = 0.018) (Additional file Table 7). Upon dividing the patients into subgroups according to renal impairment, no statistically significant difference in TLR expression was observed. It is notable that patients undergoing corticosteroid therapy exhibited heightened SLE activity, as evidenced by elevated SLEDAI scores (3.95 vs. 7.62, p = 0.022). Mann-Whitney U-test showed no significant difference between patients on and off antimalarial therapy in TLR values.

TLR family expression in CNS-SLE versus no CNS-SLE patients

The increased relative mRNA expression of TLR1 (0.088 vs. 0.048, p = 0.01), TLR4 (0.173 vs. 0.066, p = 0.0003) and TLR6 (0.087 vs. 0.036, p = 0.007) was observed in comparison to No CNS-SLE patients. (Table 2C, Additional file Fig. 4).

Correlations between individual TLR mRNA expression

Multiple significant correlations were identified between the expression of individual TLRs in the entire lupus cohort. Specifically, TLR1 correlated with TLR2, TLR4, TLR5, TLR7, TLR8, and TLR10. TLR2 also correlated with TLR5 and TLR8, TLR3 with TLR9, TLR4 with TLR6-9, and TLR5 with TLR8 and TLR10 (Table 3, Additional file Fig. 8). The multiple mutual relationships of TLRs and their correlation with disease activity demonstrate the activation of innate immunity in SLE.

Discussion

This study further supports the important role of TLR in SLE and its severe organ manifestations and offers insight into TLR activation, which is characteristic of kidney and central nervous system involvement. The specific blood profile of TLRs in the LN makes it suitable for minimally invasive diagnostics of active LN but also highlights a novel attractive treatment strategy for SLE by targeting TLR signalling cascades using antagonists or monoclonal antibodies that selectively prevent extracellular or endosomal TLR ligation [7].

The mosaic of TLR expression in SLE patients published thus far in the literature includes a number of interesting aspects. SLE is commonly associated with the overexpression of TLR7, which leads to IFNα production [4, 5]. Genetic or environmental signals, including female sex due to the presence of the TLR7 gene on the X chromosome [16], a high-fat diet and a cytokine environment, including the level of type I interferons, lead to the overexpression of TLR7 and increase the susceptibility of individuals to SLE [4]. Two pathways of B-cell activation lead to the development of autoreactive antibody-secreting cells in SLE: the extrafollicular response and the germinal centre reaction. These pathways can enhance the production of autoantibodies [17]. TLR7 plays a key role in both responses [16]. Targeted inhibition of TLR7-stimulated inflammatory cytokine production has already been tested in lupus mouse and monkey models and is a promising therapeutic option for the treatment of SLE [7].

TLR2 and TLR4 are membrane-bound receptors that recognise bacterial wall components and endogenous components of dying cells [6, 18]. Upregulated TLR2 mRNA expression has been found in SLE patients and is associated with increased production of IL-17 [19]. Overexpression of TLR4 can also lead to lupus progression and the presence of lupus nephritis [20]. TLR5 recognises bacterial flagellin, which activates both innate inflammatory responses and the development of adaptive immunity [21]. The role of TLR5 in SLE remains questionable. There are several single nucleotide polymorphisms (SNPs), such as the TLR5-rs5744168 SNP, that are significantly more frequent in SLE patients with LN [22], but another study revealed no association between this polymorphism and SLE [23]. Increased activation of TLR3 in peripheral mononuclear cells (PBMCs), B cells and T cells has also been observed in SLE patients [24] and may be associated with various SLE phenotypes [25]. TLR8 expression is upregulated in SLE patients, especially in females, due to the presence of 17β-estradiol, which is a potential factor in SLE development [26, 27]. However, some studies have indicated that TLR8 may also inhibit TLR7-mediated autoimmunity and renal inflammation [28]. In contrast, activation of TLR9 has an antagonistic effect and may be protective against lupus-like diseases due to its ability to limit the stimulatory activity of TLR7 [29].

Our study revealed a decrease in TLR10 mRNA expression in patients with SLE in comparison to healthy controls. TLR10 is a plasma membrane receptor that is highly expressed in monocytes, dendritic cells and B cells. Several reports have shown that this receptor has an anti-inflammatory function by inhibiting the production of inflammatory cytokines such as interleukin (IL)-1β, IL-6 and tumor necrosis factor α (TNF-α) [30,31,32]. All of these cytokines are elevated in patients with SLE and are associated with increased disease activity [33]. Our results confirmed a significant downregulation of TLR10 in SLE patients, as previously described [34]. This finding suggests a possible link between SLE development and its activity during TLR10 downregulation in our study group. TLR10 was also downregulated in active rheumatoid arthritis (RA) patients [13] and in patients with microscopic polyangiitis [35]. This finding supports its possible role in autoimmune disease development. Despite these findings, SNPs in the TLR10 gene have been reported to be associated with increased susceptibility to RA, hip osteoarthritis, and other diseases [36, 37]. In contrast, TLR7 and TLR9 mRNA expression, which are most often associated with the development of SLE, was only slightly upregulated in our SLE cohort, and the difference from healthy subjects did not reach statistical significance. Nevertheless, the elevation of TLR9 has been demonstrated in patients with Inactive LN.

Up to 40% of SLE patients develop severe organ manifestations, such as lupus nephritis. This is a risk factor for increased morbidity and all-cause mortality and can lead to progression to end-stage renal disease [38]. Numerous studies have investigated the involvement of TLRs in the pathogenesis of LN. The specific role of TLRs in kidney injury is complex. The available experimental and clinical data indicates that TLRs are involved in the pathogenesis of immune complex glomerulonephritis and intrarenal inflammation [6]. Inflammation of kidney tissue in LN can result from several pathways. Mesangial cells and other renal parenchymal cells express TLR1–4 and TLR6. Renal cells respond to TLR activation by secreting interleukins and chemokines [39, 40]. Another mechanism is that immune complex deposits containing TLR agonists activate glomerular mesangial cells and contribute to immune complex glomerulonephritis [6]. TLR2-mediated activation of podocytes and endothelial cells at the glomerular filtration barrier results in the aggravation of albuminuria [6]. Higher expression of TLR9 was detectable in the tubulointerstitial compartment and glomeruli of LN patients. A significant increase in TLR5 gene expression in LN biopsies has been also observed, confirmed by intense and diffuse staining for TLR5 only in LN tubules [41]. Overexpression of various polymorphisms of TLRs (TLR2, TLR3, TLR4, TLR7, TLR8, and TLR9) has been found in LN patients [6]. In our study, we confirmed the significant overexpression of two members of the TLR family (TLR2 and TLR5) in active LN patients and TLR2 and TLR9 when compared Inactive LN with patients without kidney involvement.

In our study, TLR1 and TLR2 were suggested as possible diagnostic markers for active LN. We proposed cut-off values for the mRNA expression of these receptors. The positivity of TLR1 or TLR2 predicts lupus nephritis with high sensitivity and high predictive value of a negative test. Furthermore, if both are positive, there is a high specificity for active LN and also a high predictive value for a negative test. These findings highlight the potential role of TLR2 as a biomarker of active LN in our study. However, further studies are needed to confirm the high expression of TLR1 and TLR2 as possible diagnostic or therapeutic targets in LN. The overexpression of TLR1 may be related to proinflammatory cytokine release, activation of type I IFN signalling leading to podocyte injury and proteinuria in LN patients [42]. TLR1 is a cell surface receptor that is localised on monocytes, macrophages, dendritic cells and B cells. TLR1 recognises bacterial lipoproteins and glycolipids in complex with TLR2, and its overexpression has been also detected in Lyme disease, sepsis, candidaemia and leprosy [43].

Little is known about the role of TLR6 in the literature. TLR6 is a cell surface TLR that mainly recognises microbial membrane components to induce an inflammatory response [44]. TLR6 forms a heterodimer with TLR2 to mediate the cellular response to gram-positive bacteria, mycoplasma, fungi, some viruses and even protozoa [45]. According to one study, two specific coding and regulatory regions of TLR6 were enriched in LN, although this did not reach statistical significance [41]. Apart from this difference, our study confirmed the results of previous research.

It is already established that antimalarial treatments have an immunomodulatory effect. However, our study did not identify any significant difference in the expression of specific TLRs in patients undergoing antimalarial therapy and those who were not. The impact of corticosteroid treatment on TLRs expression in patients with autoimmune diseases remains poorly understood. The results of our study indicated an elevated expression of TLR7 and TLR8 in the subgroup that had been treated with corticosteroids. One potential explanation for this observation is that the patients in this subgroup exhibited a more active disease state, as indicated by their SLEDAI scores. This is consistent with the literature, which predominantly associates TLR7 and TLR8 with SLE [4, 5, 16, 26, 27]. However, our data on the impact of therapy on TLR expression are limited by the fact that the blood samples were not collected before and during treatment.

Neuropsychiatric systemic lupus erythematosus (NPSLE) is one of the least understood and most complex subtypes of SLE. NPSLE can affect both the central and peripheral nervous systems with a wide range of symptoms. The most common symptoms observed in our patients were organic brain syndrome, seizures, and lupus headache, with a predominance of diffuse CNS manifestations. According to previous studies, two pathological mechanisms contribute to NPSLE development. The ischemic or thrombotic pathway is associated with immune complexes, the complement system and antiphospholipid antibody positivity, while the autoimmune or inflammatory pathway is characterised by autoantibody, proinflammatory cytokine, chemokine, microglial and C1q production [46, 47]. The specific role of TLRs in the development of NPSLE is poorly understood due to the complexity of the problem. The concentration of high-mobility group box protein 1 (HMGB1) in serum and cerebrospinal fluid was elevated in patients with seizure disorders, accompanied by the overexpression of TLR4 mRNA in PBMCs from these patients [48]. Our study confirmed that the overexpression of TLR4 in these patients was accompanied by increased expression of TLR1 and TLR6.

The aim of our study was to investigate the relationship between disease activity and TLRs. We did not find any significant correlation between TLRs and scoring systems such as the SLEDAI-2 K, SLICC/ACR DI, or major antibody and complement levels. The only significant association observed was between TLR9 mRNA relative expression and active lupus nephritis. Additionally, an inverse correlation was found between TLR9 mRNA expression and C4 complement levels. Low C4 complement levels have been shown to be a useful marker for diagnosing and measuring activity in patients with SLE and may also predict renal flares in some patients [49].

Several significant correlations were found between the expression of individual TLRs across the lupus cohort, demonstrating the interconnectedness and complexity of innate immune activation leading to the development of SLE.

The strengths of this study are the detailed analysis of the expression of TLRs in a real-world clinical setting in a group of lupus patients and the statistical evaluation of the mutual relationships of TLRs. The authors of the study are aware of some constraints of their results. First, the study was performed as a cross-sectional analysis in a real-world setting of patients at different stages of the disease. Secondly, the study was performed in an ethnically homogeneous population of the Central Moravia region. It would be of interest to compare the results with those of an analysis of patients of different ethnic backgrounds or of multicentric cohorts. Finally, our study group focused mainly on patients with active LN and NPSLE. Patients with other SLE phenotypes and other different subsets of disease were not included. The authors of this study are aware that the utilisation of transcriptomics in smaller clinical settings may be constrained by elevated costs and the necessity for specific equipment.

Conclusion

The TLR family plays a key role in the innate immune response as well as in the pathogenesis of autoimmune rheumatic diseases. We found heterogeneity in TLR expression according to the main clinical manifestations of SLE. In our study, overexpression of TLR2 was found in both the active and inactive LN groups. ROC curves indicate that the expression of TLR1 and TLR2 is indicative of the presence of active LN and the discriminative values between active LN and other patients with inactive or no LN were established. These findings highlight the potential role of TLR2 as a biomarker of active LN in our study. In addition, the expression and role of TLRs in NPSLE are not well described in the literature. Although our study cannot provide a detailed analysis of the possible pathogenetic link between the expression of TLRs and the development of CNS involvement, it does provide some insight into this issue. Given the small number of patients with this severe manifestation in SLE, further investigation is needed.

In this context, an analysis of innate gene expression patterns can provide useful information for understanding the molecular mechanisms of SLE, which may impact patient prognosis and treatment outcomes. However, the practical use of determining TLR expression requires additional data.

Table 1 Demographic and clinical characteristics of the enrolled patients
Table 2 Relative TLR mRNA expression according to SLE subgroup and healthy control
Table 3 Correlations of Toll-like receptor (TLR) mRNA expression levels
Fig. 1
figure 1

Toll-like receptor mRNA relative expression– SLE patients vs. healthy controls. The relative mRNA expression of the TLR family members was compared between patients with SLE and healthy controls. Compared to that in healthy controls, the relative mRNA expression of TLR4 was upregulated in the SLE group (Table 2A). Conversely, TLR10 expression was significantly lower in the SLE group compared to healthy control group (Table 2A). Furthermore, increased expression of TLR6 and TLR7 and decreased expression of TLR3 were observed, although the differences did not reach statistical significance

Fig. 2
figure 2

TLR family expression in Active LN-SLE, Inactive LN-SLE and No LN-SLE patients. The objective of this study was to investigate the relationship between innate immune gene expression and clinical manifestations in patients with SLE without renal impairment and in patients with active and inactive LN according to the renal SLEDAI. Our findings revealed a significant upregulation of the TLR2 and TLR5 (Table 2B). Furthermore, TLR2 and TLR9 expression was increased when the No LN-SLE and Inactive LN-SLE groups were compared. Additionally, a significant distinction was observed in TLR4 expression when the Active and Inactive LN-SLE group were compared

Fig. 3
figure 3

ROC analysis of TLR1 and TLR2 mRNA expression. ROC analysis indicated that TLR1 and TLR2 were potential predictive markers for active LN among the studied TLRs. The area under the curve (AUC) of TLR1 was 0.707, with a 95% confidence interval (CI) of 0.556–0.857 and a p value of 0.020. The AUC of TLR2 was 0.715, with a 95% CI of 0.564–0.866 and a p value of 0.015. AUC values above 0.7 were regarded as a reasonable discriminating value for predicting diagnostic tests

Data availability

The datasets used and/or analysed during the current study are available from the corresponding author upon reasonable request.

Abbreviations

Anti-dsDNA:

Anti-double-stranded DNA antibodies

ACR:

American College of Rheumatology

CNS:

Central nervous system

C3:

Complement component 3

C4:

Complement component 4

DAMPs:

Damage-associated molecular patterns

HMGB1:

High-mobility group box protein 1

IFNα:

Interferonα

IL:

Interleukin

HPF:

High-power field

LN:

Lupus nephritis

mRNA:

Messenger ribonucleic acid

NPSLE:

Neuropsychiatric systemic lupus erythematosus

PAMPs:

Pathogen-associated molecular patterns

PBMCs:

Peripheral mononuclear cells

SLE:

systemic lupus erythematosus

SLEDAI-2 K:

Systemic Lupus Eryhematosus Disease Activity Index 2000

SLICC/ACR DI:

The Systemic Lupus International Collaborating Clinics American College of Rheumatology Damage index

SNP:

Single nucleotide polymorphism

TLR:

Toll-like receptor

TNF-α:

Tumor necrosis factor α

References

  1. Akhil A, Bansal R, Anupam K, et al. Systemic lupus erythematosus: latest insight into etiopathogenesis. Rheumatol Int. 2023;43(8):1381–93. https://doiorg.publicaciones.saludcastillayleon.es/10.1007/s00296-023-05346-x. Epub 2023 May 24. PMID: 37226016; PMCID: PMC10208557.

    Article  PubMed  PubMed Central  Google Scholar 

  2. Pan L, Lu MP, Wang JH, et al. Immunological pathogenesis and treatment of systemic lupus erythematosus. World J Pediatr. 2020;16(1):19–30. https://doiorg.publicaciones.saludcastillayleon.es/10.1007/s12519-019-00229-3. Epub 2019 Feb 22. PMID: 30796732; PMCID: PMC7040062.

    Article  PubMed  Google Scholar 

  3. Takeuchi O, Akira S. Toll-like receptors; their physiological role and signal transduction system. Int Immunopharmacol. 2001;1(4):625– 35. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/s1567-5769(01)00010-8. PMID: 11357875.

  4. Fillatreau S, Manfroi B, Dörner T. Toll-like receptor signalling in B cells during systemic lupus erythematosus. Nat Rev Rheumatol. 2021;17:98–108. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/s41584-020-00544-4.

    Article  CAS  PubMed  Google Scholar 

  5. García-Ortiz H, Velázquez-Cruz R, Espinosa-Rosales F, et al. Association of TLR7 copy number variation with susceptibility to childhood-onset systemic lupus erythematosus in Mexican population. Ann Rheum Dis. 2010;69:1861–5.

    Article  PubMed  Google Scholar 

  6. Devarapu SK, Anders HJ. Toll-like receptors in lupus nephritis. J Biomed Sci. 2018;25(1):35. https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s12929-018-0436-2. PMID: 29650017; PMCID: PMC5898010.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Manno A, Honda T, Kuwata C et al. Targeting Toll-Like receptor 7 with DS-7011a, a promising novel antagonistic antibody for the treatment of systemic lupus erythematosus [abstract]. Arthritis Rheumatol 2023; 75 (suppl 9).

  8. Aringer M. EULAR/ACR classification criteria for SLE. Semin Arthritis Rheum. 2019;49(3S):S14-S17. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.semarthrit.2019.09.009. PMID: 31779843.

  9. Gladman DD, Ibañez D, Urowitz MB. Systemic lupus erythematosus disease activity index 2000. J Rheumatol. 2002;29(2):288–91. PMID: 11838846.

    PubMed  Google Scholar 

  10. Gladman DD, Goldsmith CH, Urowitz MB, et al. The systemic lupus international collaborating clinics/american college of rheumatology (SLICC/ACR) damage index for systemic lupus erythematosus international comparison. J Rheumatol. 2000;27(2):373–6. PMID: 10685799.

    CAS  PubMed  Google Scholar 

  11. Mina R, Abulaban K, Klein-Gitelman MS, et al. Validation of the lupus nephritis clinical indices in Childhood-Onset systemic lupus erythematosus. Arthritis Care Res (Hoboken). 2016;68(2):195–202. https://doiorg.publicaciones.saludcastillayleon.es/10.1002/acr.22651. PMID: 26213987; PMCID: PMC4720587.

    Article  PubMed  Google Scholar 

  12. The American College of Rheumatology nomenclature and case definitions for neuropsychiatric lupus syndromes. Arthritis Rheum. 1999;42(4):599–608. https://doiorg.publicaciones.saludcastillayleon.es/10.1002/1529-0131(199904)42:4<599::AID-ANR2>3.0.CO;2-F. PMID: 10211873.

  13. Petrackova A, Horak P, Radvansky M et al. Revealed heterogeneity in rheumatoid arthritis based on multivariate innate signature analysis. Clin Exp Rheumatol 2020 Mar-Apr;38(2):289–98. https://doiorg.publicaciones.saludcastillayleon.es/10.55563/clinexprheumatol/qb2ha3. Epub 2019 Aug 3. PMID: 31376255.

  14. Youden WJ. Index for rating diagnostic tests. Cancer. 1950;3(1):32– 5. doi: 10.1002/1097-0142(1950)3:1<32::aid-cncr2820030106>3.0.co;2-3. PMID: 15405679.

  15. Zhou XH, Obuchowski NA. McClish DK Statistical methods in diagnostic medicine. Wiley-Interscience. 2002;(1). Print ISBN:9780471347729|Online ISBN:9780470317082|DOI:10.1002/9780470317082.

  16. Souyris M, Cenac C, Azar P et al. JC. TLR7 escapes X chromosome inactivation in immune cells. Sci Immunol. 2018;3(19):eaap8855. https://doiorg.publicaciones.saludcastillayleon.es/10.1126/sciimmunol.aap8855. PMID: 29374079.

  17. Jenks SA, Cashman KS, Zumaquero E et al. Distinct Effector B Cells Induced by Unregulated Toll-like Receptor 7 Contribute to Pathogenic Responses in Systemic Lupus Erythematosus. Immunity. 2018;49(4):725–739.e6. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.immuni.2018.08.015. Epub 2018 Oct 9. Erratum in: Immunity. 2020;52(1):203. PMID: 30314758; PMCID: PMC6217820.

  18. Anders HJ. Toll-like receptors and danger signaling in kidney injury. J Am Soc Nephrol. 2010;21(8):1270-4. https://doiorg.publicaciones.saludcastillayleon.es/10.1681/ASN.2010030233. Epub 2010 Jul 22. PMID: 20651159.

  19. Liu Y, Liao J, Zhao M, et al. Increased expression of TLR2 in CD4(+) T cells from SLE patients enhances immune reactivity and promotes IL-17 expression through histone modifications. Eur J Immunol. 2015;45(9):2683–93. https://doiorg.publicaciones.saludcastillayleon.es/10.1002/eji.201445219. Epub 2015 Jun 26. PMID: 26079624.

    Article  CAS  PubMed  Google Scholar 

  20. Lee TP, Huang JC, Liu CJ et al. Interactions of surface-expressed TLR-4 and endosomal TLR-9 accelerate lupus progression in anti-dsDNA antibody transgenic mice. Exp Biol Med (Maywood). 2014;239(6):715– 23. https://doiorg.publicaciones.saludcastillayleon.es/10.1177/1535370214525299. PMID: 24719374.

  21. Hayashi F, Smith KD, Ozinsky A et al. The innate immune response to bacterial flagellin is mediated by Toll-like receptor 5. Nature. 2001;410(6832):1099– 103. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/35074106. PMID: 11323673.

  22. Gharbia OM, Bassiouni SAR, Zaki MES, et al. Toll-like receptor 5 and Toll-like receptor 9 single nucleotide polymorphisms and risk of systemic lupus erythematosus and nephritis in Egyptian patients. Egypt Rheumatol Rehabil. 2021;48:49. https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s43166-021-00093-y.

    Article  Google Scholar 

  23. Lee YH, Choi SJ, Ji JD, et al. Association between toll-like receptor polymorphisms and systemic lupus erythematosus: a meta-analysis update. Lupus. 2016;25(6):593–601. https://doiorg.publicaciones.saludcastillayleon.es/10.1177/0961203315622823.

    Article  CAS  PubMed  Google Scholar 

  24. Klonowska-Szymczyk A, Wolska A, Robak T, et al. Expression of toll-like receptors 3, 7, and 9 in peripheral blood mononuclear cells from patients with systemic lupus erythematosus. Mediators Inflamm. 2014;2014:381418. https://doiorg.publicaciones.saludcastillayleon.es/10.1155/2014/381418. Epub 2014 Feb 25. PMID: 24692849; PMCID: PMC3955595.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Wang CM, Chang SW, Wu YJ, et al. Genetic variations in Toll-like receptors (TLRs 3/7/8) are associated with systemic lupus erythematosus in a Taiwanese population. Sci Rep. 2014;4:3792. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/srep03792. PMID: 24445780; PMCID: PMC3896912.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Young NA, Wu LC, Burd CJ, et al. Estrogen modulation of endosome-associated toll-like receptor 8: an IFNα-independent mechanism of sex-bias in systemic lupus erythematosus. Clin Immunol. 2014;151(1):66–77. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.clim.2014.01.006. Epub 2014 Jan 24. PMID: 24525049; PMCID: PMC4066385.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Young NA, Valiente GR, Hampton JM, et al. Estrogen-regulated STAT1 activation promotes TLR8 expression to facilitate signaling via microRNA-21 in systemic lupus erythematosus. Clin Immunol. 2017;176:12–22. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.clim.2016.12.005. Epub 2016 Dec 27. PMID: 28039018; PMCID: PMC5815376.

    Article  CAS  PubMed  Google Scholar 

  28. Desnues B, Macedo AB, Roussel-Queval A, et al. TLR8 on dendritic cells and TLR9 on B cells restrain TLR7-mediated spontaneous autoimmunity in C57BL/6 mice. Proc Natl Acad Sci U S A. 2014;111(4):1497–502. https://doiorg.publicaciones.saludcastillayleon.es/10.1073/pnas.1314121111. Epub 2014 Jan 13. PMID: 24474776; PMCID: PMC3910605.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Nündel K, Green NM, Shaffer AL, et al. Cell-intrinsic expression of TLR9 in autoreactive B cells constrains BCR/TLR7-dependent responses. J Immunol. 2015;194(6):2504–12. https://doiorg.publicaciones.saludcastillayleon.es/10.4049/jimmunol.1402425. Epub 2015 Feb 13. PMID: 25681333; PMCID: PMC4382804.

    Article  CAS  PubMed  Google Scholar 

  30. Hess NJ, Felicelli C, Grage J, et al. TLR10 suppresses the activation and differentiation of monocytes with effects on DC-mediated adaptive immune responses. J Leukoc Biol. 2017;101(5):1245–52. https://doiorg.publicaciones.saludcastillayleon.es/10.1189/jlb.3A1116-492R. Epub 2017 Feb 24. PMID: 28235773; PMCID: PMC5380377.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Oosting M, Cheng SC, Bolscher JM, et., Al. Human TLR10 is an anti-inflammatory pattern-recognition receptor. Proc Natl Acad Sci U S A. 2014;111(42):E4478–84. https://doiorg.publicaciones.saludcastillayleon.es/10.1073/pnas.1410293111. Epub 2014 Oct 6. PMID: 25288745; PMCID: PMC4210319.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Fore F, Indriputri C, Mamutse J, et al. TLR10 and its unique Anti-Inflammatory properties and potential use as a target in therapeutics. Immune Netw. 2020;20(3):e21. https://doiorg.publicaciones.saludcastillayleon.es/10.4110/in.2020.20.e21. PMID: 32655969; PMCID: PMC7327153.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Umare V, Pradhan V, Nadkar M, et al. Effect of Proinflammatory cytokines (IL-6, TNF-α, and IL-1β) on clinical manifestations in Indian SLE patients. Mediators Inflamm. 2014;2014:385297. https://doiorg.publicaciones.saludcastillayleon.es/10.1155/2014/385297. Epub 2014 Dec 7. PMID: 25548434; PMCID: PMC4273527.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Petrackova A, Horak P, Radvansky M, et al. Cross-Disease innate gene signature: emerging diversity and abundance in RA comparing to SLE and SSc. J Immunol Res. 2019;2019:3575803. https://doiorg.publicaciones.saludcastillayleon.es/10.1155/2019/3575803. PMID: 31396542; PMCID: PMC6664489.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Lai Y, Xue C, Liao Y et al. Differential Expression of Toll-Like Receptor Signaling Pathway Is Associated with Microscopic Polyangiitis in Peripheral Blood Neutrophils. Immunol Invest. 2017;46(4):375–384. https://doiorg.publicaciones.saludcastillayleon.es/10.1080/08820139.2017.1288236. Epub 2017 Apr 4. PMID: 28375762.

  36. Zhang Y, Cao R, Ying H et al. Increased Expression of TLR10 in B Cell Subsets Correlates with Disease Activity in Rheumatoid Arthritis. Mediators Inflamm. 2018;2018:9372436. https://doiorg.publicaciones.saludcastillayleon.es/10.1155/2018/9372436. Erratum in: Mediators Inflamm. 2019;2019:8419439. PMID: 30686934; PMCID: PMC6327257.

  37. Tang H, Cheng Z, Ma W, et al. TLR10 and NFKBIA contributed to the risk of hip osteoarthritis: systematic evaluation based on Han Chinese population. Sci Rep. 2018;8:10243. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/s41598-018-28597-2.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Pryor KP, Barbhaiya M, Costenbader KH, et al. Disparities in lupus and lupus nephritis care and outcomes among US medicaid beneficiaries. Rheum Dis Clin North Am. 2021;47(1):41–53. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.rdc.2020.09.004. Epub 2020 Oct 29. PMID: 34042053; PMCID: PMC8171807.

    Article  PubMed  Google Scholar 

  39. Patole PS, Pawar RD, Lech M et al. Expression and regulation of Toll-like receptors in lupus-like immune complex glomerulonephritis of MRL-Fas(lpr) mice. Nephrol Dial Transplant. 2006;21(11):3062-73. https://doiorg.publicaciones.saludcastillayleon.es/10.1093/ndt/gfl336. Epub 2006 Sep 5. PMID: 16954173.

  40. Imaizumi T, Aizawa T, Segawa C, et al. Toll-like receptor 3 signaling contributes to the expression of a neutrophil chemoattractant, CXCL1 in human mesangial cells. Clin Exp Nephrol. 2015;19(5):761–70. https://doiorg.publicaciones.saludcastillayleon.es/10.1007/s10157-014-1060-4. Epub 2014 Dec 4. PMID: 25471749.

    Article  CAS  PubMed  Google Scholar 

  41. Elloumi N, Fakhfakh R, Abida O, Ayadi L, Marzouk S, Hachicha H, Fourati M, Bahloul Z, Mhiri MN, Kammoun K, Masmoudi H. Relevant genetic polymorphisms and kidney expression of Toll-like receptor (TLR)-5 and TLR-9 in lupus nephritis. Clin Exp Immunol. 2017;190(3):328–39. https://doiorg.publicaciones.saludcastillayleon.es/10.1111/cei.13022. Epub 2017 Aug 30. PMID: 28763101; PMCID: PMC5680057.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Yavuz S, Bianchi M, Kozyrev S, et al. Toll-like receptors revisited; a possible role for TLR1 in lupus nephritis. Ann Rheum Dis. 2021;80(3):404–6. https://doiorg.publicaciones.saludcastillayleon.es/10.1136/annrheumdis-2020-218373. Epub 2020 Sep 29. PMID: 32994161; PMCID: PMC7892377.

    Article  PubMed  Google Scholar 

  43. Rahman S, Shering M, Ogden NH, et al. Toll-like receptor cascade and gene polymorphism in host-pathogen interaction in Lyme disease. J Inflamm Res. 2016;9:91–102. https://doiorg.publicaciones.saludcastillayleon.es/10.2147/JIR.S104790. PMID: 27330321; PMCID: PMC4898433.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Kawai T, Akira S. The role of pattern-recognition receptors in innate immunity: update on Toll-like receptors. Nat Immunol. 2010;11(5):373–84. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/ni.1863. Epub 2010 Apr 20. PMID: 20404851.

    Article  CAS  PubMed  Google Scholar 

  45. Oliveira-Nascimento L, Massari P, Wetzler LM. The role of TLR2 in infection and immunity. Front Immunol. 2012;3:79. https://doiorg.publicaciones.saludcastillayleon.es/10.3389/fimmu.2012.00079. PMID: 22566960; PMCID: PMC3342043.

    Article  PubMed  PubMed Central  Google Scholar 

  46. Sarwar S, Mohamed AS, Rogers S, et al. Neuropsychiatric systemic lupus erythematosus: A 2021 update on diagnosis, management, and current challenges. Cureus. 2021;13(9):e17969. https://doiorg.publicaciones.saludcastillayleon.es/10.7759/cureus.17969. PMID: 34667659; PMCID: PMC8516357.

    Article  PubMed  PubMed Central  Google Scholar 

  47. Lindblom J, Mohan C, Parodis I. Biomarkers in neuropsychiatric systemic lupus erythematosus: A systematic literature review of the last decade. Brain Sci. 2022;12:192. https://doiorg.publicaciones.saludcastillayleon.es/10.3390/brainsci12020192.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Huang Q, Shen S, Qu H, et al. Expression of HMGB1 and TLR4 in neuropsychiatric systemic lupus erythematosus patients with seizure disorders. Ann Transl Med. 2020;8(1):9. https://doiorg.publicaciones.saludcastillayleon.es/10.21037/atm.2019.12.44. PMID: 32055600; PMCID: PMC6995744.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Birmingham DJ, Irshaid F, Zou X, et al. The complex nature of serum C3 and C4 as biomarkers of renal flare. Lupus. 2010;19:1272–80. https://doiorg.publicaciones.saludcastillayleon.es/10.1177/0961203310371154.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

Supported by MH CZ– DRO (FNOl, 00098892), IGA_LF_2025_003, Faculty of Medicine and Dentistry, Palacky University Olomouc and University Hospital, Olomouc, Czech Republic.

Author information

Authors and Affiliations

Authors

Contributions

D. M.: Patient selection, obtaining informed consent, assessment of disease activity, major contributor in writing the manuscript, S. A., V. J., S. M.: Patient selection, obtaining informed consent, assessment of disease activity. P. A., M. J., S. K. V., M.F.: Immunological analysis. P.A., R. M., K. M., L. K.: Statistical analysis, images. K. E., H. P.: Study design, secondary contribution to the writing of the manuscript. H.P.: Corresponding author. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Pavel Horak.

Ethics declarations

Ethics approval and consent to participate

The depersonalised and anonymous data and samples were collected and stored after signed written informed consent for data collection was obtained from all participants and healthy controls. The study was approved by the local ethical committee (Multicenter Research Ethics Committee of the University Hospital Olomouc and Faculty of Medicine and Dentistry, Palacky University Olomouc, Czech Republic). The study was performed in accordance with the Helsinki Declaration of 1964 and its later amendments. No animal data were used in this study.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Footnote

The findings in this manuscript were presented in virtual poster format at the Annual European Congress of Rheumatology, EULAR 2024, Vienna, Austria, is available in the EULAR Abstract Archive.

Additional information

Publisher’s note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary Material 1

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Dudkova, M., Petrackova, A., Radvansky, M. et al. Blood gene expression of Toll-like receptors in SLE patients with lupus nephritis or neuropsychiatric systemic lupus erythematosus. Arthritis Res Ther 27, 41 (2025). https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s13075-025-03512-z

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s13075-025-03512-z

Keywords