- Research article
- Open Access
De-novo assembly and characterization of the transcriptome of Metschnikowia fructicola reveals differences in gene expression following interaction with Penicillium digitatumand grapefruit peel
BMC Genomics volume 14, Article number: 168 (2013)
The yeast Metschnikowia fructicola is an antagonist with biological control activity against postharvest diseases of several fruits. We performed a transcriptome analysis, using RNA-Seq technology, to examine the response of M. fructicola with citrus fruit and with the postharvest pathogen, Penicillium digitatum.
More than 26 million sequencing reads were assembled into 9,674 unigenes. Approximately 50% of the unigenes could be annotated based on homology matches in the NCBI database. Based on homology, sequences were annotated with a gene description, gene ontology (GO term), and clustered into functional groups. An analysis of differential expression when the yeast was interacting with the fruit vs. the pathogen revealed more than 250 genes with specific expression responses. In the antagonist-pathogen interaction, genes related to transmembrane, multidrug transport and to amino acid metabolism were induced. In the antagonist-fruit interaction, expression of genes involved in oxidative stress, iron homeostasis, zinc homeostasis, and lipid metabolism were induced. Patterns of gene expression in the two interactions were examined at the individual transcript level by quantitative real-time PCR analysis (RT-qPCR).
This study provides new insight into the biology of the tritrophic interactions that occur in a biocontrol system such as the use of the yeast, M. fructicola for the control of green mold on citrus caused by P. digitatum.
Considerable research has been devoted in the past two decades to investigate the use of naturally occurring yeasts for managing postharvest diseases on a variety of fruits and vegetables [1, 2]. At present, only two commercial products are available for postharvest use, each having a very small market share of the collective technologies used to manage postharvest diseases. Biosave (Pseudomonas syringae Van Hall) used for the control of pome fruit, sweet potato, and potato diseases , and “Shemer” (Metschnikowia fructicola) used commercially for the control of sweet potato and carrot storage diseases , strawberry , apple  and citrus fruit .
Knowledge about mechanisms of action of yeast biocontrol agents used to manage postharvest diseases is relatively superficial because appropriate methods to study microbial interactions in wounds of fruit has been lacking. Various mechanisms, however, have been described, including antibiosis, production of lytic enzymes, parasitism, induction of host resistance, and competition for limiting nutrients and space [2, 8–10]. Results described in these studies are mostly correlative and provide indirect evidence for possible involvement of one or more of these components in the mechanism of action. This may be due to the lack of an understanding of the key biochemical and molecular processes occurring within tritrophic interactions (host-pathogen-biocontrol agent) that define the effectiveness of a biocontrol system.
A characteristic response of fruit to pathogen attack or application of a biocontrol agent is the activation of defence responses, including production of inhibitors of cell wall-degrading enzymes of the pathogen, activity of antifungal compounds (such as phenolic compounds and phytoalexins), active oxygen species, and reinforcement of the cell wall of the host . In relation to yeast biocontrol agents, induction of resistance in fruit tissue has been reported to involve, to a greater or lesser extent, the active elicitation of all of the above responses. Hershkovitz et al.  recently reported the first molecular study of the response of citrus fruit to the application of the yeast biocontrol agent, M. fructicola. Results indicated that significant changes in gene expression occur in grapefruit wounded tissue in response to the application of M. fructicola including changes in a variety of genes involved in response to biotic and abiotic stresses, signaling, defence and secondary metabolism. These findings imply that complex biochemical and molecular processes are involved in the reaction of fruit host tissue to the introduction of yeast cells, which may have the potential to influence the efficacy of the biocontrol agent.
Castoria et al.  indicated that the ability to tolerate high levels of reactive oxygen species (ROS) produced in fruit tissue in response to wounding is an essential characteristic of effective yeast antagonists. Macarisin et al.  demonstrated that the antagonistic yeasts, Candida oleophila and M. fructicola, actually induce defence-related oxidative responses in intact fruit exocarp tissues and within wounds of apple and citrus, generating greater levels of super oxide anion (O2−) and hydrogen peroxide (H2O2). Torres et al.  reported that, in addition to nutrient competition, the biocontrol activity of Pantoea agglomerans could be attributed to the ability of the antagonist to induce an oxidative response in orange fruit peel tissue by triggering H2O2 production, superoxide dismutase, and catalase activities. In another study, transient exposure of the yeast M. fructicola to a sublethal heat or oxidative stress increased overall abiotic stress tolerance, providing cross-protection to a range of other more severe abiotic stresses . The study showed that higher levels of trehalose, induced by exposure to a sub-lethal abiotic stress, was associated with an increase in ROS scavenging, general stress tolerance, and more rapid growth in apple wounds, collectively resulting in improved biocontrol activity.
Larralde-Corona et al.  investigated the effect of the exposure of the yeast Pichia guilliermondii to Penicillium digitatum on yeast gene expression by analyzing expressed sequencing tags (ESTs) obtained by subtractive suppressive hybridization (SSH). The yeast was exposed to starvation by carbon source, actively growing mycelium of P. digitatum separated from yeast cells with membrane and induction by fungal cell walls. Only one EST, associated with energy metabolism, was obtained in response to the limited carbon source, seven ESTs in the separated membrane system, all relating to metabolic networks such as: energy, nitrogen, cell cycle, ABC transporters, response to stress and one unknown sequence. The exposure of the yeast cells to fungal cell walls produced the highest number of ESTs, a total of 22, including all the metabolic networks mentioned above plus ESTs associated with signal transduction.
Compared to the information available about the response of host tissue to the introduction of yeast biocontrol agents, little is known about global molecular changes taking place in yeast cells in response to an interaction with either fruit tissue or a pathogen's mycelia. In the current study, we characterized transcriptional changes in M. fructicola cells following their interaction with grapefruit peel tissue and mycelium of P. digitatum using next generation sequencing technology for RNA (RNA-Seq) on the Illumina platform. Along with de-novo assembly and characterization of the transcriptome, analysis of global patterns of gene expression and functional categorization were performed.
Illumina sequencing, assembly, and annotation
Transcriptome analysis of M. fructicola utilizing Illumina sequencing generated a total of 26,829,178 reads. The raw reads were assembled into 15,803 contigs by Trinity software with an average size of 490.2 bp and median of 349 bp. The sequences were determined to represent 9684 unigenes with an average length of 546.2 bp and median of 385 bp. The size distribution is shown in Figure 1. To annotate these genes, reference sequences were searched using BLASTX within the non-redundant (nr) NCBI nucleotide database using a cut-off E-value of 10-3 (the sequences of the unigenes are found in Additional file 1 in fasta format). A total of 2452 genes (25.3% of all gene sequences) could not be annotated based on matches in the NCBI database (the full annotations of M. fructicola unigenes is found in Additional file 2). A high level (75%) of transcripts from the individual libraries could be mapped to the assembled transcriptome using bowtie software. These transcripts were blasted against the M. fructicola draft genome (available under genbank accession number ANFW01000000) with cutoff e-value of 1e-50 and identity >85%. According to this analysis 9631 unigenes (99.4%) could be aligned to the genome, indicating high accuracy of our assembly (mapping information is found in Additional file 3). The remaining 53 unigenes could be a result of erroneous assembly.
As shown in Figure 2, 42% of sequences had significant homology matches (< 1.0E10− 50) in the NCBI database, while 58% of the sequences had matches ranging between 1.0E10-3 and 1.0E10-50. The similarity distribution showed that 11% of the sequences had > 90% homology, followed by 86% of the sequences ranging between 50 to 90% homology. Only 3% of the sequences had homology lower than 50% (Figure 2B). Species specific distribution indicated that 69% of the sequences matched gene sequences of the yeast Clavispora lustiniae, followed by other yeasts Debar-yomyces hansenii (8%), Scheffersomyces stiptis (5%), Meyerozyma guilliermondii (4%), Spathaspora passa-lidarum, Millerozyma farinose and Candida tenuis (2% each) (Figure 2C). The phylogenetic analyses presented in Figure 2D shown close sequence similarity between M. fructicola and C. lustiniae (Figure 2D).
Classification of gene ontology (GO)
To classify the functions of M. fructicola unigenes, we searched the Gene Ontology (GO) database and used the WEGO tool (web gene ontology annotation Plotting) for plotting GO annotation results (http://wego.genomics.org.cn). A total of 4126 sequences out of 9674 (57.6% of annotated sequences) were categorized into 43 functional groups (Figure 3). In each of the three main categories (biological process, cellular component and molecular function) of the GO classification the 'Cell and Cell part', 'Antioxidant and Binding' and 'Metabolic and cellular processes' terms were the most dominant. In contrast, 'metallochapreon', and 'rhythmic process' were least abundant.
RNA-seq technology is a sensitive and reliable tool to study transcriptional changes. Using this technique we identified 153 genes differentially expressed (DE) with at least a 1.8-fold change in yeast cells in response to yeast 24 h after interaction with P. digitatum (antagonist-pathogen: Mf-Pdig). Among these, 105 were upregulated and 45 were downregulated. In the Mf-fruit interaction, 143 genes had differential expression where 108 were upregulated and 36 were downregulated, while 42 genes were regulated in both interactions. DE genes in either Mf-Pdig and Mf-fruit interactions or specific to one of these interactions are shown in a two-dimensional hierarchical cluster heat map and Venn diagram (Figure 4) Additional file 4.
Functional analysis of differentially expressed genes
All differentially regulated genes were organized in functional categories based on their gene ontology (GO) terms (http://www.blast2go.com/b2glaunch) and were members of at least one of the three principal GO branches: biological process, molecular function, and cellular component (Figure 5). All 16 functional groups were found to be over-expressed in Mf-Pdig interaction compared to all transcribed genes of M. fructicola controls (Figure 5A). We focused our analyses on “biological process” categories. Those categories were organized in three groups: response to stress and stimulus transport (Figure 6A), transport (Figure 6B) and metabolic processes (Figure 6C).
The induced genes in the category of response to stress in Mf-Pdig interaction are mainly involved in multidrug resistance and resistance to singlet oxygen species (YHK, YOR1, SNQ2 and STB5) and in oxidation-reduction activity in the Mf-fruit interaction (MET10, MET5, RNR2, HEM13, JLP1, FET3, AFG1 and PRX1) (Figure 6A; Additional file 5). The genes encoding enzymes involved in response to heat stress were, such as heat shock proteins HSP78, HSP40, HSP104 and SGT2, were repressed, in both of interactions.
Transport genes involved in iron transport and uptake (FTR1, FET3 and SIT2), zinc transport (ZTR2 and ZTR3) and SulP anion (SUL1) were upregulated in the Mf-fruit interaction while in the Mf-Pdig interaction, ammonium permease (MEP2), inorganic phosphate (Pi) transporter and low-affinity manganese transporter (PHO84) were induced (Figure 6B; Additional file 6).
A glucose sensor (SNF3) and two glucose transporter (HXT 6 and 4) genes were induced in the Mf-Pdig interaction. A group of genes associated with transmembrane transport, such as general amino acid permease (GAP1, myo-inositol transporter (ITP2), oligopeptide transporter (OPT2), a permease (GIT1), which mediates uptake of glycerophosphoinositol and glycerophosphocholine as sources for inositol and phosphate, an a sugar transporter (HXT3), were induced in both interactions.
Functional analysis of the metabolic processes category included lipid, cofactor and cellular amino acid metabolic process annotations. Many genes involved in the ergosterol-biosynthesis pathway (ERG1, ERG 11, ERG5, PSD and DAP1) were up-regulated in the Mf-fruit interaction (Additional file 7).
Many genes involved in amino acid metabolism (AGP3, GAP1, UGA1, CPA2, MAE1, CAR2 and FMS1) were up-regulated in the Mf-Pdig interaction. In the Mf-fruit interaction, two groups of genes involved in cellular amino acid metabolism were found to be enriched: one up-regulated (CAP1, CHI, LAP3, ARO1, LIA1 and PDC1) and the second, down-regulated (MIS1, OXP1, ALD3, CAR1 and GCV1) (Additional file 7).
Real time PCR (qRT-PCR) analysis
To confirm the expression profiles obtained from the RNA-seq expression data, RT-qPCR analysis was carried out for 8 genes selected on the basis of their biological significance: genes involved in the degradation of fungal cell wall, such as glucanse (GLU) and chitinase (CHI), were selected as were genes related to response to stresses, such as an ABC transporter, a transcription factor (STB5), a heat shock protein (HSP-40), superoxide dismutase (SOD1), a serine threonine (Ser-Thr) kinase, and a signal transducing MEK kinase. Although minor variation of transcription levels were obtained, 7 of the tested 8 genes correlated well with the fold-change, either up or down- obtained in the RNA-Seq analysis (analyzed the 24 h time point) (Figure 7). Time points included in the RT-qPCR analysis were 5, 12, 24 and 48 h after interaction of the yeast cells with either grapefruit peel or P. digitatum mycelia (Figure 7). Expression level of the SOD1 gene in RT-qPCR was two-fold higher compare to the RNA-seq data set. Transcript levels of GLU were induced in M. fructicola cells after 24 and 48 h in the Mf-Pdig interaction (Figure 8A). CHI levels were up-regulated in both Mf-Pdig and Mf-fruit interactions after 24 h and 48 h (Figure 8B). STB5 expression level displayed transient induction in both Mf-fruit and Mf-Pdig interactions at 24 h (Figure 8C), while the Ser-Thr kinase was transiently repressed in both Mf-fruit and Mf-Pdig interactions at 24 h (Figure 8H). The ABC transporter was induced in the Mf-Pdig interaction at 24 and 48 h, while in the Mf-fruit interaction it was induced only at 24 h (Figure 8D). Expression levels of SOD1 and MEKK were repressed in the Mf-Pdig interaction at all-time points tested compared to control cells grown in NYDB media. In the Mf-fruit interaction these genes showed a transient induction at 12 and 24 h (Figure 8E, G). HSP transcript levels increased following contact with either the fruit peel or P. digitatum mycelia but were down-regulated in comparison to the control, M. fructicola cells grown in NYDB (Figure 8F).
Information about the interactions between yeast biocontrol agents, pathogens, and plant hosts (mainly harvested fruit) is needed to better understand biochemical and molecular processes involved in the biocontrol system [10, 15]. In a previous study, we identified global changes in gene expression in grapefruit peel tissue following the application of the yeast Metschnikowia fructicola. Marked changes in the profiles of genes expressed in fruit tissue were observed in response to application of yeast into wounded fruit tissue. In the present study, the main objective was to compare changes in gene expression in M. fructicola cells following the interaction of the yeast with either grapefruit peel tissue or P. digitatum mycelia.
Although the genome of M. fructicola has not been sequenced, de-novo assembly of the transcriptome of M. fructicola resulted in the identification of a total of 9674 unigenes, half of which could be annotated based on homology to genes in the NCBI database (Figures 1 and 2). As shown in Figure 2C, 69% of the unigene sequences identified in M. fructicola showed highest homology to Clavispora lustiniae genes. Thus the RNA-Seq based transcriptome analysis of M. fructicola generated a large number of newly identified yeast genes, providing a substantial contribution to the existing sequence database.
Using RNA-seq, we compared transcriptomes of M. fructicola after the yeast were exposed to either grapefruit peel tissue (Mf-Pdig) or after coming into direct contact with the fungal mycelia of P. digitatum (Mf-Pdig). Our analysis identified more than 250 genes that were differentially expressed and that may potentially affect the biocontrol activity of M. fructicola (Figure 4). The number of genes regulated in Mf-fruit interaction approximately same that regulated in Mf-Pdig and the overlap between the two interaction is only 16.5% of genes co-regulated.
Stress related genes
Activation of stress-related genes in yeast biocontrol agents induced in response to the pathogen or plant host has been observed in several biocontrol systems [14, 16, 17]. In the current study, several genes involved in oxidation – reduction activity (GO: 0055114) and implicated in oxidative stress response, such as peroxidases and reductases, were up-regulated in the Mf-fruit interaction (Additional file 5). These results could be associated with the generation and detoxification of reactive oxygen species (ROS).
M. fructicola cells produce ROS as in response to contact with fruit tissue [7, 11]. These studies indicate that ROS plays a major role in the initial response of yeast cells to host tissue and act as a signal, inducing an oxidative burst in the fruit tissue which subsequently leads to the induction of resistance mechanisms in the host tissue.
The time course study of superoxide dismutase (SOD1) gene expression showed a transient induction in Mf-fruit samples at 12 and 24 h (Figure 8E). These results provide further evidence that yeast cells are undergoing oxidative stress while in contact with fruit tissue and suggest the possibility of an adaptive response. In this regard, Liu et al.  demonstrated that a mild heat-shock pretreatment (30 min at 40°C) improved the tolerance of M. fructicola to subsequent high temperature (45°C, 20–30 min) and oxidative stress (0.4 mol L-1) and increased the biocontrol efficacy of the yeast. On the other hand, lower or unaltered expression levels of SOD (Figure 8E) and low levels of ROS in Mf-Pdig interaction (data not shown) indicate that in this type of interaction, yeast cells are not under oxidative stress and hence activation of protective systems against superoxide radicals is not necessary.
High levels of oxidative stress in microorganisms is accompanied by the induction of heat-shock proteins (HSP) in an attempt to ameliorate damage to proteins caused by ROS . Induction of HSPs in response to biotic and abiotic stimuli is well documented where they assist in the folding and unfolding, and general stability, of various proteins. Production of HSPs in microorganisms is associated with various abiotic stresses, including oxidative stress . In our present study, genes encoding HSP78 and HSP104 were down regulated in both types of interactions compared to the control (yeast grown in NYDB) (Additional file 5). However, the time course study of HSP78 expression by RT-qPCR revealed significantly higher levels of transcript after 12 h, 24 and 48h in both interactions compared to its level at 5 h (Figure 8F). This result is surprising since it was anticipated that HSPs would be induced in yeast cells that are in contact with fruit peel tissue or with fungal mycelia considering the stressful environment it represents, such as desiccation due to the dryness of the fruit surface, ROS, and phytoalexins.
Response to chemical stress
Genes belonged to plasma membrane ATP-binding cassette (ABC) transporters, such as YHK8, YOR1 and SNQ2, were upregulated in the MF-Pdig interaction (Additional file 6). YOR1 and SNG2 genes were induced in S. cerevisiae following exposure to the monoterpene thymol . The function of these transporters in S. cerevisiae and their possible role in yeasts tolerance to monoterpene have yet to be elucidated . In Pichia guilliermondii, two ABC transporter related proteins involved in the secretion of toxic compounds were induced during its interaction with P. digitatum. Many of the ABC transporters are reported to be tightly regulated by transcription factors . In this regard, the current study showed that STB5, a transcription factor involved in regulating multidrug resistance and oxidative stress response, was transiently upregulated in the Mf-Pdig interaction at 24 h (Figure 8C). Induction of ABC transporters in M. fructicola cells may be needed for the yeast to tolerate exposure to fruit surface wounds or to mycelia where many potentially toxic substances (e.g. ROS, phytoalexines, toxins) may be present. This explanation may also be relevant to the ability of yeasts to tolerate high levels of ROS in apple fruit wounds which was suggested to be an essential characteristic of a successful biocontrol agent .
Lytic enzymes, such as glucanases and chitinases, involved in degradation of fungal cell walls have been considered to play an important role in the biological control of postharvest pathogens [23–25]. Jijakli and Lepoivre  indicated that P. guillermondii, a yeast biocontrol agent, exhibits a high level of β-1,3-glucanase activity that could function to degrade the cell walls of postharvest pathogens. P. membranefaciens had higher β-1,3-glucanase and exochitinase activity but less endo-chitinase activity than C. albidus in the presence of fungal cell walls .
Several studies have demonstrated the ability of yeasts to produce and secrete glucanases and chitinases in culture medium amended with fungal cells walls [16, 25, 28]. The ability of yeasts to produce exo-ß-1,3-glucanase and chitinase has been suggested to play a role in the firm attachment of yeast cells to fungal hyphae and the partial degradation of fungal mycelia . The current study clearly indicated that glucanase (GLU) and chitinase (CHI) genes are up regulated when yeast cells were in contact with P. digitatum mycelia (Mf-Pdig). When yeast cells were in contact with grapefruit peel tissue (Mf-fruit), however, only CHI was up-regulated. The fact that CHI expression is induced in both exposures to fruit peel tissues and to fungal mycelia (Figure 8B) may indicate that chitinases are involved in several biological interactions and not specifically in pathogen cell wall degradation.
Many pathogens utilize secreted siderophores and/or high-affinity uptake systems to acquire iron, and many are able to utilize ferritin, transferrin, lactoferrin, heme and heme-containing proteins for cellular processes related to virulence . Iron competition was reported as the main mode of action of M. pulcherrima inhibition of Botrytis cinerea, Alternaria alternata and P. expansum. Iron is essential for fungal growth and pathogenesis. Iron metabolism has been well characterized in the model yeast Saccharomyces cerevisiae. There are two different pathways that allow iron uptake from the environment: the reductive iron-uptake pathway mediated by Fet3p, a multicopper oxidase, and Ftr1p, an iron permease localized on the plasma membrane; and the siderophore-mediated uptake pathway. In our transcriptional analysis, we identified homologs of transmembrane ferric reductases (FET3), iron permease (FTR1), and an iron transporter (SIT1) all of which were up-regulated in the Mf-fruit interaction (Additional file 6). This information provides a clue to the potential important role of the reductive iron uptake system in the interaction of M. fructicola with fruit tissues. Iron competition was reported as the main mode of action in the ability of M. pulcherrima to inhibit Botrytis cinerea, Alternaria alternata and P.expansum infection of apples stored at 1°C for 8 months under controlled atmosphere (2% O2 and 3% CO2) . M. fructicola, similar to M. pulcherrima, produces a red pigment that is believed to be associated with the binding of iron (unpublished data). Iron permease (CaFTR1), in Candida albicans, was shown to mediate iron acquisition from transferrin and was required for systemic infection . In the plant pathogenic fungus, Fusarium graminearum, iron permeases (FgFtr1 and FgFtr2) function in the reductive iron uptake pathway but do not play a major role in pathogenicity .
Metal ion homeostasis is interdependent, linking iron availability to that of other metals, like zinc. Zinc itself was reported to be a determinant in cell-cell signaling in C. albicans biofilm formation and regulated by the transcription factor, Zap1 . Zap1 has been identified as a regulator of yeast-hypha balance in biofilms through intercellular signaling. Zap1 promotes accumulation of farnesol, an inhibitor of hypha formation in yeasts. Zap1 promotes accumulation of a postulated diffusible yeast cell inhibitor via the control of the zinc transporter, ZTR2. Zap1 activates the expression of the zinc transporter homologs, ZRT1 and ZRT2 (zinc regulated transporter). Overexpression of ZRT2, improves the growth of the zap1_/zap1_ mutant on a low-zinc medium .
In our study two transporters, ZTR3 and ZTR2, were up-regulated in the Mf-fruit interaction and may be related to low zinc levels on fruit peels. Interestingly, ZTR2 was down-regulated in cells interacting with fungal mycelia. Zap1 is required for full production of farnesol. Farnesol is believed to be generated from the ergosterol biosynthetic intermediate, farnesyl pyrophosphate [32, 33]. C. albicans biofilm growth is associated with the overall upregulation of ergosterol biosynthesis as well as increased resistance to antifungal compounds that target ergosterol . Ergosterol biosynthetic genes (ERG) are oppositely regulated by Zap1 in C. albicans and S. cerevisiae. ERG is positively regulated by ScZAP1 in S. cerevisiae and negatively regulated by CaZAP1 in C. albicans. The apparently opposite roles of Zap1 in ERG gene regulation in the two organisms may arise from the difference in growth conditions. In our study we observed up-regulation of a group of genes involved in ergosterol metabolism and signaling; ERG1, ERG 11, ERG5, PSD and DAP1 (Additional file 7). Because of the relationship between drug resistance and membrane ergosterol composition, many genes of the ergosterol biosynthetic pathway have been analyzed in C. albicans. Pasrija et al.  showed that suppressing the expression of an squalene epoxidase ERG1, leads to increased sensitivity to a number of drugs, including fluconazole and cycloheximide.
The results presented in this study provide the first use of RNA-seq, and de novo assembly, to conduct a comprehensive analysis of the transcriptome of M. fructicola cells that are in contact with fruit peel tissue or P. digitatum mycelia. De-novo assembly of the Illumina sequences generated a large number of new hypothetical proteins, which were categorized into 43 functional groups, in M. fructicola. The identification of M. fructicola genes associated with biocontrol interactions will form the basis of understanding the biology of biocontrol systems and mechanisms associated with effective biocontrol agents.
Yeast and pathogen cultures
Metschnikowia fructicola, Strain 277,  was grown in NYDP (nutrient broth (8 g l-1), yeast extract (5 g l-1), D-glucose (10 gl-1) and chloramphenicol (250 mg l-1). One ml of the yeast cell suspension was aseptically transferred from 24 h old starter culture to 250 ml. Ehrlenmeyer flasks and place on an orbital shaker at 160 rpm for 24 h at 26°C. Yeast cells were pelleted by centrifugation at 6000 rpm, washed twice with sterile distilled water, re-suspended in sterile water to its initial volume and the cell suspension concentration was adjusted to 1 × 108 cells ml-1.
Penicillium digitatum was isolated from decayed citrus fruit and maintained on potato dextrose agar (PDA) (Difco, Sparks, MD, USA). Fungal spore suspension was obtained from 2-week-old cultures grown on PDA at 26°C. The number of spores was calculated using a Cellometer (Neubauer, Germany), and the spore concentration was adjusted to 1 × 104 spores ml-1.
Red grapefruit ('Star Ruby') were obtained from a local orchard and used within 24 h after harvest. Before use, fruits were washed with tap water and detergent, rinsed, dried and surface disinfected by wiping with technical ethanol (70%).
Metschnikowia fructicola -Penicillium digitatum (Mf-Pdig) interaction
To obtain an even mycelial mat of P. digitatum on PDA, 1 ml of fungal spore suspension (1 × 104 spores ml-1) was thoroughly mixed with 9 ml molten 1% PDA (cooled down to about 40°C) in a petri dish and allowed to solidify at room temperature. After 5 days incubation at 25°C, 1 ml of yeast cells in water (1 × 108 cells ml-1) was placed on the mycelial mat that formed on the surface of the PDA. Yeast cells were carefully washed several times from the surface of mycelial mats with 1 ml of sterilized distilled water using a pipette after 5, 12, 24 and 48 h of incubation. Washings from 10 plates were collected at each time point, filtered through 4 layers of sterile cheese cloth to remove any mycelia, and transferred into sterile 50 ml tubes. Yeast cells were then pelleted by centrifugation at 6000 rpm, washed twice with sterile distilled water and pelleted again by centrifugation. The pellet was frozen and stored at −80°C for subsequent RNA extraction. Yeast cells gown in NYDB for 5, 12, 24 and 48 h were used as controls following pelleting and washing with water.
Metschnikowia fructicola-fruit (Mf-fruit) interaction
M. fructicola cells (1 × 108 cell ml-1) in water were prepared as indicated above. Surface washed and disinfected grapefruits were peeled using hand peeler to make 3 mm thick strips. Alternatively, discs were cut using a 9 mm diameter cork borer under aseptic conditions. Disks were immersed for 1 min in M. fructicola cell suspension (1 × 108 cells ml-1) and then transferred into Petri dishes lined with moist, sterile filter paper. After 5, 12, 24 and 48 h of incubation at 25°C under moist conditions, 10 disks from each treatment were transferred to 50 ml sterile conical tubes containing 10 ml sterile water and shaken on an orbital shaker at 180 rpm for 10 min. Fruit disks were then discarded and dislodged yeast cells were collected by centrifugation at 6000 rpm for 10 min, washed twice with sterile distilled water, centrifuged between rinses, and re-suspended in sterile water. Concentration of the cell suspension was adjusted to 1 × 107 cells ml-1 and then pelleted by centrifugation. Yeast pellets were frozen and stored at −80°C for subsequent RNA preparation.
RNA extraction and cDNA synthesis
Total RNA from frozen each yeast pellet was extracted using a MasterPure™ Yeast RNA Purification Kit (Epicentre® Biotechnologies, Madison, WI, USA) accor-ding to the manufacturer`s instructions. Quality and concentration of total RNA was analyzed by gel electrophoresis and a spectrophotometer ND-1000, respectively (NanoDrop, Wilmington, DE, USA). First-strand cDNA was synthesized with a Verso cDNA Kit (Thermo Fisher Scientific, Epson, UK), from 1 μg of total RNA that had been pretreated with 1.5 units of RNase I (Epicentre® Biotechnologies, Madison, WI, USA).
RNA samples isolated from control cells of M. fructicola grown in media (Mf), M. fructicola – fruit interaction (Mf-fruit), and M. fructicola – P. digitatum (Mf-Pdig) after 24 h of incubation were further prepared for RNA sequencing using a MasterPure™ Yeast RNA Purification Kit as described above. cDNA preparation, cDNA library construction, and sequencing were performed at The Genome High-Throughput Sequencing Laboratory, Sackler Faculty of Medicine and George Wise Faculty of Life Sciences, Tel Aviv University. According to the Illumina manufacturer's instructions, poly(A) + RNA was purified from 20 μg of total RNA using oligo(dT) magnetic beads and fragmented into short sequences in the presence of divalent cations at 94°C for 5 min. The cleaved poly(A) + RNA was transcribed and second-strand cDNA synthesis was performed. After end-repair and the ligation of adaptors, cDNA products were amplified by PCR and purified using the QIAquick PCR Purification Kit to create a specific cDNA library. Cluster formation, primer hybridization, and singled-end 76 cycle sequencing were performed on an Illumina Genome Analyzer IIx using proprietary reagents according to the manufacturer’s instructions.
Bioinformatic analysis of the expression RNA-seq libraries
De-novo assembly of the M. fructicola transcriptome without a reference genome using the four libraries obtained from the different treatments of the yeast was accomplished using Trinity de-novo transcriptome assembly software . After assembling the M. fructicola transcriptome, every RNA-seq library was separately aligned to the generated transcriptome assembly using Bowtie . The counting of alignments was done using RSEM . The differential expression Statistical analysis was done using the statistical method described in the R package DESeq . Annotation of the de-novo transcript contigs was done using blast2go software . Using blast2GO software, blastx  was performed against the NCBI non-redundant protein database  with an E-value cut off of 1e-5. Mapped records were associated with matching Gene Ontology (GO) terms in the GO database  used the WEGO tool (web gene ontology annotation Plotting) for plotting GO annotation results (http://wego.genomics.org.cn) . Finally, InterProScan [45, 46] was used to identify additional GO annotations for our sequences . The transcriptome datasets are available at the NCBI Sequence Read Archive (SRA), under accession number SRA054245.
Mapping of Metschnikowia fructicolaunigenes to the draft genome
The Unigene transcript were blasted against the Metschnikowia fructicola draft genome available under gene bank accession number ANFW01000000 with cutoff e-value of 1e-50 and identity >85%.
Phylogenetic tree construction
Sequences of the ribosomal RNA large subunit were downloaded from the SILVA database . Sequences were aligned using MUSCLE program . Then, phylogenetic trees were reconstructed for each dataset based on maximum likelihood (ML) framework using the phyml software with 100 bootstrap replicates .
Real-time quantitative reverse transcriptase PCR
Real-time quantitative reverse transcriptase PCR (RT-qPCR) was performed to validate transcription levels determined by RNA-Seq and to evaluate transcription levels at additional time points (5, 12, 24 and 48 h). A ten μl reaction contained specific primers (200 nM final concentration) Additional file 8, absolute QPCR SYBR Green ROX Mix (ABgene, Epson, UK) and cDNA templates (1:10 dilution), and was performed on a GENE 3000 Real-Time PCR system (Corbett Life Science, Sydney, Australia). The cycle threshold values (CT) were determined and the relative fold differences were calculated by the 2_ΔΔCT method  using actin (AJ745127) as the endogenous reference gene . Samples were run in triplicate, and each experiment was repeated twice. Two experiments of independently cultivated yeast were performed to confirm the reproducibility of the results.
Wisniewski M, Biles C, Droby S, McLaughlin R, Wilson C, Chalutz E: Mode of action of the postharvest biocontrol yeast, pichia Pichia-guilliermondii.1. Characterization of attachment to botrytis Botrytis-cinerea. Physiol Mol Plant Pathol. 1991, 39: 245-258. 10.1016/0885-5765(91)90033-E.
Wilson C, Wisniewski M, Droby S, Chalutz E: A selection strategy for microbial antagonists to control postharvest diseases of fruits and vegetables. Sci Hortic. 1993, 53: 183-189. 10.1016/0304-4238(93)90066-Y.
Stockwell VO, Stack JP: Using Pseudomonas spp. for integrated biological control. Phytopathology. 2007, 97: 244-249. 10.1094/PHYTO-97-2-0244.
Kurtzman CP, Droby S: Metschnikowia fructicola, a new ascosporic yeast with potential for biocontrol of postharvest fruit rots. Syst Appl Microbiol. 2001, 24: 395-399. 10.1078/0723-2020-00045.
Karabulut OA, Tezcan H, Daus A, Cohen L, Wiess B, Droby S: Control of preharvest and postharvest fruit rot in strawberry by Metschnikowia fructicola. Biocontrol Sci Technol. 2004, 14: 513-521. 10.1080/09583150410001682287.
Liu J, Wisniewski M, Droby S, Tian S, Hershkovitz VTT: Effect of heat shock treatment on stress tolerance and biocontrol efficacy of Metschnikowia fructicola. FEMS Microbiol Ecol. 2011, 76: 145-155. 10.1111/j.1574-6941.2010.01037.x.
Macarisin D, Droby S, Bauchan G, Wisniewski M: Superoxide anion and hydrogen peroxide in the yeast antagonist-fruit interaction: a new role for reactive oxygen species in postharvest biocontrol?. Postharvest Biol Technol. 2010, 58: 194-202. 10.1016/j.postharvbio.2010.07.008.
Castoria R, Wright SAI: Host responses to biological control agents postharvest pathology. Post-harvest pathology. Edited by: Prusky D, Gullino ML. 2010, Netherlands: Springer, 171-181. vol. 2
Jamalizadeh M, Etebarian HR, Aminian H, Alizadeh A: A review of mechanisms of action of biological control organisms against post-harvest fruit spoilage. EPPO Bulletin. 2011, 41 (1): 65-71. 10.1111/j.1365-2338.2011.02438.x.
Nunes C: Biological control of postharvest diseases of fruit. Eur J Plant Pathol. 2012, 133: 181-196. 10.1007/s10658-011-9919-7.
Hershkovitz V, Ben-Dayan C, Raphael G, Pasmanik-Chor M, Liu J, Belausov E, Aly R, Wisniewski M, Droby S: Global changes in gene expression of grapefruit peel tissue in response to the yeast biocontrol agent Metschnikowia fructicola. Mol Plant Pathol. 2012, 13: 338-349. 10.1111/j.1364-3703.2011.00750.x.
Castoria R, Caputo L, De Curtis F, De Cicco V: Resistance of postharvest biocontrol yeasts to oxidative stress: a possible new mechanism of action. Phytopathology. 2003, 93: 564-572. 10.1094/PHYTO.2003.93.5.564.
Torres R, Teixido N, Usall J, Abadias M, Mir N, Larrigaudiere C, Vinas I: Anti-oxidant activity of oranges after infection with the pathogen Penicillium digitatum or treatment with the biocontrol agent Pantoea agglomerans CPA-2. Biol Control. 2011, 57: 103-109. 10.1016/j.biocontrol.2011.01.006.
Larralde-Corona CP, Ramirez-Gonzalez MD, Perez-Sanchez G, Oliva-Hernandez AA, Narvaez-Zapata JA: Identification of differentially expressed genes in the citrus epiphytic-yeast Pichia guilliermondii during interaction with Penicillium digitatum. Biol Control. 2011, 57: 208-214. 10.1016/j.biocontrol.2011.02.012.
Droby S, Wisniewski M, Macarisin D, Wilson C: Twenty years of postharvest biocontrol research: is it time for a new paradigm?. Postharvest Biol Technol. 2009, 52: 137-145. 10.1016/j.postharvbio.2008.11.009.
Haissam JM: Pichia anomala in biocontrol for apples: 20 years of fundamental research and practical applications. Antonie Van Leeuwenhoek. 2011, 99: 93-105. 10.1007/s10482-010-9541-2.
Liu J, Wisniewski M, Droby S, Vero S, Tian SP, Hershkovitz V: Glycine betaine improves oxidative stress tolerance and biocontrol efficacy of the antagonistic yeast Cystofilobasidium infirmominiatum. Int J Food Microbiol. 2011, 146: 76-83. 10.1016/j.ijfoodmicro.2011.02.007.
Gasch AP, Werner-Washburne M: The genomics of yeast responses to environmental stress and starvation. Funct Integr Genomics. 2002, 2: 181-192. 10.1007/s10142-002-0058-2.
Gasch AP, Spellman PT, Kao CM, Carmel-Harel O, Eisen MB, Storz G, Botstein D, Brown PO: Genomic expression programs in the response of yeast cells to environmental changes. Mol Biol Cell. 2000, 11: 4241-4257.
Bi X, Guo N, Jin J, Liu J, Feng H, Shi J, Xiang H, Wu X, Dong J, Hu H: The global gene expression profile of the model fungus Saccharomyces cerevisiae induced by thymol. J Appl Microbiol. 2010, 108: 712-722. 10.1111/j.1365-2672.2009.04470.x.
Hu F, Liu J, Du G, Hua Z, Zhou J, Chen J: Key cytomembrane ABC transporters of Saccharomyces cerevisiae fail to improve the tolerance to D-limonene. Biotechnol Lett. 2012, 10.1007/s10529-10012-10931-10526.
Paumi CM, Chuk M, Snider J, Stagljar I, Michaelis S: ABC transporters in Saccharomyces cerevisiae and their interactors: new technology advances the biology of the ABCC (MRP) subfamily. Microbiol Mol Biol Rev. 2009, 73: 577-593. 10.1128/MMBR.00020-09.
Friel D, Pessoa NMG, Vandenbol M, Jijakli MH: Separate and combined disruptions of two exo-beta-1,3-glucanase genes decrease the efficiency of Pichia anomala (strain K) biocontrol against Botrytis cinerea on apple. Mol Plant Microbe Interact. 2007, 20: 371-379. 10.1094/MPMI-20-4-0371.
Saravanakumar D, Clavorella A, Spadaro D, Garibaldi A, Gullino ML: Metschnikowia pulcherrima strain MACH1 outcompetes Botrytis cinerea, Alternaria alternata and Penicillium expansum in apples through iron depletion. Postharvest Biol Technol. 2008, 49: 121-128. 10.1016/j.postharvbio.2007.11.006.
Saravanakumar D, Spadaro D, Garibaldi A, Gullino M: Detection of enzymatic activity and partial sequence of a chitinase gene in Metschnikowia pulcherrima strain MACH1 used as post-harvest biocontrol agent. Eur J Plant Pathol. 2009, 123: 183-193. 10.1007/s10658-008-9355-5.
Jijakli H, Lepoivre L: Characterization of an exobeta-1,3-glucanase produced by Pichia anomala strain K, antagonist of Botrytis cinerea on apples. Phytopathology. 1998, 88: 335-343. 10.1094/PHYTO.1922.214.171.1245.
Chan ZL, Tian SP: Interaction of antagonistic yeasts against postharvest pathogens of apple fruit and possible mode of action. Postharvest Biol Technol. 2005, 36: 215-223. 10.1016/j.postharvbio.2005.01.001.
Zhang DP, Spadaro D, Garibaldi A, Gullino ML: Potential biocontrol activity of a strain of Pichia guilliermondii against grey mold of apples and its possible modes of action. Biol Control. 2011, 57: 193-201. 10.1016/j.biocontrol.2011.02.011.
Heymann P, Gerads M, Schaller M, Dromer F, Winkelmann G, Ernst JF: The siderophore iron transporter of Candida albicans (Sit1p/Arn1p) mediates uptake of ferrichrome-type siderophores and is required for epithelial invasion. Infect Immun. 2002, 70: 5246-5255. 10.1128/IAI.70.9.5246-5255.2002.
Park YS, Choi ID, Kang CM, Ham MS, Kim JH, Kim TH, Yun SH, Lee YW, Chang HI, Sung HC: Functional identification of high-affinity iron permeases from Fusarium graminearum. Fungal Genet Biol. 2006, 43: 273-282. 10.1016/j.fgb.2005.12.005.
Ganguly S, Bishop AC, Xu WJ, Ghosh S, Nickerson KW, Lanni F, Patton-Vogt J, Mitchell AP: Zap1 control of cell-cell signaling in Candida albicans biofilms. Eukaryot Cell. 2011, 10: 1448-1454. 10.1128/EC.05196-11.
Nobile CJ, Nett JE, Hernday AD, Homann OR, Deneault JS, Nantel A, Andes DR, Johnson AD, Mitchell AP: Biofilm matrix regulation by Candida albicans Zap1. PLoS Biol. 2009, 7: e1000133-10.1371/journal.pbio.1000133.
Shea JM, Del Poeta M: Lipid signaling in pathogenic fungi. Curr Opin Microbiol. 2006, 9: 352-358. 10.1016/j.mib.2006.06.003.
Garcia-Sanchez S, Aubert S, Iraqui I, Janbon G, Ghigo JM, d'Enfert C: Candida albicans biofilms: a developmental state associated with specific and stable gene expression patterns. Eukaryot Cell. 2004, 3: 536-545. 10.1128/EC.3.2.536-545.2004.
Pasrija R, Krishnamurthy S, Prasad T, Ernst JF, Prasad R: Squalene epoxidase encoded by ERG1 affects morphogenesis and drug susceptibilities of Candida albicans. J Antimicrob Chemother. 2005, 55: 905-913. 10.1093/jac/dki112.
Grabherr M, Haas B, Yassour M, Levin J, Thompson D, Amit I, Adiconis X, Fan L, Raychowdhury R, Zeng Q: Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011, 29: 644-652. 10.1038/nbt.1883.
Langmead B, Trapnell C, Pop M, Salzberg S: Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2008, 10: R25-
Li B, Dewey CN: RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinforma. 2011, 12: 323-10.1186/1471-2105-12-323.
Anders S, Huber W: Differential expression analysis for sequence count data. Genome Biol. 2010, 11: R106-10.1186/gb-2010-11-10-r106.
Conesa A, Götz S, García-Gómez JM, Terol J, Talón 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.
Altschul S, Gish W, Miller W, Myers E, Lipman D: Basic local alignment search tool. J Mol Biol. 1990, 215: 403-410.
Pruitt KD, Tatusova T, Maglott DR: NCBI reference sequences (RefSeq): a curated non-redundant sequence database of genomes, transcripts and proteins. Nucleic Acids Res. 2007, 35: 61-65.
Harris M, Clark J, Ireland A, Lomax J, Ashburner M, Foulger R, Eilbeck K, Lewis S, Marshall B, Mungall C: The Gene Ontology GO database and informatics resource. Nucleic Acids Res. 2004, 32: 258-261. 10.1093/nar/gkh036.
Ye J, Fang L, Zheng H, Zhang Y, Chen J, Zhang Z, Wang J, Li S, Li R, Bolung L: WEGO: a web tool for plotting GO annotations. Nucleic Acids Res. 2006, 34: W293-W297. 10.1093/nar/gkl031.
Hunter S, Jones P, Mitchell A, Apweiler R, Attwood TK, Bateman A, Bernard T, Binns D, Bork P, Burge S: InterPro in 2011: new developments in the family and domain prediction database. Nucleic Acids Res. 2011, 40: D306-D312.
Quevillon E, Silventoinen V, Pillai S, Harte N, Mulder N, Apweiler RRL: InterProScan: protein domains identifier. Nucleic Acids Res. 2005, 1: W116-W120.
Burge S, Kelly E, Lonsdale D, Mutowo-Muellenet P, McAnulla C, Mitchell A, Sangrador-Vegas A, Yong S, Mulder N, Hunter S: Manual GO annotation of predictive protein signatures: the InterPro approach to GO curation. Database (Oxford). 2012, bar068. doi: 10.1093/database/bar068.
Pruesse E, Quast C, Knittel K, Fuchs B, Ludwig W, Peplies J, Glöckner FO: SILVA: a comprehensive online resource for quality checked and aligned ribosomal RNA sequence data compatible with ARB. Nucleic Acids Res. 2007, 35 (21): 7188-7196. 10.1093/nar/gkm864.
Edgar RC: MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32: 1792-1797. 10.1093/nar/gkh340.
Guindon S, Gascuel OA: Simple, fast and accurate algorithm to estimatelarge phylogenies by maximum likelihood. Syst Biol. 2003, 52: 696-704. 10.1080/10635150390235520.
Livak K, Schmittgen T: Analysis of relative gene expression data using real-time quantitative PCR and the 2-ΔΔCT method. Methods. 2001, 25: 402-408. 10.1006/meth.2001.1262.
This study was supported in part by a grant (IS-4268-09) from the U.S. - Israel Binational Agricultural Research and Development (BARD) Foundation.
Authors declare that they have no competing interests.
VH did the experiments and drafted the manuscript. NS did all the bioinformatic analysis and drafted the manuscript. LTS, JL, GR, CK, RA, ML helped with the analysis and experiments or contributed additional analysis. MW and SD conceived the project, supervised the project and drafted the manuscript. All authors have read and approved the final manuscript.
Vera Hershkovitz, Noa Sela contributed equally to this work.
Electronic supplementary material
Additional file 1: A fasta file containing all Unigene sequences of Metchnikowia fructicola transcriptome. (TXT 5 MB)
Additional file 2: Annotation of all Metschnikowia fructicola unigenes: the gene description, go terms associated and proscan analysis of every transcript. (XLSX 887 KB)
Additional file 3: Mapping of Metschnikowia fructicola de-novo assembled transcript to Metschnikowia fructicola genome draft genbank accession number ANFW01000000. (XLSX 785 KB)
Additional file 4: The differentially expressed genes in Metschnikowia fructicola interaction with Penicilium digitatum (Mf-Pdig) and interaction with fruit (Mf-fruit) and with fruit (p < 0.05). (DOCX 62 KB)
Additional file 5: Summary of differential expressed genes in Metschnikowia fructicola interaction with Penicillium digitatum and interaction with fruit (p < 0.05) involved in the response to stresses (chemical GO:042221), oxidative (GO:006979), osmotic (GO:006970), heat (GO:009408), starvation (GO:042594), DNA damage stimulus (GO:006974). (DOCX 24 KB)
Additional file 6: Summary of differential expressed genes in Metschnikowia fructicola interaction with Penicillium digitatum and interaction with fruit (p < 0.05) involved in transport (transmembrane GO:055085), ion (GO:006811), carbohydrate (GO:008643), lipid (GO:006869), amino acid transport (GO:006865). (DOCX 28 KB)
Additional file 7: Summary of differential expressed genes in Metschnikowia fructicola interaction with Penicillium digitatum and interaction with fruit (p < 0.05) involved in metabolic (lipid GO:006629), (vitamin GO:006766), (cofactor GO:051186), (cellular amino acid GO:006520) processes. (DOCX 24 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
About this article
Cite this article
Hershkovitz, V., Sela, N., Taha-Salaime, L. et al. De-novo assembly and characterization of the transcriptome of Metschnikowia fructicola reveals differences in gene expression following interaction with Penicillium digitatumand grapefruit peel. BMC Genomics 14, 168 (2013). https://doi.org/10.1186/1471-2164-14-168
- Antagonist-fruit interaction
- Antagonist-pathogen interaction
- Biological agent