Purifying selection and birth-and-death evolution in the class II hydrophobin gene families of the ascomycete Trichoderma/Hypocrea
© Kubicek et al; licensee BioMed Central Ltd. 2008
Received: 15 August 2007
Accepted: 10 January 2008
Published: 10 January 2008
Hydrophobins are proteins containing eight conserved cysteine residues that occur uniquely in mycelial fungi. Their main function is to confer hydrophobicity to fungal surfaces in contact with air or during attachment of hyphae to hydrophobic surfaces of hosts, symbiotic partners or themselves resulting in morphogenetic signals. Based on their hydropathy patterns and solubility characteristics, hydrophobins are divided into two classes (I and II), the latter being found only in ascomycetes.
We have investigated the mechanisms driving the evolution of the class II hydrophobins in nine species of the mycoparasitic ascomycetous genus Trichoderma/Hypocrea, using three draft sequenced genomes (H. jecorina = T. reesei, H. atroviridis = T. atroviride; H. virens = T. virens) an additional 14,000 ESTs from six other Trichoderma spp. (T. asperellum, H. lixii = T. harzianum, T. aggressivum var. europeae, T. longibrachiatum, T. cf. viride). The former three contained six, ten and nine members, respectively. Ten is the highest number found in any ascomycete so far. All the hydrophobins we examined had the conserved four beta-strands/one helix structure, which is stabilized by four disulfide bonds. In addition, a small number of these hydrophobins (HFBs)contained an extended N-terminus rich in either proline and aspartate, or glycine-asparagine. Phylogenetic analysis reveals a mosaic of terminal clades containing duplicated genes and shows only three reasonably supported clades. Calculation of the ratio of differences in synonymous vs. non-synonymous nucleotide substitutions provides evidence for strong purifying selection (K S /K a >> 1). A genome database search for class II HFBs from other ascomycetes retrieved a much smaller number of hydrophobins (2–4) from each species, and most were from Sordariomycetes. A combined phylogeny of these sequences with those of Trichoderma showed that the Trichoderma HFBs mostly formed their own clades, whereas those of other Sordariomycetes occurred in shared clades.
Our study shows that the genus Trichoderma/Hypocrea has a proliferated arsenal of class II hydrophobins which arose by birth-and-death evolution followed by purifying selection.
Hydrophobins are small proteins that are unique for mycelial fungi [1–4]. Their core structure consists of four beta strands, crosslinked by four disulfide bridges [5, 6], creating a structure enabling the self-assembly at a hydrophilic/hydrophobic interface between the hydrophilic cell wall and a hydrophobic environment (such as air or the hydrophobic surface of living and non-living material). They are classified on the basis of chemical properties (hydrophobicity, solubility) into class I or class II hydrophobins , with class I hydrophobins having been identified in both ascomycetes and basidiomycetes, and class II hydrophobins having so far been only detected in ascomycetes . Hydrophobins can fulfill a plethora of biological functions ranging from the formation of aerial structures, elicitation of morphogenesis and interaction with potential hosts or symbiotic partners [2, 7–10].
In the primary sequence, the most important feature common to all hydrophobins is the characteristic pattern of eight Cys-residues, which gives rise to a common disulfide network [5, 11]. Besides these conserved Cys-residues and similar hydropathy patterns, however, the hydrophobins share only a few conserved residues . The poor amino acid sequence conservation of hydrophobins raises the question as to the evolutionary mechanism driving the rapid differentiation of hydrophobin gene sequences. Other genes involved in the response of organisms to their immediate environment have sometimes been shown to be driven by positive selection [12–16], so called "arms races" . However, concerted evolution and birth-and-death evolution under strong purifying selection have also been reported [18–24].
The fungal genus Trichoderma/Hypocrea contains a large number of mycoparasitic species [25, 26], and some of them (e.g. H. lixii = T. harzianum, H. virens = T. virens, H. atroviridis = T. atroviride, T. asperellum) are also commercially used as biological fungicides . Interestingly, there are only a few reports describing the occurrence of hydrophobins in mycoparasitic strains of Trichoderma/Hypocrea [28–30] and a characterization of the function or biological role for each has not been determined. The class II hydrophobin genes hfb1 and hfb2 of the weakly mycoparasitic species H. jecorina (T. reesei) have been studied [31, 32] and shown to serve different functions during vegetative development. Viterbo and Chet  recently showed that a class I hydrophobin from T. asperellum is involved in root colonization. Interestingly, this hydrophobin is the only class I hydrophobin identified in any Trichoderma/Hypocrea sp. so far.
The assignment of a function to individual members of large gene families like the hydrophobins is complicated by the possibility that several of them may have overlapping functions [34, 35], which in turn is dependent on the selective pressures acting on the organism. Understanding the evolution of such genes and identifying stable clusters within the phylogeny may therefore help illustrate members with a potentially critical function.
Protein structure of the Trichoderma class II hydrophobin proteins
Trichoderma/Hypocrea class II hydrophobin genes*
GenBank: EU053448] [GenBank: EU053449];
HA_2a = SRH1
T. cf. Koningiopsis
T. cf. Koningiopsis
T. cf. Koningiopsis
T. aggressivum var. europeae
T. aggressivum var. europeae
T. aggressivum var. europeae
Genomic organisation of the Trichoderma hydrophobin genes
In order to identify the mechanisms acting on the evolution of the Trichoderma hydrophobin genes, we first looked at their genomic organisation and exon structures. The six T. reesei hydrophobins are located on five different, large scaffolds, and even the two which are located on the same scaffold (hfb1 and hfb6) are separated by over 700,000 bp, and are therefore unlinked. Similar, although several of the hfb genes of H. atroviridis and H. virens were located on the same scaffold, they were separated by over 100,000 bp. In order to analyse whether any of these loci would be syntenic across Trichoderma species, we subjected each of the H. jecorina hfb genes plus 5 kb of its up- and downstream nt-sequences to a TBLAST search in the genome sequences of H. virens and H. atroviridis. The loci flanking the six hfb genes in H. jecorina did not flank any of the H. virens and H. atroviridis genes found by this analysis, although the flanking genes alone were sometimes located in the same region on a different scaffold (data not shown). Together, these data suggest that the genome regions containing hydrophobin gene loci have undergone extensive recombination during their evolution.
Intron/Exon structure of the Trichoderma hydrophobin genes
All of the chromosomal hfb genes of H. jecorina, H. virens and H. atroviridis contain two introns, which are very similar although not identical in size, and are positionally conserved. Interestingly, the length of the second exon encodes the aa sequence which folds exactly into the single α-helix of the hydrophobin and its third beta-sheet (cf. ) is absolutely conserved in all genes.
Phylogeny of the Trichoderma class II hydrophobin proteins
Evidence for gene duplications within the Trichoderma HFB proteins
One feature, which became obvious from the phylogenetic analysis and which is unaffected by the low internal branching support, is the high number of paralogous proteins. Examples for this are: Ha_1b and Ha_1c; Ha_6a and Ha_2c; Hv_21a and Hv_22a; TCK1 and TCK2; Ha_6b and Ha_6c; and HFB1 and HFB2. Most of these twins form a terminal branch, or are connected by a single node, indicating that they arose by gene duplication. The Trichoderma class II hydrophobins thus display a significant pattern of gene duplications in their evolutionary history.
Nucleotide Sequence Divergence of the hydrophobin genes
The lack of recombination within the hfb genes, together with the observation of gene duplications suggests that the Trichoderma hydrophobins undergo purifying selection without concerted evolution. In such a case, the member genes would evolve independently and display a birth-and-death evolution . This model of evolution assumes that new genes are created by repeated gene duplication and that some of the duplicate genes are maintained in the genome for a long time whereas others are deleted or become non-functional [18–20]. Thus nucleotide sequence differences between genes will primarily occur at synonymous sites, thus resulting in K S >> K a .
Nucleotide diversity of the Trichoderma/Hypocrea class II hydrophobin genes
all three exons
no of sites
no syn sites
no nonsyn sites
Eta, no of mutations
Pi nt diversity
0.38773 ± 0.0164
0.42532 ± 0.0139
Theta per site (from Eta)
The Trichoderma HFB proteins form unique clades within the ascomycetous hydrophobins
Class II hydrophobin genes from other ascomycetes used in this study
Accession number *
[B. fuckeliana genome database: BC1G_03994.19]
[B. fuckeliana genome database: BC1G_01012.1 ]
[M. graminicola genome database: FGENESH2_PG.C_SCAFFOLD_8000534]
[M. graminicola genome database: FGENESH2_PG.C_SCAFFOLD_2000556]
[M. graminicola genome database: FGENESH2_PG.C_SCAFFOLD_11000390]
[N. haematococca genome database: e_gw.18.104.22.168]
[GenBank: XM 364289]
Expression analysis of the "HFB4"-clade hydrophobin genes from H. atroviridis
In this paper, we have investigated the evolutionary processes which give rise to the biodiversity of the fungal class II hydrophobin genes and proteins. This class of hydrophobins has so far been reported to be restricted to Ascomycetes only. Results from this study, however, suggest that the distribution of these genes may be even more restricted, i.e. the majority of the members of this class was actually found in the Sordariomycetes, and only few were found in Leotiomycetes and Eurotiomycetes. This picture may however be biased by the fact that Sordariomycetes are overrepresented sequenced genomes, and the six genes which we retrieved from the two species of Dothidiomycetes (Passalora, Mycosphaerella) suggest that this subphylum may also be rich in class II hydrophobins. We cannot completely rule out that the low number retrieved for Leotio- and Eurotiomycetes could be due to a failure to identify these genes by BLAST search. However, our approach also identified several class I hydrophobin genes from all these fungi (data not shown), and we would therefore assume that our screening was broad enough to identify all class II genes. Also, our results are in agreement with the results of manual annotation of several fungal genomes (Aspergillus spp., M. grisea, N. crassa, G. zeae, N. haematococcae). Therefore, while it is possible that a potential HFB encoding gene has been overlooked, our data indicate that while most species contained only 1 or 2 genes (e.g. Gibberella, Nectria, Botryotinia, Aspergillus spp.), species of Trichoderma/Hypocrea clearly exceed this with their gene number (i.e. 6 genes in H. jecorina, 9 in H. virens and 10 in H. atroviridis). The reason for this remains obscure: neither the morphology of the hyphae, the conidia or of the perithecium of Hypocrea/Trichoderma show microscopic differences which may necessitate new or multiple hydrophobins to support these structures. What differentiates this fungal genus from others, however, is its mycoparasitic and necrotrophic lifestyle . While completely speculative at this moment, it is nevertheless possible that a versatile arsenal of class II hydrophobins may help the fungus to attach to the hyphae of a broad range of asco- and basidiomycetes. An amplified spectrum of genes has also been found for the chitinases of H. jecorina, which undoubtedly also aid to its mycoparasitic abilities . With the availability of the hydrophobin gene sequences now in hand for two strongly mycoparasitic species – H. atroviridis and H. virens – this work lays a strong phylogenetic foundation to investigate this possibility by means of respective knock-out strains and expression analysis.
The results from this paper show that the class II hydrophobin genes of Trichoderm/Hypocrea contain a high number of duplicated genes, and at least two cases of pseudogenes. This suggests that the class II hydrophobins evolve by a death-and-birth mechanism , a term which has been created for a process in which genes undergo gene duplications, resulting in the maintenance of some of the copies for a considerable period of time whereas other copies are rapidly lost or converted to pseudogenes. Our data render the operation of concerted evolution unlikely, because of the high sequence divergence and also by the absence of recombination at the hydrophobin loci. (in concerted evolution, member genes evolve together as a unit by mechanisms such as gene conversion or unequal crossing-over). The fact that most of the duplicated genes occupy terminal branches in phylogenetic trees and that only few obvious pseudogenes were found, indicates that the rate of evolution of the class II hydrophobins in Trichoderma is relatively fast. This rapid evolution, and equally rapid loss of some genes is also reflected in the findings that the clades leading to the hfb genes in Trichoderma seldom contain members of other fungi, and their evolution thus took place after formation of the genus Trichoderma.
In addition, the numbers of synonymous differences of nucleotide sequences between genes from the same species are very large and frequently close to the saturation level. This high level of synonymous differences further supports the claim of a birth-and-death evolution at the DNA level, and supports the long time persistence of these genes in the genome. On the other hand, genes from different species (e.g. H. atroviridis and H. virens) but belonging to the same phylogenetic clade are highly similar (cf. Fig. 3). Such a long-term conservation of amino acid sequence is best explained by strong purifying selection. Interestingly, and in contrast to Rajashekar et al. , we found only a few individual cases where the K a /Ks ratio was >1 and would reveal a history of accelerated evolution. If such a period of accelerated evolution occurred, most of the gene duplicates from this time apparently have not been maintained and the Trichoderma/Hypocrea hydrophobin genes characterized in this study are therefore mostly of recent origin.
The present study also expands our knowledge on the structure of class II hydrophobins. While most of them are small, compact proteins, which contain little other structures than the four beta-sheets and the single helix [5, 6], we have detected several proteins which display a long extended N-terminus (ENT). With respect to class II HFBs, such structures have so far only been found in H. jecorina HFB6 , and in the pseudohydrophobin QID3 . Interestingly, an ENT was recently also identified in the class I HFB Hum3 from Ustilago mayidis . Lora et al.  hypothesized that the ENT of QID3 mediates cell wall binding because it resembles a module which is also present in plant bimolecular proteins [47–49]. Interestingly, our work reveals that there are at least two types of these ENTs: a major one, typified by H. virens HV_21a, H. atroviridis HA_2c, H. lixii QID3, and also in P. fulva HCF6, and in the spacers between the hydrophobin units in the multipartite genes of C. paspali and C. fusii, which are characterized by a conserved repeat of glycine and asparagines; and second type, shown by e.g. H. atroviridis HA_6a, H. virens HV_1d, and H. jecorina HFB6, in which the repeated motif is replaced by several PG/PD repeats, a P-rich stretch or a D-rich stretch, respectively. These proteins did not cluster together, indicating that these proteins do not show a common ancestry of the cysteine-containing core domain. Among the proteins with this terminus, one (HV_13a) is intriguing as its ENT is very short, which gives rise to the speculation that this extended N-terminus may arise by segment duplication. Support for this hypothesis would also come from the multipartite hydrophobins found in Claviceps spp. [50, 51], wherein paralogous hydrophobin gene copies are connected by P, G and N-rich loops, and which may have been trapped in the stage of gene duplications at the extended N-terminus before recombining individual copies into new loci. It is thereby intriguing to observe the similarity of the nucleotide sequence of the "GN" repeat (GGTAAT) to that found to act as a recombination hot-spot in Penicillium chrysogenum (TGTAA [A/T]; ). Therefore, the occurrence of the Claviceps multipartite hydrophobins would be due to multiplication of some of the class II hydrophobins by tandem duplication [53, 54], for which these sequences could act as recombination targets.
Nevertheless, it may still be likely that these ENTs are not only evolutionary artefacts: the [GN] repeats are reminiscent of S. cerevisiae Ure2p, a regulator of nitrogen catabolism, which can become transformed into a prion form by polymerization into filaments . These filaments have an amyloid fibril backbone formed by an N-rich sequence which form a parallel superpleated beta-structure. The prion domain is thereby divided into nine seven-residue segments, each with a four-residue strand and a three-residue turn, that zig-zag in a planar serpentine arrangement, the interior of the filament being stabilized by H-bond networks generated by the stacking of N side chains. Interestingly, hydrophobins themselves are known to form amyloid-like structures [56–58], and we consider it therefore possible that the ENTs form defined structures which additionally contribute to the structural rigidity of the hydrophobins.
During phylogenetic analysis of the amino acid sequence, most hydrophobins from Trichoderma/Hypocrea did not group into strongly supported clades. However, a few exceptions were noted, notably the clades containing H. jecorina HFB1, HFB2 and HFB4, respectively. Clade "HFB4" is intriguing as its members – in contrast to HFB1 and HFB2  were not expressed in submerged culture but only found in surface cultivation. This clade may thus contain hydrophobins relevant for hyphal growth. Unfortunately, the differences in aa-sequence with that of the other Trichoderma hydrophobins do not provide a clear clue as to the understanding of its function. One notable change is the substitution of the phenylalanine residue in the middle of the single helix, which is otherwise conserved in the class II hydrophobins from all other fungi (with the exception of A. niger and V. lecanii which contain an L and M, respectively), by a leucine, which may give rise to weaker hydrophobic interaction within the protein (hydrophobicity index F = 2.24; L = 1.99). However, Linder et al.  speculated that the aromatic ring of F39 is inserted between two Pro-residues (P11 and P50) from the two β hairpin structures into the protein and may serve to stabilize the fold through hydrogen bonds. Interestingly, members of the "HFB 4 clade" consistently have P11 replaced by an A which may hydrophobically interact with this L. In addition, the first beta-strand contains a conserved motif of two asparagines which provide it with a positive charge. Hydropathy plots show that members of the "HFB4 clade" have almost no hydrophobicity in the area between aa20 and aa40, and their helix is in contrast positively charged. While the consequence of these changes on the structure and function is however unclear our phylogenetic analysis sets the stage for future functional studies that may include transcript and gene deletion analysis.
Summarizing, this study offers a model of evolution which gave rise to the biodiversity of class II hydrophobins. The more than 70 members identified in this study enabled us to delineate the consensus for both essential aa sequence parts as well as for the tolerance to aa modifications to these small compact proteins. Our phylogenetic analysis will inform future functional genomic studies aimed at determination of more specific functions for each of the Trichoderma/Hypocrea class II hydrophobins. In view of the strong potential of the hydrophobins in "white" biotechnology, this information may offer their further improvement by molecular evolution [59–61].
Conditions of fungal growth
Hypocrea atroviridis (anamorph: Trichoderma atroviride) strain P1 (ATCC 74058) was maintained on potato dextrose agar (PDA) (Difco, Franklin Lakes, NJ, USA). The following medium was used for its cultivation (g·l-1; ): D-glucose, 10; peptone, 0.35; Tween 80, 0.175; KH2PO4, 0.68; K2HPO4, 0.87; (NH4)2SO4, 1.7; KCl, 0.2; CaCl2, 0.2; MgSO4·7H2O, 0.2; FeSO4·7H2O, 0.02; ZnSO4·7H2O, 0.02; MnSO4·7H2O, 0.02. The fungus was grown either in 1 L shake-flasks containing 200 ml of the medium at 25°C (250 RPM) or on plates (in this case the medium was solidified by the addition of 15 g·l-1 agar).
Analysis of hydrophobin gene expression
Mycelia were withdrawn at selected time points as indicated and total RNA isolated by the method of Chomzynski and Sacchi . The RNA extract was treated with DNAse I (Fermentas, St. Leon-Rot, Germany) and purified using the RNeasy MinElute Cleanup Kit (Qiagen, Hilden, Germany). The purified RNA was reverse transcribed using the RevertAid H Minus First Strand cDNA Synthesis Kit (Fermentas) with the oligo(dT)18 primer supplied by the manufacturer.
Primers used for RT-PCR
5' -> 3' sequence
In silico screening for hydrophobin sequences from Trichoderma and other fungi
To obtain the Trichoderma class II hydrophobins, we used the sequences of the six class II hydrophobins of H. jecorina  as a tool to retrieve genes encoding proteins with similarity from the genome sequence databases of Hypocrea virens and Hypocrea atroviridis and the TrichoEST EST database which includes EST sequences from 9 different Trichoderma spp.  using TBLAST (protein vs translated nucleotide). In the case of duplicates, the genomic sequence was given preference. Table 1 summarizes the genes, proteins, accession numbers or locations of the sequences thereby retrieved.
To obtain class II hydrophobin genes from other ascomycetes, the genes compiled by Linder et al.  were used as a starting point for a TBLAST search of the NCBI data base, using the filtering option turned off, and sequences which had not yet been included by Linder et al.  were retrieved. Apart of genes deposited from specific research, this database contains genome sequences from Neurospora crassa, Magnaporthe grisea, Gibberella zeae, and several Aspergillus spp. In addition, the same procedure was used to mine the genome databases of Nectria cinnabarina , Mycophaerella graminicola  and Botryotinia fuckeliana (= Botrytis cinerea, ), fungi whose sequences have not yet been included in the NCBI database. Genes thereby retrieved were then themselves used as a query in BLAST search as described above, and the procedure repeated until no new gene/protein was detected. The amino acid sequences of the retrieved proteins were aligned, and class I hydrophobins (identified according to the criteria described by Linder et al. ; e.g. by the difference in the number of amino acids between the conserved cysteins and their hydropathy profile) removed.
DNA and protein sequences were visually aligned using Genedoc 2.6 . Phylogenetic trees were constructed by the neighbor-joining (NJ) method , using the computer program MEGA, Version 4.0 , and by Bayesian analysis (MrBayes v3.0B4 program). The model of evolution and prior settings for individual loci was GTR + I + Γ. Metropolis-coupled Markov chain Monte Carlo (MCMCMC) sampling was performed with four incrementally heated chains that were simultaneously run for 1 and 3 millions of generations. To check for potentially poor mixing of MCMCMC, each analysis was repeated three times. The convergence of MCMCMC was monitored by examining the value of the marginal likelihood through generations. Convergence of substitution rate and rate heterogeneity model parameters was also checked. Bayesian posterior probabilities (PP) were obtained from the 50% majority rule consensus of trees sampled every 100 generations after removing the 500 first trees using the "burn" command. PP values lower then 0.95 were not considered significant.
Test for recombination
Two different procedures were applied to detect recombination by comparing adjacent sequence windows and to detect significant departures from a single phylogenetic history within the same alignment. First, we used difference of sums of squares (DSS, ) to compare the fit of genetic distance matrices for two adjacent windows to the same tree topology to produce the DSS statistic, the significance of which was determined through parametric bootstrapping in TOPALi .
In addition we tested whether recombination would be apparent from the phylogeny of hydrophobin genes. This was done by the split decomposition method in SplitsTree, version 2.4 , using pairwise distances under the Kimura 3-ST model . This method visualizes recombination events by depicting the shortest pathway linking sequences, rather than forcing them into a bifurcating tree [39, 40].
Tests for evolutionary mechanisms
To test the fit of the sequences to the model of neutral evolution, the D test statistic proposed by Tajima and Nei  was computed with the DnaSP program . To this end, the Genedoc alignment was exported as a PHYLIP interleaved format. Only coding sequences, after removal of the preprosequence-encoding nt areas, were used for the analyses. Introns were removed from chromosomal nt-sequences by comparison with available cDNA (EST) sequences, or by relying on prediction of consensus splicing sites [75, 76]. In the case of the multipartite hydrophobins from Claviceps spp. [50, 51], each hydrophobin-encoding nt-area was treated as a separate entity. The extent of nucleotide divergence was estimated by using the uncorrected p distance . The proportions of synonymous (were calculated pS) and nonsynonymous (pN) differences per site by the modified Nei-Gojobori method implemented in DNASp .
The authors are grateful to Heather Wilkinson for discussions on the manuscript. This work was supported by the EC-funded project TrichoEST (QLK3-2002-02032). CPK wants to acknowledge the important contribution of their colleagues from the TrichoEST consortium to the generation of the EST database. The genomic sequences for H. atroviridis and H. virens were provided by the DOE Joint Genome Institute, Walnut Creek, CA. This work was performed under the auspices of the US Department of Energy's Office of Science, Biological and Environmental Research Program and the by the University of California, Lawrence Livermore National Laboratory under Contract No. W-7405-Eng-48, Lawrence Berkeley National Laboratory under contract No. DE-AC03-76SF00098 and Los Alamos National Laboratory under contract No. W-7405-ENG-36.
- Wessels JG: Fungi in their own right. Fungal Genet Biol. 1999, 27: 134-145. 10.1006/fgbi.1999.1125.View ArticlePubMed
- Talbot NJ: Fungal biology. Coming up for air and sporulation. Nature. 1999, 398: 295-296. 10.1038/18575.View ArticlePubMed
- Linder MB, Szilvay GR, Nakari-Setala T, Penttila ME: Hydrophobins: the protein-amphiphiles of filamentous fungi. FEMS Microbiol Rev. 2005, 29: 877-896. 10.1016/j.femsre.2005.01.004.View ArticlePubMed
- Wösten HA: Hydrophobins: multipurpose proteins. Annu Rev Microbiol. 2001, 55: 625-646. 10.1146/annurev.micro.55.1.625.View ArticlePubMed
- Hakanpaa J, Paananen A, Askolin S, Nakari-Setala T, Parkkinen T, Penttila ME, Linder MB, Rouvinen J: Atomic resolution structure of the HFBII hydrophobin, a self-assembling amphiphile. J Biol Chem. 2004, 279: 534-539. 10.1074/jbc.M309650200.View ArticlePubMed
- Hakanpaa 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. Protein Sci. 2006, 15: 2129-2140. 10.1110/ps.062326706.PubMed CentralView ArticlePubMed
- 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.View ArticlePubMed
- Tucker S, Talbot NJ: Surface attachment and pre-penetration stage development by plant pathogenic fungi. Annu Rev Phytopathol. 2001, 39: 385-417. 10.1146/annurev.phyto.39.1.385.View ArticlePubMed
- Klimes A, Dobinson KF: A hydrophobin gene, VDH1, is involved in microsclerotial development and spore viability in the plant pathogen Verticillium dahliae. Fungal Genet Biol. 2006, 43: 283-294. 10.1016/j.fgb.2005.12.006.View ArticlePubMed
- Sakamoto Y, Ando A, Tamai Y, Yajima T: Pileus differentiation and pileus-specific protein expression in Flammulina velutipes. Fungal Genet Biol. 2007, 44: 14-24. 10.1016/j.fgb.2006.06.002.View ArticlePubMed
- Kwan AH, Winefield RD, Sunde M, Matthews JM, Haverkamp RG, Templeton MD, Mackay JP: Structural basis for rodlet assembly in fungal hydrophobins. Proc Natl Acad Sci USA. 2004, 103: 3621-3626. 10.1073/pnas.0505704103.View Article
- Zhang J, Rosenberg HF, Nei M: Positive Darwinian selection after gene duplication in primate ribonuclease genes. Proc Natl Acad Sci USA. 1998, 95: 3708-3713. 10.1073/pnas.95.7.3708.PubMed CentralView ArticlePubMed
- Weber E, Köbnik R: Positive selection of the Hrp pilin HrpE of the plant pathogen Xanthomonas. J Bacteriol. 2006, 188: 1405-1410. 10.1128/JB.188.4.1405-1410.2006.PubMed CentralView ArticlePubMed
- Benderoth M, Textor AJ, Windsor T, Mitchell-Olds J, Gershenzon J, Kroymann J: Positive selection driving diversification in plant secondary metabolism. Annu Rev Phytopathol. 2006, 44: 469-487. 10.1146/annurev.phyto.44.032905.091147.View Article
- Brown NF, Wickham ME, Coombes BK, Finlay BB: Crossing the line: selection and evolution of virulence traits. PLoS Patho. 2006, 2: e42-10.1371/journal.ppat.0020042. Review.View Article
- De Mita S, Santoni S, Hochu I, Ronfort J, Bataillon T: Molecular evolution and positive selection of the symbiotic gene NORK in Medicago truncatula. J Mol Evol. 2006, 62: 234-244. 10.1007/s00239-004-0367-2.View ArticlePubMed
- Haraguchi Y, Sasaki A: Host-parasite arms race in mutation modifications: indefinite escalation despite a heavy load?. J Theor Biol. 1996, 183: 121-137. 10.1006/jtbi.1996.9999.View ArticlePubMed
- Nei M, Hughes AL: Balanced polymorphism and evolution by the birth-and-death process in the MHC loci. HLA. Edited by: Tsuji K, Aizawa M, Sasazuli T. 1991, Oxford University Press, 2: 27-38. . Proc 11th Histocompatibility Workshop and Conference
- Ota T, Nei M: Divergent evolution and evolution by the birth-and-death process in the immunoglobulin VH gene family. Mol Biol Evol. 1994, 15: 469-482.
- Nei M, Gu X, Sitnikova T: Evolution by a birth-and-death process in multigene families of the vertebrate immune system. Proc Natl Acad Sci USA. 1997, 94: 7799-7806. 10.1073/pnas.94.15.7799.PubMed CentralView ArticlePubMed
- Nei M, Rogozin IB, Piontkivska H: Purifying selection and birth-and-death evolution in the ubiquitin gene family. Proc Natl Acad Sci USA. 2000, 97: 10866-10871. 10.1073/pnas.97.20.10866.PubMed CentralView ArticlePubMed
- Nei M, Rooney AP: Concerted and birth-and-death evolution of multigene families. Annu Rev Genet. 2005, 39: 121-152. 10.1146/annurev.genet.39.073003.112240.PubMed CentralView ArticlePubMed
- Jiang RH, Tyler BM, Whisson SC, Hardham AR, Govers F: Ancient origin of elicitin gene clusters in Phytophthora genomes. Mol Biol Evol. 2006, 23: 338-351. 10.1093/molbev/msj039.View ArticlePubMed
- Loria R, Kers J, Joshi M: Evolution of plant pathogenicity in Streptomyces. Proc Natl Acad Sci USA. 2006, 103: 9118-9123. 10.1073/pnas.0601738103.View Article
- Klein D, Eveleigh DE: Ecology of Trichoderma. Trichoderma and Gliocladium. Basic biology, taxonomy and genetics. Edited by: Kubicek CP, Harman GE. 1998, Taylor & Francis Ltd, London, U, 1: 57-74.
- Harman GE, Howell CR, Viterbo A, Chet I, Lorito M: Trichoderma species – opportunistic, avirulent plant symbionts. Nature Rev Microbiol. 2004, 2: 43-56. 10.1038/nrmicro797.View Article
- Hjeljord L, Tronsmo A: Trichoderma and Gliocladium in biological control: an overview. Trichoderma and Gliocladium. Enzymes, Biological Control and Commercial Applications. Edited by: Harman GE, Kubicek CP. 1998, Taylor and Francis, London, 2: 131-151.
- Munoz G, Nakari-Setala T, Agosin E, Penttila ME: Hydrophobin gene srh1, expressed during sporulation of the biocontrol agent Trichoderma harzianum. Curr Genet. 1997, 32: 225-230. 10.1007/s002940050270.View ArticlePubMed
- Vizcaino JA, Gonzalez JF, Suarez MB, Redondo J, Heinrich J, Delgado-Jarana J, Hermosa R, Gutierrez S, Monte E, Llobell A, Rey M: Generation, annotation and analysis of ESTs from Trichoderma harzianum CECT 2413. BMC Genomics. 2006, 27;7: 193-10.1186/1471-2164-7-193.View Article
- Vizcaino JA, Redondo J, Suarez MB, Cardoza RE, Hermosa R, Gonzalez FJ, Rey M, Monte E: Generation, annotation, and analysis of ESTs from four different Trichoderma strains grown under conditions related to biocontrol. Appl Microbiol Biotechnol. 2007, 75: 853-862. 10.1007/s00253-007-0885-0.View ArticlePubMed
- Askolin S, Penttila ME, Wösten HA, Nakari-Setala T: The Trichoderma reesei hydrophobin genes hfb1 and hfb2 have diverse functions in fungal development. FEMS Microbiol Letts. 2005, 253: 281-288. 10.1016/j.femsle.2005.09.047.View Article
- Askolin S, Linder M, Scholtmeijer K, Tenkanen M, Penttila M, de Vocht ML, Wosten HA: Interaction and comparison of a class I hydrophobin from Schizophyllum commune and class II hydrophobins from Trichoderma reesei. Biomacromolecules. 2006, 7: 1295-1301. 10.1021/bm050676s.View ArticlePubMed
- Viterbo A, Chet I: TasHyd1, anew hydrophobin gene from the biocontrol agent Trichoderma asperellum, is involved in plant root colonization. Mol Plant Pathol. 2006, 7: 249-258. 10.1111/j.1364-3703.2006.00335.x.View ArticlePubMed
- van Wetter MA, Wosten HA, Wessels JG: SC3 and SC4 hydrophobins have distinct roles in formation of aerial structures in dikaryons of Schizophyllum commune. Mol Microbiol. 2000, 36: 201-210. 10.1046/j.1365-2958.2000.01848.x.View ArticlePubMed
- Kershaw MJ, Wakley G, Talbot NJ: Complementation of the mpg1 mutant phenotype in Magnaporthe grisea reveals functional relationships between fungal hydrophobins. EMBO J. 1998, 17: 3838-3849. 10.1093/emboj/17.14.3838.PubMed CentralView ArticlePubMed
- Lora JM, de la Cruz J, Benitez T, Llobell A, Pintor-Toro JH: A putative carbon catabolite repressed cell wall protein from the mycoparasitic fungus Trichoderma harzianum. Mol Gen Genet. 1994, 242: 461-466. 10.1007/BF00281797.View ArticlePubMed
- Kullnig C, Krupica T, Woo SL, Mach RL, Rey M, Benitez T, Lorito M, Kubicek CP: Confusion abounds over identities of Trichoderma biocontrol isolates. Mycol Res. 2001, 105: 770-772. 10.1017/S0953756201229967.View Article
- The TrichoEST database. [http://www.trichoderma.org]
- Bandelt HJ, Dress AW: Split decomposition: a new and useful approach to phylogenetic analysis of distance data. Mol Phylogenet Evol. 1992, 1: 242-252. 10.1016/1055-7903(92)90021-8.View ArticlePubMed
- Huson DH: SplitsTree: a program for analyzing and visualizing evolutionary data. Bioinformatics. 1998, 14: 68-73. 10.1093/bioinformatics/14.1.68.View ArticlePubMed
- Nei M, Gojobori T: Simple methods for estimating the numbers of synonymous and nonsynonymous nucleotide substitutions. Mol Biol Evol. 1986, 3: 418-426.PubMed
- Seidl V, Huemer B, Seiboth B, Kubicek CP: A complete survey of Trichoderma chitinases reveals a new family 18 subgroup. The FEBS J. 2005, 272: 5923-5939. 10.1111/j.1742-4658.2005.04994.x.View ArticlePubMed
- Rajashekar B, Samson P, Johansson T, Tunlid A: Evolution of nucleotide sequences and expression patterns of hydrophobin genes in the ectomycorrhizal fungus Paxillus involutus. New Phytol. 2007, 174: 399-411. 10.1111/j.1469-8137.2007.02022.x.View ArticlePubMed
- Neuhof T, Dieckmann R, Druzhinina IS, Kubicek CP, Nakari-Setala T, Penttilä M, von Döhren H: Direct identification of hydrophobins and their processing in Trichoderma using intact-cell MALDI-TOF MS. FEBS J. 2007, 274: 841-852. 10.1111/j.1742-4658.2007.05636.x.View ArticlePubMed
- Teertstra WR, Deelstra HJ, Vranes M, Bohlmann R, Kahmann R, Kämper J, Wösten HA: Repellents have functionally replaced hydrophobins in mediating attachment to a hydrophobic surface and in formation of hydrophobic aerial hyphae in Ustilago maydis. Microbiology. 2006, 152: 3607-12. 10.1099/mic.0.29034-0.View ArticlePubMed
- Lora JM, Pintor-Toro JA, Benitez T, Romero LC: Qid3 protein links plant bimodular proteins with fungal hydrophobins. Mol Microbiol. 1995, 18: 380-382. 10.1111/j.1365-2958.1995.mmi_18020380.x.View ArticlePubMed
- Coupe SA, Taylor JE, Isaac PG, Roberts JA: Identification and characterization of a proline-rich mRNA that accumulates during pod development in oilseed rape (Brassica napus L.). Plant Mol Biol. 1993, 23: 1223-1232. 10.1007/BF00042355.View ArticlePubMed
- Castonguay Y, Laberge S, Nadeau P, Vezina LP: A cold-induced gene from Medicago sativa encodes a bimodular protein similar to developmentally regulated proteins. Plant Mol Biol. 1994, 24: 799-804. 10.1007/BF00029861.View ArticlePubMed
- He CY, Zhang JS, Chen SY: A soybean gene encoding a proline-rich protein is regulated by salicylic acid, an endogenous circadian rhythm and by various stresses. Theor Appl Genet. 2002, 104: 1125-1131. 10.1007/s00122-001-0853-5.View ArticlePubMed
- Arntz C, Tudzynski P: Identification of genes induced in alkaloid-producing cultures of Claviceps sp. Curr Genet. 1997, 31: 357-360. 10.1007/s002940050216.View ArticlePubMed
- De Vries OM, Moore S, Arntz C, Wessels JG, Tudzynski P: Identification and characterization of a tri-partite hydrophobin from Claviceps fusiformis. A novel type of class II hydrophobin. Eur J Biochem. 1999, 262: 377-385. 10.1046/j.1432-1327.1999.00387.x.View ArticlePubMed
- Fierro F, Garcia-Estrada C, Castillo NI, Rodriguez R, Velasco-Conde T, Martin JF: Transcriptional and bioinformatic analysis of the 56.8 kb DNA region amplified in tandem repeats containing the penicillin gene cluster in Penicillium chrysogenum. Fungal Genet Biol. 2006, 43: 618-29. 10.1016/j.fgb.2006.03.001.View ArticlePubMed
- Kong H, Landherr LL, Fröhlich MW, Leebens-Mack J, Ma H, Depamphilis CW: Patterns of gene duplication in the plant SKP1 gene family in angiosperms: evidence for multiple mechanisms of rapid gene birth. Plant J. 2007, Apr 23,
- Lynch M, Conery JS: The evolutionary fate and consequences of duplicate genes. Science. 2001, 290: 1151-1155. 10.1126/science.290.5494.1151.View Article
- Kajava AV, Baxa U, Wickner RB, Steven AC: A model for Ure2p prion filaments and other amyloids: the parallel superpleated beta-structure. Proc Natl Acad Sci USA. 2004, 101: 7885-7890. 10.1073/pnas.0402427101.PubMed CentralView ArticlePubMed
- Wosten HA, de Vocht ML: Hydrophobins, the fungal coat unravelled. Biochim Biophys Acta. 2000, 1469: 79-86.View ArticlePubMed
- Butko P, Buford JP, Goodwin JS, Stroud PA, McCormick CL, Cannon GC: Spectroscopic evidence for amyloid-like interfacial self-assembly of hydrophobin Sc3. Biochem Biophys Res Comm. 2001, 280: 212-215. 10.1006/bbrc.2000.4098.View ArticlePubMed
- 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. 2001, 9: 83-91. 10.1016/S0969-2126(00)00559-1.View ArticlePubMed
- Hektor HJ, Scholtmeijer K: Hydrophobins: proteins with potential. Curr Opin Biotechnol. 2005, 16: 434-439. 10.1016/j.copbio.2005.05.004.View ArticlePubMed
- Scholtmeijer K, Janssen MI, Gerssen B, de Vocht ML, van Leeuwen BM, van Kooten TG, Wösten HA, Wessels JG: Surface modifications created by using engineered hydrophobins. Appl Environ Microbiol. 2002, 68: 1367-1373. 10.1128/AEM.68.3.1367-1373.2002.PubMed CentralView ArticlePubMed
- Scholtmeijer K, Janssen MI, van Leeuwen MB, van Kooten TG, Hektor H, Wösten HA: The use of hydrophobins to functionalize surfaces. Biomed Mater Engin. 2004, 14: 447-454.
- Chomczynski P, Sacchi N: Single-step method of RNA isolation by acid guanidinium thiocyanate-phenol-chloroform extraction. Anal Biochem. 1987, 162: 156-159. 10.1016/0003-2697(87)90021-2.View ArticlePubMed
- Nakari T, Alatalo E, Penttila ME: Isolation of Trichoderma reesei genes highly expressed on glucose-containing media: characterization of the tef1 gene encoding translation elongation factor 1 alpha. Gene. 1993, 136: 313-318. 10.1016/0378-1119(93)90486-M.View ArticlePubMed
- DOE Joint Genome Institute. The Nectria cinnabarina genome database v 1.0. 2006, [http://genome.jgi-psf.org/Necha2/Necha2.home.html]
- DOE Joint genome Institute. The Mycophaerella graminicola genome database v 1.0. 2007, [http://genome.jgi-psf.org/Mycgr1/Mycgr1.home.html]
- Broad Institute and Syngenta AG. The Botrytis cinerea assembly. 2006, [http://www.broad.mit.edu/annotation/genome/botrytis_cinerea/Home.html]
- Nicholas KB, Nicholas HB: Genedoc: a tool for editing and annotating multiple sequence alignments. 1997, [http://www.psc.edu/biomed/genedoc]
- Saitou N, Saitou M: The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987, 4: 406-425.PubMed
- Kumar S, Tamura K, Nei M: MEGA3: Integrated software for Molecular Evolutionary Genetics Analysis and sequence alignment. Brief Bioinform. 2004, 5: 150-63. 10.1093/bib/5.2.150.View ArticlePubMed
- McGuire G, Wright F: TOPAL 2.0: improved detection of mosaic sequences within multiple alignments. Bioinformatics. 2000, 16: 130-134. 10.1093/bioinformatics/16.2.130.View ArticlePubMed
- Milne I, Wright F, Rowe G, Marshall DF, Husmeier D, McGuire G: TOPALi: software for automatic identification of recombinant sequences within DNA multiple alignments. Bioinformatics. 2004, 20: 1806-1807. 10.1093/bioinformatics/bth155.View ArticlePubMed
- Kimura M: Estimation of evolutionary distances between homologous nucleotide distances. Proc Natl Acad Sci USA. 1981, 78: 454-458. 10.1073/pnas.78.1.454.PubMed CentralView ArticlePubMed
- Tajima F, Nei M: Estimation of evolutionary distance between nucleotide sequences. Mol Biol Evol. 1984, 1: 269-285.PubMed
- Rozas J, Sanchez-DelBarrio JC, Messeguer X, Rozas R: DnaSP, DNA polymorphism analyses by the coalescent and other methods. Bioinformatics. 2003, 19: 2496-2497. 10.1093/bioinformatics/btg359.View ArticlePubMed
- Kupfer DM, Drabenstot SD, Buchanan KL, Lai H, Zhu H, Dyer DW, Roe BA, Murphy JW: Introns and splicing elements of five diverse fungi. Eukaryot Cell. 2004, 3: 1088-1100. 10.1128/EC.3.5.1088-1100.2004.PubMed CentralView ArticlePubMed
- Rep M, Duyvesteijn RG, Gale L, Usgaard T, Cornelissen BJ, Ma LJ, Ward TJ: The presence of GC-AG introns in Neurospora crassa and other euascomycetes determined from analyses of complete genomes: implications for automated gene prediction. Genomics. 2006, 87: 338-347. 10.1016/j.ygeno.2005.11.014.View ArticlePubMed
- Nei M, Kumar S: Molecular evolution and phylogenetics. 2000, Oxford University Press, New York
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.