Skip to main content

ORIGINAL RESEARCH article

Front. Immunol., 22 February 2023
Sec. Autoimmune and Autoinflammatory Disorders : Autoimmune Disorders
This article is part of the Research Topic Immunoglobulin A in Infection and Inflammation View all 8 articles

Salivary IgA subtypes as novel disease biomarkers in systemic lupus erythematosus

  • 1Red de Apoyo a la Investigación, Universidad Nacional Autónoma de México e Instituto Nacional de Ciencias Médicas y Nutrición Salvador Zubirán, Mexico City, Mexico
  • 2Facultad de Medicina, Universidad Nacional Autónoma de México, Mexico City, Mexico
  • 3Departamento de Biomedicina Molecular, Centro de Investigación y de Estudios Avanzados del Instituto Politécnico Nacional, Mexico City, Mexico
  • 4Departamento de Medicina Interna, Instituto Nacional de Ciencias Médicas y Nutrición Salvador Zubirán, Mexico City, Mexico
  • 5Departamento de Inmunología y Reumatología, Instituto Nacional de Ciencias Médicas y Nutrición Salvador Zubirán, Mexico City, Mexico

Introduction: Immunoglobulin A (IgA) is the main antibody isotype in body fluids such as tears, intestinal mucous, colostrum, and saliva. There are two subtypes of IgA in humans: IgA1, mainly present in blood and mucosal sites, and IgA2, preferentially expressed in mucosal sites like the colon. In clinical practice, immunoglobulins are typically measured in venous or capillary blood; however, alternative samples, including saliva, are now being considered, given their non-invasive and easy collection nature. Several autoimmune diseases have been related to diverse abnormalities in oral mucosal immunity, such as rheumatoid arthritis, Sjogren’s syndrome, and systemic lupus erythematosus (SLE).

Methods: We decided to evaluate the levels of both IgA subtypes in the saliva of SLE patients. A light chain capture-based ELISA measured specific IgA1 and IgA2 levels in a cohort of SLE patients compared with age and gender-matched healthy volunteers.

Results: Surprisingly, our results indicated that in the saliva of SLE patients, total IgA and IgA1 subtype were significantly elevated; we also found that salivary IgA levels, particularly IgA2, positively correlate with anti-dsDNA IgG antibody titers. Strikingly, we also detected the presence of salivary anti-nucleosome IgA antibodies in SLE patients, a feature not previously reported elsewhere.

Conclusions: According to our results and upon necessary validation, IgA characterization in saliva could represent a potentially helpful tool in the clinical care of SLE patients with the advantage of being a more straightforward, faster, and safer method than manipulating blood samples.

Introduction

Systemic lupus erythematosus (SLE) is an autoimmune disease characterized by a deleterious immune response affecting several organs and tissues (1), challenging diagnostic and treatment approaches. SLE is typically defined by the presence of high titers of circulating autoantibodies (2) with abnormal numbers of T and B lymphocytes (3). Currently, the diagnostic criteria for SLE are predominantly based on the presence of clinical manifestations as well as the results of laboratory tests, such as low levels of C3/C4 complement components and the presence of anti-dsDNA or anti-Sm antibodies (4). However, these currently available markers for SLE diagnosis remain suboptimal regarding either their sensitivity or specificity (5). Consequently, more robust biomarkers for SLE are still needed to accurately diagnose patients, monitor disease progression and treatment effectiveness or predict future flares.

The common consensus in systemic lupus erythematosus (SLE) is that there is a general breakdown in lymphocyte tolerance. Consequently, autoreactive B cells are usually considered one of the central effector cell subsets responsible for maintaining inflammatory status in patients. Thus, B cells are considered among the main therapeutic targets of SLE treatment. Furthermore, B cell-targeted therapies have been applied to treat several autoimmune diseases such as pemphigus, multiple sclerosis, ANCA-associated vasculitis, and rheumatoid arthritis; therefore, hypogammaglobulinemia is recognized as a potential complication of these treatments (6, 7).

The European League Against Rheumatism (EULAR) and The American Academy of Asthma, Allergy, and Immunology (AAAAI) guidelines recommend assessing baseline immune function by testing serum immunoglobulins before or even after rituximab treatment in autoimmune disease (811). Nevertheless, only a small number of reports have emerged regarding aberrant immunoglobulin levels prior to starting B cell-depleting approaches in most autoimmune disorders; for example, polyclonal hypergammaglobulinemia is well documented in SLE. However, hypogammaglobulinemia has also been lupus-associated in patients exhibiting selective-isotype deficiencies (12).

Humoral immunity is essentially assessed by measuring titers of IgG-dominated serum antibodies, thus neglecting the contribution of IgA as the major immunoglobulin isotype in humans. To dimension that, it has been documented that the daily IgA production rates around 70 mg/kg of body weight exceed that of all other antibody isotypes combined (1315). Moreover, IgA is the predominant antibody isotype in external secretions, including tears, intestinal mucous, colostrum, milk, and saliva (1315), thus being well-known as the main mucosal immunoglobulin, playing a fundamental role as an immunological barrier that recognizes and excludes human pathogens (1315). Beyond that, there are two human IgA subclasses: IgA1 and IgA2, both ubiquitously present but displaying a differential distribution in the body: IgA1 is dominant in serum or saliva, while IgA2 is most abundant in some other secretion fluids and the colon (1620). In all body fluids, both IgA1 and IgA2 are mainly present as dimeric secretory IgA (SIgA) (15).

Numerous reports have documented a link between alterations in SIgA levels and inflammatory or autoimmune entities, nearly all using saliva for sampling purposes. In this way, patients with SLE, type I diabetes, oral lichen planus, overweight/obesity, oral submucous fibrosis, ankylosing spondylitis, rheumatoid arthritis, mixed connective tissue disease, and Sjögren’s syndrome have displayed significant higher total or antigen-specific salivary SIgA content than their healthy individuals’ counterparts (2127). Although showing differences in SIgA titers, most of these studies have excluded salivary IgA1 and IgA2 subtypes assessment.

Information regarding the levels of salivary IgA subclasses in SLE is still lacking across existing reports. Most of the available efforts have been limited to the measurement of total IgA either in serum or saliva. Consequently, this study was undertaken to investigate the possible alterations in salivary IgA1 and IgA2 in a small cohort of patients with SLE and to characterize their association with disease features.

Methods

Patients and Healthy individuals

The study population comprised 14 healthy individuals and 38 SLE subjects divided into the following groups: 27 inactive lupus (SLEDAI ≤6), 11 with active disease (SLEDAI >6) were recruited from the department of Immunology and Rheumatology of the Instituto Nacional de Ciencias Médicas y Nutrición Salvador Zubirán. All SLE patients fulfilled ACR/SLICC 2012 classification criteria (28), and disease activity was addressed by the SLE disease activity index (SLEDAI). We excluded subjects with ongoing acute or chronic infections (i.e., HIV or viral hepatitis), pregnancy, and patients with a diagnosis of other concomitant autoimmune diseases except for antiphospholipid (aPL) syndrome. None of the study participants received any B cell-depleting or other biological therapies. This study was approved by the Institutional Ethics and Research Committees of the Instituto Nacional de Ciencias Médicas y Nutrición Salvador Zubirán (Ref. 2306). The patients/participants provided their written informed consent to participate prior to inclusion in the study. Demographic and Clinical characteristics of the study population are depicted in Table 1.

TABLE 1
www.frontiersin.org

Table 1 Demographics, clinical and laboratory features of the cohort.

Saliva sampling

All recruited subjects briefly rinsed their mouth with purified water before were asked to deposit a total volume of around 2 mL of unstimulated whole saliva by passive drooling (letting the saliva drop) into a sterile 15 mL polyethylene centrifuge tube. Immediately after collection, the tubes were transported to the laboratory into an ice bucket and were centrifuged for 15 min at 10,000 g and 4°C to remove debris and cells. The supernatants were then separated and added with 2x protease inhibitors (Pierce, Protease Inhibitor Tablets, Thermo Scientific) and finally stored at -70°C until used. All samples (from patients and healthy donors) were treated with this same procedure, not delaying more than 5 min from collection to centrifugation to minimize the risk of protein degradation.

IgA1 and IgA2 quantification in saliva

Detection of IgA1 and IgA2 in saliva was performed by developing a sandwich ELISA previously reported (29). We used flat bottom microtiter plates (Thermo Scientific) that were coated overnight at 4°C with 1 µg/mL of anti-human Ig light chain antibody in 0.2 M Na2CO3/NaHCO3, pH: 9.4. Blocking was performed by using phosphate-buffered saline 1x containing 0.05% Tween-20 at room temperature for 2 h. Fifty microliters of saliva samples (in triplicate) and standard samples were pipetted into the microtiter wells and incubated for 2 h at 37°C. Among the various incubation steps, the wells were washed five times with phosphate-buffered saline 1x containing 0.05% Tween-20. Then, we added biotinylated anti-human IgA1 antibody (monoclonal mouse antibody, Abcam, Cat. Num. ab99796) diluted 1:2000 or biotinylated anti-human IgA2 antibody (monoclonal mouse antibody, Abcam, Cat. Num. ab128731) diluted 1:1000 in phosphate-buffered saline 1x containing 0.05% Tween-20 and incubated for 1 h at 37°C. After washing, we added Streptavidin-HRP diluted 1:5000 in phosphate-buffered saline 1x and incubated for 1h at 37°C. Then we washed six times with phosphate-buffered saline 1x containing 0.5% Tween-20. After a final wash, 50µL of the substrate solution, tetramethylbenzidine, was added, and the plates were incubated for 5min (to IgA1) and 10 min (to IgA2), and we used to stop enzyme reaction 50µL of 1N HCl. The absorbance was measured at 450nm using a microplate spectrophotometer. When indicated, the total IgA (tIgA) amount was obtained as the sum of measured IgA1 and IgA2 concentrations.

Detection of salivary IgA anti-nucleosome and IgA anti-double stranded DNA antibodies

Salivary IgA anti-dsDNA antibodies and IgA anti-nucleosome were detected by ELISA kits QUANTA Lite ® HA dsDNA and HA Nucleosome, respectively, employing an HRP anti-tIgA antibody for detection and following manufacturer instructions. The plates were then read at 450nm on a Bio-Rad xMark spectrophotometer.

Statistical analysis

As indicated, differences between groups were analyzed using two-way ANOVA followed by Tukey’s post hoc test, Kruskal-Wallis tests followed by Dunn’s multiple comparisons tests, or Mann-Whitney U tests. The correlation of IgA1, IgA2, tIgA, laboratory, and clinical features, were evaluated by calculating the Spearman’s rank correlation coefficient. A p value less than 0.05 was considered statistically significant. ROC curve analysis was performed to distinguish IgA between healthy individuals and SLE patients or between inactive and active patients. All statistical analysis was performed using GraphPad Prism 9 software. Additionally, sensibility and specificity parameters for ROC curves were assessed using MedCalc 20.215 software.

Results

IgA1 is the predominant subtype in SLE saliva

We evaluated the levels of IgA antibodies in saliva, considering that the levels of the immunoglobulin subtypes vary in the different mucosal surfaces (30). A total of 52 individuals: 38 patients with SLE and 14 healthy individuals, were then evaluated to assess their concentration of SIgA salivary subtypes. As shown in Figure 1, the most dominant isotype in saliva was IgA1, but most importantly, healthy individuals showed significantly lower levels of IgA1 and IgA2 than patients with SLE.

FIGURE 1
www.frontiersin.org

Figure 1 Differences in salivary IgA subtype levels in healthy individuals and SLE patients. Data were analyzed by an ordinary two-way ANOVA followed by Tukey’s post hoc test. ***p<0.001, ****p<0.0001.

SLE patients exhibit higher levels of salivary SIgA subtypes than healthy controls

Given the clinical differences between inactive and active patients, we segregated our patients into these two groups as described in Methods. Then we analyzed their salivary IgA1, IgA2, and total IgA (tIgA) levels, as shown in Figure 2. Interestingly, both IgA subtypes and tIgA are significantly elevated in inactive and active patients compared to healthy individuals. Nonetheless, only increased tIgA levels were significantly different according to disease activity.

FIGURE 2
www.frontiersin.org

Figure 2 Increase of salivary IgA subtypes in inactive and active SLE patients. (A) Levels of IgA1 between healthy and inactive and active SLE patients. (B) Levels of IgA2 between healthy and inactive and active SLE patients. (C) Levels of tIgA (add of IgA1+IgA2) between healthy and inactive and active SLE patients. Data were assessed by Kruskall-Wallis tests followed by Dunn’s multiple comparisons tests. ns, not statistically significant, *p<0.05, **p<0.01, ***p<0.001, ****p<0.0001.

Increasing concentrations in salivary subtypes of IgA correlate with titers of circulating anti-dsDNA antibodies in SLE

Trying to assess the importance of salivary IgA1 and IgA2 levels in SLE, we performed correlation analyses between these immunoglobulin isotype concentrations in saliva and different clinical and laboratory features of SLE activity (Figure 3A).

FIGURE 3
www.frontiersin.org

Figure 3 Salivary IgA1, IgA2, and tIgA correlated with different clinical variables in SLE. (A) Correlation matrix showing a graphical representation of calculations between IgA1, IgA2, tIgA, and laboratory and clinical variables of SLE patients. The underlying color scale indicates Spearman’s coefficient values. Correlation analysis between antibodies anti-dsDNA with salivary IgA1 (B), salivary IgA2 (C), and salivary tIgA (D). Red slopes present a positive correlation. Correlations were assessed by calculating Spearman’s rank correlation coefficient (r) and p values (p) depicted. ***p<0.001, ****p<0.0001.

Remarkably, we only found positive and highly significant correlations between circulating anti-dsDNA with the salivary concentration of IgA1 (Figure 3B), IgA2 (Figure 3C), and tIgA (Figure 3D) in patients with SLE.

Anti-nucleosome IgA autoantibodies are present and increased in the saliva of patients with SLE

Since the clinical diagnosis of SLE usually involves the detection of different blood-circulating autoantibodies, including levels of anti-dsDNA (double-stranded DNA) IgG and anti-nucleosome IgG in serum, we become interested in finding analog autoantibodies of IgA isotype in other types of samples such as saliva. Therefore, we evaluate the presence of salivary IgA-ANAs with nuclear or cytoplasmic staining patterns by indirect immunofluorescence. As shown in Supplementary Fig. 1, we detected typical autoimmune recognition patterns over HEp-2 cells when an anti-IgA was used as a detection antibody.

To address if salivary IgA fraction recognizes specific nuclear antigens like nucleosomes or dsDNA in our cohort, we performed commercially available ELISA assays employing an anti-human IgA detection antibody instead of the regular anti-human IgG used. As depicted by Figure 4A, we did not detect differences in salivary anti-dsDNA IgA antibodies between samples of healthy donors and SLE patients. In contrast, when IgA anti-nucleosome were measured, we found significantly higher titers of these autoantibodies in SLE patients compared with healthy individuals (Figure 4B). To verify if IgA anti-nuclear antigens were also present in the circulation of these patients, we measured these same autoantibodies in their serum samples. As expected, we found both anti-dsDNA IgA and anti-nucleosome IgA significantly elevated in SLE patients (Supplementary Fig. 2).

FIGURE 4
www.frontiersin.org

Figure 4 Levels of salivary IgA anti-dsDNA and salivary IgA anti-nucleosome in patients with SLE. (A) Salivary IgA anti-dsDNA in healthy individuals and SLE patients. (B) Salivary IgA anti-nucleosome in healthy individuals and SLE patients. Data were assessed by Mann-Whitney U tests. ns, not statistically significant, ****p<0.0001.

IgA subtypes as potential salivary biomarkers of SLE

Finally, to assess the usefulness of salivary IgA measurement as a potential SLE biomarker, we generated receiver operating characteristic (ROC) curves to determine the discriminative capacity of salivary IgA subtypes (IgA1 or IgA2) concentration in SLE patients vs. healthy donors (Figure 5A). We determined the corresponding areas under the curve (AUC), and we observed that IgA1 displayed an outstanding discriminative value (AUC of 0.855), but even IgA2 displayed a good value (AUC = 0.761).

FIGURE 5
www.frontiersin.org

Figure 5 Receiver operator characteristics (ROC) curves of subtypes of salivary IgA for SLE discrimination and SLE activity prediction. (A) ROC curves of IgA subtypes for SLE discrimination. Data of n= 38 SLE patients and n=14 healthy individuals. (B) ROC curves of IgA subtypes for SLE activity prediction. Data of n= 27 inactive SLE patients and n=11 active SLE patients. Area under the curve (AUC) values and p values (p) are depicted.

As mentioned, SLE diagnosis relies on composite data from clinical manifestations and biomarkers such as serum autoantibodies (anti-dsDNA or anti-nucleosome) or complement (C3/C4) levels. Interestingly, when we compared the diagnostic robustness of salivary SIgA concentrations through their displayed AUC or sensitivity at 95% of specificity values with those from anti-dsDNA, anti-nucleosome, and complement levels in independent SLE cohorts reported previously (31, 32) as seen in Supplementary Table 1, these sensitivity/specificity values from SIgA are very similar to those displayed by traditional biomarkers, being IgA1 concentration the more robust.

Besides that, trying to distinguish between SLE inactive and active disease state contingent salivary IgA subclasses levels, we constructed an additional pair of ROC curves segregating these groups of patients. As shown in (Figure 5B) , the discriminative capacity of both IgA subtypes in saliva was satisfactory but lower than the diagnostic approach, according to their displayed AUC values (above 0.6 in both cases).

Discussion

Humans generate about 1.5 L of saliva, which is secreted by the three major types of salivary glands (parotid, submandibular, and sublingual glands) and the minor salivary glands. Among the most recognized components of this fluid, secretory IgA constitutes the predominant immunoglobulin isotype in this and other body secretions.

Since IgA is also the predominant antibody isotype produced by the human body (33), it possesses an essential role in health and several diseases. Interestingly, increases in serum IgA and IgA autoantibodies have been reported in different autoimmune disorders like Sjogren’s syndrome, rheumatoid arthritis, IgA nephropathy, inflammatory bowel disease, and SLE (3437).

On the other hand, the evidence on salivary IgA in those autoimmune diseases is limited and mostly inconclusive (3437). Most studies are focused on reporting individuals with low IgA levels since IgA deficiency’s mortality and morbidity rates correlate with SLE activity, mainly due to recurrent infections in these patients (38). Conversely, only one previous report informs about significantly higher levels of salivary (total) IgA in patients with SLE compared to healthy individuals (39), a reason why we wanted to go further with this observation evaluating the IgA subclasses.

Beyond total IgA fraction, this antibody isotype consists of two subclasses in humans where their proportions vary depending on the mucosal site: in saliva, we usually detect around 60% of IgA1 and close to 40% of IgA2 in healthy individuals (30); however, there was no previous data on salivary IgA subtypes in autoimmune diseases such as SLE.

Interestingly, we found that both salivary IgA1 and IgA2 are elevated in SLE patients compared to healthy individuals, being IgA1 levels significantly higher than those from IgA2. Besides that, when the SLE group was segregated into inactive and active patients, we observed a clear trend towards increased IgA subtype levels in active SLE individuals. Although these elevations did not exhibit significant differences when each IgA subtype was evaluated, the analysis of the levels of total IgA (tIgA=IgA1+IgA2) displayed a significant increase according to global disease activity, an observation that was not previously reported. Accordingly, increases in the antibody fraction (IgA+IgG) in the saliva of lupus patients were observed in a previous study (39), but it does not refer to any individual subtypes’ characterization. Regarding subtypes, the only prior published evidence comes from Roos Ljungberg et al., which found salivary IgA1 and IgA2 anti-citrullinated protein antibodies in rheumatoid arthritis patients (26). The authors found a strong association between salivary IgA and disease activity, even better than serum IgA, suggesting effector mechanisms in this disease pathogenesis due to oral mucosal immune responses to citrullinated proteins (26). However, the significance of altered salivary IgA subclasses was not yet elucidated in this or other autoimmune diseases.

Our results suggest that increased salivary IgA could be associated with disease activity in SLE. Interestingly, observing the absolute amounts of each antibody subclass, and although IgA2 is clearly elevated in both groups of patients, it becomes evident that the most dramatic increase is in the IgA1 fraction that displays a mean value below 2 μg/mL in healthy controls and close to 10 μg/mL in active SLE patients. This last observation is remarkable since, in the oral cavity, IgA1 would be prone to degradation mediated by bacterial proteases, given its particular structure with a larger hinge region (40, 41). This rise in salivary IgA1 concentrations must result from an overproduction of this antibody that could probably be related to the increments in systemic IL-10 levels as previously reported in lupus patients (42, 43). IL-10 is an anti-inflammatory cytokine that promotes B cell responses and plays a pathogenic role in SLE (44). Beyond that, serum levels of this cytokine have been reported to correlate with lupus disease activity (45). Furthermore, it has been proposed that the class switch to IgA1 is mediated by IL-10 and TGF-β (46, 47). As IL-10 was also previously reported elevated in the saliva of SLE patients (48), the dysregulation exerted by that cytokine in the oral microenvironment could be a key element that supports our data.

Interestingly, both subclasses and the total amount of IgA showed robust and highly significant positive correlations with the circulating anti-dsDNA autoantibodies (defined as IgG isotype in serum). Anti-dsDNA antibodies represent a hallmark of SLE and constitute an effective parameter for diagnosing and classifying patients. Additionally, their fluctuating titers during the progression of the disease reflect its activity in many patients and even may predict disease relapse (49). Still limited by our cohort size, these strong correlation values raise the feasibility of employing salivary IgA measuring, either total or subtypes, as a surrogate marker of SLE activity. This possibility needs further exploration in more significant and prospective cohorts.

So far, we have discussed changes in the amount of IgA regardless of its specificity. Thus, with these salivary antibodies correlating with a well-established biomarker as serum IgG anti-dsDNA, we determined if saliva could contain these and other IgA-isotype anti-nuclear antibodies (ANAs). So, beyond the qualitative findings depicted in our autoreactive IgA detection by immunofluorescence employing HEp-2 cells, we decided to perform a semi-quantification of specific salivary IgA ANAs. Surprisingly, we could not find any difference regarding IgA anti-dsDNA. However, we detected a highly significant increase of IgA anti-nucleosome antibody titers in the saliva of patients with SLE versus healthy individuals. To demonstrate that these IgA autoantibodies were also present in our patient’s blood, we measured their levels with the same approach in serum samples. In this regard, we found that circulating IgA anti-dsDNA or IgA anti-nucleosome levels were highly increased in SLE patients, as expected, since they have been reported as elevated previously (50, 51).

It is important to clarify that the method used to quantify these autoantibodies is limited by the unavailability of IgA standards, thus not making it possible to measure them beyond a semi-quantitative approach. Besides that, saliva constitution (in contrast with the serum used to develop and validate the commercially available tests employed here) could interfere with the adequate determination of these autoantibodies and partially explain why our measurements in this fluid seem to present such a high fluorescence background. Besides this problem, which requires an independent standardization of a new ELISA method for saliva, the significantly increased levels of salivary IgA anti-nucleosome antibodies are still of great interest. Several reports mention those anti-nucleosome antibodies as a disease activity marker in patients with SLE (52) and even correlate better than other conventional biomarkers (C3/C4 or anti-dsDNA) with SLE disease activity over time (53), making of the measurement of salivary IgA anti-nucleosome an attractive possibility for the routine monitoring of SLE patients in clinical practice due to the easy access and availability of the type of sample contingent on the validation of these observations in larger/multicenter cohorts and the development of specific antibody tests for this purpose.

To gain insight into the potential of IgA antibodies as diagnostic biomarkers in SLE, we constructed ROC curves for IgA subtypes in the saliva of these patients. Our results indicate outstanding discriminative values for IgA1, and a lower but still good for IgA2 when distinguishing patients with SLE from healthy individuals.

As depicted, the sensitivity/specificity values given by measuring IgA isotypes in saliva regarding SLE diagnostic are comparable to those exhibited by other clinically relevant markers such as serum anti-dsDNA, anti-nucleosome, or complement levels in independent cohorts. This observation makes these secretory antibodies (particularly salivary IgA1 concentration) appealing candidates to be included as other criteria in SLE diagnostic current approaches upon prior validation in a larger multiethnic cohort.

Interestingly, when we performed the same analysis for assessing salivary IgA subtypes’ predictive ability when distinguishing inactive from active disease states, we could only obtain satisfactory AUC values for both antibody subtypes concentrations, making these two variables not as robust for activity discrimination. However, it becomes evident that including more patients, particularly those with active disease, could improve predictive values.

Again, as the measurement of these immunoglobulins supports the discrimination of SLE patients, their potential employment as a clinically useful biomarker for this disease becomes evident. One more time, our study remains limited due to the cross-sectional nature of our approaches; hence, we propose to conduct longitudinal monitoring of SLE patients to assess and/or validate the utility of IgA subtypes as new clinical biomarkers. Beyond that, these measurements should also be performed considering their predictive potential for different SLE implications, such as different types of flares.

Finally, as SLE patients can be treated with B cell-depleting biological therapeutics that consequently could diminish immunoglobulin levels, our proposed measurement of SIgA towards establishing a useful biomarker would be limited. The potential development of hypogammaglobulinemia is among the main concerns after administering biological treatments for SLE; for example, one of the most common drugs of this type employed for lupus treatment, rituximab, depletes CD20+ peripheral B cells for an average of 6-12 months, including naive and unswitched B cells, both of which are direct precursors for IgM production, thus further reducing the circulating levels of this isotype (54).

Different studies have informed that rituximab administration apparently does not affect the circulating IgA baseline levels in patients with SLE in the long term (55, 56). Conversely, other reports indicate only a slight initial decrease of IgA median levels that started recovering as early as two months after rituximab (54), a decrease that becomes significant upon cumulative cycles of treatment but only in a small percentage (around 3%) of patients (56). Supporting these data, it has previously been described that circulating IgA+ plasmablasts can persist early after rituximab, suggesting resistance to depletion of switched IgA+ precursor B cells, likely in the mucosal microenvironment and/or due to an early replenishment (57).

Currently, there is little information about secretory IgA levels upon biological treatment administration in contexts of autoimmune disease. However, one case report (58) indicates that rituximab treatment leads to an unexpected increase in the percentage of IgA+ plasmablasts in parotid salivary glands, compared with a biopsy before treatment.

All the previously mentioned data support the idea that biological treatments do not significantly affect the levels of IgA, even the secretory one in saliva; however, this hypothesis needs to be corroborated, including these treated patients in an independent study.

In conclusion, as the saliva sample is easily accessible and non-invasive for the patients, its IgA measurement emerges as an attractive alternative to being proposed as a novel biomarker of SLE. Given our results, salivary IgA subtypes correlate with specific autoantibodies related to disease diagnosis or activity but also may allow us to differentiate healthy individuals and SLE patients through the measurement of total salivary IgA or individual salivary subclasses, independently of antibody specificity.

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics statement

The studies involving human participants were reviewed and approved by Institutional Ethics and Research Committees of the Instituto Nacional de Ciencias Médicas y Nutrición Salvador Zubirán (Ref. 2306). The patients/participants provided their written informed consent to participate in this study.

Author contributions

SR-R, VS-H, and RC-D contributed to the design and performance of experiments, analysis, and interpretation of data. SR-R, VS-H, RC-D, DM-S, and CN-A performed experiments and analyzed data. DC-V and JT-R assisted in the processing and preservation of patient samples, collected patient data, and generated and organized our clinical database. JT-R and DG-M assisted in writing and editing the manuscript. SR-R, DG-M, and JM-M designed and performed experiments, supervised general work, wrote, and edited the manuscript. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by CONACyT [Grant A3-S-36875] and UNAM-DGAPA-PAPIIT Program [Grant IN212122].

Acknowledgments

SR-R is a doctoral student from the Programa de Doctorado en Ciencias Biomédicas, Universidad Nacional Autónoma de México (UNAM) and has received the Consejo Nacional de Ciencia y Tecnología (CONACyT) fellowship 736578/CVU 818733.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

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

Supplementary material

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

Abbreviations

AAAAI, The American Academy of Asthma, Allergy, and Immunology; ACR, American College of Rheumatology; ANAs, Anti-nuclear antibodies; ANCA, Antineutrophil cytoplasmic antibody; Anti-dsDNA, anti-double stranded DNA; aPL, anti-phospholipid; AUC, Areas under the curve; EULAR, European League Against Rheumatism; IgA, Immunoglobulin A; ROC, Receiver Operating Characteristic; SIgA, Secretory Immunoglobulin A; SLE, Systemic Lupus Erythematosus; SLEDAI, Systemic Lupus Erythematosus Activity Index; SLICC, Systemic Lupus International Collaborating Clinics; tIgA, total Immunoglobulin A.

References

1. Kett K, Brandtzaeg P, Radl J, Haaijman JJ. Different subclass distribution of IgA-producing cells in human lymphoid organs and various secretory tissues. J Immunol (1986) 136(10):3631–5. doi: 10.4049/jimmunol.136.10.3631

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Yaniv G, Twig G, Shor DB, Furer A, Sherer Y, Mozes O, et al. A volcanic explosion of autoantibodies in systemic lupus erythematosus: A diversity of 180 different antibodies found in SLE patients. Autoimmun Rev (2015) 14(1):75–9. doi: 10.1016/j.autrev.2014.10.003

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Erkeller-Yusel F, Hulstaart F, Hannet I, Isenberg D, Lydyard P. Lymphocyte subsets in a large cohort of patients with systemic lupus erythematosus. Lupus (1993) 2(4):227–31. doi: 10.1177/0961203393002001081

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Aringer M, Costenbader K, Daikh D, Brinks R, Mosca M, Ramsey-Goldman R, et al. 2019 European League against Rheumatism/American college of rheumatology classification criteria for systemic lupus erythematosus. Ann Rheum Dis (2019) 78(9):1151–9. doi: 10.1136/annrheumdis-2018-214819

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Tan G, Baby B, Zhou Y, Wu T. Emerging molecular markers towards potential diagnostic panels for lupus. Front Immunol (2021) 12:808839. doi: 10.3389/fimmu.2021.808839

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Barmettler S, Ong MS, Farmer JR, Choi H, Walter J. Association of immunoglobulin levels, infectious risk, and mortality with rituximab and hypogammaglobulinemia. JAMA Netw Open (2018) 1(7):e184169. doi: 10.1001/jamanetworkopen.2018.4169

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Perez EE, Orange JS, Bonilla F, Chinen J, Chinn IK, Dorsey M, et al. Update on the use of immunoglobulin in human disease: A review of evidence. J Allergy Clin Immunol (2017) 139(3S):S1–S46. doi: 10.1016/j.jaci.2016.09.023

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Buch MH, Smolen JS, Betteridge N, Breedveld FC, Burmester G, Dorner T, et al. Updated consensus statement on the use of rituximab in patients with rheumatoid arthritis. Ann Rheum Dis (2011) 70(6):909–20. doi: 10.1136/ard.2010.144998

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Bukhari M, Abernethy R, Deighton C, Ding T, Hyrich K, Lunt M, et al. BSR and BHPR guidelines on the use of rituximab in rheumatoid arthritis. Rheumatol (Oxford) (2011) 50(12):2311–3. doi: 10.1093/rheumatology/ker106a

CrossRef Full Text | Google Scholar

10. Smilek DE, Lim N, Ding L, Murray SG, Diamond B, Wofsy D. Correlation of hypogammaglobulinaemia with proteinuria, and the relationship between hypogammaglobulinaemia and infection in active lupus nephritis. Lupus Sci Med (2017) 4(1):e000229. doi: 10.1136/lupus-2017-000229

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Yates M, Watts RA, Bajema IM, Cid MC, Crestani B, Hauser T, et al. EULAR/ERA-EDTA recommendations for the management of ANCA-associated vasculitis. Ann Rheum Dis (2016) 75(9):1583–94. doi: 10.1136/annrheumdis-2016-209133

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Yong PF, Aslam L, Karim MY, Khamashta MA. Management of hypogammaglobulinaemia occurring in patients with systemic lupus erythematosus. Rheumatol (Oxford) (2008) 47(9):1400–5. doi: 10.1093/rheumatology/ken255

CrossRef Full Text | Google Scholar

13. Conley ME, Delacroix DL. Intravascular and mucosal immunoglobulin a: two separate but related systems of immune defense? Ann Intern Med (1987) 106(6):892–9. doi: 10.7326/0003-4819-106-6-892

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Fagarasan S, Honjo T. Intestinal IgA synthesis: regulation of front-line body defences. Nat Rev Immunol (2003) 3(1):63–72. doi: 10.1038/nri982

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Woof JM, Mestecky J. Mucosal immunoglobulins. Immunol Rev (2005) 206:64–82. doi: 10.1111/j.0105-2896.2005.00290.x

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Brandtzaeg P, Johansen FE. Mucosal b cells: Phenotypic characteristics, transcriptional regulation, and homing properties. Immunol Rev (2005) 206:32–63. doi: 10.1111/j.0105-2896.2005.00283.x

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Delacroix DL, Dive C, Rambaud JC, Vaerman JP. IgA subclasses in various secretions and in serum. Immunology (1982) 47(2):383–5.

PubMed Abstract | Google Scholar

18. Kutteh WH, Prince SJ, Hammond KR, Kutteh CC, Mestecky J. Variations in immunoglobulins and IgA subclasses of human uterine cervical secretions around the time of ovulation. Clin Exp Immunol (1996) 104(3):538–42. doi: 10.1046/j.1365-2249.1996.36742.x

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Ladjeva I, Peterman JH, Mestecky J. IgA subclasses of human colostral antibodies specific for microbial and food antigens. Clin Exp Immunol (1989) 78(1):85–90.

PubMed Abstract | Google Scholar

20. Mestecky J, Russell MW. IgA subclasses. Monogr Allergy (1986) 19:277–301.

PubMed Abstract | Google Scholar

21. Ahmadiafshar A, Mohsenifard MR, Mazloomzadeh S. Evaluation of serum & salivary IgA in patients with type 1 diabetes. PloS One (2015) 10(4):e0122757. doi: 10.1371/journal.pone.0122757

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Ghaleyani P, Sardari F, Akbari M. Salivary IgA and IgG in oral lichen planus and oral lichenoid reactions diseases. Adv BioMed Res (2012) 1:73. doi: 10.4103/2277-9175.102977

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Hung YH, Lee YH, Chen PP, Lin YZ, Lin CH, Yen JH. Role of salivary immune parameters in patients with primary sjogren’s syndrome. Ann Lab Med (2019) 39(1):76–80. doi: 10.3343/alm.2019.39.1.76

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Patidar KA, Parwani RN, Wanjari SP. Correlation of salivary and serum IgG, IgA levels with total protein in oral submucous fibrosis. J Oral Sci (2011) 53(1):97–102. doi: 10.2334/josnusd.53.97

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Perez MM, Pessoa JS, Ciamponi AL, Diniz MB, Santos M, Alves HHO, et al. Correlation of salivary immunoglobulin a with body mass index and fat percentage in overweight/obese children. J Appl Oral Sci (2018) 27:e20180088. doi: 10.1590/1678-7757-2018-0088

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Roos Ljungberg K, Borjesson E, Martinsson K, Wettero J, Kastbom A, Svard A. Presence of salivary IgA anti-citrullinated protein antibodies associate with higher disease activity in patients with rheumatoid arthritis. Arthritis Res Ther (2020) 22(1):274. doi: 10.1186/s13075-020-02363-0

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Trull AK, Panayi GS. Serum and secretory IgA immune response to klebsiella pneumoniae in ankylosing spondylitis. Clin Rheumatol (1983) 2(4):331–7. doi: 10.1007/BF02041551

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Petri M, Orbai AM, Alarcon GS, Gordon C, Merrill JT, Fortin PR, et al. Derivation and validation of the systemic lupus international collaborating clinics classification criteria for systemic lupus erythematosus. Arthritis Rheumatol (2012) 64(8):2677–86. doi: 10.1002/art.34473

CrossRef Full Text | Google Scholar

29. Sanchez-Salguero E, Mondragon-Ramirez GK, Alcantara-Montiel JC, Cerbulo-Vazquez A, Villegas-Dominguez X, Contreras-Vargas VM, et al. Infectious episodes during pregnancy, at particular mucosal sites, increase specific IgA1 or IgA2 subtype levels in human colostrum. Matern Health Neonatol Perinatol. (2019) 5:9. doi: 10.1186/s40748-019-0104-x

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Brandtzaeg P. Secretory immunity with special reference to the oral cavity. J Oral Microbiol (2013) 5:1–24. doi: 10.3402/jom.v5i0.20401

CrossRef Full Text | Google Scholar

31. Putterman C, Furie R, Ramsey-Goldman R, Askanase A, Buyon J, Kalunian K, et al. Cell-bound complement activation products in systemic lupus erythematosus: comparison with anti-double-stranded DNA and standard complement measurements. Lupus Sci Med (2014) 1(1):e000056. doi: 10.1136/lupus-2014-000056

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Brunekreef T, Limper M, Melchers R, Mathsson-Alm L, Dias J, Hoefer I, et al. Microarray testing in patients with systemic lupus erythematosus identifies a high prevalence of CpG DNA-binding antibodies. Lupus Sci Med (2021) 8(1):e000531. doi: 10.1136/lupus-2021-000531

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Kerr MA. The structure and function of human IgA. Biochem J (1990) 271(2):285–96. doi: 10.1042/bj2710285

PubMed Abstract | CrossRef Full Text | Google Scholar

34. da Silva LS, Almeida BL, de Melo AK, de Brito DC, Braz AS, Freire EA. IgA nephropathy in systemic lupus erythematosus patients: case report and literature review. Rev Bras Reumatol Engl Ed (2016) 56(3):270–3. doi: 10.1016/j.rbr.2014.10.011

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Peen E, Mellbye OJ, Haga HJ. IgA rheumatoid factor in primary sjogren’s syndrome. Scand J Rheumatol (2009) 38(1):46–9. doi: 10.1080/03009740802366043

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Van Hoovels L, Vander Cruyssen B, Sieghart D, Bonroy C, Nagy E, Pullerits R, et al. IgA rheumatoid factor in rheumatoid arthritis. Clin Chem Lab Med (2022) 60(10):1617–26. doi: 10.1515/cclm-2022-0244

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Song DJ, Shen J, Chen MH, Liu ZJ, Cao Q, Hu PJ, et al. Association of serum immunoglobulins levels with specific disease phenotypes of crohn’s disease: A multicenter analysis in China. Front Med (Lausanne) (2021) 8:621337. doi: 10.3389/fmed.2021.621337

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Cassidy JT, Kitson RK, Selby CL. Selective IgA deficiency in children and adults with systemic lupus erythematosus. Lupus (2007) 16(8):647–50. doi: 10.1177/0961203307077543

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Ben-Aryeh H, Gordon N, Szargel R, Toubi E, Laufer D. Whole saliva in systemic lupus erythematosus patients. Oral Surg Oral Med Oral Pathol (1993) 75(6):696–9. doi: 10.1016/0030-4220(93)90425-4

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Male CJ. Immunoglobulin A1 protease production by haemophilus influenzae and streptococcus pneumoniae. Infect Immun (1979) 26(1):254–61. doi: 10.1128/iai.26.1.254-261.1979

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Janoff EN, Rubins JB, Fasching C, Charboneau D, Rahkola JT, Plaut AG, et al. Pneumococcal IgA1 protease subverts specific protection by human IgA1. Mucosal Immunol (2014) 7(2):249–56. doi: 10.1038/mi.2013.41

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Jin S, Yu C, Yu B. Changes of serum IL-6, IL-10 and TNF-alpha levels in patients with systemic lupus erythematosus and their clinical value. Am J Transl Res (2021) 13(4):2867–74.

PubMed Abstract | Google Scholar

43. Godsell J, Rudloff I, Kandane-Rathnayake R, Hoi A, Nold MF, Morand EF, et al. Clinical associations of IL-10 and IL-37 in systemic lupus erythematosus. Sci Rep (2016) 6:34604. doi: 10.1038/srep34604

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Facciotti F, Larghi P, Bosotti R, Vasco C, Gagliani N, Cordiglieri C, et al. Evidence for a pathogenic role of extrafollicular, IL-10-producing CCR6(+)B helper T cells in systemic lupus erythematosus. Proc Natl Acad Sci U S A. (2020) 117(13):7305–16. doi: 10.1073/pnas.1917834117

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Park YB, Lee SK, Kim DS, Lee J, Lee CH, Song CH. Elevated interleukin-10 levels correlated with disease activity in systemic lupus erythematosus. Clin Exp Rheumatol (1998) 16(3):283–8.

PubMed Abstract | Google Scholar

46. He B, Xu W, Santini PA, Polydorides AD, Chiu A, Estrella J, et al. Intestinal bacteria trigger T cell-independent immunoglobulin A(2) class switching by inducing epithelial-cell secretion of the cytokine APRIL. Immunity (2007) 26(6):812–26. doi: 10.1016/j.immuni.2007.04.014

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Lin M, Du L, Brandtzaeg P, Pan-Hammarstrom Q. IgA subclass switch recombination in human mucosal and systemic immune compartments. Mucosal Immunol (2014) 7(3):511–20. doi: 10.1038/mi.2013.68

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Marques CP, Victor EC, Franco MM, Fernandes JM, Maor Y, de Andrade MS, et al. Salivary levels of inflammatory cytokines and their association to periodontal disease in systemic lupus erythematosus patients. A case-control study. Cytokine (2016) 85:165–70. doi: 10.1016/j.cyto.2016.06.025

CrossRef Full Text | Google Scholar

49. Linnik MD, Hu JZ, Heilbrunn KR, Strand V, Hurley FL, Joh T, et al. Relationship between anti-double-stranded DNA antibodies and exacerbation of renal disease in patients with systemic lupus erythematosus. Arthritis Rheumatol (2005) 52(4):1129–37. doi: 10.1002/art.20980

CrossRef Full Text | Google Scholar

50. Zeng Y, Xiao Y, Zeng F, Jiang L, Yan S, Wang X, et al. Assessment of anti-nucleosome antibody (ANuA) isotypes for the diagnosis and prediction of systemic lupus erythematosus and lupus nephritis activity. Clin Exp Med (2022). doi: 10.1007/s10238-022-00942-w

CrossRef Full Text | Google Scholar

51. Villalta D, Bizzaro N, Bassi N, Zen M, Gatto M, Ghirardello A, et al. Anti-dsDNA antibody isotypes in systemic lupus erythematosus: IgA in addition to IgG anti-dsDNA help to identify glomerulonephritis and active disease. PloS One (2013) 8(8):e71458. doi: 10.1371/journal.pone.0071458

PubMed Abstract | CrossRef Full Text | Google Scholar

52. Pradhan VD, Patwardhan MM, Ghosh K. Anti-nucleosome antibodies as a disease marker in systemic lupus erythematosus and its correlation with disease activity and other autoantibodies. Indian J Dermatol Venereol Leprol. (2010) 76(2):145–9. doi: 10.4103/0378-6323.60558

PubMed Abstract | CrossRef Full Text | Google Scholar

53. Li T, Prokopec SD, Morrison S, Lou W, Reich H, Gladman D, et al. Anti-nucleosome antibodies outperform traditional biomarkers as longitudinal indicators of disease activity in systemic lupus erythematosus. Rheumatol (Oxford) (2015) 54(3):449–57. doi: 10.1093/rheumatology/keu326

CrossRef Full Text | Google Scholar

54. Marco H, Smith RM, Jones RB, Guerry MJ, Catapano F, Burns S, et al. The effect of rituximab therapy on immunoglobulin levels in patients with multisystem autoimmune disease. BMC Musculoskelet Disord (2014) 15:178. doi: 10.1186/1471-2474-15-178

PubMed Abstract | CrossRef Full Text | Google Scholar

55. Aguiar R, Araujo C, Martins-Coelho G, Isenberg D. Use of rituximab in systemic lupus erythematosus: A single center experience over 14 years. Arthritis Care Res (Hoboken) (2017) 69(2):257–62. doi: 10.1002/acr.22921

PubMed Abstract | CrossRef Full Text | Google Scholar

56. Almaghlouth I, Johnson SR, Pullenayegum E, Gladman D, Urowitz M. Immunoglobulin levels in systemic lupus erythematosus: A narrative review. Lupus (2021) 30(6):867–75. doi: 10.1177/09612033211004714

PubMed Abstract | CrossRef Full Text | Google Scholar

57. Cambridge G, Isenberg DA, Edwards JC, Leandro MJ, Migone TS, Teodorescu M, et al. B cell depletion therapy in systemic lupus erythematosus: Relationships among serum b lymphocyte stimulator levels, autoantibody profile and clinical response. Ann Rheum Dis (2008) 67(7):1011–6. doi: 10.1136/ard.2007.079418

PubMed Abstract | CrossRef Full Text | Google Scholar

58. Pijpe J, van Imhoff GW, Vissink A, van der Wal JE, Kluin PM, Spijkervet FK, et al. Changes in salivary gland immunohistology and function after rituximab monotherapy in a patient with sjogren’s syndrome and associated MALT lymphoma. Ann Rheum Dis (2005) 64(6):958–60. doi: 10.1136/ard.2004.030684

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: IgA, IgA1, IgA2, systemic lupus erythematosus, saliva

Citation: Romero-Ramírez S, Sosa-Hernández VA, Cervantes-Díaz R, Carrillo-Vázquez DA, Meza-Sánchez DE, Núñez-Álvarez C, Torres-Ruiz J, Gómez-Martín D and Maravillas-Montero JL (2023) Salivary IgA subtypes as novel disease biomarkers in systemic lupus erythematosus. Front. Immunol. 14:1080154. doi: 10.3389/fimmu.2023.1080154

Received: 25 October 2022; Accepted: 13 February 2023;
Published: 22 February 2023.

Edited by:

Renato C. Monteiro, Université de Paris, France

Reviewed by:

Klara Martinsson, Linköping University, Sweden
Zhiqiu Xia, University of Nebraska Medical Center, United States

Copyright © 2023 Romero-Ramírez, Sosa-Hernández, Cervantes-Díaz, Carrillo-Vázquez, Meza-Sánchez, Núñez-Álvarez, Torres-Ruiz, Gómez-Martín and Maravillas-Montero. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: José L. Maravillas-Montero, maravillas@cic.unam.mx

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