Determining the fungal diversity within an environment is a method utilized in occupational health studies to identify health hazards. This protocol describes DNA extraction from occupational air samples for amplification and sequencing of fungal ITS regions. This approach detects many fungal species that can be overlooked by traditional assessment methods.
Traditional methods of identifying fungal exposures in occupational environments, such as culture and microscopy-based approaches, have several limitations that have resulted in the exclusion of many species. Advances in the field over the last two decades have led occupational health researchers to turn to molecular-based approaches for identifying fungal hazards. These methods have resulted in the detection of many species within indoor and occupational environments that have not been detected using traditional methods. This protocol details an approach for determining fungal diversity within air samples through genomic DNA extraction, amplification, sequencing, and taxonomic identification of fungal internal transcribed spacer (ITS) regions. ITS sequencing results in the detection of many fungal species that are either not detected or difficult to identify to species level using culture or microscopy. While these methods do not provide quantitative measures of fungal burden, they offer a new approach to hazard identification and can be used to determine overall species richness and diversity within an occupational environment.
Fungal exposures in indoor and occupational environments can result in respiratory morbidities, including allergic sensitization and asthma1. Identification of fungal hazards is important for assessing risk and preventing worker exposure. These fungal hazards may be a result of indoor contamination, outdoor air intrusion, or environmental disturbances that result in the transport of fungal materials into areas where workers are present2. Methods to assess fungal exposure have included viable culture sampling as well as microscopic identification of fungal spores. These approaches have several limitations and often overlook many fungal species that could be contributing to the overall fungal burden3. Culture-based approaches can only differentiate those viable fungal organisms that can be cultivated on nutrient media. Identifying fungal spores to species level via microscopy can be confounded by spores sharing similar morphologies. Both methods are highly dependent on mycologists to analyze and identify the fungal species, with many remaining unidentified.
To improve upon existing methodologies used in occupational hazard identification and exposure assessments, many researchers have turned to molecular-based technologies. Sequencing-based approaches for assessing microbial diversity within indoor and occupational environments have revealed a broader spectrum of fungal species encountered compared to methods such as microscopy and viable culture3,4,5. The method presented here describes the air sampling of occupational environments and extraction of genomic DNA for the identification of potential fungal hazards. Hazard identification is accomplished by sequencing the nuclear ribosomal internal transcribed spacer, or ITS, regions that are highly variable among fungi and have been commonly used to differentiate fungal species6,7,8,9. Many species found in occupational settings, such as some species belonging to the phylum Basidiomycota, are not identifiable in viable culture and are difficult to differentiate microscopically. These fungi have been observed in high relative abundance within indoor and occupational environments assessed by sequencing fungal ITS regions3,4,10. ITS sequencing has provided greater knowledge into the diversity of fungi encountered within indoor and occupational environments.
The protocol described here details the methods used to collect, extract, and amplify fungal ITS regions from bioaerosols for sequence analysis. This approach utilizes the National Institute for Occupational Safety and Health (NIOSH) two-stage cyclone aerosol sampler to collect particulates in the air. This sampler was developed to collect bioaerosols and separate respirable (≤4 µm aerodynamic diameter) and non-respirable (>4 µm aerodynamic diameter) particles, which allows for identification of fungal organisms within indoor environments that are most likely to be inhaled by a worker11. Other air samplers, including cyclone samplers, are available on the market that have the ability to collect particles within the respirable range (<4 µm) using filters12,13. In contrast, the NIOSH two-stage cyclone aerosol sampler separates fungal species based on their aerodynamic diameter into disposable, polypropylene tubes that can be immediately processed for downstream applications14.
The processes of extracting genomic DNA and amplifying the fungal ITS regions are detailed in this protocol. The extraction methodologies presented have been developed specifically for the extraction of genomic DNA from fungi and bacteria, as many commercial kits target mammalian cells, bacteria, or specifically yeasts15. The primers used in this study are selected based on their overall coverage of both the fungal ITS 1 and ITS 2 regions4,5. Sequencing of these regions allows for the comparison of many banked ITS sequences, including those that sequence the ITS 1 region, the ITS 2 region, or both the ITS 1 and ITS 2 regions. The fungal diversity of air samples collected in an indoor setting using these methods are shown, revealing a substantial number of sequences placed in the phyla Ascomycota and Basidiomycota as well as other sequences belonging to less dominant fungal phyla, such as Zygomycota. The broad diversity of fungal sequences identified using this approach would not be captured using traditional hazard identification methodologies like cultivation or microscopy. Sequencing of fungal ITS regions provides an enhanced method to identify fungal hazards and allow for a better understanding of indoor and occupational fungal exposures.
1. Preparing the NIOSH aerosol sampler
NOTE: The NIOSH aerosol sampler is a two-stage cyclone aerosol sampler that collects bioaerosols using two sampling tubes and a polytetrafluoroethylene (PTFE) filter.
2. Static and personal aerosol sampling
3. Extraction of genomic DNA from air samples
4. Amplification of fungal ribosomal DNA
5. Verification of fungal ITS amplification using agarose gel electrophoresis
6. Sequencing and analysis of fungal ITS regions
The species distribution within an environment can be assessed using relative abundance by determining the number of clones of each OTU identified in the air samples. Figure 6 is a Krona chart representing the taxonomically placed species within an indoor environment following 60 min of air sampling. It can be observed that the environment contains a variety of species within two major fungal phyla, Ascomycota and Basidiomycota, as well as species belonging to the phylum Zygomycota (Rhizopus microsporus). Traditional methods of assessment are biased toward Ascomycota species, as many Basidiomycota are not culturable or cannot be differentiated using microscopy-based methods. Sequencing of fungal ITS regions allows for the detection of numerous fungal species, specifically some Basidiomycota, that often go undetected. Analysis of this environment shows 68% of the fungal sequences identified belonged to the Basidiomycota with most of them placed in the order Polyporales.
The taxonomic data obtained using sequencing-based approaches can be used to determine the taxonomic diversity within an environment. Diversity indices, such as those presented in Table 1, determine how rich the environment is, which describes the number of species identified in each sample, as well as how diverse the environment is, which takes into account the number of species and the abundance of each. Table 1 shows the Chao 1 richness indices and Shannon diversity indices for two indoor environments. The indices indicate higher richness and diversity in air-conditioned environments compared to evaporative cooler environments. Also included are Bray-Curtis dissimilarity coefficients. These describe how dissimilar samples within the environment are regarding the species identified within each sample. The samples within both of the environments described in Table 1 are highly dissimilar at 98% and 97% for air conditioner and evaporative cooler environments, respectively.
Figure 1: Schematic representation of a three-piece filter cassette assembly. The filter cassette is made up of three pieces: a top or cap piece (inlet), an extension cowl and a base piece (outlet). A support pad and filter are placed on the gridded surface of the base piece. The extension cowl is then sealed onto the base using a manual or pneumatic press. For use with the NIOSH bioaerosol cyclone sampler, the top piece is left off during sampling and is then replaced for storage and transport. Please click here to view a larger version of this figure.
Figure 2: Schematic representation of the NIOSH bioaerosol cyclone sampler. The NIOSH cyclone aerosol sampler collects airborne particulate by drawing air into the sampler through the inlet, and producing a cyclone that deposits particles on the wall of the sample tubes. Air is first drawn into a 15 mL polypropylene tube, where large particles (>4 µm) collect on the walls of the tube. Air then flows out of the first tube and is drawn into a 1.5 mL microcentrifuge tube, where smaller particles (1 to 4 µm) collect. The remaining particles (<1 µm) collect on a PTFE filter as the air is drawn out of the sampler. Please click here to view a larger version of this figure.
Figure 3: Example of a static sampler set-up. The image demonstrates samplers set up for static area sampling. The NIOSH samplers are set up to collect 40 inches (102 cm) and 60 inches (152 cm) from the floor, representing heights of a sitting and standing adult, respectively. Please click here to view a larger version of this figure.
Figure 4: Example of a personal sampler assembly. The image illustrates a NIOSH sampler worn on a backpack by an individual. The sampler and power switch are located on the straps of the pack while the sampling pump is located within the backpack itself. Please click here to view a larger version of this figure.
Figure 5: Example of an agarose gel visualizing amplified ITS DNA from air samples. The image shows DNA bands between 750 and 1,000 base pairs. These bands represent amplified fungal ITS regions from extracted air samples. ITS regions vary in size depending on the species origin. Please click here to view a larger version of this figure.
Figure 6: Example Krona chart presenting taxonomically placed fungal species identified within an indoor environment. The Krona chart depicts the relative abundance of fungal species found in 4 samples collected from homes following a 60 min air sampling period with the NIOSH bioaerosol cyclone sampler. 68% of the species identified belong to the phylum Basidiomycota with the remaining belonging to the Ascomycota (33%) and Zygomycota (1%). The most abundant ITS sequences in this environment belonged to the order Polyporales and were identified as Phanerodontia chrysosporium and Gelatoporia dichroa. Please click here to view a larger version of this figure.
Chao-1 Richness | Shannon Diversity | Bray-Curtis Distance | |
mean (min, max) | mean (min, max) | mean (min, max) | |
Air Conditioner (n = 9) | 33.86 (15.0, 102.0) | 1.39 (0.68, 2.51) | 0.9778 (0.8313, 1.00) |
Evaporative Cooler (n = 10) | 25.46 (12.5, 49.5) | 1.10 (0.50, 1.72) | 0.9624 (0.3804, 1.00) |
Table 1: Example data presenting diversity and species richness indices comparing indoor environments. The table shows Chao 1 richness indices (number of species per sample), Shannon diversity indices (number of species and abundance of each), and Bray-Curtis dissimilarity coefficients (how dissimilar the species distribution is between samples on a scale of 0 – 1). These data demonstrate higher richness and diversity in air-conditioned environments and high dissimilarity among samples in both environments. This table has been modified from Lemons et al.10 with permission from the Royal Society of Chemistry.
Determining the fungal diversity within an occupational environment using sequencing-based approaches has improved fungal hazard identification and exposure assessment. Using this approach has allowed for the detection of many additional fungal species that are often not detected using culture or microscopy-based methods of assessment. A method for sampling bioaerosols from occupational and indoor environments and the extraction of genomic DNA from air samples for ITS amplification and sequencing is presented here. Determinations of fungal diversity using these methods is highly dependent on: (1) successful and complete extraction of genomic DNA from the air samples, (2) amplification of the gDNA with primers that allow for comparison of the amplified ITS sequences to banked sequences in the database and limit amplification biases, and (3) correct taxonomic identification of the sequences within the databases.
Successful genomic DNA extraction is reliant on the extraction methodologies and reagents employed in a study. As is mentioned in the protocol, a lot of the bioaerosols and other biological and non-biological particulate collected in the tube phases of the NIOSH two-stage cyclone sampler collect at the top of the tube. It is very important that the tubes be carefully opened when adding the lysis buffer as not to lose any particulate and that they be vortexed in a manner that allows for all the particulate to fall into the lysis buffer (right-side up and upside down). It is also important that the filter be carefully handled using aseptic methods as not to disrupt any bioaerosols that have collected on the surface prior to extraction or introduce contaminating DNA. It has been demonstrated that there is a large amount of variation between many of the commercially available genomic DNA extraction kits15. Some extraction procedures bias toward specific fungal species and result in varying DNA yields20. If the DNA yield following extraction is low, higher PCR reaction replicates can be performed. Not only do the extraction yields vary, but many of the kits contribute a substantial amount of microbial DNA contamination. For this reason, it is important to include proper controls throughout the procedure to identify and remove potential reagent contaminants from the analysis.
ITS amplification is another critical step of the protocol. The extraction methodologies used can vary in their ability to remove PCR inhibitors from the environmental samples. This is particularly concerning with environmental dust samples21. PCR amplification of samples extracted using the method presented here exhibited some PCR inhibition after the addition of 10 mg of dust15. Although PCR inhibition is of greatest concern in environmental dust samples, it should still be a consideration when amplifying DNA extracted from air samples for determining fungal diversity. The PCR primer set used in this method provides coverage of both the ITS 1 and ITS 2 regions. This allows for comparison to many of the sequences placed in the sequence databases. As with the extraction procedure, many amplification and primer biases have been previously described22. The fungal genome size and gene copy number could also influence ITS amplification23. These biases should be taken into consideration when interpreting the resultant sequence data. Taxonomic identification of ITS sequences is perhaps the most critical component of this method. It is important to keep identification criteria consistent. The ability to identify ITS sequences is dependent on the sequences banked in the databases. Many sequences cannot be identified to the species level. Having specific criteria when making taxonomic identifications based on banked sequences, like percentage identity cutoffs, is critical for keeping datasets consistent from study to study.
Compared to traditional assessment approaches, extracting and sequencing fungal DNA from occupational air samples provides a more complete representation of fungal diversity within an environment. In addition to the limitations discussed above, it must be noted that the results of these types of analyses are semi-quantitative unlike culture, spore counts, or quantitative PCR24 that can yield quantitative data. Sequencing data can be used to determine relative abundance per sample as well as calculate fungal diversity metrics. These methods can be used in conjunction with other methods, like qPCR, to get more quantifiable data. The datasets created using these approaches can be used in occupational health studies to gain a better understanding of fungal hazards in specific occupational environments. Developing more standardized approaches of extraction and primer selection are necessary to better characterize and compare sequencing-based studies of fungal diversity.
The authors have nothing to disclose.
This work was supported in part by an interagency agreement between NIOSH and NIEHS (AES12007001-1-0-6).
NIOSH BC251 bioaerosol cyclone sampler | NIOSH | BC251 | The NIOSH sampler is not yet commercially available. Please contact William Lindsley, PhD (wlindsley@cdc.gov) for information on obtaining the NIOSH sampler |
Fisherbrand Sterile Microcentrifuge Tubes with Screw Caps | Fisher Scientific | 02-681-373 | 1.5 mL polypropylene microcentrifuge tubes for air sampling; screw top threading must match the threading of the NIOSH sampler |
Falcon 15 mL Conical Centrifuge Tubes | Corning | 352096 | 15 mL polypropylene tubes for air sampling |
Clean Room Vinyl Tape, Easy-Remove, 1/4" Width | McMaster-Carr | 76505A1 | sealing tape |
Filter Cassette, Clear Styrene, 37 mm | SKC Inc. | 225-3LF | 3-piece sampling cassette (no filter). Contains: cassette base, extension cowl, cassette cap and inlet/outlet plugs |
PTFE hydrophobic fluoropore membrane filters, 3.0 µm, 37 mm | EMD Millipore | FSLW03700 | Contains: 37 mm, 3.0 µm PTFE filters and support pads |
Fisherbrand filter forceps | Fisher Scientific | 09-753-50 | filter forceps |
Model 502 Precision PanaPress | PanaVise | 502 | pneumatic cassette press is constructed from this precision arbor press |
Scotch Super 33+ vinyl electrical tape | McMaster-Carr | 76455A21 | 19 mm tape |
Multi-purpose Calibration Jar, Large | SKC Inc. | 225-112 | calibration jar |
Universal PCXR4 Sample Pump | SKC Inc. | 224-PCXR4 | sampling pump |
Mass Flowmeter 4140 | TSI Inc. | 4140 | flow meter |
Roche High Pure PCR Template Kit | Roche Diagnostics | 11796828001 | Kit used for genomic DNA extraction. Contains: Lysis buffer, Binding buffer, Proteinase K, Inhibitor removal buffer, Wash buffer, Elution buffer, Glass fiber filter tubes and 2 ml collection tubes |
Fisherbrand 2 mL Reinforced Polypropylene Screw Cap Tubes with Caps | Fisher Scientific | 15340162 | 2 mL reinforced tubes for bead homogenization |
Glass beads, acid washed, 212-300 µm | Sigma-Aldrich | G1277 | glass beads |
Fisher Scientific Bead Mill 24 Homogenizer | Fisher Scientific | 15-340-163 | bead homogenizer |
CelLytic B Cell Lysis Reagent, 10X | Sigma-Aldrich | C8740 | lysis reagent |
Platinum Taq polymerase | Invitrogen | 10966-018 | Contains: Platinum Taq polymerase, 10X PCR buffer (no MgCl2), 50 mM MgCl2, KB Extender |
dNTP Mix | Invitrogen | 18427-088 | 10 mM dNTP mix |
QIAquick PCR Purification Kit | Qiagen | 28106 | Kit used to purify fungal amplicons. Contains: Buffer PB (binding buffer), Buffer PE (washing buffer), Buffer EB (elution buffer), pH Indicator dye (optional), and GelPilot loading dye |
Owl EasyCast Mini Gel Electrophoresis System | Thermo Fisher | B1 or B2 | |
TrackIt 1 KB Plus DNA Ladder | Thermo Fisher | 10488-085 | DNA ladder |