Insights on the Evolution of Mycoparasitism from the Genome of Clonostachys rosea

1Agricultural chemicals for pest control and fertilization have played an important role for the substantial increase in the world’s food production during the past 50 years. The global human population is estimated to peak at around 9 billion people in the middle of this century, which will require further increases in food production. This must be achieved by producing more food from the same area of available land while reducing environmental impact. Integrated pest management (IPM) is one approach to reduce the amount of pesticides used, through preventive cultural practices, the use of disease resistant plant cultivars, and the use of mechanical and biological control of pathogen populations. Biological control of plant pathogens by microbial antagonists is one promising component in future disease control strategies that is compatible with both organic agriculture and IPM.

2One group of biological control agents (BCAs) includes mycotrophic fungi such as Trichoderma spp. and Clonostachys spp., which have the ability to parasitize and kill other fungi (mycoparasitism) and utilize dead fungal biomass (saprotrophy). These fungi can antagonize plant pathogenic fungi directly by secretion of cell wall degrading enzymes such as chitinases, β-1,3-glucanases, β-1,6-glucanases and proteases, and antibiotics such as peptaibols, gliotoxin, viridin and 6-pentyl-2H-pyran-2-one, but also through competition for space and nutrients. Certain Trichoderma and Clonostachys species can also promote plant growth and elicit induced resistance that can protect plants from pathogen attack. The exact mechanisms behind the disease control exerted by these BCAs may vary depending on the agricultural setting, including plant species, pathogen species, and environmental conditions. Therefore, a better understanding of the biological mechanisms that determine the outcome of biocontrol interactions is crucial for improving biological control in agricultural production systems.

3Clonostachys rosea (Link: Fr.) Schroers, Samuels, Seifert & W. Gams, comb. nov. is the anamorph stage, and preferred name, of the teleomorph Bionectria ochroleuca (Schw.) Schroers & Samuels. Clonostachys rosea belongs to the order Hypocreales and the family Bionectriaceae, but its taxonomic position in relation with other families within Hypocreales is debated. reported Bionectriaceae to be sister taxa with the family Nectriaceae, with Hypocreaceae in a basal position, whereas reported Bionectriaceae as basal to both Nectriaceae and Hypocreaceae. The more studied mycotrophic Trichoderma spp. also belong to the order Hypocreales, but to the family Hypocreaceae. Hence, an understanding the phylogenetic relationships of Bionectriaceae, Nectriaceae, and Hypocreaceae within Hypocreales is a key for understanding evolution of mycotrophism as a fungal lifestyle. Comparative genomics of mycotrophic species from different families enables us to identify key similarities and differences between their respective life strategies, which has direct implications for the implementation of biocontrol in agriculture.

4Recently, the sequenced and annotated genomes of the mycoparasitic species Trichoderma atroviride, Trichoderma hamatum, Trichoderma longibrachiatum and Trichoderma virens, and the saprotrophic Trichoderma reesei, were published. Comparative genomics revealed that the mycoparasitic lifestyle in Trichoderma was associated with gene copy number expansions of gene families involved in fungal cell wall degradation and secondary metabolite biosynthesis.

5Here, we report the genome sequencing and analysis of the first species from the family Bionectriaceae; the mycoparasite C. rosea. A comparative analysis of gene family evolution is performed under the hypothesis that evolution of mycoparasitism in Bionectriaceae and Hypocreaceae results in selection for similar interaction mechanisms. However, the majority of expanded gene families in C. rosea does not evolve in the same manner in Trichoderma spp. Clonostachys rosea contain high numbers of polyketide synthases (PKSs) and ATP-binding cassette (ABC) transporters predicted to be involved in drug resistance, which emphasizes the role of secondary metabolites in C. rosea biology. We further show that several ABC transporter genes are induced by xenobiotic substances and illustrate the agro-industrial potential for C. rosea to be applied together with low dose fungicide treatments or other BCAs with complementing mode of action, to achieve additive disease control effects.

6A complete description of Materials and Methods is found in supplementary file S1, Supplementary Material online.

7The C. rosea isolate IK726, originally isolated from barley roots in Denmark and evaluated extensively regarding its biocontrol efficiency, was used for genome sequencing. The isolate was grown in potato dextrose broth (PDB; Thermo Scientific Oxoid, UK) for DNA extraction and on solid Vogel’s minimal medium with 1% (w/v) sucrose for RNA extraction. DNA was extracted according to established methods and RNA was extracted using the RNeasy Plant Mini Kit (Qiagen, Hilden, Germany) according to manufacturer’s instructions.

8Base coverage of the C. rosea genome was generated using Illumina HiSeq paired end sequencing with an insert length of 0.5 kb and read length of 100 bp using standard library preparation kits. This datum was complemented with a mate pair library with 4 kb inserts sequenced on the Life Technologies SOLiD instrument with Exact Call Chemistry to generate sequence reads of 61 and 53 bp, forward and reverse, respectively. Illumina reads were quality trimmed using Nesoni clip (www.vicbioinformatics.com/software.nesoni, last accessed March 2014). The genome was then assembled with ABySS v. 1.3.3 with a k-mer length of 41.

9Protein-coding genes in the C. rosea genome were annotated using MAKER 2. We configured MAKER 2 to use SNAP, Augustus, and GenemarkES for ab initio gene calls. Two RNA libraries of pure fungal culture were sequenced on the Illumina HiSeq, de novo assembled with Trinity and then provided to MAKER 2 as expressed sequence tag evidence. The evidence set was also complemented with the predicted proteomes of T. virens, T. atroviride, T. reseei, and Neurospora crassa. Assembly sequences and annotation data for C. rosea IK726 is available from www.slu.se/Clonostachys-rosea-IK726 (last accessed January 13, 2015).

10The program Composition Vector Tree (CVTree) version 4.2.1 was used for construction of a phylogenomic tree by using whole predicted proteome sequences. Ten species were included for the topology construction, including C. rosea (isolate IK726), Fusarium graminearum (isolate PH-1), Fusarium oxysporum forma specialis (f. sp.) lycopersici (isolate FOL 4287), Fusarium solani (isolate 77-13-4), Fusarium verticillioides (isolate 7600), T. atroviride (isolate IMI 206040), T. reesei (isolate QM6a), T. virens (isolate Gv29-8), N. crassa (isolate OR74A), and Magnaporthe oryzae (isolate 70-15). Except for C. rosea, the predicted proteome sequences of the included species were obtained from the Comparative Fungal Genomics Platform. A bootstrap test was performed by randomly resampling 63.21% of the proteomes for each species 100 times, as described previously. The distance matrices generated by CVTree were converted into neighbor-joining trees by using “Neighbor” in the PHYLIP v. 3.6 package. A consensus tree of the resulting 100 CVTrees was calculated by using MEGA5 to evaluate bootstrap support of the topology.

11Branch lengths of the species topology generated by CVTree was calculated with PAML v. 4.4, based on a five-gene alignment including actin, calmodulin, glyceraldehyde 3-phosphate dehydrogenase, DNA-directed RNA polymerase II subunit B, and translation elongation factor 1 alpha. Coding gene sequences were retrieved from the respective genome sequences. Each gene was aligned individually using ClustalW, and all alignments were then concatenated. The resulting alignment, containing all five genes, was used to calculate branch lengths using codeml in the PAML package (using a global model for dN/dS). The resulting species phylogeny was calibrated to the fossil record by setting the base of the Sordariomycetes to 335 Ma, according to the fossil recalibration by.

12For manually annotated gene families, the predicted proteomes from individual fungal genome sequences were screened using BLASTP in an iterative process, using reference proteins that covered the diversity of the respective gene family (supplementary file S1, Supplementary Material online). OrthoMCL v. 2.0.8 was used to cluster protein sequences from the included fungal proteomes. Gene family evolution analysis was performed on families that contained ≥2 genes in at least one species and were present in ≥2 species. The program CAFE (Computational Analysis of gene Family Evolution) v. 3 was used to test whether gene family sizes were compatible with a stochastic birth and death model to estimate gene family size in extinct species and to identify lineages with accelerated rates of gene gain or loss. Mutation rate (λ) was estimated from the data and was 0.0024. A separate analysis on ABC transporter subgroups was performed and included data retrieved from T. hamatum, T. longibrachiatum, and. Branch lengths of the species topology were based on a five-gene alignment as described above, and λ was 0.001.

13ABC transporter amino acid sequences were aligned using ClustalW and phylogenetic analyses performed using Neighbor-Joining implemented in MEGA6. The JTT (Jones, Taylor, and Thorton) amino acid substitution model was used with uniform rates among sites, and pairwise deletion of gaps. Statistical support for phylogenetic grouping was assessed by 100 bootstrap resamplings.

14The in vitro effect of the fungicide Boscalid (BASF Hellas, Greece) on mycelial growth of C. rosea on PDA (Thermo Scientific Oxoid) was studied in triplicates in order to assess half maximal effective concentration (EC50). In dual culture assays, Pseudomonas chlororaphis strain MA342 (or sterile water in controls) was streaked at a 6-cm distance from an agar plug of C. rosea or Microdochium nivale (isolate MG1) in a 9-cm-diameter PDA or vegetative peptone agar (VPA, Thermo Scientific Oxoid) plate. Fungal growth was measured in triplicates as the distance between the inoculation point and the hyphal front after 8 days. Percentage fungal growth rate reduction by secreted P. chlororaphis metabolites was calculated as: (1 − [mean growth in P. chlororaphis plates/mean growth in control plates]) × 100.

15Phenazine-1-carboxamide (PCN) was isolated from P. chlororaphis strain PCL1391 using Thin Layer Chromatography (TLC) as described previously. The TLC plate was covered with a thin layer of PDA containing 1 × 106 C. rosea or F. oxysporum f. sp. radicis lycopersici conidia/ml, and assessment of growth on the PCN-containing area was done in triplicates. Biocontrol efficiency was evaluated using a gnotobiotic sand system, with minor modifications. Plant nutrient solution previously inoculated with 5 × 102 F. oxysporum f. sp. radicis lycopersici conidia/ml was mixed with sand to give a final concentration of 50 conidia/g sand. Application of BCAs was done by dipping surface-sterilized tomato pregerminated seeds in a suspension of 1 × 107 P. chlororaphis strain PCL1391 cells/ml, 1 × 106 C. rosea conidia/ml, or both BCAs mixed together, prior to sowing into the sand. Ten days after sowing, assessment of disease severity was carried out in 20 replicates according to the following disease index scale: 1 = no visible symptoms, 2 = mild symptoms on roots, 3 = severe symptoms on roots and wilting, 4 = dead plants. Data were analyzed statistically by analysis of variance and compared using the Duncan test at P = 0.05, using the program SPSS version 16 (SPSS, IL).

16Colonization of tomato roots by a green fluorescent protein (gfp)-expressing C. rosea IK726 isolate was monitored using a Nikon D-Eclipse C1 confocal microscope, using the default filter set. Digital images were acquired with the manufacturer’s software. Fusarium oxysporum f. sp. radicis lycopersici was transformed to express red fluorescent protein (rfp) as described previously, and interactions with gfp-expressing C. rosea were studied in vitro and in planta.

17Gene expression of 11 ABC transporters (supplementary file S2, Supplementary Material online) was measured using quantitative reverse transcription polymerase chain reaction (RT-qPCR). Two Fusarium mycotoxins deoxynivalenol (DON) and zearalenone (ZEN), the fungicide Boscalid, and a 2 days postinoculation (dpi) culture filtrate of P. chlororaphis strain MA342 were used as exogenous toxic substances to generate environmental stresses with relevance for biocontrol conditions. Clonostachys rosea strain IK726 was cultured in 20 ml liquid Czapek-Dox media (Sigma-Aldrich, MO) in 250-ml Erlenmeyer flasks. After 3 dpi at 25 °C, the fungus was subjected to treatment with DON, ZEN, or Boscalid at a final concentration of 200 ppm. In the control treatment, DON and ZEN were replaced with an equal volume of methanol whereas Boscalid was replaced by water. For the P. chlororaphis culture filtrate treatment, C. rosea was cultured in 100 ml PDB for 7 dpi, harvested by filtration and placed in 100 ml P. chlororaphis culture filtrate or fresh VPB as control treatment. Fungal mycelia were subsequently collected after 2-h incubation by vacuum filtration. The harvested mycelia were frozen with liquid nitrogen and kept at −80 °C.

18Total RNA was extracted from mycelia using RNeasy Plant Mini kit (Qiagen) according to the manufacturer’s protocol. To remove DNA impurities, treatment of DNase I and Ribolock RNAase inhibitor (Fermentas, Germany) was performed prior to cDNA synthesis using iScript cDNA synthesis kit (Bio-Rad, CA) following the manufacturer’s instructions. RT-qPCR was carried out as described previously, except that each treatment was measured in a minimum of three biological replicates, each based on two technical replicates. Gene expression data were analyzed by Student’s t-test with 95% confidential interval implemented in Statistica v.10 (StatSoft, OK).

19The C. rosea genome (EMBL: PRJEB4200) was deep sequenced using a short read, whole-genome shotgun approach to at least 150× coverage and assembled with ABySS. The total assembly length was 58.3 Mb spread across 591 scaffolds longer than 2 kb (table 1). A total of 14,268 gene models coding for proteins were annotated using MAKER 2. RNA sequencing data as well as proteomes of related species were used for gene model prediction in MAKER 2. Descriptive names were assigned to 7,385 gene models based on homology to the Uniprot database where the best hit with E value below 1 × 10−20 was used. Finally, all gene models were functionally annotated using InterProScan5 to assign Gene Ontology (GO) terms and Pfam domains (table 1, fig. 1). Pfam domains were found in 11,063 gene models, whereas GO terms were assigned to 8,918 genes.

20In order to resolve the phylogenetic relationship between Bionectriaceae, Nectriaceae, and Hypocreaceae within the order Hypocreales, a whole-proteome analysis was conducted with CVTree. The phylogenetic analysis included C. rosea (Bionectriaceae), F. graminearum, F. oxysporum f. sp. lycopersici, F. solani and F. verticillioides (Nectriaceae), T. atroviride, T. reesei and T. virens (Hypocreaceae) and the outgroup species N. crassa (order: Sordariales) and M. oryzae (order: Magnaporthales). Based on the whole-proteome analysis, C. rosea clustered as a sister taxa to Fusarium spp. with 97% bootstrap support, with Trichoderma spp. in a basal position (fig. 2). Among Fusarium spp., F. oxysporum f. sp. lycopersici and F. verticillioides clustered as sister taxa, whereas F. solani formed the most basal lineage. Trichoderma reesei and T. virens were identified as sister taxa, with T. atroviride in a basal position.

21As C. rosea is reported to be homothallic, the structure of the mating type locus was analyzed. The genome of C. rosea contained two mating type (MAT) 1-2-1 orthologs (BN8691_T00013892_1 on scaffold 161 and BN8691_T000141771_1 on scaffold 347). The high sequence similarity of the genes (100%) and scaffolds (99%) suggests that the apparent duplication was an artifact due to the assembly. The predicted mat1-2-1 gene encoded a predicted protein of 242 amino acids with a conserved high mobility group (HMG) domain, with high similarity (E value = 2 × 10−77) to the mating type 2 HMG-protein in Epichloe festucae (GenBank: AEI72619). No clear orthologs to the MAT1-1 proteins were found in C. rosea. The APN2 DNA lyase and the SLA2 cytoskeleton assembly control factor have been reported neighboring the MAT loci in many other ascomycetes. Both genes, apn2 (BN8691_T00013667_1) and sla2 (BN8691_T00013890_1) were predicted in C. rosea. The sla2 gene was located on scaffold 161, upstream from mat1-2-1, and therefore we considered scaffold 161 to contain the MAT1-2 locus in C. rosea. The apn2 gene was located in the end of scaffold 113 and may be adjacent to mat1-2-1, as mat1-2-1 was at the end of scaffold 161. One additional gene (BN8691_T00013891_1) was predicted to be located between sla2 and mat1-2-1, and displayed some sequence similarities with the MAT1-2-3 proteins of Fusarium species, but the highest similarity (E value = 3 × 10−11) was toward a hypothetical protein from T. virens (GenBank: EFY88586).

22Gene family expansions and contractions in C. rosea were compared with nine other Sordariomycete genomes, with varying life strategies (fig. 2). Gene family data in extant species were either manually annotated (154 families), or predicted with OrthoMCL (17,138 families). The program CAFE was used to test whether gene family sizes were compatible with a stochastic birth and death model, for estimation of gene family size in extinct species and for identification of lineages with accelerated rates of gene gains or losses (supplementary file S3, Supplementary Material online).

23In total, 91 gene families were identified to evolve nonrandomly (P ≤ 0.05), although changes in 62 families were restricted to F. oxysporum f. sp. lycopersici, F. verticillioides, and their ancestor (table 2). The genome of C. rosea contained significantly more polysaccharide lyase family 1 (PL1, pectin lyase), auxiliary activity family 3 (AA3, glucose-methanol-choline oxidoreductases), AA9 (lytic polysaccharide monooxygenase), ABC transporter, PKS, cytochrome P450 monooxygenase (CYP) and OrthoMCL family 1001 (PKS), OrthoMCL family 1003 (ABC transporter), and OrthoMCL family 1044 (Ankyrin-repeat protein) genes than predicted under a random model, but significantly fewer carbohydrate-binding family 18 (CBM18, chitin-binding) module-containing genes (table 2). This was in sharp contrast to Trichoderma spp. as no gene family identified as expanded in C. rosea was expanded in any Trichoderma species, with the exception for CYP that was expanded in T. virens and contracted in T. atroviride and T. reesei (table 2). When compared with Fusarium spp., AA3 was identified as expanded in F. solani, ABC transporter was expanded in F. oxysporum f. sp. lycopersici, whereas CYP was expanded in F. oxysporum f. sp. lycopersici and F. solani. The CBM18 and CYP gene families evolved nonrandomly in a majority of the ten included species (seven and ten species, respectively). Fusarium oxysporum f. sp. lycopersici contained the highest number of expanded gene families (79), whereas F. verticillioides contained the highest number of contracted gene families (22).

24The PL1 gene family, which was significantly expanded in C. rosea (17 genes, supplementary file S4, Supplementary Material online) compared with all other included species, contained enzymes with pectin degrading activity (pectate lyase, EC 4.2.2.2; exo-pectate lyase, EC 4.2.2.9; pectin lyase, EC 4.2.2.10) as listed in the CAZy (carbohydrate-active enzymes) database. Furthermore, a more detailed analysis of the AA3 gene family was conducted, as this family consists of four defined subfamilies (AA3_1 through AA3_4; ). This analysis showed that the expansion in C. rosea and F. solani and the contraction in T. reesei and N. crassa took place in AA3_2 (supplementary file S4, Supplementary Material online). The AA3_2 subfamily included two closely related H2O2-producing enzyme classes, aryl-alcohol oxidase (AO, EC 1.1.3.7) and glucose 1-oxidase (GOX, EC 1.1.3.4, ). The significantly lower number of chitin-binding CBM18s, and the low number of peptidoglycan/chitin-binding CBM50s, in C. rosea was related to the low number of GH18 chitinases (EC 3.2.1.14) in this species. More specifically, CBM18s are found exclusively in GH18 subgroup C (killer toxin-like chitinases), but not in GH18 subgroup A (predicted to contain mostly exo-acting chitinases) or GH18 subgroup B (endo-acting chitinases). A more detailed analysis confirmed that C. rosea only contained two B group and two C group GH18 chitinases, whereas the number of group A GH18s (eight genes) was comparable to the other included species (supplementary file S4, Supplementary Material online).

25Clonostachys rosea harbored an impressive repertoire of genes coding for PKS and nonribosomal peptide synthetase (NRPS) secondary metabolite assembly lines, with 31 PKS, 17 NRPS, and 1 PKS–NRPS hybrid (supplementary file S5, Supplementary Material online). Furthermore, eight terpenoid synthase genes were identified. In order to assign the PKS genes to functional categories, a phylogenetic network comprising keto synthase domains of all deduced PKSs from C. rosea, combined with functionally characterized PKSs of other ascomycetes, was performed (fig. 3). This functional clustering revealed that most PKS genes code for highly reducing PKSs, whereas only four PKSs (PKS1–PKS4) clustered with nonreducing (aromatic) PKSs. The phylogenetic analysis placed PKS1 with anthraquinone-type PKSs, whereas PKS2 and PKS4 represented the citrinin-type of tetraketide synthases. PKS3 resembled a dimethyl orsellinic acid synthase.

26The product spectrum of NRPSs encoded in the C. rosea genome was predicted by domain analysis and extraction of amino acid specificity codes (supplementary file S5, Supplementary Material online). These results suggested that C. rosea possess one putative epipolythiodiketopiperazine synthetase (NPS1), which is presumably involved in biosynthesis of the antimicrobial compound glioperazine or a related metabolite. Furthermore, several putative peptaibol assembly lines were identified. The largest peptaibol synthetase in C. rosea (NPS9) comprised 20 biosynthetic modules, whose specificity signatures corresponded well to the structure of antimicrobial peptaibols isolated from C. rosea. However, the majority of NRPS genes belonged to orphan pathways.

27Clonostachys rosea contained not only the highest number of predicted ABC transporters (86 genes, supplementary file S6, Supplementary Material online) among the included species but also the highest number of predicted major facilitator superfamily (MFS) transporters (620 genes, supplementary file S6, Supplementary Material online). The OrthoMCL analysis identified six MFS transporter genes in two families (OrthoMCL12606 and OrthoMCL12615) that were unique for C. rosea (supplementary file S7, Supplementary Material online).

28The ABC transporter gene family is shown to be divided into groups with various functions. As mycoparasites need to defend themselves against toxic metabolites produced by the fungal prey, we hypothesized that the increased number of ABC transporters in C. rosea was related to drug resistance. Therefore, a more detailed evolutionary analysis of ABC transporters in filamentous ascomycetes was performed. The analysis of gene family evolution estimated the number of ABC transporter genes in ancestral species (fig. 4) and number of gene gains and losses for each ABC transporter group in each branch of the species tree (supplementary file S6, Supplementary Material online).

29The analyses showed that the expansion of ABC transporters in C. rosea took place in group B (multidrug resistance proteins) and in group G (pleiotropic drug resistance proteins) (P ≤ 0.001, table 3). ABC transporters in groups B and G are further divided into subgroups, and phylogenetic analyses showed that most group B ABC transporters in C. rosea were found in subgroup B-III (fig. 5), which include the M. oryzae B-III member Abc3 that is partly required for resistance toward valinomycin, actinomycin D, and reactive oxygen species. The expansion in group G involved subgroups G-I and G-V (fig. 6), which include many members with established roles in drug resistance.

30In Trichoderma, a significant increase (P ≤ 0.001) of ABC transporter genes was identified for T. virens, which could be traced to group ABC C (multidrug resistance-associated proteins), subgroup C-V (supplementary file S8, Supplementary Material online) that contain members hypothesized to be involved in export of endogenous toxins. In contrast, T. reesei and the ancestor to Trichoderma experienced a significant (P ≤ 0.032) loss of group G ABC transporters (table 3). Significant losses (P ≤ 0.054) of ABC transporter genes were also detected in Coccidioides immitis, F. verticillioides and N. crassa, whereas an expansion (P ≤ 0.003) was evident in Aspergillus oryzae and F. oxysporum f. sp. lycopersici (table 3).

31The expansion of ABC transporters associated with drug resistance in C. rosea prompted us to investigate whether expression of ABC transporter genes was induced in response to mycotoxins DON and ZEN, the fungicide Boscalid and metabolites from the BCA bacterium P. chlororaphis strain MA342. First, the ability of C. rosea to tolerate Boscalid was evaluated by measuring in vitro growth on agar plates, and the resulting EC50 value was greater than 500 µg/ml. Second, the ability of C. rosea to tolerate metabolites produced by P. chlororaphis strain MA342 was evaluated using a dual culture assay. The growth inhibition of the snow mould pathogen Mi. nivale by secreted P. chlororaphis MA342 metabolites was significantly (P ≤ 0.05) higher than the growth reduction of C. rosea under the same conditions, on both PDA and VPA (supplementary file S9, Supplementary Material online).

32Expression of a subset of 11 genes representing groups ABC B, ABC C, and ABC G was analyzed (supplementary file S2, Supplementary Material online). Two ABC transporter genes were significantly (P ≤ 0.043) induced by the ZEN treatment; abcB26 and abcG8 (fig. 7). Gene abcB3 was induced (P = 0.012) by the Boscalid treatment whereas two genes (abcC12 and abcC14) were repressed (P ≤ 0.048) (fig. 7). Five genes were induced (P ≤ 0.037) by exposure to a culture filtrate from P. chlororaphis; abcB1, abcB4, abcB18, abcB26, and abcG8 (fig. 7). No gene expression changes were induced by the DON treatment.

33The high tolerance of C. rosea toward metabolites from P. chlororaphis suggested the potential for combination treatments of these two BCAs for additive biocontrol effects. For this experiment, P. chlororaphis strain PCL1391 was used. The ability of C. rosea to tolerate the main metabolite from PCL1391, PCN, was evaluated. PCN completely inhibited in vitro growth of the tomato foot and root rot pathogen F. oxysporum f. sp. radicis lycopersici, whereas growth of C. rosea was observed under the same condition (supplementary file S9, Supplementary Material online). A gnotobiotic sand system was used to evaluate the effect of coinoculating C. rosea together with P. chlororaphis PCL1391 on tomato seeds to protect against F. oxysporum f. sp. radicis lycopersici foot and root rot disease. The combination treatment of the two BCAs reduced the disease index to a significantly (P ≤ 0.05) lower level than either the single C. rosea or the single P. chlororaphis treatments (supplementary file S9, Supplementary Material online).

34Colonization of tomato seedling roots by C. rosea IK726 in the gnotobiotic sand system was studied by fluorescence microscopy using a strain expressing gfp. Clonostachys rosea colonized the plant’s roots by growing along the junctions of epidermal cells, formed a net of hyphae around the main root and conidiated 6 dpi (supplementary file S9, Supplementary Material online). Hyphae of C. rosea were frequently observed intermingled, attached and coiled around root hairs, and hyphal swellings putatively representing penetration structures were repeatedly observed (supplementary file S9, Supplementary Material online), although no disease symptoms were observed. In the presence of F. oxysporum f. sp. radicis lycopersici both fungi were observed on the root surface, and close contact of hyphae indicating mycoparasitism was present at several sites. A potential for mycoparasitism was further suggested by in vitro interactions where hyphae of C. rosea IK726 were found to tightly coil around hyphae of F. oxysporum f. sp. radicis lycopersici (supplementary file S9, Supplementary Material online).

35Our phylogenomic analysis shows that the Bionectriaceae family (represented by C. rosea) is a sister taxon with the family Nectriaceae (Fusarium spp.), with Hypocreaceae (Trichoderma spp.) in a basal position, thereby supporting the results from. Our data also support the close relationship between Trichoderma sections Longibrachiatum (T. longibrachiatum and T. reesei) and Pachybasium (T. virens). Mycoparasitism was previously shown to be the ancestral life style of Trichoderma spp., but the mycoparasitic lifestyle within Hypocreales may be even older. In fact, certain fossils from the Early Devonian (400 Myr) Rhynie chert sediments depict coiling hyphae surrounding partially degraded perithecia, which are interpreted as remnants of mycoparasites. Whether mycoparasitism in C. rosea represents an ancestral trait shared by Clonostachys spp. and Trichoderma spp., or represents a case of convergent evolution, requires better understanding of the mycoparasitic potential of species within Bionectriaceae.

36We hypothesized that the mycoparasitic lifestyle in Bionectriaceae and Hypocreaceae results in selection for similar interaction mechanisms. However, the majority of expanded gene families in C. rosea did not evolve in the same manner in Trichoderma spp., indicating differences in the underlying mechanisms of antagonism. One major difference concerns the ABC transporter gene family, where the lifestyle of C. rosea results in selection for increased numbers of genes in groups B-III, G-I, and G-V, presumably reflecting the need for tolerance toward toxic metabolites produced by the mycoparasite itself, the fungal prey, or plants. Subgroup ABC G-I includes several members with characterized roles in drug resistance, such as Pdr5p, Pdr10p, and Pdr15p in Saccharomyces cerevisiae. Interestingly, the TAABC2 transporter from T. atroviride belongs to the ABC G-V group, and is responsible for resistance toward exogenous toxins from other fungi. The ABC transporter genes abcG5 and abcG29 (groups G-I and G-V, respectively) from C. rosea are induced by ZEN and deletion of abcG5 results in reduced ZEN tolerance and partial loss of the ability to protect barley against mycotoxigenic F. graminearum.

37Deletion of the group B-III gene abc3 in M. oryzae results in mutants sensitive to valinomycin and actinomycin D, but highly sensitive to oxidative stress. The M. oryzae Δabc3 mutants also contained higher endogenous levels of H2O2, and if this function can be proven also for C. rosea B-III ABC transporters it is interesting to speculate about a connection between the expansions of ABC B-III genes and the AA3_2 subfamily genes encoding H2O2-producing enzymes. Group C-V ABC transporters are hypothesized to be involved in secretion of endogenous toxins due to their close association with secondary metabolite biosynthesis clusters. The increased numbers of ABC C-V genes in T. virens, but not C. rosea, may therefore be connected with the dependence on gliotoxin production in this species.

38The induced expression of several group B and G ABC transporters in C. rosea during exposure to ZEN, Boscalid, and bacterial metabolites indicates that membrane transporters may contribute to the high tolerance against these compounds in C. rosea. We show that C. rosea can tolerate at least 10-fold higher levels of Boscalid than the gray mould pathogen Botrytis cinerea, and that C. rosea can tolerate toxic compounds produced by different strains of the BCA bacterium P. chlororaphis to a substantially higher degree than certain plant pathogenic species such as Mi.nivale and F. oxysporum f. sp. radicis lycopersici. These results may find important agro-industrial application areas, and we show that the combined application of both C. rosea and P. chlororaphis PCL1391 on tomato seeds protects the emerging seedlings better against root rot caused by F. oxysporum f. sp. radicis lycopersici, compared with single application of either BCA.

39Mycoparasitic Trichoderma species contain high numbers of PKS and NRPS genes, and transcriptome data support the view that production of secondary metabolites that poison the fungal prey is important for both T. virens and T. atroviride, although T. atroviride also induces genes involved in hydrolytic degradation of the cell wall of the fungal prey. In this context, the high numbers of PKS and NRPS genes in the C. rosea genome may suggest a prolific capacity of C. rosea to produce, hitherto mostly unknown, polyketide natural products and peptide metabolites. The putative NPS9 peptaibol synthetase in C. rosea displays a specificity signature that corresponds well with peptaibols isolated from C. rosea strain BAFC3874, which contribute to the antifungal activity of C. rosea culture filtrates toward Sclerotinia sclerotiorum.

40However, very few metabolites are in fact isolated from C. rosea. Besides the mentioned peptaibols, Clonostachys sp. produce some epipolysulfanyldioxopiperazines with nematicidal activity, polyterpenoid glisoprenins that inhibit appressorium formation of phytopathogenic fungi, and cyclic peptides with chitinase inhibitor activity that may be involved in protecting the fungus from chitinases produced by itself or by the fungal prey. Certain GH18 subgroup C killer toxin-like chitinases are hypothesized to make the cell wall of the fungal prey more permeable to enhance penetration of toxins, and the high subgroup C gene numbers in T. virens and T. atroviride (15 and 9 genes, respectively) are hypothesized to represent an adaptation to the dependence on toxic metabolites during the mycoparasitic attack. Consequently, the low number (two genes) of subgroup C chitinase genes in C. rosea suggests that the products of the orphan secondary metabolite assembly lines in C. rosea are not primarily involved in chemical attack.

41We hypothesize that certain PKSs are involved in biosynthesis of the pink, yellow, orange, and brown pigments frequently observed in C. rosea cultures as various chromophores, for example, polyenes or aromatic systems can be produced by reducing or nonreducing PKSs, respectively. Deletion of the pks4 gene in T. reesei results not only in loss of green conidial pigmentation but also in lower resistance toward toxic metabolites produced by Alternaria alternata, Rhizoctonia solani, and S. sclerotiorum. Many PKS genes in T. reesei are highly expressed during vegetative growth, which suggest more diverse functional roles than merely in biosynthesis of secondary metabolites. Two PKS genes of the reducing, lovastatin/citrinin-like type are induced in T. atroviride during interaction with R. solani, but the identity and function of the resulting products are not known.

42In relation to T. atroviride and T. virens, and even the saprotrophic T. reesei, the fungal cell wall degrading capacity of C. rosea appears to be modest. The three Trichoderma spp. contain high numbers of GH18 chitinases, GH75 chitosanases, and β-1,3-glucanases (GH17, GH55, GH64 and GH81), whereas the numbers in C. rosea are in most cases lower. The mycoparasitic lifestyle of T. atroviride and T. virens is correlated with selection for high GH18 subgroup B endochitinase gene numbers. It is therefore significant to observe that C. rosea only possesses two B group chitinases, from which one is the ortholog to a conserved ascomycete cell wall bound chitinase (Chi18-18 in T. reesei, chit-1 in N. crassa) with a cell wall plasticizing role during growth. This suggests that cell wall degradation of the fungal prey may not be a prominent biocontrol trait in C. rosea. We must recognize, however, that a low number of genes in a particular isozyme gene family does not necessarily imply a low capacity to degrade a particular substrate, but may suggest a narrower substrate range. For example, the genome of the well-known cellulase and hemicellulase producer T. reesei contains fewer cellulase and hemicellulase-encoding genes than most species. Deletion of the C. rosea GH18 chitinase genes ech37, ech42 (also referred to as chi1), and ech58 results in lower in vitro antagonistic activity toward Fusarium culmorum, but has no effect on its biocontrol ability toward F. culmorum on barley or Alternaria radicina on carrots.

43Sequencing of the C. rosea genome and establishment of its taxonomic relationship with mycoparasitic Trichoderma spp. have provided a unique opportunity to study key similarities and differences between their respective life strategies, which have direct implications for the implementation of biocontrol in agricultural production systems. A unique feature in C. rosea is the high numbers of ABC transporters in subgroups predicted to be involved in drug resistance, which may be exploited in agro-industrial applications to reduce the amount of used pesticides.

44Supplementary files S1–S9 are available at Genome Biology and Evolution online (http://www.gbe.oxfordjournals.org/).