Genome sequence of the insect pathogenic fungus Cordyceps militaris, a valued traditional chinese medicine
- Peng Zheng†1,
- Yongliang Xia†1,
- Guohua Xiao†1,
- Chenghui Xiong1,
- Xiao Hu1,
- Siwei Zhang1,
- Huajun Zheng2,
- Yin Huang2,
- Yan Zhou2,
- Shengyue Wang2,
- Guo-Ping Zhao1, 2,
- Xingzhong Liu3,
- Raymond J St Leger4 and
- Chengshu Wang1Email author
© Zheng et al.; licensee BioMed Central Ltd 2011
Received: 4 July 2011
Accepted: 23 November 2011
Published: 23 November 2011
Species in the ascomycete fungal genus Cordyceps have been proposed to be the teleomorphs of Metarhizium species. The latter have been widely used as insect biocontrol agents. Cordyceps species are highly prized for use in traditional Chinese medicines, but the genes responsible for biosynthesis of bioactive components, insect pathogenicity and the control of sexuality and fruiting have not been determined.
Here, we report the genome sequence of the type species Cordyceps militaris. Phylogenomic analysis suggests that different species in the Cordyceps/Metarhizium genera have evolved into insect pathogens independently of each other, and that their similar large secretomes and gene family expansions are due to convergent evolution. However, relative to other fungi, including Metarhizium spp., many protein families are reduced in C. militaris, which suggests a more restricted ecology. Consistent with its long track record of safe usage as a medicine, the Cordyceps genome does not contain genes for known human mycotoxins. We establish that C. militaris is sexually heterothallic but, very unusually, fruiting can occur without an opposite mating-type partner. Transcriptional profiling indicates that fruiting involves induction of the Zn2Cys6-type transcription factors and MAPK pathway; unlike other fungi, however, the PKA pathway is not activated.
The data offer a better understanding of Cordyceps biology and will facilitate the exploitation of medicinal compounds produced by the fungus.
The Ascomycete genus Cordyceps includes over 500 species that are pathogens of arthropods. Cordyceps militaris (CCM) is the type species and occurs throughout much of the Northern Hemisphere as a pathogen of lepidopteran insect pupae . C. militaris is readily characterized by the sexual fruiting bodies forming on mycosed pupae, the structures giving the fungus its common name of 'pupa grass' in China. Anamorphic Cordyceps species, such as Beauveria spp., Metarhizium spp. and Paecilomyces spp., have been developed as insect biocontrol agents [2, 3]. Although C. militaris and Cordyceps sinensis (syn. Ophiocordyceps sinensis) are best known as traditional Chinese medicines, they are also increasingly being studied and used in the West [4, 5]. An array of pharmacologically active components has been identified, including cordycepin, cordycepic acids, polysaccharides and macrolides . Cordycepin (3'-deoxyadenosine) has so far only been reported in C. militaris and is a broad spectrum antimicrobial  and polyadenylation inhibitor that is currently undergoing clinical trials against cancers . The biosynthetic pathway of cordycepin production has not been determined.
There is also much to be learnt about the nature and evolution of interactions of Cordyceps spp. with their hosts and with the wider environment. As entomopathogenicity appears to have evolved independently in Cordyceps and two Metarhizium species , comparative genomics will provide independent assessments of what is required to be entomopathogenic, identify the degree to which evolution between these fungi has been convergent, and identify the genomic basis of their differing physiologies and host-specificity. Last but not least, genomic sequencing of C. militaris will enable a systematic exploration of the biology and pharmaceuticals underlying the widespread medical impact of Cordyceps spp., and identify potential safety hazards, including genes for known human mycotoxins.
Genome sequencing and general features
Comparison of genome features among three insect pathogens
Percentage G+C content
Percentage repeat rate
Gene density (genes per Mb)
Exons per gene
Percentage secreted proteins
An InterproScan analysis identified 2,736 conserved protein families in C. militaris (containing 6,725 proteins), fewer than those in M. anisopliae (7,556 proteins in 2,796 families) or M. acridum (6,948 proteins in 2,746 families) . In particular, the number of transposases is much fewer in C. militaris (4) than in Metarhizium spp. (148 in M. anisopliae and 20 in M. acridum) or other sequenced ascomycetes (15 to 426) (Table S1 in Additional file 1). The C. militaris genome lacks retrotransposase (Table S2 in Additional file 1), and has more than three-fold fewer pseudogenes than Metarhizium spp. (Table S3 in Additional file 1). About 16% of the predicted C. militaris genes (1,547) are putatively involved in pathogen-host interactions; this proportion is slightly lower than for Metarhizium spp. (17.3% in MAA and 16.5% in MAC) but is higher than four plant pathogens (10.8 to 15.5%; P = 0.0476; false discovery rate (FDR) = 0.0152) (Table S4 in Additional file 1).
Genome-wide analysis of C. militaris gene sets.
Number of genes
Mean gene length (bp)
Mean number of introns per gene
Percentage genes without introns
Percentage GC content (excluding introns)
Number of InterproScan protein families
Number of secreted proteins
Number of PHI genesd
Number of TSA proteases
Number of MFS genes
Number of cytochrome P450s
Number of Pth11-like GPCRs
Number of protein kinases
Number of transcription factors
Number of glycoside hydrolases
Number of SM backbone genes
Number of horizontally transferred genes
Number of orthologs in M. anisopliae
Number of orthologs in M. acridum
Number of orthologs in F. graminearum
Number of orthologs in M. oryzae
Percentage identity to M. anisopliae orthologs
Percentage identity to M. acridum orthologs
Percentage identity to F. graminearum orthologs
Percentage identity to M. oryzae orthologs
The regions containing at least three contiguous open reading frames that are not present in the reference genome are designated as genomic islands (GIs) . Whole genome reciprocal analysis of three insect pathogens demonstrated that, in comparison to Metarhizium spp., C. militaris has 52 GIs (2% coverage of its genome, harboring 21% of its species-specific genes), which is many more than M. anisopliae (8 GIs, 0.3%) or M. acridum (5 GIs, 0.2%) when referenced to C. militaris. As in aspergilli , many C. militaris species-specific gene-encoding proteins do not have conserved domains and the genes are clustered together to form GIs (Table 2). A phylogenomic analysis established that the Cordyceps lineage is more closely related to the wheat pathogen F. graminearum (divergence time of 200 to 260 million years ago (MYA)) than it is to Metarhizium spp. (26 to 34 MYA) (Figure 3c). Thus, the lineage leading to C. militaris appears to have diverged from plant pathogens around the Triassic-Jurassic boundary (200 MYA), while M. anisopliae and M. acridum diverged after the Cretaceous Extinction Event (65 MYA) . Analysis of paralogous genes found only one pair of C. militaris genes with > 90% nucleotide sequence similarities (Figure 3d), which is similar to Neurospora crassa (one pair)  and F. graminearum (two pairs) . Analysis of 24 paired C. militaris genes showing > 70% nucleotide identities found a strong overall C:G to T:A mutation bias (Figure S1 in Additional file 2), consistent with repeat-induced point mutations, the DNA methylation-linked processes that cause mutations of repeated fungal sequences [15, 20].
Protein family analysis
We identified gene family expansions for proteases, chitinases, lipases and protein kinases in C. militaris when compared with phytopathogenic fungi, whereas gene family contractions occurred for glycoside hydrolases (GHs; P = 0.0144; FDR = 0.02), cutinases (P = 0.0065; FDR = 0.0226) and pectin lyases (P = 0.0245; FDR = 0.0284) (Table S1 in Additional file 1). The largest family expansions were for proteases. The C. militaris genome contains 61 families of proteases but most of them were included in families of serine proteases (180/381) and metallopeptidases (108/381) (Table S5 in Additional file 1). Gene expansions within the subtilisin (P = 0.0109; FDR = 0.0189) and trypsin (P = 0.0077; FDR = 0.0178) families are consistent with their being virulence factors in insect pathogens . However, different families of proteases are expanded in Metarhizium spp. and C. militaris, consistent with each lineage 'reinventing the wheel' during the evolution of entomopathogenicity. Thus, relative to Metarhizium spp., the S01 trypsin and S08 subtilisin subfamilies are smaller and the S53 subfamily is larger (Table S6 in Additional file 1). The C. militaris genome has 12 trypsin genes compared to 4 or less in plant pathogens. It lacks four subfamilies of trypsins present in M. anisopliae. Interestingly, the bacterial-like chymotrypsin identified in M. anisopliae  is absent in M. acridum  but present as two copies in C. militaris. The A01 aspartyl proteases are virulence factors of both mammalian and plant pathogens because of their ability to cleave an array of host proteins . Compared to phytopathogenic fungi (average 17), their number is significantly (P = 0.0059; FDR = 0.0057) expanded in the three insect pathogens (average 24) (Table S4 in Additional file 1).
Compared to many plant pathogens, Metarhizium spp. and C. militaris have fewer cutinases for degrading plant cell walls (Table S1 in Additional file 1). They also have fewer (average 137, P < 0.05) GHs than plant pathogens (average 199), including the lack of 20 GH families used by most plant pathogens and saprobes to target plant cell walls - for example, GH6, GH7 and GH61 cellulases, GH10 and GH11 xylanases, GH28 pectinases and GH78 rhamnosidases (Table S7 in Additional file 1). There are also significant differences in the spectrum of enzymes produced by the entomopathogens. For example, compared to M. anisopliae, C. militaris has few xyloglucosyl transferases (GH16) for xyloglucan catabolism and lacks α-glucuronidases (GH115) active on xylan oligomers or polymeric xylan . Consistent with this, C. militaris grows very poorly on xylose when compared with M. anisopliae (Figure S2 in Additional file 2). A phosphoketolase MPK1 involved in pentose metabolism is required for full virulence of M. anisopliae , but the homolog is absent in C. militaris. However, GH18 chitinases similar to those used by Metarhizium to degrade insect cuticles  are well represented in the C. militaris genome (20 in CCM versus 30 in MAA and 19 in MAC) relative to plant pathogens (average 11) (Table S7 in Additional file 1).
Cytochrome P450s (CYPs) play essential roles in fungal physiologies, including detoxification, degradation of xenobiotics and the biosynthesis of secondary metabolites . C. militaris has only about half as many CYPs as Metarhizium spp., and most other fungi (Table S8 in Additional file 1). Seventy CYP subfamilies present in M. anisopliae and/or M. acridum are absent in C. militaris. Of particular interest, C. militaris lacks CYP55, CYP58 and CYP65. CYP55 is a nitric oxide reductase required for denitrification . Thus, unlike most filamentous fungi, C. militaris may not respond to hypoxia through the bacterial ammonia fermentation mechanism. The absence of CYP58 (trichodiene oxygenase) and CYP65 (trichothecene C-15 hydroxylase) suggests that C. militaris will not produce the mycotoxin trichothecene . M. anisopliae can efficiently metabolize insect epicuticle alkanes . The CYP52 subfamily for alkane hydroxylation  is well represented in Cordyceps.
The major facilitator superfamily (MFS) and ATP-binding cassette (ABC) transporters are the two biggest families of fungal transporters. Members of the former typically function as nutrient symporters and drug antiporters, whereas the latter are more often implicated in defense against toxic metabolites . C. militaris has approximately half (123) the number of these transporters as Metarhizium (269 in MAA and 236 in MAC) (Table S9 in Additional file 1). The MFS transporters that are underrepresented in Cordyceps include the carbohydrate symporters (37 in CCM versus 48 in MAA, 51 in MAC and an average of 58 in plant pathogens), vitamin B2 (riboflavin) transporters (2 in CCM versus 17 each in Metarhizium species and an average of 4 in plant pathogens) and multidrug antiporters (23 in CCM versus 110 in MAA, 77 in MAC and an average of 10 in plant pathogens). Consistent with their having many multidrug transporters, Metarhizium spp. are resistant to diverse antibiotics and fungicides . Cordyceps has more ABC-type drug and metal resistant proteins than Metarhizium and plant pathogens (63 in CCM, 56 in MAA, 51 in MAC and an average of 54 in plant pathogens). The amino acid and dipeptide transporters are similarly represented in the three insect pathogens and other fungi (46 in CCM versus 53 in MAA, 49 in MAC and an average of 45 in plant pathogens).
Fungal G-protein coupled receptors (GPCRs) are required for pheromone/nutrient sensing and host recognition . Thus, the Pth11-like GPCR of Magnaporthe mediates cell differentiation in responses to plant inductive cues . C. militaris has fewer GPCRs than Metarhizium spp. and is particularly impoverished in Pth11-like GPCRs (Table S10 in Additional file 1). C. militaris has a similar number (167) of protein kinases as M. anisopliae (161) but less than M. acridum (192) (Table S11 in Additional file 1). Like other fungi, fungal specific transcription factors (TFs) and zinc finger TFs represent the two largest classes of TFs in C. militaris and their numbers are similar to those of other fungi (Table S1 in Additional file 1).
Mating-type and sexuality analysis
Metabolism of medically active components and mycotoxins
Numbers of core genes involved in the biosynthesis of secondary metabolites in different fungi
The mycotoxin ergot alkaloids have a wide range of biological activities and are important in pharmaceuticals and agriculture . Dimethylallyl tryptophan synthase (DMAT) catalyzes the alkylation of L-tryptophan, the first committed step in the ergot alkaloid biosynthetic pathway . C. militaris has one putative DMAT gene (CCM_04410), in contrast to five in M. anisopliae and three in M. acridum (Table 3). A phylogenetic analysis showed that CCM_04410 is not clustered with the Claviceps DMAT clade involved in ergot alkaloid production (Figure S5 in Additional file 2). The trichothecenes T-2 toxin and deoxynivalenol (type B trichothecene) are natural fungal products that are toxic to both animals and plants . Consistent with lacking CYP58 and CYP65, the C. militaris genome also lacks trichodiene synthase (Table S15 in Additional file 1). Thus, unlike Fusarium , C. militaris is not predicted to produce trichothecene mycotoxins. The presence of terpenoid cyclase, terpenoid synthase, fatty-acid synthase and geranylgeranyl diphosphate synthase genes in the C. militaris genome suggests that the fungus is capable of producing an array of metabolites, but the identity of these and their biological activities remain to be determined.
Transcriptional regulation of fruiting body development
Unlike other fungi, C. militaris can fruit sterilely in the absence of a sexual partner (Figure 5a-c). Perhaps because of this, 31 of the 42 C. militaris orthologs of sex-related genes identified in other ascomycetes were not expressed or transcribed at low levels (< 10 transcripts per million tags (TPM)) in sterile fruiting bodies (Table S18 in Additional file 1). However, in some cases, C. militaris expresses paralogous genes to those employed by other fungi, suggesting they have co-opted different components of the same signal transduction pathways to fulfill similar functions. For example, GATA-type TFs are important for fruiting in both A. nidulans and N. crassa , but C. militaris fruiting structures expressed orthologs of these genes at very low levels or not at all (Table S19 in Additional file 1). In contrast, the Zn2Cys6-type TFs were highly transcribed during fruiting but not in undifferentiated fungal mycelia - for example, CCM_01809 and CCM_09644 (Figure 9b) - indicating that Zn2Cys6 type TFs are predominately involved in the major developmental switch of production of fruiting structures.
We report here the first genome analysis of a Cordyceps species, the medicinal lepidopteran pathogen C. militaris, and show that the fungus is capable of fruiting without an opposite mating-type partner. We also show that it lacks genes known to be involved in production of human mycotoxins. Being an insect pathogen, the C. militaris genome contains thousands of genes putatively involved in interactions with insect hosts. Cordyceps resembles Metarhizium spp. in having a very high percentage of secreted proteins relative to plant pathogens and saprophytes and expanded families of proteases and chitinases with targets in insect hosts. However, insect-killing strategies may differ between Cordyceps and Metarhizium due to differences in gene content. Mating-type analysis indicated that sexual reproduction in C. militaris is heterothallic. Transcriptional profiling indicated that fruiting of the MAT1-1 C. militaris strain involves induction of the MAPK pathway, but unlike other homothallic or heterothallic fungi, the PKA pathway was not up-regulated. It remains to be determined whether this reflects the very unusual ability of C. militaris to produce fruiting bodies without an opposite mating-type partner.
Aside from knowing that C. militaris infects lepidopteran pupae , the life cycle of C. militaris in nature is poorly understood . Following disease, survival in soil may depend on the sexual stage of Cordyceps providing resilient long-lived ascospores as described in other fungi . Micro-cycle conidiation from germinated ascospores could adapt the fungus to nutrient poor niches (Figure 1). Metarhizium does not produce ascospores but flourishes in plant rhizospheres, which thus provide an alternative habitat in the absence of insect hosts. C. militaris can grow on germinated soybeans , suggesting a potential for an association with plants. However, relative to Metarhizium and most other ascomycetes, many protein families are smaller in the C. militaris genome, especially serine proteases, GHs, CYPs, MFS transporters and signal transduction factors. These families would be involved in scavenging for nutrients, avoidance of host defenses and toxins and other processes related to pathogenicity and a saprobic lifestyle. Around two-thirds of these protein families include pathogen-host interaction genes in plant-associated fungi (Table S4 in Additional file 1). Further studies on the ecology of Cordyceps spp. will shed more light on the relevance of the C. militaris genome to the evolution of gene families in relation to acquisition/loss of capability for dual plant/insect colonization and host range specialization.
The phylogenomic analysis demonstrated that the lineage leading to Cordyceps spp. diverged after most well known plant pathogens, including F. graminearum, but 130 MYA before Metarhizium diverged from the grass endophyte Epichloë festucae. The estimate of a Triassic-Jurassic boundary origin for the Cordyceps lineage and the post-Cretaceous origin of Metarhizium spp. is consistent with the hypocrealean fungi of Cordycipitaceae (includes Cordyceps spp.), Clavicipitaceae (includes Metarhizium spp.) and Ophiocordycipitaceae splitting about the same time as insects and angiosperms were diversifying . Families of proteases and chitinases are not expanded or lost in the E. festucae genome as they are in Cordyceps and Metarhizium, exemplifying convergent evolution to insect pathogenicity. Besides proteases and chitinases, experimentally verified Metarhizium virulence-associated genes with homologs in the C. militaris genome include a perilipin-like protein (CCM_06103 versus MAA_08819) to control cellular lipid storage and appressorium penetration , and an osmosensor (CCM_04885 versus MAA_01551) to mediate adaptation to the insect hemocoel . Homologs of these genes are broadly distributed among ascomycetes, indicative of an ancient origin. However, Cordyceps lacks other Metarhizium pathogenicity-related genes, including a collagen-like protein to evade the host immune system , a phosphoketolase for pentose metabolism  and the adhesins to mediate spore adhesions to insect and plant surfaces . The absence of key components of the Metarhizium entomopathogenicity 'toolkit' from C. militaris indicates that it has evolved different determinants to mediate its interactions with insects.
Like N. crassa and F. graminearum, C. miltaris lacks highly similar paralogs, a hallmark of the repeat-induced point mutation (RIP) mechanism . C. militaris has an ortholog (CCM_03609) of the N. crassa RIP defective gene (NCU02034), a cytosine methyltranserase essential for RIP . The high C→T and G→A mutation bias in the C. militaris genome and the readiness of C. miltaris to undergo the sexual cycle suggests that RIP is commonplace in C. militaris like many ascomycetes [10, 15, 19]. Since RIP can function effectively against selfish DNAs , it likely contributes, at least in part, to C. militaris having few DNA type transposon encoded genes, that is, transposases .
There are many more orphan genes in Cordyceps than in Metarhizium spp., underscoring that much about the proteome of Cordyceps spp. remains unknown. It is speculated that orphan genes arise from gene duplication, shuffling of gene fragments, mobile element effects, mutation of existing sequences, horizontal gene transfer and de novo origination from non-coding DNAs . De novo creation of new genes is probably rare . A role for mobile element effects is also unlikely given how few putative transposase genes are present in the C. militaris genome. Putative horizontal gene transfer genes are even fewer in Cordyceps than in Metarhizium. Thus, the numerous orphans in the C. militaris genome most likely arose from frequent mutations caused by RIP in existing (duplicated) sequences. Just as the Metarhizium-specific collagen-like protein is essentially required to camouflage cells from host immune recognition . The transcriptome data showed that 428 of the 1,329 orphan genes were transcribed during fruiting. Of the 100 most highly expressed genes in developing C. militaris fruiting bodies, about one-third are orphans (Table S17 in Additional file 1), underscoring the potential of orphans to have specific functions. Likewise, genes that are apparently unique to the mushroom Schizophyllum commune are more likely to be expressed during mushroom formation .
Concern has been raised about the possibility of harmful side effects of traditional Chinese medicines, including Cordyceps . Consistent with genotoxicity and cytotoxicity assays that show Cordyceps products to be safe for consumption , there is no evidence in the C. militaris genome for genes involved in the production of known mycotoxins. However, safety could only be completely verified by meticulous profiling of the metabolites produced by the fungus under diverse growth conditions. The C. militaris genome data will facilitate these processes as well as help with elucidation of the biosynthetic pathways of different metabolites.
The analysis of C. militaris genome indicates that it is sexually heterothallic, but strikingly, both the MAT1-1 single mating-type and MAT1-1/MAT1-2 hybrid strains can form fruiting bodies, which means that C. militaris is capable of fruiting without a partner. Single mating-type (haploid) fruiting has also been observed in the human pathogens Cryptococcus neoformans and Candida albicans . Given that perithecia and ascospores are not produced by MAT1-1 fruiting bodies, C. militaris haploid fruiting is different from the same-sex mating and fruiting of C. neoformans, in which diploidization and meiosis can occur. In budding yeast, meiotic recombination is initiated by the formation of double-strand breaks catalyzed by SPO11, a meiosis-specific endonuclease . The meiosis-specific recombinase DMC1 and the DNA repair enzyme RAD51 then co-localize to double-strand breaks and function together for meiotic recombination . The C. militaris homologue of yeast SPO11 (CCM_09527) was up-regulated more than five-fold during fruiting body maturation (the TPM ratio of FB2/FB1 = 8.1; FB3/FB2 = 5.7). Intriguingly, the C. militaris genome lacks a yeast RAD51 ortholog, but its DMC1 ortholog (CCM_06822) contains a RAD51 domain. CCM_06822 was not expressed by C. militaris during fruiting, which may explain why the C. militaris MAT1-1 strain forms fruiting bodies without meiosis. Consistent with this, a putative cyclin dependent kinase 7 (CDK7; CCM_03900) was up-regulated during fungal fruiting (Table S19 in Additional file 1). Orthologs of CDK7 initiate DNA synthesis and facilitate mitosis instead of meiosis .
In conclusion, we report on the genome sequencing, comparative genome analysis and transcriptional regulation of fruiting body development in the medicinal fungus C. militaris. The sequence data should markedly enhance the pace of molecular research on Cordyceps biology, fungal sex and pathogenicity, and will have impacts on the commercial production of fruiting structures. The genomic sequence will also be an essential tool to unravel the mechanisms by which C. militaris produces medicinal compounds and so further their exploitation.
Materials and methods
C. militaris strain Cm01 (CGMCC 3.14242) was selected for genome sequencing as it is culturally stable and commercialized in China. The culture was maintained either on artificial medium or silkworm pupae as previously described . Several different C. militaris strains were included in this study for PCR genotyping of mating-type genes (Figure 5h).
Genome sequencing and assembly
The genome of C. militaris strain Cm01 was shotgun sequenced using a Roche 454 GS FLX system for massively parallel pyrosequencing for 2.25 runs at the Chinese National Human Genome Center (Shanghai, China). This resulted in 951 Mb of sequence data (29.6 × coverage) with an average read length of 385 bp. Assembly was performed using the Newbler software (v2.3) within the Roche 454 suite package , which produced 597 contigs with a total size of 32.2 Mb. For sequence scaffolding, a DNA library of 2- to 5-kb inserts was generated and sequenced with an ABI SOLiD system (Carlsbad, California, USA). This resulted in 3.8 Gb of mate-pair reads (117.4 × coverage) to improve sequence quality and construct scaffolds. By mapping the reads to contigs, 578 contigs were assembled into 13 scaffolds and 19 contigs less than 2 kb left outside. The raw data of 454 and SOLiD reads have been deposited at NCBI's Sequence Read Archive under accession number SRA047932 and the whole project has been deposited at DDBJ/EMBL/GenBank under accession number AEVU00000000.
To maximize gene prediction accuracy, the gene structures of Cordyceps were predicted with a combination of different algorithms plus manual inspections [11, 57]. The inconsistent open reading frames were individually subject to Blast searches against the NCBI curated refseq_protein database. The prediction with the best hit was selected. Pseudogene identification was conducted with the pipeline of PseudoPipe with default settings . The potential secreted proteins of C. militaris and other fungal species included for comparison were predicted by SignaIP 3.0 analysis using a hidden Markov model . Genome repetitive elements were analyzed by Blast against the RepeatMasker library (Open 3.2.9)  and with the Tandem Repeats Finder . The transposases/retrotransposases were classified by Blastp analysis against the Repbase  plus manual inspections.
Blast score ratio test
BSR tests  were conducted to compare the differences between C. militaris and the sequenced Metarhizum genomes and the plant pathogen F. graminearum, respectively. The BSR index for each reference protein is calculated by dividing the query bit score by the reference score and normalized from 0 to 1. A score of 1 indicates a perfect match while a score of 0 indicates no Blast match of a query protein in the reference proteome. The normalized pairs of BSR indices were then plotted using the Matlab (v7.0) program (Natick, Massachusetts, USA). The same analysis was conducted for the three genomes of C. militaris, M. anisopliae and F. graminearum.
Orthology and phylogenomic analysis
In total, 2,106 orthologous proteins were acquired by a reciprocal Blast method with a cutoff E value of 1e-20 and a Blast alignment length greater than 60% of the query sequence. Corresponding orthologous gene protein sequences were aligned with Clustal X 2.0 and the concatenated amino acid sequences were used for the generation of a maximum likelihood phylogenomic tree with the program TREE-PUZZLE  using a Dayhoff model. The divergence time between species was estimated with the program r8s using a Langley-Fitch model  by calibration with the origin of the Ascomycota at 500 to 650 MYA .
Protein family classifications
Whole genome protein families were classified by InterproScan  and Pfam  analysis. The families of proteases were identified by Blastp searching against the MEROPS peptidase database release 9.4 with a cutoff E value of 1e-20 . The CYPs were named according to the classifications collected at the P450 database . Transporters were classified based on the Transport Classification Database . Kinases were classified by Blastp analysis against the KinBase database with a cutoff E value of 1e-10 . Carbohydrate-active enzymes were classified by local Blastp searching against a library of catalytic and carbohydrate-binding module enzymes . G-protein-coupled receptors were selected from the best hits to GPCRDB sequences  and by confirmation that they contained seven transmembrane helices with the amino terminus outside and the carboxyl terminus inside the plasma membrane. Homologs of the Magnaporthe Pth11-like GPCRs  were identified by local Blastp analysis with a cutoff E value of 1e-10. Putative Cordyceps virulence factors were identified by searching against the pathogen-host interaction database  with a cutoff E value of 1e-5, plus additional searches of known virulence genes reported in entomopathogenic fungi. Two sample t-tests were conducted to compare the differences in protein family sizes between insect and plant pathogens. Estimation of FDR of P-values was conducted using the program mafdr (Matlab 220.127.116.117(R2009a)).
Analysis of genes involved in purine synthesis and secondary metabolism
To model the biosynthesis of cordycepin, the purine metabolic pathway in C. militaris was constructed based on the KEGG (Kyoto Encyclopedia of Genes and Genomes) annotations . To identify NRPS, PKS or NRPS-PKS hybrid genes and gene clusters, the whole genome data set was subjected to analysis with the program SMURF with default settings . Modulation analysis and domain extraction of different NRPS or PKS proteins were conducted by Blast searching against the SBSPKS database . For phylogenetic analysis, the domain sequences were aligned with Clustal X 2.0 and the tree was generated using a Poisson model with 1,000 bootstrap replications and pair-wise deletions for gaps or missing data. The mycotoxin-encoding PKSs used in the analysis include Gibberella zeae PKS4 (ABB90283) and PKS13 (ABB90282) for the biosynthesis of zearalenones , Aspergillus ochraceus PKS (AAT92023) for ochratoxin , A. clavatus PatK (ACLA_093660) for patulin , Gibberella moniliformis Fum1p (AAD43562) for fumonisin , Monascus purpureus PKS (BAD44749) for citrinin , Aspergillus flavus PksA (AAS90093) for aflatoxin  and A. nidulans StcA (Q12397) for sterigmatocystin . The mycotoxin-encoding NRPSs included in the analysis are A. fumigatus Glip (AAW03307) for gliotoxin , Fusarium equiseti NRPS (CAA79245) for enniatin , Cochliobolus carbonum NRPS (AAA33023) for HC-toxin  and Tolypocladium inflatum NRPS (CAA82227) for cyclosporin .
Conidia of C. militaris from day 14 potato dextrose agar were inoculated into SDB and the undifferentiated mycelia harvested after a 72-hour incubation at 25°C, 180 rpm. The transcriptome of the mycelia provided a control for comparison with the transcriptomes of fruiting bodies. Chinese Tussah silkmoth (Antheraea pernyi) pupae were injected with 50 μl of a conidial suspension (5 × 106 conidia/ml) and incubated at 22°C in a 12:12 hour light:dark cycle for up to 14 days to allow emergence of fruiting bodies (a stage designated as FB1), 29 days for half-grown fruiting bodies (FB2) and 59 days, by which time fruiting bodies were mature (FB3) [11, 12]. RNA was extracted with a Qiagen RNeasy kit plus on-column treatment with RNase-free DNase I (Germantown, Maryland, USA). Messenger RNA was purified and after reverse transcription into cDNA, the libraries were constructed for tag preparation according to the massively parallel signature sequencing protocol . The tags were sequenced with an Illumina technique. We omitted tags from further analysis if only one copy was detected or it could be mapped to a different transcript. Other tags were mapped to the genome or annotated genes if they possessed no more than one nucleotide mismatch [11, 88]. The abundance of each tag was converted to the value of transcripts per million (TPM) for each mapped gene for expressional comparison between samples. The significance of differential gene expression between samples and the FDR of P-values were estimated for each individual gene with a cutoff of P ≤ 0.05 and FDR ≤0.001 [11, 89]. The RNA_seq expression dataset is available at the NCBI's Gene Expression Omnibus under the accession code GSE28001.
Blast score ratio
dimethylallyl tryptophan synthase
false discovery rate
G-protein coupled receptor
mitogen-activated protein kinase
major facilitator superfamily
million years ago
non-ribosomal peptide synthetase
protein kinase A
repeat-induced point mutation
ribonucleotide trisphosphate reductase
Sabouraud dextrose broth
transcripts per million tags.
CSW was supported by the Ministry of Science and Technology of China (grant number 2009CB118904), the Ministry of Agriculture of China (grant number 2009ZX08009-035B), the Science and Technology Commission of Shanghai Municipality (grant number 08DZ1970200) and the Chinese Academy of Sciences (KSCX2-EW-N-06 and KSCX2-EW-G-16).
- Sung GH, Hywel-Jones NL, Sung JM, Luangsa-Ard JJ, Shrestha B, Spatafora JW: Phylogenetic classification of Cordyceps and the clavicipitaceous fungi. Stud Mycol. 2007, 57: 5-59. 10.3114/sim.2007.57.01.PubMedPubMed CentralView ArticleGoogle Scholar
- de Faria MR, Wraight SP: Mycoinsecticides and Mycoacaricides: a comprehensive list with worldwide coverage and international classification of formulation types. Biol Control. 2007, 43: 237-256. 10.1016/j.biocontrol.2007.08.001.View ArticleGoogle Scholar
- St Leger RJ, Wang C: Genetic engineering of fungal biocontrol agents to achieve greater efficacy against insect pests. Appl Microbiol Biotechnol. 2010, 85: 901-907. 10.1007/s00253-009-2306-z.PubMedView ArticleGoogle Scholar
- Paterson RR: Cordyceps: a traditional Chinese medicine and another fungal therapeutic biofactory?. Phytochemistry. 2008, 69: 1469-1495. 10.1016/j.phytochem.2008.01.027.PubMedView ArticleGoogle Scholar
- Zhou X, Gong Z, Su Y, Lin J, Tang K: Cordyceps fungi: natural products, pharmacological functions and developmental products. J Pharm Pharmacol. 2009, 61: 279-291.PubMedView ArticleGoogle Scholar
- Xiao JH, Zhong JJ: Secondary metabolites from Cordyceps species and their antitumor activity studies. Recent Pat Biotechnol. 2007, 1: 123-137. 10.2174/187220807780809454.PubMedView ArticleGoogle Scholar
- ClinicalTrials.gov. [http://clinicaltrials.gov/show/NCT00709215]
- Stone R: Last stand for the body snatcher of the Himalayas?. Science. 2008, 322: 1182-10.1126/science.322.5905.1182.PubMedView ArticleGoogle Scholar
- Lee SC, Ni M, Li W, Shertz C, Heitman J: The evolution of sex: a perspective from the fungal kingdom. Microbiol Mol Biol Rev. 2010, 74: 298-340. 10.1128/MMBR.00005-10.PubMedPubMed CentralView ArticleGoogle Scholar
- Galagan JE, Calvo SE, Cuomo C, Ma LJ, Wortman JR, Batzoglou S, Lee SI, Bastürkmen M, Spevak CC, Clutterbuck J, Kapitonov V, Jurka J, Scazzocchio C, Farman M, Butler J, Purcell S, Harris S, Braus GH, Draht O, Busch S, D'Enfert C, Bouchier C, Goldman GH, Bell-Pedersen D, Griffiths-Jones S, Doonan JH, Yu J, Vienken K, Pain A, Freitag M, et al: Sequencing of Aspergillus nidulans and comparative analysis with A. fumigatus and A. oryzae. Nature. 2005, 438: 1105-1115. 10.1038/nature04341.PubMedView ArticleGoogle Scholar
- Gao Q, Jin K, Ying SH, Zhang Y, Xiao G, Shang Y, Duan Z, Hu X, Xie XQ, Zhou G, Peng G, Luo Z, Huang W, Wang B, Fang W, Wang S, Zhong Y, Ma LJ, St Leger RJ, Zhao GP, Pei Y, Feng MG, Xia Y, Wang CS: Genome sequencing and comparative transcriptomics of the model entomopathogenic fungi Metarhizium anisopliae and M. acridum. PLoS Genet. 2011, 7: e1001264-10.1371/journal.pgen.1001264.PubMedPubMed CentralView ArticleGoogle Scholar
- Xiong CH, Xia YL, Zheng P, Shi SH, Wang CS: Developmental stage-specific gene expression profiling for a medicinal fungus Cordyceps militaris. Mycology. 2010, 1: 25-66.View ArticleGoogle Scholar
- Galagan JE, Calvo SE, Borkovich KA, Selker EU, Read ND, Jaffe D, FitzHugh W, Ma LJ, Smirnov S, Purcell S, Rehman B, Elkins T, Engels R, Wang S, Nielsen CB, Butler J, Endrizzi M, Qui D, Ianakiev P, Bell-Pedersen D, Nelson MA, Werner-Washburne M, Selitrennikoff CP, Kinsey JA, Braun EL, Zelter A, Schulte U, Kothe GO, Jedd G, Mewes W, et al: The genome sequence of the filamentous fungus Neurospora crassa. Nature. 2003, 422: 859-868. 10.1038/nature01554.PubMedView ArticleGoogle Scholar
- Ma LJ, van der Does HC, Borkovich KA, Coleman JJ, Daboussi MJ, Di Pietro A, Dufresne M, Freitag M, Grabherr M, Henrissat B, Houterman PM, Kang S, Shim WB, Woloshuk C, Xie X, Xu JR, Antoniw J, Baker SE, Bluhm BH, Breakspear A, Brown DW, Butchko RA, Chapman S, Coulson R, Coutinho PM, Danchin EG, Diener A, Gale LR, Gardiner DM, Goff S, et al: Comparative genomics reveals mobile pathogenicity chromosomes in Fusarium. Nature. 2010, 464: 367-373. 10.1038/nature08850.PubMedPubMed CentralView ArticleGoogle Scholar
- Kubicek CP, Herrera-Estrella A, Seidl-Seiboth V, Martinez DA, Druzhinina IS, Thon M, Zeilinger S, Casas-Flores S, Horwitz BA, Mukherjee PK, Mukherjee M, Kredics L, Alcaraz LD, Aerts A, Antal Z, Atanasova L, Cervantes-Badillo MG, Challacombe J, Chertkov O, McCluskey K, Coulpier F, Deshpande N, von Döhren H, Ebbole DJ, Esquivel-Naranjo EU, Fekete E, Flipphi M, Glaser F, Gómez-Rodríguez EY, Gruber S, et al: Comparative genome sequence analysis underscores mycoparasitism as the ancestral life style of Trichoderma. Genome Biol. 2011, 12: R40-10.1186/gb-2011-12-4-r40.PubMedPubMed CentralView ArticleGoogle Scholar
- Rajashekara G, Glasner JD, Glover DA, Splitter GA: Comparative whole-genome hybridization reveals genomic islands in Brucella species. J Bacteriol. 2004, 186: 5040-5051. 10.1128/JB.186.15.5040-5051.2004.PubMedPubMed CentralView ArticleGoogle Scholar
- Fedorova ND, Khaldi N, Joardar VS, Maiti R, Amedeo P, Anderson MJ, Crabtree J, Silva JC, Badger JH, Albarraq A, Angiuoli S, Bussey H, Bowyer P, Cotty PJ, Dyer PS, Egan A, Galens K, Fraser-Liggett CM, Haas BJ, Inman JM, Kent R, Lemieux S, Malavazi I, Orvis J, Roemer T, Ronning CM, Sundaram JP, Sutton G, Turner G, Venter JC, et al: Genomic islands in the pathogenic filamentous fungus Aspergillus fumigatus. PLoS Genet. 2008, 4: e1000046-10.1371/journal.pgen.1000046.PubMedPubMed CentralView ArticleGoogle Scholar
- McElwain JC, Punyasena SW: Mass extinction events and the plant fossil record. Trends Ecol Evol. 2007, 22: 548-557. 10.1016/j.tree.2007.09.003.PubMedView ArticleGoogle Scholar
- Cuomo CA, Güldener U, Xu JR, Trail F, Turgeon BG, Di Pietro A, Walton JD, Ma LJ, Baker SE, Rep M, Adam G, Antoniw J, Baldwin T, Calvo S, Chang YL, Decaprio D, Gale LR, Gnerre S, Goswami RS, Hammond-Kosack K, Harris LJ, Hilburn K, Kennell JC, Kroken S, Magnuson JK, Mannhaupt G, Mauceli E, Mewes HW, Mitterbauer R, Muehlbauer G, et al: The Fusarium graminearum genome reveals a link between localized polymorphism and pathogen specialization. Science. 2007, 317: 1400-1402. 10.1126/science.1143708.PubMedView ArticleGoogle Scholar
- Galagan JE, Selker EU: RIP: The evolutionary cost of genome defense. Trends Genet. 2004, 20: 417-423. 10.1016/j.tig.2004.07.007.PubMedView ArticleGoogle Scholar
- Screen SE, St Leger RJ: Cloning, expression, and substrate specificity of a fungal chymotrypsin. Evidence for lateral gene transfer from an actinomycete bacterium. J Biol Chem. 2000, 275: 6689-6694. 10.1074/jbc.275.9.6689.PubMedView ArticleGoogle Scholar
- Coetzer TH, Goldring JP, Huson LE: Oligopeptidase B: a processing peptidase involved in pathogenesis. Biochimie. 2008, 90: 336-344. 10.1016/j.biochi.2007.10.011.PubMedView ArticleGoogle Scholar
- Chong SL, Battaglia E, Coutinho PM, Henrissat B, Tenkanen M, de Vries RP: The α-glucuronidase Agu1 from Schizophyllum commune is a member of a novel glycoside hydrolase family (GH115). Appl Microbiol Biotechnol. 2011, 90: 1323-1332. 10.1007/s00253-011-3157-y.PubMedView ArticleGoogle Scholar
- Duan ZB, Shang Y, Gao Q, Zheng P, Wang CS: A phosphoketolase Mpk1 of bacterial origin is adaptively required for full virulence in the insect-pathogenic fungus Metarhizium anisopliae. Environ Microbiol. 2009, 11: 2351-2360. 10.1111/j.1462-2920.2009.01961.x.PubMedView ArticleGoogle Scholar
- Crešnar B, Petrič S: Cytochrome P450 enzymes in the fungal kingdom. Biochim Biophys Acta. 2011, 1814: 29-35.PubMedView ArticleGoogle Scholar
- Keller NP, Turner G, Bennett JW: Fungal secondary metabolism - from biochemistry to genomics. Nat Rev Microbiol. 2005, 3: 937-947. 10.1038/nrmicro1286.PubMedView ArticleGoogle Scholar
- Jarrold SL, Moore D, Potter U, Charnley AK: The contribution of surface waxes to pre-penetration growth of an entomopathogenic fungus on host cuticle. Mycol Res. 2007, 111: 240-249. 10.1016/j.mycres.2006.10.007.PubMedView ArticleGoogle Scholar
- Morschhauser J: Regulation of multidrug resistance in pathogenic fungi. Fungal Genet Biol. 2010, 47: 94-106. 10.1016/j.fgb.2009.08.002.PubMedView ArticleGoogle Scholar
- Luz C, Netto MC, Rocha LF: In vitro susceptibility to fungicides by invertebrate-pathogenic and saprobic fungi. Mycopathologia. 2007, 164: 39-47. 10.1007/s11046-007-9020-0.PubMedView ArticleGoogle Scholar
- Kulkarni RD, Thon MR, Pan H, Dean RA: Novel G-protein-coupled receptor-like proteins in the plant pathogenic fungus Magnaporthe grisea. Genome Biol. 2005, 6: R24-10.1186/gb-2005-6-3-r24.PubMedPubMed CentralView ArticleGoogle Scholar
- Hejtmanek M, Hejtmankova N: Teleomorphs and mating types in Trichophyton mentagrophytes complex. Acta Univ Palacki Olomuc Fac Med. 1989, 123: 11-33.PubMedGoogle Scholar
- Lennon MB, Suhadolnik RJ: Biosynthesis of 3'-deoxyadenosine by Cordyceps militaris. Mechanism of reduction. Biochim Biophys Acta. 1976, 425: 532-536.PubMedView ArticleGoogle Scholar
- Reichard P: Ribonucleotide reductases: substrate specificity by allostery. Biochem Biophys Res Commun. 2010, 396: 19-23. 10.1016/j.bbrc.2010.02.108.PubMedView ArticleGoogle Scholar
- Liu M, Panaccione DG, Schardl CL: Phylogenetic analyses reveal monophyletic origin of the ergot alkaloid gene dmaW in fungi. Evol Bioinform. 2009, 5: 15-30.Google Scholar
- Brown DW, McCormick SP, Alexander NJ, Proctor RH, Desjardins AE: A genetic and biochemical approach to study trichothecene diversity in Fusarium sporotrichioides and Fusarium graminearum. Fungal Genet Biol. 2001, 32: 121-133. 10.1006/fgbi.2001.1256.PubMedView ArticleGoogle Scholar
- Pöggeler S, Nowrousian M, Ringelberg C, Loros JJ, Dunlap JC, Kuck U: Microarray and real-time PCR analyses reveal mating type-dependent gene expression in a homothallic fungus. Mol Genet Genomics. 2006, 275: 492-503. 10.1007/s00438-006-0107-y.PubMedView ArticleGoogle Scholar
- Kamata N, Sato H, Shimazu M: Seasonal changes in the infection of pupae of the beech caterpillar, Quadricalcarifera punctatella (Motsch.) (Lep., Notodontidae), by Cordyceps militaris Link (Clavicipitales, Clavicipitaceae) in the soil of the Japanese beech forest. J Appl Entomol. 1997, 121: 17-21. 10.1111/j.1439-0418.1997.tb01364.x.View ArticleGoogle Scholar
- Clarkson JP, Staveley J, Phelps K, Young CS, Whipps JM: Ascospore release and survival in Sclerotinia sclerotiorum. Mycol Res. 2003, 107: 213-222. 10.1017/S0953756203007159.PubMedView ArticleGoogle Scholar
- Ohta Y, Lee JB, Hayashi K, Fujita A, Park DK, Hayashi T: In vivo anti-influenza virus activity of an immunomodulatory acidic polysaccharide isolated from Cordyceps militaris grown on germinated soybeans. J Agric Food Chem. 2007, 55: 10194-10199. 10.1021/jf0721287.PubMedView ArticleGoogle Scholar
- Sung GH, Poinar GO, Spatafora JW: The oldest fossil evidence of animal parasitism by fungi supports a Cretaceous diversification of fungal-arthropod symbioses. Mol Phylogenet Evol. 2008, 49: 495-502. 10.1016/j.ympev.2008.08.028.PubMedView ArticleGoogle Scholar
- Wang CS, St Leger RJ: The Metarhizium anisopliae perilipin homolog MPL1 regulates lipid metabolism, appressorial turgor pressure, and virulence. J Biol Chem. 2007, 282: 21110-21115. 10.1074/jbc.M609592200.PubMedView ArticleGoogle Scholar
- Wang CS, Duan ZB, St Leger RJ: MOS1 osmosensor of Metarhizium anisopliae is required for adaptation to insect host hemolymph. Eukaryot Cell. 2008, 7: 302-309. 10.1128/EC.00310-07.PubMedPubMed CentralView ArticleGoogle Scholar
- Wang CS, St Leger RJ: A collagenous protective coat enables Metarhizium anisopliae to evade insect immune responses. Proc Natl Acad Sci USA. 2006, 103: 6647-6652. 10.1073/pnas.0601951103.PubMedPubMed CentralView ArticleGoogle Scholar
- Wang CS, St Leger RJ: The MAD1 adhesin of Metarhizium anisopliae links adhesion with blastospore production and virulence to insects, and the MAD2 adhesin enables attachment to plants. Eukaryot Cell. 2007, 6: 808-816. 10.1128/EC.00409-06.PubMedPubMed CentralView ArticleGoogle Scholar
- Freitag M, Williams RL, Kothe GO, Selker EU: A cytosine methyltransferase homologue is essential for repeat-induced point mutation in Neurospora crassa. Proc Natl Acad Sci USA. 2002, 99: 8802-8807. 10.1073/pnas.132212899.PubMedPubMed CentralView ArticleGoogle Scholar
- Montiel MD, Lee HA, Archer DB: Evidence of RIP (repeat-induced point mutation) in transposase sequences of Aspergillus oryzae. Fungal Genet Biol. 2006, 43: 439-445. 10.1016/j.fgb.2006.01.011.PubMedView ArticleGoogle Scholar
- Long M, Betran E, Thornton K, Wang W: The origin of new genes: glimpses from the young and old. Nat Rev Genet. 2003, 4: 865-875.PubMedView ArticleGoogle Scholar
- Ekman D, Elofsson A: Identifying and quantifying orphan protein sequences in fungi. J Mol Biol. 2010, 396: 396-405. 10.1016/j.jmb.2009.11.053.PubMedView ArticleGoogle Scholar
- Ohm RA, de Jong JF, Lugones LG, Aerts A, Kothe E, Stajich JE, de Vries RP, Record E, Levasseur A, Baker SE, Bartholomew KA, Coutinho PM, Erdmann S, Fowler TJ, Gathman AC, Lombard V, Henrissat B, Knabe N, Kües U, Lilly WW, Lindquist E, Lucas S, Magnuson JK, Piumi F, Raudaskoski M, Salamov A, Schmutz J, Schwarze FW, vanKuyk PA, Horton JS, Grigoriev IV, et al: Genome sequence of the model mushroom Schizophyllum commune. Nat Biotechnol. 2010, 28: 957-963. 10.1038/nbt.1643.PubMedView ArticleGoogle Scholar
- Corson TW, Crews CM: Molecular understanding and modern application of traditional medicines: triumphs and trials. Cell. 2007, 130: 769-774. 10.1016/j.cell.2007.08.021.PubMedPubMed CentralView ArticleGoogle Scholar
- Yan WJ, Li TH, Lin QY, Song B, Jiang ZD: Safety assessment of Cordyceps guangdongensis. Food Chem Toxicol. 2010, 48: 3080-3084. 10.1016/j.fct.2010.08.001.PubMedView ArticleGoogle Scholar
- Heitman J: Evolution of eukaryotic microbial pathogens via covert sexual reproduction. Cell Host Microbe. 2010, 8: 86-99. 10.1016/j.chom.2010.06.011.PubMedPubMed CentralView ArticleGoogle Scholar
- Keeney S: Mechanism and control of meiotic recombination initiation. Curr Top Dev Biol. 2001, 52: 1-53.PubMedView ArticleGoogle Scholar
- Masson JY, West SC: The Rad51 and Dmc1 recombinases: a non-identical twin relationship. Trends Biochem Sci. 2001, 26: 131-136. 10.1016/S0968-0004(00)01742-4.PubMedView ArticleGoogle Scholar
- Fisher RP: Secrets of a double agent: CDK7 in cell-cycle control and transcription. J Cell Sci. 2005, 118: 5171-5180. 10.1242/jcs.02718.PubMedView ArticleGoogle Scholar
- Schatz MC, Delcher AL, Salzberg SL: Assembly of large genomes using second-generation sequencing. Genome Res. 2010, 20: 1165-1173. 10.1101/gr.101360.109.PubMedPubMed CentralView ArticleGoogle Scholar
- Ter-Hovhannisyan V, Lomsadze A, Chernoff YO, Borodovsky M: Gene prediction in novel fungal genomes using an ab initio algorithm with unsupervised training. Genome Res. 2008, 18: 1979-1990. 10.1101/gr.081612.108.PubMedPubMed CentralView ArticleGoogle Scholar
- Zhang Z, Carriero N, Zheng D, Karro J, Harrison PM, Gerstein M: PseudoPipe: an automated pseudogene identification pipeline. Bioinformatics. 2006, 22: 1437-1439. 10.1093/bioinformatics/btl116.PubMedView ArticleGoogle Scholar
- SignalP 3.0 Server. [http://www.cbs.dtu.dk/services/SignalP/]
- RepeastMasker Server. [http://www.repeatmasker.org/]
- Benson G: Tandem repeats finder: a program to analyze DNA sequences. Nucleic Acids Res. 1999, 27: 573-580. 10.1093/nar/27.2.573.PubMedPubMed CentralView ArticleGoogle Scholar
- Kohany O, Gentles AJ, Hankus L, Jurka J: Annotation, submission and screening of repetitive elements in Repbase: RepbaseSubmitter and Censor. BMC Bioinformatics. 2006, 7: 474-10.1186/1471-2105-7-474.PubMedPubMed CentralView ArticleGoogle Scholar
- Schmidt HA, Strimmer K, Vingron M, von Haeseler A: TREE-PUZZLE: maximum likelihood phylogenetic analysis using quartets and parallel computing. Bioinformatics. 2002, 18: 502-504. 10.1093/bioinformatics/18.3.502.PubMedView ArticleGoogle Scholar
- Taylor JW, Berbee M: Dating divergences in the Fungal Tree of Life: review and new analyses. Mycologia. 2006, 98: 838-849. 10.3852/mycologia.98.6.838.PubMedView ArticleGoogle Scholar
- Lucking R, Huhndorf S, Pfister DH, Plata ER, Lumbsch HT: Fungi evolved right on track. Mycologia. 2009, 101: 810-822. 10.3852/09-016.PubMedView ArticleGoogle Scholar
- InterPro Protein Sequence Analysis & Classification. [http://www.ebi.ac.uk/interpro/]
- Pfam Database. [http://pfam.sanger.ac.uk/]
- Schomburg D, Schomburg IL: Enzyme databases. Methods Mol Biol. 2010, 609: 113-128. 10.1007/978-1-60327-241-4_7.PubMedView ArticleGoogle Scholar
- Cytochrome P450 Server. [http://drnelson.uthsc.edu/CytochromeP450.html]
- Saier MH, Tran CV, Barabote RD: TCDB: the Transporter Classification Database for membrane transport protein analyses and information. Nucleic Acids Res. 2006, 34: D181-186. 10.1093/nar/gkj001.PubMedPubMed CentralView ArticleGoogle Scholar
- GPCR Database. [http://www.gpcr.org/7tm/]
- PHI-base: Pathogen-Host Interactions. [http://www.phi-base.org]
- Kanehisa M, Goto S, Kawashima S, Nakaya A: The KEGG databases at GenomeNet. Nucleic Acids Res. 2002, 30: 42-46. 10.1093/nar/30.1.42.PubMedPubMed CentralView ArticleGoogle Scholar
- SMURF. [http://www.jcvi.org/smurf/index.php]
- Anand S, Prasad MV, Yadav G, Kumar N, Shehara J, Ansari MZ, Mohanty D: SBSPKS: structure based sequence analysis of polyketide synthases. Nucleic Acids Res. 2010, 38: W487-496. 10.1093/nar/gkq340.PubMedPubMed CentralView ArticleGoogle Scholar
- Kim YT, Lee YR, Jin J, Han KH, Kim H, Kim JC, Lee T, Yun SH, Lee YW: Two different polyketide synthase genes are required for synthesis of zearalenone in Gibberella zeae. Mol Microbiol. 2005, 58: 1102-1113. 10.1111/j.1365-2958.2005.04884.x.PubMedView ArticleGoogle Scholar
- Dao HP, Mathieu F, Lebrihi A: Two primer pairs to detect OTA producers by PCR method. Int J Food Microbiol. 2005, 104: 61-67. 10.1016/j.ijfoodmicro.2005.02.004.PubMedView ArticleGoogle Scholar
- Artigot MP, Loiseau N, Laffitte J, Mas-Reguieg L, Tadrist S, Oswald IP, Puel O: Molecular cloning and functional characterization of two CYP619 cytochrome P450s involved in biosynthesis of patulin in Aspergillus clavatus. Microbiology. 2009, 155: 1738-1747. 10.1099/mic.0.024836-0.PubMedPubMed CentralView ArticleGoogle Scholar
- Proctor RH, Desjardins AE, Plattner RD, Hohn TM: A polyketide synthase gene required for biosynthesis of fumonisin mycotoxins in Gibberella fujikuroi mating population A. Fungal Genet Biol. 1999, 27: 100-112. 10.1006/fgbi.1999.1141.PubMedView ArticleGoogle Scholar
- Shimizu T, Kinoshita H, Ishihara S, Sakai K, Nagai S, Nihira T: Polyketide synthase gene responsible for citrinin biosynthesis in Monascus purpureus. Appl Environ Microbiol. 2005, 71: 3453-3457. 10.1128/AEM.71.7.3453-3457.2005.PubMedPubMed CentralView ArticleGoogle Scholar
- Ehrlich KC, Chang PK, Yu J, Cotty PJ: Aflatoxin biosynthesis cluster gene cypA is required for G aflatoxin formation. Appl Environ Microbiol. 2004, 70: 6518-6524. 10.1128/AEM.70.11.6518-6524.2004.PubMedPubMed CentralView ArticleGoogle Scholar
- Yu JH, Leonard TJ: Sterigmatocystin biosynthesis in Aspergillus nidulans requires a novel type I polyketide synthase. J Bacteriol. 1995, 177: 4792-4800.PubMedPubMed CentralGoogle Scholar
- Gardiner DM, Howlett BJ: Bioinformatic and expression analysis of the putative gliotoxin biosynthetic gene cluster of Aspergillus fumigatus. FEMS Microbiol Lett. 2005, 248: 241-248. 10.1016/j.femsle.2005.05.046.PubMedView ArticleGoogle Scholar
- Haese A, Schubert M, Herrmann M, Zocher R: Molecular characterization of the enniatin synthetase gene encoding a multifunctional enzyme catalyzing N-methyldepsipeptide formation in Fusarium scirpi. Mol Microbiol. 1993, 7: 905-914. 10.1111/j.1365-2958.1993.tb01181.x.PubMedView ArticleGoogle Scholar
- Scott-Craig JS, Panaccione DG, Pocard JA, Walton JD: The cyclic peptide synthetase catalyzing HC-toxin production in the filamentous fungus Cochliobolus carbonum is encoded by a 15.7-kilobase open reading frame. J Biol Chem. 1992, 267: 26044-26049.PubMedGoogle Scholar
- Weber G, Schörgendorfer K, Schneider-Scherzer E, Leitner E: The peptide synthetase catalyzing cyclosporine production in Tolypocladium niveum is encoded by a giant 45.8-kilobase open reading frame. Curr Genet. 1994, 26: 120-125. 10.1007/BF00313798.PubMedView ArticleGoogle Scholar
- Brenner S, Johnson M, Bridgham J, Golda G, Lloyd DH, Johnson D, Luo S, McCurdy S, Foy M, Ewan M, Roth R, George D, Eletr S, Albrecht G, Vermaas E, Williams SR, Moon K, Burcham T, Pallas M, DuBridge RB, Kirchner J, Fearon K, Mao J, Corcoran K: Gene expression analysis by massively parallel signature sequencing (MPSS) on microbead arrays. Nat Biotechnol. 2000, 18: 630-634. 10.1038/76469.PubMedView ArticleGoogle Scholar
- Bruno VM, Wang Z, Marjani SL, Euskirchen GM, Martin J, Sherlock G, Snyder M: Comprehensive annotation of the transcriptome of the human fungal pathogen Candida albicans using RNA-seq. Genome Res. 2010, 20: 1451-1458. 10.1101/gr.109553.110.PubMedPubMed CentralView ArticleGoogle Scholar
- Benjamini Y, Yekutieli D: The control of the false discovery rate in multiple testing under dependency. Ann Stat. 2001, 29: 1165-1188. 10.1214/aos/1013699998.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.