Evolutionary analysis of hydrophobin gene family in two wood-degrading basidiomycetes, Phlebia brevispora and Heterobasidion annosums.l.
- Anthony C Mgbeahuruike†1Email author,
- Andriy Kovalchuk†1,
- Hongxin Chen1,
- Wimal Ubhayasekera2 and
- Fred O Asiegbu1
© Mgbeahuruike et al.; licensee BioMed Central Ltd. 2013
Received: 23 April 2013
Accepted: 11 October 2013
Published: 4 November 2013
Hydrophobins are small secreted cysteine-rich proteins that play diverse roles during different phases of fungal life cycle. In basidiomycetes, hydrophobin-encoding genes often form large multigene families with up to 40 members. The evolutionary forces driving hydrophobin gene expansion and diversification in basidiomycetes are poorly understood. The functional roles of individual genes within such gene families also remain unclear. The relationship between the hydrophobin gene number, the genome size and the lifestyle of respective fungal species has not yet been thoroughly investigated. Here, we present results of our survey of hydrophobin gene families in two species of wood-degrading basidiomycetes, Phlebia brevispora and Heterobasidion annosum s.l. We have also investigated the regulatory pattern of hydrophobin-encoding genes from H. annosum s.s. during saprotrophic growth on pine wood as well as on culture filtrate from Phlebiopsis gigantea using micro-arrays. These data are supplemented by results of the protein structure modeling for a representative set of hydrophobins.
We have identified hydrophobin genes from the genomes of two wood-degrading species of basidiomycetes, Heterobasidion irregulare, representing one of the microspecies within the aggregate H. annosum s.l., and Phlebia brevispora. Although a high number of hydrophobin-encoding genes were observed in H. irregulare (16 copies), a remarkable expansion of these genes was recorded in P. brevispora (26 copies). A significant expansion of hydrophobin-encoding genes in other analyzed basidiomycetes was also documented (1–40 copies), whereas contraction through gene loss was observed among the analyzed ascomycetes (1–11 copies). Our phylogenetic analysis confirmed the important role of gene duplication events in the evolution of hydrophobins in basidiomycetes. Increased number of hydrophobin-encoding genes appears to have been linked to the species’ ecological strategy, with the non-pathogenic fungi having increased numbers of hydrophobins compared with their pathogenic counterparts. However, there was no significant relationship between the number of hydrophobin-encoding genes and genome size. Furthermore, our results revealed significant differences in the expression levels of the 16 H. annosum s.s. hydrophobin-encoding genes which suggest possible differences in their regulatory patterns.
A considerable expansion of the hydrophobin-encoding genes in basidiomycetes has been observed. The distribution and number of hydrophobin-encoding genes in the analyzed species may be connected to their ecological preferences. Results of our analysis also have shown that H. annosum s.l. hydrophobin-encoding genes may be under positive selection. Our gene expression analysis revealed differential expression of H. annosum s.s. hydrophobin genes under different growth conditions, indicating their possible functional diversification.
KeywordsHydrophobins Gene evolution Basidiomycetes Ascomycetes Pathogen
Hydrophobins are surface-active proteins produced by filamentous fungi [1, 2]. They are small secreted proteins with eight cysteine residues arranged in a strictly conserved motif [3, 4]. The cysteine residues form four disulfide bridges connecting beta strands and stabilizing the protein structure. Another important feature of hydrophobins is their ability to self-assemble into amphiphilic films at hydrophilic/hydrophobic interfaces [5–7]. Hydrophobins have been reported in filamentous fungi belonging to the phyla Ascomycota and Basidiomycota . They are often secreted extracellularly but can also be found inside fungal structures such as fruiting bodies and hyphae . Based on solubility and sequence characteristics, hydrophobins can be classified into two major classes: Class I and Class II [1, 10]. Although the two classes are similar in many ways, class I hydrophobins have larger size and more diversity in amino acid sequence than class II [1, 10]. Hydrophobins are expressed at different stages of fungal life cycle: sporulation, fruiting body formation and during growth of vegetative hyphae . Studies have shown that hydrophobins play important role in fungal pathogenesis where they act as virulence factors to enhance fungal infection [11–15]. They have also been reported to be involved in the attachment of fungal structures [16, 17] and the emergence of aerial hyphae from submerged conditions [18–20]. Hydrophobins have been implicated in diverse fungal interactions such as symbiosis , mycorrhiza formation , and antagonistic interactions [20, 23]. Evidence of hydrophobin involvement in cell wall assembly during pathogenic interactions where the monomers act as elicitors and toxins have been reported . Recent gene expression studies have shown that P. gigantea hydrophobin encoding genes 1 and 2 (Pgh1 and Pgh2) are highly transcribed in the interaction zone between the biological control agent P. gigantea and the tree pathogen H. annosum s.l. [20, 23]. However, the actual roles of these genes in the interaction are not known. Previous studies have demonstrated a high level of sequence divergence in P. gigantea hydrophobin-encoding genes Pgh1 and Pgh2 and the hydrophobins from the pathogenic fungus H. irregulare Hah1 and Hah2. These observations raised crucial questions on the evolutionary forces driving the rapid differentiation of this gene family. Available data also indicate a considerable amount of variation in the numbers of hydrophobin-encoding genes across other fungal taxa, ranging from 1 gene in Acremonium alcalophilum to 40 genes in Trametes versicolor. This may suggest that hydrophobin genes could be under a dynamic evolutionary process across most fungal taxa.
Gene duplication is an important evolutionary process that plays a crucial role in an organism’s complexity, adaptation and diversification to closely related strains and species [25, 26]. Reports have shown that duplication of genes results in functional diversification and gene expression patterns observed in different fungi and other organisms . Paralogous genes resulting from duplication events create genetic redundancy, which may be vulnerable to selection pressure [25, 27]. Genome-wide analysis of gene duplication has shown that this biological phenomenon occurs at a very high rate. However, the fate of duplicated genes and the forces driving their fixation and divergence still remain unknown . There has been increasing number of evidences showing that mutated genes with deleterious effects are purged from the genome through purifying selection whereas copies with enhanced functions are fixed in the population through positive or diversifying selection [28, 29]. Furthermore, evidence of duplication and losses among stress-related genes has been documented in the literature whereas growth-related genes have been shown to be selected against change in the copy number . Gene contraction through gene loss and expansion through duplication are common processes in gene family evolution and it has been documented in chitinases  and glycosyl hydrolases (GH28)  gene families. However, the mechanisms driving these evolutionary processes in fungal hydrophobins are poorly understood.
In this study, we investigated hydrophobin gene family evolution in several diverse fungal groups, we also analyzed recombination events in H. annosum s.l. by examining the ratio of non-synonymous (dN) to synonymous substitutions per site (dS). We tested correlation between number of hydrophobin-encoding genes, overall genome size and their ecological strategy. The transcript abundance of sixteen hydrophobin genes from H. annosum s.s. during saprotrophic growth on pine wood as well as on culture filtrate from P. gigantea was further evaluated using micro-array. The results from this study have further highlighted the possible involvement of hydrophobin genes in fungal ecological lifestyle.
Hydrophobin sequence identification and alignment
Sequence alignment of selected fungal hydrophobins
Spacing between conserved cysteine residues in different classes of hydrophobins
Numbers of amino acid residues between conserved cysteine residues
Class I (basidiomycetes)
Class I (ascomycetes)
T. terrestris protein 159967
U. maydis protein 5010
Domain structure, hydropathy pattern and homology modeling of selected hydrophobins
Phylogenetic analysis of the fungal hydrophobins and ecological strategy
Basidiomycetes and ascomycetes
In order to understand the evolutionary relationships of the fungal hydrophobins from both ascomycetes and basidiomycetes screened in this study (Additional file 1: Table S1), species based trees were reconstructed using the NJ method with JTT matrix-based model . Two distinct separations along the two classes of hydrophobins were evident, all the class I hydrophobins clustered in clade A whereas class II members grouped together in clade B (Additional file 2: Figure S1). Other subclades such as C, D, E, F and G were evident. Subclade C is a mixed group containing class I proteins from both ascomycetes and basidiomycetes. Most of the basidiomycetes in this group are from the order Polyporales except few species like Auricularia delicata and Schizophyllum commune which were from the orders Auriculariales and Agaricales, respectively. It was also interesting to see that all the hydrophobins from A. delicata, a basidiomycete, are closer related to the class I proteins from the ascomycetes than to the remaining members of class I hydrophobins from basidiomycetes (Additional file 2: Figure S1, subclade C). Hydrophobins from U. maydis and T. terrestris which showed a deviation from the general consensus of classes I and II (49 aas at C3/C4 and 5 aa at C3/C4) respectively, were nested within the Class I hydrophobins (Additional file 2: Figure S1) in our phylogenetic analysis and thus, most likely, they are just two extreme examples of the variation in length of C3/C4 region in the Class I hydrophobins. Furthermore, the hydrophobin from T. terrestris with unusually short stretch of amino acids between cysteine residues C3 and C4 appeared within group C. Subclade D contains white rot fungi of the order Polyporales except Wolfiporia cocos, a brown rot fungus. Most of the hydrophobins from P. brevispora clustered in this group. However sequences of hydrophobins from H. irregulare could not be resolved into a distinct clade. Groups E, F and G consisted of class I hydrophobins from basidiomycetes of different systematic positions and ecological lifestyles (brown rots, white rots, mycorrhiza-formers and saprotrophs). Another interesting observation from the phylogenetic grouping is that many ascomycetes such as Magnaporthe grisea, Trichoderma virens, Neurospora discreta, Neurospora crassa and Neurospora tetrasperma appeared to have both Classes I and II proteins and the two classes separated into different clades (Additional file 2: Figure S1).
A similar trend was observed for the phylogenetic tree involving only hydrophobins from basidiomycetes. Seven main clades and two smaller clades were formed (Additional file 3: Figure S2). Members of clades H, I and J were mostly hydrophobins from fungi of the order Polyporales, all the fungi in this group except W. cocos and Punctularia strigosozonata are white rotters. Most of the sequences of hydrophobins from P. brevispora were resolved in clade J. However, sequences of hydrophobins from H. irregulare have not formed a distinct clade. Clade K consists of hydrophobins from fungi of diverse systematic positions (orders Agaricales, Polyporales, Russulales, Corticiales and Boletales) and ecological strategies. Clades L and P are relatively smaller clades with hydrophobins from fungi representing different orders and lifestyles. In clade M, sequences of hydrophobins from S. commune, a white rot fungus dominated the group with hydrophobins from A. delicata. Group N is a relatively small group comprising of hydrophobins from fungal species of the orders Agaricales, Boletales and Dacrymycetales. Fungal species in this group have mixed lifestyle, while some are saprotrophic in nature (Coprinopsis cinerea), others are brown rotters. In addition, some sequences of hydrophobins from the mycorrhizal fungus, Laccaria bicolor were found in this clade. In clade O, all the hydrophobin sequences separated into fungal species with brown rotting habits. The fungal species (Coniophora puteana and Serpula lacrymans) in this group are from the order Boletales.
Distribution, genome size and hydrophobin gene family evolution in ascomycetes and basidiomycetes
Gene clusters and inventory of hydrophobins in P. brevispora and H. irregularehydrophobins
The 26 class I hydrophobin-encoding-genes found in the genome of P. brevispora, were arranged in a relatively clustered pattern (Additional file 6: Table S2). A good number of the proteins (6) were located in scaffold 12, whereas 4 proteins each where found in scaffolds 11 and 14 respectively. Scaffold 23 had 3 proteins, while scaffolds 9, 14 and 19 contained 2 proteins each. Other scaffolds, 75, 30 and 38 had one protein each.
In H. irregulare, 24 putative class I hydrophobin encoding gene sequences were identified, 16 appeared to be functional proteins whereas 8 were probable pseudogenes. A similar clustering pattern of hydrophobin-encoding genes was observed in the genome of H. irregulare (Additional file 7: Table S3), most of the genes (15) clustered in scaffold 11. Scaffold 6 had 3 genes, whereas 2 genes were found in scaffold 3. Scaffolds 4, 8 and 9 had 1 gene each.
Selection tests for H. irregularehydrophobins
Selection test for hydrophobins genes from H. irregulare
Number of hydrophobins
Tajima’s D test
Non-synonymous substitutions (dN)
Synonymous substitutions (dS)
Microarray analysis of hydrophobin expression in H. annosums.s
Microarray expression a analysis of hydrophobin encoding genes from H. annosum s.s. during growth on wood and growth on culture filtrate from P. gigantea
Pg/Ha fold changesb
Bark/control fold changesb
Sapwood/control fold changesb
Heartwood/control fold changesb
Hydrophobins are surface-active proteins with diverse roles in fungal life cycle [9, 11–15]. Hydrophobins are grouped into class I and class II based on the following criteria: minimum of eight cysteines arranged in a conserved pattern, solubility of the formed aggregates, hydropathy pattern and the differences in the spacing of amino acids between the cysteine residues [1, 36]. Although we did not carry out solubility test in this study, a clear separation of the hydrophobins screened in this study into class I and II was observed. Although H. irregulare and P. brevispora have different ecological lifestyles, the domain structure of hydrophobins from these two fungi did not differ considerably, an indication that the differences in nutritional lifestyle could not be explained by their domain structures. Separation of hydrophobins into various classes based on hydropathy pattern has been reported in other studies .
The most prominent feature of the phylogenetic grouping of hydrophobins from all the fungal species is the separation of the sequences along the two classes of hydrophobins (Additional file 2: Figure S1). Similar results based on phylogenetic grouping of Class I and II into the two major fungal taxa have been reported [3, 37]. Some fungal species such as M. grisea, T. virens, N. discreta, N. crassa and N. tetrasperma have both class I and class II proteins which were resolved into distinct clades (Additional file 2: Figure S1). Similar findings have been reported in other studies where both class I and class II hydrophobins were identified in a single fungal species [1, 38]. Class I hydrophobins were common among the members of the basidiomycetes, while both class I and II were present in the ascomycetes. This is in line with findings from other studies [3, 37]. The absence of class II proteins in basidiomycetes may indicate that only the class I proteins may be important in basidiomycetes for fruit body formation. However, despite the high level of sequence conservation between the hydrophobins from P. brevispora and H. annosum s.l., the two fungi appeared to be phylogenetically distant from each other (Additional file 2: Figure S1 and Additional file 3: Figure S2).
Numerous clusters of hydrophobins were found in our analysis; paralogous genes with high bootstrap support were also nested together in the generated trees. The presence of these clusters provides additional evidence in support of our hypothesis about important role of gene duplication in the evolution of hydrophobin-encoding genes in fungi. Paralogous sequences of hydrophobins from P. brevispora clustered together and they appeared to be phylogenetically closer to hydrophobins from W. cocos, B. adusta and F. mediterranea (Additional file 2: Figure S1 and Additional file 3: Figure S2). Sequences of hydrophobins from H. irregulare were not clearly resolved into distinct clades. It was interesting to note that hydrophobins from A. delicata, a basidiomycete were closer relatives to the class I proteins from ascomycetes than members of the same class basidiomycetes (Additional file 2: Figure S1). Hydrophobins of A. delicata differ from hydrophobins of most remaining basidiomycetes in our analysis by having shorter region between the conserved cysteines C3 and C4 (26–29 amino residues versus 32–33). This structural feature might influence the results of the phylogenetic analysis. At the same time, Auricularia is the most basal member of basidiomycetes in our analysis, and it is possible that hydrophobins of Auricularia retained the similarity to the hypothetical ‘ancestral’ type of hydrophobins that were present in the common ancestor of ascomycetes and basidiomycetes, whereas in evolutionary more advanced species of basidiomycetes they are more diverged from that type, e.g. by having longer stretch of amino acids between C3 and C4 positions. The phylogenetic relationships between class I hydrophobins from ascomycetes and the same class of proteins from the basidiomycetes may suggest that the ancestor of these genes may have been formed very early during evolution. Furthermore, it is possible that the gene products may have a role in the life style or ecology of the fungal species. Based on this result, a complete survey of genome sizes and copy numbers of hydrophobin-encoding genes was made on several fungal species with varied ecological strategies (Additional file 4: Figure S3 and Additional file 5: Figure S4). A considerable variation in the number of hydrophobin-encoding genes exists in all the fungi screened, ranging from 1 in A. alcalophilum to 40 in T. versicolor. In H. irregulare 13 hydrophobins have previously been reported . In our analysis, 16 functional class I proteins and 8 putative pseudogene were identified. This difference in the number of hydrophobin-encoding genes observed in our study could be a result of automatic annotation problems. However, our microarray studies identified only the 16 functional proteins. In P. brevispora, 26 class I hydrophobin-encoding genes have been identified. Accumulation of several copies of duplicate hydrophobin genes from a single copy of ancestral sequence may have resulted in the variation in numbers of hydrophobin genes in the different fungi. An interesting finding in this work is the absence of hydrophobins among the members of hemiascomycetous yeasts, an indication that species like S. cerevisiae, Pichia stipitis, Hansenula polymorpha and Wickerhamomyces anomalus with yeast-like or monocentric growth pattern (non-filamentous growth pattern) may not need hydrophobins. In addition, the complete absence of hydrophobins in the Pucciniales may suggest that hydrophobin genes may have undergone relaxed selection after evolution. It is also possible that hydrophobins were completely lost during evolution. Among the basidiomycetes screened in this study, a significant expansion of the gene was observed (P = 0.002) whereas ascomycetes appeared to have witnessed a massive contraction of the gene during evolution. The expansion of hydrophobin-encoding genes among the members of Basidiomycota may not be unconnected with fruit body formation . It is also possible that the increase in the number of hydrophobin-encoding genes may have been positively selected for in basidiomycetes, but not in ascomycetes. Relationship between the number of hydrophobin-encoding genes and ecological strategy (pathogenic and non-pathogenic life style) was further investigated. There was a well-supported link between the number of hydrophobin genes and ecological strategy with the non-pathogenic fungi having higher numbers of hydrophobins than the pathogenic ones (P = 0.0001), this may suggest that although hydrophobins may be needed in fungal pathogenesis [11–15], higher numbers of the gene may have more ecological role in non-pathogenic conditions such as symbiosis , mycorrhiza formation  and interspecific fungi-fungi interactions [20, 23], fruit body formation  and emergence of hyphal structures [18–20]. This result is in contrast with other studies that have shown massive expansion of gene families in pathogens as compared with non-pathogenic relatives [31, 32]. Differences between ecological strategies and genome size were also tested. Genome size had no significant effect on ecological strategy (P = 0.1).
Evolutionary forces operating at a genomic scale may have some influence on gene family expansion or contractions, implying that the fungal species with larger genome sizes may have a correspondingly higher number of hydrophobins. However, there was no significant relationship between genome size and the number of hydrophobin-encoding genes in this study. To understand the evolutionary forces driving hydrophobin gene family evolution in H. irregulare, a selection test was carried out. A high dN/dS ratio (> 1) was observed in H. irregulare hydrophobins, an indication that hydrophobin genes may be under positive selection. It is also possible that the evolution of these genes could be a result of recombination and duplication events. Evidence of duplication events have been reported in hydrophobins from P. gigantea[20, 23] and H. irregulare. These finding may therefore suggest that this group of proteins may have evolved through the so called birth and death model . In the birth and death model, new gene copies evolve through duplications followed by diversification due to accumulation of spontaneous mutations; new duplicates with vital functions are retained in the genome while those with deleterious effects are purged from the genome through purifying selection. We identified 16 transcripts of class I hydrophobin encoding genes that were differentially regulated during growth on bark, heartwood and sapwood as well as during growth on culture filtrate produced by P. gigantea. Furthermore, the differences in expression pattern of these class I proteins during saprotrophic growth may suggest that different hydrophobin genes are employed by H. annosum s.s. during growth on diverse wood components. In addition, the low transcript abundance observed during growth in the submerged medium containing culture filtrate from P. gigantea may suggest that the induction of these genes could be inhibited due to the presence of the secreted molecules from the biocontrol agent (P. gigantea). P. gigantea is the biological control fungus for the control of H. annosum s.l. infection. It is possible that one of the mechanisms for action of the biocontrol fungus against the pathogen is by repression of genes such as hydrophobins through its secreted molecules or metabolites. Other studies have reported the upregulation of hydrophobin encoding genes at the zone of interaction between the biocontrol fungus and the pathogen [20, 23]. Due to lack of protein structure of hydrophobins from H. annosum s.l., hydrophobins from selected fungal species were modeled alongside hydrophobins from P. gigantea. The models revealed the surface patches of hydrophobic residues (Figure 3), which are possibly important in the formation of amphiphilic membranes as reported earlier . The conserved cysteine residues confirm probably a preserved structural feature of the protein.
We have surveyed the distribution and evolution of hydrophobin genes in P. brevispora and H. irregulare as well as in other fungi. We have also examined the relationships between the number of hydrophobin-encoding genes and ecological strategy in the examined fungal species. From our results, hydrophobin genes have witnessed a considerable expansion in P. brevispora as well as in other basidiomycetes while contraction of the same gene family has occurred in the ascomycetes. In addition, although different numbers of hydrophobins have been reported in P. brevispora, H. annosum s.l. and other fungal species, each gene has different regulatory pattern in the pathogenic fungus (H. annosum s.s.) during growth on wood and in culture filtrate produced by the antagonistic P. gigantea. The results from this study have also given some insights on some of the factors underlying the ecological habits of the fungi screened in this study.
Annotated sequences of hydrophobins from H. irregulare and P. brevispora were used as queries to search for other hydrophobin sequences of some randomly selected fungal species, using BlastP . Blast searches were performed at the fungal genomic platform of the Joint Genome Institute (JGI) (http://www.jgi.doe.gov/). The identified ORFs were used as queries to search for all possible hydrophobin proteins in the selected fungal species. Repeated blast searches were carried out until no more hydrophobin proteins were found. Sequences with E-values below 10-5 were selected for further analysis. Due to the reasonably lower number of hits generated by the blast results, filtered model was further used as a criterion to download all the hydrophobin proteins of the selected fungal taxa. Sequences without N- or C-terminal parts were corrected. Also sequences with regions of unspliced introns were corrected and aligned using MUSCLE alignment tool implemented in Molecular Evolutionary Genetic Analysis software (MEGA 5.0). MUSCLE was used because it gives a better accuracy and is computationally stronger than CLUSTAL alignment . The aligned protein sequences were viewed with the Biological sequence alignment editor (Bioedit), Windows 95/98/NT/2 K/XP . Alignments were curated manually, and all ambiguous positions were removed. The sequences were further filtered to remove pseudogenes, or sequences shorter than 50 residues in length, or missing the hydrophobin domains as defined by the InterproScan Tool  or having a different gene ontology (GO) from hydrophobins and realigned for further analysis. The signal peptides were predicted using SignalP 3.0 software. The protein sequences were classified into class I and class II using Kyte-Doolittle hydropathy plot, Version 2.0u66 in Windows 9.0  and published consensus sequence information for class I  and class II  respectively.
Species tree reconstruction
A total of 335 protein sequences from 41 fungal species (Additional file 1: Table S1) were used for the analysis. The protein sequences were further divided into the two major fungal taxa, basidiomycetes and ascomycetes. Species-based trees were reconstructed for both groups using the Neighbor-joining (NJ) method in MEGA. The evolutionary distances were calculated using the Jones Taylor–Thorton (JTT) + gamma matrix-based  taking into account rate heterogeneity among sites. The rate variation among sites was calculated with a gamma distribution with a default parameter of 5. The reliability of internal branches was evaluated using 1000 bootstrap replications .
Species ecological strategy, genome size and hydrophobin distribution
Information on the ecological strategy of each fungal species was obtained from already published data. In parallel, a complete survey of hydrophobin distribution and genome sizes of all the fungi species used in the study was obtained from JGI (http://www.jgi.doe.gov/) (Additional file 1: Table S1). The fungal species were divided into basidiomycetes and ascomycetes, and were further subdivided according to their ecological strategies, pathogens and non-pathogens (Additional file 1: Table S1). The selected fungal species (Additional file 1: Table S1) were analyzed for gene contraction through gene loss and gene expansion through duplication. The relationship between ecological strategy and the copy number of hydrophobin-encoding genes within the selected set of species was also examined (Additional file 1: Table S1). Further analysis to determine the relationship between genome size and hydrophobin gene distribution was evaluated among the fungi in this group.
Microarray expression analysis of hydrophobin genes from the pathogenic fungus, H. annosum s.s. during growth on culture filtrate of P. giganteaand saprotrophic growth on pine wood
Saprotrophic growth of H. annosum s.s.on wood bark, sapwood and heartwood
Wood discs from Scots Pine (Pinus sylvestris) were separated into bark, sapwood and heartwood. Each wood component was grounded into small particles of sizes 0.5-1 mm for 15 min at 590 rpm using a ball grinding mill (Fritsch Pulverisette, Germany) and 8 g of each wood material was weighed in a flask and autoclaved for 20 min. The wood tissues were allowed to cool for 20 min, 8 ml sterile low nitrogen medium (NH4NO3 0.6 g/L, K2HPO4 0.4 g/L, KH2PO4 0.5 g/L, MgSO4·7H20 0.4 g/L) was added to each flask followed by the addition of 8 ml of sterile distilled water ensuring that comparable moisture levels was maintained. This was followed by inoculation with three pieces of 1 × 1 cm agar plugs of malt extract agar containing H. annosum s.s. hyphae (isolate FP5, obtained from Kari Korhonen, Finish Forest Research Institute (METLA), Vantaa Finland). The plugs were put into each flask, mixed gently to allow the agar plugs to be covered by the wood tissues. Cultures were incubated at 20°C and harvested after 3 months. Harvested mycelia and wood tissues were frozen in liquid nitrogen and stored at −80°C until RNA extraction. There were 3 biological replicates for each sample.
Growth of H. annosum s.s. in culture filtrates of P. gigantea
Three Erlenmayer flask (300 ml) containing liquid malt extract (100 ml) each, were inoculated with 3 mm agar plugs of the commercial isolate of P. gigantea (Rotstop®), courtesy of Kari Korhonen (Finish Forest Research Institute (METLA), Vantaa Finland). Cultures were incubated at 20°C for 10 days. The cultures were filtered to get rid of the fungal mycelia using sterile flask and filter paper. The liquid filtrate was incubated overnight at 70°C to kill any remaining particle of P. gigantea mycelia. Freshly growing agar plugs (3 mm) of H. annosum s.s. were inoculated in each flask and cultures were incubated at 20°C. In parallel, fresh liquid medium of malt extract were inoculated with 3 mm freshly growing agar plugs from H. annosum s.s. to act as positive control and incubated at 20°C. Cultures were harvested after 10 days post inoculation (d.p.i) and the harvested mycelia were frozen in liquid nitrogen until further processing. There were 3 biological replicates for each experiment.
RNA processing and microarray
RNA was extracted from triplicate cultures of each sample using the method by Chang  with some modifications. RNA was purified by using the RNeasy® MinElute Cleanup kit (QIAGEN) according to the protocol. RNA integrity was assessed with RNA 6000 Nano kit using an Agilent Bioanalyzer (Agilent, CA). RNA concentration was measured using NanoDrop ND-1000 Spectrophotometer and the purity of the samples was estimated by the OD ratios (A260/A280, ranging within 1.8–2.2). The RNA samples were DNase treated to remove any potential DNA contamination using DNaseI according to the recommendations from the manufacturer (Fermentas, Canada). cDNA was synthesized by using the TransPlex® Complete Whole Transcriptome Amplification Kit according to the manufacturer’s protocol (SIGMA). Microarray analysis was carried out using the Nimblegene protocol (http://www.nimblegen.com/).
Pearson’s partial correlation coefficient was used to test if the number of hydrophobins in P. brevispora or H. irregulare and hydrophobins from other fungal species correlate with their genome sizes. General linear model (GLM) procedure in SAS was used to test the hypothesis that the number of hydrophobin-encoding genes and genome size differed between basidiomycetes and ascomycetes, and between pathogens and non-phytopathogens. To understand the evolutionary forces driving hydrophobin gene evolution in H. irregulare, recombinations events were tested using the method described in ; also the ratio of non-synonymous substitution per site (dN) to synonymous substitution per site (dS) in 13 paralogous sequences of hydrophobin-encoding genes was calculated. Tajima’s D test statistics  was also applied on the hydrophobin-encoding genes from H. irregulare. For the microarray data, the mean expression and fold changes were calculated with FDR (Benjamini Hochberg) multiple testing corrections using ArraySTAR software (3801 Regent Street Madison, WI53705, USA). Student t-test was also used to determine differences in mean between samples.
Homology modeling of hydrophobins
Similar sequences to the hydrophobins from P. gigantea and H. irregulare as well as hydrophobins from Verticillium dahliae, T. terrestris, U. maydis and L. bicolor were located in the entries of Protein Data Bank (PDB) using PSI-BLAST search  and aligned using CLUSTAL W . The best predictions were selected as the templates and these structures were obtained from the PDB , then superimposed and compared with the programs LSQMAN  and O . The best pair-wise alignments with the relevant parts of the structures were used to generate homology models of the hydrophobins from P. gigantea, H. irregulare, V. dahliae, T. terrestris, U. maydis and L. bicolor with thioredoxin reductase from Drosophila melanogaster (PDB entry 3DGH; identity 35%), transcriptional regulator BT_p548217 from Bacteroides thetaiotaomicron (PDB entry 2K9Q; identity 35%), Hydrophobin from Hypocrea jecorina (PDB entry 2FZ6; identity 44%) , hydrophobin from N. crassa (PDB entry 2K6A ; identity 48%), Dsl1p subunit from Saccharomyces cerevisiae (PDB entry 3ETU ; identity 31%) and glyceraldehyde-3-phosphate dehydrogenase from Trypanosoma cruzi (PDB entry 3DMT ; identity 30%) structures respectively as templates in the program SOD . The models were adjusted in O, using rotamers that would improve packing in the interior of the protein, and accounting for insertions and deletions in loop regions. The models are available upon request from the authors. The figure was prepared using http://www.pymol.org.
Availability of supporting data
The microarray data obtained in this work were deposited at the Gene Expression Omnibus (GEO, http://www.ncbi.nlm.nih.gov/geo/) database (accession number GSE39805 (http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE39805) and GSE41301 (http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE41301). The alignments of hydrophobin sequences used for the phylogenetic reconstructions were deposited at the TreeBASE database (http://www.treebase.org/) under accession number 14520 (http://purl.org/phylo/treebase/phylows/study/TB2:S14520).
This work was supported with funding from Swedish Research Council for Environment, Agricultural sciences and Spatial planning (FORMAS). Financial support from Academy of Finland is also gratefully acknowledged.
- Linder MB, Szilvay GR, Nakari-Setälä T, Penttila ME: Hydrophobins: the protein-amphiphiles of filamentous fungi. FEMS Microbiol Rev. 2005, 29: 877-896. 10.1016/j.femsre.2005.01.004.PubMedView ArticleGoogle Scholar
- Garbe LA, Schwarz P, Ehmer A: Beer Gushing. Handbook of Alcoholic Beverages Series, Beer A Quality perspective. Edited by: Bamforth, Charles W, Russell I, Stewart G. 2009, San Diego, London: Elsevier Ltd, 185-212.Google Scholar
- Wosten HA: Hydrophobins: multipurpose proteins. Ann Rev Microbiol. 2001, 55: 625-646. 10.1146/annurev.micro.55.1.625.View ArticleGoogle Scholar
- Elliot MA, Talbot NJ: Building filaments in the air: aerial morphogenesis in bacteria and fungi. Curr Opin Microbiol. 2004, 7: 594-601. 10.1016/j.mib.2004.10.013.PubMedView ArticleGoogle Scholar
- Hakanpaa J, Paananen A, Askolin S, Nakari-Setala T, Parkkinen T, Penttila M, Linder MB, Rouvinen J: Atomic resolution structure of the HFBII hydrophobin, a self-assembling amphiphile. J Biol Chem. 2004, 279: 534-539.PubMedView ArticleGoogle Scholar
- Hakanpää J, Szilvay GR, Kaljunen H, Maksimainen M, Linder M, Rouvinen J: Two crystal structures of Trichoderma reesei hydrophobin HFBI--the structure of a protein amphiphile with and without detergent interaction. Proc Natl Acad Sci U S A. 2006, 15 (9): 2129-2140.Google Scholar
- Kallio JM, Linder MB, Rouvinen J: Crystal structures of hydrophobin HFBII in the presence of detergent implicate the formation of fibrils and monolayer films. J Biol Chem. 2007, 282 (39): 28733-28739. 10.1074/jbc.M704238200.PubMedView ArticleGoogle Scholar
- Scholtmeijer K, Wessels JGH, Wösten HAB: Fungal hydrophobins in medical and technical applications. Appl Microbiol Biotechnol. 2001, 56: 1-8. 10.1007/s002530100632.PubMedView ArticleGoogle Scholar
- Linder MB: Hydrophobins: proteins that self assemble at interfaces. Curr Opin Coll Inter Sci. 2009, 14: 356-363. 10.1016/j.cocis.2009.04.001.View ArticleGoogle Scholar
- Sarlin T, Nakari-Setälä T, Linder M, Penttilä M, Haikara A: Fungal hydrophobins as predictors of the gushing activity of malt. J Inst Brew. 2005, 111 (2): 105-111. 10.1002/j.2050-0416.2005.tb00655.x.View ArticleGoogle Scholar
- Felipe MS, Andrade RV, Petrofeza SS, Maranhao AQ, Torres FA, Albuquerque P, Arraes FB, Arruda M, Azevedo MO, Baptista AJ, Bataus LA, Borges CL, Campos EG, Cruz MR, Daher BS, Dantas A, Ferreira MA, Ghil GV, Brigid MM: Transcriptome characterization of the dimorphic and pathogenic fungus Paracoccidioides brasiliensis by EST analysis. Yeast. 2003, 20 (3): 263-271. 10.1002/yea.964.PubMedView ArticleGoogle Scholar
- Karlsson M, Stenlid J, Olson A: Two hydrophobin genes from the conifer pathogen Heterobasidion annosum are expressed in aerial hyphae. Mycologia. 2007, 99: 227-231. 10.3852/mycologia.99.2.227. 10.3852/mycologia.99.2.227PubMedView 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, Nierman WC: Genomic islands in the pathogenic filamentous fungus Aspergillus fumigatus. PLoS Genet. 2008, 4: E1000046-10.1371/journal.pgen.1000046.PubMed CentralPubMedView ArticleGoogle Scholar
- Kubicek CP, Baker S, Gamauf C, Kenerley CM, Druzhinina IS: Purifying selection and birth-and-death evolution in the class II hydrophobin gene families of the ascomycete Trichoderma/Hypocrea. BMC Evol Biol. 2008, 8: 4-10.1186/1471-2148-8-4. doi:10.1186/1471-2148-8-4PubMed CentralPubMedView ArticleGoogle Scholar
- Sharpton TJ, Stajich JE, Rounsley SD, Gardner MJ, Wortman JR, Jordar VS, Maiti R, Kodira CD, Neafsey DE, Zeng Q, Hung CY, McMahan C, Muszewska A, Grynberg M, Mandel MA, Kellner EM, Barker BM, Galgiani JN, Orbach MJ, Kirkland TN, Cole GT, Henn MR, Birren BW, Taylor JW: Comparative genomic analyses of the human fungal pathogens Coccidioides and their relatives. Genome Res. 2009, 19: 1722-1731. 10.1101/gr.087551.108.PubMed CentralPubMedView ArticleGoogle Scholar
- Tasaki Y, Ohata K, Hara T, Joh T: Three genes specifically expressed during phosphate deficiency in Pholiota nameko strain N2 encode hydrophobins. Curr Genets. 2004, 45: 19-27. 10.1007/s00294-003-0457-5.View ArticleGoogle Scholar
- Kim S, Park J, Park SY, Mitchell TK, Lee YH: Identification and analysis of in planta expressed genes of Magnaporthe oryzae. BMC Genomics. 2010, 11: 104-10.1186/1471-2164-11-104.PubMed CentralPubMedView ArticleGoogle Scholar
- Mankel A, Krause K, Kothe E: Identification of a hydrophobin gene that is developmentally regulated in the ectomycorrhizal fungus Tricholoma terreum. Appl Environ Microbiol. 2002, 68: 1408-1413. 10.1128/AEM.68.3.1408-1413.2002.PubMed CentralPubMedView ArticleGoogle Scholar
- Cho EM, Kirkland B, Holder D, Keyhani N: Phage display cDNA cloning and expression analysis of hydrophobins from the entomopathogenic fungus, Beauveria (Cordyceps) bassiana. Microbiol. 2007, 153: 3438-3447. 10.1099/mic.0.2007/008532-0.View ArticleGoogle Scholar
- Mgbeahuruike AC, Karlsson M, Asigbu FO: Differential expression of two hydrophobins (Pgh1 and Pgh2) from the biocontrol fungus Phlebiopsis gigantea. Fungal Biol. 2012, 116 (5): 620-629. 10.1016/j.funbio.2012.03.001.PubMedView ArticleGoogle Scholar
- Trembley ML, Ringli C, Honegger R: Hydrophobins DGH1, DGH2, and DGH3 in the lichen-forming basidiomycete Dictyonema glabratum. Fungal Genets Biol. 2002, 35: 247-259. 10.1006/fgbi.2001.1325. DOI:10.1006/fgbi.2001.1325View ArticleGoogle Scholar
- Martin F, Aerts A, Ahren D, Brun A, Danchin EGJ, Duchaussoy F, Gibon J, Kohler A, Lindquist E, Pereda V, Salamov A, Shapiro HJ, Wuyts J, Blaudez D, Buee M, Brokstein P, Canbaeck B, Cohen D, Grigoriev IV: The genome of Laccaria bicolor provides insights into mycorrhizal symbiosis. Nature. 2008, 452: 88-92. 10.1038/nature06556.PubMedView ArticleGoogle Scholar
- Adomas A, Eklund M, Johansson M, Asiegbu FO: Identification and analysis of differentially expressed cDNAs during non-self-competitive interaction between Phlebiopsis gigantea and Heterobasidion parviporum. FEMS Microbiol Ecol. 2006, 57: 26-39. 10.1111/j.1574-6941.2006.00094.x.PubMedView ArticleGoogle Scholar
- Tucker SL, Talbot NJ: Surface attachment and pre-penetration stage development by plant pathogenic fungi. Ann Rev Phytopathol. 2001, 39: 384-415.View ArticleGoogle Scholar
- Prince VE, Pickett FB: Splitting pairs: the diverging fates of duplicated genes. Nature Rev Genets. 2002, 11: 827-837.View ArticleGoogle Scholar
- Long M, Betran E, Thornton K, Wang W: The origin of new genes: glimpses from the young and old. Nature Rev Genets. 2003, 4: 865-875.View ArticleGoogle Scholar
- Ohno S: Evolution by Gene Duplication. 1970, Berlin: Springer-VerlagView ArticleGoogle Scholar
- Hughes AL: Adaptive evolution after gene duplication. Trends Genetics. 2002, 18: 433-434. 10.1016/S0168-9525(02)02755-5.View ArticleGoogle Scholar
- Lynch M, Conery JS: The evolutionary fate and consequences of duplicate genes. Science. 2000, 290: 1151-1155. 10.1126/science.290.5494.1151.PubMedView ArticleGoogle Scholar
- Pfeffer A, Friedman N, Regev A, WapinskiI: Natural history and evolutionary principles of gene duplication in fungi. Nature. 2007, 449: 54-61. 10.1038/nature06107.PubMedView ArticleGoogle Scholar
- Karlsson M, Stenlid J: Comparative evolutionary histories of the fungal chitinase gene family reveal non-random size expansions and contractions due to adaptive natural selection. Evol Bioinformatics. 2008, 4: 47-60.Google Scholar
- Daniel DS, Helen P, Christopher BB: Evolutionary analysis of glycosyl hydrolase family 28 (GH28) suggests lineage-specific expansions in necrotrophic fungal pathogens. Gene. 2011, 479 (1–2): 29-36.Google Scholar
- Kershaw MJ, Talbot NJ: Hydrophobins and repellents: proteins with fundamental roles in fungal morphogenesis. Fungal GenetBiol. 1998, 23: 18-33.View ArticleGoogle Scholar
- Jones DT, Taylor WR, Thornton JM: The rapid generation of mutation data matrices from protein sequences. Comp Appl Biosci. 1992, 8: 275-282.PubMedGoogle Scholar
- Olson Å, Aerts A, Asiegbu F, Belbahri L, Bouzid O, Broberg A, Canbäck B, Coutinho PM, Cullen D, Dalman K, Deflorio G, van Diepen LTA, Dunand C, Duplessis S, Durling M, Gonthier P, Grimwood J, Fossdal CG, Hansson D, et al: Insight into trade-off between wood decay and parasitism from the genome of a fungal forest pathogen. New Phytol. 2012, 194: 1001-1013. 10.1111/j.1469-8137.2012.04128.x.PubMedView ArticleGoogle Scholar
- Sunde M, Kwan AH, Templeton MD, Beever RE, Mackay JP: Structural analysis of hydrophobins. Micron. 2008, 39: 773-784. 10.1016/j.micron.2007.08.003.PubMedView ArticleGoogle Scholar
- Seidle-Seiboth V, Gruba S, Sezerman U, Schwecke T, Albayrak A, Neuhof T, Döhren HV, Baker SE, Kubicek CP: Novel hydrophobins from Trichoderma defines a new hydrophobin subclass: protein properties, evolution, regulation and processing. J Mol Evol. 2011, 72: 339-351. 10.1007/s00239-011-9438-3.View ArticleGoogle Scholar
- Britt GJ, Mikael RA, Mona HP, Jens CF, Ib S: Hydrophobins from Aspergillus species cannot be clearly divided into two classes. BMC Res Notes. 2010, 3: 344-10.1186/1756-0500-3-344.View ArticleGoogle Scholar
- Yamada M, Sakuraba S, Shibata K, Inatomi S, Okazaki M: Cloning and characterization of a gene coding for a hydrophobin, Fv-hyd1, specifically expressed during fruiting body development in the basidiomycete Flammulina velutipes. Appl Microbiol Biotechnol. 2005, 67: 240-246. 10.1007/s00253-004-1776-2.PubMedView ArticleGoogle Scholar
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215 (3): 403-410.PubMedView ArticleGoogle Scholar
- Edgar RC: MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32 (5): 1792-1797. 10.1093/nar/gkh340.PubMed CentralPubMedView ArticleGoogle Scholar
- Hall TA: Nucleic Acids Symposium. BioEdit: a user friendly biological sequence alignment editor and analysis program for Windows 95/98NT, Volume 41. 1990, Oxford: Oxford University PressGoogle Scholar
- Zdobnov EM, Apweiler R: InterProScan – an integration platform for the signature-recognition methods in InterPro. Bioinformatics. 2001, 17 (9): 847-848. 10.1093/bioinformatics/17.9.847.PubMedView ArticleGoogle Scholar
- Kyte J, Doolittle RF: A simple method for displaying the hydropathic character of a protein. J Mol Biol. 1982, 5 (1): 105-132.View ArticleGoogle Scholar
- Felsenstein J: Confidence limits on phylogenies : an approach using the bootstrap. Evol. 1985, 39: 783-791. 10.2307/2408678.View ArticleGoogle Scholar
- Chang S, Puryear J, Cairney J: A simple and efficient method for isolating RNA from pine trees. Plant Mol Biol Rep. 1993, 11: 113-116. 10.1007/BF02670468.View ArticleGoogle Scholar
- Hudson RR, Kreitman M, Aguade M: A test of neutral molecular evolution based on nucleotide data. Genetics. 1987, 116: 153-159.PubMed CentralPubMedGoogle Scholar
- Tajima F: Statistical methods to test for nucleotide mutation hypothesis by DNA polymorphism. Genets. 1989, 123: 585-595.Google Scholar
- Altschul SF, Madden TL, Schäffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25 (17): 3389-3402. 10.1093/nar/25.17.3389.PubMed CentralPubMedView ArticleGoogle Scholar
- Thompson JD, Higgins DG, Gibson TJ: CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994, 22 (22): 4673-4680. 10.1093/nar/22.22.4673.PubMed CentralPubMedView ArticleGoogle Scholar
- Berman HM, Westbrook J, Feng Z, Gilliland G, Bhat TN, Weissig H, Shindyalov IN, Bourne PE: The Protein Data Bank. Nucleic Acids Res. 2000, 28 (1): 235-242. 10.1093/nar/28.1.235.PubMed CentralPubMedView ArticleGoogle Scholar
- Jones TA, Zou JY, Cowan SW, Kjeldgaard M: Improved methods for building protein models in electron density maps and the location of errors in these models. Acta Crystallogr A. 1991, 47 (Pt 2): 110-119.PubMedView ArticleGoogle Scholar
- Kleywegt GJ, Zou JY, et al: Around O. International Tables for Crystallography. Crystallography of Biological Macromolecules. Edited by: Rossmann MG, Arnold E. 2001, Dordrecht: Kluwer Academic, 353-356.Google Scholar
- Mackay JP, Matthews JM, Winefield RD, Mackay LG, Haverkamp RG, Templeton MD: The hydrophobin EAS is largely unstructured in solution and functions by forming amyloid-like structures. Structure. 9 (2): 83-91.
- Tripathi A, Ren Y, Jeffrey PD, Hughson FM: Structural characterization of Tip20p and Dsl1p, subunits of the Dsl1p vesicle tethering complex. Nature Struct Mol Biol. 2009, 16 (2): 114-123. 10.1038/nsmb.1548.View ArticleGoogle Scholar
- Guido RVC, Balliano TL, Andricopulo AD, Oliva G: Kinetic and Crystallographic Studies on Glyceraldehyde-3-Phosphate Dehydrogenase from Trypanosoma cruzi in Complex with Iodoacetate. Lett Drug Design & Disc. 2009, 6 (3): 210-214. 10.2174/157018009787847774.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.