BACKGROUND
huanglongbing disease currently threatens areas where citrus cultivation is important in the agricultural economy such as east asia, the middle east, and the americas.
huanglongbing disease is caused by three species of candidatus liberibacter asiaticus , americanus and africanus  <cit> .
the pathogen is transmitted by two species of psyllids: diaphorina citri and trioza erytreae.
recently, trioza erytreae was found for the first time in europe .
typical symptoms of huanglongbing disease in leaves include shoot yellowing and blotchy, mottled leaves.
although most of the fruits are still of commercial quality, fruits from severely affected branches are unmarketable: small, lopsided, green, and acidic, with many aborted seeds.
leaves accumulate starch, phloem is damaged and cell wall lamellae swell during calas infection  <cit> .
candidatus liberibacter spp.
belong to the alpha subdivision of the proteobacteria based on ribosomal region sequence data  <cit> .
the bacterium has not yet been definitively cultured despite attempts to do so  <cit> .
koch’s postulates have not been fulfilled for this disease, so possible interactions with other microrganisms cannot be ruled out.
the pathogen lives in the insect and in the phloem of citrus trees.
once acquired, it typically persists for the rest of the life of the host.
insecticides can decrease psyllid populations, but since the pathogen remains in the vector, disease spread can occur with the presence of just a few infected psyllids in the orchard.
all genotypes within the genus citrus are susceptible to hlb to varying degrees although species of other close-related genera showed some sort of resistance  <cit> .
there is variability in disease severity and symptoms among citrus genotypes  <cit> .
murraya paniculata , an ornamental citrus closely-related plant, showed fewer symptoms of the disease  <cit> .
a study examining the responses of  <dig> genotypes to hlb disease grouped them based on phenotypic analysis of induced symptoms  <cit> .
another recent study has evaluated  <dig> citrus accessions and  <dig> accessions belonging to other closely related genera.
resistance was reported in accessions not belonging to citrus genera  <cit> .
another work have screened citrus germplasm susceptibility to hlb analyzing sixteen citrus genotypes  <cit> .
results showed that citrus macrophylla and c. medica were the most susceptible while complex genetic hybrids ‘us 1-4-59’ and ‘fallglo’ were the least susceptible.
a metabolomic investigation was conducted comparing five different tolerant hybrids and a highly susceptible cultivar to identify potential metabolites linked with diverse response  <cit> .
the causes of the disease have been studied using different “omic” approaches to identify which genes, proteins and metabolites may be targeted by innovative diagnostic and therapeutic methods.
the genome of the pathogen was sequenced using a metagenomic approach, both from infected plants  <cit>  and the insect vector  <cit> .
no toxins or other secreted proteins have been linked to the disease and the mechanisms of its pathology are still unclear.
large scale microarray analysis revealed significant modulation of genes involved in transport, cell defense and carbohydrate metabolism  <cit> .
photosynthesis is diminished in both young and mature leaves, but it is upregulated in infected fruits  <cit> .
starch accumulation was linked to the upregulation of genes involved in glucose import into the chloroplast and starch biosynthesis  <cit> .
a modulated jasmonic -salicilic acid  crosstalk of innate responses may lead to a misdirected defence response.
an integrated approach of 2-de and mass spectrometry showed that changes in levels of several proteins involved in photosynthesis and protein synthesis were linked to reduced concentrations of ca, mg, fe, zn, mn and cu in infected grapefruit leaves  <cit> .
proteins upregulated in infected samples were involved in redox stage and cell defense such as cu/zn superoxide dismutase, peroxidases, chitinases and lectin-related proteins  <cit> .
‘madam vinous’ sweet orange plants infected by calas showed increased miraculin-like proteins, chitinase, cu/zn superoxide dismutase and lipoxygenase.
some key metabolites modulated by hlb include proline, β-elemene, -trans-caryophyllene, and α-humulene  <cit> .
increased accumulation of some amino acids  and organic acids was linked to greater susceptibility of ‘madam vinous’ sweet orange compared to carrizo citrange  <cit> .
however, it is worthy to notice that these trees were highly infected with many secondary effects so it will be necessary to confirm these results with newly infected trees.
an increased amount of most amino acids, involved in plant defense to pathogens was observed in tolerant varieties in such as phenylalanine, tyrosine, tryptophan, lysine, and asparagine  <cit> .
this study examines proteomic changes in fully photosynthesizing leaves to determine how disease mechanisms and susceptibility vary between two citrus genotypes, using an integrated approach of principal component analysis , gene set and pathway enrichment analysis.
the purpose was to characterize key proteins and post-transcriptionally modulated pathways in different citrus genotypes at a late symptomatic stage of hlb.
RESULTS
navel orange ) is an hlb-sensitive cultivar while volkameriana is moderately tolerant  <cit> .
different techniques  have been used to identify key proteins differentially regulated by hlb in citrus leaves  <cit> .
over  <dig> proteins were analyzed using isobaric tags for relative and absolute quantitation  for both genotypes .
in navel orange,  <dig> proteins were differentially regulated between infected and healthy tissue  .
in volkameriana,  <dig> differentially regulated proteins were found between infected and healthy tissue .
pca analysis
pca was used to visualize differences between the four analyzed genotypes x disease status, subdividing the entire proteomic profile into three important subcategories: biotic stress responses, overall cell metabolism and transcriptional regulation pathways .
all identified proteins belonging to these important functional categories were used for the pca plots.
the four categories of samples  were clearly separated in all three pca plots, implying significant protein changes in all three gene categories related to species and health status.
for biotic stress-related proteins, pc <dig> and pc <dig> accounted for  <dig> and 24 % of the data variability, respectively.
important proteins associated with each functional category that contributed to separation between sample types  are listed in fig.
<dig>  some key proteins involved in redox state significantly contributed in the separation of infected v from the rest of the other categories.
in the general metabolism pca plot, the pc <dig> and pc <dig> accounted for  <dig> and 24 % of the data variability, respectively.
key proteins involved in primary metabolism associated with the separation of healthy v from the infected v include malate dehydrogenase and pyruvate dehydrogenase , sucrose synthase and agpase .
the third pca plot was generated based on the expression of proteins involved in transcriptional regulation, signaling, hormone and redox state.
pc <dig> and pc <dig> accounted for  <dig> and 27 % of data variability, respectively.
interestingly, the regulation of few important proteins seems to specifically characterize the infected vs. healthy state of navel orange.
these include map kinase  <dig>  udp-glucosyl transferase, and aspartyl protease.
no changes were observed for mapk <dig>  proteins that appear in the three pca plots were highly regulated in the comparison between infected and control in both citrus genotypes.
indeed, they may be considered as putative candidate biomarkers of a clear symptomatic status in citrus at proteomic level.
their hlb-regulated pattern of expression greatly contributed in distinguishing the four different leaf sample types.
further analysis will be needed to validate these data and confirm the role in the pathogenesis of hlb disease.fig.
<dig> principal component analysis of differentially regulated proteins of four types of leaf tissues .
proteins that contribute highly to the separation of the the four samples are numbered and listed next to each graph
gene set and pathway enrichment analysis
the pageman web tool highlighted which gene categories were up- or down-regulated in each pairwise comparison .
both infected genotypes exhibited repressed amino acid biosynthesis and protein synthesis.
in navel orange, photosynthesis, isoflavone pathways, tetrapyrrole synthesis, and galactose metabolism were significantly inhibited by hlb.
in volkameriana, s-assimilation, isoprenoids, rna binding and amino acid activation were specifically diminished.
in both genotypes, hlb enhanced starch-related pathways, biotic stress-related proteins, beta  <dig>  glucan hydrolases, and protein degradation pathways.
some distinct differences were observed between genotypes in hlb response.
in volkameriana, cell wall modifications, galactose metabolism, and heat shock proteins were upregulated.
in navel orange, amino acid degradation, lipid metabolism, jasmonates, and pr-proteins were upregulated.fig.
<dig> gene set enrichment analysis using pageman web-tool.
upregulated and downregulated pathways at proteomic level in infected volkameriana  and infected navel orange  in comparison to respective healthy controls
a pathway enrichment analysis was performed using the david bioinformatic resource to determine which metabolic pathways were commonly or specifically affected by hlb disease in both genotypes .
some metabolic pathways were altered by hlb in both species.
amino acid metabolism  was significantly downregulated.
other key inhibited pathways include biosynthesis of plant hormones, terpenoids, and phenylpropanoids.
on the other hand, tyrosine metabolism was upregulated in both genotypes.
in navel orange, fatty acid biosynthesis and nitrogen metabolism were diminished while alpha-linolenic acid metabolism was enhanced.
in volkameriana, alkaloids and pyruvate metabolism-related proteins were repressed while galactose metabolism and fatty acid metabolism were upregulated.table  <dig> pathway enrichment analysis using david bioinformatics resources  <dig> .
pathways that were upregulated and downregulated for the healthy/infected comparison for each genotype are indicated with the corresponding p-value
primary metabolism
the integration of the two pairwise comparisons into a unified mapman visualization allowed us to identify proteins that were commonly or specifically regulated in response to hlb in the two species.
some proteins involved in cell wall modifications were upregulated by hlb only in volkameriana: expansin a <dig>  expansin a <dig>  expansin-like b <dig>  xyloglucan endotransglycosylase  <dig>  and xyloglucanxyloglucosyltransferase  .
in navel orange, proteins involved in fatty acid biosynthesis and elongation were generally repressed while several key proteins involved in amino acid degradation were upregulated: arginase, pyrroline-5-carboxylase reductase, lactoylglutathione lyase, and 3-hydroxylmethylglutaryl-coa lyase.
the increase of protein degradation indicates that senescence processes may be more highly activated in hlb-diseased navel orange than in volkameriana.
sucrose metabolism was only slightly affected by hlb at the protein level; only sucrose synthase was repressed.
starch metabolism was more altered .
the first enzyme of starch biosynthesis, adp-glucose pyrophosphorylase, was inhibited by hlb in both genotypes.
in volkameriana, starch synthase was upregulated while  <dig> -alpha-glucan starch branching enzyme was slightly downregulated.
among starch degradation enzymes, glucan phosphorylase and heterogycan glucosidase  <dig> were upregulated in both genotypes.
alpha-amylase was upregulated in volkameriana while beta-amylase  <dig> was enhanced in navel orange.
taken together these findings showed that starch metabolism was highly affected in both genotypes at the protein level.
no clear association between differing susceptibility to hlb is evident from starch pathway regulation alone.fig.
<dig> proteins involved in primary metabolic responses to hlb.
a) sucrose and starch metabolism; b) raffinose metabolism; c) tca cycle; d) oxidative pentose phosphate.
proteins that were differentially expressed between control and infected trees are indicated by colors, based on their pattern of expression in the two genotypes .
each colored square represents expression change in a protein catalyzing a step in the pathway.
more than one square grouped together indicates different members of the same protein family found to be differentially regulated
a significant downregulation in hlb-infected samples was observed for proteins involved in photosynthetic reactions  <cit> .
the altered transcription of sugar and starch metabolism genes caused by hlb  <cit>  mostly agreed with the corresponding protein changes presented in the present work.
the observed changes in starch-related pathways were consistent with the transcriptomic analysis  <cit> .
increased starch degradation was observed probably due to the increased starch concentrations in infected leaf tissues.
adp glucose-pyrophospholyase  was repressed in both infected species.
this enzyme is rate-limiting for starch biosynthesis, catalyzing the conversion of glucose-1-phosphate to adp-glucose that is polymerized into amylopectin or amylose  <cit> .
starch accumulation in infected leaves is a typical symptom of hlb  <cit> .
however, the greatest occurrence of this process may occur at an early, asymptomatic stage.
indeed we may speculate that the accumulation of starch may be a secondary effect of the disease instead of being the cause of symptoms.
when symptoms are already evident and yellowing is present, starch biosynthesis may slow down and starch degradation is expected to be activated as a response of the plant to limit damage to cell structures.
however, some differences in primary metabolism were observed between transcriptomic and proteomic approaches.
changes in invertase gene expression were observed  <cit> , but not its protein levels.
sucrose, the substrate of invertase, is produced by photosynthesis in leaves and then transported to sink tissues  through the phloem.
involvement of sucrose in signaling of innate responses has been described recently  <cit> .
invertase plays an important role in plant stress responses, possibly serving as an extracellular signal of pathogen attacks  <cit> .
it is possible that expression of this gene in response to hlb depends on genotype, plant physiological conditions and age, tissue developmental conditions, type of infection  and/or environmental conditions .
taken together all these primary metabolism results concur with previous findings showing significant modification of transcript abundance in minor carbohydrate metabolism  <cit>  but do not suggest a clear link to the well-known difference in tolerance between the two genotypes.
gpt <dig> has been linked with hlb disease: this gene is responsible for glucose import into the chloroplast and consequently for starch accumulation  <cit> .
the protein was not found among those extracted and characterized by itraq, therefore no conclusions can be made about changes in protein levels due to hlb.
however, the protocol used to analyze the citrus proteome favored detection of soluble cell proteins over membrane proteins such as gpt <dig>
expression of key proteins involved in the tca cycle and pep metabolism were affected by hlb.
atp-citrate lyase subunit b <dig> was repressed in both species .
alcohol dehydrogenase  <dig> which converts aceltaldehyde to ethanol, was upregulated in both genotypes.
in infected navel orange, isopropyl malate isomerase, which converts citrate to isocitrate and pyruvate decarboxylase involved in fermentation, was more abundant than in healthy tissue.
in volkameriana, pyruvate decarboxylase was repressed while aconitate hydratase involved in tca cycle was upregulated.
raffinose metabolism was drastically altered by the disease .
in navel orange, expression of alpha-galactosidase  <dig>  udp-glucose-4-epimerase, glucose-6-phosphate isomerase, and alpha-galctosidase  <dig> were enhanced by hlb.
in volkameriana, raffinose synthase, phosphofructokinase  <dig> and phosphoglucomutase were upregulated.
sucrose synthase  <dig> and phosphoglycerate mutase were donwregulated in both genotypes.
a transketolase that converts the xylulose-5-p in sedoheptulose-7-p was repressed in infected navel orange .
significant repression of aspartate biosynthesis and serine metabolism was observed in infected navel orange leaves.
this agrees with a previously described downregulation of serine-type peptidases at both asymptomatic and symptomatic stages  <cit> .
in infected navel orange, seventeen proteins involved in amino acid biosynthesis were downregulated.
fewer proteins were downregulated in infected volkameriana.
on the other hand, some upregulated proteins involved in amino acid degradation were identified only in infected navel orange.
taken together, these findings suggest that amino acid metabolism in navel orange is more sensitive to degradation during hlb infection than in volkameriana.
secondary metabolism
a general repression of key proteins involved in biosynthesis of secondary metabolites was observed in both cultivars in response to calas infection .
geranylgeranyl reductase and 1-deoxy-d-xylulose 5-phosphate reductoisomerase, involved in the non-mva pathway, were repressed.
other commonly hlb-downregulated proteins involved in the shikimate pathway included 2-dehydro-3-deoxyphosphoheptonate aldolase, 3-dehydroquinate synthase, 3-phosphoshikimate-1-carboxylvinyltransferase and mevalonate diphosphate decarboxylase.
two key proteins involved in phenylpropanoids, phenylalanine ammonia lyase and aryl-alcohol dehydrogenase, were more abundant in both species in response to hlb, while cinnamyl alcohol dehydrogenase  <dig> was repressed.
two proteins involved in alkaloid biosynthesis, tropinone reductase and strictosidine synthase-like  <dig>  were upregulated in infected volkameriana.
findings related to citrus activated defense responses against calas infection are shown in fig.
4).
three proteins involved in auxin signal transduction were activated in infected volkameriana: auxin resistant  <dig> and two aldo/keto reductases.
one aldo/keto reductase was hlb-regulated in both genotypes.
three proteins involved in jasmonic and salicylic acid responses were induced in navel orange but not in volkameriana.
lipoxygenase  <dig> was upregulated in both genotypes.. some key proteins involved in cell wall modifications were commonly regulated by both genotypes: udp-glucose-6-dehydrogenase, uxs <dig>  uxs <dig>  and rhm <dig>  proteolytic-related proteins were altered in both genotypes.
taken together these findings do not suggest any clear association between the two genotypes and proteomic changes in hormonal crosstalk, cell wall and proteolytic pathways.fig.
<dig> biotic stress–related proteins altered in response to hlb disease in both citrus genotypes.
proteins that were differentially expressed between control and infected trees are indicated by colors, based on their pattern of expression in the two genotypes .
each colored square represents expression change in a protein associated with regulatory and enzymatic functions
some proteins synthesizing volatiles via the phenylpropanoid and carotenoid pathways were affected in the present study .
the marked differences between the two species suggests that to be effective, any innovative hlb-detection system based on induced volatiles must be cultivar-specific.fig.
<dig> global view of proteomic changes in citrus leaves  in response to calas infections.
proteins, pathways, and cell functions that were differentially expressed are indicated by colors, based on their pattern of expression in the two genotypes.
each colored square represents expression change in a protein associated with regulatory and enzymatic functions.
more than one rectangle grouped together indicates different members of the same protein family were found to be differentially regulated
plant phenols not only counteract reactive oxygen species and pathogen-secreted toxins, but also play roles in transport and signal transduction pathways.
polyphenol chemistry is critical to adapting plants to environmental stresses, including pathogens  <cit> .
phenylalanine ammonia-lyase  was downregulated by hlb infection in both species .
this enzyme is a key regulatory point for the entire phenylpropanoid pathway.
enhancement of its transcript abundance is linked to phytopathogen attacks  <cit> .
another important protein encoding isoflavone reductase and involved in antioxidant reactions was repressed in both the infected species, consistent with previous proteomic analysis  <cit> .
in general, rna-seq analysis showed that phenylpropanoid pathways were transcriptionally affected by hlb in both young and mature infected leaves in field grown mature trees  <cit> .
these previous findings are not completely consistent with the present proteomic data.
gene set enrichment analysis showed a general downregulation of secondary metabolism in both genotypes, although some key proteins were upregulated in response to the disease.
these contrasting findings may reflect differences in developmental and physiological stages of the plants analyzed in the two studies or differences in environmental and agronomic conditions.
although phenylpropanoids may be activated at early stages of infection, their repression when symptoms are severe is expected.
the large number and complexity of metabolites belonging to phenylpropanoid pathways makes the clarification of their many roles in the host-pathogen battle difficult.
hormonal crosstalk
plant innate responses are finely controlled by hormonal crosstalk, particularly between ja and sa signaling.
the induction of allene oxide cyclase in infected navel orange  is intriguing because jasmonic acid signal transduction is a key pathway activated in response to necrotroph and herbivore attacks.
ja-responsive proteins, gibberellin signaling , and auxin signaling  were upregulated in mature infected leaves of valencia orange  <cit> .
because calas is a biotroph, the activation of ja may be deleterious for the host as previously suggested  <cit> .
brassinosteroids affect disease resistance in plants  <cit>  probably due to the induction of bak <dig>  which interacts with pamp receptors such as bacterial flagellin to activate immune responses  <cit> .
the st <dig> gene was repressed in hlb-infected leaves  <cit> .
however, the connections of brassinosteroids with the sa-et-ja crosstalk and plant immunity remain elusive.
on the other hand, auxin-related proteins were upregulated by hlb in both citrus genotypes.
because of the antagonist role of auxins toward sa response  <cit> , we may speculate that these effects are deleterious to the infected host.
proteomic studies have revealed post-transcriptional regulation of genes involved in key pathways which may be responsible for variations in phenotypic responses to hlb.
the regulation of carbohydrate metabolism  is clearly altered at symptomatic stage at both transcript and protein levels.
differences in signaling mechanisms and hormone-mediated defense responses may also contribute to the range of tolerance to hlb within the citrus germplasm.
the most compelling finding was that proteins involved in redox pathways and defense against xenobiotics  were more abundant in volkameriana than in navel orange, and this may be linked with the former’s greater tolerance.
four gsts were significantly upregulated in infected volkameriana and not in infected navel orange: gst <dig>  gst <dig>  ghst <dig>  ln2- <dig>  these effects on the proteome may explain the greater susceptibility of navel orange compared to volkameriana.
while these findings regarding citrus responses are valuable in the ongoing efforts to combat this deadly disease, further studies are still needed to validate these findings and deliver effective targets to develop new therapeutic strategies.
signaling and defense response pathways
some proteins encoding receptor kinases of the lrr type were hlb-regulated: one was upregulated in volkameriana, one was repressed in navel orange and two different proteins of lrr-viii were upregulated .
leucine-rich repeat receptor kinases are the largest category of receptor kinases and mediate signaling of defense responses in plants.
other receptor kinases  were upregulated in both citrus genotypes.
the proteins belonging to domain of unknown function  <dig> were upregulated in response to hlb: one in volkameriana, one in navel orange and one commonly regulated.
duf <dig> is one of the largest classes of receptor-like kinases .
these proteins play important roles in regulating pathogen defense and programmed cell death  <cit> .
based on the proteomic results in this study, we speculate that diverse signaling mechanisms occur depending on the citrus genotype.
it is possible that variability in susceptibility of citrus may result from pathogen perception due to activation of different receptors which in turn activate defense responses which vary in speed and intensity.
much remains to be learned regarding which receptor family is involved in susceptible or resistant responses to hlb.
several proteins involved in calcium regulation were less abundant in infected leaves of navel orange: calmodulin-binding and calcium-transporting atpase.
this is consistent with the significant drop in calcium concentration observed in symptomatic leaves  <cit> .
three 14-3- <dig> proteins were also hlb-downregulated in navel orange and not in volkameriana.
these are a large family of proteins present in all eukaryotic organisms that aid signaling by binding other proteins such as kinases, phosphatases, and receptors   <cit> .
volkameriana showed significant stimulation of respiratory burst and consequent redox state, a prerequisite of the upregulation of pathogenesis-related proteins.
enzymes involved in the control of reactive oxygen species were generally enhanced in response to calas infection in both genotypes although volkameriana showed a higher activation of glutathione-s-transferases .
the upregulation of these important detoxification proteins may be linked with the increased tolerance of volkameriana in comparison to navel orange.
an upregulation of enzymes involved in the biosynthesis of peroxiredoxins, cu/zn superoxide dismutase and 2fe-2s ferredoxin-like protein, occurred at both asymptomatic and symptomatic stages  <cit> .
glutathione s-transferase family proteins include several isozymes that help detoxify xenobiotic compounds  <cit> .
plant gsts add glutathione to electrophilic xenobiotic molecules pushing them into the cell vacuole.
regulation of these proteins by environmental stress stimuli suggests a role in protection against any harmful event  <cit> .
gsts have been linked with hormone homeostasis and their high affinity for auxins suggests their upregulation is a general signal of responses to stress  <cit> .
stress-inducible gsts conjugate deleterious metabolites caused by oxidative damage.
inducible gsts may play the important role of detoxifying exogenous molecules such as phytotoxins produced by pathogen attacks.
higher levels of gsts in volkameriana strengthens the hypothesis that they protect against dangerous molecules generated by calas attack.
indeed, differential activation of gsts may explain some of the variability of citrus responses to hlb.
some peroxidases were also upregulated in both genotypes.
these enzymes detoxify excess h2o <dig>  <cit>  and they are grouped based on their subcellular localization  <cit> .
as previously suggested  <cit> , greater citrus susceptibility to hlb may be linked to a failure to rapidly induce antioxidant components to alleviate the devastating effects of ros produced by calas.
it has been suggested  <cit>  that this category of proteins may be considered candidate markers in field-grown for the detection of the devastating disease “esca” in grapevine.
in the same way high expression of gsts are potential candidate markers for genotypes with useful tolerance to hlb.
further investigations will need to be conducted for a large number of genotypes.
pathogenesis-related  proteins are plant defensive proteins against biotic attacks  <cit> .
more pr-proteins were induced in volkameriana than in navel orange, consistent with the differing tolerance.
resistance  genes specifically activate a resistance reaction to a particular pathogen.
nbs-lrr proteins are the most numerous r-gene class.
nbs-lrr genes are finely controlled by regulatory mechanisms that allow their expression only when a biotic attack occurs, and limiting their metabolic cost when they are not required  <cit> .
it is possible that the two nbs-lrr proteins upregulated in volkameriana may contribute to enhanced tolerance of hlb disease.
heat shock proteins  are molecular chaperones with important functions in non-covalent protein folding or unfolding, assembly, and modifications.
genes encoding hsp <dig>  hsp <dig> and other small heat shock proteins were expressed at lower levels in hlb disease in both fruit and leaf tissues  <cit> .
down-regulation of hsp <dig>  chaperon-60kd and chaperonin-60alpha was also seen in infected grapefruit  <cit> .
in the present study, hsp <dig>  hsp <dig> and hsp <dig> were induced in volkameriana while several hsp proteins were inhibited in navel orange.
taken together, we conclude that the observed upregulation of some hsps in volkameriana may contribute to increased tolerance to hlb disease.
overall metabolism
the repression of key proteins involved in photosynthetic light reactions was linked with the upregulation of starch-related pathways.
infected volkameriana exhibited up-regulated nitrilases, oxidases, glutathione-s-transferases and other proteins involved in redox state.
infected volkameriana also exhibited enhanced production of expansins and xyloglucan endotransglycosylases.
proteins commonly altered by hlb in both genotypes strengthen the data at proteomic level.
the wd <dig> repeat-like protein was upregulated along with some enzymes involved in protein targeting, degradation, and glycosylation such as cysteine peptidase  <dig>  proteinase a <dig>  and cysteine protease.
additional file 5: figure s <dig> presents a complete list of commonly regulated proteins and their abundance in the four examined citrus sample categories.
CONCLUSIONS
forty-six of the  <dig> proteins were encoded by genes that were found to be transcriptionally regulated by hlb by previous studies in field conditions.
any comparison between transcriptomic and proteomic studies should take into consideration that experiments were performed with different plant material grown under differing environmental, developmental and physiological conditions.
we have applied an integrated approach using principal component analysis , gene set enrichment analysis and functional data mining to identify specific key proteomic changes in response to huanglongbing disease in these two citrus genotypes.
the analysis of post-transcriptional mechanisms is an essential step to link molecular regulatory networks to observed phenotypic changes.
interestingly, the clearest differences between the two genotypes were observed for proteins involved in redox state and detoxification pathways such as glutathione-s-tranferases.
hlb disease strongly affected genes and metabolites of the terpenoid, carotenoid, and jasmonic acid pathways  <cit> .

 huanglongbing disease currently threatens areas where citrus cultivation is important in the agricultural economy such as east asia, the middle east, and the americas. huanglongbing disease is caused by three species of candidatus liberibacter asiaticus, americanus and africanus. the pathogen is transmitted by two species of psyllids: diaphorina citri and trioza erytreae. recently, trioza erytreae was found for the first time in europe. typical symptoms of huanglongbing disease in leaves include shoot yellowing and blotchy, mottled leaves. although most of the fruits are still of commercial quality, fruits from severely affected branches are unmarketable: small, lopsided, green, and acidic, with many aborted seeds. leaves accumulate starch, phloem is damaged and cell wall lamellae swell during calas infection. candidatus liberibacter spp. belong to the alpha subdivision of the proteobacteria based on ribosomal region sequence data. the bacterium has not yet been definitively cultured despite attempts to do so. koch’s postulates have not been fulfilled for this disease, so possible interactions with other microrganisms cannot be ruled out. the pathogen lives in the insect and in the phloem of citrus trees. once acquired, it typically persists for the rest of the life of the host. insecticides can decrease psyllid populations, but since the pathogen remains in the vector, disease spread can occur with the presence of just a few infected psyllids in the orchard. all genotypes within the genus citrus are susceptible to hlb to varying degrees although species of other close-related genera showed some sort of resistance. there is variability in disease severity and symptoms among citrus genotypes. murraya paniculata, an ornamental citrus closely-related plant, showed fewer symptoms of the disease. a study examining the responses of genotypes to hlb disease grouped them based on phenotypic analysis of induced symptoms. another recent study has evaluated citrus accessions and accessions belonging to other closely related genera. resistance was reported in accessions not belonging to citrus genera. another work have screened citrus germplasm susceptibility to hlb analyzing sixteen citrus genotypes. results showed that citrus macrophylla and c. medica were the most susceptible while complex genetic hybrids ‘us 1-4-59’ and ‘fallglo’ were the least susceptible. a metabolomic investigation was conducted comparing five different tolerant hybrids and a highly susceptible cultivar to identify potential metabolites linked with diverse response. the causes of the disease have been studied using different “omic” approaches to identify which genes, proteins and metabolites may be targeted by innovative diagnostic and therapeutic methods. the genome of the pathogen was sequenced using a metagenomic approach, both from infected plants and the insect vector. no toxins or other secreted proteins have been linked to the disease and the mechanisms of its pathology are still unclear. large scale microarray analysis revealed significant modulation of genes involved in transport, cell defense and carbohydrate metabolism. photosynthesis is diminished in both young and mature leaves, but it is upregulated in infected fruits. starch accumulation was linked to the upregulation of genes involved in glucose import into the chloroplast and starch biosynthesis. a modulated jasmonic -salicilic acid crosstalk of innate responses may lead to a misdirected defence response. an integrated approach of 2-de and mass spectrometry showed that changes in levels of several proteins involved in photosynthesis and protein synthesis were linked to reduced concentrations of ca, mg, fe, zn, mn and cu in infected grapefruit leaves. proteins upregulated in infected samples were involved in redox stage and cell defense such as cu/zn superoxide dismutase, peroxidases, chitinases and lectin-related proteins. ‘madam vinous’ sweet orange plants infected by calas showed increased miraculin-like proteins, chitinase, cu/zn superoxide dismutase and lipoxygenase. some key metabolites modulated by hlb include proline, β-elemene, -trans-caryophyllene, and α-humulene. increased accumulation of some amino acids and organic acids was linked to greater susceptibility of ‘madam vinous’ sweet orange compared to carrizo citrange. however, it is worthy to notice that these trees were highly infected with many secondary effects so it will be necessary to confirm these results with newly infected trees. an increased amount of most amino acids, involved in plant defense to pathogens was observed in tolerant varieties in such as phenylalanine, tyrosine, tryptophan, lysine, and asparagine. this study examines proteomic changes in fully photosynthesizing leaves to determine how disease mechanisms and susceptibility vary between two citrus genotypes, using an integrated approach of principal component analysis, gene set and pathway enrichment analysis. the purpose was to characterize key proteins and post-transcriptionally modulated pathways in different citrus genotypes at a late symptomatic stage of hlb.  navel orange ) is an hlb-sensitive cultivar while volkameriana is moderately tolerant. different techniques have been used to identify key proteins differentially regulated by hlb in citrus leaves. over proteins were analyzed using isobaric tags for relative and absolute quantitation for both genotypes. in navel orange, proteins were differentially regulated between infected and healthy tissue. in volkameriana, differentially regulated proteins were found between infected and healthy tissue. pca analysis pca was used to visualize differences between the four analyzed genotypes x disease status, subdividing the entire proteomic profile into three important subcategories: biotic stress responses, overall cell metabolism and transcriptional regulation pathways. all identified proteins belonging to these important functional categories were used for the pca plots. the four categories of samples were clearly separated in all three pca plots, implying significant protein changes in all three gene categories related to species and health status. for biotic stress-related proteins, pc and pc accounted for and 24 % of the data variability, respectively. important proteins associated with each functional category that contributed to separation between sample types are listed in fig. some key proteins involved in redox state significantly contributed in the separation of infected v from the rest of the other categories. in the general metabolism pca plot, the pc and pc accounted for and 24 % of the data variability, respectively. key proteins involved in primary metabolism associated with the separation of healthy v from the infected v include malate dehydrogenase and pyruvate dehydrogenase, sucrose synthase and agpase. the third pca plot was generated based on the expression of proteins involved in transcriptional regulation, signaling, hormone and redox state. pc and pc accounted for and 27 % of data variability, respectively. interestingly, the regulation of few important proteins seems to specifically characterize the infected vs. healthy state of navel orange. these include map kinase udp-glucosyl transferase, and aspartyl protease. no changes were observed for mapk proteins that appear in the three pca plots were highly regulated in the comparison between infected and control in both citrus genotypes. indeed, they may be considered as putative candidate biomarkers of a clear symptomatic status in citrus at proteomic level. their hlb-regulated pattern of expression greatly contributed in distinguishing the four different leaf sample types. further analysis will be needed to validate these data and confirm the role in the pathogenesis of hlb disease.fig. principal component analysis of differentially regulated proteins of four types of leaf tissues. proteins that contribute highly to the separation of the the four samples are numbered and listed next to each graph gene set and pathway enrichment analysis the pageman web tool highlighted which gene categories were up- or down-regulated in each pairwise comparison. both infected genotypes exhibited repressed amino acid biosynthesis and protein synthesis. in navel orange, photosynthesis, isoflavone pathways, tetrapyrrole synthesis, and galactose metabolism were significantly inhibited by hlb. in volkameriana, s-assimilation, isoprenoids, rna binding and amino acid activation were specifically diminished. in both genotypes, hlb enhanced starch-related pathways, biotic stress-related proteins, beta glucan hydrolases, and protein degradation pathways. some distinct differences were observed between genotypes in hlb response. in volkameriana, cell wall modifications, galactose metabolism, and heat shock proteins were upregulated. in navel orange, amino acid degradation, lipid metabolism, jasmonates, and pr-proteins were upregulated.fig. gene set enrichment analysis using pageman web-tool. upregulated and downregulated pathways at proteomic level in infected volkameriana and infected navel orange in comparison to respective healthy controls a pathway enrichment analysis was performed using the david bioinformatic resource to determine which metabolic pathways were commonly or specifically affected by hlb disease in both genotypes. some metabolic pathways were altered by hlb in both species. amino acid metabolism was significantly downregulated. other key inhibited pathways include biosynthesis of plant hormones, terpenoids, and phenylpropanoids. on the other hand, tyrosine metabolism was upregulated in both genotypes. in navel orange, fatty acid biosynthesis and nitrogen metabolism were diminished while alpha-linolenic acid metabolism was enhanced. in volkameriana, alkaloids and pyruvate metabolism-related proteins were repressed while galactose metabolism and fatty acid metabolism were upregulated.table pathway enrichment analysis using david bioinformatics resources. pathways that were upregulated and downregulated for the healthy/infected comparison for each genotype are indicated with the corresponding p-value primary metabolism the integration of the two pairwise comparisons into a unified mapman visualization allowed us to identify proteins that were commonly or specifically regulated in response to hlb in the two species. some proteins involved in cell wall modifications were upregulated by hlb only in volkameriana: expansin a expansin a expansin-like b xyloglucan endotransglycosylase and xyloglucanxyloglucosyltransferase. in navel orange, proteins involved in fatty acid biosynthesis and elongation were generally repressed while several key proteins involved in amino acid degradation were upregulated: arginase, pyrroline-5-carboxylase reductase, lactoylglutathione lyase, and 3-hydroxylmethylglutaryl-coa lyase. the increase of protein degradation indicates that senescence processes may be more highly activated in hlb-diseased navel orange than in volkameriana. sucrose metabolism was only slightly affected by hlb at the protein level; only sucrose synthase was repressed. starch metabolism was more altered. the first enzyme of starch biosynthesis, adp-glucose pyrophosphorylase, was inhibited by hlb in both genotypes. in volkameriana, starch synthase was upregulated while -alpha-glucan starch branching enzyme was slightly downregulated. among starch degradation enzymes, glucan phosphorylase and heterogycan glucosidase were upregulated in both genotypes. alpha-amylase was upregulated in volkameriana while beta-amylase was enhanced in navel orange. taken together these findings showed that starch metabolism was highly affected in both genotypes at the protein level. no clear association between differing susceptibility to hlb is evident from starch pathway regulation alone.fig. proteins involved in primary metabolic responses to hlb. a) sucrose and starch metabolism; b) raffinose metabolism; c) tca cycle; d) oxidative pentose phosphate. proteins that were differentially expressed between control and infected trees are indicated by colors, based on their pattern of expression in the two genotypes. each colored square represents expression change in a protein catalyzing a step in the pathway. more than one square grouped together indicates different members of the same protein family found to be differentially regulated a significant downregulation in hlb-infected samples was observed for proteins involved in photosynthetic reactions. the altered transcription of sugar and starch metabolism genes caused by hlb mostly agreed with the corresponding protein changes presented in the present work. the observed changes in starch-related pathways were consistent with the transcriptomic analysis. increased starch degradation was observed probably due to the increased starch concentrations in infected leaf tissues. adp glucose-pyrophospholyase was repressed in both infected species. this enzyme is rate-limiting for starch biosynthesis, catalyzing the conversion of glucose-1-phosphate to adp-glucose that is polymerized into amylopectin or amylose. starch accumulation in infected leaves is a typical symptom of hlb. however, the greatest occurrence of this process may occur at an early, asymptomatic stage. indeed we may speculate that the accumulation of starch may be a secondary effect of the disease instead of being the cause of symptoms. when symptoms are already evident and yellowing is present, starch biosynthesis may slow down and starch degradation is expected to be activated as a response of the plant to limit damage to cell structures. however, some differences in primary metabolism were observed between transcriptomic and proteomic approaches. changes in invertase gene expression were observed, but not its protein levels. sucrose, the substrate of invertase, is produced by photosynthesis in leaves and then transported to sink tissues through the phloem. involvement of sucrose in signaling of innate responses has been described recently. invertase plays an important role in plant stress responses, possibly serving as an extracellular signal of pathogen attacks. it is possible that expression of this gene in response to hlb depends on genotype, plant physiological conditions and age, tissue developmental conditions, type of infection and/or environmental conditions. taken together all these primary metabolism results concur with previous findings showing significant modification of transcript abundance in minor carbohydrate metabolism but do not suggest a clear link to the well-known difference in tolerance between the two genotypes. gpt has been linked with hlb disease: this gene is responsible for glucose import into the chloroplast and consequently for starch accumulation. the protein was not found among those extracted and characterized by itraq, therefore no conclusions can be made about changes in protein levels due to hlb. however, the protocol used to analyze the citrus proteome favored detection of soluble cell proteins over membrane proteins such as gpt expression of key proteins involved in the tca cycle and pep metabolism were affected by hlb. atp-citrate lyase subunit b was repressed in both species. alcohol dehydrogenase which converts aceltaldehyde to ethanol, was upregulated in both genotypes. in infected navel orange, isopropyl malate isomerase, which converts citrate to isocitrate and pyruvate decarboxylase involved in fermentation, was more abundant than in healthy tissue. in volkameriana, pyruvate decarboxylase was repressed while aconitate hydratase involved in tca cycle was upregulated. raffinose metabolism was drastically altered by the disease. in navel orange, expression of alpha-galactosidase udp-glucose-4-epimerase, glucose-6-phosphate isomerase, and alpha-galctosidase were enhanced by hlb. in volkameriana, raffinose synthase, phosphofructokinase and phosphoglucomutase were upregulated. sucrose synthase and phosphoglycerate mutase were donwregulated in both genotypes. a transketolase that converts the xylulose-5-p in sedoheptulose-7-p was repressed in infected navel orange. significant repression of aspartate biosynthesis and serine metabolism was observed in infected navel orange leaves. this agrees with a previously described downregulation of serine-type peptidases at both asymptomatic and symptomatic stages. in infected navel orange, seventeen proteins involved in amino acid biosynthesis were downregulated. fewer proteins were downregulated in infected volkameriana. on the other hand, some upregulated proteins involved in amino acid degradation were identified only in infected navel orange. taken together, these findings suggest that amino acid metabolism in navel orange is more sensitive to degradation during hlb infection than in volkameriana. secondary metabolism a general repression of key proteins involved in biosynthesis of secondary metabolites was observed in both cultivars in response to calas infection. geranylgeranyl reductase and 1-deoxy-d-xylulose 5-phosphate reductoisomerase, involved in the non-mva pathway, were repressed. other commonly hlb-downregulated proteins involved in the shikimate pathway included 2-dehydro-3-deoxyphosphoheptonate aldolase, 3-dehydroquinate synthase, 3-phosphoshikimate-1-carboxylvinyltransferase and mevalonate diphosphate decarboxylase. two key proteins involved in phenylpropanoids, phenylalanine ammonia lyase and aryl-alcohol dehydrogenase, were more abundant in both species in response to hlb, while cinnamyl alcohol dehydrogenase was repressed. two proteins involved in alkaloid biosynthesis, tropinone reductase and strictosidine synthase-like were upregulated in infected volkameriana. findings related to citrus activated defense responses against calas infection are shown in fig. 4). three proteins involved in auxin signal transduction were activated in infected volkameriana: auxin resistant and two aldo/keto reductases. one aldo/keto reductase was hlb-regulated in both genotypes. three proteins involved in jasmonic and salicylic acid responses were induced in navel orange but not in volkameriana. lipoxygenase was upregulated in both genotypes.. some key proteins involved in cell wall modifications were commonly regulated by both genotypes: udp-glucose-6-dehydrogenase, uxs uxs and rhm proteolytic-related proteins were altered in both genotypes. taken together these findings do not suggest any clear association between the two genotypes and proteomic changes in hormonal crosstalk, cell wall and proteolytic pathways.fig. biotic stress–related proteins altered in response to hlb disease in both citrus genotypes. proteins that were differentially expressed between control and infected trees are indicated by colors, based on their pattern of expression in the two genotypes. each colored square represents expression change in a protein associated with regulatory and enzymatic functions some proteins synthesizing volatiles via the phenylpropanoid and carotenoid pathways were affected in the present study. the marked differences between the two species suggests that to be effective, any innovative hlb-detection system based on induced volatiles must be cultivar-specific.fig. global view of proteomic changes in citrus leaves in response to calas infections. proteins, pathways, and cell functions that were differentially expressed are indicated by colors, based on their pattern of expression in the two genotypes. each colored square represents expression change in a protein associated with regulatory and enzymatic functions. more than one rectangle grouped together indicates different members of the same protein family were found to be differentially regulated plant phenols not only counteract reactive oxygen species and pathogen-secreted toxins, but also play roles in transport and signal transduction pathways. polyphenol chemistry is critical to adapting plants to environmental stresses, including pathogens. phenylalanine ammonia-lyase was downregulated by hlb infection in both species. this enzyme is a key regulatory point for the entire phenylpropanoid pathway. enhancement of its transcript abundance is linked to phytopathogen attacks. another important protein encoding isoflavone reductase and involved in antioxidant reactions was repressed in both the infected species, consistent with previous proteomic analysis. in general, rna-seq analysis showed that phenylpropanoid pathways were transcriptionally affected by hlb in both young and mature infected leaves in field grown mature trees. these previous findings are not completely consistent with the present proteomic data. gene set enrichment analysis showed a general downregulation of secondary metabolism in both genotypes, although some key proteins were upregulated in response to the disease. these contrasting findings may reflect differences in developmental and physiological stages of the plants analyzed in the two studies or differences in environmental and agronomic conditions. although phenylpropanoids may be activated at early stages of infection, their repression when symptoms are severe is expected. the large number and complexity of metabolites belonging to phenylpropanoid pathways makes the clarification of their many roles in the host-pathogen battle difficult. hormonal crosstalk plant innate responses are finely controlled by hormonal crosstalk, particularly between ja and sa signaling. the induction of allene oxide cyclase in infected navel orange is intriguing because jasmonic acid signal transduction is a key pathway activated in response to necrotroph and herbivore attacks. ja-responsive proteins, gibberellin signaling, and auxin signaling were upregulated in mature infected leaves of valencia orange. because calas is a biotroph, the activation of ja may be deleterious for the host as previously suggested. brassinosteroids affect disease resistance in plants probably due to the induction of bak which interacts with pamp receptors such as bacterial flagellin to activate immune responses. the st gene was repressed in hlb-infected leaves. however, the connections of brassinosteroids with the sa-et-ja crosstalk and plant immunity remain elusive. on the other hand, auxin-related proteins were upregulated by hlb in both citrus genotypes. because of the antagonist role of auxins toward sa response, we may speculate that these effects are deleterious to the infected host. proteomic studies have revealed post-transcriptional regulation of genes involved in key pathways which may be responsible for variations in phenotypic responses to hlb. the regulation of carbohydrate metabolism is clearly altered at symptomatic stage at both transcript and protein levels. differences in signaling mechanisms and hormone-mediated defense responses may also contribute to the range of tolerance to hlb within the citrus germplasm. the most compelling finding was that proteins involved in redox pathways and defense against xenobiotics were more abundant in volkameriana than in navel orange, and this may be linked with the former’s greater tolerance. four gsts were significantly upregulated in infected volkameriana and not in infected navel orange: gst gst ghst ln2- these effects on the proteome may explain the greater susceptibility of navel orange compared to volkameriana. while these findings regarding citrus responses are valuable in the ongoing efforts to combat this deadly disease, further studies are still needed to validate these findings and deliver effective targets to develop new therapeutic strategies. signaling and defense response pathways some proteins encoding receptor kinases of the lrr type were hlb-regulated: one was upregulated in volkameriana, one was repressed in navel orange and two different proteins of lrr-viii were upregulated. leucine-rich repeat receptor kinases are the largest category of receptor kinases and mediate signaling of defense responses in plants. other receptor kinases were upregulated in both citrus genotypes. the proteins belonging to domain of unknown function were upregulated in response to hlb: one in volkameriana, one in navel orange and one commonly regulated. duf is one of the largest classes of receptor-like kinases. these proteins play important roles in regulating pathogen defense and programmed cell death. based on the proteomic results in this study, we speculate that diverse signaling mechanisms occur depending on the citrus genotype. it is possible that variability in susceptibility of citrus may result from pathogen perception due to activation of different receptors which in turn activate defense responses which vary in speed and intensity. much remains to be learned regarding which receptor family is involved in susceptible or resistant responses to hlb. several proteins involved in calcium regulation were less abundant in infected leaves of navel orange: calmodulin-binding and calcium-transporting atpase. this is consistent with the significant drop in calcium concentration observed in symptomatic leaves. three 14-3- proteins were also hlb-downregulated in navel orange and not in volkameriana. these are a large family of proteins present in all eukaryotic organisms that aid signaling by binding other proteins such as kinases, phosphatases, and receptors. volkameriana showed significant stimulation of respiratory burst and consequent redox state, a prerequisite of the upregulation of pathogenesis-related proteins. enzymes involved in the control of reactive oxygen species were generally enhanced in response to calas infection in both genotypes although volkameriana showed a higher activation of glutathione-s-transferases. the upregulation of these important detoxification proteins may be linked with the increased tolerance of volkameriana in comparison to navel orange. an upregulation of enzymes involved in the biosynthesis of peroxiredoxins, cu/zn superoxide dismutase and 2fe-2s ferredoxin-like protein, occurred at both asymptomatic and symptomatic stages. glutathione s-transferase family proteins include several isozymes that help detoxify xenobiotic compounds. plant gsts add glutathione to electrophilic xenobiotic molecules pushing them into the cell vacuole. regulation of these proteins by environmental stress stimuli suggests a role in protection against any harmful event. gsts have been linked with hormone homeostasis and their high affinity for auxins suggests their upregulation is a general signal of responses to stress. stress-inducible gsts conjugate deleterious metabolites caused by oxidative damage. inducible gsts may play the important role of detoxifying exogenous molecules such as phytotoxins produced by pathogen attacks. higher levels of gsts in volkameriana strengthens the hypothesis that they protect against dangerous molecules generated by calas attack. indeed, differential activation of gsts may explain some of the variability of citrus responses to hlb. some peroxidases were also upregulated in both genotypes. these enzymes detoxify excess h2o and they are grouped based on their subcellular localization. as previously suggested, greater citrus susceptibility to hlb may be linked to a failure to rapidly induce antioxidant components to alleviate the devastating effects of ros produced by calas. it has been suggested that this category of proteins may be considered candidate markers in field-grown for the detection of the devastating disease “esca” in grapevine. in the same way high expression of gsts are potential candidate markers for genotypes with useful tolerance to hlb. further investigations will need to be conducted for a large number of genotypes. pathogenesis-related proteins are plant defensive proteins against biotic attacks. more pr-proteins were induced in volkameriana than in navel orange, consistent with the differing tolerance. resistance genes specifically activate a resistance reaction to a particular pathogen. nbs-lrr proteins are the most numerous r-gene class. nbs-lrr genes are finely controlled by regulatory mechanisms that allow their expression only when a biotic attack occurs, and limiting their metabolic cost when they are not required. it is possible that the two nbs-lrr proteins upregulated in volkameriana may contribute to enhanced tolerance of hlb disease. heat shock proteins are molecular chaperones with important functions in non-covalent protein folding or unfolding, assembly, and modifications. genes encoding hsp hsp and other small heat shock proteins were expressed at lower levels in hlb disease in both fruit and leaf tissues. down-regulation of hsp chaperon-60kd and chaperonin-60alpha was also seen in infected grapefruit. in the present study, hsp hsp and hsp were induced in volkameriana while several hsp proteins were inhibited in navel orange. taken together, we conclude that the observed upregulation of some hsps in volkameriana may contribute to increased tolerance to hlb disease. overall metabolism the repression of key proteins involved in photosynthetic light reactions was linked with the upregulation of starch-related pathways. infected volkameriana exhibited up-regulated nitrilases, oxidases, glutathione-s-transferases and other proteins involved in redox state. infected volkameriana also exhibited enhanced production of expansins and xyloglucan endotransglycosylases. proteins commonly altered by hlb in both genotypes strengthen the data at proteomic level. the wd repeat-like protein was upregulated along with some enzymes involved in protein targeting, degradation, and glycosylation such as cysteine peptidase proteinase a and cysteine protease. additional file 5: figure s presents a complete list of commonly regulated proteins and their abundance in the four examined citrus sample categories.  forty-six of the proteins were encoded by genes that were found to be transcriptionally regulated by hlb by previous studies in field conditions. any comparison between transcriptomic and proteomic studies should take into consideration that experiments were performed with different plant material grown under differing environmental, developmental and physiological conditions. we have applied an integrated approach using principal component analysis, gene set enrichment analysis and functional data mining to identify specific key proteomic changes in response to huanglongbing disease in these two citrus genotypes. the analysis of post-transcriptional mechanisms is an essential step to link molecular regulatory networks to observed phenotypic changes. interestingly, the clearest differences between the two genotypes were observed for proteins involved in redox state and detoxification pathways such as glutathione-s-tranferases. hlb disease strongly affected genes and metabolites of the terpenoid, carotenoid, and jasmonic acid pathways. 
in this work we have conducted a proteomic analysis of mature leaves of two citrus genotypes with well-known differing tolerances to hlb: navel orange and volkameriana. although there is not any resistance source in the citrus germplasm, a certain level of moderated tolerance is present. pathway enrichment analysis showed that amino acid degradation processes occurred to a larger degree in the navel orange.
 huanglongbing disease is still the greatest threat to citriculture worldwide. although there is not any resistance source in the citrus germplasm, a certain level of moderated tolerance is present. a large-scale analysis of proteomic responses of citrus may help: 1) clarifying physiological and molecular effects of disease progression, 2) validating previous data at transcriptomic level, and 3) identifying biomarkers for development of early diagnostics, short-term therapeutics and long-term genetic resistance.  in this work we have conducted a proteomic analysis of mature leaves of two citrus genotypes with well-known differing tolerances to hlb: navel orange and volkameriana. pathway enrichment analysis showed that amino acid degradation processes occurred to a larger degree in the navel orange. no clear differences between the two genotypes were observed for primary metabolic pathways. the most important finding was that four glutathione-s-transferases were upregulated in volkameriana and not in navel orange. these proteins are involved in radical ion detoxification.  upregulation of proteins involved in radical ion detoxification should be considered as an important mechanism of increased tolerance to hlb. electronic supplementary material the online version of this article contains supplementary material, which is available to authorized users. keywords citrushuanglongbingcandidatus liberibacteritraqproteomeproteomichttp://dx.doi.org/http://dx.doi.org/ /100007565citrus research and development foundationissue-copyright-statement© the author 2016 