- 1Immunology & Microbiology, National Institute for Occupational Health, National Health Laboratory Service, Johannesburg, South Africa
- 2Department of Environmental Health, Faculty of Health Sciences, University of Johannesburg, Johannesburg, South Africa
- 3Department of Clinical Microbiology and Infectious Diseases, School of Pathology, University of the Witwatersrand, Johannesburg, South Africa
- 4Immunology & Microbiology Department, Centre for Statistical Analysis and Research, Johannesburg, South Africa
- 5Division of Occupational Medicine and Centre for Environmental & Occupational Health Research, School of Public Health and Family Medicine, University of Cape Town, Cape Town, South Africa
Cytokines elicit a pro-inflammatory response by modifying the airway microenvironment in patients with acute or chronic asthma. The expression pattern of several distinct cytokines could be a useful discriminator in asthma. This study aimed to identify asthma subject groupings based on common inflammatory patterns and to determine the relationship between these identified patterns and asthma-associated clinical indices. A sub-group of 76 dental healthcare workers (HCWs) identified from a larger cross-sectional study of 454 dental HCWs in five dental institutions were evaluated further. A self-administered questionnaire elicited the health and employment history of subjects. Sera were analyzed for atopic status, latex sensitization, and 12 cytokines (IL-1β, 3, 4, 5, 6, 7, 8, 10, 12p70, eotaxin, GM-CSF, TNF-α). Pre and post-bronchodilator spirometry was performed on all HCWs. Data clustering and factor analysis were used to identify inflammatory cluster patterns of cytokines. Associations between the cytokine cluster groupings and relevant asthma-associated clinical indices were determined using multivariate logistic regression. The classification of asthma subtype based on cytokine patterns demonstrated both eosinophilic and neutrophilic inflammatory responses. Four phenotypically distinct subgroups relating to the severity of inflammation (acute or chronic) of the cell types were identified. Cytokine determinants for the neutrophilic subtype included IL-1β, 6, 8, 10, 12p70, and TNF-α whereas for the eosinophilic subtype these included IL-3, 4, 5, 7, eotaxin, and GM-CSF. The multivariate models showed a significant association between work-related chest symptoms and all four inflammatory patterns. However, stronger associations were observed for the acute neutrophilic (OR = 6.00, p < 0.05) compared to acute and chronic eosinophilic responses (OR = 4.30, p < 0.05; OR = 4.93, p < 0.05), respectively. Subjects with airway obstruction were more likely to have a mixed cellular infiltrate. The odds of work-exacerbated asthma were increased in acute or chronic eosinophilia (OR = 7.75 and 8.12; p < 0.05), respectively as well as with acute neutrophilia (OR = 6) sub-type. This study demonstrated that neutrophilic inflammatory cell asthma phenotypes coexist with eosinophilic inflammatory phenotypes suggesting a possible dual pathway for asthma in dental health workers, probably due to mixed exposures to high molecular weight (e.g., latex) and low molecular weight (e.g., acrylates) agents.
Introduction
Cytokines play an important role in the inflammatory response of asthmatic airways since they are capable of inducing many pro-inflammatory effects characteristic of the disease (1). These molecules can modify the airway microenvironment directly or by activating resident cells (2). The implementation of tailored medicine requires the identification of endotype-specific markers in biological matrices. In this regard, allergy and asthmatic conditions are well-suited, as they present clinically with similar visible properties (phenotypes) and partially share immune mechanisms (endotypes) (3), thus providing a better understanding of risk prediction, personalized or potentially more effective treatment selection or mechanism-specific prevention strategies (4, 5). Although two major groups of cytokines have been identified based on the pathogenesis of asthma (3), these cytokines cannot be easily categorized in acute and chronic asthma due to their pleiotropic nature and overlapping properties (6). Moreover, cytokine expression patterns consisting of several distinct cytokines are more useful discriminators in multifactorial diseases like asthma compared to single cytokine levels (3, 7).
Asthma has historically been regarded as entirely an atopic disease involving allergic (IgE-mediated) sensitization resulting in interleukin (IL)-5 mediated eosinophilic airways inflammation, which was thought to be a cardinal feature (8). However, the role of IL-5 in asthma has been under scrutiny as studies have demonstrated that only 50% of asthma cases can be attributable to “allergic asthma” (9). This attribution is based on inconsistent associations between different asthma phenotypes possibly due to exposure to high molecular weight (HMW) and low molecular weight (LMW) agents (10). In addition, treatment regimens aimed at neutralizing immunoglobulin E (IgE) (anti-IgE) have only been partially effective in asthma cases. The phenomenon of non-eosinophilic inflammation (referred to as neutrophilic inflammation) is prevalent in certain occupational populations (9). These populations include those with exposure to metal-working fluids (11), LMW sensitizers usually chemical agents (12) such as isocyanates and cobalt (13), acrylates in the dentistry (10, 14, 15) and endotoxin exposure in farming (16). Although not demonstrated in working populations, bakery flour dust has also been shown to elicit a neutrophilic inflammation in mouse models (17). Asthma phenotypes can therefore be classified based on clinical or physiological categories (i.e., severity, age at onset, and chronic airflow obstruction), asthma triggers (i.e., exercise, allergens, occupational allergens or irritants) or partially sharing immune mechanisms (endotypes) of predominant inflammatory cell infiltrates (i.e., eosinophilic or neutrophilic asthma) (18). The latter pathobiological phenotypes include eosinophilic asthma, neutrophilic asthma, granulocytic asthma with concomitant eosinophilic, and neutrophilic inflammation, and a group, paucigranulocytic asthma, in which levels of both neutrophils and eosinophils are normal (9, 19, 20).
The measurement of sputum eosinophilic cationic protein (ECP) has been suggested as a predictive marker for the presence of eosinophilic inflammation in the lung (21) and changes in sputum eosinophilia also reflect fluctuations in asthma control (3). Furthermore, ECP is more predictive of sputum eosinophilia than serum IgE (3). Hancox et al., demonstrated that high blood eosinophil counts were associated with airflow obstruction (lower forced expiratory volume in one second [FEV1] % predicted values for both pre- and post-bronchodialator spirometry, p ≤ 0.048) (22). The pathobiological effects of eosinophils occur through a myriad of mediators including granulocyte-macrophage colony-stimulating factor (GM-CSF), IL-3, IL-5, and certain chemokines such as eotaxin (3, 6, 23, 24). Neutrophils also play a role in the development of well-defined clinical sub-phenotypes (25–30) and are a potential source of inflammatory mediators, enzymes (e.g., neutrophil elastase, myeloperoxidase [MPO]), reactive oxygen species, and cytokines such as IL-8 and tumor necrosis factor alpha (TNF-α). Interleukin-8 is considered a marker of severe asthma and has been linked to allergic rhinitis (27, 31).
The use of biomarkers in the assessment of airway inflammation is important for exploring the underlying mechanisms of asthma and may reveal diseased airways that are not detectable by parameters commonly used in current practice (e.g., symptom experience, lung function, airway hyperresponsiveness) (32). Airway-infiltrating cells like eosinophils modify airway responses and have profound biologic effects in the micro-environment of their release as well as at distant sites. However, they may no longer be present at the time that the effects of the cytokine/s are apparent, which has probably contributed to discordance in results reported in previous studies. The discordance may also be attributed to the heterogeneity of asthma (7, 19, 33), highlighting the need to make inferences about the relationship between airway hyperresponsiveness and airway inflammation to be made with circumspection (34). A potential advantage of using inflammatory cell markers is the ability to distinguish between heterogeneous subtypes in subjects with similar clinical characteristics (7). This study hypothesized that different cytokine groupings could be used to recognize different inflammatory cell endotypes in asthmatics and these groupings could be used to examine the heterogeneity of airway inflammation in dental healthcare workers (HCWs). The aim of the study was to identify asthma subject groupings based on common inflammatory patterns and to determine the relationship between these identified patterns and asthma-associated clinical indices.
Materials and Methods
Study Population
This study was part of a larger cross-sectional study of 454 employed dental HCWs (34). Subjects were assessed using a standardized protocol that included an interviewer-administered questionnaire, pulmonary function tests, immunological assessment for atopy and sensitization to occupational allergens (latex, chloramines). Further detailed information as to the methods employed in identifying individuals with asthma are described in this communication. In brief, the questionnaire collected a history of previous medical illnesses and other symptom variables: respiratory (wheeze and/or tight chest); ocular (itchy eyes, red eyes); and nasal (runny nose, blocked nose/stuffy nose) symptoms. The current study population comprised 76 asthmatic subjects, that were identified using epidemiological definitions of asthma as either having doctor-diagnosed asthma, probable asthma or those with current asthma-related symptoms (35). Probable asthma was defined as having reported either an attack of shortness of breath, attack of asthma or using asthma medication. Those with current asthma-related symptoms included those with a positive response to one or more of the following: shortness of breath, chest tightness, coughing, and wheezing in the last 12 months. Health outcome information was obtained from a self-administered respiratory health questionnaire based on the European Community Respiratory Health Survey (ECRHS) and spirometry indices obtained according to the American Thoracic Society (ATS) guidelines (36) as described previously (34). A positive bronchodilator test entailed an increase in FEV1 and/or forced vital capacity (FVC) ≥12% and ≥200 ml increase from baseline and is described below. Control subjects (n = 9) were classified as those with a negative response to either doctor-diagnosed asthma, probable asthma or current asthma-related symptoms.
Serum Samples
A 6 ml blood sample was taken from all workers using a Becton Dickinson Vacutainer CAT tube (with clot activator). Samples were left to coagulate for 3–6 h at room temperature and were then centrifuged at 1,500 g for 10 min. The serum was transferred to 5 ml plastic cryotubes (Nunc, Denmark) and stored at −70°C until analysis.
Measurement of Atopy and Allergen-Specific IgE
The serum was analyzed for atopy and probable latex allergy using the Phadiatop® test for common inhalants and latex-specific IgE ImmunoCAP (k82), respectively. Tests were done according to the manufacturer's instructions using the ImmunoCAP system 100 (Phadia AB, Uppsala, Sweden). Specific IgE to cross-reacting carbohydrate determinants (CCDs) and horse-radish peroxidase (HRP) and bromelain (MUXF3) were also performed. A result of >0.35 kU/l was regarded as positive for allergen-specific IgE determinations.
Measurement of Eosinophilic Cationic Protein (ECP)
Serum ECP was measured using the Phadia ImmunoCAP system ECP FEIA (Phadia AB, Uppsala, Sweden), according to the manufacturer's instructions. Briefly, anti-ECP covalently coupled to ImmunoCAP, binds with the ECP in the patient's serum. After washing, enzyme-labeled antibodies against ECP were added to form a complex. After incubation, unbound enzyme-anti-ECP was washed away and the bound complex was then incubated with a developing solution. The reaction was stopped, the fluorescence of the eluate was measured and the response transformed to concentration levels using a calibration curve. The normal range for serum ECP in the patient population that was established in the laboratory was 2.00–15.7 μg/l. Therefore, values above this were considered elevated for the study population. The maximum value was similar to that of a previous study (13.3 μg/l) (37). To ensure quality control the laboratory participated in an external quality assessment program (Phadia: Quality Club ECP [52-5213-99/05]).
Measurement of Myeloperoxidase (MPO)
MPO concentration in serum was measured using the MPO enzyme immunoassay (MPO-EIA) (MPO-EIA, Oxford Biomedical Research), following the manufacturer's instructions (Oxford Biomedical Research). Briefly, antigen captured by a solid phase monoclonal antibody was detected with a biotin-labeled goat polyclonal anti-MPO. An avidin alkaline phosphatase conjugate was bound to the biotinylated antibody. The alkaline phosphatase substrate ρ-nitrophenyl phosphate (pNPP) together with the yellow product (ρ-nitrophenol) was read at 405 nm using a spectrophotometer (Biotek Instruments, United States of America). A standard curve was obtained by plotting the absorbance values as a function of the standard MPO concentrations. Standard curves were fitted using a four-parameter data reduction method. Assays in which the standard curve had a correlation coefficient (R2) ≥ 0.98 and samples with a coefficient of variation (CV) ≤ 10% were accepted. The normal range for serum in the patient population that was established in the laboratory was 8.43–77.35 ng/ml and values above this were considered elevated for the study population.
Measurement of Non-allergic and Allergic Inflammatory Mediators (Cytokines)
The BDTM CBA Human Inflammation kit was used to quantitatively measure non-allergic inflammatory cytokines (interleukin (IL)-8, IL-1β, IL-6, IL-10, TNF-α, and IL-12p70) in each serum sample. The assay allowed for multiplexed analysis of multiple proteins from a single sample. Six bead populations with distinct fluorescence intensities coated with the specific ILs were multiplexed and resolved in the red channel of the flow cytometer (BD FacsArrayTM) (Becton Dickson, Oxford, United Kingdom). The determination of the serum inflammatory ILs was done according to the kit's instructions. Standard curves of the standard serial dilutions, expressed by a four-parameter logistic model (log CC = D + (A – D)/(1 + (log I/C)∧B)), where A = minimum asymptote, B = slope factor, C = inflection point, D = maximum asymptote, were used to determine the limit of detection (LOD) for each specific analyte based on the average median fluorescent intensity (MFI) of the negative control. The assay detection limit for each analyte was determined by the average fitted concentration of the negative control (0 pg/mL) + 2 standard deviations for each analyte. The concentration of analyte was expressed as picograms per milliliter (pg/mL). A regression coefficient (R2) for the standard curve for each analyte was accepted if ≥0.98. The results were interpreted by the level of expression of the cytokine with categories of low expression to highly expressed cytokines (0–2.9, >2.9–4.5, >4.5–6.5, >6.5–9.4, >9.4–13.2, >13.2–28, >28–80.4, >80.4–429 pg/ml). The allergy-related inflammatory cytokines (IL-7, IL-3, IL-4, IL-5, eotaxin, and GM-CSF) were measured using the BDTM CBA Human Soluble Protein Flex Set assay. The multiplexed bead populations were resolved in the near-infrared (NIR) and red channels of the flow cytometer (BD FacsArrayTM) (Becton Dickson, Oxford, UK). The procedure was done according to the kit's instructions. The results were interpreted as described for the non-allergic cytokine mediators above.
Spirometry and Bronchodilator Challenge (BDR)
Spirometry was conducted according to the American Thoracic Society (ATS) guidelines using a Jaeger Masterscope spirometer (Höchberg, Germany), calibrated at least twice a day with a 3-L syringe. Reproducibility criteria included the two best tracings for both FEV1 and FVC varying by no more than 0.2 L (38). Measurements obtained by spirometry were adjusted according to the temperature and atmospheric pressure measured throughout the day. Special instructions were given to workers to ensure that they did not smoke tobacco (at least 1 h) and did not use any anti-asthmatic inhalers (4 h) or take oral asthma medications (8 h) before the test. Workers' heights were recorded for calculating predicted lung function indices using reference values of the European Community for Coal and Steel with lower limits corresponding to the 95th percentile (39). The post-bronchodilator test was done at least 10 min after administration of the bronchodilator (inhaled beta-2 agonists, salbutamol, 400 mg). An improvement of at least 12% and 200 ml after administration of salbutamol was regarded as a significant response. Both pre- and post bronchodilator FEV1 had to be reproducible for calculating the degree of the bronchodilator response.
Statistical Analysis
Variables were log-transformed to meet the model assumptions of normality and parametric analyses were then applied. However, where certain variables followed a non-guassian distribution after transformation, non-parametric analyses were used. The results of the continuous variables were expressed as means and standard deviations. Chi-squared statistics were applied for binary comparisons. The grouping of subjects based on their cytokine profiles was achieved by both agglomerative hierarchical complete linkage clustering (40) and partitional clustering to compare how well the algorithms performed on the cytokine measurements. For the former technique, the identification was based on the largest stopping rule index (Calinski/Harabasz pseudo-F = 118.67). The study sample was also divided into four clusters using partition clustering (kmedians), based on previously known inflammatory cell asthma phenotypes (neutrophilic, eosinophilic, granulocytic and paucigranulocytic) (18). Factor analysis was used to determine inflammatory patterns based on the 12 cytokines (41, 42). Principal factor analysis (PFA), as well as principal component factor analysis (PCA), was done to determine the inflammatory cell phenotypes. The appropriateness of the correlation matrix (0.015) for factor analysis was assessed by the Bartlett's test of sphericity (p < 0.001) and the Kaiser-Meyer-Olkin test (0.75). The number of factors was determined by the Eigenvalue criterion, whose magnitude should be ≥1, as well as the percentage of total variance explained by the factors. To increase the interpretability, the factors were rotated orthogonally through the Varimax rotation procedure. The following criteria were used in labeling the factors: i) a factor loading had to be >0.45; ii) variables that loaded >0.45 on more than one factor were excluded, and iii) single factor loadings were not accepted (43). The naming of factors was based on the general association between the various cytokines and severity of asthma cited in the literature (44–46). To test for the difference between cluster groupings by phenotypic characteristics, analysis of variance (ANOVA) (with Bonferroni correction) was used. This method is similar to bi-clustering where not only the objects are clustered but also the features of the objects (i.e., the rows and columns are clustered simultaneously) (47). To further determine the association between the cluster groupings and cytokine patterns with asthma-associated clinical indices, multivariate logistic regression analysis, adjusting for age and gender, was applied. Statistical analysis was performed using STATA9 computer software (StataCorp, 2007, Texas, USA).
Results
Demographic Characteristics
The age distribution between the asthmatic subjects and controls was similar (Table 1). There was a larger proportion of females than males in both groups. The asthmatic workers had a marginally higher average body mass index (BMI), placing them in the overweight category, compared to the control group that was borderline normal. The majority of subjects from both groups were non-smokers. The control group had relatively better lung function compared to the asthmatic group, although this was not statistically significant. Just over half of the asthmatic subjects were atopic (59%) and ~ a tenth (9%) were sensitized to latex. The two cross-reacting carbohydrate determinants (CCDs), horse-radish peroxidase (HRP) and bromelain (MUXF3) were measured in all latex sensitized asthmatic subjects to exclude cross-reactivity. Furthermore, both serum inflammatory biomarkers (ECP, MPO) were relatively higher in the asthmatic group whereas the cytokine levels were variable between the two groups.
Cytokine Groupings Defined by Cluster Analysis
The expressions of 12 cytokines are shown in the heat map, which illustrates the variability of the measured cytokines in the asthmatic study population (n = 76) (Figure 1). The agglomerative clustering analysis procedure identified four potential clusters within the population. However, the agglomerative clustering produced background noise making it difficult to distinguish between groups. The population was therefore partitioned into four discernable groups using partition clustering namely cluster-1, cluster-2, cluster-3, and cluster-4. This provided some optimality criterion and therefore the clusters identified by this method were used for further analysis in this investigation.
Figure 1. Heat map showing the molecular profile of 12 inflammatory markers concentrations in asthmatic dental workers (cases). Each value is represented by a colored rectangular block. The high expressing cytokines (pg/ml) are depicted by the darker color shades and the low expressing cytokines by the lighter shades. Each row is an individual dental worker and the individual cytokine values are shown in each column. To lhe right of the heatmap, is the dendrogram showing the distinction between groups.
Identification of Inflammatory Patterns
The correlation matrix of cytokines was first subjected to PFA for the 76 subjects (Figure 2). A 2-factor model was confirmed based on the Eigen values for factor-1 (F1) and factor-2 (F2) with 2.83 and 1.95, respectively. The following cytokines met the factor selection (correlation (ρ) > 0.45), F1: IL-5, eotaxin, GM-CSF, IL-3 and F2: IL-1β, IL-10, TNF-α, and IL-12p70. The exception was eotaxin which featured in F1 but did not meet the selection criteria when Varimax rotation was applied (ρ = 0.43). Further analysis was done to extract the highest possible variance in the dataset using PCA and Varimax rotation. A four-factor (F1, F2, F3, and F4) solution was found to represent the data based on Eigenvalue criteria, explaining ~ 67% of the variance in the data set (p < 0.001). The Varimax rotation demonstrated unique factor loadings in F3 for IL-4 and IL-7, which were excluded from the initial PFA (Table 2). The labeling of the four factors identified by Varimax rotation was based on the factor selection criteria and the relation of the inflammatory patterns in asthma (Table 2).
Figure 2. Factor loading plot for principal factor analysis showing factor groupings with and without Varimax rotation. Solid line denotes variable loadings included in factors without rotations and dashed lines denote variable loadings included in factors with Varimax rotation.
Table 2. Principal component factor loading with Varimax rotation and factor labeling according to the inflammatory response of the cytokine groupings in dental healthcare workers.
Relationship Between Cluster Groupings and Inflammatory Patterns
To determine whether there were differences between the four inflammatory patterns and the four cluster groupings an ANOVA test was conducted (Table 3). All the cytokines that scored high in F1 (chronic eosinophilia) demonstrated the highest mean levels in subjects from cluster-1. Interleukin-6 and IL-8, which were highly expressed in F4 (chronic neutrophilia) showed the highest mean levels for subjects in cluster-3. There was a significant variance in the cytokine measurements featured in F2 and F4 and the cluster groupings. The majority of cytokines that scored high in F2 (acute neutrophilia) were highly expressed in cluster-4. These results demonstrated that there are distinct differences between the clusters and the majority of the factors.
Table 3. Analysis of variance between cytokine measurements of the identified factors, by the four cluster groupings in dental healthcare workers.
Relationship Between Serum Inflammatory Markers and Individual Cytokines
The relationship between serum inflammatory markers (ECP and MPO) and individual cytokines determined are presented in Table 4. A significant positive correlation was demonstrated between IL-8 and both ECP (r = 0.42, p < 0.001) and MPO (r = 0.47, p < 0.05). Similarly, a modest correlation was found between allergic inflammatory cytokines (IL-5, IL-7) and MPO. An inverse correlation was demonstrated between IL-1β and MPO (r = −0.40, p < 0.05).
Table 4. Correlation between log transformed serum inflammatory markers (ECP and MPO) and cytokine concentrations for the various factors identified in dental healthcare workers.
Relationship Between Inflammatory Patterns and Asthma-Associated Clinical Indices
Examination of the asthma-associated clinical indices and cytokine inflammatory patterns revealed a significant association between having work-related chest symptoms and all four inflammatory patterns (Table 5). Subjects with airway obstruction were significantly (p < 0.05) more likely to have a mixed cellular infiltrate of acute and chronic eosinophilia and acute neutrophilia. The odds of having increased MPO were higher in subjects with either acute and chronic neutrophilic inflammation, although this was not significant. The frequency of asthma attacks appeared to be more common in subjects with acute or chronic eosinophilia (OR range: 7.42–9.49) in the unadjusted model (data not shown). The odds of having atopic asthma were 10-fold higher in subjects with a chronic eosinophilic response and appeared to be similar in both acute eosinophilic and neutrophilic responses (OR = ~5). Furthermore, the odds of WEA were similar in subjects with acute or chronic eosinophilia (OR = 7.75 and 8.12; p < 0.05), respectively as well as in acute neutrophilia (OR = 6). Similar observations were demonstrated in other unadjusted logistic models (data not shown).
Table 5. Association between asthma related clinical indices and inflammatory phenotypes identified in dental healthcare workers (n = 76) in multivariate logistic regression models.
Discussion
This study demonstrated four cluster groupings based on differing cytokine patterns in dental health workers with asthma-associated clinical indices, resulting in multiple inflammatory phenotypes. The asthma subtypes based on cytokine measurements entailed two related to an eosinophilic and the other two to a neutrophilic response. Each cellular response also appeared to be related further to the degree of inflammation (i.e., acute or chronic) present. Similar results have been shown by Brasier et al., in which at least four phenotypically distinct subgroups were identified based on bronchoalveolar lavage-derived cytokines for severe asthma (6). Shared variance among clusters and cross-loadings of factors also suggested that some cytokines may overlap between groups or represent aetiological pathways that have smaller effects (48–50). Lemiere et al., identified three clusters in occupational asthmatics based on increases in fractional exhaled nitric oxide (FeNO) following challenge to causal agents, two clusters showing an increase post challenge to high molecular weight (HMW) agents, the two differing in degrees of severity of response, and a third cluster presenting with late reactions due to LMW agents (50). Vandenplas et al., also demonstrated differences in work-related asthma symptom patterns since subjects with OA due to LMW agents were more likely to experience chest tightness at work and daily sputum production (10). The results of the current study confirm the utility of the factor analysis (43) approach in identifying the various inflammatory patterns of asthma in dental HCWs.
Soluble markers such as ECP and MPO are generally used to assess the presence of eosinophils and neutrophils in body fluids, respectively (32, 51, 52). In the current study, a significant correlation was observed between cytokines generally associated with chronic neutrophilia with both ECP and MPO levels (53). Furthermore, a significant correlation was found between acute as well as chronic eosinophilic cytokines and MPO. A combination of cellular types (54) generally reflects an inflammatory response to multiple agents characterized by dual inflammatory pathways as a result of mixed exposures. The dual response could either indicate independent pathways caused by concomitant exposure to both HMW and LMW agents in this setting or synergistic effects between exposures that may be augmenting the inflammatory response. The primary HMW agent in this context is considered to be bacterial endotoxin (9), which has been reported previously to act as an adjuvant for certain allergens (55). An exposure characterization study conducted by our group also confirmed higher endotoxin exposure levels in dental health workers (56). It has also been observed that neutrophils can also be a source of ECP in IgE-mediated responses (57), possibly due to HMW and LMW concomitant exposures. This suggests that the nature of airway inflammation observed in dental health workers may be due to mixed (HMW and LMW) exposures. Furthermore, Suojalehto et al., have recently shown that phenotypic characteristics of acrylate-induced occupational asthma, when compared to other LMW agents, were more similar to occupational asthma caused by HMW agents (15). The role that multiple workplace exposures plays in the manifestation of different inflammatory types of asthma is further elaborated by other investigators (58, 59). The differences between acrylate exposure and other LMW agents common in dental health care settings, has therefore challenged the historical practice of pooling exposures into a single group, on the assumption that they all share similar pathophysiologic mechanisms (15).
In this study, a strong correlation was observed between acute eosinophilic inflammation and asthma-associated clinical indices for work-relatedness including work-related chest symptoms (acute eosinophilia: OR = 4.30, chronic eosinophilia: OR = 4.93), airway obstruction (acute eosinophilia: OR = 5.79, chronic eosinophilia: OR = 8.61), atopic asthma and work-exacerbated asthma (WEA) (acute eosinophilia: OR = 7.75, chronic eosinophilia: OR = 8.12). These clinical indices suggest a dominant allergic component due to its eosinophilic nature caused by agents such as common inhalants, latex and acrylates among others (60). These results are also consistent with finding from the data reported by Hancox et al., that higher eosinophils are associated with lower FEV1 values (22). Interestingly, ECP and MPO did not demonstrate a significant association with any of the asthma-associated clinical indices, except MPO which was positively associated with non-atopic asthma (OR = 3.06, 95% CI: 1.119–8.363, p < 0.05) (data not shown). This lack of association suggests that these two serum biomarkers do not appear to be useful as predictive markers as they may not be present when the effects of these cytokines manifest (7, 19, 33). It would also suggest that cytokine production, is mainly responsible for bronchial hyperresponsiveness (29, 61).
Cytokines commonly associated with eosinophilic inflammation in asthma include IL-3, IL-5, eotaxin, and GM-CSF (6, 23, 24). These cytokines were particularly prominent in the factor named chronic eosinophilic inflammation. Furthermore, subjects with chronic eosinophilia also had an increased odds of having atopic asthma. The other group of cytokines associated with acute eosinophilic inflammation was IL-4 and IL-7. These cytokines have been reported to synergize the induction of IgE and IgG4 production in a T-cell-dependent manner (62). Overexpression of IL-4 in the lungs is associated with a lymphocytic and eosinophilic inflammation, without necessarily causing airway hyperreactivity (23, 63). Hancox et al., have demonstrated that eosinophilia is a risk factor for airflow obstruction even among those without symptoms (22). Furthermore, IL-4 is generally associated with atopic diseases (e.g., allergic rhinitis and atopic dermatitis) (64). This was also observed in the current study, in that the odds of having atopic asthma were increased in both acute and chronic eosinophilic groups relative to the other groups. A previous study has shown that bronchoalveolar lavage samples with elevated CD4+ cells that expressed both IL-4 and IL-17 predicted greater asthma severity (65). Although not investigated in this current study, IL-17 appears to be an independent risk factor for severe asthma and warrants further investigation in occupational populations.
The cytokines associated with the neutrophilic responses in this study included IL-1β, IL-10, TNF-α, IL-12p70 and IL-8, IL-6, which may be potentially caused by agents such as endotoxin and mercury among others (12, 15). The results demonstrated that both acute and chronic neutrophilic groups were more likely to be associated with work-related factors (although less so than the eosinophilic responses) in the absence of HMW agents playing a major role in asthma. Furthermore, subjects with chronic neutrophilic inflammation were more likely to have work-related chest symptoms and non-atopic asthma. Taylor et al., have also shown that neutrophilic asthma is associated with an airway microbiome that is significantly different compared to patients with other inflammatory phenotypes (66). Since endotoxins, a by-product of Gram-negative microorganisms, are likely to result in a neutrophilic response. These findings are therefore consistent with other studies that non-allergic asthma associated with endotoxin exposure, possibly be due to contamination of dental unit waterlines (56, 67), is mediated by an inflammatory response involving several cytokines, including IL-1, IL-6, IL-8, and TNF-α, IL-12 (25, 68).
This study has highlighted the use of biomarkers in subjects with similar clinical characteristics of asthma to ascertain the pathophysiological basis for the symptoms in dental health care settings. While previous studies have used factor analysis to identify distinct asthma phenotypes focusing on clinical measures and quality of life, duration of the asthmatic attack and steroid use (32), this study has further extended the utility of this approach. It has demonstrated the presence of a complex cytokine network and various inflammatory subtypes in asthmatic dental HCWs. This suggests that asthma in this occupational setting is not exclusively due to an eosinophilic inflammatory response. The mixed inflammatory profile may be the result of responses to concomitant exposures to multiple different allergens or may also represent synergistic effects between allergens and agents such as endotoxin. Factors expressed as quantitative traits in the form of cytokine inflammatory patterns to characterize asthma phenotypes can therefore be useful in the exploration of asthma causation in working populations. Furthermore, it is also probable that the level of exposure to multiple agents may also result in the manifestation of a gradation of immune responses. Potential limitations of this study is the cross-sectional nature of the study design and the small study population (especially of the control group), the latter producing some unstable effect estimates that could not be interpreted optimally. Furthermore, the broad definition for asthma may have resulted in disease misclassification, although it could have also contributed to improved sensitivity of identifying subjects with asthma.
In conclusion, this study has demonstrated that dental health workers with similar clinical characteristics of asthma can have heterogeneous subtypes, denoting inflammatory patterns associated with different causative agents, especially when exposed to multiple agents co-exists. Aside from improved diagnosis, this can also contribute toward effective prevention strategies at a population level should the causative agents in this setting be identified.
Data Availability Statement
The datasets presented in this article are not readily available because of the Protection of Personal Information Act and restrictions apply to the availability of these data, which were used for the ethics application. Data are however available from the authors upon reasonable request and with permission of the University of the Witwatersrand Ethics Committee. Requests to access the datasets should be directed to dGFudXNoYXMmI3gwMDA0MDtuaW9oLmFjLnph.
Ethics Statement
The studies involving human participants were reviewed and approved by University of the Witswatersrand, Human Research Ethics Committee, Medical (Ref: R14/49), University of Kwazulu-Natal, University of the Western Cape, University of Pretoria (Ref: 32/2005), University of Michigan Medical School Institutional Review Board (H04-00006341-I). The dental institutions granted permission to conduct the study and the patients/participants provided their written informed consent to participate in this study.
Author Contributions
All authors listed have made a substantial, direct and intellectual contribution to the work, and approved it for publication.
Funding
This project was supported by an NIH Research Grant #2 D43 TW00812-06 through the Fogarty International Center ITREOH Programme awarded to the University of Michigan. Other sources of funding include the Allergy Society of South Africa-Glaxo Smith Kline Research Award and the National Institute for Occupational Health, NHLS.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We would like to acknowledge the special contributions made by the study team and thank them for their support, commitment, and endeavor to complete this work. A special thanks to the research assistants, O. Matuka, L. Majaja, and O. Rammapudu from the NIOH for conducting the fieldwork. We would also like to acknowledge the management of the dental schools for granting permission to conduct this study and the employees and students of the dental schools that participated in this study. We are grateful for the technical contributions from Prof. Andreas Lopata and Dr. Ashraf Kajee.
References
2. Haley KJ, Drazen JM. Inflammation and airway function in asthma: what you see in not necessarily what you get. Am J Respir Crit Care Med. (1998) 157:1–3. doi: 10.1164/ajrccm.157.1.ed-19
3. Eguiluz-Gracia I, Tay TR, Hew M, Escribese MM, Barber D, O'Hehir RE, et al. Recent developments and highlights in biomarkers in allergic diseases and asthma. Allergy. (2018) 73:2290–305. doi: 10.1111/all.13628
4. Agache I, Akdis CA. Endotypes of allergic diseases and asthma: an important step inbuilding blocks for the future of precision medicine. Allergol Int. (2016) 65:243–52. doi: 10.1016/j.alit.2016.04.011
5. Pavord ID, Agusti A. Blood eosinophil count: a biomarker of an important treatable trait in patients with airway disease. Eur Respir J. (2016) 47:1299–303. doi: 10.1183/13993003.00055-2016
6. Brasier AR, Victor S, Boetticher G, Ju H, Lee C, Bleecker ER, et al. Molecular phenotyping of severe asthma using pattern recognition of bronchoalveolar lavage-derived cytokine. J Allergy Clin Immunol. (2008) 121:30–7. doi: 10.1016/j.jaci.2007.10.015
7. de Magalhaes Simoes S, dos Santos MA, da Silva Oliveira M, Fontes ES, Fernezlian S, Garippo AL, et al. Inflammatory cell mapping of the respiratory tract in fatal asthma. Clin Exp Allergy. (2005) 35:602–11. doi: 10.1111/j.1365-2222.2005.02235.x
8. Douwes J, Pearce N, Heederik D. Does environmental endotoxin exposure prevent asthma? Thorax. (2002) 57:86–90. doi: 10.1136/thorax.57.1.86
9. Leigh R, Hargreave FE. Occupational neutrophilic asthma. Can Respir J. (1999) 6:194–6. doi: 10.1155/1999/858549
10. Vandenplas O, Godet J, Hurdubaea L, Rifflart C, Suojalehto H, Wiszniewska M, et al. Are high- and low-molecular-weight sensitizing agents associated with different clinical phenotypes of occupational asthma? Allergy. (2019) 74:261–72. doi: 10.1111/all.13542
11. Anees W, Huggins V, Pavord ID. Occupational asthma due to low molecular weight agents: eosinophilic and non-eosinophilic variants. Thorax. (2002) 57:231–6. doi: 10.1136/thorax.57.3.231
12. Tarlo SM, Quirce S. Impact of identification of clinical phenotypes in occupational asthma. J Allergy Clin Immunol Pract. (2020) 8:3277–82. doi: 10.1016/j.jaip.2020.06.003
13. Portengen L, Preller L, Tielen M, Doekes G, Heederik D. Endotoxin exposure and atopic sensitization in adult pig farmers. J Allergy Clin Immunol. (2005) 115:797–802. doi: 10.1016/j.jaci.2004.11.046
14. Syed M, Chopra R, Sachdev V. Allergic reactions to dental materials-a systematic review. J Clin Diagn Res. (2015) 9:ZE04–9. doi: 10.7860/JCDR/2015/15640.6589
15. Suojalehto H, Suuronen K, Cullinan P, Lindstrom I, Sastre J, Walusiak-Skorupa J, et al. Phenotyping occupational asthma caused by acrylates in a multicenter cohort study. J Allergy Clin Immunol Pract. (2020) 8:971–979.e1. doi: 10.1016/j.jaip.2019.10.017
16. Marraccini P, Brass DM, Hollingsworth JW, Maruoka S, Garantziotis S, Schwartz DA. Bakery flour dust exposure causes non-allergic inflammation and enhances allergic airway inflammation in mice. Clin Exp Allergy. (2008) 38:1526–35. doi: 10.1111/j.1365-2222.2008.03038.x
17. Siroux V, Basagana X, Boudier A, Pin I, Garcia-Aymerich J, Vesin A, et al. Identifying adult asthma phenotypes using a clustering approach. Eur Respir J. (2011) 38:310–7. doi: 10.1183/09031936.00120810
18. Simpson JL, Scott R, Boyle MJ, Gibson PG. Inflammatory subtypes in asthma: assessment and identification using induced sputum. Respirology. (2006) 11:54–61. doi: 10.1111/j.1440-1843.2006.00784.x
19. Watt AP, Schock BC, Ennis M. Neutrophils and eosinophils: clinical implications of their appearance, presence and disappearance in asthma and COPD. Curr Drug Targets Inflamm Allergy. (2005) 4:415–23. doi: 10.2174/1568010054526313
20. Wenzel SE. Asthma: defining of the persistent adult phenotypes. Lancet. (2006) 368:804–13. doi: 10.1016/S0140-6736(06)69290-8
21. Alam R, Busse WW. The eosinophil–quo vadis? J Allergy Clin Immunol. (2004) 113:38–42. doi: 10.1016/j.jaci.2003.10.054
22. Hancox RJ, Pavord ID, Sears MR. Associations between blood eosinophils and decline in lung function among adults with and without asthma. Eur Respir J. (2018) 51:1702536. doi: 10.1183/13993003.02536-2017
23. Diamant Z, Boot JD, Virchow JC. Summing up 100 years of asthma. Respir Med. (2007) 101:378–88. doi: 10.1016/j.rmed.2006.12.004
24. Mahajan S, Mehta AA. Role of cytokines in pathophysiology of asthma. Iran J Pharmacol Ther. (2006) 5:1–14.
25. Haldar P, Pavord ID. Noneosinophilic asthma: a distinct clinical and pathologic phenotype. J Allergy Clin Immunol. (2007) 119:1043–52; quiz 53–4. doi: 10.1016/j.jaci.2007.02.042
26. Mann BS, Chung KF. Blood neutrophil activation markers in severe asthma: lack of inhibition by prednisolone therapy. Respir Res. (2006) 7:59. doi: 10.1186/1465-9921-7-59
27. Mapp CE, Boschetto P, Maestrelli P, Fabbri LM. Occupational asthma. Am J Respir Crit Care Med. (2005) 172:280–305. doi: 10.1164/rccm.200311-1575SO
28. Ordonez CL, Shaughnessy TE, Matthay MA, Fahy JV. Increased neutrophil numbers and IL-8 levels in airway secretions in acute severe asthma: clinical and biologic significance. Am J Respir Crit Care Med. (2000) 161:1185–90. doi: 10.1164/ajrccm.161.4.9812061
29. Wenzel SE, Schwartz LB, Langmack EL, Halliday JL, Trudeau JB, Gibbs RL, et al. Evidence that severe asthma can be divided pathologically into two inflammatory subtypes with distinct physiologic and clinical characteristics. Am J Respir Crit Care Med. (1999) 160:1001–8. doi: 10.1164/ajrccm.160.3.9812110
30. Cicutto LC, Downey GP. Biological markers in diagnosing, monitoring, and treating asthma: a focus on noninvasive measurements. AACN Clin Issues. (2004) 15:97–111. doi: 10.1097/00044067-200401000-00008
31. Doshi U, Salat P, Parikh V. Cytokine modulators in asthma: clinical perspectives. Indian J Pharmacol. (2002) 34:16−25.
32. Lemiere C, Ernst P, Olivenstein R. Airway inflammation assessed by invasive and noninvasive means in severe asthma: eosinophil and noneosinophil phenotypes. J Allergy Clin Immunology. (2006) 118:1033–9. doi: 10.1016/j.jaci.2006.08.003
33. Brusasco V, Crimi E, Pellegrino R. Airway hyperresponsiveness in asthma: not just a matter of airway inflammation. Thorax. (1998) 53:992–8. doi: 10.1136/thx.53.11.992
34. Singh T, Bello B, Jeebhay MF. Risk factors associated with asthma phenotypes in dental healthcare workers. Am J Ind Med. (2013) 56:90–9. doi: 10.1002/ajim.22041
35. Pekkanen J, Sunyer J, Anto JM, Burney P, European community respiratory health s. operational definitions of asthma in studies on its aetiology. Eur Respir J. (2005) 26:28–35. doi: 10.1183/09031936.05.00120104
36. American Thoracic Society (ATS). Standardisation of spirometry - 1994 update. Am J Respir Crit Care Med. (1995) 152:1107–36. doi: 10.1164/ajrccm.152.3.7663792
37. Arm JP, Lee TH. The pathobiology of bronchial asthma. Adv Immunol. (1992) 51:323–82. doi: 10.1016/S0065-2776(08)60491-5
38. van Schalkwyk EM, Schultz C, Joubert JR, White NW. Guidelines for office spirometry in adults, 2004. S Afr Med J. (2004) 94:576–87.
39. Quanjer PH, Tammeling GJ, Cotes JE, Pedersen OF, Peslin R, Yernault JC. Lung volumes and forced ventilatory flows. Eur Respir J. (1993) 16(Suppl):5–40. doi: 10.1183/09041950.005s1693
40. Nriagu J, Robins T, Gary L, Liggans G, Davila R, Supuwood K, et al. Prevalence of asthma and respiratory symptoms in south-central Durban, South Africa. Eur J Epidemiol. (1999) 15:747–55. doi: 10.1023/A:1007653709188
41. Lapperre TS, Snoeck-Stroband JB, Gosman MME, Stolk J, Sont JK, Jansen DF, et al. Dissociation of lung function and airway inflammation in chronic obstructive pulmonary disease. Am J Respir Crit Care Med. (2004) 170:499–504. doi: 10.1164/rccm.200401-112OC
42. Rodrigo G, Rodrigo C. Assessment of the patient with acate asthma in the emergency department: a factor analytic study. Chest. (1993) 104:1325–8. doi: 10.1378/chest.104.5.1325
43. Riekert KA, Eakin M. Factor analysis: a primer for asthma researchers. J Allergy Clin Immunol. (2008) 121:1181–3. doi: 10.1016/j.jaci.2008.03.023
44. Rogerio AP, Cardoso CR, Fontanari C, Souza MA, Afonso-Cardoso SR, Silva EVG, et al. Anti-asthmatic potential of a D-galactose-binding lectin from Synadenium carinatum latex. Glycobiology. (2007) 17:795–804. doi: 10.1093/glycob/cwm053
45. Kasahara DI, Poynter ME, Othman Z, Hemenway D, van der Vliet A. Acrolein inhalation suppresses lipopolysaccharide-induced inflammatory cytokine production but does not affect acute airways neutrophilia. J Immunol. (2008) 181:736–45. doi: 10.4049/jimmunol.181.1.736
46. Lichtman MA, Beutler E, Kipps TJ, Seligsohn U, Kaushansky K, Prchal JT. Williams Hematology. 7th ed. New York, NY: McGraw-Hill Medical (2005).
47. Goldstein D. Statistics for Microarrays: Cluster Analysis. IPAM Group UCLA (2001). Available online at: http://statwww.epfl.ch/davison/teaching/Microarrays/ETHZ/
48. Pillai SG, Tang Y, van den Oord E, Klotsman M, Barnes K, Carlsen K, et al. Factor analysis in the genetics of asthma international network family study identifies five major quantitative asthma phenotypes. Clin Exp Allergy. (2008) 38:421–9. doi: 10.1111/j.1365-2222.2007.02918.x
49. Amelink M, de Nijs SB, de Groot JC, van Tilburg PM, van Spiegel PI, Krouwels FH, et al. Three phenotypes of adult-onset asthma. Allergy. (2013) 68:674–80. doi: 10.1111/all.12136
50. Lemiere C, S NG, Sava F, D'Alpaos V, Huaux F, Vandenplas O. Occupational asthma phenotypes identified by increased fractional exhaled nitric oxide after exposure to causal agents. J Allergy Clin Immunol. (2014) 134:1063–7. doi: 10.1016/j.jaci.2014.08.017
51. Monteseirin J, Bonilla I, Camacho J, Conde J, Sobrino F. Elevated secretion of myeloperoxidase by neutrophils from asthmatic patients: the effect of immunotherapy. J Allergy Clin Immunol. (2001) 107:623–6. doi: 10.1067/mai.2001.113566
52. Jatakanon A, Uasuf C, Maziak W, Lim S, Chung KF, Barnes PJ. Neutrophilic inflammation in severe persistent asthma. Am J Respir Crit Care Med. (1999) 160:1532–9. doi: 10.1164/ajrccm.160.5.9806170
53. Liebers V, Bruning T, Raulf-Heimsoth M. Occupational endotoxin-exposure and possible health effects on humans. Am J Ind Med. (2006) 49:474–91. doi: 10.1002/ajim.20310
54. DunnGalvin A, Hourihane JO, Frewer L, Knibb RC, Oude Elberink JN, Klinge I. Incorporating a gender dimension in food allergy research: a review. Allergy. (2006) 61:1336–43. doi: 10.1111/j.1398-9995.2006.01181.x
55. Liu AH. Endotoxin exposure in allergy and asthma: reconciling a paradox. J Allergy Clin Immunol. (2002) 109:379–92. doi: 10.1067/mai.2002.122157
56. Singh TS, Bello B, Mabe DO, Renton K, Jeebhay MF. Workplace determinants of endotoxin exposure in dental healthcare facilities in South Africa. Ann Occup Hyg. (2010) 54:299–308. doi: 10.1093/annhyg/mep095
57. Monteseirin J, Vega A, Chacon P, Camacho MJ, El Bekay R, Asturias JA, et al. Neutrophils as a novel source of eosinophil cationic protein in IgE-mediated processes. J Immunol. (2007) 179:2634–41. doi: 10.4049/jimmunol.179.4.2634
59. Jeannin P, Delneste Y, Lecoanet-Henchoz S, Gretener D, Bonnefoy JY. Interleukin-7 (IL-7) enhances class switching to IgE and IgG4 in the presence of T cells via IL-9 and sCD23. Blood. (1998) 91:1355–61. doi: 10.1182/blood.V91.4.1355
60. Mabe DO, Singh T, Bello B, Jeebhay MF, Lopata AL, Wadee A. Allergenicity of latex rubber products used in South African dental schools. S Afr Med J. (2009) 99:672–4.
61. Newby C, Heaney LG, Menzies-Gow A, Niven RM, Mansur A, Bucknall C, et al. Statistical cluster analysis of the british thoracic society severe refractory asthma registry: clinical outcomes and phenotype stability. PLoS ONE. (2014) 9:e102987. doi: 10.1371/journal.pone.0102987
62. Renauld JC. New insights into the role of cytokines in asthma. J Clin Pathol. (2001) 54:577–89. doi: 10.1136/jcp.54.8.577
63. Holberg CJ, Halonen M, Solomon S, Graves PE, Baldini M, Erickson RP, et al. Factor analysis of asthma and atopy traits shows 2 major components, one of which is linked to markers on chromosome 5q. J Allergy Clin Immunol. (2001) 108:772–80. doi: 10.1067/mai.2001.119158
64. Moore WC, Hastie AT, Li X, Li H, Busse WW, Jarjour NN, et al. Sputum neutrophil counts are associated with more severe asthma phenotypes using cluster analysis. J Allergy Clin Immunol. (2014) 133:1557–63.e5. doi: 10.1016/j.jaci.2013.10.011
65. Agache I, Ciobanu C, Agache C, Anghel M. Increased serum IL-17 is an independent risk factor for severe asthma. Respir Med. (2010) 104:1131–7. doi: 10.1016/j.rmed.2010.02.018
66. Taylor SL, Leong LEX, Choo JM, Wesselingh S, Yang IA, Upham JW, et al. Inflammatory phenotypes in patients with severe asthma are associated with distinct airway microbiology. J Allergy Clin Immunol. (2018) 141:94–103.e15. doi: 10.1016/j.jaci.2017.03.044
67. Huntington MK, Williams JF, Mackenzie CD. Endotoxin contamination in the dental surgery. J Med Microbiol. (2007) 56(Pt 9):1230–4. doi: 10.1099/jmm.0.47231-0
Keywords: neutrophilic asthma, eosinophilic asthma, cytokines, work-related asthma, phenotypes, endotypes
Citation: Singh T, Bello B and Jeebhay MF (2021) Characterizing Inflammatory Cell Asthma Associated Phenotypes in Dental Health Workers Using Cytokine Profiling. Front. Allergy 2:747591. doi: 10.3389/falgy.2021.747591
Received: 26 July 2021; Accepted: 11 October 2021;
Published: 18 November 2021.
Edited by:
Helena Ribeiro, University of Porto, PortugalReviewed by:
Cristine Secco Rosario, Federal University of Paraná, BrazilArtemii Bogomolov, National Pirogov Memorial Medical University, Ukraine
Copyright © 2021 Singh, Bello and Jeebhay. 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: Tanusha Singh, dGFudXNoYXMmI3gwMDA0MDtuaW9oLmFjLnph
†These authors have contributed equally to this work and share first authorship
‡This author has contributed equally to this work and shares senior authorship