- Research article
- Open Access
De novo genome assembly of the soil-borne fungus and tomato pathogen Pyrenochaeta lycopersici
© Aragona et al.; licensee BioMed Central Ltd. 2014
- Received: 20 June 2013
- Accepted: 22 April 2014
- Published: 27 April 2014
Pyrenochaeta lycopersici is a soil-dwelling ascomycete pathogen that causes corky root rot disease in tomato (Solanum lycopersicum) and other Solanaceous crops, reducing fruit yields by up to 75%. Fungal pathogens that infect roots receive less attention than those infecting the aerial parts of crops despite their significant impact on plant growth and fruit production.
We assembled a 54.9Mb P. lycopersici draft genome sequence based on Illumina short reads, and annotated approximately 17,000 genes. The P. lycopersici genome is closely related to hemibiotrophs and necrotrophs, in agreement with the phenotypic characteristics of the fungus and its lifestyle. Several gene families related to host–pathogen interactions are strongly represented, including those responsible for nutrient absorption, the detoxification of fungicides and plant cell wall degradation, the latter confirming that much of the genome is devoted to the pathogenic activity of the fungus. We did not find a MAT gene, which is consistent with the classification of P. lycopersici as an imperfect fungus, but we observed a significant expansion of the gene families associated with heterokaryon incompatibility (HI).
The P. lycopersici draft genome sequence provided insight into the molecular and genetic basis of the fungal lifestyle, characterizing previously unknown pathogenic behaviors and defining strategies that allow this asexual fungus to increase genetic diversity and to acquire new pathogenic traits.
- Pyrenochaeta lycopersici
- Genome assembly
- Next Generation Sequencing technologies (NGS)
Genome sequencing and assembly
Genomic DNA obtained from the virulent P. lycopersici isolate CRA-PAV_ER 1211 was sequenced using Illumina 100-bp paired-end reads. Two libraries were prepared with mean fragment sizes of 460 and 560 bp respectively and sequenced obtaining 85 millions of fragments for a total of 17Gb, corresponding to approximately 300-fold coverage of the final assembly.
The P. lycopersici genome was assembled de novo using a de Bruijn graph-based (DBG) assembly strategy with a k-mer of 55 (Additional file 1: Figure S1) before reassembly using the overlap-layout-consensus (OLC) method. We chose to use this strategy to take advantage of both the higher efficiency of DBG method in assembling millions of reads and the higher sensitivity of OLC method [22, 23]. In fact, the OLC re-assembly of sequences previously assembled with DBG method reduced the total number of sequences by a 40% (from 11,617 to 7,079) while increasing the average contig length by a 60% (from 4,834 to 7,747 bp) and not reducing significantly the total number of assembled bases (54.8 Mbps vs. 56.2 Mbps).
Pyrenochaeta lycopersici genome statistics
Number of sequences
Maximum sequence length (bp)
Average length (bp)
Sequences > 500bp
Number of sequences
Average length (bp)
Sequences > 1Kb
Number of sequences
Average length (bp)
Sequences > 5Kb
Number of sequences
Average length (bp)
Sequences > 10Kb
Number of sequences
Average length (bp)
Total number of assembled bases (without Ns)
54,841,665 bp (53,572,698 bp)
Estimated x-fold genome coverage
GC percentage content
Number of annotated gene loci
Number of annotated putative transcripts
Mean gene size (bp)
Mean exon size (bp)
Mean number of exons per gene
Genes were annotated by a combination of ab initio prediction  followed by the reference annotation-based transcript (RABT) assembly of sequences obtained from the RNA-Seq analysis of two P. lycopersici samples grown in vitro and during interaction with tomato cv Moneymaker, respectively . We annotated 17,411 genes with 27,275 transcripts, among which 9,553 loci were detected by both approaches and 2,797 only by the analysis of RNA-Seq data (Table 1). Thus, while the number of genes might be slightly overestimated because of ab initio prediction limitations, at least 70.8% of the annotations were supported by experimental evidences. We confirmed that the assembly was able to capture full-length genes by searching the predictions for full open reading frames (ORFs), finding that most of the genes (80.76%) contained start and stop codons. The RNA-Seq data were then assembled de novo and the resulting contigs were mapped onto the draft genome. From the 31,746,550 reads, we obtained 27,982 putative transcripts, 27,574 (98.5%) of which could be mapped onto the assembled genome further validating the comprehensiveness of the gene space represented.
Many of the identified transcripts (75.8%) were conserved in other species as shown by hits against sequences in the NCBI NR protein database (e-value < 1E-0.6) and the Uniprot SwissProt Fungi protein database (e-value < 1E-0.7). This analysis showed that the most represented species were fungal pathogens and the top four species were phytopathogenic fungi (Additional file 1: Figure S2). Based on the identity of the most similar proteins, at least one Gene Ontology term was assigned to 8,507 transcripts. The most represented functional categories are summarized in Additional file 1: Figure S3. We found 3,962 orphan genes (22.8%) with no matches against known proteins or protein domains , which is similar to the proportion of orphan genes found in Sordaria macrospora (22%)  and Macrophomina phaseolina (29%) . RNA-Seq analysis showed that 1,936 of the 3,962 orphan genes were transcribed and were therefore likely to be functional.
We performed a pairwise comparison between P. lycopersici and Fusarium oxysporum f. sp. lycopersici genomes to assess the genomic distribution of similarity regions. We choose F. oxysporum among the top ranked species (Additional file 2: Table S2) as a complete genome anchored to chromosomes was available for this fungus. As expected the analysis showed regions of similarity at aminoacidic level throughout all the F. oxysporum chromosomes with the exception of the four F. oxysporum dispensable chromosomes 3, 6, 14 and 15 (Additional file 1: Figure S4).
Pathogenesis related genes
Highly represented protein families
Genome assembly size ( Mb)
Number of genes
DDE superfamily endonuclease
Double stranded RNA binding domain
NAD-specific glutamate dehydrogenase
Transmembrane alpha-helix domain
hAT family dimerisation domain
Helix-turn-helix, Psq domain
Major facilitator superfamily
Fungal specific transcription factor domain
Protein kinase domain
Heterokaryon incompatibility protein (HET)
The ATP-binding cassette (ABC) transporters and MFS-type membrane transporters were largely represented in the P. lycopersici genome, comprising 109 and 229 modules respectively (Figure 3b and Additional file 2: Table S2). Transcriptome analysis showed that 85% of ABC modules and 75% of MFS-type membrane transporters belonged to transcripts expressed. The P. lycopersici genome also encoded a large number (125) of cytochrome P450 proteins, as shown for other necrotrophic and hemibiotrophic pathogens such as B. cinerea (120), C. graminicola (138), F. oxysporum (157) and S. nodorum (122).
The P. lycopersici predicted genes also included 597 sequences matching 94 different subfamilies of peptidases (Additional file 2: Table S3). Almost all known peptidase families were represented in the P. lycopersici transcriptome, with DmpA aminopeptidase 1 (P1) as the only exception. The most represented clans were the serine peptidase (S clan) and metallopeptidase (M clan), a common feature of fungal pathogens. The comparative analysis of gene families and PFAM domains with several other fungi showed that the number of peptidases in P. lycopersici is comparable to other hemibiotrophic pathogens as C. graminicola and F. oxysporum. Metallopeptidase group of aminopeptidases (M1, M18, M24, M28) were highly represented, with 43 domains in P. lycopersici (shown in yellow in Additional file 2: Table S3) compared to 34 and 30 of C. graminicola and F. oxysporum respectively. Analysis of expression data confirmed that 36 of this modules (85%) corresponded to expressed transcripts. Other metallopeptidase gene families had also undergone expansion in the P. lycopersici genome, including pappalysin, Ste24 and deubiquitinating peptidase (shown in green in Additional file 2: Table S3). A summary is reported in Figure 3d.
Genes involved in carbohydrate degradation (CAZymes)
Carbohydrate degradation is an important component of fungal pathogenicity and virulence, so we examined the P. lycopersici CAZome in detail compared with nine other fungi with complete genome sequences (Additional file 2: Table S6). The P. lycopersici genome encodes 575 CAZyme modules, including 272 glycoside hydrolases, 83 glycosyltransferases, 20 polysaccharide lyases, 149 carbohydrate esterases and 51 carbohydrate-binding modules (CMBs). CE1 and CE10 families of carbohydrate esterases including acetyl xylan esterase (EC 126.96.36.199), cinnamoyl esterase (EC 3.1.1), feruloyl esterase (EC 188.8.131.52), carboxylesterase (EC 184.108.40.206), S-formylglutathione hydrolase (EC 220.127.116.11) and sterol esterases are particularly represented (53 CE1 modules out of which 77% are detected as expressed by RNASeq analysis, 54 modules for CE10 family 73,5% of which are expressed) similarly to other hemibiotrophs such as C. graminicola and F. oxysporum (Additional file 2: Table S6). Cellulose-degrading enzymes were also well represented in P. lycopersici, including seven cellobiohydrolases (GH6), ten β-1,3-glucanases (GH55) and eight GH105 glycoside hydrolases (unknown mechanism). Finally, the complex of carbohydrate cleaving enzymes was completed by 20 genes encoding polysaccharide lyases, including 13 encoding PL3.
We report here the first genome analysis of a Pyrenochaeta species, the soil-borne filamentous fungus of ascomycete clade Pyrenochaeta lycopersici, which is the etiological agent of tomato corky root rot. The P. lycopersici genome assembly shows that a genome reconstruction approach based uniquely on paired-end Illumina reads is highly effective in reconstructing contigs containing almost full length genes. To validate the genome assembly we assessed the completeness of the gene space represented and we found that, in fact, most of the core eukaryotic conserved genes and the transcripts reconstructed from RNA-Seq data of fungus grown in vitro and during interaction with the host were represented on the assembled genome. Based on published assemblies for the phytopathogens L. maculans, P. teres, P. nodorum and P. tritici-repentis the number of annotated genes in P. lycopersici is very similar (17,411 versus 12,469, 11,799, 12,382 and 12,300 respectively) [31–34]. Most of the gene annotations are supported by the presence of a full ORF (about 80%) and a large fraction of them is validated by RNA-Seq data (70%), thus representing an high quality resource for the study of functions encoded by P. lycopersici genome. It’s worth noting that we also identified more than 2,700 genes supported by transcriptomic data but not detected by the prediction algorithm. This is not completely surprising as, even if properly trained, a prediction software will not predict all the genes of an eukaryote organism. An hybrid annotation approach, based on integration of gene predictions and massive parallel sequencing of transcripts, is thus needed to perform a comprehensive annotation of genes in a fungal genome. Overall these data confirm the high quality of the assembly and annotation obtained particularly in terms of the completeness and quality of the gene catalog represented.
The P. lycopersici assembly produced constitutes an invaluable support to understand the unique phenotypical features of this pathogen and allow to investigate the molecular basis of the reproductive behavior and of the mechanisms involved in the pathogenesis. Sexual mating of P. lycopersici has never been observed in nature, leading to speculation regarding its reproductive cycle. In agreement with its reproductive behavior based on generating spores (conidia) by mitosis, P. lycopersici apparently lacks of mating-type (MAT) genes that control the choice between sexual and asexual reproduction, suggesting the species is incapable of sexual reproduction. A potential alternative source of genetic variation in P. lycopersici is vegetative hyphal fusion controlled by HET genes, allowing horizontal gene or chromosome transfer potentially followed by non-meiotic recombination. When individuals with the same HET genotype meet, they can produce a viable heterokaryon by anastomosis, whereas individuals with different HET genotypes form a fusion cell which is compartmentalized and undergoes a form of programmed cell death termed vegetative or heterokaryon incompatibility . Although the significance of heterokaryon incompatibility responses is poorly understood , it may limit the transmission of mycoviruses and other deleterious replicons between strains . Strong selective pressure is likely to be responsible for the extensive amplification of HET genes and genes with related functions in P. lycopersici, suggesting that HET genes play a key role in the transfer of genetic information in this species, and that a strict regulation of the vegetative reproduction is important for the generation of the variability necessary for the adaptation to the environment and to host defense mechanisms.
From a pathogenetic perspective P. lycopersici is considered a hemibiotrophic fungus because cell wall hydrolytic activity is not detected during the initial infection of tomato plants and few symptoms are visible, whereas later infection involves the secretion of cell wall degrading enzymes that cause the root to collapse [15, 17]. This was confirmed by our data by the comparison of P. lycopersici sequences with those of other fungi which showed a clear phylogenetic relationship with hemibiotrophic and necrotrophic plant pathogens. Moreover, the analysis of gene functions and comparison of gene sequences with those of other fungi showed that a large fraction of genome is devoted to pathogenetic activity and showed a large overlap of the gene inventory to that of other plant pathogens such as L. maculans, P. teres, P. nodorum and P. tritici-repentis.
The first and major barrier to infection by fungal pathogens in plants is the cell wall and cellulose is the main component of plant biomass. Phytopathogenic fungi therefore secrete a cocktail of hydrolytic enzymes known as carbohydrate-active enzymes (CAZymes), which are required to penetrate and then degrade the cuticle and cell wall . The analysis of CAZome of P. lycopersici showed a stronger resemblance to that of hemibiotrophic and necrotrophic plant pathogens such as C. graminicola, L. maculans and P. tritici-repentis than to that of biotrophs such as B. graminis. In fact, the P. lycopersici genome encoded a large number of cellulose-degrading enzymes, as glycoside hydrolases, required for the complete breakdown of the plant cell wall for successful infection. In particular, the CE1 and CE10 families of carbohydrate esterases, which are required to degrade hemicellulose and thus facilitate the complete hydrolysis of polysaccharides in the cell walls of a wide range of plant species, as well as pectate-lyases (which doubles the number of PLs in P. lycopersici compared to other fungi) were largely represented. This may reflect in an enhanced activity in cleaving pectic polymers, which are more abundant in dicotyledonous plants. Among the glycoside hydrolases, the GH61 family carries out cellulose hydrolysis using a synergic mechanism in concert with canonical cellulases [39, 40]. The expansion of the GH61 family in P. lycopersici mirrors the expansion of this family in the order Pleosporales. The GH61 gene family is exclusive to fungi and structure–function analysis of some enzymes belonging to this class, showed they cleave cellulose using an oxidative mechanism, and therefore they are not canonical glycoside hydrolases. The new biochemical mechanism proposed for GH61 proteins redefined this class of enzymes as polysaccharide monooxygenases (PMOs) [42, 43]. A P. lycopersici endo-β-1,4-glucanase gene (Plegl1) from the GH61 family is expressed in a manner that corresponds to disease progression [17, 44]. Plegl1 is thus far the only gene known to be induced in a fungal phytopathogen during infection, suggesting a role in pathogenesis. These data suggest the hypothesis that P. lycopersici might have evolved diverse strategies for cellulose digestion.
Other proteins which play an important role in the interaction with the host by detoxifying plant defense compounds are cytochrome P450 proteins. P. lycopersici, similarly to other hemibiotroph and necrotroph fungal pathogens, shows an high number and an high variety of genes coding for this protein family.
Peptidases are pathogenicity-related enzymes that are secreted to facilitate penetration and colonization of the host by degrading the plant cell wall [45–48] and are required for the degradation of plant defense proteins [49–51]. The P. lycopersici genome encodes a large number of metallopeptidases, covering 10 major peptidase families and 94 subfamilies, in agreement with the general properties of Dothideomycetes genomes, mostly representing hemibiotrophs . The comparison of 18 genomes revealed that Dothideomycetes have a wider range of exopeptidases and endopeptidases than other fungal phytopathogens, including the greatest number of secreted metallopeptidases, but fewer aspartic peptidases (A01) than necrotrophs, saprotrophs and ectomycorrhizal symbionts . The aminopeptidase gene family is the most highly represented among all the metallopeptidase families represented in the P. lycopersici genome and it is well documented that bacteria and fungi which secrete aminoproteases are generally pathogenic [52, 53].
Another class of proteins which play an important role during host invasion is that of transporters. Transporters import nutrients and export secondary metabolites produced by the fungal pathogens as virulence factors but have a role also in removing toxic compounds. In particular, ABC and MFS transporters are the major families involved [54–56] as they are required to export host-specific toxins (HSTs) and mycotoxins [57–59], remove inhibitory defense compounds such as phytoalexins produced by the host plant , and confer resistance to fungicides . Although the virulence of the fungus is not strictly dependent on the abundance of these transporter families  is notable that the ABC transporter family in P. lycopersici, together with F. oxysporum, is the most diverse and abundant compared to all the other fungi, and is therefore likely to be intimately associated with the virulence mechanism. Overall the analysis of P. lycopersici genome clearly shows that this pathogen evolved its pathogenetic mechanisms through an expansion both of genes involved in the penetration and degradation of the host tissues and by the expansions of gene families necessary to counteract the defense mechanisms of the host.
P. lycopersici genome reveals a significative expansion of specific genes families related both to pathogenesis and to reproduction mechanisms, which suggests that P. lycopersici has undergone to a specialization and adaptation process during its evolution. The assembly presented constitutes an important resource to understand the molecular bases of corky root rot and more in general to enrich current knowledge of plant-pathogen interaction mechanisms.
Sample and library preparation
Genomic DNA was isolated from a virulent P. lycopersici isolate (CRA-PAV_ER 1211) grown on Potato Dextrose Agar (PDA) medium, by the method described by Cenis as previously described  with the following modifications: 200 mg of mycelium was frozen in liquid nitrogen, pulverized, and incubated in 300 μl of lysis buffer (200 mM Tris-HCl pH 8.5, 250 mM NaCl, 25 mM EDTA, 0.5% SDS) for 10 min at 65°C. We then added 150 μl 3 M sodium acetate (pH 5. 2), incubated at –20°C for 10 min and centrifuged at 10000 × g for 30 min. The supernatant was transferred in a fresh tube and the DNA was precipitated by adding an equal volume of isopropanol, centrifuging at 10000 × g for 10 min and washing with 70% ethanol. Finally the DNA was purified using the NucleoSpin Extract II kit (Macherey-Nagel, Düren, Germany). Total RNA from vegetative mycelium grown on PDA medium was extracted using the RNeasy Midi kit (Qiagen, Hilden, Germany). Total RNA from infected tomato roots of cv Moneymaker was extracted using the NucleoSpin RNA Plant 2 kit (Macherey-Nagel) after 8 days post infection (dpi).
Genomic DNA (6 μg) was fragmented by nebulization at 35 psi for 6 min. DNA libraries with insert sizes of 400 bp and 560 bp were prepared from 1 μg of fragmented genomic DNA using the paired-end TruSeq DNA Sample Preparation Kit (Illumina Inc., San Diego, CA, USA). The library quality was determined using the High Sensitivity DNA Kit (Agilent, Wokingham, UK).
Total RNA samples were assessed for quality using an RNA 6000 Nano Kit (Agilent) and 2.5-μg aliquots were used to isolate poly(A) mRNA for the preparation of a non-directional Illumina RNA-Seq library using the TruSeq RNA Sample Prep Kit (Illumina). The quality of the library was checked using the High Sensitivity DNA Kit (Agilent).
Sequencing and data preprocessing
Libraries were sequenced with an Illumina GAIIx sequencer generating 100-bp paired-end sequences for DNA libraries and 130-bp paired-end sequences for RNA libraries.
The sequences were pre-processed by removing reads with a number of N >10 or with a read quality <20 using a custom script. Adapters were clipped using Scythe v0.980  and bases on both 3′ ends with a quality <20 were trimmed using Sickle v0.940 , eventually entirely removing the fragment if the length was reduced to < 50 bp.
De novo assembly and gene catalog assessment
The genome was assembled de novo using Velvet v1.1.06  with the following parameters: –exp_cov auto (automatic calculation of expected coverage), –scaffolding (scaffolding of contigs with paired-end reads) and –min_contig_lgth 200 (mimimun contig length = 200). The optimal k-mer length was determined by adjusting the k-mer length from 39 to 67 bp in 4-bp increments and using the k-mer for which the N50 and the maximum contig length reached the highest value (Additional file 1: Figure S1). The resulting contigs were then re-assembled with CAP3 v10/15/07  using standard parameters.
Core Eukaryotic Genes (CEGs) were aligned with assembled genome using BLAST and hits were considered significant when the sequence identity was >65%.
RNA-seq reads from in vitro mycelia were assembled using Trinity v r2011-11-26  with standard parameters, jaccard clip on and a minimum contig length of 200 bp. Assembled contigs were mapped using GMAP  with standard parameters.
The final assembly was processed by GeneMark.hmm-ES v2.3e  with standard parameters and no a priori information. The genome was masked for repetitive and low-complexity regions with RepeatMasker v open-3.3.0  with standard parameters and a general repeats database. The resulting annotation was refined using TopHat v1.4.1 and Cufflinks v1.2.1  on the two RNA-seq libraries in RABT mode with standard parameters. ORFs were identified for each transcript using CPC v0.9.r2 .
A phylogenetic tree based on the comparison of whole genomes of P. lycopersici and 16 other fungi was constructed using CVTree v2  at a k-mer length of 7. Rhizopus oryzae was used as the outgroup for building the unscaled tree . The P. tritici-repentis and C. graminicola proteomes were obtained from the Colletotrichum Sequencing Project , the L. maculans proteome was obtained from the L. maculans genome project , and the B. graminis proteome was obtained from the Blumeria sequencing project . Unless otherwise stated, the remaining proteomes were obtained from the CVtree inbuilt genome database .
Functional annotation was initiated by using each sequence as a BLAST query  against the NCBI Non Redundant database retrieved 2012-09-14  (e-value < 1E-0.6) and the Uniprot SwissProt Fungi protein database retrieved 2012-04-30  (e-value < 1E-0.7). The results were analyzed using Blast2GO  and integrated with InterPro results . Conserved protein domains in P. lycopersici and all the other fungi considered (AN, BC, BG, CG, FO, LM, NC, PT-R, SN) were identified using HMMer v3.0  to identify homology with proteins in the Pfam-A database (v26.0, 2011-11) . Sequence conservation was considered significant at an e-value threshold <1e-6 for both the entire sequence match and for the independent E-value of the single domain match. CAZymes (v2.0)  homology was also inspected using HMMer. Alignments were considered significant at an alignment length > 80 residues, E-value < 1e-5 and HMM profile coverage > 30% or alignment length < 80 residues and E-value < 1e-3 and HMM profile coverage > 30%. BLASTX was used to identify sequences homologous to known pathogenic genes (PHI-base ver 3.2) , peptidases (MEROPS ver 9.8) , zinc fingers (C2H2 ZNF db, ver. 2007-10-03) , MAP kinase sequences from NCBI NR (retrieved 2013-01-30)  and membrane transport proteins (TCDB, ver. 2011-July-15)  (e-value <1e-10). Sequences with significant hits against membrane transport proteins were also used as BLASTX queries against a G protein-coupled receptors database (GPCRDB, retrieved 2013-01-30)  and fungi major facilitator superfamily sequences from NCBI NR .
Significance of protein families abundance differences between P. lycopersici and all the other fungal plant pathogens (BC, BG, CG, FO, LM, PT-R, SN) was assessed by a 1-sample t-test. Variance was estimated based on protein families abundance data of all the fungal pathogens taken into account (P. lycopersici excluded). P-values were corrected according to Benjamini and Hochberg  on the full dataset of comparisons.
Genome coverage has been estimated by mapping the reads on the assembled genome using BWA v. 0.6.2-r126  using default parameters and calculating coverage on a panel of 140 single copy genes.
RNASeq reads and transcriptome assemblies have been deposited at the NCBI Sequence Reads Archive (SRA) and NCBI Transcriptome shotgun assembly (TSA) databases respectively and are available under BioProject number PRJNA202292. Genomic reads and genome assembly have been deposited at the NCBI Sequence Reads Archive (SRA) and are available under BioProject number PRJNA202288. Assemblies have been deposited as Whole Genome Shotgun project at DDBJ/EMBL/GenBank under the accession ASRS00000000. The version described in this paper is version ASRS01000000.
This research was supported by the Italian national project ‘Identificazione di geni implicati nella resistenza e nella patogenicità in interazioni tra piante di interesse agrario e patogeni fungini, batterici e virali’ (‘RESPAT’) funded by MiPAAF and by Fondazione Cariverona (Completamento e attività del Centro di Genomica Funzionale Vegetale), Verona, Italy.
- Termohlen GP: On corky root of tomato and the corky root fungus. Tijdschr Plantenziekten. 1962, 68: 295-367.Google Scholar
- Gerlach W, Schneider R: Nachweis eines Pyrenochaeta Stadiums bei Stammen des Korkwurzelerregers der Tomate. Phytopath Z. 1964, 50: 262-269. 10.1111/j.1439-0434.1964.tb02924.x.View ArticleGoogle Scholar
- Grove GG, Campbell RN: Host range and survival in soil of Pyrenochaeta lycopersici. Plant Dis. 1987, 71: 806-809. 10.1094/PD-71-0806.View ArticleGoogle Scholar
- Infantino A, Di Giambattista G, Porta-Puglia A: First report of Pyrenochaeta lycopersici on melon in Italy. Petria. 2000, 10: 195-198.Google Scholar
- Pohronezny KL, Volin RB: Corky Root Rot. Compendium of Tomato Diseases. Edited by: Jones JB, Jones JP, Stall RE, Zitter TA. 1991, Minnesota: The American Phytopathological Society, 12-13.Google Scholar
- Aragona M, Infantino A, Papacchini M: Developing a molecular method for screening the resistance to a pathogen of tomato to contribute to limit the use of toxic chemicals in soil. WIT Trans Ecol Envir. 2009, 120: 519-524.View ArticleGoogle Scholar
- Campbell RN, Hall DH, Schweers VH: Corky root of tomato in California caused by Pyrenochaeta lycopersici and control by soil fumigation. Plant Dis. 1982, 66: 657-661. 10.1094/PD-66-657.View ArticleGoogle Scholar
- Ekengren SK: Cutting the Gordian knot: taking a stab at corky root rot of tomato. Plant Biotechnol (Tsukuba). 2008, 25: 265-10.5511/plantbiotechnology.25.265.View ArticleGoogle Scholar
- de Gruyter J, Woudenberg JHC, Aveskamp MM, Verkley GJM, Groenewald JZ, Crous PW: Systematic reappraisal of species in Phoma section Paraphoma. Pyrenochaeta Pleurophoma Mycologia. 2010, 102 (5): 1066-1108. 10.3852/09-240.View ArticleGoogle Scholar
- Infantino A, Aragona M, Brunetti A, Lahoz E, Oliva A, Porta Puglia A: Molecular and physiological characterization of Italian isolates of Pyrenochaeta lycopersici. Mycol Res. 2003, 107 (6): 707-716. 10.1017/S0953756203007962.PubMedView ArticleGoogle Scholar
- Bayraktar H, Oksal E: Molecular, physiological and pathogenic variability of Pyrenochaeta lycopersici associated with corky rot disease of tomato plants in Turkey. Phytoparasitica. 2011, 39: 165-174. 10.1007/s12600-011-0150-z.View ArticleGoogle Scholar
- Pucci N, Ferrante M, Infantino A: Study of genetic structure of Italian populations of Pyrenochaeta lycopersici by AFLP analysis. Acta Hortic. 2011, 914: 121-124.View ArticleGoogle Scholar
- White JG, Scott AC: Formation and ultrastructure of microsclerotia of Pyrenochaeta lycopersici. Ann Appl Biol. 1973, 73: 163-166. 10.1111/j.1744-7348.1973.tb01321.x.View ArticleGoogle Scholar
- Ball SFL: Morphogenesis and structure of microsclerotia of Pyrenochaeta lycopersici. T Brit Mycol Soc. 1979, 73: 366-368. 10.1016/S0007-1536(79)80129-1.View ArticleGoogle Scholar
- Goodenough PW, Kempton RJ: The activity of cell wall degrading enzymes in tomato roots infected with Pyrenochaeta lycopersici and the effect of sugar concentrations in these roots on disease development. Physiol Plant Pathol. 1976, 9: 313-320. 10.1016/0048-4059(76)90064-3.View ArticleGoogle Scholar
- Goodenough PW, Kempton RJ, Maw GA: Studies on the root rotting fungus Pyrenochaeta lycopersici: extracellular enzyme secretion by the fungus grown on cell wall material from susceptible and tolerant tomato plants. Physiol Plant Pathol. 1976, 8: 243-251. 10.1016/0048-4059(76)90019-9.View ArticleGoogle Scholar
- Valente MT, Infantino A, Aragona M: Molecular and functional characterization of an endoglucanase in the phytopathogenic fungus Pyrenochaeta lycopersici. Curr Genet. 2011, 57: 241-251. 10.1007/s00294-011-0343-5.PubMedView ArticleGoogle Scholar
- Shishkoff N: Pyrenochaeta. Methods for Research in Soilborne Phytopathogenic Fungi. Edited by: Singleton L, Mihail JD, Ryush CM. 1992, St Paul, Minnesota: APS press, 153-156.Google Scholar
- Aragona M, Infantino A: Expression profiling of tomato response to Pyrenochaeta lycopersici infection. Acta Hortic. 2008, 789: 257-262.View ArticleGoogle Scholar
- Milc J, Infantino A, Pecchioni N, Aragona M: Identification of tomato genes differentially expressed during compatible interaction with Pyrenochaeta lycopersici. J Plant Pathol. 2012, 94 (2): 283-296.Google Scholar
- Clergeot P-H, Schuler H, Mørtz E, Brus M, Vintila S, Ekengren S: The corky root rot pathogen Pyrenochaeta lycopersici secretes a proteinaceous inducer of cell death affecting host plants differentially. Phytopathology. 2012, 102 (9): 878-891. 10.1094/PHYTO-01-12-0004.PubMedView ArticleGoogle Scholar
- Miller JR, Koren S, Sutton G: Assembly algorithms for next-generation sequencing data. Genomics. 2010, 95: 315-327. 10.1016/j.ygeno.2010.03.001.PubMed CentralPubMedView ArticleGoogle Scholar
- Rawat A, Elasri MO, Gust KA, George G, Pham D, Scanlan LD, Vulpe C, Perkins EJ: CAPRG: sequence assembling pipeline for next generation sequencing of non-model organisms. PLoS One. 2012, 7 (2): e30370-10.1371/journal.pone.0030370.PubMed CentralPubMedView ArticleGoogle Scholar
- Parra G, Bradnam K, Ning Z, Keane T, Korf I: Assessing the gene space in draft genomes. Nucleic Acids Res. 2009, 37: 289-297. 10.1093/nar/gkn916.PubMed CentralPubMedView ArticleGoogle Scholar
- Ter-Hovhannisyan V, Lomsadze A, Chernoff YO, Borodovsky M: Gene prediction in novel fungal genomes using an ab initio algorithm with unsupervised training. Genome Res. 2008, 18: 1979-1990. 10.1101/gr.081612.108.PubMed CentralPubMedView ArticleGoogle Scholar
- Roberts A, Pimentel H, Trapnell C, Pachter L: Identification of novel transcripts in annotated genomes using RNA-Seq. Bioinformatics. 2011, 27: 2325-2329. 10.1093/bioinformatics/btr355.PubMedView ArticleGoogle Scholar
- Ekman D, Elofsson A: Identifying and quantifying orphan protein sequences in fungi. J Mol Biol. 2010, 396: 396-405. 10.1016/j.jmb.2009.11.053.PubMedView ArticleGoogle Scholar
- Nowrousian M, Stajich JE, Chu M, Engh I, Espagne E, Halliday K, Kamerewerd J, Kempken F, Knab B, Kuo H-C, Osiewacz HD, Pöggeler S, Read ND, Seiler S, Smith KM, Zickler D, Kück U, Freitag M: De novo assembly of a 40 Mb eukaryotic genome from short sequence reads: Sordaria macrospora, a model organism for fungal morphogenesis. PLoS Genet. 2010, 6: e1000891-10.1371/journal.pgen.1000891.PubMed CentralPubMedView ArticleGoogle Scholar
- Islam MS, Haque MS, Islam MM, Emdad EM, Halim A, Hossen QMM, Hossain MZ, Ahmed B, Rahim S, Rahman MS, Alam MM, Hou S, Wan X, Saito J a, Alam M: Tools to kill: genome of one of the most destructive plant pathogenic fungi Macrophomina phaseolina. BMC Genomics. 2012, 13: 493-10.1186/1471-2164-13-493.PubMed CentralPubMedView ArticleGoogle Scholar
- Daskalov A, Paoletti M, Ness F, Saupe SJ: Genomic clustering and homology between HET-S and the NWD2 STAND protein in various fungal genomes. PLoS One. 2012, 7 (4): e34854-10.1371/journal.pone.0034854.PubMed CentralPubMedView ArticleGoogle Scholar
- Rouxel T, Grandaubert J, Hane JK, Hoede C, Van De Wouw P, Couloux A, Dominguez V, Anthouard V, Bally P, Bourras S, Cozijnsen AJ, Ciuffetti LM, Dilmaghani A, Duret L, Fudal I, Goodwin SB, Gout L, Glaser N, Linglin J, Kema GHJ, Lapalu N, Lawrence CB, May K, Meyer M, Ollivier B, Schoch CL, Simon A, Spatafora JW, Turgeon BG, Tyler BM, et al: Effector diversification within compartments of the Leptosphaeria maculans genome affected by repeat-induced point mutations. Nat Commun. 2011, 2: 20210.1038-View ArticleGoogle Scholar
- Ellwood SR, Liu Z, Syme RA, Lai Z, Hane JK, Keiper F, Moffat CS, Oliver RP, Friesen TL: A first genome assembly of the barley fungal pathogen Pyrenophora teres f. teres. Genome Biol. 2010, 11: R10910.1186-View ArticleGoogle Scholar
- Hane JK, Williams A, Oliver RP: Genomic and comparative analysis of the class Dothideomycetes. The Mycota. Edited by: Poggeler S, Wostemeyer J. 2011, Berlin: Springer-Verlag, 14: 205-226.Google Scholar
- Pyrenophora tritici-repentis database. http://www.broadinstitute.org/,
- Hall C, Welch J, Kowbel DJ, Glass NL: Evolution and diversity of a fungal self/nonself recognition locus. PLoS One. 2010, 5: e14055-10.1371/journal.pone.0014055.PubMed CentralPubMedView ArticleGoogle Scholar
- Milgroom MG, Sotirovski K, Risteski M, Brewer MT: Heterokaryons and parasexual recombinants of Cryphonectria parasitica in two clonal populations in southeastern Europe. Fungal Genet Biol. 2009, 46: 849-854. 10.1016/j.fgb.2009.07.007.PubMedView ArticleGoogle Scholar
- Tuite MF, Serio TR: The prion hypothesis: from biological anomaly to basic regulatory mechanism. Nat Rev Mol Cell Biol. 2010, 11: 823-833. 10.1038/nrm3007.PubMed CentralPubMedView ArticleGoogle Scholar
- Knogge W: Fungal infection of plants. Cell. 1996, 8: 1711-1722.Google Scholar
- Harris PV, Welner D, McFarland KC, Re E, Navarro Poulsen JC, Brown K, Salbo R, Ding H, Vlasenko E, Merino S, Xu F, Cherry J, Larsen S, Lo Leggio L: Stimulation of lignocellulosic biomass hydrolysis by proteins of glycoside hydrolase family 61: structure and function of a large, enigmatic family. Biochemistry. 2010, 49: 3305-3316. 10.1021/bi100009p.PubMedView ArticleGoogle Scholar
- Kostylev M, Wilson D: Synergistic interaction in cellulose hydrolysis. Biofuels. 2012, 3 (1): 61-70. 10.4155/bfs.11.150.View ArticleGoogle Scholar
- Ohm RA, Feau N, Henrissat B, Schoch CL, Horwitz BA, Barry KW, Condon BJ, Copeland AC, Dhillon B, Glaser F, Hesse CN, Kosti I, LaButti K, Lindquist EA, Lucas S, Salamov AA, Bradshaw RE, Ciuffetti L, Hamelin RC, Kema GH, Lawrence C, Scott JA, Spatafora JW, Turgeon BG, de Wit PJ, Zhong S, Goodwin SB, Grigoriev IV: Diverse lifestyles and strategies of plant pathogenesis encoded in the genomes of eighteen dothideomycetes fungi. PLoS Pathog. 2012, 8 (12): e1003037-10.1371/journal.ppat.1003037.PubMed CentralPubMedView ArticleGoogle Scholar
- Quinlan RJ, Sweeney MD, Lo Leggio L, Otten H, Poulsen JCN, Johansen KS, Krogh KBRM, Jorgensen CI, Tovborg M, Anthonsen A, Tryfona T, Walter CP, Dupree P, Xu F, Davies GJ, Walton PH: Insights into the oxidative degradation of cellulose by a copper metalloenzyme that exploits biomass components. P Natl Acad Sci USA. 2011, 108 (37): 15079-15084. 10.1073/pnas.1105776108.View ArticleGoogle Scholar
- Phillips CM, Beeson WT, Cate JH, Marletta MA: Cellobiose dehydrogenase and a copper-dependent polysaccharide monooxygenase potentiate cellulose degradation by Neurospora crassa. ACS Chem Biol. 2011, 6: 1399-1406. 10.1021/cb200351y.PubMedView ArticleGoogle Scholar
- Aragona M, Valente MT: Endoglucanase expression and virulence in plant fungal pathogens. The Fungal Cell Wall. Edited by: Mora-Montes HM. 2013, New York: Nova Publishers, 253-274.Google Scholar
- Flores A, Chet I, Herrera-Estrella A: Improved biocontrol activity of Trichoderma harzianum strains by overexpression of the proteinase encoding gene prb1. Curr Genet. 1997, 31: 30-37. 10.1007/s002940050173.PubMedView ArticleGoogle Scholar
- Pozo MJ, Baek JM, Garcia JM, Kenerley CM: Functional analysis of tvsp1, a serine protease-encoding gene in the biocontrol agent Trichoderma virens. Fungal Genet Biol. 2004, 41: 336-348. 10.1016/j.fgb.2003.11.002.PubMedView ArticleGoogle Scholar
- Suárez B, Rey M, Castillo P, Monte E, Llobell A: Isolation and characterization of PRA1, a trypsin-like protease from the biocontrol agent Trichoderma harzianum CECT 2413 displaying nematicidal activity. Appl Microbiol Biotech. 2004, 65: 46-55.View ArticleGoogle Scholar
- Viterbo A, Harel M, Chet I: Isolation of two aspartyl proteases from Trichoderma asperellum expressed during colonization of cucumber roots. FEMS Microbiol Lett. 2004, 238: 151-158.PubMedGoogle Scholar
- Carlile AJ, Bindschedler LV, Bailey AM, Bowyer P, Clarkson JM, Cooper RM: Characterization of SNP1, a cell wall-degrading trypsin, produced during infection by Stagonospora nodorum. Mol Plant-Microbe In. 2000, 13: 538-550. 10.1094/MPMI.2000.13.5.538.View ArticleGoogle Scholar
- Plummer KM, Clark SJ, Ellis LM, Loganathan A, Al-Samarrai TH, Rikkerink EHA, Sullivan PA, Templeton MD, Farley PC: Analysis of a secreted aspartic peptidase disruption mutant of glomerella cingulata. Eur J Plant Pathol. 2004, 110: 265-274.View ArticleGoogle Scholar
- Thon MR, Nuckles EM, Takach JE, Vaillancourt LJ: CPR1: a gene encoding a putative signal peptidase that functions in pathogenicity of colletotrichum graminicola to maize. Mol Plant-Microbe In. 2002, 15: 120-128. 10.1094/MPMI.2002.15.2.120.View ArticleGoogle Scholar
- Goodwin SB, M’barek SB, Dhillon B, Wittenberg AH, Crane CF, Hane JK, Foster AJ, Van der Lee TA, Grimwood J, Aerts A, Antoniw J, Bailey A, Bluhm B, Bowler J, Bristow J, van der Burgt A, Canto-Canché B, Churchill AC, Conde-Ferràez L, Cools HJ, Coutinho PM, Csukai M, Dehal P, De Wit P, Donzelli B, van de Geest HC, Van Ham RC, Hammond-Kosack KE, Henrissat B: Finished genome of the fungal wheat pathogen mycosphaerella graminicola reveals dispensome structure, chromosome plasticity, and stealth pathogenesis. PLoS Genet. 2011, 7 (6): e1002070-10.1371/journal.pgen.1002070.PubMed CentralPubMedView ArticleGoogle Scholar
- Duplessisa S, Cuomob CA, Linc Y-C, Aertsd A, Tisseranta E, Veneault-Fourreya C, Jolye DL, Hacquarda S, Amselemf J, Cantarelg BL, Chiuh R, Coutinhog PM, Feaue N, Fieldh M, Freya P, Gelhayea E, Goldbergb J, Grabherrb MG, Kodirab CD, Kohlera A, Küesi U, Lindquistd EA, Lucasd SM, Magoj R, Maucelib E, Morina E, Murata C, Pangilinand JL, Parkk R, Pearsonb M, et al: Obligate biotrophy features unraveled by the genomic analysis of rust fungi. P Natl Acad Sci USA. 2011, 108 (229): 1669171-Google Scholar
- Coleman JJ, Mylonakis E: Efflux in fungi: la piece de resistance. PLoS Pathog. 2009, 5: e1000486-10.1371/journal.ppat.1000486.PubMed CentralPubMedView ArticleGoogle Scholar
- Morschhauser J: Regulation of multidrug resistance in pathogenic fungi. Fungal Genet Biol. 2010, 47: 94-106. 10.1016/j.fgb.2009.08.002.PubMedView ArticleGoogle Scholar
- Ren Q, Chen K, Paulsen IT: TransportDB: a comprehensive database resource for cytoplasmic membrane transport systems and outer membrane channels. Nucleic Acids Res. 2007, 35: D274-D279. 10.1093/nar/gkl925.PubMed CentralPubMedView ArticleGoogle Scholar
- Keller NP, Turner G, Bennett JW: Fungal secondary metabolism - from biochemistry to genomics. Nat Rev Microbiol. 2005, 3: 937-947. 10.1038/nrmicro1286.PubMedView ArticleGoogle Scholar
- Friesen TL, Faris JD, Solomon PS, Oliver RP: Host-specific toxins: effectors of necrotrophic pathogenicity. Cell Microbiol. 2008, 10: 1421-1428. 10.1111/j.1462-5822.2008.01153.x.PubMedView ArticleGoogle Scholar
- Walton JD: HC-toxin. Phytochem. 2006, 67: 1406-1413. 10.1016/j.phytochem.2006.05.033.View ArticleGoogle Scholar
- Urban M, Bhargava T, Hamer JE: An ATP-driven efflux pump is a novel pathogenicity factor in rice blast disease. EMBO J. 1999, 18: 512-521. 10.1093/emboj/18.3.512.PubMed CentralPubMedView ArticleGoogle Scholar
- de Waard MA, Andrade AC, Hayashi K, Schoonbeek HJ, Stergiopoulos I, Zwiers LH: Impact of fungal drug transporters on fungicide sensitivity, multidrug resistance and virulence. Pest Manag Sci. 2006, 62: 195-207. 10.1002/ps.1150.PubMedView ArticleGoogle Scholar
- Cenis JL: Rapid extraction of fungal DNA for PCR amplification. Nucleic Acids Res. 1992, 20: 2380-10.1093/nar/20.9.2380.PubMed CentralPubMedView ArticleGoogle Scholar
- Scythe homepage. [https://github.com/vsbuffalo/scythe]
- Sickle homepage. [https://github.com/najoshi/sickle]
- Zerbino DR, Birney E: Velvet: algorithms for de novo short read assembly using de Bruijn graphs. Genome Res. 2008, 18 (5): 821-829. 10.1101/gr.074492.107.PubMed CentralPubMedView ArticleGoogle Scholar
- Huang X, Madan A: CAP3: a DNA sequence assembly program. Genome Res. 1999, 9: 868-877. 10.1101/gr.9.9.868.PubMed CentralPubMedView ArticleGoogle Scholar
- Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, Adiconis X, Fan L, Raychowdhury R, Zeng Q, Chen Z, Mauceli E, Hacohen N, Gnirke A, Rhind N, di Palma F, Birren BW, Nusbaum C, Lindblad-Toh K, Friedman N, Regev A: Full-length transcriptome assembly from RNA-seq data without a reference genome. Nat Biotechnol. 2011, 29: 1-17.View 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
- Smit AFA, Hubley R, Green P: RepeatMasker Open-3.0. 1996-2010. [http://www.repeatmasker.org]
- Trapnell C, Roberts A, Goff L, Pertea G, Kim D, Kelley DR, Pimentel H, Salzberg SL, Rinn JL, Pachter L: Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat Protoc. 2012, 7 (3): 562-578. 10.1038/nprot.2012.016.PubMed CentralPubMedView ArticleGoogle Scholar
- Kong L, Zhang Y, Ye ZQ, Liu XQ, Zhao SQ, Wei L, Gao G: CPC: assess the protein-coding potential of transcripts using sequence features and support vector machine. Nucleic Acids Res. 2007, 35: W345-W349. 10.1093/nar/gkm391.PubMed CentralPubMedView ArticleGoogle Scholar
- Wang H, Xu Z, Gao L, Hao B: A fungal phylogeny based on 82 complete genomes using the composition vector method. BMC Evol Biol. 2009, 9: 195-10.1186/1471-2148-9-195.PubMed CentralPubMedView ArticleGoogle Scholar
- O’Connell RJ, Thon MR, Hacquard S, Amyotte SG, Kleemann J, Torres MF, Damm U, Buiate E, Epstein L, Alkan N, Altmüller J, Alvarado-Balderrama L, Bauser C, Becker C, Birren BW, Chen Z, Choi J, Crouch JA, Duvick JP, Farman M, Gan P, Heiman D, Henrissat B, Howard RJ, Kabbage M, Koch C, Kracher B, Kubo Y, Law AD, Lebrun M-H, et al: Lifestyle transitions in plant pathogenic Colletotrichum fungi deciphered by genome and transcriptome analyses. Nat Genet. 2012, 44: 1060-1065. 10.1038/ng.2372.PubMedView ArticleGoogle Scholar
- Broad Institute of Harvard and MIT. [http://www.broadinstitute.org/]
- Unitè de Recherche Genomique info. [urgi.versailles.inra.fr]Google Scholar
- The Blumeria Sequencing Project. [http://www.blugen.org]
- LTR Finder. [http://tlife.fudan.edu.cn/cvtree/]
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215: 403-410. 10.1016/S0022-2836(05)80360-2.PubMedView ArticleGoogle Scholar
- NCBI Non redundant protein database. [ftp://ftp.ncbi.nih.gov/blast/db]
- Uniprot database. [http://www.uniprot.org]
- Conesa A, Götz S, Garcia-Gomez JM, Terol J, Talon M, Robles M: Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005, 21: 3674-3676. 10.1093/bioinformatics/bti610.PubMedView ArticleGoogle Scholar
- Interpro. [http://www.ebi.ac.uk/interpro/]
- Eddy SR: Accelerated profile HMM searches. PLoS Comput Biol. 2011, 7 (10): e1002195-10.1371/journal.pcbi.1002195.PubMed CentralPubMedView ArticleGoogle Scholar
- PFam. [http://pfam.sanger.ac.uk]
- dbCAN. [http://csbl.bmb.uga.edu/dbCAN/]
- PHI-base. [http://www.phi-base.org]
- MEROPS. [http://merops.sanger.ac.uk]
- C2H2 ZNF db. [http://kzfgd.pzr.uni-rostock.de:8080/KZGD2007]
- TCDB. [http://www.tcdb.org]
- GPCRDB. [http://www.gpcr.org]
- Benjamini Y, Hochberg Y: Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B. 1995, 57: 289-300.Google Scholar
- Li H, Durbin R: Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009, 25: 1754-1760. 10.1093/bioinformatics/btp324.PubMed CentralPubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited.