- Research article
- Open Access
Genome analysis of Daldinia eschscholtzii strains UM 1400 and UM 1020, wood-decaying fungi isolated from human hosts
© Chan et al. 2015
- Received: 28 August 2014
- Accepted: 10 November 2015
- Published: 18 November 2015
Daldinia eschscholtzii is a wood-inhabiting fungus that causes wood decay under certain conditions. It has a broad host range and produces a large repertoire of potentially bioactive compounds. However, there is no extensive genome analysis on this fungal species.
Two fungal isolates (UM 1400 and UM 1020) from human specimens were identified as Daldinia eschscholtzii by morphological features and ITS-based phylogenetic analysis. Both genomes were similar in size with 10,822 predicted genes in UM 1400 (35.8 Mb) and 11,120 predicted genes in UM 1020 (35.5 Mb). A total of 751 gene families were shared among both UM isolates, including gene families associated with fungus-host interactions. In the CAZyme comparative analysis, both genomes were found to contain arrays of CAZyme related to plant cell wall degradation. Genes encoding secreted peptidases were found in the genomes, which encode for the peptidases involved in the degradation of structural proteins in plant cell wall. In addition, arrays of secondary metabolite backbone genes were identified in both genomes, indicating of their potential to produce bioactive secondary metabolites. Both genomes also contained an abundance of gene encoding signaling components, with three proposed MAPK cascades involved in cell wall integrity, osmoregulation, and mating/filamentation. Besides genomic evidence for degrading capability, both isolates also harbored an array of genes encoding stress response proteins that are potentially significant for adaptation to living in the hostile environments.
Our genomic studies provide further information for the biological understanding of the D. eschscholtzii and suggest that these wood-decaying fungi are also equipped for adaptation to adverse environments in the human host.
- Genome sequencing
- Daldinia eschscholtzii
- Wood-inhabiting endophyte
- Wood-decaying fungi
Daldinia spp. belonging to the phylum Ascomycota and class Sordariomycetes, are known as endophytes or latent pathogens which inhabit woody host plants and remain in a dormant phase until the decay of wood or formation of perithecial stromata . The ascospores and conidia of these fungi are spread to neighboring trees via wind movements or fungivorous insects . Daldinia eschscholtzii has been isolated from dead trees , marine alga , insect  and recently also from human specimens [6, 7], displaying the great adaptation ability of this organism in a diverse host range. To date, there is no extensive analysis of the D. eschscholtzii genome, although this fungal species has been shown to produce potential bioactive compounds.
Secondary metabolites produced by D. eschscholtzii have potential medical and industrial applications. Zhang et al. [5, 8] reported the isolation of immunosuppressive compounds, including dalesconol A, B and C, daeschol A, 2, 16-dihydroxyl-benzo[j]fluoranthene and dalmanol A, from mantis-associated D. eschscholtzii. Helicascolide C, a new lactone with fungistatic activity against Cladosporium cucumerinum was isolated together with helicascolide A from an Indonesian marine algicolous-associated D. eschscholtzii . Daldinia spp. have been reported to produce volatile organic compounds (VOCs) [9–12] which can be developed in industrial applications for biofuel, biocontrol, and mycofumigation.
D. eschscholtzii has previously shown a typical feature of wood-decaying fungi, which is the production of enzymes for the degradation of lignocellulosic biomass, such as endoglucanase and β-glucosidase [3, 13]. As previously described, Daldinia spp. are known as type II soft-rot fungi that cause erosive degradation of lignocelluloses ; one sp. D. concentrica has been shown to be able to degrade the recalcitrant non-phenolic structures of lignin . This indicates that Daldinia spp. may have the ability to convert lignocellulosic biomass into different value-added products including biofuel, chemicals, and cheap carbon sources for fermentation, improved animal feeds, and human nutrients.
In this study, we present the genome of D. eschscholtzii UM 1400, an isolate obtained from human skin scraping, that enabled us to perform detailed analysis with the previously published genome of D. eschscholtzii UM 1020  for shared and common biological features. The genetic information of D. eschscholtzii UM 1400, combined with that of D. eschscholtzii UM 1020, will provide the knowledge for a deeper understanding of the biological nature of D. eschscholtzii.
Morphological and molecular identification
Genome sequencing and assembly
Summary of the assembly and annotation features of Daldinia eschscholtzii UM 1400 and UM 1020
Reads from small-insert library (Gb)
Reads from large-insert library (Gb)
Total Reads (Gb)
Assembly size (Mb)
Number of contigs (≥200 bp)
Contigs size (N50) (bp)
GC content of contigs (%)
Number of scaffolds (≥1000 bp)
Scaffolds size (N50) (bp)
GC content of scaffolds (%)
Number of predicted genes (≥99 bp)
Average gene length (bp)
Average number of exons per gene
CEGMA completeness score calculated from the complete gene set (%)
CEGMA completeness score calculated from both complete and partial gene sets (%)
Repetitive sequence (%)
The completeness of genome assembly was assessed using the CEGMA (Core Eukaryotic Genes Mapping Approach) software that evaluates the presence and completeness of a widely conserved set of 248 core eukaryotic genes [17, 18]. The standard CEGMA pipeline identified 235 out of the 248 core eukaryotic genes (94.76 %) in the UM 1400 genome assembly as complete, with an additional five core eukaryotic genes detected as partial (2.02 %). Similarly, out of the 248 core eukaryotic genes, a total of 239 (96.37 %) complete and two (0.81 %) partial core eukaryotic genes were detected when assessing the genome assembly of UM 1020 isolate. These results from CEGMA indicate that both assembled genomes cover most of the eukaryote’s gene space with many of genes complete and not fragmented onto multiple contigs.
The percentage of the repetitive sequences in both genomes (1.02 % in UM 1400 and 1.42 % in UM 1020) was lower than that reported for other Sordariomycetes genomes, for instance, 9.77 % in Magnaporthe grisea  and 10 % in Neurospora crassa . Of the repetitive sequences, transposable elements comprised 0.12 and 0.14 % in the UM 1400 and UM 1020 genomes, respectively. The transposable elements were classified into eight (UM 1400) and 12 (UM 1020) families with the subclass of Ty1_Copia most abundant in UM 1400 and the subclass of ISC1316 most abundant in UM 1020 (Additional file 1: Table S1). These data suggest that the D. eschscholtzii genomes are poor in repetitive sequences. However, it has been reported that repeat contents in Illumina-sequenced genomes are likely to be underestimated owing to a difficulty with the assembly of short repetitive reads into long repeat regions . Hence, the low repeat content of UM genomes is probably due to the Illumina technology that generates short reads that are prone to errors in the estimation of repetitive sequences, especially when the repetitive sequences are longer than the length of the sequencing reads [22–25].
The genome sequence of D. eschscholtzii UM 1400 has been deposited in the European Nucleotide Archive (ENA) under the accession numbers [ENA: CCED01000001-CCED01001944 and LK023387-LK023490]. The version described in this paper is the first version of this genome sequence.
The other shared gene families associated with fungus-host interactions included genes encoding the CFEM domain-containing protein (family SORD10851), FAS1 domain-containing protein (family SORD11403), and polysaccharide lyase family 4 protein (family SORD11288). The CFEM and FAS1 domains are present in fungal membrane proteins that could function as cell surface receptors or adhesion molecules in interactions with the host [27, 28]. The conservation of these genes suggests that both UM isolates may encode specific cell surface proteins with important roles in the interaction with its specific host. The polysaccharide lyase family 4 protein (rhamnogalacturonan lyase) cleaves the backbone of rhamnogalacturonan-I, which is a major component of plant cell wall polysaccharide pectin. A previous study showed that no endophytes tested have the ability to degrade pectin, and suggested that an endophyte is likely to be a latent pathogen if it can degrade pectic substances . Hence, the presence of gene encoding rhamnogalacturonan lyase indicates the ability of both isolates to produce pectic enzyme for pectin degradation. In line with this, it is implied that both UM isolates are likely to be latent pathogens, a lifestyle of D. eschscholtzii as described previously by another research group .
Both UM isolates shared gene families associated with stress response, for instance, genes encoding acid trehalase (family SORD14461) and ClpB protein (family SORD11060). The acid trehalases are involved in the assimilation of extracellular trehalose as a carbon source under nutrient limitation, as previously revealed in acid trehalase-deficient mutants of Saccharomyces cerevisiae and Aspergillus nidulans [30, 31]. The ClpB protein is an ATP-dependent molecular chaperone that plays an essential role in disaggregation and reactivation of the aggregated proteins in response to heat stress . These stress responses are necessary for fungal survival and adaptation in harsh environmental conditions.
Plant cell wall degrading enzymes
Comparison of total number of CAZymes with enzymatic activity for plant polysaccharide degradation
Endophyte/wood decaying fungi
Comparison of total number of CAZymes with ligninolytic auxiliary activities
Class II peroxidase
Vanillyl alcohol oxidase
Iron reductase domain
The presence of CAZymes with enzymatic activities for plant cell wall degradation implies that both human host-isolated D. eschscholtzii have once lived in the environment as wood-decaying fungi with degrading ability on plant biomass. These CAZymes are suggested to be required to degrade the wood cell consisting of the primary cell wall, secondary cell wall, and middle lamella, with each cell component containing different ratios of cellulose, hemicellulose, pectin and lignin. In addition, we identified six functional domains of cutinase (CAZyme family CE5) in the UM 1400 genome and five in the UM 1020 genome. Cutinases are critical for the initial fungal penetration through the cuticular barrier attached to the epidermal cell walls in aerial parts of plants, such as leaves, flowers, fruits and young stems . This indicates that both UM isolates have the potential ability to penetrate through not only lignified woody cell walls but also plant cuticle and epidermal cell walls as well.
The wood-inhabiting endophyte D. eschscholtzii has been reported to produce arrays of secondary metabolites that have potential applications in medical and biofuel industries, such as immunosuppressive polyketides and volatile organic compounds [5, 8–10]. In the UM 1400 and UM 1020 genomes, we identified 47 and 45 secondary metabolite backbone genes respectively, including those encoding lovastatin nonaketide synthase, conidial pigment polyketide synthase Alb1, dimethylallyl tryptophan synthase (DMATS) and citrinin polyketide synthase (Additional file 1: Table S6).
Lovastatin nonaketide synthase is involved in the biosynthesis of lovastatin, a cholesterol-lowering drug . The presence of this encoding gene suggests that both UM isolates may produce essential enzyme needed to manufacture the potent drug lovastatin for lowering blood cholesterol. Another polyketide synthase, Alb1, is responsible for the heptaketide naphtopyrene YWA1 synthesis in conidial pigmentation. Tsai et al.  reported that Aspergillus fumigatus produces Alb1 protein to synthesize the conidial pigment via the pentaketide pathway. This indicates that heptaketide synthase may be involved in the initiation of pentaketide melanin biosynthesis in D. eschscholtzii. Melanin protects fungal spores and mycelium against environmental stresses, including desiccation, oxidizing agents and ultraviolet (UV) light. Thus, melanin production may be a protective trait that allows D. eschscholtzii to survive in harsh conditions such as drought that triggers desiccation and osmotic stress.
Dimethylallyl tryptophan synthase (DMATS) and citrinin polyketide synthase are involved in the synthesis of ergot alkaloids and antibiotic citrinin, respectively. As previously reported, ergot alkaloids were shown to be poisonous to herbivores , while citrinin had antimicrobial activity against pathogens . These bioactive compounds may play a similar role in both UM isolates to confer beneficial protection to its host plant from the attacks of herbivores and pathogens subjected to further confirmation. A previous study reported that an endophytic Daldinia eschscholtzii EC12 produces volatile organic compounds that are active against a broad range of plant pathogens [9, 10].
Secreted peptidases facilitate fungal penetration and colonization of the host plant by degrading plant cell wall structural proteins and plant defense-related proteins [40, 41]. Examples are subtilisin-like peptidases (MEROPS subfamily S08A) and metallopeptidases (MEROPS families M35 and M36).
Subtilisin-like peptidases are serine peptidases that have been found to be associated with colonization of the host by endophytes  and plant pathogenic fungi . We identified 15 genes (eight in the UM 1400 genome and seven in the UM 1020 genome) encoding subtilisin-like peptidases (Additional file 1: Table S7). The metallopeptidases are known to cleave the glycoproteins of extracellular matrix that have been implicated in host resistance mechanisms against pathogen invasion [43, 44]. For instance, the fungalysin of M36 family was shown to truncate non-structural host resistance proteins . The presence of genes encoding penicillolysin of M35 family and fungalysin of M36 family in both genomes indicates the ability of UM isolates to inactivate proteinaceous components from the plant defense response.
A protein blast analysis against the pathogen-host interaction database (PHI-base) revealed 602 and 606 putative PHI genes (>50 % identity; >70 % subject coverage) in the genomes of UM 1400 and UM 1020 respectively. With the Eukaryotic Orthologous Group (KOG) functional classification, the putative PHI genes were distributed into 22 functional categories with a higher number assigned to the category of signal transduction mechanisms (Additional file 2: Figure S1).
The UM genomes contained arrays of putative genes encoding signaling components, and here, we discuss those involved in mitogen-activated protein kinase (MAPK) signaling pathways. The MAPK signaling pathways are commonly found in all eukaryotes and are known to be involved in cell growth, differentiation, and stress response. From the genomic analysis, three putative MAPK signaling pathways were proposed to be present in D. eschscholtzii UM 1400 and UM 1020 isolates, including the cell wall integrity pathway mediated by the Ssk2/Ssk22-Pbs2-Osm1 cascade, the osmoregulation pathway mediated by the Bck1-Mkk1/Mkk2-Mps1 cascade, and the mating/filamentation pathway mediated by the Mst11-Ste7-Gpmk1 cascade (Additional file 1: Table S8).
Numerous homologs of components of the osmoregulation pathway were identified in the UM genomes, including Sln1, Hik1, Sho1, Cdc42, Mst20 (Ste20 homolog), Mst50 (Ste50 homolog), Mst11 (Ste11 homolog), Pbs2, Osm1 (Hog1 homolog), and Ssk2/Ssk22. The osm1 gene was shown to encode a functional homolog of MAPK Hog1 and to be required in response to osmotic stress . As referring to the high-osmolarity glycerol (HOG) pathway in Saccharomyces cerevisiae , two upstream branches were predicted to activate the MAPK Pbs2-MAPK Osm1 module in both UM isolates. One branch consisted of MAPKK Ssk2/Ssk22 and a two-component histidine kinase phospho-relay system Sln1-Ypd1-Ssk1, with lacks homologs Ypd1 and Ssk1; another branch consisted of a putative membrane protein Sho1, Cdc42, Mst11, Mst20, and Mst50. Besides osmotic stress, the osmoregulation pathway is also required for adaptation to oxidative stress, thermal stress, cellular morphogenesis regulation, and cell wall functionality .
The UM genomes contained several homologs of components of the cell wall integrity pathway, including the GTP-binding protein Rho1, MAPKKK Bck1, MAPKK Mkk1/Mkk2, and MAPK Mps1 (Slt2 homolog). Generally, this MAPK pathway is essential for cell wall integrity and pathogenesis . Some MAPK Slt2 homologs are also involved in other roles, such as conidium germination and polarized growth in Aspergillus nidulans , and response to various stresses, including oxidative and osmotic stresses, in Candida albicans .
The homolog of MAPK Gmpk1 of the mating/filamentation pathway was identified in both UM genomes. In Fusarium graminearum, Gmpk1 is required to regulate mating, conidial production, and pathogenicity as well as the early induction of extracellular endoglucanase, xylanolytic, proteolytic and lipolytic activities [50, 51]. Other identified homologous components of this pathway included a MAPKKK Mst11 (Ste11 homolog), a MAPKK Ste7, an adaptor protein Ste50, a PAK kinase Mst20 (Ste20 homolog), two small GTP-binding proteins Ras2 and Cdc42, a Gα subunit Gpa2. These components were previously reported to be involved in the activation of the mating/filamentation pathway in well-characterized Saccharomyces cerevisiae . However, the homolog of G protein-coupled, seven-transmembrane receptor Gpr1 was not found in the UM genomes. One Ste12-like transcription factor, Cst1 homolog was identified in the UM genomes. As previously reported, this downstream transcription factor regulates genes involved in penetration and infectious growth in Colletotrichum lagenarium .
Although the MAPKKK-MAPKK-MAPK cascades are generally conserved in eukaryotes, the UM isolates seemed to lack significant homologs of upstream protein kinases, such as Ypd1 and Ssk1 in the osmoregulation pathway, and Gpr1 in the mating/filamentation pathway. This suggests that the upstream components in our UM isolates may be different from those in other well-characterized organisms, like S. cerevisiae , and may be novel receptor kinases for sensing environmental signals.
Adaptation-associated stress response proteins
Daldinia eschscholtzii has been isolated from diverse environments [3–7] where it may be subjected to many extreme conditions. The transition from a moderate environment to a hostile environment causes drastic changes in various parameters, including osmotic changes, pH changes, thermal changes, nutrient deprivations, as well as oxidative and nitrosative stresses. In both UM genomes, we identified numerous stress-responsive genes as listed in Table S9 in Additional file 1.
Daldinia spp. appear to be adapted to survive during periods of drought in the natural environment, and even when their woody host plant has been fire-damaged . These harsh conditions result in osmotic and thermal stresses to Daldinia spp. To maintain cellular turgor and prevent water loss, high concentrations of osmolytes, like glycerol, erythritol, mannitol, or trehalose are generated. The gene encoding osmotic stress-responsive proteins were identified in the UM genomes, including os-4 orthologue, os-1/nik-1 orthologue, and pbsA orthologue that are involved in osmolytes accumulation, tpsA and orlA that are involved in osmolyte trehalose biosynthesis, and gfdB that is involved in osmolyte glycerol biosynthesis. The increased production of osmolytes could induce the formation of vegetative structures conferring resistance to drought condition , likes the stromatic structures formed by Daldinia spp. for survival in drought . Numerous genes encoding thermal stress-responsive proteins were found in the genomes, for instance, genes encoding heat shock proteins (hsp70, hsp78, hsp104, hsf1) which are induced to refold or degrade damaged proteins, to unfold aggregated proteins, and also to help in stabilizing proteins and membranes . Unceasing wood decay will change the chemical composition and physical structure of wood which will, in turn, lead to nutrient deprivation stress. To tolerate this stressful condition, the alternate nutrient sources may be assimilated by expressing the genes associated with sources metabolism and nutrients uptake, such as treB, mep1, mep2, prnB, and prnC. The treB gene encodes a neutral trehalase that partially contributes to the energy requirements of spore germination under carbon limitation, as shown in the tre mutant of Aspergillus nidulans . The mep1 and mep2 genes are the examples of genes encoding proteins involved in nitrogen assimilation, and are predominantly expressed at low concentrations of ammonium or on poor nitrogen sources. The ammonium permease encoded by the mep2 gene has been shown to control nitrogen starvation-induced filamentous growth in Candida albicans via interaction with Ras1 . In addition, fungi are able to utilize amino acids as sole nitrogen and/or carbon sources. In response to amino acid starvation, the transcription of the genes involved in amino acid biosynthesis are activated, such as prnB and prnC genes , both of which were found in the UM genomes.
Reactive oxygen species (ROS) and reactive nitrogen species (RNS) produced by hosts are harmful to fungi by causing damage to their proteins, lipid membranes, and deoxyribonucleic acid . In order to survive in this harsh environment, fungi must have mechanisms to detoxify these reactive molecular species and repair the cellular damages triggered by the oxidative and nitrosative stresses. The UM isolates were found to contain genes encoding antioxidant enzymes (sod1, sod2, sodA, cat1, tsa1, tsa3, grx5, gpxA, msrA, msrB) and enzymes involved in the production of secondary metabolites with antioxidant function (tpsA, orlA) to handle ROS, as well as the gene encoding nitrosative stress-responsive proteins (fhbA) to cope with RNS. The plant pathogenic fungus Botrytis cinerea has been reported to produce these proteins to thrive against the oxidative and nitrosative environments generated by host plant cells [60, 61]. These enzymes have also been implicated in the defense of opportunistic fungal pathogens (Candida albicans, Cryptococcus neoformans and Aspergillus fumigatus) against the ROS and RNS produced by human phagocytes . In the case of wood-decaying fungi, the extracellular hydrogen peroxide provided by oxidative enzymes is involved in the generation of highly reactive oxidants or hydroxyl radicals via the Fenton reaction with the presence of iron cofactor . These radicals are involved in the degradation process. The high level of generated ROS is coped by the intracellular antioxidant enzymes to prevent fungal cell damage, as investigated in the previous study on the endogenous oxidative stress response of Coriolus versicolor . Our CAZymes analysis identified oxidative enzymes in both UM genomes, with GMC oxidoreductases (family AA3) present in a high number (Table 3). This GMC oxidase (family AA3) has been thought to play an important role in peroxide production in the wood-decaying fungus Gloeophyllum trabeum . Other enzymes, the copper radical oxidase (family AA5), FAD-linked oxidoreductase (family AA3) and glucose oxidase-like protein (family AA3) have been demonstrated to be potentially involved in extracellular peroxide production in Postia placenta .
The Fenton reaction requires iron as the cofactor of peroxidase enzymes for degradation activity. However, iron is sequestered by high-affinity iron binding proteins; thus, the iron acquisition system is required for wood-decaying fungi under iron starvation condition . The UM genomes featured genes involved in iron acquisition, namely genes encoding iron permease (Ftr1 orthologue) and mitochondrial ornithine carrier (AmcA orthologue). The Ftr1 protein is required for high-affinity iron uptake in the reductive iron uptake system , while the AmcA protein is involved in the supply of ornithine for siderophore biosynthesis . The Ftr1 protein was shown to be up-regulated during the growth of P. placenta on cellulose medium .
Changes in pH can be encountered upon environmental transition. The gene encoding Pac1 ortholog of Fusarium graminearum was identified in the UM genomes. This gene has been previously reported to encode a pH regulator factor regulating the production of secondary metabolite in F. graminearum . Overall, the UM isolates harbored many genes encoding stress response proteins that cope with triggered stresses under adverse conditions in their natural habitats. This feature could also serve to their advantage in surviving the adverse microenvironments of human niches.
The genomic analysis of both UM isolates revealed a common set of putative domains or genes that improves our understanding of the biological nature of D. eschscholtzii. The environmental origin of these isolates is suggested by the identification of putative CAZyme arrays and genes encoding secreted peptidases related to plant cell wall degradation. As D. eschscholtzii has hitherto never been associated with human infections, our UM isolates might have been entering into humans via the exposure of open wounds to the decaying wood material containing this organism and have been surviving in the human without causing any disease. Both UM genomes displayed a wide range of adaptation-associated stress response genes that are required by fungi for adaptation to hostile conditions in their natural habitat. These genes most likely also confer a selective advantage for survival and adaptation in adverse microenvironments in the human host. Our genomic analysis also revealed other biological features, such as the identified genes encoding MAPK signaling pathway components that suggest three MAPK signaling cascades, and the identified secondary metabolite backbone genes that indicate the potential of the UM isolates to produce various bioactive secondary metabolites. The biological functions of predicted genes have to be validated by further studies using appropriate approaches such as insertional mutagenesis, serial analysis of gene expression, microarray analysis, proteomics, and metabolomics.
As no patient information is disclosed, it was considered unnecessary to apply for ethical approval from the University Malaya Medical Centre (UMMC) Medical Ethics Committee (http://umresearch.um.edu.my/doc/File/UMREC/6_CODE%20OF%20RESEARCH%20ETHICS%20%20IN%20UNIVERSITY%20OF%20MALAYA.pdf).
Both UM 1400 and UM 1020 isolates were recovered from a collection of fungi routinely cultured and archived in the Mycology diagnostic laboratory, UMMC. These isolates were grown on SDA plates at 30 °C and maintained on SDA slants at 4 °C until required for research use.
Fungal cultures on SDA were observed for cultural characteristics. Slide cultures mounted with lactophenol cotton blue stain were observed under the light microscope for anamorphic structures. For SEM examination, the cultures were fixed, dried, mounted on a specimen stub using electrically conductive double-sided adhesive tape, and sputter-coated with gold before observing under the XL-30 ESEM microscope (Philips, Netherlands) for the surface topography of conidia and conidiophores.
Fungal DNA was extracted using ZR Fungal/Bacterial DNA MiniPrepTM (Zymo Research, USA), according to the manufacturer’s protocol. The specific primer pair ITS1 and ITS4 was used to amplify the region of ITS1-5.8S-ITS2 rDNA, as previously described . PCR products were visualized by gel electrophoresis and purified prior to Sanger sequencing. ITS sequences were searched against the NCBI nucleotide database to determine fungal identity. For phylogenetic tree analysis, the complete ITS1-5.8S-ITS2 sequence was collected from each Daldinia species available in the GenBank database. Multiple sequence alignments of all data-mined ITS sequences were generated using M-Coffee  which uses other packages to compute the alignments and uses T-Coffee to combine all of these alignments into one unique final alignment. Phylogenetic analysis was then performed using MrBayes version 3.2.1 . Bayesian Markov Chain Monte Carlo (MCMC) analysis was performed by sampling across the entire general time reversible (GTR) model space. A total of 250,000 generations were run with a sampling frequency of 100, and diagnostics were calculated for every 1000 generations. A burn-in setting of 25 % was used to discard the first 625 trees. Convergence was assessed as suggested in the manual , with a standard deviation of split frequencies below 0.01, no obvious trend for the plot of the generation versus the log probability of the data, and the potential scale reduction factor (PSRF) close to 1.0 for all parameters.
Genome sequencing and assembly
Genomic DNA was extracted as previously described . The genomes of UM 1020 and UM 1400 were sequenced separately at different time periods. The sequencing and analysis of UM 1020 were reported previously as a genome announcement . UM 1400 was sequenced with two libraries of 500 bp and 5 kb insert size, using Illumina HiSeq 2000 sequencer. The genome was assembled using Velvet version 1.2.07  and the primary scaffolds from the Velvet assembly were further scaffolded with SSPACE Basic version 2.0  and gap filled using GapFiller version 1.10  utilizing paired-end information from both libraries. The genome completeness of both assemblies was accessed using CEGMA version 2.4 . Whole genome alignments and comparison between the assembled scaffolds of both UM isolates were performed using NUCmer version 3.1 from MUMmer package version 3.23  with default parameters. Genome assemblies of both UM isolates were then aligned and ordered using MAUVE version 2.3.1 . Dot-plot was generated using mummerplot with the –color parameter and ordered set of UM 1020 genome sequence.
Gene prediction and annotation
Gene prediction and functional annotation of both genome assemblies were carried out using the same methods to allow direct comparison. Prior to gene prediction, the repetitive elements, including interspersed repeats and low complexity DNA sequences were masked throughout the genomes using RepeatMasker version open-3.3.0 with the Repbase fungal library version rm-20,120,418. The Ab initio gene prediction was then performed on repeats-masked genomes using GeneMark-ES version 2.3 . Predicted proteins were functionally annotated by local BLAST similarity searches against NCBI NR and SwissProt databases. Results from blast searches were analyzed for Gene Ontology (GO) and KEGG pathway using local Blast2GO tools . Functional classification of the predicted proteins was performed using KOG . Interpro analysis with Pfam database was used to annotate the predicted proteins based on protein domain families using InterProScan 5 . The rRNA and tRNA of the genomes were predicted using RNAmmer version 1.2  and tRNAscan-SE version 1.3.1  respectively.
Predicted protein models were subjected to dbcan which is a web server and database for automated carbohydrate-active enzyme annotation . Peptidases were identified by querying against MEROPS database release 9.9 using the batch blast service . Secretome analysis was performed using the method as previously described , in which the predictions of cleavage sites and the signal peptide/non-signal peptide were carried out using SignalP version 4.1 . Only secreted proteins without transmembrane (TM) domains and those with single TM present at the N-terminal 40 amino acids corresponding to secretion signals were selected. The presence of TM domains was identified using TMHMM version 2.0 . Secondary metabolite backbone genes in both genomes were predicted using web-based software SMURF . Pathogenicity-associated genes and stress response genes were identified by BLASTP search against local databases, which were downloaded and built from the pathogen-host interaction database (PHI-base)  and fungal stress response database (FSRD)  respectively, with criteria of BLASTP e-value threshold of less than or equal to 1e-5, percentage identity of more than 50 % and subject coverage of more than 70 %. Putative transposon elements were predicted using Transposon-PSI  which performed PSI-TBLASTN search of our genomes with a collection of (retro-)transposon ORF homology profiles to identify statistically significant alignments encoded by diverse families of transposable elements.
To identify orthologs of D. eschscholtzii within the class Sordariomycetes, protein sequences of all currently available Sordariomycetes genomes were downloaded from the databases. Proteins sequences (≥33 amino acids) from both UM isolates and the data-mined reference genomes were clustered using OrthoMCL software version 2.02  which performs all-against-all BLASTP searches of all proteins, with reciprocal best blast hits from distinct genomes recognized as orthologs. OrthoMCL applied Markov Cluster algorithm  and was run with BLAST e-value cut-off of 1e-5 using an inflation parameter of 1.5.
For the phylogenomic tree construction, proteomes from Sordariomycetes fungi (used in gene family analysis) including our UM isolates were clustered using OrthoMCL software version 2.02 for within-class comparisons, with two Dothideomycetes fungi as outgroups. Individual multiple sequence alignments of 332 single-copy orthologs from the gene family analysis was generated using ClustalW version 2.0 . The spurious sequences and poorly aligned regions were removed by using trimAL, and the multiple alignments were concatenated into a superalignment with 174,267 characters. Subsequently, ProtTest version 3.2 was run with AIC calculation on the alignment to select the best-fit substitution model . A phylogenomic tree was constructed by using both RAxML version 7.7.9  and MrBayes version 3.2.1 . RAxML was run with model PROTGAMMAJTTF to search for the best-scoring Maximum Likelihood tree followed by 100 bootstrap replicates. Convergence was observed after 50 replicates using -I autoMRE option in RAxML. MrBayes was run using Jones amino acid model with gamma-distributed rate variation across sites and a proportion of invariable sites. The MCMC was run with a sampling frequency of 100 for 250,000 generations, and burn-in setting of 25 %.
This study was supported by High Impact Research MoE Grant UM.C/625/1/HIR/MOHE/MED/31 (Account no. H-20001-00-E000070) from the Ministry of Education Malaysia and Postgraduate Research Grant (PPP) PV051/2012A from the University of Malaya. The genomic sequence data of Daldinia eschscholtzii EC12, Hypoxylon sp. CI-4A, Hypoxylon sp. CO27-5, and Hypoxylon sp. EC38 were produced by the US Department of Energy Joint Genome Institute http://www.jgi.doe.gov/ in collaboration with the user community, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Srutka P, Pazoutova S, Kolarik M. Daldinia decipiens and Entonaema cinnabarina as fungal symbionts of Xiphydria wood wasps. Mycol Res. 2007;111(Pt 2):224–31.View ArticlePubMedGoogle Scholar
- Johannesson H, Gustafsson M, Stenlid J. Local population structure of the wood decay Ascomycete Daldinia Loculata. Mycologia. 2001;93:440–6.View ArticleGoogle Scholar
- Karnchanatat A, Petsom A, Sangvanich P, Piapukiew J, Whalley AJS, Reynolds CD, et al. A novel thermostable endoglucanase from the wood-decaying fungus Daldinia eschscholzii (Ehrenb.:Fr.) Rehm. Enzyme Microb Technol. 2008;42:404–13.View ArticleGoogle Scholar
- Tarman K, Palm GJ, Porzel A, Merzweiler K, Arnold N, Wessjohann LA, et al. Helicascolide C, a new lactone from an Indonesian marine algicolous strain of Daldinia eschscholzii (Xylariaceae, Ascomycota). Phytochem Lett. 2012;5:83–6.View ArticleGoogle Scholar
- Zhang YL, Ge HM, Zhao W, Dong H, Xu Q, Li SH, et al. Unprecedented immunosuppressive polyketides from Daldinia eschscholzii, a mantis-associated fungus. Angew Chemie Int Ed. 2008;47:5823–6.View ArticleGoogle Scholar
- Ng KP, Ngeow YF, Yew SM, Hassan H, Soo-Hoo TS, Na SL, et al. Draft genome sequence of Daldinia eschscholzii isolated from blood culture. Eukaryot Cell. 2012;11:703–4.PubMed CentralView ArticlePubMedGoogle Scholar
- Yew SM, Chan CL, Lee KW, Na SL, Tan R, Hoh C-C, et al. A five-year survey of dematiaceous fungi in a tropical hospital reveals potential opportunistic species. PLoS One. 2014;9:e104352.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhang YL, Zhang J, Jiang N, Lu YH, Wang L, Xu SH, et al. Immunosuppressive polyketides from mantis-associated Daldinia eschscholzii. J Am Chem Soc. 2011;133:5931–40.View ArticlePubMedGoogle Scholar
- Mends MT, Yu E, Strobel GA, Riyaz-Ul-Hassan S, Booth E, Geary B, et al. An endophytic Nodulisporium sp. producing volatile organic compounds having bioactivity and fuel potential. J Pet Environ Biotechnol. 2012;03:117.Google Scholar
- Strobel GA. Methods of discovery and techniques to study endophytic fungi producing fuel-related hydrocarbons. Nat Prod Rep. 2014;31:259–72.View ArticlePubMedGoogle Scholar
- Pažoutová S, Follert S, Bitzer J, Keck M, Surup F, Šrůtka P, et al. A new endophytic insect-associated Daldinia species, recognised from a comparison of secondary metabolite profiles and molecular phylogeny. Fungal Divers. 2013;60:107–23.View ArticleGoogle Scholar
- Yu ET, Tran-Gyamfi M, Strobel GA, Taatjes C, Hadi MZ. VOC profile of endophytic fungi is altered by nature of lignocellulosic biomass feedstock. NASAreport, 2013, USA- National Aeronautics and Space Administration, Washington, D. C.Google Scholar
- Karnchanatat A, Petsom A, Sangvanich P, Piaphukiew J, Whalley AJS, Reynolds CD, et al. Purification and biochemical characterization of an extracellular beta-glucosidase from the wood-decaying fungus Daldinia eschscholzii (Ehrenb.:Fr.) Rehm. FEMS Microbiol Lett. 2007;270:162–70.View ArticlePubMedGoogle Scholar
- Kuhad RC, Kuhar S, Sharma KK, Shrivastava B. Microorganisms and enzymes involved in lignin degradation vis-à-vis production of nutritionally rich animal feed: an overview. In: Kuhad RC, Singh A, editors. Biotechnol environ manag resour recover. India: Springer India; 2013. p. 44.View ArticleGoogle Scholar
- Shary S, Ralph SA, Hammel KE. New insights into the ligninolytic capability of a wood decay ascomycete. Appl Environ Microbiol. 2007;73:6691–4.PubMed CentralView ArticlePubMedGoogle Scholar
- Ju YM, Rogers JD, San Martin F. A revision of the genus Daldinia. Mycotaxon. 1997;61:243–93.Google Scholar
- Parra G, Bradnam K, Korf I. CEGMA: a pipeline to accurately annotate core genes in eukaryotic genomes. Bioinformatics. 2007;23:1061–7.View ArticlePubMedGoogle Scholar
- Parra G, Bradnam K, Ning Z, Keane T, Korf I. Assessing the gene space in draft genomes. Nucleic Acids Res. 2009;37:289–97.PubMed CentralView ArticlePubMedGoogle Scholar
- Dean RA, Talbot NJ, Ebbole DJ, Farman ML, Mitchell TK, Orbach MJ, et al. The genome sequence of the rice blast fungus Magnaporthe grisea. Nature. 2005;434:980–6.View ArticlePubMedGoogle Scholar
- Galagan JE, Calvo SE, Borkovich KA, Selker EU, Read ND, Jaffe D, et al. The genome sequence of the filamentous fungus Neurospora crassa. Nature. 2003;422:859–68.View ArticlePubMedGoogle Scholar
- Gnerre S, Maccallum I, Przybylski D, Ribeiro FJ, Burton JN, Walker BJ, et al. High-quality draft assemblies of mammalian genomes from massively parallel sequence data. Proc Natl Acad Sci U S A. 2011;108:1513–8.PubMed CentralView ArticlePubMedGoogle Scholar
- Ohm RA, Feau N, Henrissat B, Schoch CL, Horwitz BA, Barry KW, et al. Diverse lifestyles and strategies of plant pathogenesis encoded in the genomes of eighteen Dothideomycetes fungi. PLoS Pathog. 2012;8:e1003037.PubMed CentralView ArticlePubMedGoogle Scholar
- Li R, Zhu H, Ruan J, Qian W, Fang X, Shi Z, et al. De novo assembly of human genomes with massively parallel short read sequencing. Genome Res. 2010;20:265–72.PubMed CentralView ArticlePubMedGoogle Scholar
- Schatz MC, Delcher AL, Salzberg SL. Assembly of large genomes using second-generation sequencing. Genome Res. 2010;20:1165–73.PubMed CentralView ArticlePubMedGoogle Scholar
- Tsai IJ, Otto TD, Berriman M. Improving draft assemblies by iterative mapping and assembly of short reads to eliminate gaps. Genome Biol. 2010;11:R41.PubMed CentralView ArticlePubMedGoogle Scholar
- Kurtz S, Phillippy A, Delcher AL, Smoot M, Shumway M, Antonescu C, et al. Versatile and open software for comparing large genomes. Genome Biol. 2004;5:R12.PubMed CentralView ArticlePubMedGoogle Scholar
- Kulkarni RD, Kelkar HS, Dean RA. An eight-cysteine-containing CFEM domain unique to a group of fungal membrane proteins. Trends Biochem Sci. 2003;28:118–21.View ArticlePubMedGoogle Scholar
- Liu T, Chen G, Min H, Lin F. MoFLP1, encoding a novel fungal fasciclin-like protein, is involved in conidiation and pathogenicity in Magnaporthe oryzae. J Zhejiang Univ Sci B. 2009;10:434–44.PubMed CentralView ArticlePubMedGoogle Scholar
- Choi YW, Hodgkiss IJ, Hyde KD. Enzyme production by endophytes of Brucea javanica. J Agric Technol. 2005;1:55–66.Google Scholar
- Nwaka S, Mechler B, Holzer H. Deletion of the ATH1 gene in Saccharomyces cerevisiae prevents growth on trehalose. FEBS Lett. 1996;386:235–8.View ArticlePubMedGoogle Scholar
- d’Enfert C, Fontaine T. Molecular characterization of the Aspergillus nidulans treA gene encoding an acid trehalase required for growth on trehalose. Mol Microbiol. 1997;24:203–16.View ArticlePubMedGoogle Scholar
- Zolkiewski M, Zhang T, Nagy M. Aggregate reactivation mediated by the Hsp100 chaperones. Arch Biochem Biophys. 2012;520:1–6.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhao Z, Liu H, Wang C, Xu J-R. Correction: comparative analysis of fungal genomes reveals different plant cell wall degrading capacity in fungi. BMC Genomics. 2014;15:6.PubMed CentralView ArticlePubMedGoogle Scholar
- Levasseur A, Drula E, Lombard V, Coutinho PM, Henrissat B. Expansion of the enzymatic repertoire of the CAZy database to integrate auxiliary redox enzymes. Biotechnol Biofuels. 2013;6:41.PubMed CentralView ArticlePubMedGoogle Scholar
- Kolattukudy PE. Enzymatic penetration of the plant cuticle by fungal pathogens. Annu Rev Phytopathol. 1985;23:223–50.View ArticleGoogle Scholar
- Campbell CD, Vederas JC. Biosynthesis of lovastatin and related metabolites formed by fungal iterative PKS enzymes. Biopolymers. 2010;93:755–63.View ArticlePubMedGoogle Scholar
- Tsai HF, Fujii I, Watanabe A, Wheeler MH, Chang YC, Yasuoka Y, et al. Pentaketide melanin biosynthesis in Aspergillus fumigatus requires chain-length shortening of a heptaketide precursor. J Biol Chem. 2001;276:29292–8.View ArticlePubMedGoogle Scholar
- Bacetty AA, Snook ME, Glenn AE, Noe JP, Hill N, Culbreath A, et al. Toxicity of endophyte-infected tall fescue alkaloids and grass metabolites on Pratylenchus scribneri. Phytopathology. 2009;99:1336–45.View ArticlePubMedGoogle 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–7.PubMed CentralView ArticlePubMedGoogle Scholar
- Sreedhar L, Kobayashi DY, Bunting TE, Hillman BI, Belanger FC. Fungal proteinase expression in the interaction of the plant pathogen Magnaporthe poae with its host. Gene. 1999;235:121–9.View ArticlePubMedGoogle Scholar
- Naumann TA, Wicklow DT, Price NPJ. Identification of a chitinase-modifying protein from Fusarium verticillioides: truncation of a host resistance protein by a fungalysin metalloprotease. J Biol Chem. 2011;286:35358–66.PubMed CentralView ArticlePubMedGoogle Scholar
- Reddy PV, Lam CK, Belanger FC. Mutualistic fungal endophytes express a proteinase that is homologous to proteases suspected to be important in fungal pathogenicity. Plant Physiol. 1996;111:1209–18.PubMed CentralView ArticlePubMedGoogle Scholar
- Bowles DJ. Defense-related proteins in higher plants. Annu Rev Biochem. 1990;59:873–907.View ArticlePubMedGoogle Scholar
- Dow JM, Davies HA, Daniels MJ. A metalloprotease from Xanthomonas campestris that specifically degrades proline/hydroxyproline-rich glycoproteins of the plant extracellular matrix. Mol plant-microbe Interact. 1998;11:1085–93.View ArticlePubMedGoogle Scholar
- Dixon KP, Xu JR, Smirnoff N, Talbot NJ. Independent signaling pathways regulate cellular turgor during hyperosmotic stress and appressorium-mediated plant infection by Magnaporthe grisea. Plant Cell. 1999;11:2045–58.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhao X, Mehrabi R, Xu J-R. Mitogen-activated protein kinase pathways and fungal pathogenesis. Eukaryot Cell. 2007;6:1701–14.PubMed CentralView ArticlePubMedGoogle Scholar
- Brown AJP, Budge S, Kaloriti D, Tillmann A, Jacobsen MD, Yin Z, et al. Stress adaptation in a pathogenic fungus. J Exp Biol. 2014;217:144–55.PubMed CentralView ArticlePubMedGoogle Scholar
- Bussink H. A mitogen-activated protein kinase (MPKA) is involved in polarized growth in the filamentous fungus, Aspergillus nidulans. FEMS Microbiol Lett. 1999;173:117–25.View ArticlePubMedGoogle Scholar
- Navarro-Garcia F. The MAP kinase Mkc1p is activated under different stress conditions in Candida albicans. Microbiology. 2005;151:2737–49.View ArticlePubMedGoogle Scholar
- Jenczmionka NJ, Maier FJ, Lösch AP, Schäfer W. Mating, conidiation and pathogenicity of Fusarium graminearum, the main causal agent of the head-blight disease of wheat, are regulated by the MAP kinase gpmk1. Curr Genet. 2003;43:87–95.PubMedGoogle Scholar
- Jenczmionka NJ, Schäfer W. The Gpmk1 MAP kinase of Fusarium graminearum regulates the induction of specific secreted enzymes. Curr Genet. 2005;47:29–36.View ArticlePubMedGoogle Scholar
- Gustin MC, Albertyn J, Alexander M, Davenport K. MAP kinase pathways in the yeast Saccharomyces cerevisiae. Microbiol Mol Biol Rev. 1998;62:1264–300.PubMed CentralPubMedGoogle Scholar
- Tsuji G, Fujii S, Tsuge S, Shiraishi T, Kubo Y. The Colletotrichum lagenarium Ste12-like gene CST1 is essential for appressorium penetration. Mol Plant-Microbe Interact. 2003;16:315–25.View ArticlePubMedGoogle Scholar
- Stadler M, Læssøe T, Fournier J, Decock C, Schmieschek B, Tichy H-V, et al. A polyphasic taxonomy of Daldinia (Xylariaceae). Stud Mycol. 2014;77:1–143.PubMed CentralView ArticlePubMedGoogle Scholar
- Son H, Lee J, Lee YW. Mannitol induces the conversion of conidia to chlamydospore-like structures that confer enhanced tolerance to heat, drought, and UV in Gibberella zeae. Microbiol Res. 2012;167:608–15.View ArticlePubMedGoogle Scholar
- D’Enfert C, Bonini BM, Zapella PD, Fontaine T, da Silva AM, Terenzi HF. Neutral trehalases catalyse intracellular trehalose breakdown in the filamentous fungi Aspergillus nidulans and Neurospora crassa. Mol Microbiol. 1999;32:471–83.View ArticlePubMedGoogle Scholar
- Biswas K, Morschhäuser J. The Mep2p ammonium permease controls nitrogen starvation-induced filamentous growth in Candida albicans. Mol Microbiol. 2005;56:649–69.View ArticlePubMedGoogle Scholar
- Tazebay UH, Sophianopoulou V, Scazzocchio C, Diallinas G. The gene encoding the major proline transporter of Aspergillus nidulans is upregulated during conidiospore germination and in response to proline induction and amino acid starvation. Mol Microbiol. 1997;24:105–17.View ArticlePubMedGoogle Scholar
- Brown AJ, Haynes K, Quinn J. Nitrosative and oxidative stress responses in fungal pathogenicity. Curr Opin Microbiol. 2009;12:384–91.PubMed CentralView ArticlePubMedGoogle Scholar
- Mayer AM, Staples RC, Gil-ad NL. Mechanisms of survival of necrotrophic fungal plant pathogens in hosts expressing the hypersensitive response. Phytochemistry. 2001;58(1):33–41.View ArticlePubMedGoogle Scholar
- Turrion-Gomez JL, Eslava AP, Benito EP. The flavohemoglobin BCFHG1 is the main NO detoxification system and confers protection against nitrosative conditions but is not a virulence factor in the fungal necrotroph Botrytis cinerea. Fungal Genet Biol. 2010;47:484–96.View ArticlePubMedGoogle Scholar
- Aguirre J, Hansberg W, Navarro R. Fungal responses to reactive oxygen species. Med Mycol. 2006;44 Suppl 1:101–7.View ArticleGoogle Scholar
- Martinez D, Challacombe J, Morgenstern I, Hibbett D, Schmoll M, Kubicek CP, et al. Genome, transcriptome, and secretome analysis of wood decay fungus Postia placenta supports unique mechanisms of lignocellulose conversion. Proc Natl Acad Sci U S A. 2009;106:1954–9.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhao Y, Li J, Chen Y, Hang H. Response to oxidative stress of Coriolus versicolor induced by exogenous hydrogen peroxide and paraquat. Ann Microbiol. 2009;59:221–7.View ArticleGoogle Scholar
- Daniel G, Volc J, Filonova L, Plíhal O, Kubátová E, Halada P. Characteristics of Gloeophyllum trabeum alcohol oxidase, an extracellular source of H2O2 in brown rot decay of wood. Appl Environ Microbiol. 2007;73:6241–53.PubMed CentralView ArticlePubMedGoogle Scholar
- Vanden Wymelenberg A, Gaskell J, Mozuch M, Sabat G, Ralph J, Skyba O, et al. Comparative transcriptome and secretome analysis of wood decay fungi Postia placenta and Phanerochaete chrysosporium. Appl Environ Microbiol. 2010;76:3599–610.PubMed CentralView ArticlePubMedGoogle Scholar
- Fekete FA, Chandhoke V, Jellison J. Iron-binding compounds produced by wood-decaying basidiomycetes. Appl Environ Microbiol. 1989;55:2720–2.PubMed CentralPubMedGoogle Scholar
- Knight SAB, Vilaire G, Lesuisse E, Dancis A. Iron acquisition from transferrin by Candida albicans depends on the reductive pathway. Infect Immun. 2005;73:5482–92.PubMed CentralView ArticlePubMedGoogle Scholar
- Oberegger H, Schoeser M, Zadra I, Abt B, Haas H. SREA is involved in regulation of siderophore biosynthesis, utilization and uptake in Aspergillus nidulans. Mol Microbiol. 2001;41:1077–89.View ArticlePubMedGoogle Scholar
- Merhej J, Richard-Forget F, Barreau C. The pH regulatory factor Pac1 regulates Tri gene expression and trichothecene production in Fusarium graminearum. Fungal Genet Biol. 2011;48:275–84.View ArticlePubMedGoogle Scholar
- White TJ, Bruns TD, Lee S, Taylor J. Amplification and direct sequencing of fungal ribosomal DNA for phylogenetics. In: Innis MA, Gelfand DH, Sninsky JJ, White TJ, editors. PCR Protoc a Guid to methods Appl. New York: Academic; 1990. p. 315–22.Google Scholar
- Moretti S, Armougom F, Wallace IM, Higgins DG, Jongeneel CV, Notredame C. The M-Coffee web server: a meta-method for computing multiple sequence alignments by combining alternative alignment methods. Nucleic Acids Res. 2007;35(Web Server issue):W645–8.PubMed CentralView ArticlePubMedGoogle Scholar
- Huelsenbeck JP, Ronquist F. MRBAYES: bayesian inference of phylogenetic trees. Bioinformatics. 2001;17:754–5.View ArticlePubMedGoogle Scholar
- Ronquist F, Huelsenbeck J, Teslenko M. Draft MrBayes Version 3.2 Manual: Tutorials and Model Summaries. 2011. p. 103.Google Scholar
- Moslem MA, Bahkali AH, Abd-Elsalam KA, Wit PJGM. An efficient method for DNA extraction from Cladosporioid fungi. Genet Mol Res. 2010;9:2283–91.View ArticlePubMedGoogle Scholar
- Zerbino DR, Birney E. Velvet: algorithms for de novo short read assembly using de Bruijn graphs. Genome Res. 2008;18:821–9.PubMed CentralView ArticlePubMedGoogle Scholar
- Boetzer M, Henkel CV, Jansen HJ, Butler D, Pirovano W. Scaffolding pre-assembled contigs using SSPACE. Bioinformatics. 2011;27:578–9.View ArticlePubMedGoogle Scholar
- Boetzer M, Pirovano W. Toward almost closed genomes with GapFiller. Genome Biol. 2012;13:R56.PubMed CentralView ArticlePubMedGoogle Scholar
- Darling ACE, Mau B, Blattner FR, Perna NT. Mauve: multiple alignment of conserved genomic sequence with rearrangements. Genome Res. 2004;14:1394–403.PubMed CentralView ArticlePubMedGoogle Scholar
- Lomsadze A, Ter-Hovhannisyan V, Chernoff YO, Borodovsky M. Gene identification in novel eukaryotic genomes by self-training algorithm. Nucleic Acids Res. 2005;33:6494–506.PubMed CentralView ArticlePubMedGoogle Scholar
- Conesa A, Götz S, García-Gómez JM, Terol J, Talón M, Robles M. Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005;21:3674–6.View ArticlePubMedGoogle Scholar
- Tatusov RL, Fedorova ND, Jackson JD, Jacobs AR, Kiryutin B, Koonin EV, et al. The COG database: an updated version includes eukaryotes. BMC Bioinformatics. 2003;4:41.PubMed CentralView ArticlePubMedGoogle Scholar
- Quevillon E, Silventoinen V, Pillai S, Harte N, Mulder N, Apweiler R, et al. InterProScan: protein domains identifier. Nucleic Acids Res. 2005;33(Web Server issue):W116–20.PubMed CentralView ArticlePubMedGoogle Scholar
- Lagesen K, Hallin P, Rødland EA, Staerfeldt H-H, Rognes T, Ussery DW. RNAmmer: consistent and rapid annotation of ribosomal RNA genes. Nucleic Acids Res. 2007;35:3100–8.PubMed CentralView ArticlePubMedGoogle Scholar
- Lowe TM, Eddy SR. tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 1997;25:955–64.PubMed CentralView ArticlePubMedGoogle Scholar
- Yin Y, Mao X, Yang J, Chen X, Mao F, Xu Y. dbCAN: a web resource for automated carbohydrate-active enzyme annotation. Nucleic Acids Res. 2012;40(Web Server issue):W445–51.PubMed CentralView ArticlePubMedGoogle Scholar
- Rawlings ND, Barrett AJ, Bateman A. MEROPS: the database of proteolytic enzymes, their substrates and inhibitors. Nucleic Acids Res. 2012;40(Database issue):D343–50.PubMed CentralView ArticlePubMedGoogle Scholar
- Petersen TN, Brunak S, von Heijne G, Nielsen H. SignalP 4.0: discriminating signal peptides from transmembrane regions. Nat Methods. 2011;8:785–6.View ArticlePubMedGoogle Scholar
- Krogh A, Larsson B, von Heijne G, Sonnhammer EL. Predicting transmembrane protein topology with a hidden Markov model: application to complete genomes. J Mol Biol. 2001;305:567–80.View ArticlePubMedGoogle Scholar
- Khaldi N, Seifuddin FT, Turner G, Haft D, Nierman WC, Wolfe KH, et al. SMURF: genomic mapping of fungal secondary metabolite clusters. Fungal Genet Biol. 2010;47:736–41.PubMed CentralView ArticlePubMedGoogle Scholar
- Winnenburg R, Urban M, Beacham A, Baldwin TK, Holland S, Lindeberg M, et al. PHI-base update: additions to the pathogen host interaction database. Nucleic Acids Res. 2008;36(Database issue):D572–6.PubMed CentralPubMedGoogle Scholar
- Karányi Z, Holb I, Hornok L, Pócsi I, Miskei M. FSRD: fungal stress response database. Database (Oxford). 2013;2013:bat037.View ArticleGoogle Scholar
- TransposonPSI: An Application of PSI-Blast to Mine (Retro-)Transposon ORF Homologies [http://transposonpsi.sourceforge.net] Accessed 23 Aug 2013.
- Li L, Stoeckert CJ, Roos DS. OrthoMCL: identification of ortholog groups for eukaryotic genomes. Genome Res. 2003;13:2178–89.PubMed CentralView ArticlePubMedGoogle Scholar
- Van Dongen S: Graph Clustering by Flow Simulation. University of Utrecht: The Netherlands; 2000.Google Scholar
- Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, et al. Clustal W and Clustal X version 2.0. Bioinformatics. 2007;23:2947–8.View ArticlePubMedGoogle Scholar
- Abascal F, Zardoya R, Posada D. ProtTest: selection of best-fit models of protein evolution. Bioinformatics. 2005;21:2104–5.View ArticlePubMedGoogle Scholar
- Stamatakis A. RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics. 2006;22:2688–90.View ArticlePubMedGoogle Scholar