Genome sequence of the necrotrophic fungus Penicillium digitatum, the main postharvest pathogen of citrus
- Marina Marcet-Houben†1, 2,
- Ana-Rosa Ballester†3,
- Beatriz de la Fuente3,
- Eleonora Harries3,
- Jose F Marcos3,
- Luis González-Candelas3Email author and
- Toni Gabaldón1, 2Email author
© Marcet-Houben et al.; licensee BioMed Central Ltd. 2012
Received: 15 May 2012
Accepted: 9 November 2012
Published: 21 November 2012
Penicillium digitatum is a fungal necrotroph causing a common citrus postharvest disease known as green mold. In order to gain insight into the genetic bases of its virulence mechanisms and its high degree of host-specificity, the genomes of two P. digitatum strains that differ in their antifungal resistance traits have been sequenced and compared with those of 28 other Pezizomycotina.
The two sequenced genomes are highly similar, but important differences between them include the presence of a unique gene cluster in the resistant strain, and mutations previously shown to confer fungicide resistance. The two strains, which were isolated in Spain, and another isolated in China have identical mitochondrial genome sequences suggesting a recent worldwide expansion of the species. Comparison with the closely-related but non-phytopathogenic P. chrysogenum reveals a much smaller gene content in P. digitatum, consistent with a more specialized lifestyle. We show that large regions of the P. chrysogenum genome, including entire supercontigs, are absent from P. digitatum, and that this is the result of large gene family expansions rather than acquisition through horizontal gene transfer. Our analysis of the P. digitatum genome is indicative of heterothallic sexual reproduction and reveals the molecular basis for the inability of this species to assimilate nitrate or produce the metabolites patulin and penicillin. Finally, we identify the predicted secretome, which provides a first approximation to the protein repertoire used during invasive growth.
The complete genome of P. digitatum, the first of a phytopathogenic Penicillium species, is a valuable tool for understanding the virulence mechanisms and host-specificity of this economically important pest.
Postharvest losses can have a significant impact on crops, ranging from 15 to as high as 50% of total production . Citrus is one of the most economically important fruit crops in the world and it is particularly susceptible to postharvest damage because harvested fruits are usually stored before they reach the market for fresh consumption, a period in which they are subject to both biotic and abiotic stress conditions. Fungi, which are particularly adapted to a saprophytic lifestyle, are among the main biological agents causing crop deterioration during transport and storage. In citrus, a common postharvest disease known as green mold and caused by Penicillium digitatum can account for up to 90% of the total losses, especially in arid and sub-tropical climates . Nowadays, the control of post-harvest decay is achieved by the massive application of fungicides to marketed fruits. This however, has serious implications regarding toxicity, consumer acceptance and environmental risk. In this sense, reports of selection for fungicide resistant strains are increasing and represent a significant obstacle to postharvest conservation and commerce . P. digitatum is a necrotrophic wound pathogen that requires pre-existing injured fruit peel to penetrate the plant tissue, and it colonizes mostly through the deployment of maceration enzymes. Remarkably, and despite this rather unspecific infection mechanism, P. digitatum exhibits a high degree of host specificity and has not been described as ‘naturally-occurring’ in other pathosystems outside citrus fruits . Despite its considerable economic interest, the molecular bases of P. digitatum infection and host specificity remain largely unknown. We have previously analyzed the fruit’s transcriptomic response to P. digitatum infection and also to elicitors that trigger induced resistance [5, 6]. Here, we report the complete genome sequence of two P. digitatum strains that differ in their resistance to common fungicides and a comprehensive comparison across 28 other sequenced Pezizomycotina. This resource will enable progress to be made towards understanding the physiology and virulence mechanisms of this economically important phytopathogen.
Results and discussion
Summary of the main assembly and annotation features of the genomes of the two sequenced Penicillium digitatum strains
P. digitatum strain
Illumina (HiSeq 2000)
Number of contigs
Number of scaffolds
Number of Large Scaffolds (>100 Kb)
N50 (base pairs)
Number of indetermined sites (Ns)
GC content (%)
Number of genes
Mean gene length (nucleotides)
Percentage of genes with introns
Mean number of exons per gene
Sexual locus, nitrate utilization and iron uptake
As in many other ascomycetous fungi, no teleomorfic form or sexual reproduction cycle has been described for P. digitatum. We nevertheless found a conserved alpha-box mating-type protein (MAT1) locus in both strains (Additional file 1: Figure S1) in an area of conserved synteny with other Penicillium and Aspergillus species suggesting the possibility of sexual reproduction in these species. In addition, P. digitatum shows the organization typical of heterothallic Pezizomycotina, with the presence of a single MAT locus .
Comparison of the two sequenced strains
The two sequenced strains are highly similar with pairs of orthologs sharing on average 99.96% identity at the protein level (99.95%, at the nucleotide level). This similarity extends to the chromosomal organization between the two strains, with 99.8% of neighboring gene pairs being conserved. We could only detect two putative chromosomal re-arrangements (Additional file 1: Figure S2). When aligning the reads of one strain to the other strain several coding regions remained unmapped suggesting that they could have been either lost or gained specifically in only one of the strains. More specifically, two genes were unique for strain Pd1 and eight for PHI26 (Additional file 1: Table S2), the two unique sequences in the resistant strain being contiguous in the assembly. The differences in the number of predicted genes for the two strains (Table 1) are actually larger than those shown in Additional file 1: Table S2. This is mostly due to genes not predicted in Pd1 as a consequence of assembly gaps, since 52% of the PHI26-unique genes mapped to regions containing Ns and 48% appeared truncated in the Pd1 assembly. The prediction of unique genes based on the mapping of reads from the other strain (Additional file 1: Table S2) is independent of differences in the assembly or gene prediction quality and thus should be considered to be more accurate.
Few single nucleotide polymorphisms (SNPs) were noted between the two strains. In total 1,441 high confidence SNPs were found, of which 553 affected coding sequences. Remarkably, there is an excess of non-synonymous (358) vs synonymous (195) SNPs, suggesting the possibility that part of the differences are the result of positive selection. Among the 22 genes with more than one non-synonymous SNPs (Additional file 1: Table S3) we found one possibly having a regulatory function (PDIG_27530) and two putative membrane transporters (PDIG_44900 and PDIG_41990) that could be involved in the extrusion of toxic compounds. In this regard it is interesting to note that P. digitatum Pd1 is significantly more resistant than PHI26 to several additional fungicides imazalil and thiabendazole (see below).
Among the detected SNPs we confirmed the presence of the Phe-to-Tyr mutation at codon 200 in the β-tubulin gene which has been associated with thiabendazole resistance in P. digitatum. However, a second mutation affecting this gene, and related to benomyl resistance, was not found in either strain . Of particular note, we found two novel mutations in Pd1 both of which cause coding sequence truncations, one in a gene (PDIP_45840) encoding a putative member of a phosphotransferase family (APH, PFAM: PF01636) implicated in antibiotic resistance in bacteria , and another in a gene (PDIG_06580) predicted to encode a putative hydrolase. The possible involvement in fungicide resistance of these mutations and some of the remaining differences between the two strains will require experimental validation. Going one step in this direction, we specifically tested for this in the two abovementioned Pd1-specific genes by constructing the corresponding knock-out strains and exposing them to a battery of six widely-used fungicides (See Materials and Methods and Additional file 1). Neither of the two deleted strains showed significant differences in either fungicide resistance or their pathogenicity towards citrus fruits with respect to the wild-type Pd1 control or an ectopic (non-deleted) transformant (Additional file 1: Figures S3 to S6, and Additional file 1: Table S4). We thus discount a direct involvement of these two genes in resistance to the six fungicides tested.
In addition to the differences already commented, our analyses also confirm the presence of a 199 bp insertion in the promoter of the cyp51B gene in Pd1 that has been shown to confer resistance to imazalil by increasing gene expression . No other variations were found in any of the three cyp51 genes described that could account for the differential susceptibility to imazalil observed between isolates Pd1 and PHI26. Notably, this 199 bp motif is represented 9 times in the Pd1 genome and 12 times in PHI26. The locations of eight of these motifs are common between the strains but all except the resistance-associated insertion fall distant from protein-coding genes. Comparison between the two strains revealed that the motifs absent in Pd1 correspond to genomic regions that are either unassembled or unresolved for this strain. Remarkably, the insertion site of the 199 bp motif in the cyp51B promoter is identical to that found in two other imazalil-resistant isolates from China . Given the few differences found between the two sequenced Spanish strains and the observation that only a fraction of them are directly related to the resistance phenotype indicates that the appearance of fungicide resistance in P. digitatum is likely to involve few genomic mutations.
The high degree of genomic conservation also extends to the mitochondrial genomes of Pd1 and PHI26. They exhibit very few differences and all involve small variations in the lengths of polynucleotide tracks, a phenomenon attributable to known sequencing error biases in the 454 technology . It is thus highly likely that Pd1 and PHI26 have identical mitochondrial DNA sequences. We next compared this sequence to the recently reported mitochondrial sequence from the Chinese isolate pd01 , and this also turned out to be identical. This indicates a remarkable conservation of mitochondrial DNA sequences in distant isolates, and is suggestive of a recent expansion of P. digitatum around the globe. The cultivation of citrus has its origin in Southeast Asia, and, while its presence in Southern Europe can be dated back to the times of the Roman Empire, intensive production and massive worldwide trade only developed as recently as the 1940s . Our results indicate that the P. digitatum strains isolated in Spain diverged from the Chinese strain very recently and argue in favor of a recent expansion of this pathogenic species, probably coinciding with the initiation of intensive agricultural practices and worldwide trade.
Comparison across species
Based on the phylome analysis we found that 843 P. digitatum proteins lack homologs in the closely-related but non-pathogenic species P. chrysogenum. Of these, 97 have homologs exclusively in other plant pathogenic species thus highlighting them as good candidates for possible roles in virulence (Additional file 1: Table S7). Other genes in P. digitatum have homologs in a limited number of fungal species while having widespread bacterial homologs, a pattern that could be indicative of Horizontal Gene Transfer (HGT). Previous analyses  have found that genomes from Pezizomycotina species are particularly enriched for genes acquired from prokaryotes. We have applied the same strategy to find genes of putatively prokaryotic origin in P. digitatum. Four genes fulfilled our criteria (Additional file 1: Table S8). Notably, one of these (PDIG_81140) was identified as having homologs in plant pathogenic fungi but not in the non-pathogenic Penicillium, and also a homolog in a bacterial family (DEC1) that has been associated with virulence in maize infections and the synthesis of T-toxin in Cochliobolus heterostrophus. None of the putatively transferred genes identified here in P. digitatum correspond to genes that have been described as HGT in P. chrysogenum, among which the penicillin gene cluster is included (see below). Another relevant HGT described in some fungi , but absent in P. digitatum, is arsenate reductase (arsC) which in P. chrysogenum forms part of an arsenate resistance cluster. Remarkably, we also identified the absence of homologs to the alternative arsenate reductases present in Saccharomyces cerevisiae (Acr2p) or in other distantly-related fungi . This would suggest that P. digitatum may be vulnerable to arsenic, or alternatively, that this species may be using another detoxification system not previously described in fungi. We next focused on five particular aspects of P. digitatum genomic content that are likely to provide hints for understanding this organism’s pathogenic lifestyle and constitute interesting targets for future research: namely i) secondary metabolism gene clusters, ii) the prediction of putative small cysteine-rich peptides (CRPs), iii) homologs of known effector proteins iv) carbohydrate active enzymes (CAZy), and finally iv) the repertoire of secreted proteins among which we expect to find those that potentially interact with the host and are involved in invasive growth of the fungus.
Small, cysteine-rich peptides and effector proteins
Several fungi are known to secrete small cysteine-rich proteins (CRPs) that exhibit diverse biological functions including adherence , virulence [38, 39] or antimicrobiosis [40, 41]. Some CRPs are virulence effectors that show carbohydrate binding activity that interferes with host recognition of the pathogen . It has been suggested that antimicrobial proteins could become extremely useful in medicine and agriculture to treat pathogenic infections in the near future [41, 43]. Several antifungal CRPs exhibiting antimicrobial activity have been characterized from ascomycetous fungal species from the genera Aspergillus and Penicillium. These cationic proteins share the previously described γ-motif common to small antimicrobial proteins and peptides . Two distinct antifungal CRPs have been isolated and characterized from cultures of P. chrysogenum[45, 46]. Among the predicted proteins in P. digitatum we found a CRP (PDIG_68840) showing 91% amino acid identity with the Pc12g08290 gene from P. chrysogenum which corresponds to one of the previously reported proteins . Remarkably, the ortholog to Pc24g00380 (the most studied antifungal protein from this fungus) is absent in P. digitatum as it is located in the Pc24 supercontig which is virtually devoid of P. digitatum orthologs. We identified two additional putative CRPs in P. digitatum (PDIG_31210 and PDIG_70730) by performing profile-searches with CRP motifs identified in plants . Due to their short length, CRPs are likely to be under-predicted in the annotation of sequenced genomes. This prompted us to extend the genome annotation with a specific search (see Materials and Methods) for additional small CRPs that contain at least four cysteine residues, display similarities to known CRPs, harbor the γ-motif and/or have a high content of cationic residues typical of antimicrobial proteins, under the assumption that these might have been under-predicted in the initial annotation (see Materials and Methods). We found 73 additional ORFs that potentially would code for CRPs, ranging in size from 23 to 146 residues (Additional file 1: Table S11). Two of these putative CRPs (Pdigorf_18868 and Pdigorf_93050) share the γ-motif, making them good candidates for having antifungal activity.
We next searched the pathogen-host interaction database  for P. digitatum proteins homologous to effector proteins from other pathogenic fungi . These include not only CRPs but also other types of proteins that can modulate host physiology. Again, P. digitatum encodes a smaller number of homologs of known effector proteins compared to P. chrysogenum, but abundances relative to the total proteome are comparable (Figure 5). We found eight P. digitatum homologs to Ace1 of Magnaporte grisea, six to HopI1 (Pseudomonas syringae), and single homologs to Npp1 (Hyaloperonospora parasitica), Gip1 and Gip2 (Phytophtora infestans). Additionally, five proteins in P. digitatum contain the LysM domain (Additional file 1: Table S12) of which two have no orthologs in P. chrysogenum and are closely related to A. nidulans and Aspergillus carbonarius instead. Of note, the fraction of proteins with homologs in the known effector proteins dataset is roughly similar (0.01–0.2%) across all the pezyzomycotina genomes analyzed regardless of their phytopathogenicity. This result illustrates the limitations of using mere homology to predict the functional roles of proteins.
Carbohydrate active enzymes (CAZymes)
Fungi use secreted enzymes to break-down biopolymers the resulting products of which are then transported inside cells for further metabolic processing, and in phytopathogenic fungi secreted proteins constitute the first line of interaction with the host and play important roles in virulence [51–53]. Since secreted fungal enzymes have found a wide range of applications in the food, feed, pulp and paper, bioethanol and textile industries there is a growing interest in deciphering the complement of secreted proteins in sequenced genomes (i.e. the secretome). Using an in silico pipeline (Methods) we predicted 552 secreted proteins, of which 133 (24%) code for various hydrolases such as proteases, lipases or nucleases (Additional file 1: Table S14) Among the pezizomycotina species considered, P. digitatum has the second smallest predicted proteome after Uncinocarpus reesii (Additional file 1: Table S5). Similarly, P. digitatum also has one of the smallest secretomes comprising 552 proteins (6% of the proteome), which in relative terms is somewhat smaller than the secretome of P. chrysogenum (914 proteins, 7%) and the average 8% found in Aspergillus species, and much lower than the 14% of secreted proteins found in the Magnaporthe species (Figure 5c). That P. digitatum is a necrotroph and contains one of the smallest secretomes among analyzed fungi suggests that the absolute number of secreted proteins has no direct relationship with the lifestyle of a pathogen.
Previous analyses comparing the genomic features of necrotrophic fungal pathogens  lacked the inclusion of a necrotrophic eurotiomycete. Our comparison of two P. digitatum genomes from strains that differ in their fungicide resistance traits with those of other sequenced pezizomycotina thus constitutes the first comprehensive genomic analysis of a necrotrophic organism within this important fungal clade that includes other plant pathogens such as those in the Aspergillu s genus. Our data suggest a recent global expansion of P. digitatum, consistent with the relatively modern origin of intensive cultivation and worldwide commercialization of citrus. In addition, our comparative analysis suggests the action of positive selection and that important differences in fungicide resistance seem to have been achieved through a small number of genomic changes. Our analyses also uncovered other strain-specific differences, the functional roles of which will be worth exploring. We show a striking reduction in gene content in P. digitatum compared to the close relative P. chrysogenum. Large regions in the P. chrysogenum genome, including entire supercontigs, are completely absent from P. digitatum and originated via large gene family expansions rather than through horizontal gene transfer. Although we could only detect enrichment of two functional classes among the genes lost in the P. digitatum lineage, detailed analysis revealed important physiological implications such as the loss of the ability to use nitrate as a nitrogen source. Other differences extend to secreted proteins and CAZyme-encoding genes, and may serve as a basis to understand the host-specificity of P. digitatum. With respect to pathogenicity determinants, and consistent with analysis of other phytopathogenic genomes , there seem to be no global differences that distinguish these pathogens from other pathogenic and non-pathogenic fungi. These findings reinforce the view that necrotrophic pathogenicity is a very complex trait that involves many genes that may vary in different lineages. Our data pave the way for a more systematic functional and experimental analysis of candidate genes that will eventually reveal the details of the complex interaction between P. digitatum and its host.
Penicillium digitatum (Pers.:Fr.) Sacc. PHI26 and Pd1 strains have been used in this study. PHI26 was isolated from an infected orange  and Pd1 from an infected grapefruit , both in Valencia (Spain). Both strains have been deposited with the Spanish Type Culture Collection (CECT) with accession codes CECT20796 and CECT20795, respectively.
RNA-Seq data were generated from four different libraries: a cDNA library synthesized from P. digitatum spores collected after 7 days of growth in potato dextrose agar, and three libraries from P. digitatum-infected oranges at 12, 24, and 48 h post inoculation. Fruit inoculation and RNA extraction were conducted following procedures previously described . cDNA libraries were generated using the MINT-Universal cDNA synthesis kit (Evrogen) following the manufacturer’s protocol. Sequencing was carried out using the Roche 454 GS FLX Titanium system. The raw data files are available in the Sequence read archive at the National Center for Biotechnology Information (NCBI), GenBank ID: SRA059533.
Genomic DNA extraction
P. digitatum PHI26 and Pd1 (5 × 105 conidia/mL) were grown in 500 mL Potato Dextrose Broth media (PDB, Difco) on a shaker at 180 rpm at 24°C for 4 days. Fungal mass was collected, dried and frozen in liquid nitrogen. DNA was extracted from 1 g of frozen mycelium essentially as previously described . The final DNA preparation was incubated overnight at room temperature in 490 μL of TE buffer and 10 μL of DNase-free RNase (10 μg/mL), followed by phenol-chloroform extraction and isopropanol precipitation. Finally, DNA was resuspended in TE buffer. DNA quality was assessed through spectrophotometry and CHEF electrophoresis (data not shown).
Genomic sequencing and assembly
Genome sequencing was performed at the CRG ultrasequencing facility (PHI26) and at LifeSquencing (Pd1) using Illumina HiSeq 2000 and Roche 454 GS FLX Titanium technologies, respectively. For the Illumina sequencing (PHI26) DNA was fragmented by nebulization to a size ~400-600 bp. After shearing, the ends of DNA fragments were blunted with T4 DNA polymerase and Klenow fragment (New England Biolabs). DNA was purified with a QIAquick PCR purification kit (Qiagen). Thereafter, 3′-adenylation was performed by incubation with dATP and 3′-5′-exo− Klenow fragment (New England Biolabs). DNA was purified using MinElute spin columns (Qiagen) and double-stranded Illumina paired-end adapters were ligated to the DNA using rapid T4 DNA ligase (New England Biolabs). After another purification step, adapter-ligated fragments were enriched, and adapters were extended by selective amplification in an 18-cycle PCR reaction using Phusion DNA polymerase (Finnzymes). Libraries were quantified and loaded into Illumina flow-cells at concentrations of 7–20 pM (Genome Analyzer IIx; GA), or 1.4–1.75 pM (HiSeq 2000; HS). Cluster generation was performed in an Illumina cluster station (GA) or in a cBOT (HS). Sequence runs of 2 × 50 cycles were performed on the sequencing instruments. Base calling was performed using Illumina pipeline software. In case of multiplexed libraries, we used 4 bp internal indices (5′ indexed sequences). De-convolution was performed using the CASAVA software (Illumina). A total of 3.4 Gb of pair-end (350 insert size) and 2.6 Gb of mate-pairs (5 kb inserts), were produced. Reads were trimmed at the base with a Phred quality score lower than 10, or at the first undetermined base. Reads shorter than 31 bases (and their corresponding pairs) were removed. Reads were assembled using SOAPdenovo . The best Kmer size was chosen based on optimization of the N50, assembly length, largest contig length, and number of contigs. Two suboptimal assemblies were consistent with the main results presented here. Gaps were then filled in using Gapcloser . We exploited the similarity of the two sequenced strains and used OSLay  to create scaffolds. Length of gaps were controlled and contigs joined through large stretches (>10Kb) of Ns were separated.
For the Roche 454 sequencing (Pd1), a random “shotgun” genomic library was generated via fragmentation of 10 μg of P. digitatum Pd1 genomic DNA employing the GS FLX Titanium General Library Preparation Kit following the manufacturer’s recommendations (454 Life Sciences, Branford, CT). Briefly, DNA was randomly sheared via nebulization and double stranded DNA adaptors were blunt-ligated to fragment ends following post-electrophoresis agarose gel excision of the 500–800 bp fraction. The final single-stranded DNA library was isolated via streptavidin bead binding to biotinylated adaptors followed by alkaline treatment. The library was quantitated via fluorometry using Quant-iT RiboGreen reagent (Invitrogen, Carlsbad, CA) prior to emulsion PCR amplification. Subsequently, a paired end library was constructed. Briefly, 30 μg of double stranded genomic DNA was randomly fragmented via hydrodynamic shearing to an average size of 15,000 or 20,000 bp using the HydroShear apparatus (Genomic Solutions, Ann Arbor, MI). Fragment ends were blunt ended via T4 DNA polymerase and T4 polynucleotide kinase (Roche Applied Science, Indianapolis, IN) treatment, and a double-stranded DNA adaptor (containing a loxP recognition site and internal biotin moiety) was blunt-ligated to DNA ends. Library fragments in the desired size range were isolated from excised gel fragments via electroelution following overnight 0.5% agarose gel electrophoresis. Nicks present at the 3′-junctions of each of the adaptors and the library fragment were filled in by the strand-displacement activity of the Bst DNA polymerase (New England BioLabs, Ipswich, MA). Linear library molecules were circularized via Cre recombinase (New England BioLabs, Ipswich, MA) excision. These molecules along with 1 μg carrier DNA were randomly fragmented via nebulization to an average size of approximately 450–500 bp. Fragment ends were blunted and desired library fragments (linker flanked by original DNA fragment ends) were selected via the biotinylated linker, using 1 Dynal M-270 magnetic streptavidin beads (Invitrogen, Carlsbad, CA). Double-stranded DNA adaptors were ligated to blunt fragment ends. The resulting library was amplified using 20 cycles of PCR and purified with AMPure beads (Agencourt Bioscience, Beverly, MA, USA) to remove small DNA fragments and amplification primers. Single stranded DNA library was isolated via Dynal M-270 bead binding followed by alkaline treatment and quantitation using Quant-iT RiboGreen (Invitrogen, Carlsbad, CA) prior to emulsion PCR amplification. Genomic shotgun library molecules were clonally amplified via emulsion PCR following manufacturer’s recommendations employing the GS FLX Titanium LV emPCR Kit (454 Life Sciences, Branford, CT). Paired end library was handled similarly, with the exception that 30 μL of amplification primer was used per Large Volume Emulsion cup. Following amplification, emPCR reactions were collected, and emulsions broken according to the manufacturer’s protocols. Beads containing sufficient copies of clonally amplified library fragments were selected via the specified enrichment procedure and counted with a Multisizer 3 Coulter Counter (Beckman Coulter, Fullerton, CA) prior to sequencing. Following emulsion PCR enrichment, beads were deposited into the wells of a Titanium Series PicoTiterPlate device and 454 Sequencing was performed using the GS FLX instrument according to the manufacturer’s recommendations (454 Life Sciences, Branford, CT). Image analysis, signal processing and base calling were performed using supplied system software. Standard Flowgram Format (sff) files output from base calling were employed in subsequent genome assembly. The assembly of 454 reads was performed with Newbler (454-Roche).
We used a genome annotation pipeline similar to that used for the closely-related P. chrysogenum genome . In brief, we combined predictions from three alternative methods (Additional file 1: Figure S9): i) a first annotation was obtained using Augustus trained with A. nidulans files. RNA-Seq data from P. digitatum Pd1 in different growth conditions were incorporated in the predictions as hints. ii) a second prediction was done using JIGSAW  in order to unify the predictions of several predicting programs and methods which included a) Augustus  (without RNA-Seq data), b) GeneID , and c) SNAP , which were run using training files based on either P. chrysogenum or A. nidulans, depending on availability; d) a custom python script was used to detect Open Reading Frames (ORFs) using the universal genetic code, e) homology data obtained from Blast searches for each protein encoded in the three sequenced Penicillium species (P. chrysogenum, Penicillium marneffei and Talaromyces stipitatus) and several other Pezizomycotina species (A. nidulans, A. flavus, Neurospora crassa, Histoplasma capsulatum and Gibberella zeae) run against the P. digitatum genomes; f) GMAP  was used to map the genes encoded in P. chrysogenum to the newly sequenced P. digitatum genomes. iii) A third prediction was run using GeneWise  to search for all the proteins encoded in several fungal genomes (P. chrysogenum, P. marneffei, T. stipitatus, A. nidulans, Fusarium oxysporum, A. flavus, N. crassa, G. zeae, H. capsulatum) and the set of revised proteins found in Uniprot (as of June 2011). An initial Blast was performed to locate the contig in which a protein was potentially present, then GeneWise was run to establish its most likely position within the contig. Hits with a bit score over 50 were revised and extended in case of lack of starting or stop codon. If no internal stop codons were found in the final gene, this was added to the annotation. The final annotation was build using the first prediction as a starting point. Proteins detected in the second and third annotations where added when there was no overlapping with the initial set. The final predicted gene set included 9,153 and 8,969 genes for PHI26 and Pd1, respectively. Pfam domains  were searched using HMMER v3 with default values. GO terms and E.C. codes and mapping to KEGG  pathways were transferred from several of the fungal genomes when a one-to-one orthology relationship existed. A total of 6403 (70%) proteins contained at least some kind of annotation. Tansposons and repetitive elements were predicted with RepeatMasker .
The subcellular location of each protein in the 29 genomes used in this study was predicted using PSORT . Additionally, all proteomes were scanned to look for proteins containing a signal peptide using SignalP . These proteins were then scanned for the presence of transmembrane domains using TMHMM . 552 proteins in PHI26 and 547 proteins in Pd1 were considered as secreted proteins.
Prediction of small cystein - rich proteins (CRP)
Prediction of CRPs was based according to their expected sequence characteristics . Known proteins and predicted open reading frames (ORFs) of small (20 to 150 amino acids) size, which presented a predicted signal peptide from the secretory pathway and which contained at least four cysteins were pre-selected. This resulted in 6,961 short proteins in the seven genomes analyzed: P. digitatum (both strains), P. chrysogenum, A. clavatus, A. fumigatus, A. flavus, Aspergillus niger. Subsequently, HMM-derived profiles of 516 plant CRP motives described  were used to search the above mentioned collection of proteins using HMMER3 . Hits were considered if they had an e-value below 0.1. Peptides hit by the same CRP motif and showing a level of identity greater than 20% were grouped and aligned to derive a new HMM profile, which was then used as bait in a novel search. This process was repeated iteratively until no further peptides were added. This procedure yielded 73 putative CRPs.
The phylome, meaning the complete collection of phylogenetic trees for each gene in a genome, was reconstructed for the genomes of P. digitatum (strain PHI26) and of P. chrysogenum. The phylome was reconstructed using 28 other Pezizomycotina species. Considering the high degree of similarity between the two P. digitatum strains only the PHI26 strain was used in the phylomes. The phylome was reconstructed using an automated pipeline previously described in . Briefly, for each protein in a genome a Smith-Waterman search was performed against the fungal proteome database (Additional file 1: Table S5). Results were filtered using an e-value < 1e-05 and a continuous overlapping region of 0.5. At most 150 homologous sequences for each protein were accepted. Homologous sequences were then aligned using three different programs: MUSCLE v3.7 , MAFFT v6.712b , and DIALIGN-TX . Alignments were performed in forward and reverse direction (i.e. using the Head or Tail approach ), and the six resulting alignments were combined with M-Coffee . This combined alignment was trimmed with trimAl v1.3 , (consistency-score cutoff 0.1667, gap-score cutoff 0.9). Trees were reconstructed using the best-fitting evolutionary model. The selection of the model best fitting each alignment was performed as follows: a Neighbour Joining (NJ) tree was reconstructed as implemented in BioNJ ; the likelihood of this topology was computed, allowing branch-length optimization, using seven different models (JTT, LG, WAG, Blosum62, MtREV, VT and Dayhoff), as implemented in PhyML v3.0 ; the two models best fitting the data, as determined by the AIC criterion , were used to derive ML trees. Four rate categories were used and invariant positions were inferred from the data. Branch support was computed using an aLRT (approximate likelihood ratio test) based on a chi-square distribution. Resulting trees and alignments are stored in phylomeDB  (http://phylomedb.org), with the phylome IDs 150 (P. digitatum) and 179 (P. chrysogenum). Trees were scanned automatically using ETE v2 .
Orthologs between the two P. digitatum strains were established using a best bidirectional hit approach. Orthologs between P. digitatum and the other species were based on phylogenies included in the phylome. A species-overlap algorithm, as implemented in ETE v2  was used to infer orthology and paralogy relationships. Briefly the algorithm decides whether a node in a tree is a speciation of a duplication node depending on the overlap of the species branching from the node . Overlap between those species will indicate a duplication node. Otherwise a speciation node will be considered.
Detection of horizontal gene transfers
In order to detect genes putatively acquired from prokaryotes we used the procedure described in , using a database of 102 completely-sequenced fungi, 95 non-fungal eukaryotes, and 1395 eukaryotes as downloaded from KEGG as of June 2011 . Blast results were filtered using the same parameters as in the phylome reconstruction. Only proteins present in less than 10 fungal species, no other non-fungal eukaryotes and that were over-represented in bacterial species were considered to have a putative prokaryotic origin. These families were further analyzed phylogenetically.
Species tree reconstruction
The species tree was build using a concatenation method. 592 widespread, single copy genes were selected. After concatenation, the alignment was trimmed using trimAl . Columns with more than 50% of gaps were removed. A conservation score of 50% of the alignment was used. The final alignment had 371,529 positions. The tree was reconstructed using RAxML version:7.2.6 . LG model  was selected and a 4-categories GAMMA distribution was used. Bootstrap was obtained by creating 100 random sequences using SeqBoot from the Phylip package . A tree was then reconstructed for each sequence and the consensus tree was inferred using Phylip. Additionally, we reconstructed a species tree based on the super-tree reconstruction program DupTree , which took as input the 7,723 trees obtained during phylome reconstruction. Both trees resulted in the same final topology.
Percentage of identity
Pairs of orthologous proteins were established from the phylogeny-based orthology predictions. Between P. digitatum and the other fungal species only one-to-one orthologs were considered. A pair-wise alignment between the orthologous pair was reconstructed using MUSCLE v3.7 . The percentage of identity was calculated using the sident option of trimAl . As the two strains of P. digitatum are highly similar, proteins that had a length difference superior to the average difference of all protein lengths were not considered.
Gene order conservation
The degree of gene order conservation between two species was calculated as the pairs of orthologous genes that were found consecutively in both genomes and were found in the same direction. Up to five genes without orthologs were allowed between a pair of genes.
Detection of single nucleotide polymorphisms (SNPs) and assessment of evolutionary rates
SNPs were detected by aligning the short reads produced during the sequencing of P. digitatum PHI26 to the genome of P. digitatum Pd1. Only reads mapping unambiguosly were considered. SAMtools  was then used to call variants, considering that the organism is haploid. Only SNPs in regions covered at least by 10 reads and with an agreement in the base call of >90% were considered.
Ratios of non-synonymous over synonymous substitutions (dN/dS) were calculated as follows. Alignments of the relevant one-to-one orthologs (P. digitatum vs P. chrysogenum for the nitrogen cluster and P. digitatum vs A. clavatus for patulin) were computed on the protein sequences using MUSCLE v3.7 . The amino acid alignment was then back translated to nucleotides using the backtrans option implemented in trimAl v1.3 . dN/dS was then calculated as the number of codons that codified for different amino acids in the two species versus the number of codons that contained synonymous mutations. As a control, the same analysis was performed against a set of proteins involved in the glycolysis pathway and for all one-to-one orthologs detected between these species.
Experimental assessment of a possible role in fungicide resistance of selected genes
To investigate the possible role in the differential susceptibility to fungicides of the two contiguous genes that were only present in the Pd1 strain, we transformed this strain via Agrobacterium tumefaciens-mediated transformation with plasmid pRFD519. This plasmid derives from plasmid pRFHU2  and contains the flanking regions of the two genes surrounding the hygromycin resistant cassette in the T-DNA region of the plasmid. Primers PD519O1 (5′- GGTCTTAAUGGGGAGAGCATAGGTGGAAT-3′) and PD519O2 (5′-GGCATTAAUCTTGGATGGAGAGGGAACAA-3′) were used to amplify a fragment of 1.76 kb region upstream of the genes and primers PD519A3 (5′- GGACTTAAUCGTGACATAAGGGCCAAACT-3′) and PD519A4 (5′- GGGTTTAAUCCCCTTTGGTTACCCTCATT-3′) were used to amplify a 1.85 kb DNA fragment downstream of the genes. The two flanking fragments were introduced into pRFHU2 following the USER protocol described by  and the resulting plasmid was introduced into E. coli DH5α chemical competent cells. Proper fusions were checked by DNA sequencing and then the plasmid was transferred to A. tumefaciens AGL1 electrocompetent cells. Agrobacterium-mediated transformation of P. digitatum Pd1 was conducted basically as described in  with minor modifications. Briefly, equal volumes of a conidial suspension adjusted at 105 con/mL was mixed with A. tumefaciens AGL1 cells that had been grown in liquid induction medium (IM) supplemented with 200 μM acetosyringone until they reached an A600 of ~0.75. Aliquots of the mixture were spread onto 0.45 μm nitrocellulose membranes (Abet) that were layered on solid IM plates. After three days of incubation at 24°C the membranes were transferred to Petri plates containing PDA medium supplemented with 100 μg/mL hygromycin B (Invitrogen), for selection of fungal transformants, and 200 μg/mL cefotaxim (Calbiochem), for killing the bacterial cells. Putative transformants appeared after four days of incubation at 24°C. They were transferred to fresh plates containing PDA supplemented with hygromycin B and monosporic isolates were obtained. Conidia from monosporic isolates were inoculated in PDB medium containing 100 μg/mL hygromicin B and incubated at 24°C for two days with shaking. DNA was extracted from fungal mycelia following the procedure described in  and resuspended in 100 μL of TE. Transformants were analyzed by PCR with primers PD519F (5′-GCTTTCCCGCTTTAGTGTGG-3′) and PD519R (5′-ACAGCCGCAGCCTGTTATCT-3′). In the wild type strain these two primers amplify a fragment of 2.3 kb, whereas in gene knockout transformants this fragment is replaced by a 3.0 kb band. Ectopic transformants contain both bands.
Fungicide resistance assays were conducted in 96-well microtiter plates as described in . Wells contained 105 con/mL in PDB medium supplemented with fungicide. The fungicides azoxystrobin (Quadris; Syngenta Crop Protection), imazalil (Textar I; Tecnidex), myclobutanil (Thiocur 12; Rohmand Haas Italia SRL), prochloraz (Ascurit; Tecnidex) and trifloxystrobin (Flint; Bayer CropScience) were assayed at final concentrations of active ingredient of 0, 0.16, 0.31, 0.63, 1.25, 2.5, 5.0, and 10.0 μg/mL. Thiabendazole (Textar 60T; Tecnidex) was assayed at concentrations of 0, 0.63, 1.25, 2.5, 5.0, 10.0, 20.0 and 40.0 μg/mL. There were three replicate wells for each fungicide concentration. Plates were incubated at 24°C and growth, measured as A600, was followed for up to seven days. Minimum inhibitory concentration (MIC) values were considered as those in which no fungal growth was observed by seven days of incubation.
The authors wish to thank Heinz Himmelbauer and the CRG ultrasequencing facility for their technical support, Ana Izquierdo for her technical assistance and Andrew MacCabe for critical reading of the manuscript. TG research is funded in part by a grant from the Spanish Ministry of Economy and Innovation (BFU2009-09168). JFM research is funded in part by the Spanish Ministry of Science (BIO2009-12919). LGC research is funded in part by the Spanish Ministry of Economy and Innovation (AGL-2008-04828-C03-02 and AGL2011-30519-C03-01) and the Generalitat Valenciana (PROMETEO/2010/010 and ACOMP/2011/250). ARB acknowledges the support of the JAE-Doc program from CSIC and the European FEDER funds.
- FAO: FAO Training series. vol. 17. Prevention of post-harvest food losses: food, vegetables and root crops - a training manual. 1989, Rome: FAOGoogle Scholar
- Eckert J, Eaks I: Postharvest disorders and diseases of citrus fruits. The citrus industry. Edited by: Reuther W, Calavan E, Carman G. 1989, Berkeley: University of California Press, 179-260.Google Scholar
- Barkai-Golan R: Chemical control. Postharvest diseases of fruits and vegetables. Edited by: Barkai-Golan R. 2001, Amsterdam: Elsevier, 147-188.View ArticleGoogle Scholar
- Barkai-Golan R: Each fruit or vegetable and its characteristic pathogens. Postharvest diseases of fruits and vegetables. Edited by: Barkai-Golan R. 2001, Amsterdam: Elsevier, 25-32.View ArticleGoogle Scholar
- Gonzalez-Candelas L, Alamar S, Sanchez-Torres P, Zacarias L, Marcos JF: A transcriptomic approach highlights induction of secondary metabolism in citrus fruit in response to Penicillium digitatum infection. BMC Plant Biol. 2010, 10: 194-10.1186/1471-2229-10-194.PubMed CentralView ArticlePubMedGoogle Scholar
- Ballester AR, Lafuente MT, Forment J, Gadea J, De Vos RC, Bovy AG, Gonzalez-Candelas L: Transcriptomic profiling of citrus fruit peel tissues reveals fundamental effects of phenylpropanoids and ethylene on induced resistance. Mol Plant Pathol. 2011, 12 (9): 879-897. 10.1111/j.1364-3703.2011.00721.x.View ArticlePubMedGoogle Scholar
- Lopez-Garcia B, Gonzalez-Candelas L, Perez-Paya E, Marcos JF: Identification and characterization of a hexapeptide with activity against phytopathogenic fungi that cause postharvest decay in fruits. Mol Plant Microbe Interact. 2000, 13 (8): 837-846. 10.1094/MPMI.2000.13.8.837.View ArticlePubMedGoogle Scholar
- Sánchez-Torres P, Tuset J: Molecular insights into fungicide resistance in sensitive and resistant Penicillium digitatum strains infecting citrus. Postharvest Biol Tec. 2011, 59 (2): 159-165. 10.1016/j.postharvbio.2010.08.017.View ArticleGoogle Scholar
- Tarailo-Graovac M, Chen N: Using RepeatMasker to identify repetitive elements in genomic sequences. Curr Protoc Bioinformatics. 2009, 4: 4-10.Google Scholar
- Marcet-Houben M, Gabaldón T: The tree versus the forest: the fungal tree of life and the topological diversity within the yeast phylome. PLoS One. 2009, 4 (2): e4357-10.1371/journal.pone.0004357.PubMed CentralView ArticlePubMedGoogle Scholar
- Gabaldon T: Large-scale assignment of orthology: back to phylogenetics?. Genome Biol. 2008, 9 (10): 235-10.1186/gb-2008-9-10-235.PubMed CentralView ArticlePubMedGoogle Scholar
- Huerta-Cepas J, Marcet-Houben M, Pignatelli M, Moya A, Gabaldón T: The pea aphid phylome: a complete catalogue of evolutionary histories and arthropod orthology and paralogy relationships for Acyrthosiphon pisum genes. Insect Mol Biol. 2010, 19 (Suppl 2): 13-21.View ArticlePubMedGoogle Scholar
- Debuchy R, Berteaux-Lecellier V, Silar P: Mating systems and sexual morphogenesis in Ascomycetes. Cellular and molecular biology of filamentous fungi. Edited by: Borkovich K, DJ E. 2010, Washington: ASM Press, 501-531.View ArticleGoogle Scholar
- Frisvad JC, Samson RA: Polyphasic taxonomy of Penicillium subgenus Penicillium: A guide to identification of food and air-borne terverticillate Penicillia and their mycotoxins. Stud Mycol. 2004, 49: 1-173.Google Scholar
- Johnstone IL, McCabe PC, Greaves P, Gurr SJ, Cole GE, Brow MA, Unkles SE, Clutterbuck AJ, Kinghorn JR, Innis MA: Isolation and characterisation of the crnA-niiA-niaD gene cluster for nitrate assimilation in Aspergillus nidulans. Gene. 1990, 90 (2): 181-192. 10.1016/0378-1119(90)90178-T.View ArticlePubMedGoogle Scholar
- Burger G, Tilburn J, Scazzocchio C: Molecular cloning and functional characterization of the pathway-specific regulatory gene nirA, which controls nitrate assimilation in Aspergillus nidulans. Mol Cell Biol. 1991, 11 (2): 795-802.PubMed CentralView ArticlePubMedGoogle Scholar
- van den Berg MA, Albang R, Albermann K, Badger JH, Daran JM, Driessen AJ, Garcia-Estrada C, Fedorova ND, Harris DM, Heijne WH, et al: Genome sequencing and analysis of the filamentous fungus Penicillium chrysogenum. Nat Biotechnol. 2008, 26 (10): 1161-1168. 10.1038/nbt.1498.View ArticlePubMedGoogle Scholar
- Haas H, Schoeser M, Lesuisse E, Ernst JF, Parson W, Abt B, Winkelmann G, Oberegger H: Characterization of the Aspergillus nidulans transporters for the siderophores enterobactin and triacetylfusarinine C. Biochem J. 2003, 371 (Pt 2): 505-513.PubMed CentralView ArticlePubMedGoogle Scholar
- Oberegger H, Zadra I, Schoeser M, Abt B, Parson W, Haas H: Identification of members of the Aspergillus nidulans SREA regulon: genes involved in siderophore biosynthesis and utilization. Biochem Soc Trans. 2002, 30 (4): 781-783.View ArticlePubMedGoogle Scholar
- Koenraadt H, Somerville S, Jones A: Characterization of mutations in the b-tubulin gene of benomyl-resistant field strains of Venturia inaequalis and other pathogenic fungi. Mol Plant Pathol. 1992, 82: 1348-1354.Google Scholar
- Sarwar M, Akhtar M: Cloning of aminoglycoside phosphotransferase (APH) gene from antibiotic-producing strain of Bacillus circulans into a high-expression vector, pKK223–3. Purification, properties and location of the enzyme. Biochem J. 1990, 268 (3): 671-677.PubMed CentralView ArticlePubMedGoogle Scholar
- Sun X, Wang J, Feng D, Ma Z, Li H: PdCYP51B, a new putative sterol 14alpha-demethylase gene of Penicillium digitatum involved in resistance to imazalil and other fungicides inhibiting ergosterol synthesis. Appl Microbiol Biotechnol. 2011, 91 (4): 1107-1119. 10.1007/s00253-011-3355-7.View ArticlePubMedGoogle Scholar
- Gilles A, Meglecz E, Pech N, Ferreira S, Malausa T, Martin JF: Accuracy and quality assessment of 454 GS-FLX Titanium pyrosequencing. BMC Genomics. 2011, 12: 245-10.1186/1471-2164-12-245.PubMed CentralView ArticlePubMedGoogle Scholar
- Sun X, Li H, Yu D, Dijksterhuis J: Complete mitochondrial genome sequence of the phytopathogenic fungus Penicillium digitatum and comparative analysis of closely related species. FEMS Microbiol Lett. 2011, 323 (1): 29-34. 10.1111/j.1574-6968.2011.02358.x.View ArticlePubMedGoogle Scholar
- Webber J: History and development of the citrus industry. 1967, Riverside: University of California Division of Agricultural SciencesGoogle Scholar
- Huerta-Cepas J, Capella-Gutierrez S, Pryszcz LP, Denisov I, Kormes D, Marcet-Houben M, Gabaldon T: PhylomeDB v3.0: an expanding repository of genome-wide collections of trees, alignments and phylogeny-based orthology and paralogy predictions. Nucleic Acids Res. 2011, 39 (Database issue): D556-D560.PubMed CentralView ArticlePubMedGoogle Scholar
- Dimmer EC, Huntley RP, Barrell DG, Binns D, Draghici S, Camon EB, Hubank M, Talmud PJ, Apweiler R, Lovering RC: The gene ontology - providing a functional role in proteomic studies. Proteomics. 2008, 8 (Suppl 23-24): 2-11.Google Scholar
- Marcet-Houben M, Gabaldon T: Acquisition of prokaryotic genes by fungal genomes. Trends Genet. 2010, 26 (1): 5-8. 10.1016/j.tig.2009.11.007.View ArticlePubMedGoogle Scholar
- Rose MS, Yun SH, Asvarak T, Lu SW, Yoder OC, Turgeon BG: A decarboxylase encoded at the Cochliobolus heterostrophus translocation-associated Tox1B locus is required for polyketide (T-toxin) biosynthesis and high virulence on T-cytoplasm maize. Mol Plant Microbe Interact. 2002, 15 (9): 883-893. 10.1094/MPMI.2002.15.9.883.View ArticlePubMedGoogle Scholar
- Keller NP, Bennett J, Turner G: Secondary metabolism: then, now and tomorrow. Fungal Genet Biol. 2010, 48 (1): 1-3.View ArticlePubMedGoogle Scholar
- Frisvad JC, Smedsgaard J, Larsen TO, Samson RA: Mycotoxins, drugs and other extrolites produced by species in Penicillium subgenus Penicillium. Stud Mycol. 2004, 49: 201-241.Google Scholar
- Leitao AL: Potential of Penicillium species in the bioremediation field. Int J Environ Res Public Health. 2009, 6 (4): 1393-1417. 10.3390/ijerph6041393.PubMed CentralView ArticlePubMedGoogle Scholar
- Chang PK, Horn BW, Dorner JW: Clustered genes involved in cyclopiazonic acid production are next to the aflatoxin biosynthesis gene cluster in Aspergillus flavus. Fungal Genet Biol. 2009, 46 (2): 176-182. 10.1016/j.fgb.2008.11.002.View ArticlePubMedGoogle Scholar
- Artigot MP, Loiseau N, Laffitte J, Mas-Reguieg L, Tadrist S, Oswald IP, Puel O: Molecular cloning and functional characterization of two CYP619 cytochrome P450s involved in biosynthesis of patulin in Aspergillus clavatus. Microbiology. 2009, 155 (Pt 5): 1738-1747.PubMed CentralView ArticlePubMedGoogle Scholar
- Moake MM, Padilla-Zakour OI, Worobo RW: Comprehensive review of patulin control methods in foods. Compr Rev Food Sci F. 2005, 4 (1): 8-21. 10.1111/j.1541-4337.2005.tb00068.x.View ArticleGoogle Scholar
- Gao X, Chooi YH, Ames BD, Wang P, Walsh CT, Tang Y: Fungal indole alkaloid biosynthesis: genetic and biochemical investigation of the tryptoquialanine pathway in Penicillium aethiopicum. J Am Chem Soc. 2011, 133 (8): 2729-2741. 10.1021/ja1101085.PubMed CentralView ArticlePubMedGoogle Scholar
- Wosten HA: Hydrophobins: multipurpose proteins. Annu Rev Microbiol. 2001, 55: 625-646. 10.1146/annurev.micro.55.1.625.View ArticlePubMedGoogle Scholar
- Stergiopoulos I, de Wit PJ: Fungal effector proteins. Annu Rev Phytopathol. 2009, 47: 233-263. 10.1146/annurev.phyto.112408.132637.View ArticlePubMedGoogle Scholar
- de Jonge R, van Esse HP, Kombrink A, Shinya T, Desaki Y, Bours R, van der Krol S, Shibuya N, Joosten MH, Thomma BP: Conserved fungal LysM effector Ecp6 prevents chitin-triggered immunity in plants. Science. 2010, 329 (5994): 953-955. 10.1126/science.1190859.View ArticlePubMedGoogle Scholar
- Marx F: Small, basic antifungal proteins secreted from filamentous ascomycetes: a comparative study regarding expression, structure, function and potential application. Appl Microbiol Biotechnol. 2004, 65 (2): 133-142.View ArticlePubMedGoogle Scholar
- Mygind PH, Fischer RL, Schnorr KM, Hansen MT, Sonksen CP, Ludvigsen S, Raventos D, Buskov S, Christensen B, De Maria L, et al: Plectasin is a peptide antibiotic with therapeutic potential from a saprophytic fungus. Nature. 2005, 437 (7061): 975-980. 10.1038/nature04051.View ArticlePubMedGoogle Scholar
- de Jonge R, Thomma BP: Fungal LysM effectors: extinguishers of host immunity?. Trends Microbiol. 2009, 17 (4): 151-157. 10.1016/j.tim.2009.01.002.View ArticlePubMedGoogle Scholar
- Hancock RE, Sahl HG: Antimicrobial and host-defense peptides as new anti-infective therapeutic strategies. Nat Biotechnol. 2006, 24 (12): 1551-1557. 10.1038/nbt1267.View ArticlePubMedGoogle Scholar
- Yount NY, Yeaman MR: Multidimensional signatures in antimicrobial peptides. Proc Natl Acad Sci U S A. 2004, 101 (19): 7363-7368. 10.1073/pnas.0401567101.PubMed CentralView ArticlePubMedGoogle Scholar
- Marx F, Binder U, Leiter E, Pocsi I: The Penicillium chrysogenum antifungal protein PAF, a promising tool for the development of new antifungal therapies and fungal cell biology studies. Cell Mol Life Sci. 2008, 65 (3): 445-454. 10.1007/s00018-007-7364-8.View ArticlePubMedGoogle Scholar
- Rodriguez-Martin A, Acosta R, Liddell S, Nunez F, Benito MJ, Asensio MA: Characterization of the novel antifungal chitosanase PgChP and the encoding gene from Penicillium chrysogenum. Appl Microbiol Biotechnol. 2010, 88 (2): 519-528. 10.1007/s00253-010-2767-0.View ArticlePubMedGoogle Scholar
- Silverstein KA, Moskal WA, Wu HC, Underwood BA, Graham MA, Town CD, VandenBosch KA: Small cysteine-rich peptides resembling antimicrobial peptides have been under-predicted in plants. Plant J. 2007, 51 (2): 262-280. 10.1111/j.1365-313X.2007.03136.x.View ArticlePubMedGoogle 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 (Database issue): D572-D576.PubMed CentralPubMedGoogle Scholar
- Cantarel BL, Coutinho PM, Rancurel C, Bernard T, Lombard V, Henrissat B: The Carbohydrate-Active EnZymes database (CAZy): an expert resource for Glycogenomics. Nucleic Acids Res. 2009, 37 (Database issue): D233-D238.PubMed CentralView ArticlePubMedGoogle Scholar
- Amselem J, Cuomo CA, van Kan JA, Viaud M, Benito EP, Couloux A, Coutinho PM, de Vries RP, Dyer PS, Fillinger S, et al: Genomic analysis of the necrotrophic fungal pathogens Sclerotinia sclerotiorum and Botrytis cinerea. PLoS Genet. 2011, 7 (8): e1002230-10.1371/journal.pgen.1002230.PubMed CentralView ArticlePubMedGoogle Scholar
- Gonzalez-Fernandez R, Jorrin-Novo JV: Contribution of proteomics to the study of plant pathogenic fungi. J Proteome Res. 2012, 11 (1): 3-16. 10.1021/pr200873p.View ArticlePubMedGoogle Scholar
- Schmidt SM, Panstruga R: Pathogenomics of fungal plant parasites: what have we learnt about pathogenesis?. Curr Opin Plant Biol. 2011, 14 (4): 392-399. 10.1016/j.pbi.2011.03.006.View ArticlePubMedGoogle Scholar
- Idnurm A, Howlett BJ: Pathogenicity genes of phytopathogenic fungi. Mol Plant Pathol. 2001, 2 (4): 241-255. 10.1046/j.1464-6722.2001.00070.x.View ArticlePubMedGoogle Scholar
- Ballester AR, Lafuente MT: Spatial study of antioxidant enzymes, peroxidase and phenylalanine ammonia-lyase in the citrus fruit-Penicillium digitatum interaction. Postharvest Biol Tec. 2006, 39 (2): 115-124. 10.1016/j.postharvbio.2005.10.002.View ArticleGoogle Scholar
- Moller EM, Bahnweg G, Sandermann H, Geiger HH: A simple and efficient protocol for isolation of high molecular weight DNA from filamentous fungi, fruit bodies, and infected plant tissues. Nucleic Acids Res. 1992, 20 (22): 6115-6116. 10.1093/nar/20.22.6115.PubMed CentralView ArticlePubMedGoogle Scholar
- Li R, Li Y, Kristiansen K, Wang J: SOAP: short oligonucleotide alignment program. Bioinformatics. 2008, 24 (5): 713-714. 10.1093/bioinformatics/btn025.View ArticlePubMedGoogle Scholar
- Li R, Zhu H, Ruan J, Qian W, Fang X, Shi Z, Li Y, Li S, Shan G, Kristiansen K, et al: De novo assembly of human genomes with massively parallel short read sequencing. Genome Res. 2010, 20 (2): 265-272. 10.1101/gr.097261.109.PubMed CentralView ArticlePubMedGoogle Scholar
- Richter DC, Schuster SC, Huson DH: OSLay: optimal syntenic layout of unfinished assemblies. Bioinformatics. 2007, 23 (13): 1573-1579. 10.1093/bioinformatics/btm153.View ArticlePubMedGoogle Scholar
- Allen JE, Salzberg SL: JIGSAW: integration of multiple sources of evidence for gene prediction. Bioinformatics. 2005, 21 (18): 3596-3603. 10.1093/bioinformatics/bti609.View ArticlePubMedGoogle Scholar
- Stanke M, Waack S: Gene prediction with a hidden Markov model and a new intron submodel. Bioinformatics. 2003, 19 (Suppl 2): ii215-ii225. 10.1093/bioinformatics/btg1080.View ArticlePubMedGoogle Scholar
- Guigo R, Knudsen S, Drake N, Smith T: Prediction of gene structure. J Mol Biol. 1992, 226 (1): 141-157. 10.1016/0022-2836(92)90130-C.View ArticlePubMedGoogle Scholar
- Korf I: Gene finding in novel genomes. BMC Bioinformatics. 2004, 5: 59-10.1186/1471-2105-5-59.PubMed CentralView ArticlePubMedGoogle Scholar
- Wu TD, Watanabe CK: GMAP: a genomic mapping and alignment program for mRNA and EST sequences. Bioinformatics. 2005, 21 (9): 1859-1875. 10.1093/bioinformatics/bti310.View ArticlePubMedGoogle Scholar
- Birney E, Clamp M, Durbin R: GeneWise and Genomewise. Genome Res. 2004, 14 (5): 988-995. 10.1101/gr.1865504.PubMed CentralView ArticlePubMedGoogle Scholar
- Punta M, Coggill PC, Eberhardt RY, Mistry J, Tate J, Boursnell C, Pang N, Forslund K, Ceric G, Clements J, et al: The Pfam protein families database. Nucleic Acids Res. 2011, 40 (Database issue): D290-D301.PubMed CentralPubMedGoogle Scholar
- Kanehisa M, Goto S, Furumichi M, Tanabe M, Hirakawa M: KEGG for representation and analysis of molecular networks involving diseases and drugs. Nucleic Acids Res. 2010, 38 (Database issue): D355-D360.PubMed CentralView ArticlePubMedGoogle Scholar
- Horton P, Park KJ, Obayashi T, Fujita N, Harada H, Adams-Collier CJ, Nakai K: WoLF PSORT: protein localization predictor. Nucleic Acids Res. 2007, 35 (Web Server issue): W585-W587.PubMed CentralView ArticlePubMedGoogle Scholar
- Petersen TN, Brunak S, von Heijne G, von Nielsen H: SignalP 4.0: discriminating signal peptides from transmembrane regions. Nat Methods. 2011, 8 (10): 785-786. 10.1038/nmeth.1701.View ArticlePubMedGoogle Scholar
- Krogh A, Larsson B, von Heijne G, Sonnhammer EL: Predicting transmembrane protein topology with a hidden Markov model: application to complete genomes. J Mol Biol. 2001, 305 (3): 567-580. 10.1006/jmbi.2000.4315.View ArticlePubMedGoogle Scholar
- Graham MA, Silverstein KA, Cannon SB, VandenBosch KA: Computational identification and characterization of novel genes from legumes. Plant Physiol. 2004, 135 (3): 1179-1197. 10.1104/pp.104.037531.PubMed CentralView ArticlePubMedGoogle Scholar
- Eddy SR: Accelerated Profile HMM Searches. PLoS Comput Biol. 2011, 7 (10): e1002195-10.1371/journal.pcbi.1002195.PubMed CentralView ArticlePubMedGoogle Scholar
- Edgar RC: MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32 (5): 1792-1797. 10.1093/nar/gkh340.PubMed CentralView ArticlePubMedGoogle Scholar
- Katoh K, Kuma K, Toh H, Miyata T: MAFFT version 5: improvement in accuracy of multiple sequence alignment. Nucleic Acids Res. 2005, 33 (2): 511-518. 10.1093/nar/gki198.PubMed CentralView ArticlePubMedGoogle Scholar
- Subramanian AR, Kaufmann M, Morgenstern B: DIALIGN-TX: greedy and progressive approaches for segment-based multiple sequence alignment. Algorithms Mol Biol. 2008, 3: 6-10.1186/1748-7188-3-6.PubMed CentralView ArticlePubMedGoogle Scholar
- Landan G, Graur D: Heads or tails: a simple reliability check for multiple sequence alignments. Mol Biol Evol. 2007, 24 (6): 1380-1383. 10.1093/molbev/msm060.View ArticlePubMedGoogle Scholar
- Wallace IM, O’Sullivan O, Higgins DG, Notredame C: M-Coffee: combining multiple sequence alignment methods with T-Coffee. Nucleic Acids Res. 2006, 34 (6): 1692-1699. 10.1093/nar/gkl091.PubMed CentralView ArticlePubMedGoogle Scholar
- Capella-Gutierrez S, Silla-Martinez JM, Gabaldón T: trimAl: a tool for automated alignment trimming in large-scale phylogenetic analyses. Bioinformatics. 2009, 25 (15): 1972-1973. 10.1093/bioinformatics/btp348.PubMed CentralView ArticlePubMedGoogle Scholar
- Gascuel O: BIONJ: an improved version of the NJ algorithm based on a simple model of sequence data. Mol Biol Evol. 1997, 14 (7): 685-695. 10.1093/oxfordjournals.molbev.a025808.View ArticlePubMedGoogle Scholar
- Guindon S, Dufayard JF, Lefort V, Anisimova M, Hordijk W, Gascuel O: New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Syst Biol. 2010, 59 (3): 307-321. 10.1093/sysbio/syq010.View ArticlePubMedGoogle Scholar
- Akaike H: Information theory and extension of the maximum likelihood principle. Proceedings of the 2nd international symposium on information theory. 1973, Hungary: Budapest, 267-281.Google Scholar
- Huerta-Cepas J, Dopazo J, Gabaldón T: ETE: a python environment for tree exploration. BMC Bioinformatics. 2010, 11: 24-10.1186/1471-2105-11-24.PubMed CentralView ArticlePubMedGoogle Scholar
- Stamatakis A, Ludwig T, Meier H: RAxML-III: a fast program for maximum likelihood-based inference of large phylogenetic trees. Bioinformatics. 2005, 21 (4): 456-463. 10.1093/bioinformatics/bti191.View ArticlePubMedGoogle Scholar
- Le SQ, Gascuel O: An improved general amino acid replacement matrix. Mol Biol Evol. 2008, 25 (7): 1307-1320. 10.1093/molbev/msn067.View ArticlePubMedGoogle Scholar
- Retief JD: Phylogenetic analysis using PHYLIP. Methods Mol Biol. 2000, 132: 243-258.PubMedGoogle Scholar
- Wehe A, Bansal MS, Burleigh JG, Eulenstein O: DupTree: a program for large-scale phylogenetic analyses using gene tree parsimony. Bioinformatics. 2008, 24 (13): 1540-1541. 10.1093/bioinformatics/btn230.View ArticlePubMedGoogle Scholar
- Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R: The sequence alignment/Map format and SAMtools. Bioinformatics. 2009, 25 (16): 2078-2079. 10.1093/bioinformatics/btp352.PubMed CentralView ArticlePubMedGoogle Scholar
- Frandsen RJ, Andersson JA, Kristensen MB, Giese H: Efficient four fragment cloning for the construction of vectors for targeted gene replacement in filamentous fungi. BMC Mol Biol. 2008, 9: 70-10.1186/1471-2199-9-70.PubMed CentralView ArticlePubMedGoogle Scholar
- Michielse CB, Hooykaas PJ, van den Hondel CA, Ram AF: Agrobacterium-mediated transformation of the filamentous fungus Aspergillus awamori. Nat Protoc. 2008, 3 (10): 1671-1678. 10.1038/nprot.2008.154.View ArticlePubMedGoogle Scholar
- Raeder U, Broda P: Rapid preparation of DNA from filamentous fungi. Lett Appl Microbiol. 1985, 1: 17-20. 10.1111/j.1472-765X.1985.tb01479.x.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.