Interspecific and host-related gene expression patterns in nematode-trapping fungi
© Andersson et al.; licensee BioMed Central Ltd. 2014
Received: 12 February 2014
Accepted: 24 September 2014
Published: 11 November 2014
Nematode-trapping fungi are soil-living fungi that capture and kill nematodes using special hyphal structures called traps. They display a large diversity of trapping mechanisms and differ in their host preferences. To provide insights into the genetic basis for this variation, we compared the transcriptome expressed by three species of nematode-trapping fungi (Arthrobotrys oligospora, Monacrosporium cionopagum and Arthrobotrys dactyloides, which use adhesive nets, adhesive branches or constricting rings, respectively, to trap nematodes) during infection of two different plant-pathogenic nematode hosts (the root knot nematode Meloidogyne hapla and the sugar beet cyst nematode Heterodera schachtii).
The divergence in gene expression between the fungi was significantly larger than that related to the nematode species being infected. Transcripts predicted to encode secreted proteins and proteins with unknown function (orphans) were overrepresented among the highly expressed transcripts in all fungi. Genes that were highly expressed in all fungi encoded endopeptidases, such as subtilisins and aspartic proteases; cell-surface proteins containing the carbohydrate-binding domain WSC; stress response proteins; membrane transporters; transcription factors; and transcripts containing the Ricin-B lectin domain. Differentially expressed transcripts among the fungal species encoded various lectins, such as the fungal fruit-body lectin and the D-mannose binding lectin; transcription factors; cell-signaling components; proteins containing a WSC domain; and proteins containing a DUF3129 domain. A small set of transcripts were differentially expressed in infections of different host nematodes, including peptidases, WSC domain proteins, tyrosinases, and small secreted proteins with unknown function.
This is the first study on the variation of infection-related gene expression patterns in nematode-trapping fungi infecting different host species. A better understanding of these patterns will facilitate the improvements of these fungi in biological control programs, by providing molecular markers for screening programs and candidates for genetic manipulations of virulence and host preferences.
Soil contains a diverse range of fungi that are parasites on nematodes . These fungi include the nematode-trapping fungi, which have specific hyphal structures in which the nematodes can be trapped. The interest in studying these fungi is due to their potential use as biological control agents for plant and animal parasitic nematodes . The traps of the nematode-trapping fungi develop from hyphae and can be formed spontaneously or be induced in response to signals from the environment . There is a large variation in the morphology of the traps, and the type of trap depends on the species . In some species, the traps consist of an erect branch that is covered by an adhesive material. In other species such as in the well-studied Arthrobotrys oligospora, the trap is a three-dimensional net. A third type of trap is the adhesive knob, which is a single celled structure. Finally, there are some species that capture nematodes using a mechanical trap called a constricting ring . Despite the large morphological variation in trapping structures, phylogenetic analyses inferred from molecular data has shown that the majority of the nematode-trapping fungi belong to a monophyletic group consisting of a single family of the order Orbiliales (Ascomycota) [5–8]. Furthermore, these studies have shown that the nematode-trapping fungi have evolved along two major lineages: one basal lineage leading to species with constricting rings and one lineage containing species that form adhesive traps, including three-dimensional networks, knobs and branches [5–8].
The trapping mechanisms of the species with constricting rings and those with adhesive traps are distinctively different . The constricting ring consists of three cells. When the nematode enters the ring, the cells inflate and the nematode is trapped. The closure is very rapid (0.1 s) and is triggered by pressure of the nematode on the constricting-ring cells . Ultrastructural examinations revealed that the cell wall of the constricting-ring cells is folded; when the cells inflate, the folded cell wall balloons out and forms the new cell wall [10, 11]. The adhesive trap is surrounded by a layer of fibrillar, extracellular polymers. Although the molecular mechanism has not yet been characterized, ultrastructural studies have shown that the fibrillar layer is reorganized during the attachment of the traps to the nematode cuticle . Following the trapping of nematodes, the infection mechanisms appear to be rather similar in the species with constricting rings and adhesive traps: the fungus forms a penetration tube that pierces the nematode cuticle. During penetration the nematode becomes paralyzed. Subsequently, the internal tissues are rapidly colonized and digested by fungal hyphae .
In laboratory assays, most nematode-trapping fungi can trap and infect a range of different nematode species [14, 15]. However, there are a number of studies showing that different species and even strains of nematode-trapping fungi can vary in their host preferences. For example, in vitro predacity tests of four nematode-trapping fungi showed that the constricting-ring species Arthrobotrys dactyloides was the most efficient species in capturing and killing the root knot nematode Meloidogyne graminicola. Further studies on A. dactyloides showed that even strains of this species differed in their predacity to Meloidogyne incognita. Significant differences in the susceptibility to nematode-trapping fungi have also been shown in field trials with the cyst nematode Heterodera schachtii and the root knot nematode Meloidogyne javanica.
Recently, the infection mechanism of nematode-trapping fungi has been examined using the tools of genomics, transcriptomics and proteomics. The genomes of two nematode-trapping fungi have been published; the net-forming Arthrobotrys oligospora and the adhesive knob-forming Monacrosporium haptotylum. The two genomes are similar in size and consist of ~62% core genes that are shared with other fungi, ~20% genes that are specific for the two species and ~16% genes that are unique for each genome . Comparative genome analysis showed that the genomes of nematode-trapping fungi have been expanded in a number of gene families, including extracellular peptidases such subtilisins; homologs to several virulence factors identified in plant-pathogenic fungi; and families of putative cell-surface adhesins containing carbohydrate-binding domains such as the WSC domain and the mucin domain . Transcriptome analysis showed that M. haptotylum expresses a unique set of genes during the early stages of infection of the nematode Caenorhabditis briggsae. Among these is a large proportion that belongs to gene families that are significantly expanded in the nematode-trapping fungi. Transcripts encoding small secreted proteins (SSPs) and many species-specific genes were also highly expressed during the early phase of infection. Many of them were orphans, that is, genes lacking both homologs and Pfam domains [20, 21]. Furthermore, quantitative proteomics revealed proteins that were significant upregulated in the knob compared with the vegetative mycelia in M. haptotylum. Among the upregulated proteins were peptidases, tyrosinase and proteins containing the WSC domain .
In this study, we have examined in more detail the molecular basis of the infection process in nematode-trapping fungi that have various trapping mechanisms, including adhesive nets (A. oligospora), adhesive branches (Monacrosporium cionopagum) and constricting rings (A. dactyloides) [4, 13]. Two plant-parasitic nematodes were used as hosts, the root knot nematode Meloidogyne hapla and the sugar beet cyst nematode H. schachtii. These are both sedentary endoparasites belonging to the group of nematodes that causes the most damages to crops . Comparative transcriptome analysis of the infection, including the adhesion, penetration and digestion stages, showed that the divergence in interspecific gene expression was significantly larger than that related to the nematode host used. We identified a common set of genes that were expressed by all three fungi and a more variable set that were regulated depending on either the fungal species or the nematode host.
Results and discussion
Infection of plant parasitic nematodes by nematode-trapping fungi a
Characterization of the transcriptome libraries
Characterization of the transcriptome libraries
Total number of reads
Total number of isotigs
(Average size, bp)
Isotigs with Pfam
Number of unique Pfam
Isotigs with UniRef50
Mapping to Ao genome
Number of mapped reads
Number of gene modelsd
Comparing transcriptomes between species with not yet sequenced genomes is challenging due to the difficulties in identifying one-to-one orthologs. To circumvent this problem, we here clustered the isotigs into UniRef50 clusters. UniRef50 clusters are based on pre-computed sequence clusters of the UniProt database that have at least 50 percent similarity and 80 percent coverage . Istotigs were only grouped to UniRef50 clusters if passing a given threshold value (1e-10) and only one isotig per fungal species (displaying the highest sequence similarity) were assigned to a given UniRef50 cluster. A recent study including data from seven fungal genomes revealed that the grouping of gene sequences into UniRef50 clusters using the described procedure are in close agreement with traditional ortholog clustering methods (Canbäck et al., manuscript in preparation). Furthermore, the risk of clustering non-orthologous gene duplicates into a given UniRef50 cluster is reduced in nematode-trapping due to the rapid divergence of gene duplicates generated by repeat induced point (RIP) mutations [19, 20]. In our analyses, 15 713 of the in total 17 446 isotigs matched to 6 520 unique UniRef50 protein clusters. The reads of these putative orthologs were normalized with DESeq . Based on the hypothesis that most transcripts are not differentially expressed, the analysis showed that a proper normalization was obtained using the 5% most highly expressed UniRef50 clusters (Additional file 2). This cohort (“Differentially expressed UniRef50 clusters”) contained 326 unique UniRef50 clusters.
The third data set, “Host-specific gene expression”, was generated for identifying genes that were differentially expressed due to the nematode host species. Because the genome of A. oligospora is available , the analysis focused on comparing the transcriptional response of this fungus when infecting M. hapla and H. schachtii. The reads from the Ao(Mh) library were mapped to 7 351 genes and those from the Ao(Hs) library to 6 377 genes (Table 2).
Divergence in gene expression
To compare the functional groups of genes that were expressed in the different libraries, the abundances and expression levels of Pfam domains in the “Top 500 transcripts” data set were analyzed. The number of Pfam domains found in the five libraries varied between 330 and 412 (Additional files 3 and 4). In total, 700 Pfam domains were found in at least one of the libraries.
Commonly expressed transcripts
Proportion (%) of secreted proteins and orphans a
Top 500 most expressed isotigs
Top 500 most expressed isotigs
Pfam domains expressed by all fungi during nematode infection a
Peptidase_S8 (PF00082); Aspartyl protease (PF00026); Peptidase_M3 (PF01432); Peptidase_S10 (PF00450);
CFEM (PF05730); Ricin-type lectin (PF14200)
Stress response, chaperons
Thioredoxin (PF00085); Glutathione S-transferase (PF00043); Catalase (PF00199); AhpC/TSA family (PF00578); DnaJ (PF00226); HSP20 (PF00011); HSP70 (PF00012); HSP90 (PF00183); Peptidylprolyl isomerase (PF00254, PF00160)
RHO protein GDP dissociation inhibitor (PF02115); ADP Ribosylation Factors (ARFs) (PF00025); 14-3-3- proteins (PF00244); Pkinase (PF00069); Ras (PF00071)
Calreticulin (PF00262); EF-hand motif (PF13499)
Actin (PF00022); Cofilin (PF00241); Profilin (PF00235); Tropomyosin (PF12718); Tubulin (PF00091)
Membrane protein Mpv17_PMP22 (PF04117)
Multiprotein bridging factor 1 (PF08523); Homeobox (PF00046); TATA binding protein (PF00352); bZIP Transcription factor (PF00170, PF07716); Histone (PF00125, PF00538); Nucleosome assembly protein (PF00956); Helicases (PF00270, PF00271)
ABC transporter (PF00005) Amino acids permease (PF00324); Ammonium transporter (PF00909); Major Facilitator Superfamily (PF07690); Sugar and others (PF00083); Porin (PF01459)
Syntaxin (PF05739); Synaptobrevin (PF00957); Rab GDP dissociation inhibitors (PF00996)
Enolase (PF03952); Fructose-bisphosphate aldolase (PF01116); Triose-phosphate isomerase (PF00121); Glyceraldehyde 3-phosphate dehydrogenase (PF00044); Phosphoglycerate kinase (PF00162); Pyruvate kinase (PF00224); Pyruvate carboxylase (PF00682)
Metabolism, pentose phosphate pathway
Phosphogluconate dehydrogenase (PF03446, PF00479); Transaldolase (PF00923); Transketolase (PF00456)
lactate/malate dehydrogenase (PF02866, PF00056); Succinyl coenzyme A synthetase (PF00549)
Cyt-b5 (PF00173); ATP synthase (PF00887, PF00137); mitochondrial carrier (PF00153); ATPases (PF00006, PF02874)
Metabolism, amino acids
Glutamine amidotransferases (PF00310); NAD-specific glutamate dehydrogenase (PF10712); Glutamine synthetase (PF03951)
Glycoside hydrolase family 1 (PF00232); Biotin-requiring enzyme (PF00364) aldo-keto reductase family (PF00248); Aldehyde dehydrogenase family (PF00171); Short-chain dehydrogenases/reductases family (PF00106); Pyrophosphatase (PF00719); Myo-inositol-1-phosphate synthase(PF01658, PF07994); FA desaturase (PF00487); Nucleoside diphosphate kinase (PF00334); UDP-glucose pyrophosphorylase (PF01704); Transketolase, pyrimidine binding domain (PF02779); ATPase (PF00004); CoA binding domain (PF02629) ATP-grasp_2 (PF08442); Calcineurin-like phosphoesterase (PF00149); Epimerase (PF01370)
Proteasome (PF00227, PF10584); Ubiquitination (PF00240, PF00240, PF00179)
Elongation factors (PF10587, PF01873, PF00009); tRNA synthetases (PF00587); translation, initiation factors W2 (PF02020), SUI1 (PF01253), MIF4G (PF02854); Ribosomal proteins (PF00428, PF00428, PF00687, PF00466, PF00238, PF00827, PF00252, PF14204, PF01775,PF00828, PF00861, PF01280, PF03947, PF01157, PF01776, PF01246, PF01777, PF00831, PF00297, PF01198, PF01655, PF01780, PF01907, PF01020, PF00935, PF00281, PF00347, PF01159, PF01248, PF00411, PF00164, PF00416, PF00253, PF00833, PF00203, PF01090, PF00318, PF01282, PF03297,PF01283, PF01667, PF01015, PF00163, PF00333, PF00177, PF00410, PF01201, PF00380, PF08071)
Furthermore, transcripts containing the Atg8 domain were highly expressed in all fungi. Atg8 is an essential protein in the autophagic pathway  and disruption of a homolog of this gene in A. oligospora leads to reduced trap formation . In addition, all three fungi expressed transcripts with the Pfam domain RicinB_lectin_2 (PF14200). Ricin-B lectins are ribosome-inactivating proteins (RIPs) consisting of a catalytic A-chain and a sugar-binding B-chain [31, 32]. All fungi have a RicinB_lectin_2 transcript that match to G1X3G7 in A. oligospora. G1X3G7 is a protein with a length of 134 amino-acid residues (aa), without a secretion signal and with low sequence similarity to other RicinB lectins in the UniProt database. A RicinB_lectin_2 domain-containing protein (MOA) with nematotoxic activity against Caenorhabditis elegans has been identified in the basidiomycete Marasmius oreades. The nematotoxicity was dependant on the cysteine protease activity of MOA and the binding of its lectin domain to glycosphingolipids in the worm intestine. MOA consists of 293 aa and lacks a classical secretion signal . Sclerotinia sclerotiorum agglutinin (SSA) is a RicinB_lectin_2 domain-containing protein with a length more similar to the G1X3G7 protein in A. oligospora. SSA has a length of 153 aa, lacks secretion signal and shows insecticidal properties when fed to the pea aphid Acyrthosiphon pisum.
Previous studies have shown that subtilisins (peptidase_S8) are important virulence factors in nematode-trapping fungi. In A. oligospora they have a key role in the early stages of infection, including immobilization of the captured nematode [19, 35, 36]. A. oligospora has 52 genes containing the peptidase_S8 domain . However, only one transcript containing the peptidase_S8 domain was identified among the highly expressed transcripts in each library of Ao(Mh) and Ao(Hs). BLASTX searches showed that both transcripts displayed the highest sequence homology to the A. oligospora protein G1XLL2. The other cDNA libraries contained also only one transcript with the peptidase_S8 domain among the top 500 expressed genes. The three transcripts in Ad(Mh), Ad(Hs) and Mc(Hs) all displayed the highest sequence homology to H072_8474 in M. haptotylum. Interestingly, G1XLL2 and H072_8474 are orthologs (T. Meerupati, B. Canbäck, D. Ahrén, A. Tunlid, manuscript in preparation). Furthermore, G1XLL2 was the most expressed peptidase_S8 gene in A. oligospora during early infection (6 and 10 hours) of C. briggsae, and H072_8474 was the second most expressed peptidase_S8 gene in M. haptotylum during early infection (4 hours) of C. briggsae. H072_8474 was also identified in the proteome of both the knob and the mycelia in M. haptotylum. This shows that despite the large number of peptidase_S8 genes only a few are highly expressed during infection.
Stress proteins were highly expressed in all fungi independent of trapping structure. They included heat-shock proteins and chaperones such as DnaJ, HSP70 and HSP90; gluthatione S-transferases; and antioxidant enzymes such as thioredoxin and catalase. Antioxidants are enzymes involved in the protection of the cell from oxidative damages induced by reactive oxygen species (ROS) . ROS are continuously produced in the cell as byproducts from various metabolic pathways and have an important role(s) in signaling . In the plant-pathogenic fungus Magnaporthe grisea, ROS-generating NADPH oxidases (Nox1 and Nox2) are essential for pathogenicity . The authors  suggest that the generated ROS accumulate in the appressorium to facilitate oxidative cross-linking of cell-wall proteins. This leads to a strengthening of the cell wall of the appressorium that will eventually resist high turgor pressure . Transcripts with sequence similarity to the Nox proteins in M. grisea were regulated in all fungi of our study during infection. However, none of them were found among our top 500 most expressed transcripts.
Interspecific variation in gene expression
Transcripts encoding the fungal fruit-body lectin (FB_lectin) and the D-mannose binding lectin (B_lectin) were only highly expressed in A. oligospora and not in the other two fungi (Figure 5). In the UniRef50 cluster containing the FB_lectin domain we identified the previously studied AOL lectin (Q00233) . The transcript matching to this cluster was the third most expressed of all transcripts during Ao(Mh) infection and the 25th most expressed of all transcripts during Ao(Hs) infection. Earlier studies have shown that AOL functions as a storage protein during both saprophytic and parasitic growth . However, deletion of this gene did not affect the fungus’ ability to infect nematodes .
Previous studies have shown that proteins containing the carbohydrate-binding domain WSC comprise a large and rapidly evolving gene family in M. haptotylum. Phylogenetic analysis of the 33 WSC-containing proteins in M. haptotylum revealed a clade of 15 WSC paralogs . This clade contains only one (G1X6Q5) of the 16 WSC proteins identified in A. oligospora. Thirteen of the 15 WSC paralogs of M. haptotylum were at least twofold upregulated during the infection of the nematode C. briggsae. In this study, transcripts encoding WSC proteins were highly expressed by all fungi during infection of plant-parasitic nematodes (Table 4). The largest number of transcripts encoding WSC domain proteins was expressed by M. cionopagum. In total, 19 transcripts of WSC domain proteins were identified in the Mc(Hs) library, of which seven were found among the top 500 transcripts. Four of these transcripts displayed closest sequence homology to proteins found in the expanded clade of WSC proteins of M. haptotylum. The libraries of A. oligospora and A. dactyloides contained a lower number of transcripts of WSC domain proteins. Among the top 500 transcripts, the Ao(Hs), Ao(Mh) and Ad(Mh) libraries each had two WSC domain proteins, whereas none were found in the Ad(Hs) library. Transcripts displaying highest sequence homology to the previously mentioned A. oligospora protein G1X6Q5 were identified in both the Ao(Hs) and the Ao(Mh) libraries. The deduced proteins of two highly expressed WSC transcripts in A. dactyloides did not show any sequence similarity to the proteins found in the expanded clade of paralogs in M. haptotylum. Taken together, the comparative transcriptome analysis shows that the WSC domain proteins comprise a large and divergent gene family that is highly expressed during pathogenesis in nematode-trapping fungi. The specific sets of genes that are expressed depend on the fungal species and the nematodes being infected, which suggests that the function of the WSC proteins is to contribute to the specialization of the trapping mechanisms.
Virulence associated transcripts
A BLAST search of the top 500 transcripts in each library was conducted in the pathogen–host interaction protein database (PHI-base) . PHI-base contains experimentally verified pathogenicity, virulence and effector genes from fungi, oomycetes and bacterial pathogens. In total, 97 unique PHI-base genes were identified. Genes with sequence similarity to seven of them were found in at least ten gene models in either M. haptotylum or A. oligospora. They included RBT4 from Candida albicans, which is necessary for virulence . The function of RBT4 is unknown but it contains a CAP (Cysteine-rich secretory proteins, Antigen 5 and Pathogenesis-related 1 protein) domain . Among the 97 identified PHI-base genes, 15 were highly expressed by all fungal species (Additional files 8 and 9) and 82 were highly expressed by one or two of the fungal species (Additional file 10). The PHI-base genes expressed by all fungi included stress response genes and several cell signaling genes containing the Ras domain. The PHI-base genes that differ in expression between the fungi included aspartic proteases and the gas 1 and gas 2 proteins of M. grisea that contains the DUF3129 domain.
Host-specific gene expression
Differentially expressed genes encoding secreted proteins in A. oligospora during infection of M. hapla or H. schachtii a
Upregulated in M. hapla
Histidine acid phosphatase
Patched sphingolipid transporter
Glucosidase 2 subunit beta
GPI anchored cell-wall protein
Glyco_hydro_35, BetaGal_ dom2, BetaGal_dom3, BetaGal_dom4_5
Upregulated in H. schachtii
This is the first study examining the variation in gene expression patterns among nematode-trapping fungi during infection of different host species. Comparative transcriptome analysis showed that the divergence in gene expression between the fungal species was significantly larger than that related to the nematode host. A core set of transcripts that were highly expressed by all three fungi was identified. This core set included subtilisins, aspartic proteases and proteins containing the CFEM domain. These genes were also highly expressed by A. oligospora and M. haptotylum during infection of C. briggsae. Also within this core set, a putative nematotoxic protein was identified, the Ricin-B lectin. A more variable set of transcripts being regulated depending on the fungal species was also identified. A small set of genes was identified showing differential expression depending on the host. This set was enriched in genes encoding secreted proteins and also included several gene families that were expanded in the genomes of nematode-trapping fungi . Information on the genetic basis of the interspecific variation in the trapping mechanisms and host preferences fungi will be useful for researchers who are screening for more potent biological control agents of nematode-trapping fungi.
Culture of organisms and infection experiments
Cultures of A. oligospora (ATCC 24927), M. cionopagum (CBS 220.54) and A. dactyloides (CBS 109.37) were maintained on corn meal agar 1:10. Infested soil/roots of M. hapla (Strain E 226) were obtained from Prof. Dr. Gerrit Karssen (Plant Protection Service, HC Wageningen, the Netherlands) to start the culture of this nematode. Small pieces of roots infested with M. hapla were inoculated in rhizosphere of 3-week-old tomato plants raised in the green house of Department of Biology, Lund University, for induction of root knots. After 8 weeks, the infected tomato roots with well developed knots and egg masses were gently washed under running tap water. Egg masses of M. hapla were picked by fine forceps from knots of infected roots under a stereoscopic binocular microscope and surface disinfected for 1 minute in a 0.5% NaOCl solution and rinsed three times in sterile distilled water. Egg masses were then collected in Petri dishes (30 mm) in sterile distilled water and incubated at 22 ± 1°C for 48 hours for hatching of second-stage juveniles (J2). After incubation, freshly hatched J2s were separated from the egg masses and collated in Eppendorf tubes, surface sterilized with 0.5% NaOCl for 2 minutes, and rinsed five times with sterilized distilled water, and used for infection experiments. J2s of H. schachtii were obtained from HZPC in the Netherlands (http://www.hzpc.com) and used for infection experiments after sterilization and washing as described for M. hapla.
Infection experiments were performed using a dialysis membrane assay . Briefly, conidia of A. oligospora, M. cionopagum and A. dactyloides were inoculated onto several pieces of sterilized dialysis membrane (spectra/por 4, Spectrumlabs). The membranes were placed over plates containing modified low-nutrient mineral salt (LNM) medium (KCl 1.0 g/l, MgSO4 0.2 g/l, ZnSO4.7H2O 0.88 mg/l, FeCl3.6H2O 3.0 mg/l, thiamine-HCl 0.2 mg/l, biotin 0.005 mg/l, L-phenylalanine-L-valine 0.1 g/l, agar 10 g/l, pH 6.5) [48, 49]. Infection structures (traps) were induced by adding 40–50 specimens of the nematode Panagrellus redivivus L. (Goodey) to the hyphae growing on each dialysis membrane. P. redivivus was grown axenically in a soya peptone-liver extract . After several days, when substantial amount of traps have been developed and all added nematodes have been killed and digested, the infection experiments were started by adding 75–100 surface sterilized second-stage juveniles of M. hapla or H. schachtii. The following five combinations of fungi and nematodes were examined: A. oligospora and M. hapla (designated Ao(Mh)), A. oligospora and H. schachtii (Ao(Hs)), A. dactyloides and M. hapla (Ad(Mh)), A. dactyloides and H. schachtii (Ad(Hs)), and M. cionopagum and H. schachtii (Mc(Hs)). The infection was followed under a light microscope, and the number of trapped, paralyzed (i.e. with arrested movements), and colonized (hyphae growing inside the capture nematode) were counted after various time periods. For each fungal and nematode interaction, 10 replications were used. Dialysis membranes having fungal and nematode interaction of each stage were quickly transferred into liquid nitrogen and ground. Materials were collected from all infection stages (trapped, paralyzed and infected (colonized) nematodes). The ground material was stored at -80°C until use.
RNA extraction, cDNA library construction and sequencing
Total RNA was extracted from each infection stage using the RNeasy Plant Mini kit and the RLC buffer (Qiagen) and subsequently quantified using a NanoDrop 2000C spectrophotometer (Thermo Scientific). RNA integrity was inspected using a RNA 6000 Pico kit on a 2100 BioAnalyzer (Agilent). Approximately equal amounts of RNA from the three infection stages of each fungus-nematode combination were pooled. The RNA pools were concentrated by precipitation using ammonium acetate/glycogen/ethanol as described in the MicroPoly(A)Purist Kit manual (Ambion).
From total RNA, mRNA was isolated using the PolyATtract kit (Promega) according to manufacturer’s instructions. Double-stranded cDNA was synthesized using the cDNA Synthesis System (Roche Diagnostics) according to the GS FLX Titanium cDNA Rapid Library preparation protocol (454/Roche) and using adaptors with Multiplex Identifiers (MIDs) that allow for pooling of the libraries prior sequencing. Library concentration was assessed by qPCR on a Mx3005P instrument (Stratagene) and using the KAPA Library Quantification Kit - 454 Titanium (Lib-L)/Universal (Kapa Biosystems). Based on the qPCR results all libraries were pooled to contain an equal molar amount of each library. Titration and library production (aiming at 7-16% enrichment) was performed using emulsion PCR and the Lib-L kit (454/Roche). DNA-containing beads were enriched and counted using a CASY Cell Counter DT (Roche Innovatis AG), processed using aXLR70 sequencing kit (454 Life Sciences/Roche Diagnostics), and loaded onto a picotiter plate for pyrosequencing on a 454 Life Sciences Genome Sequencer FLX machine (454 Life Sciences/Roche Diagnostics). Sequencing was conducted at the Lund University Sequencing Facility (Faculty of Science).
The reads obtained from the 454 sequencing were filtered, assembled and analysed according to the flowchart shown in Figure 2. Reads matching rRNA were removed using the BLASTN algorithm  with an E-value threshold of 1e-5 against a custom made database of rRNA sequences obtained from the 5S rRNA database  and the SILVA rRNA database . The remaining reads for each of the five libraries were assembled separately using the GS de novo assembler 2.6 (454 Life Sciences/Roche Diagnostics) with the -cdna option. The reads were assembled into 17 785 isotigs and 10 contigs with a length > 500 bp. In the following text, both categories were referred as isotigs, i.e. transcripts. Low abundance isotigs with less than five reads and isotigs with a length shorter than 100 bp were removed. Isotigs with top hits to non-fungal species in the UniProt database  (the BLASTX search)  were also removed. The filtered dataset contained 17 446 isotigs.
The filtered isotigs from the two A. oligospora samples were mapped to the A. oligospora genome using Gmap  to assess the quality of the de novo isotig assemblies. In total, 3 944 out of 3 952 isotigs (including 2 634 istotigs from the Ao(Mh) library and 1 318 isotigs from the Ao(Hs) library, Table 2)) matched the genome indicating efficient filtering of non-fungal transcripts. Only 75 of the 3 944 isotigs aligned to more than one position in the genome giving a total of 4 019 genome regions aligning to the isotigs. The low number of isotigs with multiple matches indicates a low frequency of chimeric transcripts. To further investigate the quality of the transcriptome assembly, we compared the genome sequences of isotig alignments with the 11 479 predicted genes from the A. oligospora genome using the Eval software . All of the 4 019 genome regions matched to the predicted A. oligospora genes. In total, 2 652 A. oligospora genes were matched giving on average 1.5 isotigs per predicted gene which suggest that some of the isotigs may represent alternative splicing forms. In total, 94.6 percent of the aligned genome regions (3 803 out of 4 019) contained both start and stop codon (i.e. considered complete by the Eval software). The mean length of the isotigs was 1 109 basepairs compared to 1 498 basepairs for the A. oligospora gene models. The difference in length may at least partly be due to alternative splicing forms where exon skipping will give shorter transcripts than the predicted gene models. The high proportion of successful matches to the A. oligospora genome as well as to its genes, the large number of complete isotigs and the long mean isotig length clearly indicate that most of the filtered isotigs have been correctly assembled into near full length transcripts.
The filtered isotigs were used to generate two different data sets (Figure 2). The first data set (“Highly expressed transcripts”) was normalized using two different approaches; reads per kilobase pair (kb), (the number of aligned reads per transcript was divided by the transcript length), and the reads per kilobase per million reads (RPKM) method  (Additional file 1). The isotigs were annotated based on homology using the BLASTX algorithm  (threshold values of 1e-10) to the UniProt sequence database  and proteins of M. haptotylum. The isotigs were also annotated using the pfam_scan.pl tool (ftp://ftp.sanger.ac.uk/pub/databases/Pfam/Tools/) to search the Pfam-A family protein database  with default thresholds. Secretion signals were predicted using the SignalP 4.0 algorithm . Isotigs were considered to encode putative secreted proteins if fulfilling at least one of the following three criteria: 1) Isotigs having a secretion signal in the same frame as the Pfam domain; 2) Isotigs having a secretion signal in the longest predicted open-reading frame (ORF) in the same frame as the BLASTX match (threshold value of 1e-10) to a protein in the UniProt database or protein from M. haptotylum; 3) Isotigs having a BLASTX match to a protein in the UniProt database or protein from M. haptotylum that contains a secretion signal. Orphans were identified as isotigs lacking both Pfam domains and BLASTX matches in the UniProt database and the M. haptotylum genome (threshold value of 1e-5) against a species other than itself. Orphans with a secretion signal in the longest ORF were considered to be putative secreted proteins. PCA and hierarchical clustering were performed using the Omics Explorer ver. 2.2 (Qlucore). Virulence-related genes were identified by BLASTX  similarity searches against the PHI-base database version 3.2  using a cutoff of < 1e-10.
The second data set (“Differentially expressed UniRef50 clusters”) was obtained by matching the isotig sequences to UniRef50 clusters . The UniRef50 clusters contain a representative of UniProt sequences that show 50% sequence similarity and 80% overlap with the longest sequence in the cluster. The isotigs with the highest BLASTX score to a UniRef50 sequence cluster from each library were considered as putative orthologs giving maximum one isotig from each species for any given UniRef50 cluster to take alternative splicing into account. Sometimes no isotig from one species has a significant match to a particular UniRef50 cluster (cutoff 1e-10). The transcript abundance of these putative orthologs, called “UniRef50 clusters” was normalized to correct for different library sizes using the R/Bioconductor software package DESeq 1.10.1  (Additional file 2).
To identify the third data set, “Host-specific gene expression” (Figure 2), GSMapper 2.8 (454 Life Sciences/Roche Diagnostics) was used with the -cdna and -cref parameters to map the reads from the two A. oligospora libraries (Ao(Mh) and Ao(Hs)) against the coding sequences of the 11 479 genes predicted in the A. oligospora genome . The read counts were normalized using DESeq . A homology search of the mapped A. oligospora proteins was carried out using the BLASTP algorithm  to the UniProt database and proteins of M. haptotylum (threshold values of 1e-5). Secretory A. oligospora proteins were predicted using SignalP 4.0 .
Sequence accession numbers
Sequences can be accessed from the database http://mbio-serv2.mbioekol.lu.se/NematodeTrappingFungi/. The short read pyrosequences from Ad(Mh), Ad(Hs), Ao(Mh), Ao(Hs) and Mc(Hs) are available at NCBI SRA database with the Bioproject IDs PRJNA230433, PRJNA230458, PRJNA230459, PRJNA230446 and PRJNA230448, respectively.
This work was supported by grants from the Swedish Research Council (VR). DK was supported by a DST (Department of Science and Technology, Gov. of India) BOYSCAST fellowship.
- Barron GL: The nematode-destroying fungi. 1977, Guelph: Canadian Biological PublicationsGoogle Scholar
- Tunlid A, Ahrén D: Molecular mechanisms of the interaction between nematode-trapping fungi and nematodes: Lessons from genomics. Biological control of plant-parasitic nematodes, vol 11. Edited by: Davies K, Spiegel Y. 2011, Springer: Netherlands, 145-169.View ArticleGoogle Scholar
- Dijksterhuis J, Veenhuis M, Harder W, Nordbring-Hertz B: Nematophagous fungi: physiological aspects and structure-function relationships. Adv Microb Physiol. 1994, 36: 111-143.PubMedView ArticleGoogle Scholar
- Nordbring-Hertz B, Jansson HB, Tunlid A: Nematophagous Fungi. 2011, John Wiley & Sons, Chichester: In Encyclopedia of Life SciencesView ArticleGoogle Scholar
- Liou GY, Tzean SS: Phylogeny of the genus Arthrobotrys and allied nematode-trapping fungi based on rDNA sequences. Mycologia. 1997, 89: 876-884. 10.2307/3761108.View ArticleGoogle Scholar
- Ahrén D, Ursing BM, Tunlid A: Phylogeny of nematode-trapping fungi based on 18S rDNA sequences. FEMS Microbiol Lett. 1998, 158: 179-184. 10.1016/S0378-1097(97)00519-3.PubMedView ArticleGoogle Scholar
- Li Y, Hyde KD, Jeewon R, Cai L, Vijaykrishna D, Zhang KQ: Phylogenetics and evolution of nematode-trapping fungi (Orbiliales) estimated from nuclear and protein coding genes. Mycologia. 2005, 97: 1034-1046. 10.3852/mycologia.97.5.1034.PubMedView ArticleGoogle Scholar
- Yang Y, Yang E, An ZQ, Liu XZ: Evolution of nematode-trapping cells of predatory fungi of the Orbiliaceae based on evidence from rRNA-encoding DNA and multiprotein sequences. Proc Natl Acad Sci U S A. 2007, 104: 8379-8384. 10.1073/pnas.0702770104.PubMed CentralPubMedView ArticleGoogle Scholar
- Higgins ML, Pramer D: Fungal Morphogenesis - Ring formation and closure by Arthrobotrys dactyloides. Science. 1967, 155: 345-346. 10.1126/science.155.3760.345.PubMedView ArticleGoogle Scholar
- Heintz CE, Pramer D: Ultrastructure of nematode-trapping fungi. J Bacteriol. 1972, 110: 1163-1170.PubMed CentralPubMedGoogle Scholar
- Liu KK, Tian JQ, Xiang MC, Liu XZ: How carnivorous fungi use three-celled constricting rings to trap nematodes. Protein Cell. 2012, 3: 325-328. 10.1007/s13238-012-2031-8.PubMedView ArticleGoogle Scholar
- Tunlid A, Jansson HB, Nordbring-Hertz B: Fungal attachment to nematodes. Mycol Res. 1992, 96: 401-412. 10.1016/S0953-7562(09)81082-4.View ArticleGoogle Scholar
- Nordbring-Hertz B, Jansson HB, Friman E, Persson Y, Dackman C, Trude H, Poloczek E, Feldman R: Nematophagous fungi. 1995, Göttingen, Germany: Institut für den Wissenschaftlichen Film, Film No. C 1851Google Scholar
- Jansson HB, Nordbring-Hertz B: Interactions between nematophagous fungi and plant-parasitic nematodes - Attraction, induction of trap formation and capture. Nematologica. 1980, 26: 383-389. 10.1163/187529280X00323.View ArticleGoogle Scholar
- Rosenzweig WD, Premachandran D, Pramer D: Role of trap lectins in the specificity of nematode capture by fungi. Can J Microbiol. 1985, 31: 693-695. 10.1139/m85-131.View ArticleGoogle Scholar
- Singh KP, Jaiswal RK, Kumar N, Kumar D: Nematophagous fungi associated with root galls of rice caused by Meloidogyne graminicola and its control by Arthrobotrys dactyloides and Dactylaria brochopaga. J Phytopathol. 2007, 155: 193-197. 10.1111/j.1439-0434.2007.01208.x.View ArticleGoogle Scholar
- Kumar D, Singh KP: Assessment of predacity and efficacy of Arthrobotrys dactyloides for biological control of root knot disease of tomato. J Phytopathol. 2006, 154: 1-5. 10.1111/j.1439-0434.2005.01047.x.View ArticleGoogle Scholar
- Jaffee BA, Muldoon AE: Susceptibility of root-knot and cyst nematodes to the nematode-trapping fungi Monacrosporium ellipsosporum and M. cionopagum. Soil Biol Biochem. 1995, 27: 1083-1090. 10.1016/0038-0717(94)00240-2.View ArticleGoogle Scholar
- Yang J, Wang L, Ji X, Feng Y, Li X, Zou C, Xu J, Ren Y, Mi Q, Wu J, Liu S, Liu Y, Huang X, Wang H, Niu X, Li J, Liang L, Luo Y, Ji K, Zhou W, Yu Z, Li G, Liu Y, Li L, Qiao M, Feng L, Zhang KQ: Genomic and proteomic analyses of the fungus Arthrobotrys oligospora provide insights into nematode-trap formation. PLoS Pathog. 2011, 7: e1002179-10.1371/journal.ppat.1002179.PubMed CentralPubMedView ArticleGoogle Scholar
- Meerupati T, Andersson KM, Friman E, Kumar D, Tunlid A, Ahrén D: Genomic mechanisms accounting for the adaption to parasitism in nematode-trapping fungi. PLoS Genet. 2013, 9: e1003909-10.1371/journal.pgen.1003909.PubMed CentralPubMedView ArticleGoogle Scholar
- Tautz D, Domazet-Loso T: The evolutionary origin of orphan genes. Nat Rev Genet. 2011, 12: 692-702.PubMedView ArticleGoogle Scholar
- Andersson KM, Meerupati T, Levander F, Friman E, Ahrén D, Tunlid A: Proteome of the nematode-trapping cells of the fungus Monacrosporium haptotylum. Appl Environ Microbiol. 2013, 79: 4993-5004. 10.1128/AEM.01390-13.PubMed CentralPubMedView ArticleGoogle Scholar
- Williamson VM, Gleason CA: Plant-nematode interactions. Curr Opin Plant. 2003, 6: 327-333. 10.1016/S1369-5266(03)00059-1.View ArticleGoogle Scholar
- Bird A: The development and organization of skeletal structures in nematodes. Croll NA: The organization of nematodes. 1976, London: Academic Press, 107-137.Google Scholar
- Camacho C, Coulouris G, Avagyan V, Ma N, Papadopoulos J, Bealer K, Madden TL: BLAST plus: architecture and applications. BMC Bioinformatics. 2009, 10: 421-10.1186/1471-2105-10-421.PubMed CentralPubMedView ArticleGoogle Scholar
- Suzek BE, Huang HZ, McGarvey P, Mazumder R, Wu CH: UniRef: comprehensive and non-redundant UniProt reference clusters. Bioinformatics. 2007, 23: 1282-1288. 10.1093/bioinformatics/btm098.PubMedView ArticleGoogle Scholar
- Anders S, Huber W: Differential expression analysis for sequence count data. Genome Biol. 2010, 11: R106-10.1186/gb-2010-11-10-r106.PubMed CentralPubMedView ArticleGoogle Scholar
- Kulkarni RD, Kelkar HS, Dean RA: An eight-cysteine-containing CFEM domain unique to a group of fungal membrane proteins. Trends Biochem Sci. 2003, 28: 118-121. 10.1016/S0968-0004(03)00025-2.PubMedView ArticleGoogle Scholar
- Nakatogawa H, Ichimura Y, Ohsumi Y: Atg8, a ubiquitin-like protein required for autophagosome formation, mediates membrane tethering and hemifusion. Cell. 2007, 130: 165-178. 10.1016/j.cell.2007.05.021.PubMedView ArticleGoogle Scholar
- Chen YL, Gao Y, Zhang KQ, Zou CG: Autophagy is required for trap formation in the nematode-trapping fungus Arthrobotrys oligospora. Environ Microbiol Rep. 2013, 5: 511-517. 10.1111/1758-2229.12054.PubMedView ArticleGoogle Scholar
- Endo Y, Tsurugi K: RNA N-Glycosidase activity of ricin A-chain - Mechanism of action of the toxic lectin ricin on eukaryotic ribosomes. J Biol Chem. 1987, 262: 8128-8130.PubMedGoogle Scholar
- Michiels K, Van Damme EJM, Smagghe G: Plant-insect interactions: What can we learn from plant lectins?. Arch Insect Biochem Physiol. 2010, 73: 193-212. 10.1002/arch.20351.PubMedView ArticleGoogle Scholar
- Wohlschlager T, Butschi A, Zurfluh K, Vonesch SC, Keller UAD, Gehrig P, Bleuler-Martinez S, Hengartner MO, Aebi M, Kunzler M: Nematotoxicity of Marasmius oreades agglutinin (MOA) depends on glycolipid binding and cysteine protease activity. J Biol Chem. 2011, 286: 30337-30343. 10.1074/jbc.M111.258202.PubMed CentralPubMedView ArticleGoogle Scholar
- Hamshou M, Smagghe G, Shahidi-Noghabi S, De Geyter E, Lannoo N, Van Damme EJM: Insecticidal properties of Sclerotinia sclerotiorum agglutinin and its interaction with insect tissues and cells. Insect Biochem Mol Biol. 2010, 40: 883-890. 10.1016/j.ibmb.2010.08.008.PubMedView ArticleGoogle Scholar
- Tunlid A, Jansson S: Proteases and their involvement in the infection and immobilization of nematodes by the nematophagous fungus Arthrobotrys oligospora. Appl Environ Microbiol. 1991, 57: 2868-2872.PubMed CentralPubMedGoogle Scholar
- Åhman J, Johansson T, Olsson M, Punt PJ, van den Hondel CA, Tunlid A: Improving the pathogenicity of a nematode-trapping fungus by genetic engineering of a subtilisin with nematotoxic activity. Appl Environ Microbiol. 2002, 68: 3408-3415. 10.1128/AEM.68.7.3408-3415.2002.PubMed CentralPubMedView ArticleGoogle Scholar
- Morano KA, Grant CM, Moye-Rowley WS: The response to heat shock and oxidative stress in Saccharomyces cerevisiae. Genetics. 2012, 190: 1157-1195. 10.1534/genetics.111.128033.PubMed CentralPubMedView ArticleGoogle Scholar
- Heller J, Tudzynski P: Reactive oxygen species in phytopathogenic fungi: Signaling, development, and disease. Annu Rev Phytopathol. 2011, 49: 369-390. 10.1146/annurev-phyto-072910-095355.PubMedView ArticleGoogle Scholar
- Egan MJ, Wang ZY, Jones MA, Smirnoff N, Talbot NJ: Generation of reactive oxygen species by fungal NADPH oxidases is required for rice blast disease. Proc Natl Acad Sci U S A. 2007, 104: 11772-11777. 10.1073/pnas.0700574104.PubMed CentralPubMedView ArticleGoogle Scholar
- Xue CY, Park G, Choi WB, Zheng L, Dean RA, Xu JR: Two novel fungal virulence genes specifically expressed in appressoria of the rice blast fungus. Plant Cell. 2002, 14: 2107-2119. 10.1105/tpc.003426.PubMed CentralPubMedView ArticleGoogle Scholar
- Rosén S, Ek B, Rask L, Tunlid A: Purification and characterization of a surface lectin from the nematode-trapping fungus Arthrobotrys oligospora. J Gen Microbiol. 1992, 138: 2663-2672. 10.1099/00221287-138-12-2663.PubMedView ArticleGoogle Scholar
- Rosén S, Sjollema K, Veenhuis M, Tunlid A: A cytoplasmic lectin produced by the fungus Arthrobotrys oligospora functions as a storage protein during saprophytic and parasitic growth. Microbiology. 1997, 143: 2593-2604. 10.1099/00221287-143-8-2593.View ArticleGoogle Scholar
- Balogh J, Tunlid A, Rosén S: Deletion of a lectin gene does not affect the phenotype of the nematode-trapping fungus Arthrobotrys oligospora. Fungal Genet Biol. 2003, 39: 128-135. 10.1016/S1087-1845(03)00023-9.PubMedView ArticleGoogle Scholar
- Winnenburg R, Urban M, Beacham A, Baldwin TK, Holland S, Lindeberg M, Hansen H, Rawlings C, Hammond-Kosack KE, Kohler J: PHI-base update: additions to the pathogen-host interaction database. Nucleic Acids Res. 2008, 36: 572-576.View ArticleGoogle Scholar
- Braun BR, Head WS, Wang MX, Johnson AD: Identification and characterization of TUP1-regulated genes in Candida albicans. Genetics. 2000, 156: 31-44.PubMed CentralPubMedGoogle Scholar
- Tzelepis GD, Melin P, Jensen DF, Stenlid J, Karlsson M: Functional analysis of glycoside hydrolase family 18 and 20 genes in Neurospora crassa. Fungal Genet Biol. 2012, 49: 717-730. 10.1016/j.fgb.2012.06.013.PubMedView ArticleGoogle Scholar
- Bairoch A, Apweiler R, Wu CH, Barker WC, Boeckmann B, Ferro S, Gasteiger E, Huang H, Lopez R, Magrane M, Martin MJ, Natale DA, O'Donovan C, Redaschi N, Yeh LS: UniProt: the Universal Protein knowledgebase. Nucleic Acids Res. 2004, 32: 115-119. 10.1093/nar/gkh151.View ArticleGoogle Scholar
- Nordbring-Hertz B: Dialysis membrane technique for studying microbial interaction. Appl Environ Microbiol. 1983, 45: 290-293.PubMed CentralPubMedGoogle Scholar
- Nordbring-Hertz B: Peptide-induced morphogenesis in the nematode-trapping fungus Arthrobotrys oligospora. Physiol Plant. 1973, 29: 223-233. 10.1111/j.1399-3054.1973.tb03097.x.View ArticleGoogle Scholar
- Nordbring-Hertz B: Nematode-induced morphogenesis in predacious fungus Arthrobotrys oligospora. Nematologica. 1977, 23: 443-451. 10.1163/187529277X00372.View ArticleGoogle Scholar
- Szymanski M, Barciszewska MZ, Erdmann VA, Barciszewski J: 5S ribosomal RNA database. Nucleic Acids Res. 2002, 30: 176-178. 10.1093/nar/30.1.176.PubMed CentralPubMedView ArticleGoogle Scholar
- Quast C, Pruesse E, Yilmaz P, Gerken J, Schweer T, Yarza P, Peplies J, Glockner FO: The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 2013, 41: D590-D596. 10.1093/nar/gks1219.PubMed CentralPubMedView ArticleGoogle Scholar
- Wu TD, Watanabe CK: GMAP: a genomic mapping and alignment program for mRNA and EST sequences. Bioinformatics. 2005, 21: 1859-1875. 10.1093/bioinformatics/bti310.PubMedView ArticleGoogle Scholar
- Keibler E, Brent M: Eval: a software package for analysis of genome annotations. BMC Bioinformatics. 2003, 4: 50-10.1186/1471-2105-4-50.PubMed CentralPubMedView ArticleGoogle Scholar
- Mortazavi A, Williams BA, Mccue K, Schaeffer L, Wold B: Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008, 5: 621-628. 10.1038/nmeth.1226.PubMedView ArticleGoogle Scholar
- Punta M, Coggill PC, Eberhardt RY, Mistry J, Tate J, Boursnell C, Pang N, Forslund K, Ceric G, Clements J, Heger A, Holm L, Sonnhammer EL, Eddy SR, Bateman A, Finn RD: The Pfam protein families database. Nucleic Acids Res. 2012, 40: D290-D301. 10.1093/nar/gkr1065.PubMed CentralPubMedView ArticleGoogle Scholar
- Petersen TN, Brunak S, von Heijne G, Nielsen H: SignalP 4.0: discriminating signal peptides from transmembrane regions. Nat Methods. 2011, 8: 785-786. 10.1038/nmeth.1701.PubMedView 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/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.