MycoCosm portal: gearing up for 1000 fungal genomes.

MycoCosm is a fungal genomics portal (http://jgi.doe.gov/fungi), developed by the US Department of Energy Joint Genome Institute to support integration, analysis and dissemination of fungal genome sequences and other 'omics' data by providing interactive web-based tools. MycoCosm also promotes and facilitates user community participation through the nomination of new species of fungi for sequencing, and the annotation and analysis of resulting data. By efficiently filling gaps in the Fungal Tree of Life, MycoCosm will help address important problems associated with energy and the environment, taking advantage of growing fungal genomics resources.

INTRODUCTION

Genomics is a powerful tool for understanding the genomic basis of biological processes in individual organisms and their communities. The Fungal Genomics Program (http://jgi.doe.gov/fungi) of the US Department of Energy (DOE) Joint Genome Institute (JGI) has partnered with the international scientific community to support genomic exploration of the kingdom Fungi and to help address problems related to energy and the environment. The exponential growth of data resulting from several large-scale genomics initiatives like the 1000 Fungal Genomes Project requires new approaches for data integration and analysis (1). We designed MycoCosm, an integrated fungal genomics resource, to allow researchers to easily access large amounts of genomics data and analyze them using web-based analytical tools. The MycoCosm portal enables in-depth multidimensional analysis of individual genomes and efficient comparative genomics of fungi, which may be applied to phylogenetically related fungi, or to those sharing the same lifestyle, ecological niche or phenotypic trait.

The first JGI fungal portal was created for the white rot fungus Phanerochaete chrysosporium in 2004 (2). To enable comparative genome analysis within the same analytical framework, in 2008 JGI integrated data from six genomes of Dothideomycete plant pathogens (three sequenced by JGI and three others provided by the Dothideomycetes Genome Consortium). This served as a prototype for MycoCosm, the first official version of which was released in March 2010 in response to a call for a pan-fungal genomic resource from the community of fungal biologists brought together by the Burroughs Wellcome Fund (Alexandria, VA). By 2011, MycoCosm included over 100 fungal genomes (3), and continues to grow. Here we provide an update on the latest developments of MycoCosm, which includes flexible tools for the analysis of exponentially growing genomics and other ‘omics’ data sets, and extends user community participation to the entire genome project cycle, from nomination of species for sequencing through genome annotation and comparative genomics.

Significantly larger and better-integrated genomics data

Since our last report (3), MycoCosm has grown to host >250 publicly available fungal genomes (Figure 1), 80% of which were sequenced by JGI and 20% by other groups added for comparative purposes. Individual genomes, or biologically relevant groups of them, may be searched using keywords or sequence, investigated using MycoCosm analytical tools, and downloaded for custom analyses.

Figure 1.

Growth of annotated genomes in MycoCosm. The genomes sequenced by JGI are shown in blue and those sequenced by others are shown in red.

Also in the past two years, all MycoCosm nodes have been populated with newly sequenced and annotated genomes including the recently activated nodes for Glomeromycota, Entomophthoromycotina, Kickxellomycotina, Blastocladiomycota, Neocallimastigomycota and others. Several new branches such as Xylonomycetes, Zoopagomycotina and Cryptomycota have been added for a more complete picture of the Fungal Tree of Life (Figure 2). Several model fungi have been added to MycoCosm, for which genes are linked to the corresponding records in curated databases such as the Saccharomyces Genome Database (4), the Aspergillus Genome Database (5) and others. GenBank records for JGI-sequenced fungal genomes are linked to MycoCosm using “JGIDB” qualifiers (http://www.ncbi.nlm.nih.gov/genbank/collab/db_xref/).

Figure 2.

MycoCosm front page. The MycoCosm tree illustrates relationships between major groups of fungi. Each node of the tree is linked to the corresponding group page and portals for individual genomes within the group.

Each genome in MycoCosm has its own portal where genomics and other ‘omics’ data can be explored in depth, in a single-genome context (Figure 3). Genomes in MycoCosm are also grouped for comparative genomics analyses (Figure 4), into two kinds of groups: phylogenetic and ecological. Phylogenetic groups are represented (and linked to) by the nodes of the MycoCosm tree, and facilitate comparative analyses of organisms with a shared evolutionary history. Ecological groups represent a curation effort in which fungi of similar lifestyles (e.g. endophytes, mycorrhizal symbionts, plant pathgens, etc.) have been assembled into groups that are linked to the chapters of the Genomics Encyclopedia of Fungi (6).

Figure 3.

Genome Browser for Neurospora crassa. Horizontal tracks represent genomic features of a given assembly scaffold. Red bars indicate N. crassa genes predicted by the Broad Institute and uploaded into MycoCosm. The four tracks with brown bars illustrate SNPs derived from resequencing of mutant strains of N. crassa (7). The two pink curves represent transcriptome profiling using RNA-Seq during N. crassa growth under two different conditions. Additional tracks can be displayed and configured using the toolbar.

Figure 4.

Examples of MycoCosm group pages. A PhyloGroup containing seven annotated genomes of fungi from Mucoromycotina is shown on the left. An EcoGroup with seven annotated genomes of phylogenetically diverse endophytes is shown on the right.

To encourage the research community to participate in the further growth of MycoCosm, users can now nominate diverse fungal species for sequencing (Figure 5) in the context of the 1000 Fungal Genomes Project (http://1000.fungalgenomes.org). Nominated species will be peer reviewed, and, if approved, will be sequenced, annotated and fully integrated into MycoCosm’s comparative genomics framework.

Figure 5.

1000 Fungal Genomes nomination page. Over 500 families of fungi and associated genome projects are displayed. Green and yellow rows highlight families with no, or only one, sequenced genome, where additional species can be nominated for sequencing via the use of the ‘Nominate’ button. Nomination requires user registration and commitment to provide JGI with DNA and RNA samples for sequencing.

Uniform genome annotation

For genome comparison to produce meaningful results, application of the same set of functional annotation databases to all genomes is required. For genomes sequenced and annotated outside of JGI, the existing gene models are fed into the functional annotation part of the pipeline, using the same databases used for annotation of JGI genomes, to ensure consistent references for genome comparisons.

All fungi sequenced at JGI are annotated using the JGI Annotation Pipeline, which integrates an array of tools for gene prediction, annotation and analysis (8). Using repeat-masked genome assemblies as input, protein-coding gene models are predicted using a combination of gene predictors of three general kinds: ab initio, homology-based and transcriptome-based. The ab initio Fgenesh (9) and GeneMark (10) are specifically trained for the given genome. Homology-based Fgenesh+ (9) and Genewise (11) are seeded by BLASTx (12) alignments of fungal sequences in NCBI’s nonredundant protein database. Transcriptome-based models are predicted from assembled RNA-sequencing data derived from the cultured fungus. For each genomic locus, the best representative gene model is selected based on a combination of protein homology and transcriptome support; these make up a ‘gene catalog’. In addition to protein-coding genes, tRNAs are predicted using tRNAscan-SE (13).

All predicted proteins are then functionally annotated. SignalP (14) is used to detect the sequence motifs responsible for protein localization, TMHMM (15) identifies possible transmembrane domains and InterProScan (16) predicts functional domains from Pfam (17) and other databases. Protein alignments to the NCBI’s nonredundant (http://www.ncbi.nlm.nih.gov), Swiss-Prot (18), KEGG (19) and KOG (20) databases additionally facilitate functional interpretation. Interpro, KEGG and Swiss-Prot hits are used to map gene ontology (GO) terms (21). Detailed annotation reports from these methods are provided on every coding gene’s protein page. Additionally, summaries listing numbers of genes by category in the GO, KEGG and KOG classifications are accessible from the portal menu, and can be compared with other selected genomes to explore gene family expansions and contractions across genomes.

Predicted gene models, along with different lines of supporting evidence including transcriptomics and similarity to proteins in other genomes, are displayed on a genome browser (Figure 3). The browser, which is based on the UCSC Genome Browser (22), displays genomic features as horizontal tracks, which can be customized to the needs of each user. Individual elements from these tracks are linked to the corresponding pages for additional details. For example, predicted gene models are linked to annotation reports, while genome conservation tracks are linked to VISTA Point (23). Two new types of data, single-nucleotide polymorphisms (SNPs) and gene expression profiles, have been added to the browser, to enable the visualization of resequencing and transcriptomics data, respectively (Figure 3). These are important additions to the resources developed, for example, for model fungi such as the ascomycete Neurospora crassa (7) or the basidiomycete Schizophyllum commune (24).

Predicted gene models, and their annotations, can be improved through the process of manual curation, which is open to all authorized users and coordinated by the principal investigator of each project. Such collective efforts were critical for the analysis of several fungal genomes (24–31). Additionally, users can upload custom tracks to the browser and share them, either with specified collaborators or globally. This capability to upload data for individual genes or entire genomes, makes MycoCosm a powerful and customizable data exchange resource for the community.

Grouping genomes for comparative analyses

The wealth of genomes in MycoCosm creates an opportunity to compare groups of genomes, and the need for efficient tools to enable these comparisons. Mycocosm therefore includes tools that integrate single genomes into a comparative context, such as the ability to visualize variation in gene counts in different GO, KEGG and KOG categories across a user-selected assortment of genomes. Additionally, MycoCosm includes predefined groups of fungi to enable users to ask questions for fungal genomes in a group context. Two types of genome groups are available in MycoCosm and shown in Figure 4: (i) PhyloGroups, consisting of phylogenetically related species (corresponding to the nodes of the MycoCosm tree), and (ii) EcoGroups, containing fungi of similar lifestyle or ecology (but which may be phylogenetically distant).

Both kinds of group offer consistent interfaces for comparative analyses. Genome sizes and gene counts may be compared on the Info page (Figure 4), giving a coarse-grained overview. Group-specific BLAST, search and download of genomes and annotations are available. Genes shared between genomes within a group can be explored using precomputed Markov clustering (32) of the full set of protein sequences in the group. Derived from all-versus-all BLAST alignments of the proteins, these gene clusters can be thought of as protein ‘families’. The clusters can be filtered to show genes shared between all, or some, organisms in the group, or to find genes unique to them as described in details in help pages. Each cluster is linked to a page listing all the genes of the cluster with graphical display of protein domains, intron–exon gene structure and local synteny of cluster members. A larger genome-scale synteny calculated by VISTA (23) is available from the Synteny menu of individual genome portals, where dot-plots, for example, can illustrate the full spectrum of micro-, meso- and macro-synteny between different pairs of Dothideomycete species (31). This phenomenon is limited to this class of fungi and illustrates the utility of the group-centric approach for comparative genome analysis.

User nomination of new species for sequencing

Even though all major nodes of the MycoCosm tree are represented by at least one sequenced genome, many fungal clades remain unrepresented. The 250+ annotated genomes in MycoCosm are a small fraction of the total number of fungal species, estimates of which range from 1.5 million to as high as 5.1 million (33,34). Thus, sampling of the Fungal Tree of Life by genomics is sparse, and likely has left many biologically interesting clades unexplored. This undersampling will surely limit our attempts to understand functional diversity and evolution through fungal genomes, and will hamper metagenomics studies of fungal communities, where a broad array of reference genomes is required to interpret complex sequence data.

To encourage the scientific community to participate in determining what gaps in the Fungal Tree of Life will be filled in next, every node of MycoCosm is now equipped with a ‘Nominate’ link to enable users to recommend new genomes for sequencing via the 1000 Fungal Genomes Project (http://1000.fungalgenomes.org). Following this link from a given node, families of fungi with no, or only one, sequenced genomes (Figure 5, green and yellow rows) display a ‘Nominate’ button, which invites users to describe the species, access JGI User Guidelines and provide DNA and RNA samples to JGI for sequencing (if approved by peer review). The 1000 Fungal Genomes Project aims to have at least two sequenced genomes for every family of fungi. Through this and other JGI Community Sequencing Projects (35), we aim to fill the remaining gaps in the Fungal Tree of Life.

To enable visibility of ongoing JGI genome projects, the Projects menu of MycoCosm (http://jgi.doe.gov/projects) lists each project’s status, principal investigator and other pertinent information. Combined with other genome projects listed in the GOLD database (36) and other sources, they form a ‘master’ list of fungal genome sequencing projects worldwide (http://1.usa.gov/JGI-1000-Fungi). This ‘master’ list helps us better focus our sequencing efforts on unexplored regions of the Fungal Tree of Life. Thus, we integrate a full genome project cycle (from sequencing target selection through assembly and annotation, to biology-driven data analysis) into MycoCosm enabling users to be active contributors to this process.

FUTURE PLANS

Since the release of its first fungal genome 10 years ago, the number of fungal genomes sequenced by JGI has jumped to >200. Collecting, integrating and visualizing these genomes, and those in sequencing pipeline, for efficient comparative analysis requires new approaches and tools. Especially, functional genomics will quickly reveal putative functions of a large number of genes identified only through genome sequencing (as opposed to smaller-scale biochemical and genetic studies). New interactive tools for transcriptomics, epigenomics, population genomics and metagenomics will therefore be developed and integrated into MycoCosm, in support of data analysis and utilization, and hypothesis development.

FUNDING

Director, Office of Science, Office of Biological and Environmental Research, Life Sciences Division, US Department of Energy [DE-AC02-05CH11231]. Funding for open access charge: US Department of Energy Contract [No. DE-AC02-05CH11231].

Conflict of interest statement. None declared.