Spatiotemporal Variation of the Indoor Mycobiome in Daycare Centers

doi:10.21203/rs.3.rs-147547/v1

Download PDF

Research

Spatiotemporal Variation of the Indoor Mycobiome in Daycare Centers

https://doi.org/10.21203/rs.3.rs-147547/v1

This work is licensed under a CC BY 4.0 License

You are reading this latest preprint version

Background: Children spend considerable time in daycare centers in parts of the world, and are exposed to the indoor micro- and mycobiomes of these facilities. The level of exposure to microorganisms varies within and between buildings, depending on occupancy, climate and season. In order to evaluate indoor air quality, and the effect of usage and seasonality, we investigate the spatiotemporal variation in the indoor mycobiomes of two daycare centers. We collected dust samples from different rooms throughout a year and analyzed their mycobiomes using DNA metabarcoding.

Results: The fungal community composition in rooms with limited occupancy (auxiliary rooms) was similar to the outdoor samples, and clearly different from the rooms with higher occupancy (main rooms). The main rooms had higher abundance of Ascomycota, while the auxiliary rooms contained comparably more Basidiomycota. We observed a strong seasonal pattern in the mycobiome composition, mainly structured by the outdoor climate. Most markedly, basidiomycetes of the orders Agaricales and Polyporales, mainly reflecting typical outdoor fungi, were more abundant during summer and fall. In contrast, ascomycetes of the orders Saccharomycetales and Capnodiales were dominant during winter and spring.

Conclusions: Our findings provide clear evidences that the indoor mycobiome in daycare centers are structured by occupancy as well as outdoor seasonality. We conclude that the temporal variability should be accounted for in indoor mycobiome studies and in the evaluation of indoor air quality of buildings.

General Microbiology

Fungi

seasonal variation

DNA metabarcoding

dust mycobiome

built environment

microbial ecology

fungal community

Humans spend significant amount of time indoors, in private homes, but also in workplaces, schools, daycare centers and hospitals. We share these indoor environments with a variety of microorganisms, including microscopic fungi that may affect our health in different ways. In moist conditions, fungi can propagate and act as sources of indoor pollutants leading to poor indoor air quality. This has been associated with adverse health effects, such as allergies, asthma and other respiratory symptoms [1, 2]. The indoor microorganisms originate from both indoor and outdoor sources and are potentially structured by numerous factors, including building features, building usage, the number and type of occupants, and not least, our behavior [3, 4]. The bacterial indoor microbiome is known to be highly affected by the occupants and their activities, and often directly related to the human body [5, 6]. However, indoor fungi, which can be referred to as the indoor mycobiome, are known to be highly influenced by the outdoor air and climate [5, 7, 8]. Previous studies at large geographical scale in the US and Norway, have demonstrated that the composition of the indoor mycobiomes significantly correlates with variables of the outdoor environment (i.e. climate, soil and vegetation) [9, 10]. The most important indoor sources of fungi include occupants, pets, food, waste, plants, plumbing systems, mold damages, heating, ventilation and air conditioning [11]. Different rooms in buildings may have different mycobiome composition due to different occupancy and exposure to outdoor air [12, 13]. For example, central rooms with higher activity, like the kitchen and living room, promote dust resuspension in the air that facilitate dispersal of fungi from occupants, their activities and outdoor sources. Similarly, floor dust of high activity rooms contains higher levels of skin-associated yeasts of the genera Rhodotorula, Candida, Cryptococcus, Malassezia, and Trichosporon [14].

The indoor mycobiomes may not only differ in space, but also in time. Previous culture-based studies have been reviewed by Nevalainen et al. [15], where they found a general pattern of seasonal variation with lower concentrations of airborne fungi in winter than in summer. This review included studies from different climatic regions in countries like Australia [16], Denmark [17], and Taiwan [18]. DNA-based studies have also reported a clear seasonal variation of fungal richness, diversity and community composition in indoor environments, in both dust and air samples [7, 19]. By analyzing dust samples from a university housing facility in California, Adams et al. [7] reported higher fungal richness in winter than in summer. Likewise, Weikl et al. [19] showed a drop of the fungal diversity in summer, based on floor dust samples from 286 houses in Munich. This latter observation was explained by the high prevalence of a few dominant taxa during summer [19]. Hence, observed temporal trends in indoor mycobiomes are not uniform.

In boreal and temperate climatic regions, the fungal spore diversity and composition in outdoor air are expected to vary significantly more throughout the year because of clear seasons. For example, Karlsson et al. [20] reported lowest richness of fungi and bacteria for air samples collected during winter in two climatic zones from Sweden. It can be expected that this variation influences the indoor mycobiome, due to an influx of spores into buildings. Many fungi, especially basidiomycetes, produces fruit bodies during the fall leading to a relatively higher spore abundance during this period [21]. Plant pathogens, dominated by ascomycetes, may have a wider temporal distribution since many spread asexual spores during the entire plant growth season [22]. Indoor fungi originating from indoor sources, here growing on available organic materials, can be expected to have a year-round growth and sporulation connected to human activity.

A particularly interesting environment to study the spatiotemporal variation of the mycobiome is daycare centers, where children, at least in parts of the world, spend a considerable amount of time. For example, in Norway, 92.2% of children between 1–5 years old are in daycares. This particular built environment is characterized by a high occupancy with high levels of activity, and higher fungal concentrations have been detected here compared to private homes [23]. Exploring the indoor mycobiome and revealing the factors driving this spatiotemporal variation are important not only to understand the ecological context of indoor fungi, but also to recognize the effect that some fungal species may have on children’s health. To what degree the mycobiome associated with daycares affect the children’s health is still unknown.

The overarching aim of this study is to reveal the indoor mycobiomes spatiotemporal dynamics in daycare centers in order to improve evaluations of air quality in indoor air. We expect rooms with different occupancy to differ in mycobiome composition (Hypothesis 1; H1), with frequently accessed rooms being dominated by indoor fungi derived from the occupants and their activities. Given that part of the indoor mycobiome originates from outdoor sources, we hypothesize that indoor mycobiomes fluctuate with seasons (H2). In seasons with optimal fungal growth conditions outdoors, as in summer and fall, we expect that a higher proportion of the indoor mycobiome is derived from outdoor sources, with Basidiomycota dominating during the fall season (H3). In contrast, we expect that a higher proportion of the mycobiome has an indoor origin with increased amount of time spend inside during winter and spring (H4). To test these hypotheses, we collected indoor dust and outdoor air samples from two daycare centers bi-weekly during a year and performed DNA metabarcoding of the rDNA ITS2 region. Two daycare centers located in Oslo, Norway, were selected for the study. We collected dust swab samples every second week from different rooms and stores in the daycare centers (Fig. 1), as well as outdoor air samples every week. Fungi present in the samples were surveyed through DNA metabarcoding analyses of the rDNA ITS2 region.

Mycobiome composition

A variation partitioning analysis of the indoor dust mycobiome (Fig. 2a) revealed that 37% of the compositional variation could be ascribed to assessed factors, including outdoor climate, time (i.e. the biweekly sampling point), space and occupant characteristics. Most of the explained variation was accounted for by the combined effects of occupants and room type (19%). These two factors are correlated, as the activity of both staff and children are considerably lower in the basement and loft (hereafter called the auxiliary rooms) than in the remaining indoor rooms (hereafter referred to as the main rooms). Nine percent of the variation was accounted for by time alone, likely reflecting other unmeasured environmental factors changing with time, while 8% was accounted for by the combined effect of time and climate, which again are tightly coupled.

These findings were also displayed in multivariate analyses (Fig. 2b-d), with a distinct difference in mycobiome composition between samples from the auxiliary rooms and the main rooms, the latter used more frequently by the staff and children (Fig. 2b). The dust mycobiomes from the auxiliary rooms were similar to the mycobiome obtained from the outside environment (Fig. S1). The main rooms from both daycare centers overlapped in mycobiome composition, although the compositional variation (beta-diversity) across samples was higher in some rooms, like the kitchen and staff room (Fig. 2c).

In addition to the compositional differences related to room type and occupancy, we observed a clear temporal variation in the indoor mycobiome composition (Fig. 2d). Although there was some overlap, the winter and spring samples were more similar in fungal community composition, as were the samples from summer and fall. The temporal trend in mycobiome composition correlated with the yearly variation in temperature and moisture, as could be seen from the vectors fitted in Fig. 2d. Room, time and occupants were the most important factors structuring the fungal community composition, whereas the impact of which daycare the samples were from, was low (Table 1).

Table 1

Adonis test results showing the influence of the variables on the composition of the dust mycobiome from the complete indoor dataset (auxiliary and main rooms), and from the main rooms of two daycare centers sampled in Oslo, Norway throughout a full year.
Variable	Auxiliary and main rooms		Main rooms
	R squared	p value	R squared	p value
Room	0.218	0.001	0.092	0.001
Month	0.196	0.001	0.233	0.001
Children age	0.133	0.001	0.056	0.001
Season	0.115	0.001	0.124	0.001
Nr of adults	0.114	0.001	0.032	0.001
Nr of children	0.097	0.001	0.027	0.001
Temperature	0.071	0.001	0.082	0.001
Daycare	0.07	0.001	0.033	0.001
Dust coverage	0.039	0.001	0.008	0.029
Moisture	0.02	0.001	0.022	0.001

Taxonomic variation

Overall, Ascomycota was more prominent in the main room while Basidiomycota was far more abundant in samples from the auxiliary rooms (Fig. 3a). Mucoromycota was more prevalent in the main rooms, as well as other fungal phyla (mainly Chytridiomycota). Ascomycete yeasts affiliated to Saccharomycetales were more abundant in the main rooms, while basidiomycetes from the orders Agaricales and Polyporales were dominating the samples from auxiliary rooms. We observed a clear temporal trend in the composition of fungal taxonomic groups (displayed at order level in Fig. 3a). Most markedly, the proportion of basidiomycete sequences from the orders Agaricales, Polyporales and Hymenochaetales, mainly reflecting outdoor fungi, were higher during the growth season (May-November) than in winter, a pattern far more prominent in the auxiliary rooms than in the main rooms. The Ascomycota was proportionally more abundant in colder periods. However, the Saccharomycetales, likely derived from indoor sources, was prevalent year-round in the main room. The order Eurotiales, including fungal genera with allergenic potentials, such as Penicillium and Aspergillus, was relatively more prevalent in the main rooms in the colder season. Similar trends were observed in the OTU ordination plot (Fig. 3b), where the dominant Ascomycota OTUs are associated with the main rooms, while the Basidiomycota OTUs with the auxiliary rooms. Further, the main rooms are dominated by OTUs of Saccharomycetales, Mucorales, Malasseziales and Filobasidiales.

Indicator species analyses, assessing which fungal OTUs followed a significant temporal trend on a monthly basis, revealed that numerous OTUs in the already mentioned orders of Agaricales, Polyporales and Hymenochaetales increased considerably during their expected fruiting season, independently of space (i.e., room) (Fig. 3c).

Richness and evenness trends

The main and auxiliary rooms had comparable fungal richness that largely followed a similar temporal trend, with higher richness in the summer and fall (June-November) (Fig. 4a). The richness followed the variation in annual temperature. In winter, the richness deviated more from the moisture gradients. The evenness followed a similar trend as the richness (Fig. S4).

A random forest model, which was used to assess the contribution of each factor in the observed richness patterns (Fig. 4b), revealed that month and season (both enclosing various environmental factors), as well as, temperature and moisture, accounted for much of the variation in both the main and the auxiliary rooms datasets. In addition, the factor rooms was highly important in the main rooms dataset, where presence of children and adults also contributes to the richness.

In this study, we observed that the indoor mycobiomes of two daycare centers were strongly structured by room type and occupancy (hypothesis H1) and, further, that the mycobiomes changed systematically throughout the seasons (hypothesis H2). No marked difference in mycobiome composition was observed between the two studied daycare centers, indicating a common pattern of indoor mycobiomes in daycare centers from the same local geographic region.

Spatial distribution

We observed a clear separation in mycobiome composition of the main rooms and auxiliary rooms, which likely can be explained by the number of people accessing and using the rooms. These results strongly support our hypothesis (H1) and further suggest that occupancy is an important factor shaping the indoor mycobiome, in addition to the outdoor air. The outdoor samples (air sampling) were collected as point samples in one day, while the auxiliary room samples represent a collection of dust accumulated within two weeks. These different sampling methods may influence the recovered mycobiomes. Nevertheless, the mycobiome composition of the auxiliary rooms and the outdoor mycobiome were highly similar, which supports that occupancy strongly affects the indoor mycobiome. Previous studies of indoor environments suggest that the indoor mycobiomes are highly affected by outdoor air [7, 9, 10, 24, 25]. Most of these studies have not accounted for indoors environments with different levels of activity. However, in a recent study in private homes in Norway, we demonstrated that the number of inhabitants affected the indoor mycobiome composition [10].

The highest fungal richness was found in the main rooms. This may be explained by that the indoor air of the main rooms includes outdoor taxa, in addition to the more specific indoor fungi derived from the occupants and their activities. Higher fungal richness in indoor environments than in outdoor air has also been found in private houses and schools [5, 10, 14]. It should be noted, though, that richness analyses based on DNA-metabarcoding is vulnerable to various biases. For example, if some dominant species are present, they may mask the remaining richness during the PCR process, since their DNA templates may outcompete the rarer species during PCR amplification. However, the evenness follows largely the same trend for both types of rooms and are therefore probably not causing significant biases for the richness analyses.

OTUs of the phylum Basidiomycota were overrepresented in the auxiliary rooms, whereas there were relatively more OTUs of Ascomycota in the main rooms. Likewise, previous studies have demonstrated a predominance of Ascomycota in indoor samples, while Basidiomycota prevails in outdoor samples [26, 27]. As the auxiliary rooms were more similar to the outdoor air, we expected basidiomycetes to be more prevalent in these rooms, especially the mushroom-forming Agaricales and Polyporales. In the main rooms, the high abundance of ascomycetes can be explained by their high tolerance towards environmental stressors, such as high temperature and low water availability, typical conditions in indoor environments. The orders Saccharomycetales and Capnodiales were the most abundant ascomycetes. Saccharomycetales are yeasts including the well-known genera Saccharomyces, associated with foods, and the potential human pathogen Candida. Capnodiales, with the widespread genus Cladosporium, includes both plant and human pathogens [28]. In addition, the basidiomycete orders Malasseziales and Filobasidiales, together with the order Mucorales were abundant in the main rooms. These orders includes yeasts and molds, and were also more abundant in indoor mycobiomes than outdoor air in our previous study of private homes [10].

Seasonality

We observed a clear seasonal pattern in the indoor mycobiomes, supporting our hypothesis H2. Collection month was best able to explain the variation in fungal richness in all rooms. This seasonal pattern is further supported by the evenness and richness analyses of time series, which follows the shift of temperature and moisture throughout the year. Our observed patterns mirror those found in seasonal studies on outdoor mycobiomes. For example, in northern Sweden, the outdoor fungal communities shifted throughout the season [20]. Since the outdoor fungal community has a strong impact on indoor mycobiomes, it is expected that seasonal changes in the outdoor environment also affect which fungi occur indoor. During the spring, summer and fall, with temperatures above zero, fungal activity and sporulation are clearly linked to the level of precipitation (i.e. rainfall). However, during winter, the precipitation manifests largely as snow, which has less effect on the fungal communities at sub-zero temperatures. It has been suggested that during the winter, when the ground is frozen and covered by snow, the impact of the outdoor fungal community on the indoor mycobiome is limited [29]. This can explain the drop of richness during the winter observed in our study.

Although our results demonstrate that dust sampling can be used to reveal the seasonal variation in the indoor mycobiome, there are methodological constraints that should be taken into consideration when analyzing samples with relatively small amount of DNA. All samples were treated equally in the laboratory, nevertheless, some steps could represent sources of heterogeneity in the dataset e.g. variability in DNA extraction efficiency among organisms, primer bias to different fungal taxonomic group, PCR bias and sequencing errors, which might affect the fungal community. Some of the rooms with limited occupancy, such as the auxiliary rooms and the staff room, had a considerably lower amount of dust and potentially lower amount of DNA. Likewise, the outdoor air samples collected during winter contained lower amount of DNA, most likely due to considerably lower number of fungal spores in the air compared to other seasons.

We observed higher abundances of basidiomycetes during summer and fall in all rooms, with a predominance of Agaricomycetes, confirming our hypothesis H3. Agaricomycetes cover the mushroom-forming species that typically disperse spores during the summer and fall in high latitude ecosystems. In addition, more indicator species, showing a distinct temporal pattern, were found in these two seasons, in particular from Agaricales and Polyporales. Thus, high outdoor spore production of basidiomycetes during the summer and fall affects the indoor mycobiome. A high outdoor aerial abundance of basidiomycetes during summer and fall was also observed in northern Sweden [20]. However, our findings are rather opposite to what has been found in seasonal studies in Munich (Germany) and California [7, 19]. Weikl and colleagues explained their observed decline in diversity during the summer in houses in Munich with a few highly abundant OTUs, and not necessarily of lower diversity [19]. Further, in California the summers are warm and dry, and the mushroom-forming species of the Agaricomycetes often fruit during late fall and winter. Thus, all studies may show the same pattern of higher richness of indoor mycobiome during the outdoor sporulation period of basidiomycetes.

The ascomycetes, especially prevalent in indoor conditions, were proportionally more abundant indoor during winter and spring compared to summer and fall. This confirms hypothesis H4. At sub-zero temperatures during winter and spring, fungal growth and sporulation outdoor is reduced in the study area (Oslo, Norway), which will limit the input of basidiomycetes to the indoor environment. Instead, ascomycetes of indoor origin will be more prevalent during this time-period. Similar findings were reported in a seasonal study of indoor mycobiomes of four office complexes, where ascomycete molds and basidiomycete yeasts were more common in the spring and winter [25]. In contrast, in another study monitoring airborne fungi in four daycares centers over 12 months through culturing, viable counts of major indoor fungi were significantly lower in the winter [30]. Overall, they found the ascomycetes Cladosporium, Penicillium, Alternaria, and Aspergillus to be the most dominating genera. These genera, considered to be some of the most allergenic fungi normally present indoors and outdoors, have also been reported as abundant in other studies [31, 32]. De Ana et al. investigated the seasonal distribution of these species, and found that the highest presence of Aspergillus, Cladosporium and Penicillium in the indoor environment was registered in fall, whereas Alternaria was more frequent in the summer [32]. In our study, the order Eurotiales, including Penicillium and Aspergillus, was relatively more prevalent in the main rooms in the colder season. In addition, the genera Saccharomyces, Cladosporium and Didymella, often encountered in indoor environments in other studies [31], were also especially prevalent in the winter. Numerous indoor ascomycetes are known to cause allergies and disease in humans, and it is a concern if these species have a higher prevalence during the winter when the children spend more time inside. In addition, in a previous study of school environments [14], they showed that occupancy contributed more to the allergenic fungal populations in indoor air than outdoor fungi. Understanding this spatiotemporal variation of the indoor mycobiome is important as the time spent inside during the different seasons varies, and will reflect how the children are affected by these fungal species.

In conclusion, our study demonstrates clear differences in the dust mycobiome composition in daycare centers between rooms with different occupancy. The more human activity, the more the indoor mycobiome differs from the outdoor mycobiome composition. To our knowledge, this is the first study that monitors the same rooms and buildings continuously over a full year using a DNA metabarcoding approach. Thus, our results demonstrate how the mycobiome composition follows a strong seasonal trend, mirroring outdoor weather conditions. Knowledge about the seasonal trends will have important implications for monitoring and evaluation of indoor air quality.

Sampling

Dust samples from the daycare centers were collected with floq swabs (Copan Italia spa, Brescia, Italy) and adhesive tapes (Mycotape 2, Mycoteam AS, Oslo, Norway) from 30 x 40 cm² glass plates located 1–2 m above floor level. The swab collected dust from an area of 30 x 30 cm², whereas the tapes sampled dust from 3.8 x 7.5 cm² from the remaining area to calculate the percentage of dust coverage. These samples were collected once for every sampling date. The plates were placed in different rooms and stores in the daycares (Fig. 1). Five rooms were sampled in daycare A, and four rooms in daycare B. The plates were sterilized with 85% ethanol after each harvesting, every second week throughout a year. In addition, outdoor air samples were collected every week throughout a year by processing approximately 1800 L air through a 25 mm cassette with a 0.8 µm pore diameter mixed cellulose ester filter (Zefon international, Ocala, FL, USA) by using an air pump. The 294 swab and filter samples were stored at -80 ºC until DNA extraction, whereas the adhesive tapes were directly scanned for dust coverage using Epson Perfection V850 Pro (Seiko Epson Corporation, Nagano, Japan). The percentage of dust coverage was calculated with the Olympus Stream v 1.9 software.

DNA extraction and fungal metabarcoding

DNA from swabs and filter samples were extracted using the E.Z.N.A Soil DNA kit (Omega Bio-tek, Norcross, GA, USA). The swabs and filters were placed in disruptor tubes using sterilized scissor or forceps, respectively, and 800 µL SLX-Mlus Buffer was added. The samples were homogenized for 2 x 1 min at 30 Hz using TissueLyser (Qiagen, Hilden, Germany) and stored at -20 °C until further processing. The samples were thawed at 70 °C, following an incubation of 10 min and homogenized twice for 1 min at 30 Hz using a TissueLyser. The samples were cooled on ice before 600 µL chloroform was added. Then, the samples were vortexed and centrifuged at 13 000 rpm for 5 min at RT. The aqueous phase was transferred to a new 1.5 mL tube and an equal volume of XP1 Buffer was added before vortexing. The samples were then added to the HiBind DNA Mini Column and further processed by following the manufacturer´s guidelines. The extracted genomic DNA was eluted in 50 µL of elution buffer.

The ITS2 region was targeted by using the forward primer ITS4 (5′-xCTCCGCTTATTGATATG; White et al., 1990) and the reverse primer gITS7 (5′-xGTGARTCATCGARTCTTTG; Ihrmark et al., 2012). The sample barcodes x ranged from 6–9 base pairs. The PCR reaction contained 2 µl DNA template and 23 µl master mix; 14.6 µl Milli-Q water, 2.5 µl 10x Gold buffer, 0.2 µl dNTP’s (25 nM), 1.5 µl reverse and forward primers (10 µM), 2.5 µl MgCl2 (50 mM), 1.0 µl BSA (20 mg/ml) and 0.2 µl AmpliTaq Gold polymerase (5 U/µl, Applied Biosystems, Thermo Fisher Scientific). For samples with low DNA concentration (weak gel bands), 5 µl DNA template and 20 µl master mix were used. The DNA was amplified by initial denaturation at 95 °C for 5 min, followed by 32 cycles of denaturation at 95 °C for 30 s, annealing at 55 °C for 30 s, and elongation at 72 °C for 1 min. A final elongation step was included at 72 °C for 10 min. PCR products were normalized by using the SequalPrep Normalization Plate Kit (Invitrogen, Thermo Fisher Scientific, Waltham, MA, USA) and eluted in 20 µL elution buffer.

The resulting 345 PCR products, including technical replicates, negative PCR controls and mock community (1 ng/µL equimolar DNA concentration from an artificial mix of Mycena belliarum, Pycnoporellus fulgens, Serpula similis and Pseudoinonotus dryadeus), were processed in a total of four metabarcoding libraries. The technical replicates included DNA from 12 dust samples and were included in each library. The 96 uniquely barcoded PCR products within each library were pooled, and the pools were concentrated and purified using Agencourt AMPure XP magnetic beads (Beckman Coulter, CA, USA). The quality of the purified pools was controlled using Qubit (Invitrogen, Thermo Fisher Scientific, Waltham, MA, USA). The four libraries were barcoded with Illumina adapters, spiked with PhiX and sequenced in two Illumina MiSeq (Illumina, San Diego, CA, USA) runs with 2 x 250 bp paired-end reads at Fasteris SA (Plan-les-Ouates, Switzerland). The resulting metabarcoding dataset comprised 41 126 514 sequences.

Bioinformatics

The raw forward and reverse sequences, were demultiplexed independently on a sample basis using CUTADAPT v 2.7 [33], allowing no mismatches between barcode tags and sequence primer, and sequences shorter than 100 bp where discarded. DADA2 [34] was used to filter low quality sequences, with a maximum expected error of 2.5 and to correct read errors based on a machine learning model built from the sequence data. We then merged the error-corrected forward and reverse sequences using a minimum overlap of 5 bp. Chimeras were filtered out using the bimera algorithm, with default parameters implemented in DADA2 v.12. The resulting 28 346 amplicon sequence variants were further clustered into operational taxonomic units (OTUs) using VSEARCH [35] at 97% similarity. LULU [36] was used with default settings to correct for potential OTU over-splitting. Taxonomy was assigned using BLAST [37] to the final OTU table using the UNITE database [38]. All the negative PCR controls and most of the negative DNA controls were automatically removed during the bioinformatics because the number of sequences was too low. The OTUs of the remaining controls were inspected to assess any contamination issues. The final dataset (excluding controls and replicates) contained 6 800 OTUs accounting for 18 694 392 reads from 292 retained samples. The number of reads per sample varied from 470 (from outdoor air during the winter) to 257 599 with a mean value of 65 365. The number of OTUs per sample varied from 3 to 1 259.

Environmental variables

Climatic variables were retrieved from The Norwegian Climate and Service Center (https://klimaservicesenter.no/, accessed March 11th, 2020), recorded by the Meteorological station at Blindern, Oslo, Norway. The daycare centers are located within a 500 m radius to the meteorological station. The climatic variables included: mean air temperature, mean dew point temperature, max air temperature, min air temperature, mean cloud area fraction, mean water vapor partial pressure, mean surface air pressure, mean wind speed, max relative humidity, mean relative humidity, min relative humidity, humidity mixing ratio, specific humidity, snow coverage, surface snow thickness, amount of precipitation and duration of sunshine. The variables were downloaded for each week throughout the year, and averages for every two weeks prior to sampling were calculated and used for seasonal analyses. These variables were studied with principal component analyses (Fig. S2). The results indicated that the first and second dimensions explained a total of 75.6% of the variance. The first dimension was clearly correlated with variables associated with temperature while the second dimension was associated with variables related to humidity and moisture. The coordinates of dimension 1 and 2 of the PCA analyses were designated as temperature and moisture, respectively, and used as surrogate for all the above-mentioned climatic variables in downstream analyses. Season was also included as a variable, with related data averaged accordingly. The following months were grouped in four different seasons: winter from December to February, spring from March to May, summer from June to August and fall from September to November. In addition, the number of children, age of children and number of working adults (staff) having access to each daycare center and room between two sampling dates were recorded and included as variables. Continuous variables were scaled using the scale function in R.

Statistical analyses

All statistical analyses were performed in R version 3.6.2 [39] through RStudio (version 1.3.959) unless stated otherwise. We first confirmed the similarity of the technical replicates by nonmetric multidimensional scaling (NMDS) using the metaMDS function from the vegan package version 2.4-2 [40], and visualized by ggplot2 [41] (Fig. S3). Then, the complete dataset was rarefied to 1 649 sequences sample-wise, using the function rrarefy (vegan). This led to three samples being discarded for downstream statistical analyses, because of shallow sequencing depth in these samples. We then transformed the abundance of OTU per sample table (OTU table) into Hellinger abundance, using the decostand function (vegan). The community structure was analyzed using NMDS as described above. A stable solution, for NMDS, was searched with a maximum number of 200 random starts and iterations with the convergence criteria set to stress and/or scale factor of the gradient below 1 x 10e-7, using a Bray-Curtis dissimilarity distance. The community structure was visualized using ggplot2 [41] with the axes transformed into half-change units.

The results showed a clear distinction between outdoor samples and indoor samples, with the exception of auxiliary rooms (the indoor samples belonging to rooms with very low frequency of occupancy), which showed very similar patterns to the outdoor samples. Considering that the outdoor samples were collected in a different way and time-frame and that outdoor air seasonality was not the main focus of our hypotheses, we decided to focus on the indoor space and refrained for downstream analyses regarding the outdoor samples. Since the indoor samples showed a clear segregation between auxiliary rooms and the main rooms, we decided to analyze the indoor data in two separate sets; auxiliary and main rooms together, and only the main rooms. For both datasets we rarefied all the samples to the sample with the lowest number reads in the respective dataset, 2 657 sequences for the auxiliary and main rooms dataset and 3 381 sequences in the main rooms dataset. We used the same procedure described above to analyze and visualize community composition. The function envfit (vegan) was used to regress the environmental variables onto the Bray-Curtis dissimilarity matrix. Significance of the regression was assessed using 999 permutations. The variables with significant effect were overlaid as vectors in the ordination (NMDS) graphic with arrows pointing in the increasing direction. In addition, we used the function adonis2 (vegan) with 999 permutations to perform a permutational multivariate analysis of variance to assess the contribution of each environmental variable in explaining variability in the community structure. Additionally, we performed variation partitioning analysis using varpart (vegan) to assess the interaction and total variability explained by the following groups of variables: climate (temperature and moisture), time (month and season), space (daycare and room), and occupants (number of adults and children, and age of children).

The taxonomic compositional summary was achieved by summing all the rarefied reads, at the order level, within a sample and averaged across the time period. Richness, Shannon-Weaver and evenness indices were determined using the functions, specnumber, diversity (vegan) and , respectively. Richness and Shannon-Weaver were strongly correlated; we therefore retained richness as a representative of alpha diversity. To estimate the effects of temperature, moisture, season, month, room, dust coverage, number of children and adults and children age on richness, we conducted linear models followed by analyses of variance as implemented by lm and anova functions in R [39]. Random forest models with permutations, as implemented in the R package rfPermute [42], with all predictor variables randomly sampled at each tree node, 500 trees and 999 permutations were applied to determine the significance and importance of each variable. In all models the squared-root of richness was used to normalize the response variable.

We further identified indicator OTUs on a monthly basis for indicator species analyses using the function multipatt in the R package indicspecies [43]. We then retained only OTUs with a p < 0.05 and present in at least three samples per month. The results were summarized by number of OTUs per order per month. The full lists of indicator OTUs of the auxiliary rooms and the main rooms are provided in Table S1 and Table S2, respectively.

Ethics approval and consent to participate

Personal data of the occupants in the day care centers is not collected, the addresses of the day care centers are omitted from the metadata.

Consent for publication

Not applicable.

Availability of data and materials

All sequencing data, metadata, and scripts used in the bioinformatic analyses and statistics are deposited in Dryad Digital Repository (http://dx.doi.org/) and will be made available upon publication or earlier by request.

Competing interests

The authors declare no competing interests.

Funding

The research was financially supported by the University of Oslo and the Norwegian Asthma and Allergy Association (NAAF). PMMS was funded by the European Union’s Horizon 2020 research and innovation programme (Marie Skłodowska-Curie Individual Fellowship; grant agreement MycoIndoor No 741332).

Authors´ contributions

ELFE, IBE, HK and IS planned and initiated the study, ELFE and IS collected the dust samples, ELFE did all lab and bioinformatic analyses, LM did the statistical analyses. All authors drafted the manuscript, discussed the results and agreed about the final version.

Acknowledgements

We would like to thank the daycare centers for allowing us to perform sampling and for providing metadata of the usage of the daycare centers. Mycoteam AS contributed to the sampling and provided sampling equipment.

Bornehag CG, Blomquist G, Gyntelberg F, Jarvholm B, Malmberg P, Nordvall L, Nielsen A, Pershagen G, Sundell J: Dampness in buildings and health. Nordic interdisciplinary review of the scientific evidence on associations between exposure to "dampness" in buildings and health effects (NORDDAMP). Indoor air 2001, 11:72-86.
Mendell MJ, Mirer AG, Cheung K, Tong M, Douwes J: Respiratory and allergic health effects of dampness, mold, and dampness-related agents: A review of the Epidemiologic Evidence. Environ Health Pers 2011, 119:748-756.
Adams RI, Bhangar S, Dannemiller KC, Eisen JA, Fierer N, Gilbert JA, Green JL, Marr LC, Miller SL, Siegel JA et al: Ten questions concerning the microbiomes of buildings. Building Environ 2016, 109:224-234.
Horve PF, Lloyd S, Mhuireach GA, Dietz L, Fretz M, MacCrone G, Van Den Wymelenberg K, Ishaq SL: Building upon current knowledge and techniques of indoor microbiology to construct the next era of theory into microorganisms, health, and the built environment. J Expo Sci Environ Epidemiol 2020, 30:219-235.
Barberán A, Dunn RR, Reich BJ, Pacifici K, Laber EB, Menninger HL, Morton JM, Henley JB, Leff JW, Miller SL et al: The ecology of microscopic life in household dust. Proc Biol Sci 2015, 282.
Lax S, Smith DP, Hampton-Marcell J, Owens SM, Handley KM, Scott NM, Gibbons SM, Larsen P, Shogan BD, Weiss S et al: Longitudinal analysis of microbial interaction between humans and the indoor environment. Science 2014, 345:1048-1052.
Adams RI, Miletto M, Taylor JW, Bruns TD: Dispersal in microbes: fungi in indoor air are dominated by outdoor air and show dispersal limitation at short distances. ISME 2013, 7:1262-1273.
Amend AS, Seifert KA, Samson R, Bruns TD: Indoor fungal composition is geographically patterned and more diverse in temperate zones than in the tropics. PNAS 2010, 107:13748-13753.
Barberan A, Ladau J, Leff JW, Pollard KS, Menninger HL, Dunn RR, Fierer N: Continental-scale distributions of dust-associated bacteria and fungi. PNAS 2015, 112:5756-5761.
Martin-Sanchez PM, Estensmo ELF, Morgado LN, Maurice S, Engh IB, Skrede I, Kauserud H: Analyzing indoor mycobiomes through a large-scale citizen science study of houses from Norway. Research Square 2020.
Prussin AJ, Marr LC: Sources of airborne microorganisms in the built environment. Microbiome 2015, 3:78.
Adams RI, Miletto M, Lindow SE, Taylor JW, Bruns TD: Airborne bacterial communities in residences: Similarities and differences with fungi. PloS one 2014, 9:e91283.
Hegarty B, Haverinen-Shaughnessy U, Shaughnessy RJ, Peccia J: Spatial gradients of fungal abundance and ecology throughout a damp building. Environ Sci Technol Lett 2019, 6:329-333.
Yamamoto N, Hospodsky D, Dannemiller KC, Nazaroff WW, Peccia J: Indoor emissions as a primary source of airborne allergenic fungal particles in classrooms. Environ Sci Technol 2015, 49:5098-5106.
Nevalainen A, Täubel M, Hyvärinen A: Indoor fungi: companions and contaminants. Indoor air 2015, 25:125-156.
Garrett MH, Rayment PR, Hooper MA, Abramson MJ, Hooper BM: Indoor airborne fungal spores, house dampness and associations with environmental factors and respiratory health in children. Clin Exp Allergy 1998, 28:459-467.
Frankel M, Bekö G, Timm M, Gustavsen S, Hansen EW, Madsen AM: Seasonal variations of indoor microbial exposures and their relation to temperature, relative humidity, and air exchange rate. Appl Environ Microbiol 2012, 78:8289-8297.
Su HJ, Wu PC, Chen HL, Lee FC, Lin LL: Exposure assessment of indoor allergens, endotoxin, and airborne fungi for homes in southern Taiwan. Environ Res 2001, 85:135-144.
Weikl F, Tischer C, Probst AJ, Heinrich J, Markevych I, Jochner S, Pritsch K: Fungal and bacterial communities in indoor dust follow different environmental determinants. PloS one 2016, 11:e0154131.
Karlsson E, Johansson A-M, Ahlinder J, Lundkvist MJ, Singh NJ, Brodin T, Forsman M, Stenberg P: Airborne microbial biodiversity and seasonality in Northern and Southern Sweden. PeerJ 2020, 8:e8424.
Peay KG, Bruns TD: Spore dispersal of basidiomycete fungi at the landscape scale is driven by stochastic and deterministic processes and generates variability in plant–fungal interactions. New Phytol 2014, 204:180-191.
Chen W, Hambleton S, Seifert KA, Carisse O, Diarra MS, Peters RD, Lowe C, Chapados JT, Lévesque CA: Assessing performance of spore samplers in monitoring aeromycobiota and fungal plant pathogen diversity in Canada. Appl Environ Microbiol 2018, 84:e02601-02617.
Madureira J, Paciência I, Rufo JC, Pereira C, Teixeira JP, de Oliveira Fernandes E: Assessment and determinants of airborne bacterial and fungal concentrations in different indoor environments: Homes, child day-care centres, primary schools and elderly care centres. Atmos Environ 2015, 109:139-146.
Frankel M, Bekö G, Timm M, Gustavsen S, Hansen EW, Madsen AM: Seasonal variations of indoor microbial exposures and their relation to temperature, relative humidity, and air exchange rate. Appl Environ Microbiol 2012, 78:8289.
Pitkäranta M, Meklin T, Hyvärinen A, Paulin L, Auvinen P, Nevalainen A, Rintala H: Analysis of fungal flora in indoor dust by ribosomal DNA sequence analysis, quantitative PCR, and culture. Appl Environ Microbiol 2008, 74:233-244.
Hanson B, Zhou Y, Bautista EJ, Urch B, Speck M, Silverman F, Muilenberg M, Phipatanakul W, Weinstock G, Sodergren E et al: Characterization of the bacterial and fungal microbiome in indoor dust and outdoor air samples: a pilot study. Environ Sci Process Impacts 2016, 18:713-724.
Hanson B, Zhou Y, Bautista EJ, Urch B, Speck M, Silverman F, Muilenberg M, Phipatanakul W, Weinstock G, Sodergren E et al: Characterization of the bacterial and fungal microbiome in indoor dust and outdoor air samples: a pilot study. Environ Sci Process Impacts 2016, 18:713-724.
Crous PW, Schoch CL, Hyde KD, Wood AR, Gueidan C, de Hoog GS, Groenewald JZ: Phylogenetic lineages in the Capnodiales. Stud Mycol 2009, 64:17-47S17.
Reponen T, Nevalainen A, Jantunen M, Pellikka M, Kalliokoski P: Normal range criteria for indoor air bacteria and fungal spores in a subarctic climate. Indoor air 1992, 2:26-31.
Aydogdu H, Asan A: Airborne fungi in child day care centers in Edirne City, Turkey. Environ Monit Assess 2008, 147:423-444.
Fairs A, Wardlaw A, Thompson J, Pashley CH: Guidelines on ambient intramural airborne fungal spores. J Investig Allergol Clin Immunol 2010, 20:490-498.
de Ana SG, Torres-Rodríguez JM, Ramírez EA, García SM, Belmonte-Soler J: Seasonal distribution of Alternaria, Aspergillus, Cladosporium and Penicillium species isolated in homes of fungal allergic patients. J Investig Allergol Clin Immunol 2006, 16:357-363.
Martin M: Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet 2011, 17:3.
Callahan BJ, McMurdie PJ, Rosen MJ, Han AW, Johnson AJ, Holmes SP: DADA2: High-resolution sample inference from Illumina amplicon data. Nat Methods 2016, 13:581-583.
Rognes T, Flouri T, Nichols B, Quince C, Mahé F: VSEARCH: a versatile open source tool for metagenomics. PeerJ Preprints 2016, 4:e2409v2401.
Frøslev TG, Kjøller R, Bruun HH, Ejrnæs R, Brunbjerg AK, Pietroni C, Hansen AJ: Algorithm for post-clustering curation of DNA amplicon data yields reliable biodiversity estimates. Nat Comm 2017, 8:1188.
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol 1990, 215:403-410.
Koljalg U, Larsson KH, Abarenkov K, Nilsson RH, Alexander IJ, Eberhardt U, Erland S, Hoiland K, Kjoller R, Larsson E et al: UNITE: A database providing web-based methods for the molecular identification of ectomycorrhizal fungi. New Phytol 2005, 166:1063-1068.
R Core Team: R: A language and environment for statistical computing. R Foundation for Statistical Computing. Austria: Vienna. In. https://www.R-project.org/; 2018.
Oksanen J, Blanchet FG, Friendly M, Kindt R, Legendre P, McGlinn D, Minchin PR, O’hara R, Simpson GL, Solymos P: Vegan: Community ecology package. R package version 2.4-2. Vienna: R Foundation for Statistical Computing; 2016.
Wickham H: ggplot2: elegant graphics for data analysis: Springer; 2016.
Archer E: rfPermute: Estimate permutation p-values for random forest importance metrics. In., R package version 2.1.81 edn. https://CRAN.R-project.org/package=rfPermute; 2020.
De Cáceres M, Sol D, Lapiedra O, Legendre P: A framework for estimating niche metrics using the resemblance between qualitative resources. Oikos 2011, 120:1341-1350.

Download PDF

Editorial decision: Minor revision
10 Mar, 2021
Review #3 received at journal
03 Mar, 2021
Review #1 received at journal
15 Feb, 2021
Review #2 received at journal
15 Feb, 2021
Reviewer #3 agreed at journal
06 Feb, 2021
Reviewer #2 agreed at journal
04 Feb, 2021
Reviewers invited by journal
04 Feb, 2021
Reviewer #1 agreed at journal
04 Feb, 2021
Editor assigned by journal
12 Jan, 2021
Submission checks completed at journal
12 Jan, 2021
Editor invited by journal
12 Jan, 2021
First submitted to journal
11 Jan, 2021

You are reading this latest preprint version

Spatiotemporal Variation of the Indoor Mycobiome in Daycare Centers

Status:

Version 1

Abstract

Figures

Background

Results

Mycobiome composition

Taxonomic variation

Richness and evenness trends

Discussion

Spatial distribution

Seasonality

Conclusions

Methods

Sampling

DNA extraction and fungal metabarcoding

Bioinformatics

Environmental variables

Statistical analyses

Declarations

References

Supplementary Files

Status:

Version 1