Lukas Storz1,2†
Bettina Schmid1†
Philipp Peter Bosshard1,2
Peter Schmid-Grendelmeier1,2,3
Marie-Charlotte Brüggen1,2,3
Claudia Lang1,2*- 1Department of Dermatology, University Hospital Zurich, Zurich, Switzerland
- 2Faculty of Medicine, University of Zurich, Zurich, Switzerland
- 3Medical Campus Davos, Davos, Switzerland
Background: Atopic dermatitis (AD) is a chronic relapsing inflammatory skin disease in which patients are sensitized towards a plethora of allergens. The hosts fungal microbiota, the mycobiota, that is believed to be altered in patients suffering from AD acts as such an allergen. The correlation context of specific sensitization, changes in mycobiota and its impact on disease severity however remains poorly understood.
Objectives: We aim to enhance the understanding of the specific sensitization towards the mycobiota in AD patients in relation to their fungal skin colonization.
Methods: Sensitization pattern towards the Malassezia spp. and Candida albicans of 16 AD patients and 14 healthy controls (HC) were analyzed with the newly developed multiplex-assay ALEX2® and the established singleplex-assay ImmunoCAP®. We compared these findings with the fungal skin colonization analyzed by DNA sequencing of the internal transcribed spacer region 1 (ITS1).
Results: Sensitization in general and towards Malassezia spp. and C. albicans is increased in AD patients compared to HC with a quantitative difference in severe AD when compared to mild to moderate AD. Further we saw an association between sensitization towards and skin colonization with Candida spp. yet a negative correlation between sensitization towards and skin colonization with Malassezia spp.
Conclusion: We conclude that AD in general and severe AD in particular is associated with increased sensitization towards the hosts own mycobiota. There is positive correlation in Candida spp. skin colonization and negative in Malassezia spp. skin colonization when compared to AD, AD severity as well as to specific sensitization patterns.
Introduction
Atopic dermatitis (AD) is a chronic relapsing inflammatory skin disease. It is characterized by age-dependent typical eczematous lesions, xerosis, pruritus with a chronic disease course and is often associated with allergic rhino conjunctivitis and asthma (1). In most cases onset occurs in early childhood from 3 months to 2 years, yet a substantial portion of these patients go into complete remission before the age of 2 years (2). AD has experienced an epidemic rise, i.e., a three-fold increase over the past 30 years. It now affects 15% to 30% of all children and 2 to10% of the adults (1). Having the highest burden of Disability-adjusted life years (DALY) of all skin diseases (3) AD has an important health effect at population level (4).
Common risk factors of developing AD are mutations in genes encoding for structural proteins such as the loss-of-function mutation within the FLG-gene encoding for filaggrin, epidermal proteases and protease inhibitors (5, 6). The strongest risk factor overall is a positive family history for atopic diseases with active AD in particular (7). Further, environmental properties such as small family size, living in urban settings and regions with low exposure to UV radiation, low humidity and western-diet predispose to AD (4).
With these structural and environmental hazards, the pathophysiology is not yet fully explained. It is widely accepted as a pathophysiological concept that AD evolves due to an interplay of a genetic background, impaired skin barrier as well as an aberrant skin immune response (8, 9). On one hand, the barrier dysfunction allows substances like bacteria or allergens from the environment to penetrate the skin. On the other hand, these substances can stimulate keratinocytes to produce a variety of cytokines and thus trigger a Th2-mediated immune response (1, 10). These environmental particles include the host’s fungal and bacterial skin microbiota (11).
The role of Staphylococcus aureus, a pathogenic bacteria that colonizes human skin with a much higher burden in AD patients than in healthy controls (HC) (12, 13) and acts as a driver of inflammatory skin eruption through immune modulation, is already well studied (14, 15). In contrast studies describing the host’s fungal microbiota (mycobiota) are rare (16, 17). Two of the most abundant fungal microorganisms on human skin are Malassezia spp. and Candida spp. (18, 19).
Malassezia is a lipid-dependent yeast and is by far the most abundant fungus found on Caucasian skin (19, 20). It currently includes 14 species, among which 9 species colonize human and 5 animal skin (21). Within the human colonizing species, there are differences in Th2-dependent immune response found in human keratinocytes and dendritic cells. Dependent on the species the immune pattern shows elevated levels of pro-inflammatory cytokines such as IL-4, IL-13 and IL-17 which compromise the skin’s innate immune response towards microorganisms (17, 22). The microenvironment shows susceptibility for secondary infections on top of amplifying specific AD inflammation patterns (23). While sensitization towards Malassezia spp. has been shown to be an important allergen-specific marker for AD severity in Caucasian adults (24, 25), the impact of the amount of fungal colonization in people suffering from AD remains unclear (25–27).
Another commonly detected yeast on human skin is Candida albicans (C. albicans), which contains several antigenic components capable of stimulating immediate hypersensitivity responses and sensitization towards its allergens and thus contributes to the pathogenesis of AD (11, 28). C. albicans burden seems higher in AD patients compared to HC, yet effects on pathomechanism and disease progression remain unclear (29, 30).
Treatment of AD includes pharmacological and non-pharmacological therapy and should be applied depending on clinical severity an individual’s phenotype in a stepwise manner (31). Standard pharmacological therapy comprises of basic skin care with emollients, topical anti-inflammatory therapy with corticosteroids or calcineurin inhibitors and eventually systemic therapy with short term corticosteroids, ciclosporin, biologicals and Janus kinase inhibitors (31–33). Non-pharmacological treatment includes phototherapy, stress management, psychoeducation and identification and reduction of AD triggers and (31). AD patients with head and neck phenotypes, as well as patients with IgE sensitization towards Malassezia spp. (34) may benefit from antifungal therapy, although additional studies to properly identify beneficial patient characteristics are needed (35). Azole derivates such as ketoconazole and itraconazole are among the most studied antifungal therapies and can be applied topically or systemically (31, 36).
The aim of this study was to examine differences in sensitization towards Malassezia spp. and C. albicans between AD patients and HC and to put it in perspective to the skin mycobiome and to the impact on disease severity.
Methods
Study design
We included 16 AD patients from the Prospective Longitudinal Observational Research in Atopic Dermatitis (ProRaD)-Cohort who were recruited from March 2018 to December 2019 and treated at the University hospital Zurich, Switzerland. Patients being older than 18 years, suffering from atopic dermatitis and who provided written informed consent were eligible to participate. Before enrolment, the clinical diagnosis of AD had to be confirmed by a dermatologist at the Department of Dermatology, University Hospital Zurich, Switzerland, based on the criteria of Hanifin and Rajka (37). Exclusion criteria were refusal to participate on the study or being unable to give consent.
As an age-matched control group, we included 14 healthy controls (HC). Inclusion criteria were age above 18 years, no history of AD or any other skin disease, negative skin prick test (SPT) to common allergens, no chronic medical condition or treatment and neither used antifungals nor antibiotics within 6 months prior to sampling. Exclusion criteria were refusal to participate on the study or being unable to give consent.
This study was approved by the local ethics committee (EK 2016-00301) and conducted according to the Declaration of Helsinki. All study subjects participated voluntarily and gave written informed consent.
AD clinical assessment
All AD patients underwent physical examination by a dermatologist at the Department of Dermatology, University Hospital Zurich, Switzerland. To assess disease severity, the SCORing Atopic Dermatitis (SCORAD) index (38) was used. Classification was made in mild, moderate or severe AD with a SCORAD of <25, 25–50 or >50 (with a maximum score of 103) respectively (39).
Multiple allergen testing with ALEX2®
Blood was withdrawn from all study participants, processed immediately after collection and stored at −80°C until further processing. For ALEX2® protocol, the samples were analyzed and calibrated according to the manufacturer’s instructions (Allergy Explorer-ALEX2® Macro-Array DX Wien, Austria) for the presence of total IgE (tIgE) and specific IgE (sIgE) antibodies using ALEX2®, a colorimetric enzyme multiplex assay. The test contained 183 components and 117 allergen-extracts (Supplementary Table S1). Concentrations of sIgE and tIgE were measured and classified in 5 levels (0–4 corresponding to <0.3 kUA/l, 0.3–1 kUA/l, 1–5 kUA/L, 5–15 kUA/l and >15 kUA/l).
ImmunoCAP®
An ImmunoCAP® (Thermo Fisher Scientific, Massachusetts, United States) assay was used to confirm total sensitization, sensitization to Malassezia spp. antigen and to further investigate specific sensitization towards C. albicans. The samples were analyzed and calibrated according to the manufacturer’s instructions for the presence of tIgE and sIgE using whole allergen extract m227 and m5. Concentrations of sIgE and tIgE were measured and classified in 7 levels (0–6 corresponding to <0.35 kUA/l, 0.35–0.7 kUA/l, 0.7–3.5 kUA/l, 3.5–17.5 kUA/l, 17.5–50 kUA/l, 50–100 kUA/l and >100 kUA/l).
Mycobiome sampling, DNA extraction, amplification, sequencing and taxonomic classification
The fungal skin mycobiota was assessed as in Schmid et al. (40). Shortly, four common eczema skin sites in adults (antecubital crease, dorsal neck, glabella, and vertex) were swabbed with Floqswabs or eSwabs (COPAN, Brescia, Italy) soaked in sterile NaCl (0.9%, Braun, Sempach, Switzerland) by repeatedly rubbing 4–8 cm2 of skin surface on each skin site. As negative controls, swabs without skin contact and water were used.
DNA extraction was performed by adding extraction buffer [1 M Tris-HCl (pH 8), 50 mM EDTA (Thermo Fisher Scientific, Rheinach, Switzerland) and 0.5% Tween 20 (Bio-Rad Laboratories, Cressier, Switzerland)] with proteinase K (Roche Diagnostics, Rotkreuz, Switzerland) to the sample and overnight pre-digestion at 56°C. The following day, the cell walls were mechanically disrupted with 0.5 mm beads (Qiagen, Hilden, Germany) in a Tissuelyser (Qiagen) and further with the Masterpure Yeast DNA Purification kit (Epicentre, LuBioScience GmbH, Zurich, Switzerland) according to the manufacturer’s protocol.
PCR amplification of the internal transcribed spacer region 1 (ITS1) region was performed with the primers 18S-F (5′-GTAAAAGTCGTAACAAGGTTTC-3′) and 5.8S-1R (5′-GTTCAAAGAYTCGATGATTCAC-3′) (19). PCR was executed with 1 to 5 μL DNA and the Kapa Hifi Hotstart polymerase (Roche Diagnostics). Cycling conditions were as follows: pre-incubation at 95°C for 300 s, 33 cycles of 98°C for 20 s, 61°C for 20 s, 72°C for 40 s, and final extension at 72°C for 60 s. Libraries were created following the guidelines of Illumina’s 16S rRNA metagenomics sequencing library preparation, incorporating certain modifications (41). Shortly, the DNA amplicon was evaluated on a 1.5% agarose gel after the first PCR before adding Nextera XT Indices (Illumina, San Diego, CA, United States) by a second PCR with 8 cycles. Agencourt Ampure XP beads (Beckman Coulter, Krefeld, Germany) were used for purification. Subsequently PCR amplicons were quantified with a Qubit 2.0 Fluorometer (Thermo Fisher Scientific), normalized to 4 nM, and pooled. Libraries were finalized with the Miseq Regent kit v3 (Illumina) and 2 × 300 pb paired-end reads were generated with the Miseq sequencing machine (Illumina).
Raw ITS reads were processed with PIPITS (v2.3) (42) with default parameters and the RDP Classifier (v2.1211) (43) against UNITE database (04.02.2020) (44) to generate taxonomic classification. BLAST was used for manual adjustment of the taxonomy table for interesting sequences with >1,000 contigs and >97% identity (45).
Data analysis and statistics
For data analysis and statistics for interpretation of baseline characteristics and sensitization patterns IBM® SPSS ® Statistics for Windows, version 28 (IBM Corp., Armonk, NY, United States) was used. Nonparametric statistical analysis was performed using Kruskal–Wallis test for multiple comparisons among the whole study-population with post-hoc Mann–Whitney-U test to further investigate pairwise subgroup differences. p-values <0.05 were considered as statistically significant. The effect size of Mann–Whitney-U test was calculated using Pearson’s correlation coefficient, where values <0.3 were considered small, between 0.3–0.5 medium and >0.5 large using Cohen’s classification. Correlation analysis was performed using Spearman’s rank-order correlation (R), with two-sided p-values of <0.05 considered as statistically significant.
Data analysis of the ITS sequencing data was performed with the R software (46) by using following packages: phyloseq (47), ape (48), and tidyverse (49). For decontamination, the R package decontam was used by applying the prevalence method with a threshold of 0 (50).
Results
Patient characteristics
A total of 16 patients diagnosed with AD were included in this study, with 9 patients classified as mild to moderate AD and 7 patients as severe AD. 10 (63%) of these 16 patients were males. The median age for all patients was 31 years, ranging from 18 to 65 years. Five (31%) patients had head and neck type AD and median onset of AD was 1 year, ranging from 0 to 38 years. 9 (56%) patients had history of asthma, 14 (88%) patients had history of allergies. Family history for AD was positive in 7 patients (44%). Fourteen patients (88%) used topical steroids during the last 6 months while only 2 (13%) used topical antimycotics (Table 1).
Table 1. Characteristics of study population according to disease-severity including demographic information, disease background and medication.
Gradual increase in tIgE and sensitization towards yeast antigens from HC to severe AD patients
First, we investigated the total amount of tIgE in all 16 AD patients and 14 HC. There was a significant difference between HC and AD patients altogether as well as between HC and mild to moderate and severe AD, respectively. The differences were significant both in ALEX2® (Supplementary Table S2; Mann–Whitney-U test, p < 0.001, Cohen’s effect size r = 0.84 and r = 0.77 respectively) and ImmunoCAP® (Supplementary Table S2; Mann–Whitney-U test, p = 0.002 and p < 0.001 respectively, Cohen’s effect size r = 0.66 and r = 0.77 respectively). tIgE levels did not significantly differ between mild to moderate AD and severe AD. Further, positive and statistically significant correlation with large effect size was found between disease severity and tIgE among both protocols of ImmunoCAP® and ALEX2® (Figure 1).
Figure 1. Correlation of tIgE and disease severity. tIgE using and ALEX2®. tIgE is measured in kUA/l and shown in correlation to the 3 subgroups HC, mild to moderate AD and severe AD. Results of Spearman’s correlations are illustrated in form of R and p-values.
Next, we investigated the specific sensitization patterns towards Malassezia spp., C. albicans and further molds (Alternaria alternata, Aspergillus fumigatus, Cladosporium herbarum, Penicilium chrysogenum) present in the ALEX2® macroarray. There was a significant difference in M. sympodialis Mala s 6 sIgE (ALEX2®) as well as C. albicans m5 sIgE (ImmunoCAP®) between HC and severe AD patients as well as between mild to moderate AD and severe AD (Supplementary Table S2; Mann–Whitney-U test, p = 0.002, p = 0.01 with Cohen’s effect size r = 0.67, r = 0.63, respectively for Mala s 6. p = <0.001, p = 0.63 with Cohen’s effect size r = 0.96, r = 0.63, respectively for m5). Further we saw a significant difference of Malassezia spp. m227 (ImmunoCAP®) between HC and mild to moderate AD as well as between mild to moderate AD to severe AD (Supplementary Table S2; Mann–Whitney-U test, p = 0.01 and p = 0.02 respectively, Cohen’s effect size r = 0.78 and r = 0.57 respectively). There was only one subject of HC that showed a positive yet weak sensitization towards C. albicans m5 with no other positive sensitization towards the other fungal antigens neither in ALEX2® nor in ImmunoCAP® among all HC. When looking at correlation between disease severity and sensitization pattern, we found a positive, statistically significant correlation with large effect size in M. sympodialis Mala s 6 using ALEX2®, Malassezia spp. m227 and C. albicans m5 using ImmunoCAP® respectively with increasing disease severity (Figures 2–4). Between disease severity and the other two M. sympodialis components (Mala s 5 & Mala s 11) in ALEX2®, we saw a positive, statistically significant correlation with medium effect size in Mala s 5 (R = 0.366, p-value = 0.046) and Mala s 11 (R = 0.374, p-value = 0.042) respectively.
Figure 2. Spearman correlations of M. sympodialis ALEX2® and disease severity. sIgE using ALEX2® and displayed as class in correlation to the 3 subgroups HC, mild to moderate AD and severe AD. Darker dots indicate larger number of subjects. Results of Spearman’s correlation are illustrated in form of R and p-values.
Figure 3. Spearman’s correlation of Malassezia spp. ImmunoCAP® and disease severity. sIgE using ImmunoCAP® and displayed as ImmunoCAP® class in correlation to the 3 subgroups HC, mild to moderate AD and severe AD. m227 displays allergens of M. sympodialis, M. globosa and M. restricta. Darker dots indicate larger number of subjects. Results of Spearman’s correlation are illustrated in form of R and p-values.
Figure 4. Spearman’s correlation of C. albicans ImmunoCAP® and disease severity. sIgE m5 using ImmunoCAP® and displayed as ImmunoCAP® class in correlation to the 3 subgroups HC, mild to moderate AD and severe AD. Darker dots indicate larger number of subjects. Results of Spearman’s correlation are illustrated in form of R and p-values.
Concerning further molds in the ALEX2® macroarray, only Penicilium chrysogenum showed a significant difference between HC and severe AD as well as between mild to moderate AD to severe AD (Supplementary Table S2).
Increased sIgE levels against common inhalative and food allergens in AD vs. HC
Next, to tIgE and fungal sIgE levels, the ALEX2® macroarray measures sIgE against a series of allergens such as tree-, grass and weed pollen, dander & epithelia, cockroaches & mites, molds & yeasts as well as food. Being sensitized towards M. sympodialis highly correlates with sensitization towards these common inhalative and food allergens among our study population (Supplementary Table S3). Further, there is a specific sensitization in HC (Supplementary Figure S1) with polysensitization in AD study population (Supplementary Figures S2, S3 respectively).
Comparison between ImmunoCAP® and ALEX2® protocol
Comparing the two methods to detect tIgE and sIgE, we saw a positive, statistically significant and strong correlation between tIgE ImmunoCAP® and ALEX2® (Supplementary Figure S4). Looking at the Malassezia spp. and M. sympodialis extracts and components in ImmunoCAP® and ALEX2®, we only saw medium effect size yet statistically significant correlation between m227 sIgE and Mala s 5 & 6 sIgE, respectively, with weak effect size and no significant correlation between m227 sIgE and Mala s 11 sIgE (Supplementary Figure S5).
Decrease in Malassezia spp. colonization in severe AD
We further explored differences in the skin mycobiota abundance of Malassezia spp. M. sympodialis, M. restricta, M. globosa as displayed in Figure 5.
Figure 5. Colonization pattern of Malassezia spp. for each disease severity subgroup. Whiskers box plot displaying the relative abundance of (A) Malassezia spp., (B) M. sympodialis, (C) M. globosa and (D) M. restricta at different skin site (glabella, vertex, dorsal neck, antecubital crease) for every severity subgroup. Pairwise comparison for differences among the colonization pattern were performed using Mann–Whitney-U test, whereas * and ** indicate a p-value ≤0.05 and ≤0.01, respectively.
We found an overall decrease of Malassezia spp. abundance in severe AD compared to HC and mild to moderate AD. Yet only between mild to moderate AD to severe AD at skin site “glabella” and “antecubital” crease we could show a significant decrease (Supplementary Table S2; Mann–Whitney-U test, p = 0.01 and p = 0.03 respectively, Cohen’s effect size r = 0.65 and r = 0.54 respectively).
Looking at species level, we saw a decrease of M. sympodialis abundance in AD compared to HC with significant differences between HC to AD subgroups at antecubital crease swab site (Supplementary Table S2; Mann–Whitney-U test, p = 0.01 and p = 0.02 respectively, Cohen’s effect size r = 0.60 and r = 0.59 respectively) as well as in HC to mild to moderate AD at dorsal neck swab site (Supplementary Table S2; Mann–Whitney-U test, p = 0.04, Cohen’s effect size r = 0.44). Likewise in M. restricta colonization there was a decrease of abundance in severe AD compared to HC and mild to moderate AD at vertex swab site (Supplementary Table S2; Mann–Whitney-U test, p = 0.04, Cohen’s effect size r = 0.44 and r = 0.52 respectively) and in severe compared to HC at antecubital crease swab site (Supplementary Table S2; Mann–Whitney-U test, p = 0.02, Cohen’s effect size r = 0.52). There was no difference in M. globosa colonization pattern among the 3 different groups.
Increase in Candida spp. colonization in severe AD
When looking at Candida spp. colonization patterns (Figure 6A), we found an increase of relative abundance in severe AD when compared to HC as well as mild to moderate AD. At the vertex and dorsal neck, this difference was significant between severe AD and HC as well as mild to moderate AD (Supplementary Table S2; Mann–Whitney-U test, p = 0.02 and p = 0.05 respectively, Cohen’s effect size r = 0.53 and r = 0.43 respectively). At the glabella there was a significant difference in relative abundance between severe AD and mild to moderate AD (Supplementary Table S2; Mann–Whitney-U test, p = 0.04, Cohen’s effect size r = 0.52), for C. albicans, we saw a very low relative abundance of fungal colonization on species level with values ranging between 0–0.72% and no detectable difference among the 3 subgroups (Figure 6B).
Figure 6. Colonization pattern of Candida spp. for each disease severity subgroup. Whiskers box plot displaying the relative abundance of (A) Candida spp. and (B) C. albicans at different skin site (glabella, vertex, dorsal neck, antecubital crease) for every severity subgroup. Pairwise comparison for differences among the colonization pattern were performed using Mann–Whitney-U test, whereas * and ** indicate a p-value ≤0.05 and ≤0.01, respectively.
Negative correlation between specific sensitization against Malassezia spp. and relative abundance of Malassezia sympodialis
By assessing the colonization patterns, Spearman’s rank correlation analysis showed a statistically significant, negative correlation between M. sympodialis colonization and m227 as well as Mala s 6, respectively, (Supplementary Figures S6A–C). For Malassezia spp., M. restricta and M. globosa colonization, there was a negative yet not significant correlation towards Malassezia spp. sensitization.
Positive correlation between specific sensitization against Candida albicans and relative abundance of Candida spp.
Correlation analysis for Candida spp. and C. albicans colonization showed only a significant positive correlation between Candida spp. at skin site vertex and C. albicans antigen m5 (Supplementary Figure S6D). For other skin sites, correlations were positive yet not significant. For C. albicans colonization, there was no correlation with its antigen m5.
Discussion
Increase in sensitization, decrease in Malassezia spp. colonization
The aim of this study was to examine differences in sensitization towards Malassezia spp. and C. albicans between AD patients HC and put it in perspective to the skin mycobiota and impact on disease severity. This study shows an inverse correlation between sensitization towards and colonization with Malassezia spp. with increasing disease severity in AD patients.
Relative abundance of Malassezia spp. in severe AD patients was significantly decreased (Figure 5A) when compared to mild to moderate AD patients. Further, relative abundance of M. sympodialis and M. restricta decreased significantly with increasing disease severity. Yet sensitization towards M. sympodialis (Mala s 6 in ALEX) and the mixture of components of M. restricta, M. globosa and M. sympodialis (m227 in ImmunoCAP) correlated positively with disease severity (Figure 1). When comparing sensitization and colonization, we found a negative correlation between Malassezia spp. as well as on species level (M. sympodialis, M. restricta, M. globosa) with only significant results in M. sympodialis.
The question arises, how the shown decrease in skin colonization with Malassezia spp. can be put into perspective with increasing sensitization. The overall increase of total and specific sensitization with disease severity (Supplementary Table S2; Figures 1–4) with higher levels in AD patients is already well-known (11, 23, 32, 33) and might be attributable to the fact, that our AD population is on the extrinsic end of the phenotypical AD spectrum and thus is prone to react towards common allergens with the production of IgE (2) (Supplementary Figures S2, S3). On the other hand, among all Malassezia spp., with exception of M. pachydermatis, lipid-dependence for growth is a common feature (8). Pathological features of AD patients next to the increased susceptibility towards allergens and the production of sIgE (9) are the elevation in stratum corneum pH which enhances the activity of degradatory proteases and decreases the activity of lipid synthesis enzymes (6), thus withdrawing the substrate for growth of Malassezia spp. (51, 52). Those findings show the potential use of multiplex-assays such as ALEX2® to detect multiple specific IgE to Malassezia species (amongst other allergens), making it a valuable diagnostic tool in the management of AD. Results could have a major impact on add on antifungal therapies, being another step to precision medicine.
Considering skin colonization, our findings are consistent with the results from Sandström Falk et al. (53) while other studies have described no difference in skin colonization with Malassezia spp. at all (25), or even an increased colonization with Malassezia spp. with increasing disease severity in AD patients (26). When comparing these study populations with ours, the one from Sandström et al. is also Caucasian whilst the other two studies have been conducted on an Asian population. There is evidence, that the skin mycobiota differs among different ethnicities (8, 21, 54) which might be attributed not only to geographical differences but also to different sampling techniques (53). There has even been a description of different specific sensitization towards Malassezia spp. in a sub-Saharan AD population whilst no difference in skin colonization was found (55). Thus, we assume that with increase in inflammation and withdrawal of lipids as a substrate, colonization with Malassezia spp. decreases. Further, we conclude that ethnicity has an impact on the skin mycobiota which makes it difficult to draw conclusions from our study population when compared to other ethnicities. Our findings underline, that ethnic and geographical differences among AD patients could have an impact on a therapeutic level. A distinct group could benefit from antifungal treatment, which could be easily available in wide regions off the world (compared to other systemic therapies) (3).
Increase in Candida spp. colonization
Interestingly along the reduced frequency in skin colonization with Malassezia spp. we saw an increase with Candida spp. in severe AD patients (Figure 6A).
The increase in skin colonization with Candida spp. at the expense of Malassezia spp. is consistent with other studies (30, 56). It is hypothesized that the skin mycobiota is affected by topical and systemic use of corticosteroids and antibiotics. With more frequent use in severe AD it might be attributable to the change in the composition of skin flora eventually resulting in conditions that favor growth of, i.e., Candida spp. (57). When looking at species level, relative abundance of C. albicans was very low with maximum values under 1% and with no difference among HC and AD. Looking at our study population, there is no relevant difference in intake or application of topical or systemic antibiotics or corticosteroids, respectively. Further, only two patients used topical antimycotics, which did not show a relevant effect on colonization. Another hypothesis is, that Candida spp., i.e., being one of the most common pathogenic yeasts of normal skin flora (29) could either profit from or generate an inflammatory environment leading to the above mentioned decrease in lipid synthesis (6) and thus decrease in colonization with Malassezia spp. A recent meta-analysis from He et al. (27) has shown a superiority in treatment of AD using topical antifungals compared to topical glucocorticoids. This suggests that the mycobiota is of central cause in AD progression (17) and already established topical antifungal therapy, nowadays consisting mostly of azole derivates (17, 27, 58), could reduce the amount of allergens causing a type I inflammatory response (35, 59).
Along with the increase in skin colonization the severe AD patients showed significant sensitization towards C. albicans which is in line with the above-mentioned increased specific sensitization in our polysensitized severe AD population and is consistent with findings of other studies (28).
Further investigation concerning causality between inflammatory processes, fungal skin colonization and sensitization towards skin mycobiota could have a major impact on targeted and maybe even preventive management in AD patients. With this assumption, our study could also lay a headstone for clinical investigative trials to recess antifungal treatment approaches in bigger cohorts. This was beyond the scope of our study and with our AD study population counting only 16 patients, and only two (13%) of them using topical antifungals, it would have been difficult to draw conclusion, making sample size one of the major limitations of our study.
Limitations
The sensitization pattern was primarily assessed using the relatively newly developed ALEX2®, a macroarray that allows to obtain multiplex sIgE profile as well as tIgE (60, 61). Due to the removal of C. albicans extract in ALEX2® compared to the ALEX1®, we could not compare the sensitization pattern using this multiplex assay only. Further, ALEX2® could only measure sensitization towards M. sympodialis, while lacking all the other Malassezia spp., especially M. restricta and M. globosa, who are the most abundant on human skin (62, 63), making it one of the major drawbacks of this study. To expand the specific sensitization measurement towards C. albicans, M. restricta and M. globosa, we had to use the ImmunoCAP® assay.
To reduce the influence of fungal colonization of non-measured skin sites on systemic sensitization, we selected lipid rich skin swab sites where we expected Malassezia spp. to be present in a high number. Nevertheless, the possibility that the sIgE towards Malassezia spp. is elevated due to colonization of other skin sites and colonization in total is underestimated with increasing disease severity cannot be excluded.
Another already above-mentioned limitation is the small sample size which makes our study population vulnerable to confounders such as sex, topical and systemic medication, and skin care regimens. Further, we saw differences especially concerning colonization patterns among HC and AD patients, with many of them not being significant. This might as well be due to the very small sample size.
Strengths
One of the strengths of this study is the direct comparison of specific sensitization towards the allergen and the allergen load, i.e., sIgE of Malassezia spp. and C. albicans with and their abundance on human skin. This is the first study using this design on a Caucasian study population while comparing AD with HC with only one study from Zhang et al. (25) using a similar approach on AD patients only. With assessing tIgE and sIgE using ImmunoCAP®, we were able to use the best studied sIgE single-and multiplex assay (64). It shows substantial to good accordance to ALEX® (61), and we were able to cross-validate the total and specific sensitization which correlated strongly (Supplementary Figure S4).
A further strength lies in the use of a PCR-approach to measure fungal skin colonization. Malassezia spp. are very difficult to cultivate and even if growth occurs, some species grow much slower than others and thus would be overwhelmed, leading a wrong representation of Malassezia spp. which has been described as cultivation bias (62).
Conclusion/outlook
With this study we gathered further evidence for increased specific sensitization as well as change in composition of skin mycobiota in severe AD patients. Further investigation with larger sample size to put these findings into a clinical context need to be done.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found at: https://www.ebi.ac.uk/ena/browser/view/PRJEB44167.
Ethics statement
The studies involving humans were approved by Kantonale Ethikkommission Zürich; EK 2016-00301. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.
Author contributions
LS: Writing – original draft, Writing – review & editing, Conceptualization, Investigation, Methodology, Project administration, Visualization, Formal analysis, Software. BS: Conceptualization, Investigation, Methodology, Project administration, Writing – original draft, Writing – review & editing, Formal analysis, Software. PB: Conceptualization, Supervision, Writing – review & editing. PS-G: Writing – review & editing, Conceptualization, Supervision, Funding acquisition, Project administration, Resources, Validation. M-CB: Funding acquisition, Project administration, Resources, Supervision, Validation, Writing – review & editing, Conceptualization. CL: Funding acquisition, Project administration, Resources, Supervision, Validation, Writing – review & editing, Conceptualization, Methodology.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was supported by funds of the University Zurich, Faculty of Medicine, the Christine Kühne foundation (CkCare) and the Pierre Fabre Eczema Foundation.
Acknowledgments
Special thanks go to Muriel Rentschler for her advice in planning and support in conducting the fungal sensitization pattern analysis.
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/fmed.2024.1353784/full#supplementary-material
References
2. Bieber, T, D’Erme, AM, Akdis, CA, Traidl-Hoffmann, C, Lauener, R, Schäppi, G, et al. Clinical phenotypes and endophenotypes of atopic dermatitis: where are we, and where should we go? J Allergy Clin Immunol. (2017) 139:S58–64. doi: 10.1016/j.jaci.2017.01.008
3. Laughter, MR, Maymone, MBC, Mashayekhi, S, Arents, BWM, Karimkhani, C, Langan, SM, et al. The global burden of atopic dermatitis: lessons from the global burden of disease study 1990–2017*. Br J Dermatol. (2021) 184:304–9. doi: 10.1111/bjd.19580
4. Weidinger, S, and Novak, N. Atopic dermatitis. Lancet. (2016) 387:1109–22. doi: 10.1016/S0140-6736(15)00149-X
5. Irvine, AD, McLean, WHI, and Leung, DYM. Filaggrin mutations associated with skin and allergic diseases. N Engl J Med. (2011) 365:1315–27. doi: 10.1056/NEJMra1011040
6. Cork, MJ, Danby, SG, Vasilopoulos, Y, Hadgraft, J, Lane, ME, Moustafa, M, et al. Epidermal barrier dysfunction in atopic dermatitis. J Invest Dermatol. (2009) 129:1892–908. doi: 10.1038/jid.2009.133
7. Apfelbacher, CJ, Diepgen, TL, and Schmitt, J. Determinants of eczema: population-based cross-sectional study in Germany. Allergy. (2011) 66:206–13. doi: 10.1111/j.1398-9995.2010.02464.x
8. Gaitanis, G, Magiatis, P, Hantschke, M, Bassukas, ID, and Velegraki, A. The Malassezia genus in skin and systemic diseases. Clin Microbiol Rev. (2012) 25:106–41. doi: 10.1128/CMR.00021-11
9. Ong, PY, Ohtake, T, Brandt, C, Strickland, I, Boguniewicz, M, Ganz, T, et al. Endogenous antimicrobial peptides and skin infections in atopic dermatitis. N Engl J Med. (2002) 347:1151–60. doi: 10.1056/NEJMoa021481
10. Ishibashi, Y, Sugawara, K, Sugita, T, and Nishikawa, A. Secretion of thymic stromal lymphopoietin from human keratinocytes in response to Malassezia yeasts. J Dermatol Sci. (2011) 62:134–8. doi: 10.1016/j.jdermsci.2011.02.012
11. Sonesson, A, Bartosik, J, Christiansen, J, Roscher, I, Nilsson, F, Schmidtchen, A, et al. Sensitization to skin-associated microorganisms in adult patients with atopic dermatitis is of importance for disease severity. Acta Derm Venereol. (2013) 93:340–5. doi: 10.2340/00015555-1465
12. Paller, AS, Kong, HH, Seed, P, Naik, S, Scharschmidt, TC, Gallo, RL, et al. The microbiome in patients with atopic dermatitis. J Allergy Clin Immunol. (2019) 143:26–35. doi: 10.1016/j.jaci.2018.11.015
13. Edslev, SM, Agner, T, and Andersen, PS. Skin microbiome in atopic dermatitis. Acta Derm Venereol. (2020) 100:adv00164. doi: 10.2340/00015555-3514
14. Travers, JB, Norris, DA, and Leung, DYM. The keratinocyte as a target for staphylococcal bacterial toxins. J Investig Dermatology Symp Proc. (2001) 6:225–30. doi: 10.1046/j.0022-202x.2001.00045.x
15. Inoue, Y, Aihara, M, Kirino, M, Harada, I, Komori-Yamaguchi, J, Yamaguchi, Y, et al. Interleukin-18 is elevated in the horny layer in patients with atopic dermatitis and is associated with Staphylococcus aureus colonization. Br J Dermatol. (2011) 164:560–7. doi: 10.1111/j.1365-2133.2010.10145.x
16. Zhang, E, Tanaka, T, Tajima, M, Tsuboi, R, Nishikawa, A, and Sugita, T. Characterization of the skin fungal microbiota in patients with atopic dermatitis and in healthy subjects. Microbiol Immunol. (2011) 55:625–32. doi: 10.1111/j.1348-0421.2011.00364.x
17. Glatz, M, Bosshard, P, Hoetzenecker, W, and Schmid-Grendelmeier, P. The role of Malassezia spp. in atopic dermatitis. J Clin Med. (2015) 4:1217–28. doi: 10.3390/jcm4061217
18. Thammahong, A, Kiatsurayanon, C, Edwards, SW, Rerknimitr, P, and Chiewchengchol, D. The clinical significance of fungi in atopic dermatitis. Int J Dermatol. (2020) 59:926–35. doi: 10.1111/ijd.14941
19. Findley, K, Oh, J, Yang, J, Conlan, S, Deming, C, Meyer, JA, et al. Human skin fungal diversity. Nature. (2013) 498:367–70. doi: 10.1038/nature12171
20. Cabañes, FJ. Malassezia yeasts: how many species infect humans and animals? PLoS Pathog. (2014) 10:1–4. doi: 10.1371/journal.ppat.1003892
21. Harada, K, Saito, M, Sugita, T, and Tsuboi, R. Malassezia species and their associated skin diseases. J Dermatol. (2015) 42:250–7. doi: 10.1111/1346-8138.12700
22. Sparber, F, De Gregorio, C, Steckholzer, S, Ferreira, FM, Dolowschiak, T, Ruchti, F, et al. The skin commensal yeast Malassezia triggers a type 17 response that coordinates anti-fungal immunity and exacerbates skin inflammation. Cell Host Microbe. (2019) 25:389–403.e6. doi: 10.1016/j.chom.2019.02.002
23. Ishibashi, Y, Sugita, T, and Nishikawa, A. Cytokine secretion profile of human keratinocytes exposed to Malassezia yeasts. FEMS Immunol Med Mocrobiol. (2006) 48:400–9. doi: 10.1111/j.1574-695X.2006.00163.x
24. Glatz, M, Buchner, M, von Bartenwerffer, W, Schmid-Grendelmeier, PS, Worm, M, Hedderich, J, et al. Malassezia spp. specific immunoglobulin E level is a marker for severity of atopic dermatitis in adults. Acta Derm Venereol. (2015) 95:191–6. doi: 10.2340/00015555-1864
25. Zhang, E, Tanaka, T, Tajima, M, Tsuboi, R, Kato, H, Nishikawa, A, et al. Anti-Malassezia-specific IgE antibodies production in Japanese patients with head and neck atopic dermatitis: relationship between the level of specific IgE antibody and the colonization frequency of cutaneous Malassezia species and clinical severity. J Allergy. (2011) 2011:645670. doi: 10.1155/2011/645670
26. Kaga, M, Sugita, T, Nishikawa, A, Wada, Y, Hiruma, M, and Ikeda, S. Molecular analysis of the cutaneous Malassezia microbiota from the skin of patients with atopic dermatitis of different severities. Mycoses. (2011) 54:e24–8. doi: 10.1111/j.1439-0507.2009.01821.x
27. He, D, Han, Y, Wu, H, Liu, M, Li, M, and Lu, M. Treatment of atopic dermatitis using topical antifungal drugs: a meta-analysis. Dermatol Ther. (2022) 35:1–9. doi: 10.1111/dth.15930
28. Chang, FY, Lee, JH, Yang, YH, Yu, HH, Wang, LC, Lin, YT, et al. Analysis of the serum levels of fungi-specific immunoglobulin E in patients with allergic diseases. Int Arch Allergy Immunol. (2010) 154:49–56. doi: 10.1159/000319208
29. Faergemann, J. Atopic dermatitis and fungi. Clin Microbiol Rev. (2002) 15:545–63. doi: 10.1128/CMR.15.4.545-563.2002
30. Ghaffari, J, Taheri Sarvtin, M, Hedayati, MT, Hajheydari, Z, Yazdani, J, and Shokohi, T. Erratum: evaluation of Candida colonization and specific humoral responses against Candida albicans in patients with atopic dermatitis. Biomed Res Int. (2015) 2015:142453. doi: 10.1155/2015/142453
31. Werfel, T, Heratizadeh, A, Aberer, W, Augustin, M, Biedermann, T, Bauer, A, et al. S3 guideline atopic dermatitis: part 1—general aspects, topical and non-drug therapies, special patient groups. J Dtsch Dermatol Ges. (2023) 22:137–53. doi: 10.1111/ddg.15230
32. Werfel, T, Heratizadeh, A, Aberer, W, Augustin, M, Biedermann, T, Bauer, A, et al. S3 guideline atopic dermatitis: part 2—systemic treatment. J Dtsch Dermatol Ges. 1–14 (2023). doi: 10.1111/ddg.15229
33. Patruno, C, Potestio, L, and Napolitano, M. Clinical phenotypes of adult atopic dermatitis and related therapies. Curr Opin Allergy Clin Immunol. (2022) 22:242–9. doi: 10.1097/ACI.0000000000000837
34. Brodská, P, Panzner, P, Pizinger, K, and Schmid-Grendelmeier, P. IgE-mediated sensitization to Malassezia in atopic dermatitis: more common in male patients and in head and neck type. Dermatitis. (2014) 25:120–6. doi: 10.1097/DER.0000000000000040
35. Kaffenberger, BH, Mathis, J, and Zirwas, MJ. A retrospective descriptive study of oral azole antifungal agents in patients with patch test-negative head and neck predominant atopic dermatitis. J Am Acad Dermatol. (2014) 71:480–3. doi: 10.1016/j.jaad.2014.04.045
36. Svejgaard, E, Larsen, PØ, Deleuran, M, Ternowitz, T, Roed-Petersen, J, and Nilsson, J. Treatment of head and neck dermatitis comparing itraconazole 200 mg and 400 mg daily for 1 week with placebo. J Eur Acad Dermatology Venereol. (2004) 18:445–9. doi: 10.1111/j.1468-3083.2004.00963.x
38. Severity scoring of atopic dermatitis: the SCORAD index. Consensus Report of the European Task Force on Atopic Dermatitis. Dermatology. (1993) 186:23–31. doi: 10.1159/000247298
39. Oranje, AP, Glazenburg, EJ, Wolkerstorfer, A, and de Waard-van der Spek, FB. Practical issues on interpretation of scoring atopic dermatitis: the SCORAD index, objective SCORAD and the three-item severity score. Br J Dermatol. (2007) 157:645–8. doi: 10.1111/j.1365-2133.2007.08112.x
40. Schmid, B, Künstner, A, Fähnrich, A, Bersuch, E, Schmid-Grendelmeier, P, Busch, H, et al. Dysbiosis of skin microbiota with increased fungal diversity is associated with severity of disease in atopic dermatitis. J Eur Acad Dermatology Venereol. (2022) 36:1811–9. doi: 10.1111/jdv.18347
41. Illumina. 16S metagenomic sequencing library. (2013). p. 1–28. Available at: http://support.illumina.com/content/dam/illumina-support/documents/documentation/chemistry_documentation/16s/16s-metagenomic-library-prep-guide-15044223-b.pdf
42. Gweon, HS, Oliver, A, Taylor, J, Booth, T, Gibbs, M, Read, DS, et al. PIPITS: an automated pipeline for analyses of fungal internal transcribed spacer sequences from the Illumina sequencing platform. Methods Ecol Evol. (2015) 6:973–80. doi: 10.1111/2041-210X.12399
43. Wang, Q, Garrity, GM, Tiedje, JM, and Cole, JR. Naïve Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl Environ Microbiol. (2007) 73:5261–7. doi: 10.1128/AEM.00062-07
44. Nilsson, RH, Larsson, KH, Taylor, AFS, Bengtsson-Palme, J, Jeppesen, TS, Schigel, D, et al. The UNITE database for molecular identification of fungi: handling dark taxa and parallel taxonomic classifications. Nucleic Acids Res. (2019) 47:D259–64. doi: 10.1093/nar/gky1022
45. Ciardo, DE, Schär, G, Böttger, EC, Altwegg, M, and Bosshard, PP. Internal transcribed spacer sequencing versus biochemical profiling for identification of medically important yeasts. J Clin Microbiol. (2006) 44:77–84. doi: 10.1128/JCM.44.1.77-84.2006
46. R Foundation for Computing. R Core Team (2013). R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. Available at: http://www.R-project.org/. (Accessed August 27, 2023).
47. McMurdie, PJ, and Holmes, S. Phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data. PLoS One. (2013) 8:e61217. doi: 10.1371/journal.pone.0061217
48. Paradis, E, and Schliep, K. Ape 5.0: an environment for modern phylogenetics and evolutionary analyses in R. Bioinformatics. (2019) 35:526–8. doi: 10.1093/bioinformatics/bty633
49. Wickham, H, Averick, M, Bryan, J, Chang, W, McGowan, L, François, R, et al. Welcome to the Tidyverse. J Open Source Softw. (2019) 4:1686. doi: 10.21105/joss.01686
50. Davis, NM, DiM, P, Holmes, SP, Relman, DA, and Callahan, BJ. Simple statistical identification and removal of contaminant sequences in marker-gene and metagenomics data. Microbiome. (2018) 6:226. doi: 10.1186/s40168-018-0605-2
51. Nakabayashi, A, Sei, Y, and Guillot, J. Identification of Malassezia species isolated from patients with seborrhoeic dermatitis, atopic dermatitis, pityriasis versicolor and normal subjects. Med Mycol. (2000) 38:337–41. doi: 10.1080/mmy.38.5.337.341
52. Schafer, L, and Kragballe, K. Abnormalities in epidermal lipid metabolism in patients with atopic dermatitis. J Invest Dermatol. (1991) 96:10–5. doi: 10.1111/1523-1747.ep12514648
53. Sandström Falk, MH, Linder, MT, Johansson, C, Bartosik, J, Bäck, O, Särnhult, T, et al. The prevalence of Malassezia yeasts in patients with atopic dermatitis, seborrhoeic dermatitis and healthy controls. Acta Derm Venereol. (2005) 85:17–23. doi: 10.1080/00015550410022276
54. Torrelo, A. Atopic dermatitis in different skin types. What is to know? J Eur Acad Dermatol Venereol. (2014) 28:2–4. doi: 10.1111/jdv.12480
55. Lang, CCV, Masenga, J, Semango, G, Kaderbhai, H, Li, N, Tan, G, et al. Evidence for different immune signatures and sensitization patterns in sub-Saharan African vs. Central European atopic dermatitis patients. J Eur Acad Dermatology Venereol. (2021) 35:e140–2. doi: 10.1111/jdv.16871
56. Arzumanyan, VG, Magarshak, OO, and Semenov, BF. Yeast fungi in patients with allergic diseases: species variety and sensitivity to antifungal drugs. Bull Exp Biol Med. (2000) 129:601–4. doi: 10.1007/BF02434889
57. Oliveira, ADT, Sodré, CS, Ferreira, DC, Abad, ED, Saintive, S, Ribeiro, M, et al. Oral aspects identified in atopic dermatitis patients: a literature review. Open Dent J. (2018) 12:424–34. doi: 10.2174/1874210601812010424
58. Taudorf, EH, Jemec, GBE, Hay, RJ, and Saunte, DML. Cutaneous candidiasis—an evidence-based review of topical and systemic treatments to inform clinical practice. J Eur Acad Dermatology Venereol. (2019) 33:1863–73. doi: 10.1111/jdv.15782
59. Saunte, DML, Gaitanis, G, and Hay, RJ. Malassezia-associated skin diseases, the use of diagnostics and treatment. Front Cell Infect Microbiol. (2020) 10:112. doi: 10.3389/fcimb.2020.00112
60. Heffler, E, Puggioni, F, Peveri, S, Montagni, M, Canonica, GW, and Melioli, G. Extended IgE profile based on an allergen macroarray: a novel tool for precision medicine in allergy diagnosis. World Allergy Organ J. (2018) 11:7. doi: 10.1186/s40413-018-0186-3
61. Buzzulini, F, Da Re, M, Scala, E, Martelli, P, Conte, M, Brusca, I, et al. Evaluation of a new multiplex assay for allergy diagnosis. Clin Chim Acta. (2019) 493:73–8. doi: 10.1016/j.cca.2019.02.025
62. Gemmer, CM, DeAngelis, YM, Theelen, B, Boekhout, T, and Dawson, TL. Fast, noninvasive method for molecular detection and differentiation of Malassezia yeast species on human skin and application of the method to dandruff microbiology. J Clin Microbiol. (2002) 40:3350–7. doi: 10.1128/JCM.40.9.3350-3357.2002
63. Sugita, T, Tajima, M, Tsubuku, H, Tsuboi, R, and Nishikawa, A. Quantitative analysis of cutaneous Malassezia in atopic dermatitis patients using real-time PCR. Microbiol Immunol. (2006) 50:549–52. doi: 10.1111/j.1348-0421.2006.tb03825.x
Keywords: atopic dermatitis, Malassezia spp., Candida spp., sensitization patterns, skin microbiota
Citation: Storz L, Schmid B, Bosshard PP, Schmid-Grendelmeier P, Brüggen M-C and Lang C (2024) Decreased skin colonization with Malassezia spp. and increased skin colonization with Candida spp. in patients with severe atopic dermatitis. Front. Med. 11:1353784. doi: 10.3389/fmed.2024.1353784
Edited by:
Wei Li, Zhejiang University, ChinaReviewed by:
Luca Potestio, University of Naples Federico II, ItalyVijaykumar Patra, Research and Innovation, L’Oreal, France
Copyright © 2024 Storz, Schmid, Bosshard, Schmid-Grendelmeier, Brüggen and Lang. 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: Claudia Lang, claudia.lang@usz.ch
†These authors have contributed equally to this work and share first authorship