Label-Free Quantitative Proteomic Analysis of Puccinia psidii Uredospores Reveals Differences of Fungal Populations Infecting Eucalyptus and Guava

Abstract

Puccinia psidii sensu lato (s.l.) is the causal agent of eucalyptus and guava rust, but it also attacks a wide range of plant species from the myrtle family, resulting in a significant genetic and physiological variability among populations accessed from different hosts. The uredospores are crucial to P. psidii dissemination in the field. Although they are important for the fungal pathogenesis, their molecular characterization has been poorly studied. In this work, we report the first in-depth proteomic analysis of P. psidii s.l. uredospores from two contrasting populations: guava fruits (PpGuava) and eucalyptus leaves (PpEucalyptus). NanoUPLC-MS^E was used to generate peptide spectra that were matched to the UniProt Puccinia genera sequences (UniProt database) resulting in the first proteomic analysis of the phytopathogenic fungus P. psidii. Three hundred and fourty proteins were detected and quantified using Label free proteomics. A significant number of unique proteins were found for each sample, others were significantly more or less abundant, according to the fungal populations. In PpGuava population, many proteins correlated with fungal virulence, such as malate dehydrogenase, proteossomes subunits, enolases and others were increased. On the other hand, PpEucalyptus proteins involved in biogenesis, protein folding and translocation were increased, supporting the physiological variability of the fungal populations according to their protein reservoirs and specific host interaction strategies.

About the Authors:

Maria Carolina Quecine

Contributed equally to this work with: Maria Carolina Quecine, Thiago Falda Leite

Affiliation: Departament of Genetics, Escola Superior de Agricultura “Luiz de Queiroz”, Universidade de São Paulo, Piracicaba -SP, Brazil

Thiago Falda Leite

Contributed equally to this work with: Maria Carolina Quecine, Thiago Falda Leite

Affiliation: Departament of Genetics, Escola Superior de Agricultura “Luiz de Queiroz”, Universidade de São Paulo, Piracicaba -SP, Brazil

Andressa Peres Bini

Affiliation: Departament of Genetics, Escola Superior de Agricultura “Luiz de Queiroz”, Universidade de São Paulo, Piracicaba -SP, Brazil

Thais Regiani

Affiliation: Departament of Genetics, Escola Superior de Agricultura “Luiz de Queiroz”, Universidade de São Paulo, Piracicaba -SP, Brazil

Lívia Maria Franceschini

Affiliation: Departament of Genetics, Escola Superior de Agricultura “Luiz de Queiroz”, Universidade de São Paulo, Piracicaba -SP, Brazil

Ilara Gabriela Frasson Budzinski

Affiliation: Departament of Genetics, Escola Superior de Agricultura “Luiz de Queiroz”, Universidade de São Paulo, Piracicaba -SP, Brazil

Felipe Garbelini Marques

Affiliation: Departament of Genetics, Escola Superior de Agricultura “Luiz de Queiroz”, Universidade de São Paulo, Piracicaba -SP, Brazil

Mônica Teresa Veneziano Labate

Affiliation: Departament of Genetics, Escola Superior de Agricultura “Luiz de Queiroz”, Universidade de São Paulo, Piracicaba -SP, Brazil

Simone Guidetti-Gonzalez

Affiliation: Departament of Genetics, Escola Superior de Agricultura “Luiz de Queiroz”, Universidade de São Paulo, Piracicaba -SP, Brazil

David Henry Moon

Affiliation: Departament of Genetics, Escola Superior de Agricultura “Luiz de Queiroz”, Universidade de São Paulo, Piracicaba -SP, Brazil

Carlos Alberto Labate

* E-mail: [email protected]

Affiliation: Departament of Genetics, Escola Superior de Agricultura “Luiz de Queiroz”, Universidade de São Paulo, Piracicaba -SP, Brazil

Introduction

Puccinia psidii sensu lato (s.l.) is a biotrophic rust fungi native to South America, where it was first described in guava [1], causing severe infections resulting in losses up to100% of the fruits [2]. Eucalyptus rust is apparently a specialization of P. psidii genotype evolved from rust occurring in Myrtaceae in South America, enabling it to ‘host jump’ and infect Eucalyptus, an introduced Myrtaceae species to Brazil at the end of the 19^th century [3, 4]. Likely, some isolates from guava did not infect eucalyptus and vice versa [4] indicating an occurrence in the evolution of a host-specific genotype. Since then, rust has become one of the most important eucalyptus diseases, principally attacking nurseries and young plants, decreasing productivity, and in some cases, causing the death of the highly susceptible plants [5].Guava and other Myrtaceae from South America most likely served as the sources of inoculum during the adaptation of P. psidii to eucalyptus [6].

Recently, the genetic variability in P. psidii s.l. populations infecting different host plants was evaluated and proved that, in general, populations collected from guava had a greater level of diversity than populations from Eucalyptus spp. [6]. The classification of this fungus has been deeply debated. Morphological and molecular studies indicate that Puccinia is polyphyletic group [7, 8, 9], named P. psidii s.l. complex that is able to infect a great range of hosts from the Mytaceaes family. Thus, the disease caused by the fungus currently treated as P. psidii s.l. is best referred to as Myrtle rust [10].

P. psidii s.l. is considered an autoecious species with an incomplete life cycle (the spermogonia stage may be absent).The aecial state is extremely rare [5].The uredia are the most common and principal P. psidii s.l. dissemination structure and are constantly produced during natural and artificial inoculations [4]. Moreover, the uredosporic pustules are the only visual rust symptom that can be used to diagnose infection under field conditions [5, 11]. The uredospore protein profile may play a key role in the establishment of the disease, because they are the first fungal structures that have contact with the host before infection starts, and may carry some key proteins with respect to the host, because the rust uredopores germination is strongly dependent on physical and chemical signals from the host [12].

The physiological variability of P. psidii s.l.uredospores according to the host has been described previously [5, 13, 14]. Coelho et al. [13] identified three groups of rust biotypes, each pathogenic on different host combinations: Psidium guajava only, Eucalyptus grandis and P. guajava or E. grandis and Sizygium jambos. Despite the importance of the uredospore as the primary source of fungal inoculum, little is known about the molecular features of P. psidii s.l. uredospores or how they relate to its many host organisms, for instance guava and eucalyptus.

Global proteomic analysis can provide the framework for an in-depth understanding of the cellular processes and organization of a particular organism. The complexity of the proteome exceeds that of the genome when protein isoforms, alternative splicing variants, and post-translational modifications are taken into account [15]. However, proteomic analysis of filamentous fungi is still lacking in this field [16]. Among the filamentous fungi, the proteomes of the Ascomycota have been studied more intensively than those of Basidiomycota, most likely due to the number of pathogenic fungi that are found in Ascomycota [17]. Gradually, the proteomes of filamentous fungi are being revealed, including the mycelia proteome [18, 19, 20], secretomes [20, 21] and sub-proteomes [22] of a wide range of species, such as Aspergillus spp. [23, 24, 25, 26], Botrytis cinerea [17], Neurospora crassa [22], Penicillium expansum [27], Pleurotus sapidus [28], Sclerotinia sclerotiorum [20], Terebralia palustris [29], P. triticina [30], Blumeria graminis [31], Verticillium albo-atrum [32], and others.

In the present work, we report the first label free proteomic quantification of P. psidii s.l. uredospores from two contrasting fungal populations with known physiological variability that reflects specific molecular interactions directed by the host. The results obtained allowed the identification of pathogenicity or virulence factors of this important pathogen and their correlation with the fungal physiological variation as a possible specific host interaction strategy during the infection process.

Materials and Methods

Pathogen material

P. psidii uredospores were collected from E. grandis leaves (PpEucalyptus) in the experimental area belonging to FuturaGene, Itapetininga-SP, Brazil in accordance with Dr. Esteban R. Gonzalez (FuturaGene Corporation Manager). Uredospores from P. guajava (PpGuava) where isolated from fruits of trees located in the Campus of ESALQ/USP, Piracicaba-SP, Brazil in accordance with the Chief of the Intitution Dr. José Caixeta Filho. Both fields do not involve endangered or protected species. Infected leaves and fruits from E. grandis and P. guajava tree`s were individually collected. We sampled three plants per species, a bulk of leaves or fruits from each plant considered as an independent biological replicate, totalizing three biological replicates per population sample. The uredospores samples from each replicate were divided, and half was immediately frozen at -80°C for proteomic analysis. The remaining sample was dehydrated in silica gel for 48 h at 4°C and frozen at -80°C for light microscopy analysis and for infection assay in E. grandis.

Light microscopy analysis

The morphology of the PpEucalyptus and PpGuava uredospores and the germination viability assays were performed in water-agar medium (8 g/L). The uredospores from each population were suspended in mineral oil (Invitrogen) and 400 μL of each suspension (10⁴ uredospores/mL), spread onto the medium and incubated at 22°C in the dark for 24 h. The uredospores were analyzed using a light microscope (Axiophot II; Zeiss). For morphological analysis, we measured the length and width of 50 uredospores from each population, and the viability of the uredospores was confirmed by the presence of germ tubes after 24 h.

Plant-infection assay

Eighteen highly susceptible E. grandis plantlets (monoprogenie D901) were grown under greenhouse conditions for 120 days and were then transferred to a controlled growth chamber (E15, Conviron),where they were maintained under a 12 h photoperiod (200 μmol m^-1 s^-1) at 20°C for acclimatization. After seven days, the PpEucalyptus and PpGuava uredospore suspensions containing 10⁴ uredospores mL^-1 and 0.05% Tween 20 were sprayed onto six plantlets, two plants per biological sample. To maintain a high level of humidity, the plants were enclosed in transparent plastic bags for the first 48 h, with the first 24 h occurring in complete darkness at 20°C. Thereafter, the plants were returned to the previously described growth conditions. The control plantlets were treated in the same way, but were sprayed only with the 0.05% Tween solution. The evolution of the infection was monitored every two days for two weeks.

Protein extraction and trypsin digestion

Proteins were extracted from uredospores according to the protocol described by Damerval et al. [33] with some modifications. Uredospore samples (1 g) were ground into a powder in liquid nitrogen and homogenized in 30 mL of TCA solution (10% (w/v) trichloroacetic acid and 0.07% (v/v) 2-mercaptoethanol in acetone) per biological replicate. The samples were then incubated for 1 h at -20°C, followed by centrifugation for 20 min at 16,000 g at 2°C. The supernatant was discarded, and the pellet was suspended in 30 mL of 100% cold acetone containing 0.07% (v/v) 2-mercaptoethanol. The pellets were incubated at -20°C, for one h, followed by centrifugation for 20 min at 16,000 g at 2°C. The pellet was again suspended in 30 mL of 100% cold methanol and incubated for one hour at -80°C. After centrifugation at 16,000 g for 20 min at 2°C, the supernatant was discarded, and the resulting pellet was lyophilized.

The dried pellets were solubilized in 800 μL of TCT buffer [7 M urea, 2 M thiourea, 10 mM DTT and 0.4% (v/v) Triton X-100]. Complete protein solubilization was achieved by vigorous shaking using a vortex. Protein extracts were desalinized using Amicon Ultra-0.5 mL 3K-NMWL filter devices (Millipore Corporation).

The total protein concentration was determined using an Agilent Bioanalyzer 2100 (Agilent) and the 230 kDa protein kit. Fifty micrograms of protein was denatured with 25 μL of 0.2% RapiGest (RapiGest SF, Waters), reduced with 2.5 μL of 100 mM dithiothreitol and alkylated with 2.5 μL of 300 mM iodoacetamide. Trypsin digestion was performed with sequencing Grade Modified Trypsin (Promega) at a 1:100 (w/w) enzyme: protein ratio. After digestion, 10 μL of 5% (v/v) trifluoroacetic acid was added to the digested mixture to hydrolyze the RapiGest. The peptide mixture was then desalted using PepClean C 18-columns (Thermo-Fisher Scientific). The final volume of 50 μL was obtained by the addition of 20 mM ammonium formate (pH 10) solution containing 1 pM rabbit phosphorylase (internal standard to data normalization and label free protein quantification- P00489) to the lyophilized, desalted peptide sample.

MS^E analysis

Mass spectra of the peptide fragments were acquired by reverse-phase ultraperformance liquid chromatography (2D Technology nanoACQUITY-Waters). First dimension separation was achieved in an XBridge BEH130 C18 5μm 300 μm x 50 mm column. Elution was performed using 5 different binary gradients with 20 mM pH 10 ammonium formate in acetonitrile, at a flow rate of 2 μL/min. Eluted peptides from the first dimension column were trapped in a Symmetry C18 5 μm 180 μm x 20 mm column and diluted, online, with acetonitrile containing 0.1% formic acid. Second dimension separation was performed in a HSS T3 1.8 μm 75 μm x 100 mm column, using a binary gradient from 3 to 85% of acetonitrile with 0.1% formic acid, during 52 min, at a flow rate of 350 μL/min.

Mass spectrometry acquisition was achieved in a Synapt G2 HDMS mass spectrometer equipped with an ion mobility cell and a nanolockspray source in the positive ion and ‘V’ mode. The doubly-charged ion ([M+2H]²⁺) was used for initial single-point calibration and MS/MS fragment ions of GFP [Glu 1]-Fibrinopeptide B with a [M + 2H]²⁺ = 785,84206 m/z was used as lock mass, to obtain the final instrument calibration. Data-independent scanning (MSE) experiments were performed by switching between low (3 eV) and elevated collision energies (15–50 eV), applied to the trap ‘T-wave’ cell filled with argon. Scan time of 0.8 s were used for low and high energy scans from m/z 50 to 2000 [34].

Processing Parameters and Database Search

The raw data processing, protein identification and relative quantitative analyses were all performed using the ProteinLynx Global Server (PLGS- v 2.5.1). The processing parameters included: automatic tolerance for precursor and product ions, minimum of 3 fragment ions matched per peptide, minimum of 7 fragment ions matched per protein, minimum of 2 peptides matched per protein, 1 possible trypsin missed cleavage, carbamidomethylation of cysteine as fixed modification and oxidation of methionine as variable modification, and a maximum false positive discovery rate (FDR) of 4% that was determined based on the search of a reversed database, which was generated automatically using PLGS 2.5.1 by reversing the sequence of each entry.

To identify and quantify the proteins, the intensities of the spectra were calculated by the stoichiometric method, according to the internal standard, the sequence of rabbit phosphorylase (Uniprot entry: P00489), by MSE analysis [34] and normalized using the PLGS auto normalization function. The amount and sequence of the matched peptide and protein fmols were obtained, based on the ratio of its three most abundant peptides (“High Top 3” method) determined in each individual experiment [34], considering the data from 3 biological replicates for each sample. All protein hits were identified with a confidence of > 95%. Only those proteins identified in at least two out of 3 replicates performed for each biological sample, were regarded as having undergone a significant change.

Protein identifications were obtained with the embedded ion accounting algorithm of PLGSsoftware searching into the UniProtKB/Swiss (http://www.uniprot.org/ release Version 2015_01) from Puccinia spp proteins (28,122 entries) appended in the internal standard. The final expression data was screened according to the following criteria: statistical significant fold change as well as standard deviation (SD) and their correspondent p value: p≤0.05 and p≥ 0.95 were considered whose abundance was statistically increased in the PpEucalyptus and PpGuava populations, respectively. Data labeled as PpEucalyptus and PpGuava represent proteins that were only identified in that population.

After the PLGS 2.5.1 expression analysis the data was manually inspected, and some parameters obtained using the MassPivot v 101.jar were included: (i) the average amount (fmol) of protein, (ii) average score of proteins and (iii) average amount of matched peptides to each protein.

Functional Analysis

The protein sequences that were matched with ionized peptides were functionally categorized by the Gene Ontology Consortium (GO) [35] (http://www.geneontology.org). The terms were obtained using the Blast2GO software with the default parameters specified by the program [36]. These annotations were previously simplified using the GO Slim feature. The EC numbers (Enzyme Commission) were obtained and the prediction of the subcellular localization of the proteins performed using the Signal P (http://www.cbs.dtu.dk/services/SignalP/) and TMHMM (http://www.cbs.dtu.dk/services/TMHMM-2.0/) programs by the InterProScan tool (http://www.ebi.ac.uk/interpro/).

Results

Light microscopy analysis and plant infection

The morphological analysis showed that the size and shape of the uredospores were similar for both the PpEucalyptus and PpGuava populations (Fig 1A and 1B). The average length of the uredospores were 19.5 ± 2.2 μm and 19.9 ± 0.4 μm, respectively, for the PpEucalyptus and PpGuava populations, and their respective average widths were 13.5 ± 1.9 μm and 11.7 ± 0.7 μm, suggesting that there was no relationship between the uredospore size or shape and virulence or physiological variability of the populations. The populations also showed similar uredospore viability under in vitro conditions, with approximately 25% of the spores germinating on water-agar media (Fig 1C and 1D).

[Figure omitted. See PDF.]

Fig 1. Morphological and viability analysis of Puccinia psidii uredospores.

P. psidii uredospores from E. grandis (A) and P. guajava (B) exhibit similar morphology and germination viability, respectively (C and D).The arrows indicate the fungal germ tube in both uredospore populations, 24 hours after inoculation in water-agar medium. Light microscopy images of PpEucalyptus and PpGuava uredospores are shown at 100 X (A and B) and 200 X (C and D) magnification. Scale bar: 20 μm in A and B, 50 μm in C and D.

https://doi.org/10.1371/journal.pone.0145343.g001

Two weeks after inoculating E. grandis plantlets, uredospore pustules were visible on the host leaves inoculated with PpEucalyptus (Fig 2A), whereas the plantlets inoculated with PpGuava (Fig 2B) and the control plantlets (Fig 2C) did not exhibit any rust symptoms. The absence of pustules, chlorotic flecks or leaf deformation, demonstrate the existence of physiological variability and pathogenic specificity between the sampled populations.

[Figure omitted. See PDF.]

Fig 2. Eucalyptus infection by guava and eucalyptus rust.

Symptoms induced by inoculation of P. psidii uredospores from PpEucalyptus (A) on E. grandis variety D901. This clone is rust susceptible when grown under field conditions. The white arrows indicate the fungal pustules. The leaves are shown 15 days after inoculation.The PpGuava populations (B) and control (C) did not show typical rust symptoms.

https://doi.org/10.1371/journal.pone.0145343.g002

Overview of the P. psidii uredospore proteome

In total, 340 proteins were identified using the UniProt Puccinia genera sequences (UniProt database). Using spectral counting based label free quantification, we observed changes in relative protein abundance and a greater amount of proteins up regulated in PpGuava (S1 Table). Seventy five proteins were found exclusively in PpGuava and 29 in PpEucalyptus. Among the common proteins, 25 and 120 were increased in PpEucalyptus and PpGuava respectively (Fig 3).

[Figure omitted. See PDF.]

Fig 3. Proteins identified from P. psidii uredospore populations.

The proteins exclusively found in PpGuava (right) and PpEucalyptus (left) and common to both populations (center). Of the common proteins, 25 and 120 whose abundance were increased in PpEucalyptus and PpGuava, respectively.

https://doi.org/10.1371/journal.pone.0145343.g003

Table 1 shows some significant differences between PpGuava and PpEucalyptus protein profiles. Proteins commonly found in fungal proteome studies were identified in both fungal populations, such as ribosomal subunits, elongation factors and tubulins, as well as proteins related to energy generation, such as ATP synthase subunits.

[Figure omitted. See PDF.]

Table 1. Identification of proteins that were differentially represented in the P. psidii uredospore populations from E. grandis (PpEucalyptus) and Psidium guajava (PpGuava).

https://doi.org/10.1371/journal.pone.0145343.t001

The functional analysis showed a clear predominance of proteins related to stress response, such as heat shock proteins, catalases that were more abundant in PpGuava than in PpEucalyptus. Proteins related to fungal pathogenesis such as, citrate synthase and enolase, among other, were also unique or increased in PpGuava population. Additionally, fatty acid synthase was unique in the PpEucalyptus population.

Some proteins were found to have one or more variants, with different abundance levels. Our results suggest that there are a variety of protein homologs in the uredospores; for example, there were matches for 4 different malate dehydrogenase proteins and 2 matches for pyruvate carboxylase. Moreover, 167 proteins were described as uncharacterized proteins with different degrees of abundance (S1 Table).

Functional categorization of identified proteins

Proteins were categorized according to biological process using the Blast2GO software and major categories were present in both populations, such as cellular amino acid metabolic process, carbohydrate metabolism process, generation of precursor metabolites and energy, cellular protein modification process, signal transduction, vacuolar transport and others (Fig 4). However some GO terms, such as anatomical structure formation involved in morphogenesis, growth, mitosis, reproduction, transmembrane transport, mRNA processing, homeostatic process, mitochondrion organization, were only represented by the unique or more abundant from the PpGuava than PpEucalyptus population.

[Figure omitted. See PDF.]

Fig 4. Gene ontology of biological process terms in the proteomic analysis.

Bar graph represents the ratio of % composition of term in the proteomic data.

https://doi.org/10.1371/journal.pone.0145343.g004

This information and more details about the proteins as the Enzyme Code, Signal P and the presence of transmembrane helice domain (TMHMM), as well as proteins scores and peptide sequence are listed in S2 Table.

Discussion

P. psidii s.l. is unusual rust in that it exhibits a wide host range and is known to infect over 33 plant genera and 129 species, primarily in the Myrtaceae [37, 38, 39]. Formerly, P. psidii s.l. is conspecific and has been known as guava or eucalyptus rust [4, 40, 41, 42] due to the importance of these two hosts to the pathogen.

Despite the rust symptoms in the field being more pronounced in guava, cross-infection using PpGuava uredospores did not result in visible rust symptoms on E. grandis. However, light microscopy analysis of the PpEucalyptus and PpGuava uredospores showed the ability of both to produce germ tubes under in vitro conditions and probably having the potential to infect the host. Differences in the pathogenesis of P. psidii population’s were first reported by Maclachlan [43]. Since then, several cross-infection studies involving several hosts and P.psidii, collected from different plants of the Myrtle family, have shown clear differences in the pathogenicity of the organism [4, 13, 14, 44, 45, 46]. Aparecido et al. [14] inoculated five different hosts (Corymbia citriodora, P. guajava, S. jambos, Eugenia sp., and Eugenia uvalha) with uredospores collected from eight different host species.From this experiment, the authors described four distinct pathogenic groups. In our study, the PpGuava did not infected E. grandis, or there were no visible symptoms. It suggests the existence of a host-specific interaction, as a natural physiological variability among the uredospore populations. The differences observed in their protein profiles may be explained by the natural genetic variability within the investigated fungal populations. Our data is corroborated by recent molecular analysis of P. psidii populations that clearly proved that there is a correlation between the genetic variability of the P. psidii populations and the host from which they were collected [6].

Although rust is an economically important disease in eucalyptus and in other myrtle cultures, very little is known about the changes in the proteome of different populations of P. psidii s.l. Moreover, the monitoring of the changes in the uredospores protein profile, according to original host, is a new perspective to better understanding the plant-microbe interactions.

We evaluated two populations of the same fungal species, and observed a large number of proteins whose levels changed between the samples, as well as the presence of unique proteins in both samples. We opted to perform fungal population protein analysis based on the fact that the P. psidii has a biotrophic life style. There is no report of its growth and multiplication on synthetic media [6]. Being unable to grow outside of the host, it is difficult to characterize the proteome from other fungal structure. Therefore, it is laborious to produce spores using susceptible plants for growth and purification of the isolates. During this process, artificial selection may occur potentially producing a more aggressive pathogen depending on the host.

Some studies, which used different methods and equipments, found larger number of proteins than the 340 described for P. psidii. For instance, Bindschedler et al. [47] found 441 proteins from ungerminated spore of Blumeria graminis f. sp. Hordei. Cooper et al. [48] classified approximately 400 proteins from asexual uredospores of Uromyces appendiculatus and Liang et al. [49] observed the presence of 719 proteins from the three stages of sclerotia formation of Sclerotinia sclerotiorum. The lack of a P. psidii database was a strong challenge in the proteomic studies. The protein spectrums were matched with those from Puccinia genus from Uniprot. P. psidii has a non-complete genome sequenced. There are just few proteins sequences from this organism [50]. Many of these proteins are beta tubulins and transcription factors, which points a difficulty for protein identification of P. psidii uredospores. Even though, apart from the difficulties, we could have a good view of the protein roles in the uredospores between two contrasting P. psidii populations.

Using uredospores from populations that have been previously described as different biotypes [13], we identified proteins that are potentially related to the degree of variability and aggressiveness of the P. psidii allowing us to make some inferences regarding fungal diversity. A larger number of unique and/or abundant proteins were found in the PpGuava uredospores; this is most likely to be related to the higher genetic diversity within PpGuava, resulting in a broader range of proteins produced by this population, as well as a larger number of possible isoforms, which may be justified by the probable higher number of alleles in this population. These data are in agreement with a recent study that compared rust populations from eucalyptus and other myrtle species in Brazil and revealed a higher genetic diversity in the guava populations [6].The authors suggested that the pathogen recently migrated from a native South American host. When the fungus moved to other plant species, such as the exotic Eucalyptus spp., a new selection process was introduced by the new host [6]. Graça et al. [39] did not prove that P. psidii jumped from guava to eucalyptus host, but the existence of host-associated biotypes of P. psidii in Brazil, indicates that this diversity must be considered. Morin et al. [51] investigated the host-range of the P. psidii s.l. present in Australia. The authors concluded that the development of the disease in one host species in a tribe did not mean that a related species or genus also developed the disease. Moreover, they demonstrated the inability of the Australian rust accession to develop uredinia on P. guajava compared to the accession from Florida [52], supporting the hypothesis of physiological diversity between these populations.

Interestingly, during the uredospore phase, we identified many proteins related to fungal pathogenicity, such as peptidases, proteases, and proteins that are able to modify host factors [53, 54, 55, 56, 57], that were unique or more abundant only in the PpGuava than in PpEucalyptus population, for instance, proteassomes. Proteasomes are multi-subunit and multicatalytic proteolytic complexes that are found in the cells of most organisms, and these complexes play important roles in protein turnover in both the cytosol and the nucleus. They selectively degrade intracellular proteins in eukaryotic cells, particularly via the ubiquitin-proteasome pathway [58], and are involved in multiple specific functions, such as cell cycle regulation, cell signaling, the selective elimination of abnormal proteins, and the flux of substrates through metabolic pathways [59, 60]. As previously suggested by Kim et al.[61], the wax coated hydrophobic surfaces, which naturally occur on the host surface, may activate the 20S proteasome to mobilize the storage proteins that are present in uredospores. Chaperones are responsible for the initial folding and maintenance of proteins structure compared to the conventional heat shock proteins (HSPs), which are present only under conditions of high temperature or stress. Some studies suggested that these proteins also act to prevent protein aggregation, thus controlling protein turnover, which will increase its longevity and virulence [62]. Hodgetts et al. [63] observed that the up-regulation of HSP90 in Saccharomyces cerevisiae resulted in an increase of virulence in infected mice. We identified many HSP classes, some with peptide signals for secretion, similar to the data reported by Song et al [30]. The authors found high level of HSP in P. triticina haustorium. All identified HSPs were less abundant in the PpEucalyptus than in PpGuava population, most likely because the uredospores from the guava fruit are under more stressfull conditions [47] than those on PpEucalyptus, suggesting a differential response to the host organism. Larger amounts of HSPs were also identified in the Uromyces appendiculatus uredospores, which cause bean rust [48]. More investigation into the expression of HSPs is needed to better understand the roles of these proteins in the infection process during the P. psidii s.l. life cycle.

Tubulin and actin proteins were also increased in abundance in PpGuava. Similarly, Mandelc et al. [32] also observed that these proteins were increased in abundance in the lethal phytotype of V. albo-atrum. Tubulins and actins are the major proteins involved in the formation of the microtubules and microfilaments. GTPases regulates cellular polarization and motility, which affects the cytoskeleton, and are involved in the regulation of the secretion process and endocytosis [64]. The enolases identified in this study were also higher in abundance in the PpGuava population. Enolases participate in other cellular activities and can even act as HSPs; they can bind to cytoskeletal and chromatin structures to modulate transcription, thus playing a crucial role in pathophysiological processes [65, 66]. Enolase proteins increased it abundance in a virulent isolate of Pyrenophora tritici-repentis, the causal agent of tan spot in wheat [67]. Another example is the increased abundance of malate dehydrogenase in PpGuava. This protein is related to fungal virulence, the secretion of this compound acidifies the environment, signaling the expression and secretion of virulence factors and phytotoxins [17, 32].

Another interesting protein class, not previously related to fungal virulence, is prohibitin. This protein, with homologues in animals, fungi, plants, and unicellular eukaryotes was first reported as a negative regulator of cell proliferation. However, prohibitin is expressed as two transcripts with varying lengths of the 3' untranslated region. The longer one is present at higher levels in proliferating cells, suggesting that this longer 3' untranslated region sequence may function as a trans-acting regulatory RNA [68, 69].

Acetyl CoA carboxylase biotin carboxylase is abiotin-dependent enzyme involved in the biosynthesis, as well as, in the metabolism of fatty acids [70, 71]. Fatty acid synthase catalyzes the formation of long-chain fatty acids from acetyl-CoA, malonyl-CoA and NADPH [72]. Both were found as PpEucalyptus unique proteins., maybe in PpGuava these proteins may be present in lower abundance.

Our data also support the model proposed for U. appendiculatus by Cooper et al. [48], which suggests that the uredospore exists in a suspended translational state. Uredospores can produce proteins very quickly without having to assemble de novo translational components. Moreover, the peptide chains dispersed in the ribosome are protected by heat shock proteins and can resist adverse conformational changes caused by desiccation or high temperatures that the uredospores might be subjected to, once released into the environment [48].Corroborating this hypothesis, we identified a large number of ribosomal proteins subunits with different level of abundance, which may be related to the beginning of the infection process by the pathogen.

Additionally, categorization by protein function showed increased abundance of proteins related to translation, ribosome biogenesis and protein folding from the PpEucalyptus sample. Response to stress, generation of precursor metabolites and energy and carbohydrate metabolic process were more abundant in PpGuava. The proteins involved in signal transduction were present in similar proportion in both populations.

In turn, uredospores populations have specific reservatory of proteins according to the host from supporting the physiological variability of these two populations. The virulent proteins from PpEucalyptus were less abundant when compared to PpGuava. In guava, close to 100% of the fruits can be lost due to severe infections most likely because the uredospores reservatory of proteins of this population; related with virulence and stress response. Likely, PpEucalyptus uredospores need to process of proteins of virulence to start the host infection process, its reservatory is composed by biogenesis, protein folding and translation (Fig 5). It is commonly observed that uredopores from a specific host are less aggressive when infecting a host from other species. The aggressiveness increases after successive cycles of infection in the new host. One hypothesis is the modulation of the expression of some genes according to the host, resulting in different protein reservoirs in each population.

[Figure omitted. See PDF.]

Fig 5. Protein profiles in PpGuava and PpEucalyptus and the correlation with their physiological variability.

While proteins correlated to fungal virulence and stress response had the abundance increased in PpGuava, proteins related to biogenesis, protein folding and translation had the abundance increased in PpEucalyptus.

https://doi.org/10.1371/journal.pone.0145343.g005

The proteomic information generated in the present work contributes to a better understanding of the mechanisms underlying the pathogenicity of P. psidii s.l. and its interaction with its host organism. More comprehensive proteomic analyses of the P. psidii-host interaction, combined with other studies, will be important for elucidating the molecular basis of the interactions in this intriguing pathosystem.

Supporting Information

[Figure omitted. See PDF.]

S1 Table. Complete set of proteins identified using the PLGS 2.5.1 expression analysis.

Access ID is the access identification of the protein in Uniprot database (http://www.uniprot.org/ release Version 2015_02). The score of protein expression analysis, ratio and Log (PpGuava/PpEucalyptus) of quantified proteins, standard deviation (SD) of the PpGuava/PpEucalyptus ratio log and their p-values of the Log of the PpGuava/PpEucalyptus ratio obtained from PLGS 2.5.1. A p value ≤ 0.05 identifies proteins whose abundance decrease, whereas a p value ≥ 0.95 identifies proteins whose abundance increased in PpGuava; p value identifies proteins that were unique to either the PpEucalyptus or PpGuava P. psidii uredospore populations isolated from eucalyptus and guava, respectively (XLS).

https://doi.org/10.1371/journal.pone.0145343.s001

(XLSX)

S2 Table. Additional information of set of proteins.

Organism source, GO biological process, Enzyme Code, Signal P, Trans-membrane helice domain (TMHMM), as well as protein score, average amount (fmol) of protein, average score of proteins and peptide sequences and average amount of matched peptides to each protein.

https://doi.org/10.1371/journal.pone.0145343.s002

(XLSX)

Acknowledgments

We thank Dr. Esteban R. Gonzalez (FuturaGene Corporation) for providing the uredospores from P. psidii populations.

Author Contributions

Conceived and designed the experiments: MCQ TFL CAL. Performed the experiments: MCQ TFL FGM APB SGG. Analyzed the data: MCQ TR APB IGFB. Contributed reagents/materials/analysis tools: LMF DHM TR. Wrote the paper: MCQ MTVL CAL.

Citation: Quecine MC, Leite TF, Bini AP, Regiani T, Franceschini LM, Budzinski IGF, et al. (2016) Label-Free Quantitative Proteomic Analysis of Puccinia psidii Uredospores Reveals Differences of Fungal Populations Infecting Eucalyptus and Guava. PLoS ONE 11(1): e0145343. https://doi.org/10.1371/journal.pone.0145343

References

1. Winter G. Repertorium. Rubenhorstii fungi europaei et extraeuropaei.Cent XXXI et XXXII. Hedwigin. 1884; 23, 164–172.

2. Ferrari JT, Nogueira EMC, Santos AJT. Control of rust (Puccinia psidii) in guava (Psidium guajava). Acta Hort. 1997; 452: 55–57.

3. Gonçalves S. Preliminary list of plant diseases in Espirito Santo State. Rio de Janeiro, Ministério da Agricultura, 1929; 1–12.

4. Ferreira FA. Ferrugem do Eucalipto. Revista Árvore. 1983; 7: 91–109.

5. Glen M, Alfenas AC, Zauza EAV, Wingfield MJ, Mohammed C. Puccinia psidii: A threat to the Australian environment and economy–A review. Australas Plant Path. 2007; 36: 1–16.

6. Quecine MC, Bini AP, Romagnoli EM, Andreote FD, Moon DH, Labate CA. Genetic variability of Puccinia psidii populations revealed by PCR-DGGE and T-RFLP markers. Plant Dis. 2014; 98: 16–23.

7. Simpson JA, Thomas K, Grgurinovic CA. Uredinales species pathogenic on species of Myrtaceae. Australas Plant Pathol. 2006; 35: 549–562.

8. Maier W, Begerow D, Weis M, Oberwinkler F. Phylogeny of the rust fungi: an approach using the nuclear large subunit ribosomal DNA sequences. Can J Bot. 2003; 81: 12–23.

9. Walker J. Pacific mycogeography: deficiencies and irregularities in the distribution of plant parasitic fungi. Aust J Bot. 1983; 10: 89–136.

10. Roux J, Greyling I, Coutinho TA, Verleur M, Wingfield MJ. The myrtle rust pathogen, Puccinia psidii, discovered in Africa. IMA Fungus. 2013; 4: 155–159. pmid:23898420

11. Figueiredo MB. Life cycle and ecology of Puccinia psidii. O Biol. 2001; 63: 69–71.

12. Wiethölter N, Horn S, Reisige K, Beike U, Moerschbacher BM. In vitro differentiation of haustorial mother cells of the wheat stem rust fungus, Puccinia graminis f. sp. tritici, triggered by the synergistic action of chemical and physical signals.Fungal Genet Biol. 2003; 38: 320–326. pmid:12684021

13. Coelho L, Alfenas AC, Ferreira FA. Variabilidade fisiológica de Puccinia psidii—Ferrugem do Eucalipto. Summa Phytopathol. 2001; 27: 295–300.

14. Aparecido CC, Figueiredo MB, Furtado EL. Groups of physiological variability in Puccinia psidii populations. Summa Phytopathol. 2003; 29: 234–238.

15. Gokce E, Franck WL, Oh Y, Dean RA, Muddiman DC. In-Depth analysis of the Magnaporthe oryzae conidial proteome. J. Proteome Res. 2012; 11: 5827−5835. pmid:23039028

16. Kim Y, Nandakumar MP, Marten MR. Proteomics of filamentous fungi. Trends Biotechnol. 2007; 25: 395–400. pmid:17681627

17. Doyle S. Fungal proteomics: from identification to function. FEMS Microbiol Lett. 2011; 321: 1–9. pmid:21517945

18. Ebstrup T, Saalbach G, Egsgaard HA. Proteomics study of In vitro germination and appressoria formation in Phytophthora infestans. Proteomics. 2005; 5: 2839–2848. pmid:15996011

19. Fernández-Acero FJ, Jorge I, Calvo E, Vallejo I, Carbú M, Camafeita E, et al. Proteomic analysis of phytopathogenic fungus Botrytis cinerea as a potential tool for identifying pathogenicity factors, therapeutic targets and for basic research. Arch Microbiol. 2007; 187: 207–215. pmid:17124592

20. Yajima W, Kav NNV. The proteome of the phytopathogenic fungus Sclerotinia sclerotiorum. Proteomics. 2006; 6: 5995–6007. pmid:17051649

21. Paper JM, Scott-Craig JS, Adhikari ND, Cuomo CA, Walton JD. Comparative proteomics of extracellular proteins in vitro and in planta from the pathogenic fungus Fusarium graminearum. Proteomics. 2007; 7: 3171–3183. pmid:17676664

22. Schmitt S, Prokisch H, Schlunck T, Camp DG, Ahting U, Waizenegger T, et al. Proteome analysis of mitochondrial outer membrane from Neurospora crassa. Proteomics. 2006; 6: 72–80. pmid:16294304

23. Kniemeyer O, Lessing F, Scheibner O, Hertweck C, Brakhage AA. Optimisation of a 2-D gel electrophoresis protocol for the human-pathogenic fungus Aspergillus fumigatus. Curr Genet. 2006; 49: 178–189. pmid:16362819

24. Melin P, Schnürer J, Wagner EG. Proteome analysis of Aspergillus nidulans reveals proteins associated with the response to the antibiotic concanamycin a, produced by Streptomyces species. Mol. Genet. Genomics. 2002; 267: 695–702. pmid:12207217

25. Carberry S, Neville CM, Kavanagh KA, Doyle S. Analysis of major intracellular proteins of Aspergillus fumigatus by MALDI mass spectrometry: identification and characterization of an elongation factor 1b protein with glutathione transferase activity. Biochem Biophys Res Commun. 2006; 341: 1096–1104. pmid:16455047

26. Oda K, Kakizono D, Yamada O, Iefuji H, Akita O, Iwashita K. proteomic analysis of extracellular proteins from Aspergillus oryzae grown under submerged and solid-state culture conditions. Appl. Environ. Microbiol. 2006; 72: 3448–3457. pmid:16672490

27. Qin G, Tian S, Chan Z, Li B. Crucial role of antioxidant proteins and hydrolytic enzymes in pathogenicity of Penicillium expansum: analysis based on proteomic approach. Mol. Cell. Proteomics. 2007; 6: 425–438. pmid:17194899

28. Zorn H, Peters T, Nimtz M, Berger RG. The secretome of Pleurotus sapidus. Proteomics. 2005; 5: 4832–4838. pmid:16281184

29. Shimizu M, Wariishi H. Development of a sample preparation method for fungal proteomics. FEMS Microbiol Lett. 2005; 247: 17–22. pmid:15927743

30. Song X, Rampitsch C, Soltani B, Mauthe W, Linning R, Banks T, et al. Proteome analysis of wheat leaf rust fungus, Puccinia triticina, infection structures enriched for haustoria. Proteomics. 2011; 11: 944–963. pmid:21280219

31. Godfrey D, Zhang Z, Saalbach G, Thordal-Christensen H. A proteomics study of barley powdery mildew haustoria. Proteomics. 2009; 9: 3222–3232. pmid:19562796

32. Mandelc S, Radisek S, Jamnik P, Javornik B. Comparison of mycelial proteomes of two Verticillium albo-atrum pathotypes from hop. Eur J Plant Pathol. 2009;125: 159–171.

33. Damerval C, Vienne D, Zivy M Thiellement H. Technical improvements in two-dimensional electrophoresis increase the level of genetic variation detected in wheat-seedling protein. Electrophoresis. 1986; 7: 53–54.

34. Silva JC, Gorenstein MV, Li GZ, Vissers JP, Geromanos SJ. Absolute quantification of proteins by LCMSE: a virtue of parallel MS acquisition. Mol Cell Proteomics. 2006; 5: 144–156. pmid:16219938

35. Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, et al. Gene Ontology: tool for the unification of biology. the geneontology consortium. Nat Genet. 2000; 25: 25–29. pmid:10802651

36. 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. pmid:16081474

37. Carnegie AJ, Lidbetter JR. Rapidly expanding host range for Puccinia psidii sensu lato in Australia. Australas Plant Pathol. 2012; 41: 13–29.

38. Tommerup IC, Alfenas AC, Old KM. Guava rust in Brazil–a threat to Eucalyptus and other Myrtaceae. New Zealand J Forest Sci. 2003; 33: 420–428.

39. Graça RN, Ross-Davis AL, Klopfenstein NB, Kim MS, Peever TL, Cannon PG, et al. Rust disease of eucalypts, caused by Puccinia psidii, did not originate via host jump from guava in Brazil. Mol Ecol. 2013; 22: 6033–6047. pmid:24112757

40. Joffily J. Ferrugem do eucalipto. Bragantia. 1944; 4: 475–487.

41. Coutinho TA, Wingfield MJ, Alfenas AC, Crous PW. Eucalyptus rust: A disease with the potential for serious international implications. Plant Dis. 1998; 82: 819–825.

42. Kriticos DJ, Morin L, Leriche A, Anderson RC, Caley P. Combining a climatic niche model of an invasive fungus with its host species distributions to identify risks to natural as sets: Puccinia psidii Sensu Lato in Australia. PLoS One. 2013; 8:e64479. pmid:23704988

43. Maclachlan JD. A Rust of the pimento tree in jamaica, BWI. Phytopathol. 1938; 28: 157–170.

44. Marlatt RB, Kimbrough JW. Rust (Puccinia psidii) of Allspice (Pimenta dioica) Appears in Florida. P Fl St Hort Soc. 1980; 93: 111.

45. Castro HA, Krugner TL, Idehira CHF, Cappello MSC, Marchi AB. Inoculação cruzada de eucalyptus, goiaba (Psidium Guajava) e jambeiro (Sizygium Jambos) com Puccinia psidii. Fitopatol Bras. 1983; 8: 491–497.

46. Castro HA, Bergamim FA, Krugner TL. Padrão de produção de uredosporos em mudas de eucalyptus spp inoculadas artificialmente com Puccinia psidii. Summa Phytopathol. 1984; 10: 155–170.

47. Bindschedler LV, Burgis TA, Mills DJ, Ho JTC, Cramer R, Spanus PD. In planta proteomics and proteogenomics of the biotrophic barley fungal pathogen Blumeria graminis f. sp. hordei. Mol Cell Proteomics. 2009; 8: 2368–2381. pmid:19602707

48. Cooper B, Garrett WM, Campbell KB. Shotgun identification of proteins from uredospores of the bean rust Uromyces appendiculatus. Proteomics. 2006; 6: 2477–2484. pmid:16518873

49. Liang Y, Rahman MH, Strelkov SE, Kav NNV. Developmentally induced changes in the sclerotial proteome of Sclerotinia sclerotiorum. Fungal Biol. 2010; 114: 619–627. pmid:20943173

50. Tan MK, Collins D, Chen Z, Englezou A, Wilkins MR. A brief overview of the size and composition of the myrtle rust genome and its taxonomic status. Mycology. 2014; 5: 52–63. pmid:24999437

51. Morin L, Aveyard R, Lidbetter JR, Wilson PG. Investigating the host-range of the rust fungus Puccinia psidii sensu lato across tribes of the family myrtaceae present in Australia. PLoS One. 2012; 7:e35434. pmid:22523596

52. Rayachhetry MB, Van TK, Center TD, Elliott ML. Host range of Puccinia psidii, a potential biological control agent of Melaleuca quinquenervia in Florida. Biol Control. 2001; 22: 38–45.

53. Kim HS, Desveaux D, Singer AU, Patel P, Sondek J, Dangl JL. The Pseudomonas syringae effector AvrRpt2 cleaves its C-terminally acylated target, RIN4, from Arabidopsis membranes to block RPM1 activation. P Natl Acad Sci USA. 2005; 102: 6496–6501.

54. Xia Y. Proteases in pathogenesis and plant defence. Cel. Microbiol. 2004; 6: 905–913.

55. Roden J, Eardley L, Hotson A, Cao Y, Mudgett MB. Characterization of the Xanthomonas AvrXv4 effector, a SUMO protease translocated into plant cells. Mol Plant Microbe In. 2004; 17: 633–643.

56. Jia Y, McAdams SA, Bryan GT, Hershey HP, Valent B. Direct interaction of resistance gene and avirulence gene products confers rice blast resistance. Embo J. 2000; 19: 4004–4014. pmid:10921881

57. Debono M, Goedee RS. Antibiotics that inhibit fungal cell wall development. Annu Rev Microbiol. 1994; 48: 471–497. pmid:7826015

58. Baumeister W, Walz J, Zühl F, Seemüller E. The proteasome: paradigm of a self-compartmentalizing protease. Cell. 1998; 92: 367–380. pmid:9476896

59. Monaco JJ, Nandi D. The genetics of proteasomes and antigen processing. Annu Rev Genet. 1995; 29: 729–754. pmid:8825492

60. Coux O, Tanaka K, Goldberg AL. Structure and functions of the 20S and 26S proteasomes. Annu Rev Biochem. 1996; 65: 801–847. pmid:8811196

61. Kim ST, Yu S, Kim SG, Kim HJ, Kang SY, Hwang DH, et al. Proteome analysis of rice blast fungus (Magnaporthe grisea) proteome during appressorium formation. Proteomics. 2004; 4: 3579–3587. pmid:15378734

62. Burnie JP, Carter TL, Hodgetts SJ, Matthews RC. Fungal heat-shock proteins in human disease. FEMS Microbiol Rev. 2006; 30: 53–88. pmid:16438680

63. Hodgetts S, Matthews R, Morrissey G, Mitsutake K, Piper P, Burnie J. Over-expression of Saccharomyces cerevisiae hsp90 enhances the virulence of this yeast in mice. FEMS Immunol Med Microbiol. 1996; 16: 229–234. pmid:9116640

64. Ridley AJ. Rho GTPases and actin dynamics in membrane protrusions and vesicle trafficking. Trends Cell Biol. 2006; 16: 522–529. pmid:16949823

65. Pancholi V. Multifunctional alpha-enolase: its role in diseases. Cell Mol Life Sci. 2001; 58: 902–920. pmid:11497239

66. Pancholi V, Chhatwal GS. Housekeeping enzymes as virulence factors for pathogens. Int J Med Microbiol. 2003; 293: 391–401. pmid:14760970

67. Cao T, Kim YM, Kav NN, Strelkov SE. A proteomic evaluation of Pyrenophora tritici-repentis, causal agent of tan spot of wheat, reveals major differences between virulent and avirulent isolates. Proteomics. 2009; 9: 1177–1196. pmid:19206107

68. Van Aken O, Pecenkova T, van de Cotte B, De Rycke R, Eeckhout D, Fromm H, et al. Mitochondrial type-I prohibitins of Arabidopsis thaliana are required for supporting proficient meristem development. Plant J. 2007; 52: 850–864. pmid:17883375

69. Mishra S, Murphy LC, Murphy LJ. The prohibitins: emerging roles in diverse functions. J Cell Mol Med. 2006; 10: 353–363. pmid:16796804

70. Tong L. Acetyl-coenzyme A carboxylase: crucial metabolic enzyme and attractive target for drug discovery. Cell. Mol. Life Sci. 2005; 62: 1784–803. pmid:15968460

71. Brownsey RW, Zhande R, Boone AN. Isoforms of acetyl-CoA carboxylase: structures, regulatory properties and metabolic functions. Biochem Soc Trans. 1997; 25: 1232–1238. pmid:9449982

72. Niwa H, Katayama E, Yanagida M, Morikawa K. Cloning of the fatty acid synthetase beta subunit from fission yeast, coexpression with the alpha subunit, and purification of the intact multifunctional enzyme complex. Protein Expr Purif.1998; 13:403–413. pmid:9693066