Meeting report: fungal its workshop (october 2012).

This report summarizes a meeting held in Boulder, CO USA (19-20 October 2012) on fungal community analyses using ultra-high-throughput sequencing of the internal transcribed spacer (ITS) region of the nuclear ribosomal RNA (rRNA) genes. The meeting was organized as a two-day workshop, with the primary goal of supporting collaboration among researchers for improving fungal ITS sequence resources and developing recommendations for standard ITS primers for the research community.

ITS sequence database for fungal community analyses

Sequencing-based techniques have allowed characterization of microbial communities from environmental samples without relying on cultivation. Because of improvements in these techniques and of methods required to analyze the resulting datasets, microbial ecology has been advancing rapidly: microbial diversity can be surveyed to an extent previously unimaginable. The internal transcribed spacer (ITS) and other ribosomal RNA gene sequence regions (“rRNA” in the rest of this document) have been used successfully to profile fungal communities using Sanger [1] and 454 [2] sequencing. The arrival of ultra-high-throughput sequencing platforms [3] promise to offer new insights into the diversity and ecology of fungal communities, yet few studies of fungal communities have employed this technology successfully. Progress has been hampered, in part, by the lack of a high-quality reference database of fungal sequences for the ITS region of the rRNA operon [4], which is now the most widely sequenced DNA region in fungi [5] and is the marker of choice for molecular identification of most fungal taxa [6].

For and , where the ribosomal small subunit (SSU/16S) gene is the primary marker in environmental sequencing, efforts have been made to improve the quality of the public reference sequence datasets, including GreenGenes [7] and RDP [8]. The same is true for more general SSU/large subunit (LSU) rRNA gene sequence databases, such as SILVA, which includes all three domains of life [9]. For fungi, a “hand-curated” LSU sequence reference set is currently available, and work is underway to apply similar methods to improve the ITS database [10]. Because sequencing technologies continue to improve, the number of ITS sequences in primary sequence repositories such as INSDC will steadily increase, and quality control via hand-curation for specialized, publically available rRNA gene sequence databases will not be sustainable.

Primary sequence repositories are already experiencing explosive growth in the number of unidentified environmental fungal ITS sequences [11], yet these sequences will be of limited use in improving our overall understanding of fungal diversity unless they are properly identified and can be placed within a phylogenetic context. When sufficiently closely related sequences exist, environmental sequences can be placed within a phylogenetic context today simply by aligning with related sequences and constructing trees. However, because ITS evolves rapidly, constructing phylogenies that span large taxonomic ranges remains extremely challenging. An even more important problem is that of misidentified sequences (environmental sequences included) currently in public databases. These can lead to erroneous placement of unknowns, even if tree-based approaches are used. Therefore, identifying these errors and re-annotating in an automated fashion is a critical challenge, especially for extremely large datasets where manual phylogenetic analysis is not feasible due to the presence of millions of sequence reads that correspond to “unknown” operational taxonomic units (OTUs). Consequently, clean reference databases as well as automated phylogenetic assignment and analysis methods are critical needs.

Facilitate communication, potential data exchange, and collaboration with the aim of improving fungal ITS sequence resources for the research community.

Identify suitable ITS primers for fungal community analyses using ultra-high-throughput sequencing.

Develop strategies for automated (and manual) database curation as well as the naming of environmental sequences and OTUs at various levels of resolution.

Establish a sustainable plan for reference database development and maintenance.

Participants

The meeting participants included researchers representing publicly available databases that contain microbial sequence data (e.g., GenBank, GreenGenes, RDP, SILVA) or fungi-specific resources (e.g., MycoBank, UNITE), as well as researchers currently using ultra-high-throughput sequencing to examine fungal communities or those involved in developing software, such as QIIME [12] and PhyloSift [13], to facilitate such studies.

Activities

The meeting was conducted as a two-day workshop. The first day was devoted primarily to brief presentations by participants outlining their involvement in curating public sequence databases, developing high-throughput sequencing pipelines, or using ultra-high-throughput sequencing to examine fungal diversity in environmental samples (e.g., air or soil). The presentations are available online [14]. The second day focused on discussions related to the assembly of a high-quality reference database of fungal ITS sequences, selecting ITS primers suitable for ultra-high-throughput sequencing, as well as methods to link ITS sequences to the fungal phylogeny for automated curation, quality control, and phylogeny-based community analysis methods.

Conclusions / Outcomes

Ultra-high-throughput sequence processing/analytical pipelines, such as those implemented in QIIME, rely on de-replication of large sequence datasets through clustering for the creation of reference sequence sets that can be used to assist in the recognition of OTUs from environmental samples. The meeting participants largely supported the use of the UNITE database [15,16] as a focal point for the development of high quality fungal ITS reference sequence sets. UNITE currently has implemented several desirable features for this task, including:

A comprehensive set of approximately 300,000 fungal ITS sequences extracted from public databases.

An annotation management system (PlutoF) that allows qualified third-party users to add pertinent metadata (e.g., on ecology or geography), improve the taxonomic resolution, tag problematic entries, or correct misidentifications for sequences in the UNITE database.

Global Key Annotations that permit visualization of sequences clustered at a range of similarity levels (99–97%), with the sequences representing each OTU in the cluster depicted in an alignment.

Cluster centroids selected with a preference toward using sequences that are reliable (e.g., generated from trustworthy sources such as the Assembling the Fungal Tree of Life project or hand-picked by taxonomic experts), taxonomically informative (e.g., identified to the species-level), or that have particular relevance for taxonomy/nomenclature (e.g., sequences from type specimens).

Improvements on the horizon (slated for early 2013) included labeling sequences representing cluster centroids that will allow these unique sequences to be tracked through time as clusters change, as well as options for downloading cluster centroid sequence sets for different sequence similarity levels.

With the availability of these cluster centroids, reference sequence sets for ultra-high-throughput pipelines can be directly generated from the UNITE database in a rapid manner. The meeting allowed for coordination that resulted in the creation of an alpha version of the UNITE reference set to facilitate OTU picking and taxonomic assignment for fungal ITS sequence reads generated in high-/ultra-high-throughput sequencing runs. This reference set is now publically available on the QIIME website [12,17].

UNITE currently provides taxonomic strings based on classification schema culled from fungal taxonomic resources such as Index Fungorum [18] and MycoBank [19], which are comprehensive databases for fungal names that offer expertise and resources for improving the quality and availability of fungal taxonomic information. Journals publishing novel fungal taxa now typically require authors to register new names in MycoBank, which in turn is encouraging submission of informative DNA sequences, such as the ITS region, associated with new taxa to the public databases. In addition to acting as sequence vouchers for type material, these data also have the potential to inform molecular studies examining environmental samples.

Synergistic collaboration and the flow of information between, and within, online taxonomic resources and the public sequence databases (that have expressed interest in using a global standardized taxonomy) were seen by meeting participants as being highly desirable. The integration of fungal taxonomy and phylogeny was deemed another important consideration. The Assembling the Fungal Tree of Life (AFToL) project made considerable progress toward refining our understanding of the fungal phylogeny, which informed taxonomy for the kingdom [20]. Sequences generated under AFToL represent reliable data that are desirable for cluster centroids in the fungal reference sets. New projects, such as the Open Tree of Life [21], hold potential as phylogeny-based taxonomic resources for reference databases, while others, such as Fungal Barcoding [22], promise to produce additional high-quality sequences across a wide range of fungal groups for improving the reference datasets.

Although the ITS marker allows fungal sequences to be resolved at the genus- or species-level, aligning ITS sequences across a wide taxonomic range is essentially unworkable. Thus, meeting participants also discussed strategies for anchoring fungal ITS sequences to the broader phylogeny (e.g., one constrained to the AFToL backbone) using SSU or LSU sequences having contiguous ITS reads. Given that for fungi SSU (18S) is typically uninformative below the order- or family-level, LSU (26/28S) was seen as an appropriate marker for this task. The large subunit rRNA gene was considered desirable, as it has been extensively used in fungal phylogenetic studies, it allows for accurate placement in the phylogeny both at higher and lower taxonomic ranks (e.g., phylum/class and family/genus, respectively), and many reliable AFToL sequences (spanning both LSU and the ITS region) are currently available for this purpose [10,23]. Other reliable sequence sources for full-length LSU/ITS reads were identified from complete genomes, such as fungal genomes project [24]. With the anchor tree established, fungal ITS reads produced by ultra-high-throughput sequencing can potentially be placed within a phylogenetic context for taxonomic validation, improving the taxonomy associated with unknown environmental sequences, identifying and naming novel environmental groups, as well as for other automated curation tasks. Linking ITS reads to the fungal phylogeny will also allow for phylogenetic metrics of community distances (e.g., UniFrac) [25] to be used in beta-diversity analyses.

Two talks on the first day of the workshop presented preliminary data from ultra-high-throughput sequence surveys of soil fungi targeting the ITS1 region using the primer pairs ITS1-F/ITS2 [see 26]. Due to spliceosomal inserts known to exist toward the 3′ end of SSU rRNA gene that could interfere with priming sites and cause biases against groups of fungi (e.g., Helotiales) where such inserts exits, participants expressed a preference for using primers for the ITS2 region, for which such spliceosomal inserts are not known. Additional reasons for adopting the use of the ITS2 region marker include its close proximity to LSU (e.g., for anchoring to the phylogeny, see above), less variation in read length compared to ITS1, and the availability of data on ITS2 secondary structure [27,28] that can inform sequence alignments. With continually improving read lengths of ultra-high-throughput sequencing platforms, full length ITS/LSU reads may be possible in the near future. Primers targeting the ITS2 region have been identified, and are currently being tested for use in ultra-high-throughput sequencing. Recommendations for fungal ITS2 as well as current versions of the UNITE centroid reference sequence sets will be available on the QIIME website [12] in the coming year.