Environmental Mold and Mycotoxin Exposures Elicit Specific Cytokine and Chemokine Responses

Abstract

Background

Molds can cause respiratory symptoms and asthma. We sought to use isolated peripheral blood mononuclear cells (PBMCs) to understand changes in cytokine and chemokine levels in response to mold and mycotoxin exposures and to link these levels with respiratory symptoms in humans. We did this by utilizing an ex vivo assay approach to differentiate mold-exposed patients and unexposed controls. While circulating plasma chemokine and cytokine levels from these two groups might be similar, we hypothesized that by challenging their isolated white blood cells with mold or mold extracts, we would see a differential chemokine and cytokine release.

Methods and Findings

Peripheral blood mononuclear cells (PBMCs) were isolated from blood from 33 patients with a history of mold exposures and from 17 controls. Cultured PBMCs were incubated with the most prominent Stachybotrys chartarum mycotoxin, satratoxin G, or with aqueous mold extract, ionomycin, or media, each with or without PMA. Additional PBMCs were exposed to spores of Aspergillus niger, Cladosporium herbarum and Penicillium chrysogenum. After 18 hours, cytokines and chemokines released into the culture medium were measured by multiplex assay. Clinical histories, physical examinations and pulmonary function tests were also conducted. After ex vivo PBMC exposures to molds or mycotoxins, the chemokine and cytokine profiles from patients with a history of mold exposure were significantly different from those of unexposed controls. In contrast, biomarker profiles from cells exposed to media alone showed no difference between the patients and controls.

Conclusions

These findings demonstrate that chronic mold exposures induced changes in inflammatory and immune system responses to specific mold and mycotoxin challenges. These responses can differentiate mold-exposed patients from unexposed controls. This strategy may be a powerful approach to document immune system responsiveness to molds and other inflammation-inducing environmental agents.

About the Authors:

Jamie H. Rosenblum Lichtenstein

* E-mail: [email protected]

Affiliations Molecular and Integrative Physiological Sciences Program, Department of Environmental Health, Harvard T.H. Chan School of Public Health, Boston, MA, United States of America, Institute for Liberal Arts, Emerson College, Boston, MA, United States of America

Yi-Hsiang Hsu

Affiliations Molecular and Integrative Physiological Sciences Program, Department of Environmental Health, Harvard T.H. Chan School of Public Health, Boston, MA, United States of America, Hebrew Senior Life Institute for Aging Research and Harvard Medical School, Boston, MA, United States of America, Broad Institute of MIT and Harvard, Cambridge, MA, United States of America

Igor M. Gavin

Affiliation: Research Resources Center, University of Illinois at Chicago, Chicago, IL, United States of America

Thomas C. Donaghey

Affiliation: Molecular and Integrative Physiological Sciences Program, Department of Environmental Health, Harvard T.H. Chan School of Public Health, Boston, MA, United States of America

Ramon M. Molina

Affiliation: Molecular and Integrative Physiological Sciences Program, Department of Environmental Health, Harvard T.H. Chan School of Public Health, Boston, MA, United States of America

Khristy J. Thompson

Affiliation: Molecular and Integrative Physiological Sciences Program, Department of Environmental Health, Harvard T.H. Chan School of Public Health, Boston, MA, United States of America

Chih-Lin Chi

Affiliation: Institute for Health Informatics, University of Minnesota School of Nursing, Minneapolis, MN, United States of America

Bruce S. Gillis

Affiliation: Department of Medicine, University of Illinois College of Medicine, Chicago, IL, United States of America

Joseph D. Brain

Affiliation: Molecular and Integrative Physiological Sciences Program, Department of Environmental Health, Harvard T.H. Chan School of Public Health, Boston, MA, United States of America

Introduction

Indoor environments contaminated by molds cause adverse human health effects [1]. Cellulose, when combined with moisture and warm temperatures, promotes mold growth. Chronic mold exposures at home, work or school are associated with increased upper and lower respiratory symptoms [2]. This is commonly attributed to an allergen-dependent pathway, but there is evidence that mold may also trigger asthma in an allergen-independent manner [3, 4]. Mold and mold-related odors are an important cause of atopic symptoms, allergic sensitization and asthma [5]. Mold exposures have been shown to cause a nine-fold increase in emergency room visits for asthma symptoms among asthmatics [6]. Mold exposures may increase sensitivity to commonly inhaled microorganisms and inert substances and increase risks of secondary infections [7].

Both mold hyphae and conidia induce immune responses in humans [8]; however, no reliable tests linking clinical symptoms with exposures have been reported. Most clinical studies have used self-reported symptoms and were based on subjective complaints prone to bias and confounders [9, 10]. Existing clinical tests also fail to establish a definitive link between chronic mold exposures and adverse health effects. More reliable mold-related tests are needed [11]. Common immunological tests, such as IgE measurements (RAST analysis) or skin prick tests, are poor indicators of mold exposure [12, 13].

We looked for immune system modulators that might link molds and mold-related substances with respiratory illness in humans. The mold Stachybotrys chartarum (S. chartarum) can be found in both indoor and outdoor environments [14]. Modern building materials, such as cellulose-based wallboard and ceiling tiles, have made S. chartarum more common in indoor environments [15]. Severe weather such as Hurricane Katrina can cause water intrusions in homes, offices and schools; the resulting wet substrates foster the growth of S. chartarum [16, 17]. Nevertheless, there are limited studies of humans focusing on inflammatory responses [18] or biomarkers of exposures [19] to S. chartarum [20].

S. chartarum causes symptoms such as runny nose, cough, headache, and asthma exacerbations [21]. Several studies correlated infant idiopathic pulmonary hemorrhage with S. chartarum exposures in homes [22–25], though a subsequent CDC report [26] noted some flaws in the initial reports. There exists strong evidence that S. chartarum causes acute inflammatory responses, macrophage cytotoxicity [27], pulmonary hemorrhage [28], lung inflammation [29, 30] and asthma-like responses [31] in mice. S. chartarum has also been correlated with asthma in children exposed at school [32]. Acute responses to S. chartarum are usually associated with mycotoxins, including the trichothecene satratoxin G (SG), which interferes with RNA synthesis and leads to apoptosis [33–35]. We have shown that Stachybotrys spore toxin (SST), a mixture of mycotoxins, causes pulmonary hemorrhage as well as cytokine and chemokine production in a murine model [27].

To better understand the human health risks associated with mold exposures and to identify strategies to document health consequences from indoor mold exposures, we studied responses of peripheral blood mononuclear cells (PBMCs) from individuals who were exposed to molds in their workplaces. We hypothesized that chronic exposures to molds may induce tolerance and/or sensitization to these allergens, thereby decreasing some allergen-specific immune responses while increasing others. Mouse models show that chronic allergen exposure sometimes creates both tolerance, marked by suppression of some inflammatory responses such as eosinophilia, and sensitization, marked by GM-CSF expression and dendritic and CD4+ T-cell activation [36].

Cytokines and chemokines are mediators of inflammatory and immune system responses. Many macrophage and epithelial cell lines produce cytokines and chemokines after challenges to mold spores or hyphae [28, 37–39]. These facts provided additional rationale for this study. Since tolerance and sensitization are induced at the cellular level [40], we hypothesized that repeated mold exposures might cause specific persistent alterations in the expression of cytokines and chemokines when isolated PBMCs respond to subsequent exposures to specific molds, their mycotoxins, or other immune challenges. The pattern of altered expression of these cytokines and chemokines might then represent a signature of the mold exposure. To test this hypothesis, we isolated PBMCs from mold-exposed patients and unexposed controls. We then exposed those PBMCs to mold or mycotoxins ex vivo and compared their immunologic and inflammatory responses. We also subdivided the mold-exposed patients into asthmatics and nonasthmatics to explore whether any subset of cytokines and chemokines might distinguish them.

Materials and Methods

We investigated (1) if PBMCs respond to mold and/or mycotoxin challenges by altering cytokine and chemokine production, and (2) if the magnitude of any response was affected by previous chronic exposures to mold.

Ethics Statement

Informed signed written consent for blood collection and testing was obtained from all participants using a form and procedures approved by the Institutional Review Board Services, Aurora, Ontario, Canada (FDA/OHRP IORG Registration # IORG0000456). The Harvard T.H. Chan School of Public Health Office of Human Research Administration approved the ex vivo mycotoxin exposures and associated data collection.

Human Mold Exposures

Industrial hygiene consultants (Public Health and Safety, Irvine, CA) tested two office buildings with water intrusions in the Los Angeles metro area. They reported detectable levels of culturable molds. The industrial hygiene testing was approved by and performed according to the guidelines of the State of California Department of Workers’ Compensation. Air, surface and material samples collected from 37 sites throughout the two affected buildings, including multiple portions of their ventilation systems, were tested for viable mold. The fungal colonies were cultured and identified by microscopic and macroscopic morphology.

Patient and Control Selection

All 33 patients in this study worked in one of the affected office settings. The primary jobsites of mold-exposed patients were spread throughout different parts of these two buildings. All of the mold-exposed patients had occupational exposures from 1–5 years in length to S. chartarum, A. niger, C. herbarum, and P. chrysogenum. None reported experiencing water intrusions and/or mold contamination in any other personal environments, such as their residences. The mold-exposed patients did not complain of chronic respiratory symptoms until working in the contaminated buildings and had subsequently sought medical evaluations. The onset of symptoms varied from 1–3 months to over a year after starting in these job environments. Twenty-nine mold-exposed patients had blood samples collected while they were working in a contaminated building. Four mold-exposed patients had blood samples collected between 6 months and 2 years after transferring to a different work environment. We collected additional blood samples from 15 of the mold-exposed patients 18 months after the initial collection to assess the stability of responses. We do not know how often or how intensely the patients were exposed to the molds, spores, glucans, and mycotoxins detailed in the industrial hygiene reports.

The control groups were healthy adults from the Los Angeles metro area who did not work in the same buildings. The control subjects had no known history of exposures to molds at work or home and they lacked any related symptoms. None of the control subjects complained of chronic respiratory symptoms nor did they possess any medical history of chronic respiratory ailments. The controls for the mycotoxin experiments were on average slightly younger and had a different ethnic composition than the mold-exposed patients. The gender distribution was similar between the mold-exposed patients and controls.

These experiments included a total of 33 mold-exposed patients. Cells from 27 mold-exposed patients and 17 unexposed control subjects were challenged with S. chartarum mycotoxins ex vivo. Cells from five of these mold-exposed patients, together with six additional mold-exposed patients and eleven additional unexposed controls, comprising a “mold spore group,” were also challenged with spores from three additional molds. Demographic information, including age, gender, ethnicity, smoking status, pulmonary function and BMI, was collected from the mold-exposed patients.

Mold Spores and Mycotoxins

For immunologic challenges of cultured PBMCs with molds, the following mold strains were purchased from the American Tissue Culture Collection (Manassas, VA): S. chartarum, ATCC201212; A. niger, ATCC16888; C. herbarum, ATCC11281; and P. chrysogenum, ATCC10106.

The S. chartarum strain used was an ATCC clone originally isolated from an affected home during a Cleveland outbreak of infant idiopathic pulmonary hemorrhage [41]. This strain was chosen because it is well characterized and expresses a range of mycotoxins: satratoxin G (SG), satratoxin H, stachylysin, trichoverrol B, and roridin L-2 [41].

S. chartarum spores were grown on potato dextrose agar (PDA) plates at 15°C. After approximately 35 days of growth, spores were dried for 21 days and vacuumed using a modified filter cassette with a 37-mm, 0.4-μm polycarbonate membrane filter (Poretics Corp., Livermore, CA). These methods are modified from Rao et al. [42]. Dried spores were stored in glass at room temperature. We prepared Stachybotrys Spore Toxin (SST) from the S. chartarum spores, as previously described [27]. SST is a mixture of mycotoxins and other water-soluble extracts [27]. SG was obtained from Dr. James Pestka (Michigan State University, East Lansing, MI) [43]. Both toxin preparations were suspended in D-PBS and sterile filtered at 0.1 μm (Acrodisc syringe filters, Pall Life Sciences, Port Washington, NY). SG was diluted 1:1000 in media to a final concentration of 14.3 ng/ml. SST concentrations were standardized among batches based on SG content to a final SG concentration of 26.4 ng/ml in media. These toxin concentrations were determined in a small dose-response pilot study to cause <1% cytotoxicity.

Spore samples of A. niger, C. herbarum, and P. chrysogenum were suspended to a total volume of 1.4 ml in phosphate buffered saline (PBS) to final concentrations of 800 and 4,000 spores/ml (A. niger), 2,000 and 10,000 spores/ml (C. herbarum), and 4,000 and 20,000 spores/ml (P. chrysogenum). Spore concentrations were determined by counting spores in a counting chamber (Hausser Scientific, Horsham, PA).

Collection of Human Blood Samples

Twenty-eight ml of blood were obtained by venipuncture from mold-exposed patients and unexposed controls after obtaining their consent for research. Blood was withdrawn in four 7 ml tubes containing 0.081 ml of 15% K3 EDTA solution (BD Vacutainer, BD, Franklin Lakes, NJ). The samples were coded and shipped to the processing lab in insulated containers at room temperature.

Pulmonary Function Testing

All mold-exposed patients were requested to undergo methacholine challenges according to a protocol adapted from American Thoracic Society recommendations [44]. All but four of the 33 mold-exposed patients agreed to methacholine challenges. Briefly, prior to administration of methacholine, each patient was first challenged with NaCl. If this dropped their FEV1 by 10% from their baseline, the patient response was deemed positive. Mold-exposed patients without a positive response to NaCl were given methacholine at increasing doses and a positive response was defined as an FEV1 decrease of 20% or more from baseline following any dose of methacholine up to 16 mg/ml.

Isolation and Culture of PBMCs

PBMCs were isolated by Ficoll gradient centrifugation, as previously described [45, 46]. Cells were suspended at 106 cells/ml in RPMI 1640 medium supplemented with 1% penicillin-streptomycin, 1% L-glutamine and 10% fetal bovine serum (Invitrogen, Carlsbad, CA).

Ex Vivo Mold and Mycotoxin Exposures of PBMCs

We isolated PBMCs, challenged them ex vivo, and measured the cytokines and chemokines that were produced by the cells [46, 47].

We conducted two studies: a “mycotoxin study” where the challenges were either a mycotoxin (SG), an aqueous S. chartarum extract with a mix of mycotoxins and other antigens (SST), a non-specific toxin (ionomycin [positive control, VWR, Radnor, PA]), or media alone (negative control), each with or without the adjuvant phorbol 12-myristate 13-acetate (PMA, VWR, Radnor, PA); and a “mold spore study” where the challenges were with intact mold spores, a non-specific toxin (phytohemagglutinin [PHA]), or media alone.

S. chartarum mycotoxins and aqueous extracts.

Isolated PBMCs were suspended in RPMI at a concentration of 106 cells/ml and mixed with the designated treatment in a microcentrifuge tube, then plated in 96-well plates at a density of 200,000 PBMCs per well. Isolated PBMCs from each patient and control were exposed to SST, SG, ionomycin (1μg/ml), or media, in the presence or absence of PMA (1:1000). The mixture was plated in triplicate 200 μl aliquots in 96-well plates, and then incubated at 37°C with 5% CO2. We collected cell culture supernatants 18 hours after the exposures began by centrifuging the 96-well plates and transferring supernatants to new 96-well plates, which were then frozen at -80°C.

A niger, C. herbarum, and P. chrysogenum mold spores.

PBMCs were cultured in 6-well plates at 106 cells/ml, 3 ml/well. Mold spores were added at two concentrations each: A. niger, 800 and 4,000 spores/ml; C. herbarum, 2,000 and 10,000 spores/ml; and P. chrysogenum, 4,000 and 20,000 spores/ml. These concentrations were five and twenty-five times higher, respectively, than the lowest observed effect level (LOEL): mold spore concentrations that activated extracellular expression of cytokines and chemokines by a statistically significant value (not shown). We added 10 μg/ml phytohemagglutinin (PHA-P, Sigma-Aldrich, St. Louis, MO) to separate plates as a positive control. PHA has high mitogenic activity and induces PBMC proliferation and cytokine secretion [48, 49]. Three replicate plates were used for each condition. Negative control samples contained media only and were used for determining the basal levels of cytokine and chemokine expression.

Cytokine and Chemokine Assays

At two points during the study, cell culture supernatants were sent as a batch to Eve Technologies (Calgary, Alberta, Canada) for the measurement of 41 cytokines and chemokines (EGF, eotaxin-1, FGF-2, Flt-3L, fractalkine, G-CSF, GM-CSF, GRO(pan), IFNα2, IFN-γ, IL-1α, IL-1β, IL-1ra, IL-2, IL-3, IL-4, IL-5, IL-6, IL-7, IL-8 (CXCL8), IL-9, IL-10, IL-12 p40, IL-12 p70, IL-13, IL-15, IL-17A, IP-10 (CXCL10), MCP-1, MCP-3, MDC, MIP-1α (CCL3), MIP-1β (CCL4), PDGF-AA, PDGF-AB/BB, RANTES, sCD40L, TGF-α, TNF-α, TNF-β, and VEGF-A). The first batch used Milliplex polystyrene beads (Millipore, Billerica, MA) on a Luminex platform. The second batch used a Milliplex magnetic bead panel, also manufactured by Millipore, on a Luminex platform. The change in platform was necessitated by a change in technology and kit availability. Since some of the magnetic beads were more likely to aggregate in the presence of cellular debris, the cell culture supernatants were filtered before analysis. In our statistical analyses, we excluded cytokines and chemokines that were not consistent between polystyrene and magnetic beads in reference samples (p-values < 0.01). We also excluded from the analyses cytokines and chemokines that were not stable over time based on the subset of the nine mold-exposed patients whose repeat blood samples were collected 18 months apart.

Cytokine and chemokine concentrations from cell-culture supernatants from PBMCs exposed to A. niger, C. herbarum, and P. chrysogenum were measured by Luminex xMAP bead array immunoassays using a Bio-Plex 200 fluorescence bead reader (BioRad Laboratories, Hercules, CA). Three panels of antibody-conjugated beads for measuring human inflammatory cytokines (GM-CSF, IL-1β, IL-6, IL-8, TNF-α), Th1/Th2 cytokines (IFN-γ, IL-2, IL-4, IL-5, IL-10) and chemokines (MIP-1α, MIP-1β, MCP-1, eotaxin, RANTES) (BioSource, Camarillo, CA) were used in the assay according to the manufacturers’ instructions.

Statistical Methods

We performed both single and multiple-marker analyses to identify cytokine and chemokine response signatures for each mold and mycotoxin exposure.

Single-marker analysis.

For each mold challenge, to examine associations of cytokine and chemokine response levels with case-control status, we performed single-marker analyses using multiple logistic regressions. We log-transformed cytokine response levels in the model. We adjusted for age, gender, ethnicity and smoking status in the models for mycotoxin challenges. We did not adjust for any covariates in the models for whole mold spore challenges due to missing demographic information for some of the subjects. To correct for multiple testing, we estimated false discovery rates (FDR) and defined statistical significance as FDR q-values <0.05, which indicated that no more than 5% of the significant results by p-value <0.05 were considered to be false positive findings [50, 51]. FDR is a common approach to correct for multiple testing in high-throughput screening data, such as microarrays, proteomics and large-scale genotyping studies. In mold-exposed patients, we also explored associations between lung function (methacholine responses) and biomarker (cytokine and chemokine) response levels stratified by mold challenges, as well as associations between clinical symptoms and biomarkers (cytokine and chemokine response levels) by multiple logistic regression models adjusting for age, gender, ethnicity, smoking status and BMI. All of the above analyses were done using R statistical software.

Multiple-marker analysis to identify cytokine and chemokine signatures.

Multiple-marker analyses were conducted using the support vector machine (SVM) algorithm. Analyzing multiple biomarkers together generally increases the accuracy of disease risk prediction [52]. The SVM approach is a widely used machine learning algorithm for disease prediction and classification with high dimension data [53]. The SVM approach has been used to identify gene expression patterns [54] and biomarker signatures [55] of cancers[56], diabetes [57], and major adverse cardiac events [58], as well as a recent application to classify liver cirrhosis patients from healthy controls by gut microbiota metagenomics [59]. One of the major advantages of applying SVM is that it is an effective classifier without local minima issues that only depends on parsimonious parameters. In addition to linear classification (such as the commonly used logistic regression model), SVM can efficiently perform non-linear statistical classification, which is an important feature when the relationship between predictors and outcomes are not linear.

To identify the cytokine and chemokine response signatures of mold-exposed patients, we used the heuristic search-SVM [60–63] approach. The model building for the heuristic search-SVM involves four major steps: (1) imputation of missing data, (2) feature (biomarker) selection, (3) classification and (4) cross-validation. The Area Under the Receiver Operating Characteristic (ROC) Curve (AUC) [64] was estimated in the SVM analysis and used to quantify the overall ability of the selected biomarkers to discriminate disease patients from controls. In our study, the AUC was used as the classification performance index to find the best signature (biomarker set) of cytokine and chemokine response levels to differentiate mold-exposed patients from unexposed controls. The SVM analysis was performed using the R svmpath (version 0.953) [65] and the LibSVM (version 3.18) packages [66].

Step 1: Imputation of missing data.

To minimize the influence of missing biomarkers data in a few samples, we applied a k-Nearest Neighbor (kNN) approach [67] to estimate the missing values for these cytokine response measurements. The selection of the k-nearest cytokines was based on the comparison of cytokine response levels with those of the cytokines of interest available in other samples [67]. In addition, case-control status and mold toxicity challenge information were included in the imputation model to avoid any attenuation of the estimated effects in later analyses [68].

Step 2 and Step 3: Feature selection and classification.

The feature selection and classification steps were used to identify the optimal combination of cytokine and chemokine responses that can best distinguish cases from controls. This combination of biomarkers can be considered a unique signature of cytokine and chemokine responses in mold-exposed cases. The input features in the SVM model included cytokine responses, chemokine responses, age, gender, ethnicity, and cigarette smoking. Body weight and pulmonary function were not available for unexposed controls and so were not included in the analyses. The feature selection was done by heuristic search [60–63] and performed with a simulated annealing algorithm [69]. This heuristic search algorithm was embedded within the optimization procedure of classification and cross-validation (described below), which increases prediction accuracy [61, 70]. For classification, we used a non-linear SVM classifier with the Gaussian radial basis function kernel. The non-linear SVM classifier mapped all data points into a higher dimensional feature space with non-linear transformation.

Step 4: Cross-validation.

To reduce the risk of over-fitting and to estimate AUC accuracy, we performed a 10-fold cross-validation [70] 100 times. Each 10-fold cross-validation had a different random seed in order to randomly partition the overall sample into 10 subsets. Since cross-validation produced a potentially different model for each subset of the data, the classification using all observations (i.e., without cross-validation) was displayed for purposes of describing the optimal model.

Results

Patient and Control Demographics

The demographics of the mold-exposed patients and unexposed controls for our study are shown in Table 1. The patients were 49.9 ± 8.4 years old (mean ± standard deviation) and the controls were 40.9 ± 15.9 years old. Fifty-two percent of the patients were female and 48% were male, while 29% of the unexposed controls were female and 71% were male. Two patients and four controls were current smokers. Of the patients, 70.4% were Caucasian, 14.8% were Hispanic, 11.1% were Asian, and 3.7% were African-American. Of the controls, 52.9% were Caucasian and 47.1% were African-American.

[Figure omitted. See PDF.]

Table 1. Characteristics of subjects who provided blood for ex vivo mycotoxin exposures.

https://doi.org/10.1371/journal.pone.0126926.t001

The patients in the mold spore study ranged from 37–64 years old. Sixty-four percent of the patients were female and 36% were male. One patient was a current smoker and another was a former smoker. Of the patients, 27% were Caucasian, 36% were Hispanic, and 36% were African American.

Industrial hygiene.

The industrial hygiene consultants tested multiple places within the workplaces including individual offices and commonly used spaces, such as the lobby, front desk area and lounge. They reported measurable levels of Cladosporium, Stachybotrys, A. niger, and P. chrysogenum in multiple places in the building. The full industrial hygiene report is available in S1 Table.

Mycotoxins Alter Expression of Cytokines and Chemokines by PBMCs

Single-marker analysis.

After we challenged PBMCs from mold-exposed patients and controls with media (negative control), SST, SG, or ionomycin, with or without PMA, we identified the differences in cytokine and chemokine responses between these groups by running a single-marker analysis that controlled for age, gender, ethnicity, and smoking status. The majority of the cytokine and chemokine responses with q<0.05 (false discovery rate q-value correcting for multiple testing, q<0.05 corresponding to p<0.021) were observed in response to SG or SST alone. When PMA was used as an adjuvant with SST, SG or ionomycin, it blunted the differences between mold-exposed patients and controls, and PMA alone caused marked cytokine and chemokine responses. The following cytokines and chemokines were significantly higher in response to both SG and SST in mold-exposed patients compared to controls: eotaxin, INF-α, IL-1α, IL-12 p40, IL-12 p70, IP-10, PDGF-AA, TNF-β, and VEGF (q<0.05 in Table 2). Notably, none of these cytokines or chemokines were significantly higher in cases when PBMCs were exposed to media alone.

[Figure omitted. See PDF.]

Table 2. Cytokine and chemokine concentrations in PBMC cultures exposed to Satratoxin G (SG), Stachybotrys Spore Toxin (SST), or Ionomycin ex vivo.

https://doi.org/10.1371/journal.pone.0126926.t002

Multiple-marker analysis.

To identify the cytokine and chemokine response signature in PBMCs of mold-exposed patients after different mycotoxin challenges, we applied the support vector machine (SVM) approach. We found significant response signatures in patient cells challenged with SG, SST and ionomycin (without PMA) (Fig 1).

[Figure omitted. See PDF.]

Fig 1. Cytokine and chemokine response signatures and cumulative AUCs from multiple-marker analysis.

The estimated cumulative AUCs from SVM analysis are shown to demonstrate the prediction probability of separating mold-exposed patients from controls for combinations of cytokine and chemokine measurements after challenges. The first cytokine or chemokine shown in each graph is the best cytokine or chemokine (among the tested cytokines and chemokines) to differentiate mold-exposed patients from controls after challenge with the listed cellular exposure. Each subsequent cytokine or chemokine adds additional separation (AUC), as depicted on the graph. Each graph shows the top 10 cytokines and chemokines and the resulting AUC when each additional cytokine or chemokine is added to the model. Note that specific cytokines and chemokines on the x-axis generate a toxin-specific signature. A. Satratoxin G (SG), B. SG + phorbol 12-myristate 13-acetate (PMA), C. S. chartarum spore [mixed] toxins (SST), D. SST + PMA, E. Ionomycin, F. Ionomycin + PMA, G. media (negative control), H. media + PMA.

https://doi.org/10.1371/journal.pone.0126926.g001

The area under a receiver operating characteristic (ROC) curve (AUC) ranges between 0.5 and 1.0. A truly insignificant biomarker (one no better at identifying true positives than flipping a coin) has an AUC of 0.5. A perfect set of biomarkers (zero false positives and zero false negatives) has an AUC of 1.00. The higher the AUC, the more accurately the biomarkers have identified (classified) subjects as disease patients versus non-disease controls. These cytokine and chemokine response signatures accurately classified mold-exposed patients from unexposed controls with AUC > 0.97 after controlling for age, gender, ethnicity and smoking status. The cytokine and chemokine response signature of mold-exposed patients was different between SG-exposed PBMCs and SST-exposed PBMCs. The cytokine and chemokine response signature of SG-exposed PBMCs separating mold-exposed patients from unexposed controls comprised six cytokines and chemokines: IL-12 p40, RANTES, IL-1α, MCP-3, TNF-α, and PDGF-AA. The cytokine and chemokine response signature of SST-exposed PBMCs from the mold-exposed group comprised four cytokines and chemokines: IL-12 p40, MIP-1β, G-CSF, and MCP-1. The full cytokine and chemokine data set is available in S2 Table.

Mold Spores Alter Expression of Cytokines and Chemokines by PBMCs

To evaluate cellular responses to whole mold spores, we challenged PBMCs from mold-exposed patients and unexposed controls with mold spores ex vivo and measured the expression of 15 cytokines and chemokines (GM-CSF, IFN-γ, IL1-β, IL-2, IL-4, IL-5, IL-6, IL-8, IL-10, TNF-α, MIP-1β, MCP-1, eotaxin, MIP-1α and RANTES) in cell culture supernatants. Different mold species evoked distinct cytokine and chemokine responses. Seven cytokines and chemokines had significantly different expression profiles between mold-exposed patients and controls in at least one mold challenge (eotaxin, IL1-β, IL-6, IL-8, MCP-1, MIP-1α, and MIP-1β). For ex vivo challenges with P. chrysogenum and C. herbarum, we observed a significant decrease in multiple cytokine and chemokine responses in mold-exposed patients compared to controls that were not evident with media alone (Table 3). For example, a decrease in the production of eotaxin, IL-8, MCP-1, MIP-1α, and MIP-1β was observed in patient samples challenged with P. chrysogenum (q<0.05). Lower production of IL-1β, IL-6, MCP-1, and MIP-1α was observed in patient samples challenged with C. herbarum (q<0.05). In patient samples challenged with A. niger, there were no significant changes in cytokines or chemokines that were not evident in cells incubated in media alone (q<0.05). Expression levels varied for different cytokines and chemokines but concentrations of most cytokines and chemokines were increased at higher mold doses.

[Figure omitted. See PDF.]

Table 3. Cytokine and chemokine concentrations in PBMC cultures exposed to P. chrysogenum, C. herbarum, or media alone ex vivo.

https://doi.org/10.1371/journal.pone.0126926.t003

To determine if differing responses to mold spores in mold-exposed patients and unexposed controls were mold-specific or reflected a change in responsiveness of patient cells, we compared responses to PHA and media alone from mold-exposed patients and controls. There were no significant differences between mold-exposed patients and controls in cytokine or chemokine responses to the PHA challenge (FDR q<0.05, not shown). PBMCs from mold-exposed patients and unexposed controls produced different amounts of eotaxin, MIP-1α, and RANTES in media alone. The magnitude of these differences was more pronounced after ex vivo challenges with C. herbarum and P. chrysogenum, but not with A. niger (Table 3). The full cytokine and chemokine data set for the mold spore exposures is available in S3 Table.

Stratification of Cases

In addition to comparing cases to controls, we took advantage of additional data and stratified cytokine and chemokine responses by methacholine responsiveness and associated symptoms among the mold cases.

Asthmatics had a mold-independent cytokine and chemokine response signature.

We evaluated the correlation between cytokine and chemokine responses and methacholine test results. In the regression model we used, we adjusted for BMI, age, gender, ethnicity, and smoking status. This adjustment is important because BMI (p = 0.004) and age (p = 0.003) were strongly associated with methacholine results and may directly affect sensitivity to methacholine in a cytokine- and chemokine-independent manner. Though SG and SST alone were the best separators of mold-exposed patients and controls in the single-marker analysis and were excellent separators of mold-exposed patients and controls in the multiple-marker analysis, ionomycin was a better separator of mold-exposed patients with asthma (positive methacholine) from mold-exposed patients without asthma (negative methacholine). We found that MIP-1β, IL-17, TGF-α, and IL-10 were all significantly higher in mold-exposed asthmatics than mold-exposed nonasthmatics (p = 0.007, p = 0.0026, p = 0.0056, and p = 0.0065, respectively; FDR q<0.05) in ionomycin-exposed PBMCs. These four cytokines were likewise significantly higher in mold-exposed nonasthmatics compared to unexposed control subjects (p = 4x10-5, p = 6x10-4, p = 8x10-8, and p = 1x10-5, respectively; FDR q<0.05). There were no significant differences in cytokine or chemokine production between mold-exposed nonasthmatics and asthmatics when their PBMCs were exposed to media, SG, SST, A. niger, C. herbarum or P. chrysogenum. This asthma cytokine and chemokine signature was mold-species independent since these differences only became apparent when their cells were exposed to ionomycin ex vivo.

Cytokine and chemokine responses correlate with clinical symptoms in mold-exposed patients.

The mold-exposed patients reported many symptoms consistent with respiratory allergies and asthma. Every patient reported at least one respiratory symptom. More than ⅔ of the mold-exposed patients complained of shortness of breath, coughing, sneezing, runny nose, or nasal congestion. Approximately ½ reported eye itching, burning or tearing and ½ reported wheezing. About ⅓ reported headaches, throat irritation, chronic phlegm production or sinus congestion.

Several cytokine and chemokine responses in SG- and SST-exposed PBMCs were associated with self-reported patient symptoms. Lower levels of CD40L, PDGF-AA, and RANTES after SG challenge were associated with a self-report of head, eye, ear, nose and throat symptoms including itchy, burning or watery eyes; stuffy or runny nose; sneezing; and headaches (p = 0.0003, p = 0.0009, and p = 0.004, respectively; FDR q<0.05). A higher level of human macrophage-derived cytokine (MDC) after SST challenge was also associated with a self-report of these symptoms (p = 0.002, FDR q<0.05). Similarly, a self-report of weight gain is highly correlated with several cytokines and chemokines in response to SST, including INF-α, MDC, and RANTES (p = 0.0002, p = 0.003, and p = 0.01, respectively; FDR q<0.05).

Because 100% of mold-exposed patients reported respiratory symptoms, no correlation with any specific cytokine or chemokine levels in the patient group can be made.

Discussion

To identify abnormalities associated with indoor environmental mold exposures, we challenged PBMCs from 33 patients who were chronically exposed to molds at their jobsites. When their PBMCs were exposed to mycotoxins, mold extracts, and/or molds ex vivo, we found significant differences in the chemokines and cytokines produced. Expression of these proteins differed between the mold-exposed patients and unexposed controls. All the mold-exposed patients reported respiratory symptoms. A third of the mold-exposed patients had abnormal methacholine responsiveness indicative of asthma.

Cytokine and Chemokine Responses to Mold and Mycotoxin Exposures of PBMCs

From our single-marker analysis, we identified unique cytokine and chemokine response signatures for a history of mold exposure and asthma status. In addition, we used an SVM multiple-marker analysis to identify several mold-specific cytokine and chemokine response signatures in mold-exposed patients. SVM is a popular machine learning algorithm that relies on multiple artificial intelligence algorithms to recognize the patterns/trends in experimental data. It uses recognized patterns to classify and predict in a wide range of scientific applications [53], such as gene expression profiling [54] and biomarker signatures [55] of patients with cancers [71, 72], diabetes [57], or major adverse cardiac events (MACE) [58]. It is an effective classifier without local minima issues and depends only on parsimonious parameters. In addition to linear classification (such as the commonly used logistic regression model), SVM can efficiently perform non-linear statistical classifications, an important feature where the relations between predictors and outcomes are not linear.

The multiple-marker SVM analysis gives an AUC of 0.98 for SG and 0.97 for SST for a particular group of cytokines and chemokines. We identified somewhat different sets of cytokines and chemokines using multiple-marker and the single-marker analyses. Two cytokines or chemokines with significant p-values in the single-marker analysis might be highly correlated. In the multiple-marker analysis, adding a second highly correlated cytokine or chemokine might not add additional separation between the groups versus one of those cytokines or chemokines alone. In the online supplement (S1 Text), we highlight the functions and underlying pathways of the cytokines and chemokines that were identified as potential mediators of inflammatory responses to both SG and SST.

Although each cytokine or chemokine may influence the phenotypes observed, we point out that cytokines and chemokines work together in an orchestrated fashion. It is unlikely that a single cytokine or chemokine will be a perfect biomarker for separating mold-exposed patients and controls. Instead, a pattern of cytokine and chemokine expression must be identified. Quantifying multiple biomarkers increases the accuracy of disease correlations [52]. Most common diseases are caused by the deregulation of complex processes in multiple pathways. Aggregating biomarkers reduces noise and provides more power to detect deregulation of complete functional units, and hence it produces a clearer picture of the underlying pathophysiological process.

Molds may induce persistent changes in inflammatory and immune responses. Epigenetic mechanisms can mediate changes in cytokine and chemokine gene expression and sensitization after chronic mold exposure. Regulatory regions of DNA could be methylated or histones acetylated or deacetylated in response to mold or mycotoxin exposures [73]. This could alter the expression of cytokines and chemokines in response to subsequent stimulation. In addition, memory T cells might activate signaling pathways that trigger cytokine or chemokine release in the presence of mold or mycotoxin challenges.

Mold and Asthma

We stratified mold-exposed patients by self-reported symptoms as well as by a more objective measure of asthma, methacholine responsiveness. Cytokines and chemokines identified in the initial case-control analysis and the subsequent stratification analysis are more likely to be important mediators of long-term sensitivity to molds. Higher BMI was also correlated with methacholine responsiveness. This is consistent with previously reported correlations between obesity and asthma [74].

Among the mold-exposed patients, only ionomycin challenge differentiated asthmatics from nonasthmatics. Ionomycin challenges of PBMCs may reveal cytokine and chemokine response signatures that differentiate asthmatics from nonasthmatics among mold-exposed patients who have respiratory symptoms. The asthma cytokine and chemokine response signature from ionomycin is distinct from the cytokine and chemokine response signatures from the molds and mycotoxins we tested. Since ionomycin challenges are not specific to mold exposures, ionomycin is a strong candidate for objectively separating asthmatics from nonasthmatics even in a population with other respiratory symptoms. These results should be validated in a larger clinical study.

Our data are novel in that they suggest alterations in immune and inflammatory system pathways in response to molds. The ex vivo ionomycin challenges identified four cytokines, IL-17, IL-10, TGF-α, and MIP-1β, that are elevated in mold-exposed asthmatics compared to mold-exposed nonasthmatics and unexposed controls. IL-17 is a lymphokine classically associated with asthma. Th17 cells secrete both it and IL-10 to balance protective and pathologic T-cell responses to multiple pathogens, including molds [75]. Increased levels of IL-10 limit the inflammatory response to pathogens. IL-10 may play a role in down-regulating some proinflammatory cytokines in asthmatic patients with chronic mold exposure [76]. TGF-α is a member of the epidermal growth factor family. It plays a critical role in airway remodeling and mucus production in asthma [77]. TGF-α may exacerbate the respiratory symptoms in asthmatics with chronic mold exposure. MIP-1β (CCL4) is a neutrophil chemoattractant that contributes to both early and late phase asthma responses [78].

The asthmatic response to ionomycin could be the result of differences in sub-populations of PBMCs in asthmatics and nonasthmatics or might be the result of epigenetic changes during asthma development in genes regulating cytokine responses. This result should be validated in a larger population and further mechanistic studies would be required to clarify why the asthmatic patients are more responsive to ionomycin.

Moldy environments may contribute to upper respiratory symptoms even in the absence of IgE-mediated allergic sensitivity [79]. Our results support a link between respiratory illness and human cellular immune and inflammatory responses that have been previously reported [80, 81]. When PBMCs were challenged with mycotoxins or molds, the different patterns of cytokine and chemokine expression observed in mold-exposed patients versus controls was consistent with an activation of inflammatory cells and an increase in cytokine and chemokine production both in vivo [28, 82–85] and in cell cultures [37–39, 86–88].

We observed different cytokine and chemokine response signatures after challenges with different mycotoxins or mold species. This specificity may explain differences in cytokine and chemokine expression patterns in other studies [89]. These different signatures support our view that chronic exposures to mold induce altered immune system responses to subsequent mold challenges in a mycotoxin- and mold-specific manner. Because similar cytokine and chemokine responses could be caused by other inflammatory or immune responsive agents such as other allergens or bacteria [90] and we don’t have industrial hygiene reports from the buildings where unexposed control subjects work, a larger clinical study that better characterizes the environment where control subjects work will be required to fully validate a cellular ex vivo test that can be used to characterize human health effects from exposures to individual molds. In addition, a future clinical study should include other common respiratory exposures that can lead to cytokine and chemokine release to establish whether these cytokine and chemokine signatures are unique to specific molds.

We demonstrated for the first time that chronic exposures to molds alter cytokine and chemokine responses of PBMCs to specific mold or mycotoxin challenges. Inhibited or enhanced production of cytokines in mold-exposed patients may reflect elements of both sensitization and tolerance. Inhibition of production of some cytokines and chemokines in mold-exposed patients may be analogous to the suppression of inflammatory responses observed in PBMCs from allergen-sensitized subjects [91]. Mold toxins may suppress the immune system through a balance of cytotoxicity and altered Th1/Th2 balance. The alteration of immune responses due to chronic mold exposures may also adversely affect the ability of the immune system to fight infections and other environmental challenges. This may explain patient complaints of concurrent susceptibility to infectious organisms and enhanced responses to chemical irritants. Some changes we observed with ex vivo exposure to S. chartarum mycotoxins were consistent with sensitization. Prior exposures can lead to increases in subsequent inflammatory responses [92].

Future Studies

Our data are consistent with studies that emphasize the prevention of mold exposures as one way to reduce the incidence and severity of asthma [93]. Since many molds have an adverse impact on respiratory health, it is vital to reduce mold contamination in indoor environments. Especially when there is water intrusion-related mold growth, we should aggressively remediate these contaminated environments as soon as possible to minimize associated health risks. Given the ubiquitous nature of mold spores and their ability to remain dormant and survive in isolated spaces, remediation and renovation of some contaminated buildings may be inadequate and total replacement may be required. Studies of the effectiveness of various remediation strategies are needed.

Mold-induced changes in chemokine and cytokine release may suggest a review of best practices for the diagnosis and treatment of asthma and other respiratory illnesses. For example, current management of patients with mold-related illnesses often involves oral, parenteral or inhalation administration of corticosteroid-based medications that act by suppressing the immune system. However, the underlying illness of asthma, when induced or aggravated by molds, is accompanied by both increases and decreases in immune and inflammatory pathways, as shown by our findings. Perhaps a goal of treatment should include the restoration of normal immune system functioning, in which case steroid-based therapies should be reassessed.

Four of the mold-exposed patients in this study transferred to uncontaminated worksites six months or more before their blood samples were collected. It is important to document to what extent mold-related immune system responses persist after mold exposures cease. This phenomenon may serve as a basic etiologic mechanism in the development of chronic respiratory illnesses such as asthma. Future studies should include a larger number of subjects who have histories of mold exposures but who are no longer being actively exposed.

The results of this study provide guidance for an expanded clinical study to examine these findings in larger population groups. Such a study should measure mold-specific IgE and pulmonary function in all mold-exposed patients and healthy controls. In parallel, a full industrial hygiene assessment is necessary for all patient and control workplaces and homes since we know that mold exposures are not all or nothing, but rather a continuum. A future study might also be designed to assess the contribution of mold exposures to the initiation and aggravation of pediatric and adult asthma.

Conclusions

We demonstrate for the first time that chronic exposures to environmental molds are associated with changes in cytokine and chemokine production by isolated PBMCs in response to specific mold or mycotoxin challenges. Ex vivo exposures of PBMCs to molds or mycotoxins can differentiate mold-exposed patients from unexposed controls. This strategy may be applicable to other inflammation-related environmental toxins.

Supporting Information

[Figure omitted. See PDF.]

S1 Table. Full Industrial Hygiene Report.

https://doi.org/10.1371/journal.pone.0126926.s001

(PDF)

S2 Table. Full cytokine and chemokine data set for ex vivo mycotoxin exposures.

https://doi.org/10.1371/journal.pone.0126926.s002

(XLSX)

S3 Table. Full cytokine and chemokine data set for ex vivo mold spore exposures.

https://doi.org/10.1371/journal.pone.0126926.s003

(XLSX)

S1 Text. Detailed descriptions of the possible cytokine and chemokine pathways altered by S. chartarum exposure.

https://doi.org/10.1371/journal.pone.0126926.s004

(DOCX)

Acknowledgments

We thank Dr. James Pestka at Michigan State University for providing us with purified satratoxin G, Dr. Zarema Arbieva for her help with statistical analyses, and Ms. Melissa Curran for her editorial assistance.

Author Contributions

Conceived and designed the experiments: JHRL YH IMG TCD RMM KJT BSG JDB. Performed the experiments: JHRL IMG TCD KJT BSG. Analyzed the data: JHRL YH IMG TCD RMM KJT CC JDB. Contributed reagents/materials/analysis tools: JHRL YH IMG TCD RMM KJT CC BSG JDB. Wrote the paper: JHRL YH IMG TCD RMM KJT CC BSG JDB. Responsible for mold-exposed patient and unexposed control recruitment: BSG.

Citation: Rosenblum Lichtenstein JH, Hsu Y-H, Gavin IM, Donaghey TC, Molina RM, Thompson KJ, et al. (2015) Environmental Mold and Mycotoxin Exposures Elicit Specific Cytokine and Chemokine Responses. PLoS ONE 10(5): e0126926. https://doi.org/10.1371/journal.pone.0126926

References

1. 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 Perspect. 2011;119(6):748–56. pmid:21269928; PubMed Central PMCID: PMC3114807.

2. Clark N, Ammann H, Brunekreef B, Eggleston P, Fisk W, Fullilove R, et al. Damp indoor spaces and health. Washington, DC: Institute of Medicine of the National Academies; 2004.

3. Flamant-Hulin M, Annesi-Maesano I, Caillaud D. Relationships between molds and asthma suggesting non-allergic mechanisms. A rural-urban comparison. Pediatr Allergy Immunol. 2013;24(4):345–51. pmid:23692327.

4. Weinmayr G, Gehring U, Genuneit J, Buchele G, Kleiner A, Siebers R, et al. Dampness and moulds in relation to respiratory and allergic symptoms in children: results from Phase Two of the International Study of Asthma and Allergies in Childhood (ISAAC Phase Two). Clin Exp Allergy. 2013;43(7):762–74. pmid:23786283.

5. Zock JP, Jarvis D, Luczynska C, Sunyer J, Burney P. Housing characteristics, reported mold exposure, and asthma in the European Community Respiratory Health Survey. J Allergy Clin Immunol. 2002;110(2):285–92. pmid:12170270.

6. Blatter J, Forno E, Brehm J, Acosta-Perez E, Alvarez M, Colon-Semidey A, et al. Fungal exposure, atopy, and asthma exacerbations in Puerto Rican children. Ann Am Thorac Soc. 2014;11(6):925–32. pmid:24915164; PubMed Central PMCID: PMC4213992.

7. 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(4):459–67. pmid:9641573.

8. Purokivi MK, Hirvonen MR, Randell JT, Roponen MH, Meklin TM, Nevalainen AL, et al. Changes in pro-inflammatory cytokines in association with exposure to moisture-damaged building microbes. Eur Respir J. 2001;18(6):951–8. pmid:11829101.

9. Engman LH, Bornehag CG, Sundell J. How valid are parents' questionnaire responses regarding building characteristics, mouldy odour, and signs of moisture problems in Swedish homes? Scand J Public Health. 2007;35(2):125–32. pmid:17454915.

10. Dales RE, Miller D, McMullen E. Indoor air quality and health: validity and determinants of reported home dampness and moulds. Int J Epidemiol. 1997;26(1):120–5. pmid:9126511.

11. Tischer CG, Heinrich J. Exposure assessment of residential mould, fungi and microbial components in relation to children's health: achievements and challenges. Int J Hyg Environ Health. 2013;216(2):109–14. pmid:22704485.

12. Hyvarinen A, Husman T, Laitinen S, Meklin T, Taskinen T, Korppi M, et al. Microbial exposure and mold-specific serum IgG levels among children with respiratory symptoms in 2 school buildings. Arch Environ Health. 2003;58(5):275–83. pmid:14738273.

13. Immonen J, Laitinen S, Taskinen T, Pekkanen J, Nevalainen A, Korppi M. Mould-specific immunoglobulin G antibodies in students from moisture- and mould-damaged schools: a 3-year follow-up study. Pediatr Allergy Immunol. 2002;13(2):125–8. pmid:12000485.

14. Shelton BG, Kirkland KH, Flanders WD, Morris GK. Profiles of airborne fungi in buildings and outdoor environments in the United States. Appl Environ Microbiol. 2002;68(4):1743–53. pmid:11916692; PubMed Central PMCID: PMC123871.

15. Betancourt DA, Krebs K, Moore SA, Martin SM. Microbial volatile organic compound emissions from Stachybotrys chartarum growing on gypsum wallboard and ceiling tile. BMC Microbiol. 2013;13(1):283. pmid:24308451; PubMed Central PMCID: PMC4234204.

16. Barbeau DN, Grimsley LF, White LE, El-Dahr JM, Lichtveld M. Mold exposure and health effects following hurricanes Katrina and Rita. Annu Rev Public Health. 2010;31:165–78. pmid:20070193.

17. Grimsley LF, Chulada PC, Kennedy S, White L, Wildfire J, Cohn RD, et al. Indoor environmental exposures for children with asthma enrolled in the HEAL study, post-Katrina New Orleans. Environ Health Perspect. 2012;120(11):1600–6. pmid:22894816; PubMed Central PMCID: PMC3556612.

18. Johanning E, Biagini R, Hull D, Morey P, Jarvis B, Landsbergis P. Health and immunology study following exposure to toxigenic fungi (Stachybotrys chartarum) in a water-damaged office environment. Int Arch Occup Environ Health. 1996;68(4):207–18. Epub 1996/01/01. pmid:8738349.

19. Yike I, Distler AM, Ziady AG, Dearborn DG. Mycotoxin adducts on human serum albumin: biomarkers of exposure to Stachybotrys chartarum. Environ Health Perspect. 2006;114(8):1221–6. Epub 2006/08/03. pmid:16882529; PubMed Central PMCID: PMC1552036.

20. Al-Ahmad M, Manno M, Ng V, Ribeiro M, Liss GM, Tarlo SM. Symptoms after mould exposure including Stachybotrys chartarum, and comparison with darkroom disease. Allergy. 2010;65(2):245–55. Epub 2009/10/03. pmid:19796210.

21. Crook B, Burton NC. Indoor moulds, Sick Building Syndrome and building related illness. Fungal Biol Rev. 2010;24(3–4):106–13.

22. Dearborn DG, Yike I, Sorenson WG, Miller MJ, Etzel RA. Overview of investigations into pulmonary hemorrhage among infants in Cleveland, Ohio. Environ Health Perspect. 1999;107 Suppl 3:495–9. pmid:10346998; PubMed Central PMCID: PMC1566217.

23. Dearborn DG, Smith PG, Dahms BB, Allan TM, Sorenson WG, Montana E, et al. Clinical profile of 30 infants with acute pulmonary hemorrhage in Cleveland. Pediatrics. 2002;110(3):627–37. pmid:12205270.

24. Etzel RA, Montana E, Sorenson WG, Kullman GJ, Allan TM, Dearborn DG, et al. Acute pulmonary hemorrhage in infants associated with exposure to Stachybotrys atra and other fungi. Arch Pediatr Adolesc Med. 1998;152(8):757–62. Epub 1998/08/13. pmid:9701134.

25. Update: pulmonary hemorrhage/hemosiderosis among infants—Cleveland, Ohio, 1993–1996. MMWR Morb Mortal Wkly Rep. 1997;46(2):33–5. pmid:9011781.

26. Update: Pulmonary hemorrhage/hemosiderosis among infants—Cleveland, Ohio, 1993–1996. MMWR Morb Mortal Wkly Rep. 2000;49(9):180–4. pmid:11795499.

27. Lichtenstein JH, Molina RM, Donaghey TC, Amuzie CJ, Pestka JJ, Coull BA, et al. Pulmonary responses to Stachybotrys chartarum and its toxins: mouse strain affects clearance and macrophage cytotoxicity. Toxicol Sci. 2010;116(1):113–21. Epub 2010/04/14. pmid:20385656; PubMed Central PMCID: PMC2886860.

28. Rosenblum Lichtenstein JH, Molina RM, Donaghey TC, Brain JD. Strain differences influence murine pulmonary responses to Stachybotrys chartarum. Am J Respir Cell Mol Biol. 2006;35(4):415–23. Epub 2006/05/13. pmid:16690987; PubMed Central PMCID: PMC2643262.

29. Leino M, Makela M, Reijula K, Haahtela T, Mussalo-Rauhamaa H, Tuomi T, et al. Intranasal exposure to a damp building mould, Stachybotrys chartarum, induces lung inflammation in mice by satratoxin-independent mechanisms. Clin Exp Allergy. 2003;33(11):1603–10. Epub 2003/11/18. pmid:14616875.

30. Flemming J, Hudson B, Rand TG. Comparison of inflammatory and cytotoxic lung responses in mice after intratracheal exposure to spores of two different Stachybotrys chartarum strains. Toxicol Sci. 2004;78(2):267–75. pmid:14718650.

31. Viana ME, Coates NH, Gavett SH, Selgrade MK, Vesper SJ, Ward MD. An extract of Stachybotrys chartarum causes allergic asthma-like responses in a BALB/c mouse model. Toxicol Sci. 2002;70(1):98–109. pmid:12388839.

32. Baxi SN, Muilenberg ML, Rogers CA, Sheehan WJ, Gaffin J, Permaul P, et al. Exposures to molds in school classrooms of children with asthma. Pediatr Allergy Immunol. 2013;24(7):697–703. pmid:24112429; PubMed Central PMCID: PMC3782748.

33. Bae HK, Shinozuka J, Islam Z, Pestka JJ. Satratoxin G interaction with 40S and 60S ribosomal subunits precedes apoptosis in the macrophage. Toxicol Appl Pharmacol. 2009;237(2):137–45. pmid:19306889; PubMed Central PMCID: PMC3124149.

34. Islam Z, Hegg CC, Bae HK, Pestka JJ. Satratoxin G-induced apoptosis in PC-12 neuronal cells is mediated by PKR and caspase independent. Toxicol Sci. 2008;105(1):142–52. pmid:18535002; PubMed Central PMCID: PMC2734304.

35. Pestka JJ, Yike I, Dearborn DG, Ward MD, Harkema JR. Stachybotrys chartarum, trichothecene mycotoxins, and damp building-related illness: new insights into a public health enigma. Toxicol Sci. 2008;104(1):4–26. pmid:18007011.

36. Swirski FK, Sajic D, Robbins CS, Gajewska BU, Jordana M, Stampfli MR. Chronic exposure to innocuous antigen in sensitized mice leads to suppressed airway eosinophilia that is reversed by granulocyte macrophage colony-stimulating factor. J Immunol. 2002;169(7):3499–506. pmid:12244139.

37. Huttunen K, Hyvarinen A, Nevalainen A, Komulainen H, Hirvonen MR. Production of proinflammatory mediators by indoor air bacteria and fungal spores in mouse and human cell lines. Environ Health Perspect. 2003;111(1):85–92. pmid:12515684; PubMed Central PMCID: PMC1241310.

38. Murtoniemi T, Nevalainen A, Suutari M, Toivola M, Komulainen H, Hirvonen MR. Induction of cytotoxicity and production of inflammatory mediators in raw264.7 macrophages by spores grown on six different plasterboards. Inhal Toxicol. 2001;13(3):233–47. pmid:11295859.

39. Shahan TA, Sorenson WG, Paulauskis JD, Morey R, Lewis DM. Concentration- and time-dependent upregulation and release of the cytokines MIP-2, KC, TNF, and MIP-1alpha in rat alveolar macrophages by fungal spores implicated in airway inflammation. Am J Respir Cell Mol Biol. 1998;18(3):435–40. pmid:9490662.

40. Ozdemir C. An immunological overview of allergen specific immunotherapy—subcutaneous and sublingual routes. Ther Adv Respir Dis. 2009;3(5):253–62. pmid:19880430.

41. Jarvis BB, Sorenson WG, Hintikka EL, Nikulin M, Zhou Y, Jiang J, et al. Study of toxin production by isolates of Stachybotrys chartarum and Memnoniella echinata isolated during a study of pulmonary hemosiderosis in infants. Appl Environ Microbiol. 1998;64(10):3620–5. pmid:9758776; PubMed Central PMCID: PMC106476.

42. Rao CY, Burge HA, Brain JD. The time course of responses to intratracheally instilled toxic Stachybotrys chartarum spores in rats. Mycopathologia. 2000;149(1):27–34. Epub 2001/03/03. pmid:11227851.

43. Chung YJ, Jarvis BB, Tak H, Pestka JJ. Immunochemical Assay for Satratoxin G and other Macrocyclic Trichothecenes Associated with Indoor Air Contamination by Stachybotrys chartarum. Toxicol Mech Methods. 2003;13(4):247–52. Epub 2003/01/01. pmid:20021149.

44. Crapo RO, Casaburi R, Coates AL, Enright PL, Hankinson JL, Irvin CG, et al. Guidelines for methacholine and exercise challenge testing-1999. Am J Respir Crit Care Med. 2000;161(1):309–29. pmid:10619836.

45. Gavin IM, Gillis B, Arbieva Z, Prabhakar BS. Identification of human cell responses to hexavalent chromium. Environ Mol Mutagen. 2007;48(8):650–7. pmid:17685462.

46. Gillis B, Gavin IM, Arbieva Z, King ST, Jayaraman S, Prabhakar BS. Identification of human cell responses to benzene and benzene metabolites. Genomics. 2007;90(3):324–33. Epub 2007/06/19. pmid:17572062.

47. Behm FG, Gavin IM, Karpenko O, Lindgren V, Gaitonde S, Gashkoff PA, et al. Unique immunologic patterns in fibromyalgia. BMC Clin Pathol. 2012;12(1):25. pmid:23245186; PubMed Central PMCID: PMC3534336.

48. Leavitt RD, Felsted RL, Bachur NR. Biological and biochemical properties of Phaseolus vulgaris isolectins. J Biol Chem. 1977;252(9):2961–6. pmid:853039.

49. Yachnin S, Svenson RH. The immunological and physicochemical properties of mitogenic proteins derived from Phaseolus vulgaris. Immunology. 1972;22(5):871–83. pmid:4623260; PubMed Central PMCID: PMC1407857.

50. Storey JD. The positive false discovery rate: a Bayesian interpretation and the q-value. Ann Stat. 2003;31(6):2013–35.

51. Storey JD. A direct approach to false discovery rates. J Roy Stat Soc Ser B (Stat Method). 2002;64(3):479–98.

52. Xiong M, Fang X, Zhao J. Biomarker identification by feature wrappers. Genome Res. 2001;11(11):1878–87. pmid:11691853; PubMed Central PMCID: PMC311150.

53. Vapnik VN. The Nature of Statistical Learning Theory. New York: Springer-Verlag New York Inc.; 2000.

54. Ramaswamy S, Tamayo P, Rifkin R, Mukherjee S, Yeang CH, Angelo M, et al. Multiclass cancer diagnosis using tumor gene expression signatures. Proc Natl Acad Sci U S A. 2001;98(26):15149–54. Epub 2001/12/14. pmid:11742071; PubMed Central PMCID: PMC64998.

55. Agranoff D, Fernandez-Reyes D, Papadopoulos MC, Rojas SA, Herbster M, Loosemore A, et al. Identification of diagnostic markers for tuberculosis by proteomic fingerprinting of serum. Lancet. 2006;368(9540):1012–21. Epub 2006/09/19. pmid:16980117.

56. Zhang F, Kaufman HL, Deng Y, Drabier R. Recursive SVM biomarker selection for early detection of breast cancer in peripheral blood. BMC Med Genomics. 2013;6 Suppl 1:S4. pmid:23369435; PubMed Central PMCID: PMC3552693.

57. Alkhalaf A, Zurbig P, Bakker SJ, Bilo HJ, Cerna M, Fischer C, et al. Multicentric validation of proteomic biomarkers in urine specific for diabetic nephropathy. PLoS One. 2010;5(10):e13421. Epub 2010/10/27. pmid:20975990; PubMed Central PMCID: PMC2958112.

58. Jin G, Zhou X, Wang H, Zhao H, Cui K, Zhang XS, et al. The knowledge-integrated network biomarkers discovery for major adverse cardiac events. J Proteome Res. 2008;7(9):4013–21. Epub 2008/07/31. pmid:18665624; PubMed Central PMCID: PMC2854538.

59. Qin N, Yang F, Li A, Prifti E, Chen Y, Shao L, et al. Alterations of the human gut microbiome in liver cirrhosis. Nature. 2014;513(7516):59–64. pmid:25079328.

60. Boardman M, Trappenberg T, editors. A heuristic for free parameter optimization with support vector machines. International Joint Conference on Neural Networks, 2006 IJCNN'06 2006: IEEE.

61. Kohavi R, John GH. Wrappers for feature subset selection. Artif Intell. 1997;97(1–2):273–324.

62. Peng S, Xu Q, Ling XB, Peng X, Du W, Chen L. Molecular classification of cancer types from microarray data using the combination of genetic algorithms and support vector machines. FEBS Lett. 2003;555(2):358–62. pmid:14644442.

63. Saeys Y, Inza I, Larranaga P. A review of feature selection techniques in bioinformatics. Bioinformatics. 2007;23(19):2507–17. Epub 2007/08/28. btm344 [pii] pmid:17720704.

64. Faraggi D, Reiser B. Estimation of the area under the ROC curve. Stat Med. 2002;21(20):3093–106. pmid:12369084.

65. Hastie T, Rosset S, Tibshirani R, Zhu J. The entire regularization path for the support vector machine. J Mach Learn Res. 2004:1391–415.

66. Chang C- C, Lin C- J. LIBSVM: a library for support vector machines. ACM Transactions on Intelligent Systems and Technology (TIST). 2011;2(3):27.

67. Troyanskaya O, Cantor M, Sherlock G, Brown P, Hastie T, Tibshirani R, et al. Missing value estimation methods for DNA microarrays. Bioinformatics. 2001;17(6):520–5. pmid:11395428.

68. Moons KG, Donders RA, Stijnen T, Harrell FE Jr. Using the outcome for imputation of missing predictor values was preferred. J Clin Epidemiol. 2006;59(10):1092–101. pmid:16980150.

69. Kim MS, Guest CC. Simulated annealing algorithm for binary phase only filters in pattern classification. Appl Opt. 1990;29(8):1203–8. pmid:20562979.

70. Ambroise C, McLachlan GJ. Selection bias in gene extraction on the basis of microarray gene-expression data. Proc Natl Acad Sci U S A. 2002;99(10):6562–6. pmid:11983868 PubMed Central PMCID: PMC124442.

71. Furey TS, Cristianini N, Duffy N, Bednarski DW, Schummer M, Haussler D. Support vector machine classification and validation of cancer tissue samples using microarray expression data. Bioinformatics. 2000;16(10):906–14. pmid:11120680.

72. Liu JJ, Cutler G, Li W, Pan Z, Peng S, Hoey T, et al. Multiclass cancer classification and biomarker discovery using GA-based algorithms. Bioinformatics. 2005;21(11):2691–7. pmid:15814557.

73. Fitzpatrick DR, Wilson CB. Methylation and demethylation in the regulation of genes, cells, and responses in the immune system. Clin Immunol. 2003;109(1):37–45. pmid:14585274.

74. Shore SA. Obesity and asthma: possible mechanisms. J Allergy Clin Immunol. 2008;121(5):1087–93. pmid:18405959.

75. Stumhofer JS, Silver JS, Laurence A, Porrett PM, Harris TH, Turka LA, et al. Interleukins 27 and 6 induce STAT3-mediated T cell production of interleukin 10. Nat Immunol. 2007;8(12):1363–71. pmid:17994025.

76. Saraiva M, O'Garra A. The regulation of IL-10 production by immune cells. Nat Rev Immunol. 2010;10(3):170–81. pmid:20154735.

77. Barnes PJ. The cytokine network in asthma and chronic obstructive pulmonary disease. J Clin Invest. 2008;118(11):3546–56. pmid:18982161; PubMed Central PMCID: PMC2575722.

78. Bisset LR, Schmid-Grendelmeier P. Chemokines and their receptors in the pathogenesis of allergic asthma: progress and perspective. Curr Opin Pulm Med. 2005;11(1):35–42. pmid:15591886.

79. Hope AP, Simon RA. Excess dampness and mold growth in homes: an evidence-based review of the aeroirritant effect and its potential causes. Allergy Asthma Proc. 2007;28(3):262–70. pmid:17619553.

80. Kanda N, Tani K, Enomoto U, Nakai K, Watanabe S. The skin fungus-induced Th1- and Th2-related cytokine, chemokine and prostaglandin E2 production in peripheral blood mononuclear cells from patients with atopic dermatitis and psoriasis vulgaris. Clin Exp Allergy. 2002;32(8):1243–50. pmid:12190666.

81. Kimura M, Tsuruta S, Yoshida T. Differences in cytokine production by peripheral blood mononuclear cells (PBMC) between patients with atopic dermatitis and bronchial asthma. Clin Exp Immunol. 1999;118(2):192–6. pmid:10540178; PubMed Central PMCID: PMC1905416.

82. Havaux X, Zeine A, Dits A, Denis O. A new mouse model of lung allergy induced by the spores of Alternaria alternata and Cladosporium herbarum molds. Clin Exp Immunol. 2005;139(2):179–88. pmid:15654816; PubMed Central PMCID: PMC1809297.

83. Hudson B, Flemming J, Sun G, Rand TG. Comparison of immunomodulator mRNA and protein expression in the lungs of Stachybotrys chartarum spore-exposed mice. J Toxicol Environ Health A. 2005;68(15):1321–35. Epub 2005/07/16. pmid:16020192.

84. Stark H, Roponen M, Purokivi M, Randell J, Tukiainen H, Hirvonen MR. Aspergillus fumigatus challenge increases cytokine levels in nasal lavage fluid. Inhal Toxicol. 2006;18(13):1033–9. pmid:16966303.

85. Yike I, Rand TG, Dearborn DG. Acute inflammatory responses to Stachybotrys chartarum in the lungs of infant rats: time course and possible mechanisms. Toxicol Sci. 2005;84(2):408–17. Epub 2005/01/14. pmid:15647601.

86. Inoue Y, Matsuwaki Y, Shin SH, Ponikau JU, Kita H. Nonpathogenic, environmental fungi induce activation and degranulation of human eosinophils. J Immunol. 2005;175(8):5439–47. pmid:16210651.

87. Kauffman HF, Tomee JF, van de Riet MA, Timmerman AJ, Borger P. Protease-dependent activation of epithelial cells by fungal allergens leads to morphologic changes and cytokine production. J Allergy Clin Immunol. 2000;105(6 Pt 1):1185–93. pmid:10856154.

88. Shalit I, Halperin D, Haite D, Levitov A, Romano J, Osherov N, et al. Anti-inflammatory effects of moxifloxacin on IL-8, IL-1beta and TNF-alpha secretion and NFkappaB and MAP-kinase activation in human monocytes stimulated with Aspergillus fumigatus. J Antimicrob Chemother. 2006;57(2):230–5. pmid:16352735.

89. Grazziutti ML, Rex JH, Cowart RE, Anaissie EJ, Ford A, Savary CA. Aspergillus fumigatus conidia induce a Th1-type cytokine response. J Infect Dis. 1997;176(6):1579–83. pmid:9395370.

90. Lappalainen MH, Hyvarinen A, Hirvonen MR, Rintala H, Roivainen J, Renz H, et al. High indoor microbial levels are associated with reduced Th1 cytokine secretion capacity in infancy. Int Arch Allergy Immunol. 2012;159(2):194–203. pmid:22678428.

91. de Mello LM, Bechara MI, Sole D, Rodrigues V. TH1/TH2 balance in concomitant immediate and delayed-type hypersensitivity diseases. Immunol Lett. 2009;124(2):88–94. pmid:19433108.

92. O'Driscoll BR, Hopkinson LC, Denning DW. Mold sensitization is common amongst patients with severe asthma requiring multiple hospital admissions. BMC Pulm Med. 2005;5(1):4. pmid:15720706; PubMed Central PMCID: PMC550663.

93. Quansah R, Jaakkola MS, Hugg TT, Heikkinen SA, Jaakkola JJ. Residential dampness and molds and the risk of developing asthma: a systematic review and meta-analysis. PLoS One. 2012;7(11):e47526. pmid:23144822; PubMed Central PMCID: PMC3492391.