Skip to main content

Transcriptome analysis of sweet orange trees infected with ‘Candidatus Liberibacter asiaticus’ and two strains of Citrus Tristeza Virus



Huanglongbing (HLB) and tristeza, are diseases of citrus caused by a member of the α-proteobacteria, ‘Candidatus Liberibacter asiaticus’ (CaLas), and Citrus tristeza virus (CTV) respectively. HLB is a devastating disease, but CTV strains vary from very severe to very mild. Both CaLas and CTV are phloem-restricted. The CaLas-B232 strain and CTV-B6 cause a wide range of severe and similar symptoms. The mild strain CTV-B2 doesn’t induce significant symptoms or damage to plants.


Transcriptome profiles obtained through RNA-seq revealed 611, 404 and 285 differentially expressed transcripts (DETs) after infection with CaLas-B232, CTV-B6 and CTV-B2. These DETs were components of a wide range of pathways involved in circadian rhythm, cell wall modification and cell organization, as well as transcription factors, transport, hormone response and secondary metabolism, signaling and stress response. The number of transcripts that responded to both CTV-B6 and CaLas-B232 was much larger than the number of transcripts that responded to both strains of CTV or to both CTV-B2 and CaLas-B232. A total of 38 genes were assayed by RT-qPCR and the correlation coefficients between Gfold and RT-qPCR were 0.82, 0.69, 0.81 for sweet orange plants infected with CTV-B2, CTV-B6 and CaLas-B232, respectively.


The number and composition of DETs reflected the complexity of symptoms caused by the pathogens in established infections, although the leaf tissues sampled were asymptomatic. There were greater similarities between the sweet orange in response to CTV-B6 and CaLas-B232 than between the two CTV strains, reflecting the similar physiological changes caused by both CTV-B6 and CaLas-B232. The circadian rhythm system of plants was perturbed by all three pathogens, especially by CTV-B6, and the ion balance was also disrupted by all three pathogens, especially by CaLas-B232. Defense responses related to cell wall modification, transcriptional regulation, hormones, secondary metabolites, kinases and stress were activated by all three pathogens but with different patterns. The transcriptome profiles of Citrus sinensis identified host genes whose expression is affected by the presence of a pathogen in the phloem without producing symptoms (CTV-B2), and host genes whose expression leads to induction of symptoms in the plant (CTV-B6, CaLas-B232).


Citrus is grown in more than one hundred countries and is one of the most important commercial fruit crops in the world. Citrus is affected by different diseases and different geographic areas have distinct disease problems. Tristeza and huanglongbing (HLB) are two of the most destructive and widely distributed diseases [1]. Tristeza is caused by Citrus tristeza virus (CTV), a member of the Closteroviridae that is aphid-transmitted and phloem-restricted, with a single-stranded and positive sense RNA genome. CTV can be classified into mild and severe types based on the severity of symptoms produced on citrus indicators. Mild CTV strains cause only mild or no symptoms on indicators and usually bring about no economic loss. Severe CTV strains induce different disease syndromes, including decline and death of sweet orange on sour orange rootstocks, stem pitting (SP) of sweet orange and grapefruit scions, seedling yellows (SY) in sour orange, and vein clearing in ‘Mexican’ lime seedlings [2, 3]. HLB is associated with a member of the α-subdivision of the proteobacteria: ‘Candidatus Liberibacter asiaticus’ (CaLas), which grows only in the phloem of infected plants and the hemocoel and salivary glands of infected insect vectors. Because CaLas has not been grown in pure culture and Koch’s postulates have not been completed, it has ‘Candidatus’ status in the accepted nomenclature [4, 5]. HLB affected trees may produce new shoots that are uniformly yellow but mature to present blotchy mottle symptoms, usually with symptoms of zinc deficiency. Fruit are misshaped with reduced quality and quantity, abnormal coloration and aborted seeds. HLB is considered to be one of the most serious plant diseases in the world and there are no resistant varieties. Growers attempt to maintain yield in the presence of the disease with nutritional supplementation [6], and rigorous control of the citrus psyllid vector even after trees in a grove have already been infected by CaLas [7]. Benefits in terms of yield of these treatments have been reported [6, 7].

An understanding of host responses to pathogen infection is essential to clarify the mechanisms of plant-microbe interactions and to develop novel strategies for therapy and detection. Pathogens may activate plant defense responses, such as fortification of cell walls, the production of reactive oxygen species, phytoalexins and signaling molecules and the synthesis of pathogenesis-related proteins [8]. Transcriptome studies of the host can provide valuable information on such defense responses and may lead to new plant protection strategies. Transcriptome analyses of citrus plants in response to CTV and CaLas have been reported in separate studies [912]. Many host genes and related metabolic pathways and processes were affected by CaLas [10] and severe strain CTV-T305 compared to mild strain CTV-T385 [12].

Both CTV and CaLas are phloem-restricted pathogens and can occur separately or in combination on an individual tree. Mild CTV strains don’t cause damage to citrus hosts, while severe CTV strains and CaLas have overlapping symptom patterns. Therefore, we present a comparison of transcriptome profiles of citrus plants infected with a mild CTV strain, severe CTV strain and CaLas. We hypothesized that many host genes would be regulated in a similar manner by a severe strain of CTV and CaLas, even in asymptomatic host tissues.


Thirty-eight to 44 million clean reads were obtained from infected and control plants and mapped to the Citrus sinensis ‘Valencia’ reference genome [13], with approximately 73 % for CTV-B2, 76 % for the healthy, CTV-B6 and CaLa-B232 libraries (Additional file 1: Figure S1). Dramatic differences in the transcriptome profiles of sweet orange trees were observed in response to the three pathogens. Many more transcripts identified in the citrus gene database were up-regulated than down-regulated in response to the three pathogens (Additional file 1: Figure S2). More transcripts were differentially expressed in response to CaLas-B232, followed by CTV-B6 and CTV-B2 (Additional file 1: Figure S2). Among these differentially expressed transcripts (DETs), 50 transcripts were up-regulated, while only two transcripts were down-regulated in response to all three pathogens (Fig. 1). The number of transcripts regulated in response to CTV-B6 and CaLas-B232 (105 up and two down) was approximately four times the number of transcripts regulated in response to the two strains of CTV (27 up and 0 down) or CTV-B2 and CaLas-B232 (22 up and 3 down) (Fig. 1). The number of DETs responding to CaLas-B232 only (412 up and 15 down) was much larger than the number responding only to CTV-B2 (143 up and 37 down) or CTV-B6 (146 up and 72 down) (Fig. 1). A similar pattern was also found when the DETs were identified using the Arabidopsis Information Resource (TAIR) (Fig. 2). DETs with unknown functions were also identified (Additional file 1: Table S1).

Fig. 1

Transcripts in response to infection of C. sinensis with CTV-B2, CTV-B6 and CaLas-B232. Transcripts were mapped to the C. sinensis reference genome. a Number of up-regulated transcripts (b) Number of down-regulated transcripts

Fig. 2

Differentially expressed genes (DEGs) in C. sinensis in response to infection by CTV-B2, CTV-B6 and CaLas-B232. DEGs were identified as Arabidopsis thaliana orthologs

Plant Gene Sequence Enrichment Analysis (PlantGSEA) was performed with excellent p-values and low false discovery rates to reveal perturbed pathways. The functions of DETs observed in response to all three pathogens were mainly associated with the maintenance of circadian rhythms, such as flowering and photoperiodism, response to stimulus and response to other organism and signaling (Additional file 1: Table S2). DETs observed in response to CTV-B2 and CTV-B6 were disproportionally involved in hormone related responses, signaling pathways, and water stress (Additional file 1: Table S3). DETs observed in response to CTV-B2 and CaLas-B232 disproportionally encoded binding activities, including to iron, tetrapyrrole and heme molecules, which could affect chloroplast function, in addition to genes that respond to biotic stimulus or stress (Additional file 1: Table S4). A much larger set of DETs observed in response to both CTV-B6 and CaLas-B232 disproportionally encoded biological process (BP) functions in signaling, defense, immune and stress responses and other BPs (Additional file 1: Table S5). DETs observed to respond only to CTV-B2 were mainly associated with molecular functions (MF) and cellular components (CC) (Fig. 3) and disproportionately functioned in active transport of molecules across membranes and cell junctions (Additional file 1: Table S6). In contrast, the DETs observed to respond only to CTV-B6 or CaLas-B232 were primarily involved in BPs (Fig. 3). The long lists of DETs observed in response to CTV-B6 and CaLas-B232 are the reflection of the BPs disrupted by the severe disease symptoms caused by these pathogens, although symptoms had not yet developed in the tissues at the time they were sampled. Genes whose products functioned in transcription and inositol metabolism, as well as membrane bound intracellular organelles were up-regulated but signaling pathways were down regulated in response to CTV-B6 (Additional file 1: Table S7). The effects on inositol metabolism and membrane bound organelles could represent either the cause or an effect of chloroplast dysfunction seen in citrus infected by this pathogen, and negative regulation of signaling pathways likely impedes the defense response of the host. The long list of DETs that were observed in response only to CaLas-B232 included a large group associated with response to fungal pathogens, programmed cell death, hypersensitive responses and other host defense responses, and the production of jasmonic acid (JA) catalytic activity (Additional file 1: Table S8). The substantial number of genes annotated as ‘response to fungus’ may be correlated with the root decline and rot caused by CaLas and a potential interaction with Phytophthora spp. [14, 15]. Similar defense responses could also be triggered by effectors secreted by CaLas and delivered by fungus. DETs were also assigned to several specific function groups by Mapman (Fig. 4), including cell wall, hormone metabolism, secondary metabolism, transcription factors and signaling etc. Many pathways related to biotic stress were affected by CTV-B2, CTV-B6 and CaLas-B232 (Fig. 5). For more complete information about DETs, please refer to the supplemental information (Additional file 1: Table S9).

Fig. 3

Summary of Gene Ontology categories for differentially expressed genes (DEGs). BP, biological process; MF, molecular function; CC, cellular component. The DEGs were observed in response to infection with one or more of the pathogens CTV-B2, CTV-B6 or CaLas-B232 as shown in the figure. Y-axis, number of DEGs involved in BP, MF and CC

Fig. 4

Functional categories of DEGs in response to infection by CTV-B2, CTVB6 and CaLas-B232 through Mapman. X-axis, number of differentially expressed genes involved functional categories. Y-axis, functional categories of DEGs

Fig. 5

Summary of differentially expressed genes related to biotic stress in response to infection by CTV-B2, CTV-B6 and CaLas-B232

RT-qPCR validation of RNA-seq results

Thirty-eight genes were tested by replicated RT-qPCR. Expression levels ranged from −7 to 8 (log2FC) (Additional file 1: Table S10). Spearman’s rho values, 0.821 for CTV-B2, 0.687 for CTV-B6 and 0.813 for CaLas-B232, indicated excellent correlations between Gfold analysis of the RNA Seq libraries and RT-qPCR (Fig. 6). All primers showed amplification efficiencies between 90–110 % (Additional file 1: Table S11) and melting curves showed single peaks for the evaluated genes (Additional file 1: Figure S3).

Fig. 6

Correlation of differentially expressed genes was estimated by RT-qPCR and Gfold (log2FC). FC: Fold Change

Genes involved in cell wall and cell organization

Many transcripts related to cell wall modification and cell organization were differentially regulated, especially in response to CaLas-B232 (Additional file 1: Figure S4A). Transcripts encoding xyloglucan endotransglycosylase 6 (XTR6) were up-regulated by infection with all three pathogens (Table 1). A gene annotated as encoding a plant invertase/pectin methylesterase inhibitor (PMEI), was strongly up-regulated in response to CaLas-B232 (Table 1). Inhibition of pectin methyl esterase would strengthen the cell wall and is a common defense response [16, 17]. Another set of genes encoding expansin were up-regulated in response to CaLas-B232, in accordance with another study [10]. Expansin functions to loosen the structure of cell walls. Taken together this set of genes could contribute to the thickening and stiffening of leaves that occurs as a primary symptom of HLB and distortions of phloem cell walls as CaLas multiplies [18]. Transcripts encoding cellulose synthesis, including cellulose synthase-like B3 (CSLB3) and A2 (CSLA2), were sharply up-regulated by infection with CTV-B2 (Table 1). A transcript encoding COBRA-like protein 7 precursor (COBL7), likely involved in the deposition of cellulose, were also remarkably up-regulated by infection with CTV-B2 (Table 1), all consistent with a stronger cell wall in response to CTV-B2. A pectin lyase-like superfamily protein and beta-D-xylosidase 4 (XYL4) protein related to cell wall degradation were down-regulated by infection with CTV-B6 (Table 1). Transcripts encoding phloem protein (PP2-B13 and PP2-B15) were strongly down-regulated in response to CTV-B2 (Table 1), but up-regulated by infection with CaLas-B232. The expression of PP2-B10 gene was up-regulated by CTV-B6 and CaLas-B232 (Table 1). Because these phloem proteins move with assimilate from source to sink tissues [19], they may be markers for altered source sink flow of assimilates. Transcripts encoding NPR3 and NPR4, participating in cell organization, were also up-regulated responding to CaLas-B232 (Table 1). Transcripts encoding CYP83B1, an enzyme involved in the biosynthesis of glucosinolates and defense response by callose deposition, were up-regulated by both CTV-B6 and CaLas-B232 [20] (Table 1).

Table 1 DETs involved in cell wall modification and callose deposition in Citrus sinensis responding to CTV-B2, CTV-B6 and CaLas-B232

Transcription factors (TFs)

The expression of transcription factors, including WRKYs, B-box, MYB, bHTH, heat shock transcription factors (HSF), CDR1, pseudo-response regulator (PRR), ABR1, and jumonji domain-containing protein (JMJD5), were differentially modulated in response to the three pathogens. This was especially the case in trees infected by CTV-B6 (Additional file 1: Figure S4B). Transcripts encoding PRR5 and PRR7 were up-regulated responding to all three pathogens (Table 2). A transcript encoding PRR3 was up-regulated by infection with CTV-B6 and CaLas-B232 (Table 2). Transcripts encoding JMJD5, a histone demethylase, and CCT motif family protein were up-regulated in response to both CTV strains (Table 2). JMJD5 participates in circadian systems in both plants and animals [21]. The expression level of the late elongated hypocotyl (LHY) gene was down-regulated in response to the three pathogens, especially CTV-B6 (Table 2), consistent with down-regulation in flavedo (FF) and vascular tissue (VT) of both symptomatic ‘Hamlin’ and ‘Valencia’ sweet orange infected with CaLas [22].

Table 2 DETs involved in transcription factors and transportation in Citrus sinensis infected with CTV-B2, CTV-B6 and CaLas-B232

Plant-specific WRKY transcription factors with a DNA-binding domain are involved in plant senescence, plant defense and response to various environmental stresses [23, 24]. Expression of WRKY13 and WRKY 70 were especially up-regulated by infection with CTV-B6 and CaLas-B232. WRKY27 and WRKY50 were also up-regulated in response to both CTV-B6 and CaLas-B232 and WRKY28, WRKY33 and WRKY40 also were up-regulated in response to CaLas-B232 (Table 2). The expression of constitutive disease resistance (CDR1) gene, encoding an aspartic protease, was very strongly induced in response to both CTV-B6 and CaLas-B232, but not at all to CTV-B2 (Table 2). CDR1is involved in salicylic-acid-dependent inducible resistance responses and over-expression of CDR1 enhances resistance to bacterial pathogens in Arabidopsis [25]. A Myb-like HTH (helix-turn-helix) transcription factor was very highly down-regulated by infection with CTV-B2, but up-regulated with CaLas-B232 (Table 2). Transcripts encoding an integrase-type DNA-binding protein (ABR1), responsive to ethylene stimulus, were up-regulated by infection with CTV-B2 (Table 2).


Many transporter genes were differentially regulated in response to the three pathogens, especially CTV-B2 and CaLas-B232 (Additional file 1: Figure S4C). Transcripts encoding 4 zinc transporters (ZIP1, 4, 5 and 11), especially ZIP4 and ZIP5, were up-regulated in association with CaLas infection, but down regulated responding to CTV-B2, consistent with symptoms of zinc deficiency in HLB-diseased plants (Table 2). Transcripts encoding cyclic nucleotide gated channel 1 (CNGC1), involved in cyclic nucleotide or calcium transport, were differentially up-regulated in response to the three pathogens (Table 2). Transcripts encoding copper transporter 1 (COPT1) were up-regulated in response to all three pathogens (Table 2). Transcripts encoding ABC-2 type transporters were highly up-regulated by infection with CTV-B2 (Table 2). These ATP-binding cassette transporters typically function to export substrates from eukaryotic cells [26]. The expression of phosphate transporter 1;7 (PHT1;7) was induced by infection with CaLas-B232 (Table 2), consistent with a prominent role of phosphate deficiency in the HLB disease syndrome [27]. The expression of another phosphate transporter PHT1;4 was slightly up-regulated in response to both strains of CTV (Table 2). Transcripts encoding an oligopeptide transporter 7 (OPT7), involved in transport of small peptides [28], were up-regulated by infection with CTV-B6 and CaLas-B232 (Table 2).

Hormone metabolism

Transcripts encoding jasmonic acid carboxyl methyl transferase (JMT) were highly up-regulated responding to the three pathogens (Table 3). Transcripts encoding lipoxygenase 2 (LOX2), a key enzyme in the synthesis of JA [29], and involved in the JA signaling pathway, were highly up-regulated in response to CaLas-B232, but not to the other pathogens (Table 3). Different transcripts of the senescence-related gene (SRG1), involved in ethylene metabolism, were induced in response to the three pathogens (Table 3). The expression of ethylene responsive element binding factors (ERF1 and ERF104), involved in ethylene signal transduction, was also up-regulated by infection with CaLas-B232 (Table 3). A transcript encoding a gibberellin-regulated protein was down-regulated only by infection with CTV-B6 (Table 3).

Table 3 DETs involved in hormone and secondary metabolism pathways in Citrus sinensis plants infected with CTV-B2, CTV-B6 and CaLas-B232

Secondary metabolism

Many different transcripts encoding 2-oxoglutarate and Fe (II)-dependent (2OG-Fe (II)) oxygenases were highly up-regulated in response to all three pathogens (Table 3). Transcripts encoding terpene synthase (TPS21), involved in isoprenoid and terpenoid metabolism, were induced up to eight fold by infection with CTV-B2, and to a lesser extent infection by CaLas-B232, but were repressed by CTV-B6 (Table 3). These perturbations of the terpenoid pathways likely contribute to the different volatile components detected following infection by CTV or CaLas [30, 31]. The expression of O-methyltransferase (OMT), participating in flavonoid and phenylpropanoid metabolism, was up-regulated in response to both strains of CTV (Table 3). A transcript encoding geranylgeranyl pyrophosphatesynthase1 (GGPS1), also involved in isoprenoid biosynthetic process, was strongly up-regulated in response to CTV-B2 and less so to CaLas-B232 (Table 3). Transcripts encoding homogentisate phytyltransferase (VTE2) were up-regulated by infection with CTV-B2, but were down-regulated by infection with CTV-B6 and did not respond to CaLas-B232 (Table 3). VTE2 catalyzes the synthesis of tocopherol, which plays an important role in protecting chloroplast membranes. The expression of VTE2 was induced in symptomatic and repressed in asymptomatic leaves of sweet orange in response to CaLam [32]. Transcripts encoding terpenoidcyclases/proteinprenyl transferase were also significantly up-regulated by infection with CTV-B2 (Table 3). Laccase genes LAC4, LAC5 and LAC7 participating in simple phenol metabolism, were over-expressed in response to CaLas-B232 (Table 3). Transcripts encoding UDP-glucosyl transferase 74B1 (UGT74B1) involved in plant secondary metabolism and playing a potential role in stress or defense responses were up-regulated by both CTV-B2 and CaLas-B232 [33] (Table 3).


Plant receptor-like proteins/kinases (RLPs/RLKs) regulate developmental processes and respond to pathogens and abiotic stresses [34]. A large number of transcripts encoding RLRs/RLKs, including leucine rich repeat (LRR) kinases, wall-associated kinases (WAKs) and Domain of Unknown Function 26 (DUF26) kinases, were differentially expressed in response to all three pathogens, and very notably to CaLas-B232 (Additional file 1: Figure S4F). The expression of LRR kinases were up-regulated responding to CTV-B6 and CaLas-B232 but not to CTV-B2, consistent with their role in defense response [34]. These included RLP6, RLP9, RLP13, RLP14, RLP15 and RLP45, which were consistently up-regulated by infection with CaLas-B232 and CTV-B6 (Table 4). Transcripts related to wall associated kinases were differentially expressed in response to CTV-B2 and CaLas-B232, including WAK2 and WAK7 (Table 4). Notably, the expression of two members of DUF26 RLKs, cysteine-rich RLK2 (CRK2) and receptor-like kinase1 (RLK1), was significantly induced by infection with CTV-B2 (Table 4). The expression of E1-E2 ATPase was up-regulated in response to all three pathogens. The expression of EF-hand binding protein and calmodulin (CaM) was induced responding to CTV-B6 and CaLas-B232 (Table 4). Transcripts encoding cold-circadian rhythm RNA binding-like protein (CCR-like, CCL) were up-regulated in response to all three pathogens, especially CTV-B6, approximately five fold (Table 4).

Table 4 DETs involved in signaling, stress response and miscellaneous pathways in Citrus sinensis plants infected with CTV-B2, CTV-B6 and CaLas-B232

Stress response

Many disease resistance (R) genes were differentially regulated in response to CTV-B2 and CaLas-B232, and fewer to CTV-B6 (Additional file 1: Figure S4G), including CC (TIR)-NB-LRR (coiled-coil/Toll-interleukin1 receptor-nucleotide binding-Leucine rich repeat) and NB-ARC (nucleotide binding-adaptor shared by APAF-1, certain R gene products and CED4) domain-containing disease resistance proteins, the largest and well-characterized group of R genes. Transcripts encoding heat shock proteins (HSPs) were especially regulated in response to CTV-B6 and CaLas-B232 (Table 4). The expression of Kunitz-type protease inhibitor (PI), located in the cell wall and responsive to peroxidase, was significantly up-regulated only in response to CaLase-B232 (Table 4). Remarkably, a transcript encoding accelerated cell death 2 (ACD2) was highly up-regulated only responding to CTV-B2 (Table 4).

Oxidation/reduction processes

Transcripts related to the peroxidase super family were also differentially regulated in response to the three pathogens (Additional file 1: Table S9). Peroxidase 2 (PA2) mRNA accumulated in response to CTV-B2 (Table 4). Transcripts encoding WCRKC thioredoxin1 (WCRKC1), likely involved with redox homeostasis in the chloroplast, were greatly up-regulated only responding to CTV-B6 (Table 4). Transcripts encoding oxidoreductase, a zinc-binding dehydrogenase family protein, were up-regulated in response to CTV-B6 and CaLas-B232 (Table 4).

Development, protein degradation and miscellaneous functions

Transcripts encoding thiazole biosynthetic enzyme (THI1, TZ, THI4) were up-regulated responding to all three pathogens, especially CTV-B6 (Additional file 1: Table S9). Thiamine biosynthetic genes respond to bacteria, and also act in starch and maltose metabolism [35]. These enzymes may contribute to the accumulation of starch in leaves of diseased trees. Transcripts encoding gigantea (GI) and NAC domain containing protein 2 involved in the development and strengthening of secondary cell walls were strongly up-regulated in response to all three pathogens (Table 4). A set of transcripts involved in protein degradation were also differentially regulated in response to the three pathogens (Additional file 1: Table S9), consistent with a previous result that protein degradation and misfolding processes were activated in citrus fruits infected with CaLas [36]. A larger set of transcripts encoding UDP-glycoyl transferase and cytochrome P450 were differentially regulated in response to all three pathogens (Additional file 1: Table S9). Both UDP-glycoyl transferase and cytochrome P450 are involved in secondary metabolism, responses to biotic and abiotic stress and oxidation/reduction process. A transcript encoding a concanavalin A-like (ConA-like) protein with a jacalin-like lectin domain was very highly up-regulated by infection with CTV-B2 (Table 4). Jacalin-related proteins were demonstrated as components in plant defense system [37].


Recently, global analysis of gene expression has been applied widely to study plant-pathogen interactions. Though studies have been carried out to identify DETs in citrus infected with CaLas and CTV by microarray or RNA-seq, these studies have been conducted with each pathogen separately or in the context of susceptible vs. tolerant plant genotypes [9, 10, 12, 38, 39]. In the present study, RNA-seq was used to identify DETs in susceptible ‘Valencia’ sweet orange after CaLas and CTV infection, with all plants grown and inoculated together in the same greenhouse. We found many biological processes and pathways affected by the three pathogens. The number of DETs reflected the complexity of symptoms in response to the three pathogens. The transcriptome profiles were more similar between CTV-B6 (VT genotype) and CaLas-B232 than between the two CTV genotypes CTV-B2 (T30 genotype) and CTV-B6 or between CTV-B2 and CaLas-B232. This is consistent with similar symptoms. CTV-B6 causes severe stem pitting [2, 3], which is a major disturbance in the phloem system and CaLas-B232 induces typical symptoms of HLB which include degeneration of phloem cells. Both disease syndromes also include yellowing and cupping of leaves caused by the failure of phloem transport and starch accumulation in the leaves. These symptoms are not induced by CTV-B2, although the pathogen is fully systemic in the host.

Circadian system affected by all three pathogens

The circadian system, an endogenous mechanism vital to organisms, allows them to coordinate physiological changes to predictable day/night cycles and plays extremely important roles in reproductive development, germination, growth and carbohydrate metabolism processes throughout the life cycle of the plant [40, 41]. The circadian rhythm regulates the expression of 31 % of the Arabidopsis genome and have been found to be significantly disrupted by abiotic and biotic stress [42]. The circadian clock controls the cotton plant’s (Gossypium hirsutum) sensitivity to cold [43] and also regulates the transcription levels of some pathogen response related genes in Arabidopsis [44]. Disruption of the circadian rhythms by CTV or CaLas has not been reported previously, but off season flowering is observed in citrus groves in response to CTV and HLB [4547]. In our study, the circadian system was heavily affected by the three pathogens (Fig. 7), especially CTV-B6. Circadian clock-associated components PRR5 and PRR7 have been proven to play important roles in the circadian rhythm by repressing the expression of CCA1 (circadian clock-associated 1) and LHY (late elongated hypocoty) [48], consistent with the up-regulation of PRR and down-regulation of LHY genes in our study. CCA1 and LHY are homologous to Myb-related DNA-binding proteins and the circadian clock components form an auto-regulatory feedback loop at the transcription levels [49, 50] (Fig. 7). Accumulation of PRR1/TOC1 (timing of cab expression 1) in late day and early night promotes the rise of CCA1 and LHY proteins during the early and mid day, which in turn represses transcription of PRR1/TOC1. TOC1 also acts in concert with JMJD5 to regulate the circadian clock [21]. PRR family proteins and C-terminal (CCT) motif proteins are homologous and contribute to the control of the circadian rhythm that controls photoperiodic flowering [51]. The circadian clock controls the stability of the CCL mRNA at different times of day [52] as well as the expression of the GIGANTEA (GI) gene. The CCL gene was up-regulated in juice vesicles (JV) of ‘Hamlin’ and ‘Valencia’ by infection with CaLas [22], consistent with our study. The GI gene encodes a protein with six putative membrane-spanning domains that regulates both photoperiodic flowering and the deposition of wall ingrowths in phloem parenchyma transfer cells in A. thaliana [11, 18, 53, 54]. The formation of wall ingrowths in transfer cells plays an important role in phloem biology by producing physical barriers to limit pathogen invasion [55]. GI is also necessary to maintain the appropriate period of CCA1 and LHY gene expression and the circadian amplitude in the outer feedback loop in A. thaliana [56]. Therefore, the change of gene expression induced by all three pathogens, but especially by CTV-B6, indicated the alteration of the circadian cycle in the affected plants. Because of the large number of genes coordinated by the circadian clock [42], plants infected with these pathogens may have other abnormal gene expression patterns.

Fig. 7

Circadian rhythm pathway regulated by CTV-B2, CTV-B6 and CaLas-B232. The green boxes represent genes that were differentially expressed in response to all of the three pathogens

Cell wall modification defense response

Plant cell walls are composed of cellulose, hemicelluose and pectins and play a vital defensive role in plants against pathogens. Pectin is found in primary cell walls and the middle lamella, where it contributes to cell wall structure and binds cells together [57, 58]. Extensive disassembly of pectin is associated with tissue deterioration, leaf abscission [59], as well as swelling of the middle lamella [60], which was observed in cell walls surrounding sieve elements (SEs) at early stages of infection with CaLas [18]. Pectin hydrolysis often occurs in response to soft rot bacteria and fungi and pectin methylesterase inhibitors (PMEI) play a role in defense by regulation of cell wall susceptibility to microbial endopolygalacturonases. The expression of PMEI was also up-regulated in ‘Valencia’ leaves in association with CaLas [10], as well as in present study. The over expression of PMEI may also be associated with the interaction with Phytophthora spp. [14, 15] or opportunistic fungi.

Xyloglucan endotransglycosylase (XET) is an enzyme involved in the synthesis of xyloglucan, the predominant hemicellulose in the cell walls. Transcripts encoding XET were up-regulated in C. sinensis leaves [10] but down-regulated in C. sinensis stems infected with CaLas [39]. XET plays a profound role during the process of fruit ripening [61], which is impaired by CaLas. A. thaliana XTH33, encoding a xyloglucan endotransglycosylase/hydrolase (XTH), protects plants against aphids by modifying the cell wall [62]. Tomato leXTH1 mRNA accumulated in response to dodder, Cuscuta reflexa [63], and the expression of Medicago truncatula MtXTH1 gene was induced in association with Glomus and Gigaspora in roots [64]. Genes related to cellulose synthesis, CSLC7 and CSLD4, were down-regulated in stems in association with CaLas [39]. In the current study, the expression of cellulose synthase-like genes CSLB3, CSLA2 and COBL7, was strongly up-regulated by infection with CTV-B2, but not by infection with CaLas and CTV-B6. CSLA proteins have mannan/glucomannan synthase activity and mannans have vital roles in metabolic networks that regulate cell wall structure, carbohydrate storage and other cellular processes in plants [65, 66]. AtCSLA7 is essential for synthesis of a cell wall β-linked polysaccharide throughout the plant [67]. Both CSLH from barley and CSLF from rice mediate the synthesis of cell wall (1,3 and 1,4)-β-D-glucans in transgenic A. thaliana [68, 69]. COBRA-like (COBL) genes, predicted to encode a putative glycosylphophatidylinositol (GPI)-anchored protein in A. thaliana [70], play an important role in cell wall maintenance and biosynthesis, controlling the orientation of cell expansion and modulation of cellulose deposition in the cell walls of roots [71, 72]. Another COBRA-like protein, Brittle culm1 (BC1), provides mechanical strength for rice by regulating the biosynthesis of secondary cell walls [73]. Therefore, the plant responded to CTV-B2, but not to CTV-B6 and CaLas by activating enzymes that would work together to strengthen the cell wall, which may contribute to the reinforcement of physical barriers to restrict the invasion of CTV-B2.

Phloem damage and defense response

PP2, one of the well-known sieve element (SE) proteins, was found in plugged sieve pores and was described in cucurbit phloem sap [74]. Both CTV-B6 and CaLas-B232 damage phloem tissue and similar sets of genes may be altered by both pathogens. Several studies have reported that transcripts encoding phloem proteins (PP2-B10, PP2-B14 and PP2-B15) were up-regulated in symptomatic leaves in response to CaLas [11, 22, 32]. PP2 accumulated at the sieve plate linked to symptom development in response to CaLas [11]. Cucumber csPP2 interacts with the Hop stunt viroid (HSVd) RNA and may facilitate its long-distance movement [75, 76]. PP2 also can interact with polyadenylated mRNAs [75]. The altered expression of PP2 genes was tested and confirmed by RT-qPCR (Additional file 1: Table S10). The closure of sieve pores by PP2 is normally accompanied by the accumulation of callose on sieve plates under abiotic or biotic stress [77].

Callose is widely distributed in higher plants and is deposited at cell plates during cytokinesis [78] and in plasmodesmata, and contributes to the regulation of movement of molecules from cell-to-cell [79]. PP2 and callose were deposited at the sieve plates in apple trees infected with the phloem-limited pathogen ‘Ca. Phytoplasma mali’ [80]. Accumulation of PP2 and callose proteins may provide a physical barrier against further invasion of CTV-B6 and CaLas-B232, and contribute to long-distance signaling [10, 19, 75]. The over expression of callose and PP2 genes is consistent with the phloem damage and defense mechanisms activated by CTV-B6 and CaLas-B232. The clear down regulation of PP2 by infection with CTV-B2 suggests that CTV-B2 was not perceived by the host plant in a way that triggered the defense response.

Ions and nutrient deficiency

The expression of both EF-hand calcium-binding and calmodulin-binding proteins was down-regulated in stems of C. sinensis inoculated with CaLas [39], while the expression of EF-hand calcium-binding protein was up-regulated about 5 fold in C. sinensis 13–17 weeks after inoculation with CaLas [10]. In the current study, the expression of both genes was also up-regulated (Table 4). The differences among these studies might be attributed to the different stages of infection or tissues used. CNGC (cyclic nucleotide gated channel) mRNA accumulated in flavedo (FF), vascular tissue (VT) and juice vesicles (JV) of symptomatic ‘Hamlin’ and ‘Valencia’ sweet orange infected with CaLas [22]. Our transcriptome profiles showed up-regulation of EF-hand calcium-binding protein, calmodulin-binding protein, CNGCs and ATPase E1-E2 proteins in young leaves in response to CaLas. ATPase E1-E2 accumulated in response to CTV-B2 and EF-hand calcium-binding proteins accumulated in response to CTV-B6. ATPase E1-E2 is involved in Ca2+-translocation and CNGCs are involved in Ca2+-signaling processes and both maintain [Ca2+] homeostasis in the plant [81, 82]. In Arabidopsis, CNGC1 conducts Ca2+ into cells [83] and AtCNGC2 responds to pathogen recognition by increasing [Ca2+]cyt levels [84]. AtCNGC4 is also an essential component of the signaling pathway leading to the hypersensitive response [85]. Calmodulin (CaM) is a major Ca2+ signal transducer in both animals and plants [86] and is involved in Ca2+-mediated plant disease resistance responses [87]. Soybean CaM genes SCaM4 and ScaM5 were induced within 30 min by pathogens [87]. Evidence of altered calcium metabolism is also demonstrated by the formation of crystal idioblasts in CaLas infected plant tissues [18]. These idioblasts are believed to contain calcium oxalate (CaOX), which can serve dual functions in the regulation of calcium concentration and plant defense. An influx of Ca2+ can also activate Ca2+-dependent defense responses such as PP2 aggregation and callose deposition [80]. The accumulation of EF-hand calcium-binding protein, CaM, ATPase E1-E2 and CNGC and the formation of CaOX crystals suggest imbalances in calcium metabolism after infection by the three pathogens. Ca2+-dependent defense mechanisms were activated by the three pathogens and play more significant roles in defense response to CaLas-B232.

Zinc transporters such as ZIP1, ZIP2, ZIP3, ZIP4, ZIP5 and ZIP11, are induced in plants in response to zinc deficiency. ZIP4 and ZIP5 mediate zinc uptake from the rhizosphere and are responsible for translocation of zinc within the plant [88, 89]. Citrus plants affected by CaLas present nutritional deficiency symptoms typically associated with zinc deficiency [4, 5]. Up-regulation of ZIP1 was observed in transcriptome profiles of ‘Madame Vinous’ and ‘Navel’ sweet oranges inoculated with CaLas [9], and over expression of ZIP5 was also found in the susceptible genotype ‘Cleopatra mandarin’ but not in the tolerant genotype ‘US-897’ inoculated with CaLas [38]. In another study, transcripts encoding ZIP1, ZIP4 and ZIP5 were up-regulated in ‘Pera’ sweet orange inoculated with CaLam [32], all of which is consistent with the over expression of four zinc transporters (ZIP1, ZIP4, ZIP5 and ZIP11) in our study. The up-regulation of zinc transporters in susceptible citrus genotypes infected with CaLas or CaLam in the late stages is consistent with the zinc deficiency symptoms seen on leaves, but is also observed in young leaves without symptoms. If CaLas-infected plants experience dramatic fibrous root decline before foliar symptoms appear [14, 15], this would limit the uptake of sufficient microelements such as Zn and Cu from the rhizosphere. CaLas also assimilates zinc from the plant phloem sap with assistance of zinc uptake ABC transporters [90], further contributing to the zinc deficit in the host.

Zinc deficiency leads to a wide range of ultrastructural problems with chloroplasts, mitochondria, cytoplasm and aberrant regulation of DNA-transcription and RNA-processing in plants [91], followed by foliar symptoms on leaves and delayed maturation of fruits. Transcription of Zn and Cu transporters would then be up-regulated in response to CaLas to meet the deficits of Zn, Cu and other micronutrients. This also may explain why HLB disease severity was reduced by the application of various nutritional products including micronutrients on the foliage extended the vigor and productivity of CaLas infected trees [6]. CaLas encodes a zinc-dependent metalloprotease (Accession Number: ACT57246; GI: 204040450) similar to RseP-like site-2 proteases (S2P), which need zinc to enhance stability [92]. The introduction of zinc-dependent metalloproteases into organisms can cause impairment of physiological functions that lead to disease or death [93]. The zinc-dependent metalloprotease in CaLas may explain the up-regulation level of zinc transporters to a greater degree than transporters of other microelements in response to CaLas, and also the higher expression level of zinc transporters than in plants infected with CTV-B6. In contrast, ZIP4 and ZIP5 were notably down-regulated in asymptomatic sweet orange trees infected by CTV-B2 in this study. Taken together, the differential expression of Zn, Cu and Ca transporters is evidence of impaired homeostasis in plants caused by the three pathogens. The greater expression level of zinc transporters in response to CaLas than CTV-B6 may be attributed to root decline and usage of zinc by CaLas.

Hormone-mediated defense response

The JMT gene was induced in young leaves and in mature fruits infected with CaLas [94], consistent with up-regulation in our study (Table 3). JMT catalyzes the synthesis of methyl jasmonate (MeJA), which plays an important role in JA-regulated plant defense responses [95]. In our study, the up-regulation of LOX2 encoding a lipoxygenase involved in the JA pathway responding to CaLas is in accordance with previous study [10]. The transcription level of LOX was also increased in Citrus reticulata inoculated with Xylella fastidiosa at the early stage [96]. Many WRKY transcription factors are involved in hormone-mediated signaling pathways under biotic and abiotic stress [97]. WRKY70 coordinates the cross talk of SA-and JA-mediated signals in plant defense, and over expression of WRKY70 increased resistance to pathogens and influenced expression of SA induced pathogenesis-related genes [98]. WRKY6, WRKY40 and WRKY70 were up-regulated about five fold in symptomatic leaves of sweet orange [32] and WRKY70 was expressed at a higher level in fruits than in leaves in response to CaLas [94]. The accumulation of WRKY13 transcripts affected transcription of more than 500 genes, some that functioned in disease resistance [99]. The expression of WRKYs may also link to extensive transcriptional repression and activation by its own WRKY members [97]. Hormone-mediated defense responses were activated by all three pathogens. However, hormone-mediated defense responses were stronger in response to CTV-B6 and CaLas-B232 than the mild strain CTV-B2 which did not trigger these defense responses to nearly the same degree.

Secondary metabolism related defense

Plant secondary metabolites play an important role in defense responses to herbivores, pests and pathogens [100]. 2OG-Fe (II) oxygenases, involved in the production of plant secondary metabolites such as flavonoids, were induced by the three pathogens. DMR6, encoding a 2OG-Fe (II) oxygenase, plays a role in plant defense in A. thaliana [101]. The expression of 2OG-Fe (II) oxygenase was up-regulated in the susceptible genotypes (C. sinensis and C. reticulata), and also extremely up-regulated in the tolerant genotype (US-879, C. reticulata x P. trifoliata) infected with CaLas [10, 38]. Terpene synthases and terpenoid cyclases were highly induced in response to CTV-B2 and both are involved in the production of terpenes, the largest class of secondary metabolites [102], which are important factors in resistance to several pathogens and pests [103]. Hence, the expression of secondary metabolites was triggered by all three pathogens, but to a much greater degree in response to CTV-B2.

Pathogenesis related defense responses

Receptor-like proteins/kinases (RLPs/RLKs) are important factors in hormone signaling in response to abiotic and biotic stress in plants. More RLPs/RLKs were induced in response to CTV-B6 and CaLas-B232 than to CTV-B2. However some RLPs/RLKs, such as CRK2 and RLK1, were only induced in response to CTV-B2. These are cysteine-rich RLKs (CRKs) of the DUF26 group. Three CRKs (CRK4, CRK19 and CRK20), closely related to CRK5, were also induced by SA treatment and infection with Pseudomonas syringae and all activate rapid cell death in transgenic Arabidopsis [104]. CRK genes, PRK1 to PRK4 from Potato and PK20-1 from common bean were differentially regulated in response to symbionts, elicitor molecules and pathogens [105]. More pathogenesis-related (PR) genes were up-regulated in sweet orange inoculated with CaLas-B232 than with CTV-B2 and CTV-B6.

Transcripts encoding Kunitz-type trypsin inhibitors were up-regulated in response to CaLas. The Kunitz trypsin inhibitors are responsive to wounding [106] and protease inhibitors from plants may also inhibit the growth of a variety of pathogenic bacteria and fungi. An important gene to note is ACD2, up-regulated only in response to CTV-B2. The product of ACD2 has extensive similarities to red chlorophyll catabolite reductase, which is involved in the degradation of the porphyrin component of chlorophyll (Chl) [107]. Over expression of ACD2 in plants induced many characteristics of the hypersensitive response (HR), including cell wall modification, accumulation of mRNA of PR genes and resistance to bacteria [108]. The up-regulation of ACD2 did not restrict the growth of Pseudomonas syringae but reduced disease symptoms [108, 109], which may explain symptom suppression in response to CTV-B2. Pathogenesis-related defense responses were induced by all three pathogens, and plants used more kinase signaling and stress response pathways in response to CaLas-B232 than to either strain of CTV.


More genes and pathways were regulated in response to CaLas-B232, and to a lesser extent CTV-B6, than to CTV-B2, consistent with the symptoms caused by corresponding pathogens. Gene expression patterns of trees infected with CTV-B6 and CaLasB232 were more similar than those of infected with CTV-B2 and CTV-B6. These gene expression patterns in the asymptomatic leaves sampled lead to similar physiological changes, such as phloem damage, root decline, yellowing and cupping of leaves caused by both CTV-B6 and CaLas-B232. The circadian rhythm system of sweet orange plants was heavily perturbed by all three pathogens, most notably by CTV-B6, and ion balances were also disrupted by all three pathogens, but especially by CaLas-B232. Many defense responses were triggered by all three pathogens but with different preferences: Defense responses based on cell wall modification and secondary metabolism were activated by all three pathogens, but were more prominent in response to CTV-B2. Defense responses modulated by transcription factors responded most significantly to CTV-B6. Ca2+-dependent, kinase signaling and pathogenesis-related defense responses functioned more actively in response to CaLas-B232 than to either strain of CTV. DETs regulated by CTV-B2, but not by CTV-B6 or CaLas, are likely to be important to protect the plant from CTV-B2, based on perception of the pathogen by the host and efficient triggering of systemic defense responses.


Inoculation of experimental trees with CTV and CaLas

CTV mild strain B2 (T30 genotype, Florida), severe strain B6 (VT genotype, =SY568, California) and CaLas strain B232 (Thailand) are maintained in planta as part of the Exotic Pathogens of Citrus Collection (EPCC) at the USDA-ARS Beltsville Agricultural Research Center (BARC) in Beltsville, MD. Experimental seedling trees were ‘Valencia’ sweet orange. Each pathogen was inoculated into three replicate trees by inverted T-bud grafting when the trees were actively growing. Three trees were also mock-inoculated with healthy buds as the control. The trees were maintained in a greenhouse with supplemental fertilization provided with the irrigation water as described previously [110].

Detection of CTV and CaLas with PCR

In order to confirm that the plants were successfully inoculated and infected with the phloem-limited pathogens, RNA and DNA were extracted from leaf midribs. DNA was extracted from 0.1 g of leaf midrib using the Plant DNeasy® Mini Kit (Qiagen, Valencia, CA) according to the manufacturer’s instructions. qPCR assays were performed as described [111] with a SmartCycler (Cepheid, Sunnyvale, CA). Total RNA was extracted with Trizol (Invitrogen, Carlsbad, CA) and treated with TURBO™ DNase (AM2239, Ambion) according to the manufacturer’s instructions. cDNA synthesis was done with the Thermoscript RT-PCR System for First-Strand cDNA Synthesis (Invitrogen). CTV cDNA (2 μl as template) was amplified in a 25 μl reaction volume containing 10 X reaction buffer, 10 mM dNTP, 25 mM MgCl2, 10 μM of each primer [112] and 0.3 U of Taq polymerase. Amplification was carried out with a PTC-200 thermal cycler (MJ Research) as follows: 94 °C for 3 min, then 94 °C for 30 s, 57 °C for 40 s, and 72 °C for 40 s. After 30 cycles the reactions were incubated an additional 10 min at 72 °C prior to storage at 4 °C. Products were analyzed by electrophoresis on 1.5 % agarose gels and stained with Gel Red (Biotium, Hayward, CA). Plants with confirmed infections were used for transcription analysis.

Isolation and quantification of total RNA

After infection by the pathogens had been confirmed by PCR and qPCR, but before symptoms became evident, RNA samples were collected in the morning (8–10 AM) from young, soft, not fully expanded leaves of uniform size at 32 weeks after inoculation (Additional file 1: Figure S5). Total RNA was extracted and treated as described above. A portion of the RNA samples from three PCR-positive independent biological trees of each pathogen and from mock-inoculated trees were pooled. The quantity and quality of total RNA was determined with the Qubit 2.0 Fluorometer (Invitrogen) and the Bioanalyzer 2100 (Agilent Technologies, Santa Clara, CA), respectively. A total of 30 μg of RNA (RIN ≥ 7.5) [113] from healthy ‘Valencia’ sweet orange and ‘Valencia’ sweet orange infected with CTV-B2, CTV-B6 and CaLas-B232 were sent to SeqWright, Inc. (Houston, TX, USA) for paired-end sequencing of polyA selected mRNA with the Illumina HiSeq 2500 platform. The remainder of the RNA extracts was retained as replicates for RT-qPCR assays of selected mRNA targets.

Statistical analyses

Paired-end reads (2 X 100 nucleotides) for the four citrus cDNA libraries were generated with the Illumina HiSeq 2500 and mapped to Citrus sinensis ‘Valencia’ reference genome [13] with Bowtie ( Differentially expressed transcripts (DETs) were identified for CTV-B2, CTV-B6 and CaLas-B232 infected citrus libraries against the mock-inoculated library with the Gfold program [114]. Gfold is especially useful when replicate data is not available. Gfold normalizes the expression level of genes on reads per kilobase per million reads (RPKM) and tests for significance with the Poisson distribution. P < 0.01 and log2 fold change (log2FC) ≥ │1│were set as cut-off values. DETs were assigned into GO (Gene Ontology) categories, including biological processes (BPs), molecular functions (MFs) and cellular components (CCs) based on A. thaliana orthologs by PlantGSEA [115] (Gene Set Enrichment Analysis) [115] with a hypergeometric test (p-values < 0.05). Functional annotation of DETs was also performed with Mapman [116] based on Citrus sinensis orthologs. The raw RNA-Seq data has been deposited at the sequence read archive of NCBI under the project number SRP067360.

RT-qPCR analysis

A total of 38 genes were assayed by RT-qPCR with RNA preparations from true biological replicates to validate the Gfold results. Genes were selected based on their predicted function in disease symptom development and fold changes (log2FC) estimated by Gfold. Gene-specific primers were designed with calculated melting temperatures of 60 °C ± 5 °C by primer (Additional file 1: Table S11; Integrated DNA Technologies, Coralville, IA). Melting curve analysis was performed to ensure that a single product was amplified. The 2-Ct method was applied for relative quantification of gene expression [117]. The correlation between Gfold and RT-qPCR was performed with SPSS 16.0.



biological process


Candidatus Liberibacter asiaticus’


cellular component


Citrus tristeza virus


differentially expressed transcript


fold change


Gene Ontology


Gene Set Enrichment Analysis




molecular function


Reads per KB per million


transcription factor


  1. 1.

    Moreno P, Ambros S, Albiach-Marti MR, Guerri J, Pena L. Citrus tristeza virus: A pathogen that changed the course of the citrus industry. Mol Plant Pathol. 2008;9(2):251–68.

    CAS  Article  PubMed  Google Scholar 

  2. 2.

    Ruiz-Ruiz S, Moreno P, Guerri J, Ambrós S. Discrimination between mild and severe Citrus tristeza virus isolates with a rapid and highly specific real-time reverse transcription-polymerase chain reaction method using TaqMan LNA probes. Phytopathology. 2009;99(3):307–15.

    CAS  Article  PubMed  Google Scholar 

  3. 3.

    Dawson WO, Garnsey SM, Tatineni S, Folimonova SY, Harper SJ, Gowda S. Citrus tristeza virus-host interactions. Front Microbiol. 2013;4.

  4. 4.

    Bové JM. Huanglongbing: a destructive, newly-emerging, century-old disease of citrus. J Plant Pathol. 2006;88(1):7–37.

    Google Scholar 

  5. 5.

    Wang N, Trivedi P. Citrus huanglongbing: a newly relevant disease presents unprecedented challenges. Phytopathology. 2013;103(7):652–65.

    Article  PubMed  Google Scholar 

  6. 6.

    Rouse B. Rehabilitation of HLB infected citrus trees using severe pruning and nutritional sprays. In: Proceedings of the Florida State Horticultural Society. 2013. p. 51–4.

    Google Scholar 

  7. 7.

    Stansly PA, Arevalo HA, Qureshi JA, Jones MM, Hendricks K, Roberts PD, Roka FM. Vector control and foliar nutrition to maintain economic sustainability of bearing citrus in Florida groves affected by huanglongbing. Pest Manag Sci. 2014;70(3):415–26.

    CAS  Article  PubMed  Google Scholar 

  8. 8.

    Hammond-Kosack KE, Jones J. Resistance gene-dependent plant defense responses. Plant Cell. 1996;8(10):1773–91.

    CAS  PubMed Central  Article  PubMed  Google Scholar 

  9. 9.

    Fan J, Chen C, Yu Q, Brlansky RH, Li ZG, Gmitter FG. Comparative iTRAQ proteome and transcriptome analyses of sweet orange infected by ‘Candidatus Liberibacter asiaticus’. Physiol Plant. 2011;143(3):235–45.

    CAS  Article  PubMed  Google Scholar 

  10. 10.

    Albrecht U, Bowman KD. Gene expression in Citrus sinensis (L.) Osbeck following infection with the bacterial pathogen ‘Candidatus Liberibacter asiaticus’ causing Huanglongbing in Florida. Plant Sci. 2008;175(3):291–306.

    CAS  Article  Google Scholar 

  11. 11.

    Kim J-S, Sagaram US, Burns JK, Li J-L, Wang N. Response of sweet orange (Citrus sinensis) to ‘Candidatus Liberibacter asiaticus’ infection: microscopy and microarray analyses. Phytopathology. 2009;99(1):50–7.

    Article  PubMed  Google Scholar 

  12. 12.

    Gandía M, Conesa A, Ancillo G, Gadea J, Forment J, Pallás V, Flores R, Duran-Vila N, Moreno P, Guerri J. Transcriptional response of Citrus aurantifolia to infection by Citrus tristeza virus. Virology. 2007;367(2):298–306.

    Article  CAS  PubMed  Google Scholar 

  13. 13.

    Gmitter Jr FG, Chen C, Machado MA, de Souza AA, Ollitrault P, Froehlicher Y, Shimizu T. Citrus genomics. Tree Genet Genomes. 2012;8(3):611–26.

    Article  Google Scholar 

  14. 14.

    Johnson E, Wu J, Bright D, Graham J. Association of ‘Candidatus Liberibacter asiaticus’ root infection, but not phloem plugging with root loss on Huanglongbing-affected trees prior to appearance of foliar symptoms. Plant Pathol. 2013;63:290–8.

    Article  Google Scholar 

  15. 15.

    Graham JH, Johnson EG, Gottwald T, Irey MS. Pre-symptomatic fibrous root decline in citrus trees caused by huanglongbing and potential interaction with Phytophthora spp. Plant Dis. 2013;97(9):1195–9.

    Article  Google Scholar 

  16. 16.

    Lionetti V, Raiola A, Camardella L, Giovane A, Obel N, Pauly M, Favaron F, Cervone F, Bellincampi D. Overexpression of pectin methylesterase inhibitors in Arabidopsis restricts fungal infection by Botrytis cinerea. Plant Physiol. 2007;143(4):1871–80.

    CAS  PubMed Central  Article  PubMed  Google Scholar 

  17. 17.

    An SH, Sohn KH, Choi HW, Hwang IS, Lee SC, Hwang BK. Pepper pectin methylesterase inhibitor protein CaPMEI1 is required for antifungal activity, basal disease resistance and abiotic stress tolerance. Planta. 2008;228(1):61–78.

    CAS  PubMed Central  Article  PubMed  Google Scholar 

  18. 18.

    Folimonova SY, Achor DS. Early events of Citrus Greening (Huanglongbing) disease development at the ultrastructural level. Phytopathology. 2010;100(9):949–58.

    Article  PubMed  Google Scholar 

  19. 19.

    Golecki B, Schulz A, Thompson GA. Translocation of structural P proteins in the phloem. Plant Cell Online. 1999;11(1):127–40.

    CAS  Article  Google Scholar 

  20. 20.

    Clay NK, Adio AM, Denoux C, Jander G, Ausubel FM. Glucosinolate metabolites required for an Arabidopsis innate immune response. Science. 2009;323(5910):95–101.

    CAS  PubMed Central  Article  PubMed  Google Scholar 

  21. 21.

    Jones MA, Covington MF, Ditacchio L, Vollmers C, Panda S, Harmer SL. Jumonji domain protein JMJD5 functions in both the plant and human circadian systems. Proc Natl Acad Sci. 2010;107(50):21623–8.

    CAS  PubMed Central  Article  PubMed  Google Scholar 

  22. 22.

    Liao H-L, Burns JK. Gene expression in Citrus sinensis fruit tissues harvested from huanglongbing-infected trees: comparison with girdled fruit. J Exp Bot. 2012;63(8):3307–19.

    CAS  PubMed Central  Article  PubMed  Google Scholar 

  23. 23.

    Yang P, Chen C, Wang Z, Fan B, Chen Z. A pathogen-and salicylic acid-induced WRKY DNA-binding activity recognizes the elicitor response element of the tobacco class I chitinase gene promoter. Plant J. 1999;18(2):141–9.

    CAS  Article  Google Scholar 

  24. 24.

    Chen W, Provart NJ, Glazebrook J, Katagiri F, Chang H-S, Eulgem T, Mauch F, Luan S, Zou G, Whitham SA. Expression profile matrix of Arabidopsis transcription factor genes suggests their putative functions in response to environmental stresses. Plant Cell Online. 2002;14(3):559–74.

    CAS  Article  Google Scholar 

  25. 25.

    Xia Y, Suzuki H, Borevitz J, Blount J, Guo Z, Patel K, Dixon RA, Lamb C. An extracellular aspartic protease functions in Arabidopsis disease resistance signaling. EMBO J. 2004;23(4):980–8.

    CAS  PubMed Central  Article  PubMed  Google Scholar 

  26. 26.

    Yazaki K. ABC transporters involved in the transport of plant secondary metabolites. FEBS Lett. 2006;580(4):1183–91.

    CAS  Article  PubMed  Google Scholar 

  27. 27.

    Zhao H, Sun R, Albrecht U, Padmanabhan C, Wang A, Coffey MD, Girke T, Wang Z, Close TJ, Roose M. Small RNA profiling reveals phosphorus deficiency as a contributing factor in symptom expression for Citrus Huanglongbing disease. Mol Plant. 2013;6(2):301–10.

    CAS  PubMed Central  Article  PubMed  Google Scholar 

  28. 28.

    Koh S, Wiles AM, Sharp JS, Naider FR, Becker JM, Stacey G. An oligopeptide transporter gene family in Arabidopsis. Plant Physiol. 2002;128(1):21–9.

    CAS  PubMed Central  Article  PubMed  Google Scholar 

  29. 29.

    Delker C, Stenzel I, Hause B, Miersch O, Feussner I, Wasternack C. Jasmonate biosynthesis in Arabidopsis thaliana-enzymes, products, regulation. Plant Biol. 2006;8(3):297–306.

    CAS  Article  PubMed  Google Scholar 

  30. 30.

    Aksenov AA, Pasamontes A, Peirano DJ, Zhao W, Dandekar AM, Fiehn O, Ehsani R, Davis CE. Detection of huanglongbing disease using differential mobility spectrometry. Anal Chem. 2014;86(5):2481–8.

    CAS  Article  PubMed  Google Scholar 

  31. 31.

    Cheung WH, Pasamontes A, Peirano DJ, Zhao W, Grafton-Cardwell EE, Kapaun T, Yokomi RK, Simmons J, Doll M, Fiehn O. Volatile organic compound (VOC) profiling of Citrus tristeza virus infection in sweet orange citrus varietals using thermal desorption gas chromatography time of flight mass spectrometry (TD-GC/TOF-MS). Metabolomics. 2015;11(6):1514–25.

  32. 32.

    Mafra V, Martins PK, Francisco CS, Ribeiro-Alves M, Freitas-Astúa J, Machado MA. ‘Candidatus Liberibacter americanus’ induces significant reprogramming of the transcriptome of the susceptible citrus genotype. BMC Genomics. 2013;14(1):247.

    CAS  PubMed Central  Article  PubMed  Google Scholar 

  33. 33.

    Douglas Grubb C, Zipp BJ, Ludwig‐Müller J, Masuno MN, Molinski TF, Abel S. Arabidopsis glucosyltransferase UGT74B1 functions in glucosinolate biosynthesis and auxin homeostasis. Plant J. 2004;40(6):893–908.

    Article  CAS  PubMed  Google Scholar 

  34. 34.

    McHale L, Tan X, Koehl P, Michelmore RW. Plant NBS-LRR proteins: adaptable guards. Genome Biol. 2006;7(4):212.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  35. 35.

    Goyer A. Thiamine in plants: aspects of its metabolism and functions. Phytochemistry. 2010;71(14):1615–24.

    CAS  Article  PubMed  Google Scholar 

  36. 36.

    Martinelli F, Uratsu SL, Albrecht U, Reagan RL, Phu ML, Britton M, Buffalo V, Fass J, Leicht E, Zhao W. Transcriptome profiling of citrus fruit response to Huanglongbing disease. PLoS ONE. 2012;7(5), e38039.

    CAS  PubMed Central  Article  PubMed  Google Scholar 

  37. 37.

    Xiang Y, Song M, Wei Z, Tong J, Zhang L, Xiao L, Ma Z, Wang Y. A jacalin-related lectin-like gene in wheat is a component of the plant defence system wheat is a component of the plant defence system. J Exp Bot. 2011;62(15):5471-83.

  38. 38.

    Albrecht U, Bowman KD. Transcriptional response of susceptible and tolerant citrus to infection with ‘Candidatus Liberibacter asiaticus’. Plant Sci. 2012;185:118–30.

    Article  CAS  PubMed  Google Scholar 

  39. 39.

    Aritua V, Achor D, Gmitter FG, Albrigo G, Wang N. Transcriptional and microscopic analyses of citrus stem and root responses to ‘Candidatus Liberibacter asiaticus’ infection. PLoS ONE. 2013;8(9), e73742.

    CAS  PubMed Central  Article  PubMed  Google Scholar 

  40. 40.

    McClung CR. Circadian rhythms in plants. Annu Rev Plant Biol. 2001;52(1):139–62.

    CAS  Article  Google Scholar 

  41. 41.

    McClung CR. Plant circadian rhythms. Plant Cell. 2006;18(4):792–803.

    CAS  PubMed Central  Article  PubMed  Google Scholar 

  42. 42.

    Rosa BA, Jiao Y, Oh S, Montgomery BL, Qin W, Chen J. Frequency-based time-series gene expression recomposition using PRIISM. BMC Syst Biol. 2012;6(1):69.

    PubMed Central  Article  PubMed  Google Scholar 

  43. 43.

    Rikin A, Dillwith JW, Bergman DK. Correlation between the circadian rhythm of resistance to extreme temperatures and changes in fatty acid composition in cotton seedlings. Plant Physiol. 1993;101(1):31–6.

    CAS  PubMed Central  PubMed  Google Scholar 

  44. 44.

    Schaffer R, Landgraf J, Accerbi M, Simon V, Larson M, Wisman E. Microarray analysis of diurnal and circadian-regulated genes in Arabidopsis. Plant Cell Online. 2001;13(1):113–23.

    CAS  Article  Google Scholar 

  45. 45.

    Gopal K, Gopi V, Kalyani L, Sreelatha M, Sreenivasulu B. Symptom-based diagnosis of Huanglongbing (citrus greening) disease by PCR in sweet orange (Citrus sinensis Osbeck) and acid lime (Citrus aurantifolia Swingle). Arch Phytopathol Plant Protect. 2010;43(9):863–70.

    CAS  Article  Google Scholar 

  46. 46.

    Kyriakou A, Polycarpou D, Efstathiou A, Hadjinicoli A. Citrus tristeza virus in Cyprus. In: Proc 12th Conf IOCV IOCV, Riverside: 1993. 1993. p. 69–72.

    Google Scholar 

  47. 47.

    Halbert SE, Manjunath KL. Asian citrus psyllids (Sternorrhyncha: Psyllidae) and greening disease of citrus: a literature review and assessment of risk in Florida. Fla Entomol. 2004;87(3):330–53.

    Article  Google Scholar 

  48. 48.

    Nakamichi N, Kusano M, Fukushima A, Kita M, Ito S, Yamashino T, Saito K, Sakakibara H, Mizuno T. Transcript profiling of an Arabidopsis PSEUDO RESPONSE REGULATOR arrhythmic triple mutant reveals a role for the circadian clock in cold stress response. Plant Cell Physiol. 2009;50(3):447–62.

    CAS  Article  PubMed  Google Scholar 

  49. 49.

    Alabadı́ D, Oyama T, Yanovsky MJ, Harmon FG, Mas P, Kay SA. Reciprocal regulation between TOC1 and LHY/CCA1 within the Arabidopsis circadian clock. Science. 2001;293(5531):880–3.

    Article  PubMed  Google Scholar 

  50. 50.

    Alabadı́ D, Yanovsky MJ, Más P, Harmer SL, Kay SA. Critical Role for CCA1 and LHY in Maintaining Circadian Rhythmicity in Arabidopsis. Curr Biol. 2002;12(9):757–61.

    Article  PubMed  Google Scholar 

  51. 51.

    Onouchi H, Igeño MI, Périlleux C, Graves K, Coupland G. Mutagenesis of plants overexpressing CONSTANS demonstrates novel interactions among Arabidopsis flowering-time genes. Plant Cell Online. 2000;12(6):885–900.

    CAS  Article  Google Scholar 

  52. 52.

    Lidder P, Gutiérrez RA, Salomé PA, McClung CR, Green PJ. Circadian control of messenger RNA stability. Association with a sequence-specific messenger RNA decay pathway. Plant Physiol. 2005;138(4):2374–85.

    CAS  PubMed Central  Article  PubMed  Google Scholar 

  53. 53.

    Edwards J, Martin AP, Andriunas F, Offler CE, Patrick JW, McCurdy DW. GIGANTEA is a component of a regulatory pathway determining wall ingrowth deposition in phloem parenchyma transfer cells of Arabidopsis thaliana. Plant J. 2010;63(4):651–61.

    CAS  Article  PubMed  Google Scholar 

  54. 54.

    Fowler S, Lee K, Onouchi H, Samach A, Richardson K, Morris B, Coupland G, Putterill J. GIGANTEA: a circadian clock-controlled gene that regulates photoperiodic flowering in Arabidopsis and encodes a protein with several possible membrane-spanning domains. EMBO J. 1999;18(17):4679–88.

    CAS  PubMed Central  Article  PubMed  Google Scholar 

  55. 55.

    Chinnappa KSA, Nguyen TTS, Hou J, Wu Y, McCurdy DW. Phloem parenchyma transfer cells in Arabidopsis-an experimental system to identify transcriptional regulators of wall ingrowth formation. Front Plant Sci. 2013;4.

  56. 56.

    Park DH, Somers DE, Kim YS, Choy YH, Lim HK, Soh MS, Kim HJ, Kay SA, Nam HG. Control of circadian rhythms and photoperiodic flowering by the Arabidopsis GIGANTEA gene. Science. 1999;285(5433):1579–82.

    CAS  Article  PubMed  Google Scholar 

  57. 57.

    Caffall KH, Mohnen D. The structure, function, and biosynthesis of plant cell wall pectic polysaccharides. Carbohydr Res. 2009;344(14):1879–900.

    CAS  Article  PubMed  Google Scholar 

  58. 58.

    Pérez S, Rodriguez-Carvajal M, Doco T. A complex plant cell wall polysaccharide: rhamnogalacturonan II. A structure in quest of a function. Biochimie. 2003;85(1):109–21.

    Article  CAS  PubMed  Google Scholar 

  59. 59.

    Hadfield KA, Bennett AB. Polygalacturonases: many genes in search of a function. Plant Physiol. 1998;117(2):337–43.

    CAS  PubMed Central  Article  PubMed  Google Scholar 

  60. 60.

    Redgwell RJ, MacRae E, Hallett I, Fischer M, Perry J, Harker R. In vivo and in vitro swelling of cell walls during fruit ripening. Planta. 1997;203(2):162–73.

    CAS  Article  Google Scholar 

  61. 61.

    Baumann MJ, Eklöf JM, Michel G, Kallas ÅM, Teeri TT, Czjzek M, Brumer H. Structural evidence for the evolution of xyloglucanase activity from xyloglucan endo-transglycosylases: biological implications for cell wall metabolism. Plant Cell Online. 2007;19(6):1947–63.

    CAS  Article  Google Scholar 

  62. 62.

    Divol F, Vilaine F, Thibivilliers S, Kusiak C, Sauge MH, Dinant S. Involvement of the xyloglucan endotransglycosylase/hydrolases encoded by celery XTH1 and Arabidopsis XTH33 in the phloem response to aphids. Plant Cell Environ. 2007;30(2):187–201.

    CAS  Article  PubMed  Google Scholar 

  63. 63.

    Albert M, Werner M, Proksch P, Fry S, Kaldenhoff R. The cell wall-modifying xyloglucan endotransglycosylase/hydrolase LeXTH1 is expressed during the defence reaction of tomato against the plant parasite Cuscuta reflexa. Plant Biol. 2004;6(4):402–7.

    CAS  Article  PubMed  Google Scholar 

  64. 64.

    Maldonado-Mendoza IE, Dewbre GR, Blaylock L, Harrison MJ. Expression of a xyloglucan endotransglucosylase/hydrolase gene Mt-XTH1 from Medicago truncatula is induced systemically in mycorrhizal roots. Gene. 2005;345(2):191–7.

    CAS  Article  PubMed  Google Scholar 

  65. 65.

    Dhugga KS, Barreiro R, Whitten B, Stecca K, Hazebroek J, Randhawa GS, Dolan M, Kinney AJ, Tomes D, Nichols S, et al. Guar seed beta-mannan synthase is a member of the cellulose synthase super gene family. Science. 2004;303(5656):363–6.

    CAS  Article  PubMed  Google Scholar 

  66. 66.

    Liepman AH, Nairn CJ, Willats WG, Sørensen I, Roberts AW, Keegstra K. Functional genomic analysis supports conservation of function among cellulose synthase-like A gene family members and suggests diverse roles of mannans in plants. Plant Physiol. 2007;143(4):1881–93.

    CAS  PubMed Central  Article  PubMed  Google Scholar 

  67. 67.

    Goubet F, Misrahi A, Park SK, Zhang Z, Twell D, Dupree P. AtCSLA7, a cellulose synthase-like putative glycosyltransferase, is important for pollen tube growth and embryogenesis in Arabidopsis. Plant Physiol. 2003;131(2):547–57.

    CAS  PubMed Central  Article  PubMed  Google Scholar 

  68. 68.

    Doblin MS, Pettolino FA, Wilson SM, Campbell R, Burton RA, Fincher GB, Newbigin E, Bacic A. A barley cellulose synthase-like CSLH gene mediates (1, 3; 1, 4)-β-D-glucan synthesis in transgenic Arabidopsis. Proc Natl Acad Sci. 2009;106(14):5996–6001.

    CAS  PubMed Central  Article  PubMed  Google Scholar 

  69. 69.

    Burton RA, Wilson SM, Hrmova M, Harvey AJ, Shirley NJ, Medhurst A, Stone BA, Newbigin EJ, Bacic A, Fincher GB. Cellulose synthase-like CslF genes mediate the synthesis of cell wall (1, 3; 1, 4)-ß-D-glucans. Science. 2006;311(5769):1940–2.

    CAS  Article  PubMed  Google Scholar 

  70. 70.

    Roudier F, Schindelman G, DeSalle R, Benfey PN. The COBRA family of putative GPI-anchored proteins in Arabidopsis. A new fellowship in expansion. Plant Physiol. 2002;130(2):538–48.

    CAS  PubMed Central  Article  PubMed  Google Scholar 

  71. 71.

    Schindelman G, Morikami A, Jung J, Baskin TI, Carpita NC, Derbyshire P, McCann MC, Benfey PN. COBRA encodes a putative GPI-anchored protein, which is polarly localized and necessary for oriented cell expansion in Arabidopsis. Genes Dev. 2001;15(9):1115–27.

    CAS  PubMed Central  Article  PubMed  Google Scholar 

  72. 72.

    Sedbrook JC, Carroll KL, Hung KF, Masson PH, Somerville CR. The Arabidopsis SKU5 gene encodes an extracellular glycosyl phosphatidylinositol–anchored glycoprotein involved in directional root growth. Plant Cell Online. 2002;14(7):1635–48.

    CAS  Article  Google Scholar 

  73. 73.

    Li Y, Qian Q, Zhou Y, Yan M, Sun L, Zhang M, Fu Z, Wang Y, Han B, Pang X. BRITTLE CULM1, which encodes a COBRA-like protein, affects the mechanical properties of rice plants. Plant Cell Online. 2003;15(9):2020–31.

    CAS  Article  Google Scholar 

  74. 74.

    Read S, Northcote D. Chemical and immunological similarities between the phloem proteins of three genera of the Cucurbitaceae. Planta. 1983;158(2):119–27.

    CAS  Article  PubMed  Google Scholar 

  75. 75.

    Owens R, Blackburn M, Ding B. Possible involvement of the phloem lectin in long-distance viroid movement. Mol Plant-Microbe Interact. 2001;14(7):905–9.

    CAS  Article  PubMed  Google Scholar 

  76. 76.

    Gómez G, Pallás V. A long-distance translocatable phloem protein from cucumber forms a ribonucleoprotein complex in vivo with Hop stunt viroid RNA. J Virol. 2004;78(18):10104–10.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  77. 77.

    McNairn R, Currier H. Translocation blockage by sieve plate callose. Planta. 1968;82(4):369–80.

    CAS  Article  PubMed  Google Scholar 

  78. 78.

    Hong Z, Delauney AJ, Verma DPS. A cell plate–specific callose synthase and its interaction with phragmoplastin. Plant Cell Online. 2001;13(4):755–68.

    CAS  Article  Google Scholar 

  79. 79.

    Iglesias VA, Meins F. Movement of plant viruses is delayed in a β-1, 3-glucanase-deficient mutant showing a reduced plasmodesmatal size exclusion limit and enhanced callose deposition. Plant J. 2000;21(2):157–66.

    CAS  Article  PubMed  Google Scholar 

  80. 80.

    Musetti R, Paolacci A, Ciaffi M, Tanzarella OA, Polizzotto R, Tubaro F, Mizzau M, Ermacora P, Badiani M, Osler R. Phloem cytochemical modification and gene expression following the recovery of apple plants from apple proliferation disease. Phytopathology. 2010;100(4):390–9.

    CAS  Article  PubMed  Google Scholar 

  81. 81.

    Evans DE, Briars S-A, Williams LE. Active calcium transport by plant cell membranes. J Exp Bot. 1991;42(3):285–303.

    CAS  Article  Google Scholar 

  82. 82.

    White PJ. Calcium signals in root cells: the roles of plasma membrane calcium channels. Biologia. 2004;59(S13):77–83.

    CAS  Google Scholar 

  83. 83.

    Ali R, Zielinski RE, Berkowitz GA. Expression of plant cyclic nucleotide-gated cation channels in yeast. J Exp Bot. 2006;57(1):125–38.

    CAS  Article  PubMed  Google Scholar 

  84. 84.

    Ma W, Yoshioka K, Berkowitz GA. Cyclic nucleotide gated channels and Ca2+-mediated signal transduction during plant innate immune response to pathogens. Plant Signal Behav. 2007;2:548–50.

    PubMed Central  Article  PubMed  Google Scholar 

  85. 85.

    Balagué C, Lin B, Alcon C, Flottes G, Malmström S, Köhler C, Neuhaus G, Pelletier G, Gaymard F, Roby D . HLM1, an essential signaling component in the hypersensitive response, is a member of the cyclic nucleotide–gated channel ion channel family. Plant Cell Online. 2003;15(2):365–79.

    Article  CAS  Google Scholar 

  86. 86.

    Roberts DM, Harmon AC. Calcium-modulated proteins: targets of intracellular calcium signals in higher plants. Annu Rev Plant Biol. 1992;43(1):375–414.

    CAS  Article  Google Scholar 

  87. 87.

    Do Heo W, Lee SH, Kim MC, Kim JC, Chung WS, Chun HJ, Lee KJ, Park CY, Park HC, Choi JY. Involvement of specific calmodulin isoforms in salicylic acid-independent activation of plant disease resistance responses. Proc Natl Acad Sci. 1999;96(2):766–71.

    Article  Google Scholar 

  88. 88.

    Ishimaru Y, Suzuki M, Kobayashi T, Takahashi M, Nakanishi H, Mori S, Nishizawa NK. OsZIP4, a novel zinc-regulated zinc transporter in rice. J Exp Bot. 2005;56(422):3207–14.

    CAS  Article  PubMed  Google Scholar 

  89. 89.

    Grotz N, Fox T, Connolly E, Park W, Guerinot ML, Eide D. Identification of a family of zinc transporter genes from Arabidopsis that respond to zinc deficiency. Proc Natl Acad Sci. 1998;95(12):7220–4.

    CAS  PubMed Central  Article  PubMed  Google Scholar 

  90. 90.

    Vahling-Armstrong CM, Zhou H, Benyon L, Morgan JK, Duan Y. Two plant bacteria, S. meliloti and Ca. Liberibacter asiaticus, share functional znuABC homologues that encode for a high affinity zinc uptake system. PLoS ONE. 2012;7(5):e37340.

    CAS  PubMed Central  Article  PubMed  Google Scholar 

  91. 91.

    Hänsch R, Mendel RR. Physiological functions of mineral micronutrients (Cu, Zn, Mn, Fe, Ni, Mo, B, Cl). Curr Opin Plant Biol. 2009;12(3):259–66.

    Article  CAS  PubMed  Google Scholar 

  92. 92.

    Hernick M, Fierke CA. Zinc hydrolases: the mechanisms of zinc-dependent deacetylases. Arch Biochem Biophys. 2005;433(1):71–84.

    CAS  Article  PubMed  Google Scholar 

  93. 93.

    Merhej V, Georgiades K, Raoult D. Postgenomic analysis of bacterial pathogens repertoire reveals genome reduction rather than virulence factors. Brief Funct Genomics. 2013;12:291–304.

  94. 94.

    Martinelli F, Reagan RL, Uratsu SL, Phu ML, Albrecht U, Zhao W, Davis CE, Bowman KD, Dandekar AM. Gene regulatory networks elucidating huanglongbing disease mechanisms. PLoS ONE. 2013;8(9), e74256.

    CAS  PubMed Central  Article  PubMed  Google Scholar 

  95. 95.

    Seo HS, Song JT, Cheong J-J, Lee Y-H, Lee Y-W, Hwang I, Lee JS, Do Choi Y. Jasmonic acid carboxyl methyltransferase: a key enzyme for jasmonate-regulated plant responses. Proc Natl Acad Sci. 2001;98(8):4788–93.

    CAS  PubMed Central  Article  PubMed  Google Scholar 

  96. 96.

    Rodrigues CM, de Souza AA, Takita MA, Kishi LT, Machado MA. RNA-Seq analysis of Citrus reticulata in the early stages of Xylella fastidiosa infection reveals auxin-related genes as a defense response. BMC Genomics. 2013;14(1):676.

    CAS  PubMed Central  Article  PubMed  Google Scholar 

  97. 97.

    Dong J, Chen C, Chen Z. Expression profiles of the Arabidopsis WRKY gene superfamily during plant defense response. Plant Mol Biol. 2003;51(1):21–37.

    CAS  Article  PubMed  Google Scholar 

  98. 98.

    Li J, Brader G, Palva ET. The WRKY70 transcription factor: a node of convergence for jasmonate-mediated and salicylate-mediated signals in plant defense. Plant Cell Online. 2004;16(2):319–31.

    CAS  Article  Google Scholar 

  99. 99.

    Qiu D, Xiao J, Xie W, Liu H, Li X, Xiong L, Wang S. Rice gene network inferred from expression profiling of plants overexpressing OsWRKY13, a positive regulator of disease resistance. Mol Plant. 2008;1(3):538–51.

    CAS  PubMed  Google Scholar 

  100. 100.

    Bennett RN, Wallsgrove RM. Secondary metabolites in plant defence mechanisms. New Phytol. 1994;127(4):617–33.

    CAS  Article  Google Scholar 

  101. 101.

    Van Damme M, Huibers RP, Elberse J, Van den Ackerveken G. Arabidopsis DMR6 encodes a putative 2OG-Fe (II) oxygenase that is defense-associated but required for susceptibility to downy mildew. Plant J. 2008;54(5):785–93.

    Article  CAS  PubMed  Google Scholar 

  102. 102.

    Croteau R, Kutchan TM, Lewis NG. Natural products (secondary metabolites). Biochem Mol Biol Plants. 2000;24:1250–319.

    Google Scholar 

  103. 103.

    Harborne JB. Introduction to ecological biochemistry 4th edition. Academic Press; 1994. ISBN 9780080918594.

  104. 104.

    Chen K, Fan B, Du L, Chen Z. Activation of hypersensitive cell death by pathogen-induced receptor-like protein kinases from Arabidopsis. Plant Mol Biol. 2004;56(2):271–83.

    CAS  Article  PubMed  Google Scholar 

  105. 105.

    Lange J, Xie Z-P, Broughton WJ, Vögeli-Lange R, Boller T. A gene encoding a receptor-like protein kinase in the roots of common bean is differentially regulated in response to pathogens, symbionts and nodulation factors. Plant Sci. 1999;142(2):133–45.

    CAS  Article  Google Scholar 

  106. 106.

    Bradshaw Jr HD, Hollick JB, Parsons TJ, Clarke HR, Gordon MP. Systemically wound-responsive genes in poplar trees encode proteins similar to sweet potato sporamins and legume Kunitz trypsin inhibitors. Plant Mol Biol. 1990;14(1):51–9.

    CAS  Article  PubMed  Google Scholar 

  107. 107.

    Wüthrich KL, Bovet L, Hunziker PE, Donnison IS, Hörtensteiner S. Molecular cloning, functional expression and characterisation of RCC reductase involved in chlorophyll catabolism. Plant J. 2000;21(2):189–98.

    Article  PubMed  Google Scholar 

  108. 108.

    Greenberg JT, Guo A, Klessig DF, Ausubel FM. Programmed cell death in plants: a pathogen-triggered response activated coordinately with multiple defense functions. Cell. 1994;77(4):551–63.

    CAS  Article  PubMed  Google Scholar 

  109. 109.

    Mach JM, Castillo AR, Hoogstraten R, Greenberg JT. The Arabidopsis-accelerated cell death gene ACD2 encodes red chlorophyll catabolite reductase and suppresses the spread of disease symptoms. Proc Natl Acad Sci. 2001;98(2):771–6.

    CAS  PubMed Central  Article  PubMed  Google Scholar 

  110. 110.

    Hartung JS, Paul C, Achor D, Brlansky R. Colonization of Dodder, Cuscuta indecora, by ‘Candidatus Liberibacter asiaticus’ and ‘Ca. L americanus’. Phytopathology. 2010;100(8):756–62.

    Article  PubMed  Google Scholar 

  111. 111.

    Li W, Hartung JS, Levy L. Quantitative real-time PCR for detection and identification of ‘Candidatus Liberibacter’ species associated with citrus huanglongbing. J Microbiol Methods. 2006;66(1):104–15.

    CAS  Article  PubMed  Google Scholar 

  112. 112.

    Roy A, Ananthakrishnan G, Hartung JS, Brlansky R. Development and application of a multiplex reverse-transcription polymerase chain reaction assay for screening a global collection of Citrus tristeza virus isolates. Phytopathology. 2010;100(10):1077–88.

    CAS  Article  PubMed  Google Scholar 

  113. 113.

    Schroeder A, Mueller O, Stocker S, Salowsky R, Leiber M, Gassmann M, Lightfoot S, Menzel W, Granzow M, Ragg T. The RIN: an RNA integrity number for assigning integrity values to RNA measurements. BMC Mol Biol. 2006;7(1):3.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  114. 114.

    Feng J, Meyer CA, Wang Q, Liu JS, Liu XS, Zhang Y. GFOLD: a generalized fold change for ranking differentially expressed genes from RNA-seq data. Bioinformatics. 2012;28(21):2782–8.

    CAS  Article  PubMed  Google Scholar 

  115. 115.

    Yi X, Du Z, Su Z. PlantGSEA: a gene set enrichment analysis toolkit for plant community. Nucleic Acids Res. 2013;41(W1):W98–103.

    PubMed Central  Article  PubMed  Google Scholar 

  116. 116.

    Thimm O, Bläsing O, Gibon Y, Nagel A, Meyer S, Krüger P, Selbig J, Müller LA, Rhee SY, Stitt M. mapman: a user‐driven tool to display genomics data sets onto diagrams of metabolic pathways and other biological processes. Plant J. 2004;37(6):914–39.

    CAS  Article  PubMed  Google Scholar 

  117. 117.

    Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2−ΔΔCT method. Methods. 2001;25(4):402–8.

    CAS  Article  PubMed  Google Scholar 

Download references


The research was supported by China Scholarship Council, Southwest University and United States Department of Agriculture-Agricultural Research Service. The authors thank Cristina Paul for excellent plant care without which the study would not have been possible.

Author information



Corresponding authors

Correspondence to Changyong Zhou or John S. Hartung.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

SMF participated in the design of the experiments prepared RNA for sequencing, analyzed the data and drafted the manuscript; JS participated in the design of the study and carried out all statistical analyses; JSH conceived of the study and participated in its design and coordination and helped to analyze the data and draft the manuscript; JSH and CYZ revised the manuscript. All authors read and approved the final manuscript.

Additional file

Additional file 1: Figure S1.

Summary of reads from Citrus sinensis infected with CTV-B2, CTV-B6 and CaLas-B232 and the mock-inoculated control (H). Reads were obtained through paired-end RNA sequencing and mapped to the C. sinensis genome. Figure S2. Differentially expressed transcripts in diseased compared to healthy C. sinensis. Transcripts were up- or down-regulated in response to infection by CTV-B2, CTV-B6 and CaLas-B232. Figure S3. Melt curves of gene products assayed by RT-qPCR demonstrate specificity and reproducibility of the qPCR reactions. Gene targets are described in Table S11. Figure S4. Summary of differentially expressed transcripts based on functional categorization by Mapman. A. cell wall modification and cell organization; B. transcription factors; C. transporters; D. hormone metabolism; E. secondary metabolism; F. signaling; G. disease resistance proteins; H. proteolysis; I. development. Figure S5. Leaves that were sampled for RNA-seq. A, healthy; B, CTV-B2 positive; C, CTV-B6 positive; D, CaLas-B232 positive. Table S1. Differentially expressed transcripts of genes with unknown functions in response to CTV-B2, CTV-B6 and CaLas-B232. Table S2. GO categorization of genes that were differentially expressed in response to CTV-B2, CTV-B6 and CaLas-B232. Table S3. GO categorization of genes that were differentially expressed in response to CTV-B2, CTV-B6 but not CaLas-B232. Table S4. GO categorization of genes that responded to CTV-B2 and CaLas-B232 but not CTV-B6. Table S5. GO categorization of genes that responded to CTV-B6 and CaLas-B232 but not CTV-B2. Table S6. GO categorization of genes that responded only to CTV-B2. Table S7. GO categorization of genes that responded only to CTV-B6. Table S8. GO categorization of genes that responded only to CaLas-B232. Table S9. Complete set of differentially expressed transcripts in response to the three pathogens. Table S10. Estimates of changes in gene expression (fold change) by RT-qPCR and Gfold. Table S11. Primers used for qRT-PCR validation of differentially expressed transcripts identified by Gfold. (ZIP 1765 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Fu, S., Shao, J., Zhou, C. et al. Transcriptome analysis of sweet orange trees infected with ‘Candidatus Liberibacter asiaticus’ and two strains of Citrus Tristeza Virus . BMC Genomics 17, 349 (2016).

Download citation


  • Huanglongbing
  • Gfold
  • Host-pathogen interactions