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

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-MSE 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.


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.
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 4 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 4 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 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 Rapi-Gest. 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] 2+ ) was used for initial single-point calibration and MS/MS fragment ions of GFP [Glu 1]-Fibrinopeptide B with a [M + 2H] 2+ = 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: p0.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/).

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).
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.
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). 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.  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.
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.
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 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). (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. (XLSX)