<figure><img src="image_1.png" /></figure>

HAL
open science

# Metagenomic Characterization of Indoor Dust Bacterial and Fungal Microbiota in Homes of Asthma and Non-asthma Patients Using Next Generation Sequencing

J-P Gangneux, Mohamed Sassi, Pierre Lemire, Pierre Le Cann

► To cite this version:

J-P Gangneux, Mohamed Sassi, Pierre Lemire, Pierre Le Cann. Metagenomic Characterization of In-door Dust Bacterial and Fungal Microbiota in Homes of Asthma and Non-asthma Patients Using Next Generation Sequencing. Frontiers in Microbiology, 2020, 11, pp.1671. 10.3389/fmicb.2020.01671. hal-02930076

HAL Id: hal-02930076

https://hal.science/hal-02930076v1

Submitted on 4 Sep 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
<figure><img src="image_3.png" /></figure>

Check for
updates

# Metagenomic Characterization of Indoor Dust Bacterial and Fungal Microbiota in Homes of Asthma and Non-asthma Patients Using Next Generation Sequencing

Jean-Pierre Gangneux*, Mohamed Sassi, Pierre Lemire and Pierre Le Cann

Univ Rennes, CHU Rennes, Inserm, EHESP, Institut de Recherche en Santé, Environnement et Travail (Irset) – UMR_S 1085, Rennes, France

**Background:** The exposure of house occupants to indoor air pollutants has increased in recent decades. Among microbiological contaminants, bacterial and fungal aerosols remain poorly studied and the debate on the impact of these aerosols on respiratory health is still open. This study aimed to assess the diversity of indoor microbial communities in relationship with the health of occupants.

**Methods:** Measurements were taken from dwellings of 2 cohorts in Brittany (France), one with children without any pathology and the other with children and adults with asthma. Thirty dust samples were analyzed by next generation sequencing with a 16S and 18S targeted metagenomics approach. Analysis of sequencing data was performed using qiime 2, and univariate and multivariate statistical analysis using R software and phyloseq package.

**Results:** A total of 2,637 prokaryotic (589 at genus level) and 2,153 eukaryotic taxa were identified (856 fungal taxa (39%) and 573 metazoa (26%)). The four main bacterial phyla were identified: Proteobacteria (53%), Firmicutes (27%), Actinobacteria (11%), Bacteroidetes (8%). Among Fungi, only 136 taxa were identified at genus level. Three main fungal phyla were identified: Ascomycota (84%), Basidiomycota (12%) and Mucoromycota (3%). No bacterial nor fungal phyla were significantly associated with asthma versus control group. A significant over representation in control group versus asthma was observed for Christensenellaceae family (p-value = 0.0015, adj. p-value = 0.033). Besides, a trend for over representation in control group was observed with Dermabacteraceae family (p-value = 0.0002, adj. p-value = 0.815).

**Conclusions:** Our findings provide evidence that dust samples harbor a high diversity of human-associated bacteria and fungi. Molecular methods such as next generation

OPEN ACCESS

**Edited by:**

Sara Gago,
The University of Manchester,
United Kingdom

**Reviewed by:**

Raquel Sabino,
National Institute of Health Dr. Ricardo
Jorge, Portugal
Danielle Weaver,
The University of Manchester,
United Kingdom

**\*Correspondence:**

Jean-Pierre Gangneux
jean-pierre.gangneux@univ-rennes1.fr

**Specialty section:**

This article was submitted to
Fungi and Their Interactions,
a section of the journal
Frontiers in Microbiology

**Received:** 25 February 2020

**Accepted:** 25 June 2020

**Published:** 30 July 2020

**Citation:**

Gangneux J-P, Sassi M, Lemire P and Le Cann P (2020) Metagenomic Characterization of Indoor Dust Bacterial and Fungal Microbiota in Homes of Asthma and Non-asthma Patients Using Next Generation Sequencing. *Front. Microbiol.* 11:1671. doi: 10.3389/fmicb.2020.01671
sequencing are reliable tools for identifying and tracking the bacterial and fungal diversity in dust samples, a less easy strategy for the detection of eukaryotes at least using 18S metagenomics approach. This study showed that the detection of some bacteria might be associated to indoor air of asthmatic patients. Regarding fungi, a higher number of samples and sequencing with more depth could allow reaching significant signatures.

**Keywords:** microbiota, mycobiota, asthma, indoor environment, dust, next-generation sequencing, metagenomics, fungi

INTRODUCTION

Asthma is a common chronic inflammatory airway disease in children and adults from developed countries. In France, it affects more than 10% of children and 6% of adults (Delmas and Fuhrman, 2010).

Asthma is a very complex disease characterized by enhanced bronchial reactivity leading to airway obstruction and respiratory difficulties (Martinez and Vercelli, 2013). Exposure to environmental chemical or biological factors can lead to exacerbation of asthma and prevention measures are still needed in order to reduce the exposure to these agents (Beasley et al., 2015; Shaheen, 2019). The airway microbiome in asthma patients results from numerous factors. Environmental exposure, a complex relationship with the gastrointestinal microbiome, the development of immune function, and a predisposition to allergic sensitization and asthma are among the most referred factors (Huang and Boushey, 2015).

Indoor fungal exposure has been associated to the development of asthma (Chen et al., 2014). But, according to Ege et al. (2011), children who grow up in the environments with a wide range of microbial exposures, like farming environments, are more likely protected from childhood asthma and atopy than urban children. Feng et al. (2016) reported that parental allergic diseases, atopy, diet and early life exposures might explain the higher prevalence of asthma in the urban environment. Also, there is a complex interplay between genetic predisposition and environmental exposures, including microbes and allergens (Anderson and Jackson, 2017). Regarding exposure to molds, it is estimated that in France about 125,000 individuals suffer from severe asthma with fungal sensitization (SAFS) episodes (189 cases/100,000 adults per year) (Gangneux et al., 2016).

Among biological agents, bacteria have been less studied even if endotoxins produced by gram negative bacteria have been significantly correlated with wheezing in children (Horick et al., 2006). Moreover, many studies have suggested that exposure to lipopolysaccharide (LPS), in non-rural environments, is a risk factor for increased asthma prevalence and severity of disease (Thorne et al., 2005).

Culture-based studies have been in use for a long time for the detection of various microorganisms indoor (Rintala et al., 2008, 2012; Méheust et al., 2014). However, while culture isolation focuses on the presence of particular bacteria or fungi, next generation sequencing (NGS) studies unveil microbial communities comprising thousands of uncultured microbes. Recent advances in microbiological methods show promise in determining whether the microbes recovered during indoor

sampling campaigns may in fact be causative agents as such, or only surrogates of an underlying harmful or beneficial exposure. Given the lack of knowledge on the specificities of the relevant exposures, assessments at this stage need to be kept broad and comprehensive (Mensah-Attipoe et al., 2017).

According to Singanayagam et al. (2017), who underlined the role of microbes in the exacerbation of asthma, the emergence of next generation molecular sequencing techniques to characterize the microbiota has facilitated renewed interest. They concluded that longitudinal studies that characterize the changes in lower respiratory tract microbiota from birth up to development of asthma are currently unavailable. Additionally, few studies have focused on the interactions between bacteria and the immune system in driving development of nonatopic asthma phenotypes.

Thus, despite the recognition of the importance of microbial exposure for human health, the precise role of microbes in

<figure><img src="image_4.png" /></figure>
the development and exacerbation of respiratory symptoms and allergies remains poorly understood.

The objective of our study was to use 16S and 18S rRNA sequencing approaches to characterize and identify indoor bacterial and eukaryotic communities that may be associated with the exacerbation of asthma. This pilot study compared two groups of asthmatic and non-asthmatic people's dwellings.

# MATERIALS AND METHODS

## Study Population and Home Visits

Thirty dust samples were collected from dwellings of 2 French cohorts and were analyzed by next-generation sequencing:

(i) 15 dwellings of asthma patients (ASTHMA) from The Ecenvir cohort, which aims to evaluate the clinical and the economic impact of Indoor Environment Counselor (IEC) on the symptoms of severe asthma. Patients answered

a questionnaire during an interview with the IEC on dwellings characteristics (localization in cities or rural, apartment or house, ...) and cleaning habits of patients. Asthma was graded with the GINA (Global Initiative for Asthma) scale and patients had three clinical check-ups.

(ii) 15 control dwellings of control patients (CONTROL) from the Pélagie AC cohort, which follows more than 3,500 pregnant women in Brittany from pregnancy to age of 20 of children. The Pélagie AC project aims to assess the effects on children respiratory health of chemicals and biological pollutants exposure. One of the check-up of Pélagie AC cohort occurs once the children are 6 years old. A subset of families without asthma answered a detailed questionnaire about their dwellings characteristics and life for the Pélagie AC study.

According to the French Public Health Law, protocols of this type are exempt from the requirement for formal informed consent. However, The INSERM and Rennes Ethical Committees

<figure><img src="image_5.png" /><figcaption>**FIGURE 2** | Hierarchical clustering analysis: Dendrogram of Bray-Curtis dissimilarity matrices between samples based on 16S data at feature level with UPGMA method.</figcaption></figure>
approved the protocol, patients were informed and it was possible for them to refuse environmental samples for both cohorts:

- The Pélagie AC cohort, with dwellings being part of a larger study as stated in Dallongeville et al. (2015) (French Consulting Committee for the Treatment of Information in Medical Research (no. 09.485),

- The Ecenvir cohort, with dwellings being part of a larger study as approved by the Comité de Protection des Personnes de Rennes Ouest V (8 january 2013, N°ID RCB 2012-A01414-39).

## Dust Collection

Dust samples were collected by vacuuming the floor in the child's bedroom using a Duststream Collector™ (Indoor biotechnologies, United Kingdom) sampler-fitted vacuum cleaner (40 µm mesh nylon filter, domestic vacuum cleaner) as described by Dallongeville et al. (2015). Briefly, sampling was preferentially carried out on the carpets (when available) or on hard surface floors until filling at least 2/3 of the filter-sampler, in order to ensure collection of at least 10 mg of fine dust. Collected dust was sieved at 300 µm to discard debris, and 10 mg of fine dust were resuspended in 1 mL of PBS-Tween (0.05%) buffer, shaken at 800 rpm for 1 h and stored at -20°C until analyzed.

## DNA Extraction, Sequencing, and Analysis

Total DNA was extracted from samples using the method developed by Godon et al. (1997). A portion of the 16S rRNA was amplified using the barcoded, universal primer set (515WF/918WR) (Wang et al., 2009). DNA extraction was carried out according to the guanidium thiocyanate method. Positive extraction control is included. A portion of the 18S rRNA was amplified using universal 18S rRNA primer set (574WF/952WR) (Hadziavdic et al., 2014). PCR reactions were performed using AccuStart II PCR ToughMix kit and cleaned (HighPrep PCR beads, Mokascience). The amplification program consists of 28 PCR cycles with the following steps: 94°C for 3 min, followed by 28 cycles of 94°C for 30 s, 53°C for 40 s and 72°C for 1 min, followed by a final elongation step at 72°C for 5 min. PCR positive and negative controls are included. Pools were submitted for sequencing on Illumina MiSeq instrument at GeT-PlaGe (Auzeville, France). Sequencing read length was 250 bp. Sequences were processed using Mothur (version 1.33.1) according to MiSeq SOP pipeline (Schloss et al., 2009).

Barcodes, primers, and sequences showing homopolymers of more than 8 bp have been discarded. Sequences showing 100% homology were grouped in unique sequences, then in OTUs (operational taxonomic unit, based on 97% homology). Sequences were next assigned to match to a sequence in SILVA (version 123) and PR2 databases for prokaryotes and eukaryotes respectively to identify the genus level. We considered assigned sequences with relative abundance (RA) higher than 0.002 %.

<table><caption>TABLE 1 | Abundances and percentages of the bacterial communities in the environmental samples of the control and asthma groups.</caption><thead><tr><th rowspan="2">Taxa</th><th colspan="2">Control</th><th colspan="2">Asthma</th><th colspan="2">Classical Univariate Statistical Comparisons</th><th colspan="2">Linear Discriminant Analysis (LDA) Effect Size (LEfSe)</th><th colspan="4">Deseq method</th></tr><tr><th>Abundance</th><th>Percentage</th><th>Abundance</th><th>Percentage</th><th>Pvalues</th><th>FDR</th><th>Statistics</th><th>Pvalues</th><th>FDR</th><th>Asthma</th><th>Control</th><th>LDAscore</th><th>log2FC</th><th>IfcSE</th><th>Pvalues</th><th>FDR</th></tr></thead><tbody><tr><td>Proteobacteria</td><td>39305</td><td>42.18</td><td>67401</td><td>58.84</td><td>0.086611</td><td>0.27995</td><td>1.7759</td><td>0.071185</td><td>0.23728</td><td>5335600</td><td>4209300</td><td>-5.75</td><td>-0.54907</td><td>0.36825</td><td>0.13595</td><td>0.45318</td></tr><tr><td>Firmicutes</td><td>35733</td><td>38.35</td><td>27883</td><td>24.34</td><td>0.020982</td><td>0.10491</td><td>-2.4461</td><td>0.029436</td><td>0.14718</td><td>2639100</td><td>3859700</td><td>5.79</td><td>0.17782</td><td>0.2959</td><td>0.54788</td><td>0.93211</td></tr><tr><td>Actinobacteria</td><td>12415</td><td>13.32</td><td>10543</td><td>9.2</td><td>0.25007</td><td>0.50014</td><td>-1.1745</td><td>0.2902</td><td>0.72549</td><td>1052800</td><td>1317800</td><td>5.12</td><td>0.1867</td><td>0.28297</td><td>0.5084</td><td>0.93211</td></tr><tr><td>Bacteroidetes</td><td>4421</td><td>4.74</td><td>7542</td><td>6.58</td><td>0.11198</td><td>0.27995</td><td>1.641</td><td>0.60413</td><td>0.79302</td><td>850600</td><td>472570</td><td>-5.28</td><td>-0.90264</td><td>0.389</td><td>0.02032</td><td>0.1016</td></tr><tr><td>Cyanobacteria</td><td>790</td><td>0.85</td><td>791</td><td>0.69</td><td>0.93632</td><td>0.93632</td><td>0.08062</td><td>0.94825</td><td>0.94825</td><td>89040</td><td>86065</td><td>-3.17</td><td>-0.16357</td><td>1.1595</td><td>0.89782</td><td>0.98647</td></tr><tr><td>Acidobacteria</td><td>118</td><td>0.13</td><td>173</td><td>0.15</td><td>0.90185</td><td>0.93632</td><td>0.12445</td><td>0.71613</td><td>0.7957</td><td>13756</td><td>12537</td><td>-2.79</td><td>-0.54658</td><td>2.2049</td><td>0.80422</td><td>0.98647</td></tr><tr><td>Fusobacteria</td><td>101</td><td>0.11</td><td>154</td><td>0.13</td><td>0.67713</td><td>0.84641</td><td>0.42078</td><td>0.60465</td><td>0.79302</td><td>12344</td><td>7576.6</td><td>-3.38</td><td>0</td><td>0</td><td>1</td><td>1</td></tr><tr><td>Deinococcus_Thermus</td><td>73</td><td>0.08</td><td>51</td><td>0.04</td><td>0.5005</td><td>0.73593</td><td>-0.68254</td><td>0.63442</td><td>0.79302</td><td>4478</td><td>8131.4</td><td>3.26</td><td>1.2533</td><td>2.1463</td><td>0.55927</td><td>0.93211</td></tr><tr><td>Epsilonbacteraeota</td><td>53</td><td>0.06</td><td>12</td><td>0.01</td><td>0.51515</td><td>0.73593</td><td>-0.65921</td><td>0.55026</td><td>0.79302</td><td>2283.1</td><td>6084</td><td>3.28</td><td>-0.85519</td><td>2.9528</td><td>0.77211</td><td>0.98647</td></tr><tr><td>Others</td><td>172</td><td>0.18</td><td>0</td><td>0</td><td>0.015326</td><td>0.10491</td><td>-2.5826</td><td>0.0075619</td><td>0.075619</td><td>0</td><td>20294</td><td>4.01</td><td>5.5643</td><td>2.1585</td><td>0.0099431</td><td>0.099431</td></tr></tbody></table>
<figure><img src="image_6.png" /></figure>

## Sequencing Data Analysis

Analysis of sequencing data was performed using QIIME 2¹, and univariate and multivariate statistical analysis using R software and phyloseq package.

## RESULTS

### Comparison of Microbiota in the Environmental Samples

The sequencing of the 16S rRNA allowed the identification of 2,637 bacterial taxa, 589 at genus level with a predominance of Proteobacteria (53%), Firmicutes (27%) and Actinobacteria (11%). Regarding 18S, 2,153 eukaryotic taxa were identified: 856 fungal taxa (39%) and 573 metazoa (26%). Among Fungi, only 136 taxa were identified at genus level. The 16S and 18S sequencing data have been uploaded at EBI Metagenomics under the accession numbers PRJEB37043 and PRJEB37050 respectively.

### Bacteria

The Shannon index was used to evaluate the intra-group diversity in the environmental samples of the two groups; it was 3.68 ± 0.67 in the control group and 3.94 ± 0.22 in the asthma group, showing no significant difference in diversity between the two groups (Mann-Whitney statistical *p* = 0.5393) (Figure 1). The beta-diversity results using

ANOSIM method indicated that the bacterial composition in the two groups was different as the global *R*-value was 0.109, with a *p*-value of 0.004. Hierarchical clustering analysis (Figure 2) indicated no distinct difference in the OTUs for the environmental samples of the two groups. There was no significant difference between the compositions of the bacterial communities in the environmental samples of the two groups, when analyzed at feature level (Supplementary Figure 1). When analyzed at the phylum level the abundance for the phylum Firmicutes was 38.35% in the control group, significantly higher than and 24.34% in the asthma group (24.34%, *p* = 0.02) (Table 1). Further, Proteobacteria, Actinobacteria and Bacteroides showed an abundance of 42.18, 13.32, and 4.74% in the control group and 58.84, 9.20, and 6.58% in the asthma group, respectively, with no significant difference between the compositions of the bacterial communities in the environmental samples of the two groups (Table 1, Figure 3 and Supplementary Figure 1).

### Fungi

The Shannon index was used to evaluate the intra-group diversity in the environmental samples of the two groups; it was 3.00 ± 0.77 in the control group and 2.89 ± 0.64 in the asthma group, showing no significant difference in diversity between the two groups (Mann-Whitney statistical *p* = 0.66952) (Figure 4). The beta-diversity results using ANOSIM method indicated that the fungal composition in the two groups was not different as the global *R*- value was -0.003, with a *p*-value of 0.859. There was no significant difference between the

<footnote n="1">¹https://qiime2.org</footnote>
<figure><img src="image_7.png" /><figcaption>**FIGURE 4** | Fungal alpha-diversity, measured by Shannon diversity Index is plotted for samples with asthma (red) and controls (blue). The line inside the box represents the median, while the whiskers represent the lowest and highest values within the 1.5 interquartile range (IQR). Outliers as well as individual sample values are shown as dots. Mann-Whitney statistical testing showed no significant difference in diversity between the two groups (pShannon = 0.66952).</figcaption></figure>
compositions of the fungal communities in the environmental samples of the two groups, when analyzed at feature level (Supplementary Figure 2) and phylum level (Table 2). For the phylum Ascomycota, the abundance was 82.21% in the control group and 89.18% in the asthma group. Furthermore, Basidiomycota and Mucoromycota showed an abundance of 12.09 and 3.92% in the control group and 9.03 and 1.24% in the asthma group, respectively (Figure 5).

## Differential Abundance Analysis

Bacterial and fungal taxon-levels (OTU) differences in relative abundance were examined between the asthma and control groups. A significant over representation of specific bacteria was observed in control patients' group for Christensenellaceae family (p-value = 0.0015, adj. p-value = 0.033) (Figure 6) and only a trend for Dermabacteraceae family (p-value = 0.0002, adj. p-value = 0.065) (Figure 7).

## DISCUSSION

Health effects of environmental microbes are still debated and not completely understood. A first approach is to examine relationships between genes, microbes in airways and gut, and the environment in asthma causes. This is a complex, dynamic and very heterogeneous process still unresolved (Huang and Boushey, 2015). Besides, healthy environments may have impact on the exacerbation of asthma and the usefulness of allergen avoidance on asthma control have shown various evidence of their efficacy on the clinical improvement of patients (Le Cann et al., 2017; Gangneux et al., 2020). In this work, our aim was to compare microbial communities in asthma patients dwellings compared to non-asthma control homes, and to characterize and identify indoor bacterial and eukaryotic communities that may be associated with the exacerbation of asthma. In the literature, exposure to low fungal and bacterial richness in house dust is associated with an increased risk of asthma development (Ege et al., 2011; Dannemiller et al., 2013). Air, dust and surface sampling strategies to detect bacteria and fungi have alternately advantages and limits, as well as the different methods of detection, identification and quantification (Méheust et al., 2014; Gangneux et al., 2019). Because our objective for this work was to use NGS, we decided to sample dust by vacuuming the floor rather than aerosols in order to gain in sensitivity. Another valuable technical option described in the literature relies on the use of electrostatic dust cloths for further molecular, immunological or cultural analysis (Cox et al., 2017; Kristono et al., 2019).

In our work, the alpha- and beta-diversity results indicated that the bacterial and fungal composition in the two groups were not different, without significant difference when analyzed at the phylum level. On a global point of view, our results are in accordance to others as the main phylum detected in both group is Proteobacteria, then Firmicutes, Actinobacteria, Bacteroidetes (Hewitt et al., 2012), with Firmicutes over represented in the control group compared to asthma dwellings. Besides, there was a trend of more Proteobacteria detected in asthma population as

<table><caption>TABLE 2 | Abundances and percentages of the fungal communities in the environmental samples of the control and asthma groups.</caption><thead><tr><th rowspan="2">Taxa</th><th colspan="2">Control</th><th colspan="2">Asthma</th><th colspan="2">Classical Univariate Statistical Comparisons</th><th colspan="2">Linear Discriminant Analysis (LDA) Effect Size (LEfSe)</th><th colspan="2">Deseq method</th></tr><tr><th>Abundance</th><th>Percentage</th><th>Abundance</th><th>Percentage</th><th>Pvalues</th><th>FDR</th><th>Statistics</th><th>Pvalues</th><th>FDR</th><th>Asthma</th><th>Control</th><th>LDAScore</th><th>log2FC</th><th>IfcSE</th><th>Pvalues</th><th>FDR</th></tr></thead><tbody><tr><td>Ascomycota</td><td>186492</td><td>82.21</td><td>189747</td><td>89.18</td><td>0.11207</td><td>0.31014</td><td>142</td><td>0.10635</td><td>0.29345</td><td>8744400</td><td>8042200</td><td>-5.55</td><td>-0.39748</td><td>0.40435</td><td>0.3256</td><td>0.51503</td></tr><tr><td>Basidiomycota</td><td>20078</td><td>12.09</td><td>17194</td><td>9.03</td><td>0.18608</td><td>0.31014</td><td>74</td><td>0.17607</td><td>0.29345</td><td>996370</td><td>1388100</td><td>5.29</td><td>0.22401</td><td>0.22313</td><td>0.31541</td><td>0.51503</td></tr><tr><td>Mucoromycota</td><td>6515</td><td>3.92</td><td>2357</td><td>1.24</td><td>0.91251</td><td>0.91251</td><td>108</td><td>0.89511</td><td>0.89511</td><td>201920</td><td>358010</td><td>4.89</td><td>-1.2848</td><td>0.99581</td><td>0.19698</td><td>0.51503</td></tr><tr><td>Chytridiomycota</td><td>2564</td><td>1.54</td><td>691</td><td>0.36</td><td>0.16554</td><td>0.31014</td><td>75</td><td>0.15849</td><td>0.29345</td><td>33275</td><td>185200</td><td>4.88</td><td>0.99249</td><td>1.4868</td><td>0.50444</td><td>0.51503</td></tr><tr><td>Others</td><td>376</td><td>0.23</td><td>352</td><td>0.18</td><td>0.33805</td><td>0.42257</td><td>123.5</td><td>0.32481</td><td>0.40601</td><td>24060</td><td>26473</td><td>3.08</td><td>-1.2974</td><td>1.9928</td><td>0.51503</td><td>0.51503</td></tr></tbody></table>
<figure><img src="image_8.png" /><figcaption>**FIGURE 5** | Proportion of fungal community at phylum level in the control and asthma groups at the phylum level. \*0.01 < *p* ≤ 0.05, based on Mann-Whitney *U*-test.</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>**FIGURE 6** | Box plot of differential abundance of Christensenellaceae family in the control and asthma groups using Classical Univariate Statistical Comparisons method based on Mann-Whitney *U*-test.</figcaption></figure>

previously described (Ciaccio et al., 2015). This was associated with a higher abundance of Proteobacteria in airways microbiota of asthmatic people (Hilty et al., 2010). However, it is very difficult to compare the studies because usually the populations studied are different (sometimes urban, sometimes rural, etc.). At the family level, a signature in differential abundance was observed with significantly more Christensenellaceae detected in control dwellings than in asthma homes. Regarding other genus and family, no statistical difference was observed even if much Dermatobacteriacae were also detected in control dwellings.

Some limitations of this work can be identified. A lack of statistical power is possible, linked to the limited number of dwellings sampled. However, this study opens new perspectives with interesting trends that may drive further studies on the topic. Also, it is possible that we have underestimated the fungal burden and diversity in dust samples because of a potential loss of fungal spores as shown for aerosols (Mbareche et al., 2019). The interest of concentrating particles before DNA extraction needs to be investigated on dust samples as it has been studied for air samples. Finally, we intentionally decided to amplify and
<figure><img src="image_10.png" /><figcaption>**FIGURE 7** | Box plot of differential abundance of Dermabacteraceae family in the control and asthma groups using Classical Univariate Statistical Comparisons method based on Mann-Whitney *U*-test.</figcaption></figure>

sequence rRNA of 18S gene in order to detect fungi but also eukaryotes cells that may be involved in asthma sensitization such as cockroach, house dust mites or pets. However, we were very disappointed by the low resolution for species identification of non fungal eukaryotes and were not able to analyze them. Thus, if looking only to fungi, the question of using ITS1-2 or 28S rRNA targets could have been asked, while advantages and limits of each target and even primer sets are still debated (Hanson et al., 2016; De Filippis et al., 2017; Frau et al., 2019). Another option is to use shotgun sequencing for taxon identification and abundance assessment for both prokaryotic and eukaryotic sequences. New tools such as NGS should be more standardized from one study to another and include processing steps of quality controls in order to avoid low-quality reads and sequences and to optimize abundance thresholds and profiles. However, deep sequencing really opens new perspectives on the comprehension of the microbial exposome and its impact on airway microbiome.

## DATA AVAILABILITY STATEMENT

The datasets generated for this study are available on request to the corresponding author.

## ETHICS STATEMENT

According to the French Public Health Law, protocols of this type are exempt from the requirement for formal informed consent. However, they were informed and it was possible to refuse environmental samples for both cohorts.

## AUTHOR CONTRIBUTIONS

All authors listed have made a substantial, direct and intellectual contribution to the work, and approved it for publication.

## FUNDING

We thank the French ANSES PNR EST for the financial support of this project EIIBE - EST-2015/1/148.

## ACKNOWLEDGMENTS

The authors would like to thank Cécile Chevrier and Nathalie Costet from Epidemiology team of IRSET for giving the opportunity to have access to the dwellings of the Pelagie cohort. This work was presented as a poster presentation during the 9th Trends in Medical Mycology (TIMM) congress, in Nice, France (11–14 October 2019).

## SUPPLEMENTARY MATERIAL

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

**FIGURE S1** | Heatmap based on 16S analysis. The Euclidean dissimilarity matrix and Ward's hierarchical clustering algorithms were applied in this analysis. Rows indicate feature level (OTUs) and columns indicate the different samples.

**FIGURE S2** | Heatmap based on 18S analysis. The Euclidean dissimilarity matrix and Ward's hierarchical clustering algorithms were applied in this analysis. Rows indicate feature level (OTUs) and columns indicate the different samples.
REFERENCES

Anderson, H. M., and Jackson, D. J. (2017). Microbes, allergic sensitization, and the natural history of asthma. *Curr. Opin. Allergy Clin. Immunol.* 17, 116–122. doi: 10.1097/ACI.00000000000000338 doi: 10.1097/aci.00000000000000338

Beasley, R., Semprini, A., and Mitchell, E. A. (2015). Risk factors for asthma: is prevention possible? *Lancet* 386, 1075–1085. doi: 10.1016/S0140-6736(15) 00156-157 doi: 10.1016/s0140-6736(15)00156-7

Chen, C. H., Chao, H. J., Chan, C. C., Chen, B. Y., and Guo, Y. L. (2014). Current asthma in schoolchildren is related to fungal spores in classrooms. *Chest* 146, 123–134. doi: 10.1378/chest doi: 10.1378/chest.13-2129

Ciaccio, C. E., Barnes, C., Kennedy, K., Chan, M., Portnoy, J., and Rosenwasser, L. (2015). Home dust microbiota is disordered in homes of low-income asthmatic children. *J. Asthma* 52, 873–880. doi: 10.3109/02770903.2015.1028076

Cox, J., Indugula, R., Vesper, S., Zhu, Z., Jandarov, R., and Reponen, T. (2017). Comparison of indoor air sampling and dust collection methods for fungal exposure assessment using quantitative PCR. *Environ Sci. Process Impacts* 19, 1312–1319. doi: 10.1039/c7em00257b

Dallongeville, A., Le Cann, P., Zmirou-Navier, D., Chevrier, C., Costet, N., Annesi-Maesano, I., et al. (2015). Concentration and determinants of molds and allergens in indoor air and house dust of French dwellings. *Sci. Total Environ.* 536, 964–972. doi: 10.1016/j.scitotenv.2015.06.039

Dannemiller, K. C., Murphy, J. S., Dixon, S. L., Pennell, K. G., Suuberg, E. M., Jacobs, D. E., et al. (2013). Formaldehyde concentrations in household air of asthma patients determined using colorimetric detector tubes. *Indoor Air.* 23, 285–294. doi: 10.1111/ina.12024

De Filippis, F., Laiola, M., Blaiotta, G., and Ercolini, D. (2017). Different amplicon targets for sequencing-based studies of fungal diversity. *Appl. Environ. Microbiol.* 83, e905–e917. doi: 10.1128/AEM.00905-917

Delmas, M. C., and Fuhrman, C. (2010). Asthma in France: a review of descriptive epidemiological data. *Rev. Mal. Respir.* 27, 151–159. doi: 10.1016/j.rmr.2009.09.001

Ege, M. J., Mayer, M., Normand, A. C., Genuneit, J., Cookson, W. O. C. M., Braun-Fahrländer, C., et al. (2011). Exposure to environmental microorganisms and childhood asthma. *N. Engl. J. Med.* 364, 701–709. doi: 10.1056/NEJMoa1007302

Feng, M., Yang, Z., Pan, L., Lai, X., Xian, M., Huang, X., et al. (2016). Associations of early life exposures and environmental factors with asthma among children in rural and urban areas of Guangdong, China. *Chest* 149, 1030–1041. doi: 10.1016/j.chest.2015.12.028

Frau, A., Kenny, J. G., Lenzi, L., Campbell, B. J., Ijaz, U. Z., Duckworth, C. A., et al. (2019). DNA extraction and amplicon production strategies deeply influence the outcome of gut mycobiome studies. *Sci. Rep.* 9:9328.

Gangneux, J. P., Bougnoux, M. E., Hennequin, C., Godet, C., Chandenier, J., Denning, D. W., et al. (2016). An estimation of burden of serious fungal infections in France. *J. Med. Mycol.* 26, 385–390.

Gangneux, J. P., Bouvrais, M., Frain, S., Morel, H., Deguen, S., Chevrier, C., et al. (2020). Asthma and indoor environment: usefulness of a global allergen avoidance method on asthma control and exposure to molds. *Mycopathologia* 185, 367–371. doi: 10.1007/s11046-019-00417-419

Gangneux, J. P., Guegan, H., Vandenborght, L. E., Buffet-Bataillon, S., Enaud, R., and Delhaes, L. (2019). A European ECMM-ESCMID survey on goals and practices for mycobiota characterisation using next-generation sequencing. *Mycoses* 62, 1096–1099. doi: 10.1111/myc. 12999

Godon, J. J., Zumstein, E., Dabert, P., Habouzit, F., and Moletta, M. (1997). Molecular microbial diversity of an anaerobic digestor as determined by small-subunit rDNA sequence analysis. *Appl. Environ. Microbiol.* 63, 2802–2813. doi: 10.1128/aem.63.7.2802-2813.1997

Hadziavdic, K., Lekang, K., Lanzen, A., Jonassen, I., Thompson, E. M., and Troedsson, C. (2014). Characterization of the 18S rRNA gene for designing universal eukaryote specific primers. *PLoS One* 9:e87624. doi: 10.1371/journal.pone.0087624

Hanson, B., Zhou, Y., Bautista, E. J., Urch, B., Speck, M., Silverman, F., et al. (2016). Characterization of the bacterial and fungal microbiome in indoor dust and outdoor air samples: a pilot study. *Environ. Sci. Process Impacts* 18, 713–724. doi: 10.1039/c5em00639b

Hewitt, K. M., Gerba, C. P., Maxwell, S. L., and Kelley, S. T. (2012). Office space bacterial abundance and diversity in three metropolitan areas. *PLoS One* 7:e37849. doi: 10.1371/journal.pone.0037849

Hilty, M., Burke, C., Pedro, H., Cardenas, P., Bush, A., Bossley, C., et al. (2010). Disordered microbial communities in asthmatic airways. *PLoS One* 5:e8578. doi: 10.1371/journal.pone.0008578

Horick, N., Weller, E., Milton, D. K., Gold, D. R., Li, R., and Spiegelman, D. (2006). Home endotoxin exposure and wheeze in infants: correction for bias due to exposure measurement error. *Environ. Health Perspect.* 114, 135–140. doi: 10.1289/ehp.7981

Huang, Y. J., and Boushey, H. A. (2015). The microbiome in asthma. *J Allergy Clin. Immunol.* 135, 25–30.

Kristono, G. A., Shorter, C., Pierse, N., Crane, J., and Siebers, R. (2019). Endotoxin, cat, and house dust mite allergens in electrostatic cloths and bedroom dust. *J. Occup. Environ. Hyg.* 16, 89–96. doi: 10.1080/15459624.2018.1536827

Le Cann, P., Paulus, H., Glorennec, P., Le Bot, B., Frain, S., and Gangneux, J. P. (2017). Home environmental interventions for the prevention or control of allergic and respiratory diseases: what really works. *J. Allergy Clin. Immun. Pract.* 5, 66–79. doi: 10.1016/j.jaip.2016.07.011

Martinez, F. D., and Vercelli, D. (2013). Asthma. *Lancet* 1982, 1360–1372. doi: 10.1016/S0140-6736(13)61536-61536

Mbareche, H., Veillette, M., Teertstra, W., Kegel, W., Bilodeau, G. J., Wösten, H. A. B., et al. (2019). Recovery of fungal cells from air samples: a tale of loss and gain. *Appl. Environ. Microbiol.* 85:e02941-18.

Méheust, D., Le Cann, P., Reboux, G., Millon, L., and Gangneux, J. P. (2014). Indoor fungal contamination: health risks and measurement methods in hospitals, homes and workplaces. *Crit. Rev. Microbiol.* 40, 248–260. doi: 10. 3109/1040841X.2013.777687

Mensah-Atipoe, J., Täubel, M., Hernandez, M., Pitkäranta, M., and Reponen, T. (2017). An emerging paradox: toward a better understanding of the potential benefits and adversity of microbe exposures in the indoor environment. *Indoor Air* 27, 3–5. doi: 10.1111/ina.12344

Rintala, H., Pitkäranta, M., and Täubel, M. (2012). Microbial communities associated with house dust. *Adv. Appl. Microbiol.* 78, 75–120. doi: 10.1016/ B978-0-12-394805-2.00004-X

Rintala, H., Pitkäranta, M., Toivola, M., Paulin, L., and Nevalainen, A. (2008). Diversity and seasonal dynamics of bacterial community in indoor environment. *BMC Microbiol.* 8:56. doi: 10.1186/1471-2180-8-56 doi: 10.1186/ 1471-2180-8-56

Schloss, P. D., Westcott, S. L., Ryabin, T., Hall, J. R., Hartmann, M., Hollister, E. B., et al. (2009). Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. *Appl. Environ. Microbiol.* 75, 7537–7541. doi: 10.1128/AEM. 01541-09

Shaheen, S. (2019). Elucidating the causes of asthma: how can we do better? *Lancet Respir. Med.* 7:e25. doi: 10.1016/s2213-2600(19)30225-5

Singanayagam, A., Ritchie, A. I., and Johnston, S. L. (2017). Role of microbiome in the pathophysiology and disease course of asthma. *Curr. Opin. Pulm. Med.* 23, 41–47. doi: 10.1097/mcp.00000000000000333

Thorne, P. S., Kulhánková, K., Yin, M., Cohn, R., Arbes, S. J. Jr., and Zeldin, D. C. (2005). Endotoxin exposure is a risk factor for asthma: the national survey of endotoxin in United States housing. *Am. J. Respir. Crit. Care Med.* 172, 1371–1377. doi: 10.1164/ccm.200505-758oc

Wang, Y., Morimoto, S., Ogawa, N., Oomori, T., and Fujii, T. (2009). An improved method to extract RNA from soil with efficient removal of humic acids. *J. Appl. Microbiol.* 107, 1168–1177. doi: 10.1111/j.1365-2672.2009.04298.x

**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.

Copyright © 2020 Gangneux, Sassi, Lemire and Le Cann. 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.