Skip to main content
  • Research article
  • Open access
  • Published:

Temporal and spatial detection of Candidatus Liberibacter asiaticus putative effector transcripts during interaction with Huanglongbing-susceptible, −tolerant, and -resistant citrus hosts



Citrus Huanglongbing (HLB) is a bacterial disease with high economic significance. The associated agent Candidatus Liberibacter asiaticus is a fastidious, phloem-limited, intracellular bacterium that is transmitted by an insect vector the Asian citrus psyllid (ACP). The genome of Ca. L. asiaticus contains protein secretion machinery that suggests host cell modulation capacity of this bacterium.


A total of 28 candidate effectors, an important class of secreted proteins, were predicted from the Ca. L. asiaticus genome. Sequence specific primers were designed for reverse transcription (RT) and quantitative PCR (qPCR), and expression was validated for 20 of the effector candidates in infected citrus with multiple genetic background. Using detached leaf inoculation, the mRNA of effectors was detected from 6 h to 7 days post ACP exposure. It was observed that higher bacterial titers were associated with a larger number of effectors showing amplification across all samples. The effectors’ expression were compared in citrus hosts with various levels of HLB tolerance, including susceptible Duncan grapefruit and Washington navel orange, tolerant citron and Cleopatra mandarin, and resistant Pomeroy trifoliate and Carrizo citrange. Across all genotypes relatively high expression was observed for CLIBASIA_03695, CLIBASIA_00460, CLIBASIA_00420, CLIBASIA_04580, CLIBASIA_05320, CLIBASIA_04425, CLIBASIA_00525 and CLIBASIA_05315 in either a host-specific or -nonspecific manners. The two genotypes in each HLB-response group also show effector-expression profiles that seem to be different. In a companion study, the expression of effectors was compared between leaves and roots of own-rooted citrus that had been Ca. L. asiaticus-infected for more than a year. Results indicated relatively high expression of CLIBASIA_03875, CLIBASIA_04800 and CLIBASIA_05640 in all leaf and some root tissues of citron, Duncan and Cleopatra.


This temporal and spatial expression analysis of Ca. L. asiaticus effectors identified candidates possibly critical for early bacterial colonization, host tolerance suppression and long-term survival which are all worthy of further investigation.


Huanglongbing (HLB) is currently the most devastating disease among citrus crops, significantly impacting the production in the United States and the industry worldwide. In the U.S. it is associated with the bacterium Candidatus Liberibacter asiaticus (Ca. L. asiaticus). Ca. L. asiaticus is an α-proteobacterium restricted in citrus to the phloem cells and is transmitted by the insect vector Asian citrus psyllid (ACP, Diaphorina citri). Likely associated with its intracellular nature, Ca. L. asiaticus is fastidious and it has not been maintained in sustainable culture in axenic conditions [1, 2]. Therefore, Koch’s Postulate remains incomplete for this pathogen and studies related to bacterial pathogenesis are especially challenging. However, complete genome sequences of several Ca. L. asiaticus isolates have been reported and made publicly available [3,4,5,6]. Sequence data has provided useful genetic information on the evolution of the bacterium and permitted inferences on virulence mechanisms and essential metabolism. The genomic dataset includes many putative genes that can be subject to reverse genetics for functional analysis.

Plants have evolved layered innate immune systems that function via inter- and intracellular defenses [7]. Cell surface receptors can perceive pathogen-derived molecules such as pathogen-associated molecular patterns (PAMPs) which then activate PAMP-triggered immunity (PTI). In turn, pathogens secrete effector proteins that suppress PTI or regulate plant physiology to facilitate disease development, which results in effector-triggered susceptibility (ETS). Incompatible hosts maintain resistance (R) genes capable of recognizing the effectors and then causing rapid immune responses categorized as effector-triggered immunity (ETI). In the Ca. L. asiaticus genome, a set of flagellum-associated genes have been identified [3] including the Fla gene which includes the 22-amino acid PAMP flg22 [8]. Assays with citrus hosts indicated Ca. L. asiaticus flg22 elicited plant defenses that differed in strength between susceptible and tolerant citrus suggesting this may have a role in disease tolerance [9]. In the study of Ca. L. asiaticus effectors, in silico genome searches uncovered a repertoire of candidates [10, 11]. Protein function studies revealed that the effector CLIBASIA_05315 induced ETI-like reactions in Nicotiana benthamiana [11] and contributed to excessive cellular starch accumulation, a typical physiological disorder associated with HLB [12]. Recently, a papain-like cysteine protease was identified to be the target of CLIBASIA_05315 in citrus and this uncovered an interesting aspect of virulence mechanism for HLB [13].

Effector biology is emerging as an important aspect of the investigation on plant-pathogen interactions, as secreted effector proteins play many roles in the pathogenicity that lead to initial infection and host colonization. An understanding of fundamental effector biology is key to revealing pathogen evolution to achieve virulence, including subcellular localizations, mechanisms of host cell modulation, and in planta binding targets including susceptible and R genes [14,15,16]. Importantly, characterization of effectors provides guidance on and accelerates crop resistance breeding. In combating potato late blight caused by Phytophthora infestans, increasing numbers of R gene-effector pairs have been identified [17, 18] and implemented in breeding to select for naturally-occurring R genes that provide durable resistance [19]. Alternatively, R genes identified from biochemical or genetic approaches can be cloned and used in transgenic production of resistant cultivars, a strategy that has been demonstrated in tomato [20], rice [21], potato [22] and alfalfa [23] as examples.

Genetic resistance to HLB appears to be lacking within the Citrus genus. However, disease tolerance as shown by mild symptom and less impairment on tree development have been observed in both controlled inoculations [24, 25] and field pathogen exposure [26,27,28]. In Florida, an evaluation of 83 Citrus and Citrus relatives accessions with 6-years of field exposure of HLB identified citron (Citrus medica) as a genetic source of tolerance [28]. The Citrus relatives trifoliate orange (Poncirus trifoliata) and its hybrids have been an important resource as rootstocks and showed marked HLB tolerance/resistance in multiple studies [27, 29, 30]. On the other hand, commercially important citrus types such as grapefruits (C. paradisi), sweet oranges (C. sinensis) and mandarins (C. reticulata) are all negatively impacted by HLB, but disease severity varies greatly. Several studies report potentially useful tolerance in some mandarins [24, 26, 31], which also appears to be heritable (Stover, unpublished results).

In this study, we conducted in silico analysis of the Ca. L. asiaticus genome and identified a total of 28 effector candidates. The mRNA of 20 effectors were detected in Ca. L. asiaticus-infected citrus by RT-qPCR. Using a detached leaf assay for insect-mediated bacterial transmission, effector mRNA could be detected from 6 h to 7 days after ACP infestation, and the number of effectors detected was positively correlated with the bacterial titer. Subsequently, the expression of effectors was compared in six citrus types with different HLB tolerance levels, including citron, Duncan grapefruit, Cleopatra mandarin, Pomeroy trifoliate, Washington navel orange and Carrizo citrange. Results indicated some effector candidates had relatively high expression in multiple citrus genotypes regardless of tolerance levels, while other effectors had host-specific expression patterns. The effectors expression was also compared between leaf and root tissues in own-rooted citrus which had been infected with Ca. L. asiaticus for more than a year. Several candidate effectors were found to have relatively high transcriptional level in leaf tissue and some had higher expression in roots. Taken together, these studies provide guiding information on Ca. L. asiaticus utilization/deployment of effectors at early- and late-colonization stages and different tissues and hosts, which may lead to the selection of promising candidates for functional analysis and direct citrus resistance breeding.


Ca. L. asiaticus Secretome analysis for effector candidate identification

The Ca. L. asiaticus psy62 v1 genome (Genbank accession NC_012985.3) [3] was used for bioinformatic prediction of effectors. Based on common features of bacterial secreted proteins, the 1136 protein-coding sequences were filtered using SignalP [32] to predict the presence of signal peptides (SP) (Additional file 1: Table S1A) which were then subjected to TMHMM v2 [33, 34] to eliminate proteins with transmembrane domains [35]. Subsequently, the candidates with fewer than 200 amino acids were selected, and homology was searched using the National Center for Biotechnology Information (NCBI) BLAST tool. The proteins without any annotated homologs were analyzed in this study. A group of 28 putative Ca. L. asiaticus effectors was identified and used for the following experiments (Additional file 1: Table S1B and Fig. 1).

Fig. 1
figure 1

The flow chart of the study on expression of Ca. L. asiaticus effectors in citrus hosts. (a) The bioinformatics pipeline was applied to identify secreted candidate effectors using 1136 CDS from Ca. L. asiaticus genome. (b) Preliminary detection of mRNA using RT-qPCR measured expression of 20 candidates in infected citrus. (c) Effector expression pattern analysis during early and late bacterial-host interactions, in different citrus genotypes and between different tissue types

Detection of putative effectors by RT-qPCR in the multiple citrus genotypes infected by Ca. L. asiaticus

Sequence specific primers (SSPs) were designed for both reverse transcription and qPCR for greater specificity and higher capacity in amplification of effectors at low enrichment. Primers for reverse transcription were chosen at the 3′ end of the mRNAs to direct the synthesis of the first strand of cDNA. Forward and reverse primers were designed downstream of the cDNA synthesis primer binding region (Additional file 2: Figure S1). As one of the objectives was to determine if the expression pattern of effectors differed between HLB -susceptible, −tolerant and -resistant citrus, the SSPs efficacy were tested using multiple citrus species where specificity was defined as positive amplification in the Ca. L. asiaticus-infected tissue but not when uninfected. Hence RNA was isolated from Ca. L. asiaticus positive and negative leaves of citron, Duncan grapefruit and Cleopatra mandarin, Washington navel, Pomeroy trifoliate and Carrizo citrange for cDNA synthesis and qPCR. The SSPs of 20 effector candidates showed Ca. L. asiaticus specific amplification of targets in all citrus genotypes (Fig. 1 and Additional file 3: Table S2), whereas the remaining amplified in at least one of the uninfected citrus genotypes.

Effector expression during early Ca. L. asiaticus-citrus interactions and bacterial titer dependent effector mRNA detection

Due to the fastidious nature of this organism, common bacterial inoculation methods cannot be used to study the early interactions between Ca. L. asiaticus and citrus. This study used ACP infestation with detached leaves as previously described for inoculation [36] (Fig. 2a) to study an early infection time course. To compare Ca. L. asiaticus effector expression in citrus genotypes with different HLB tolerance levels, citron, Duncan, and Cleopatra citrus were used. Leaves from each genotype were exposed to ACP feeding for 6 h (h), 1, 3 and 7 days (d). Midrib DNA analysis indicated some successful bacterial transmissions at all of the time points and citrus genotypes, which provided multiple Ca. L. asiaticus positive samples for effector expression study (Fig. 2b). Three randomly selected infected samples at each time point and genotype were used for RT-qPCR, and showed detection of the mRNA of effectors at all the time points and genotypes studied (Additional file 4: Table S3). The number of effectors detected by RT-qPCR was generally low in most of samples which made it difficult to apply mean calculations and statistics; although some effector showed expression across more than two time points such as ‘4025’ in citron and ‘5315’ in Duncan and Cleopatra (From here on the last 4 digits of Ca. L. asiaticus gene locus name are used for simplicity). Noticeably, two samples with higher bacterial titers at 6 h and 3d in citron provided detectable mRNA of larger number of effectors (Additional file 4: Table S3), which led to a hypothesis that bacterial titer and number of effectors detected by RT-qPCR are positively correlated.

Fig. 2
figure 2

Photographs of detached leaf inoculation using Ca. L. asiaticus-infected ACP (a) and the distribution of bacterial titer measured from leaf samples collected at 6 h, 1d, 3d and 7d post ACP exposure (b-d). The experiments were conducted using leaves from citron (b), Duncan (c) and Cleopatra (d). The bacterium was quantified by qPCR amplifying the 16 s rDNA with Las-long primer from 100 ng of citrus DNA template. The cycle threshold values (Ct) were used to generate the dot plots for each time point

To test this hypothesis, additional Ca. L. asiaticus-infected samples at a later infection stage were included for a larger dataset and more diverse sample pools. These RNA samples were generated from a previous study on citrus transcriptome response to Ca. L. asiaticus at 2, 4, 6, 8 and 10 weeks (wk) after inoculation [9], from which multiple samples from each time point were used for effector mRNA detections. From the RT-qPCR data, the number of effectors expressed was counted and plotted against the Ct values quantifying Ca. L. asiaticus 16 s in each sample and each genotype. The results indicated a significant linear relationship between the two groups of variables in the three citrus genotypes tested, where the liner model fit was highest in Cleopatra (R2 = 0.80) and lowest in Duncan (R2 = 0.33) (Fig. 3a-c).

Fig. 3
figure 3

A linear relationship between the numbers of effectors detected by RT-qPCR and Ca. L. asiaticus titer in citron (a), Duncan (b) and Cleopatra (c). The calculation used samples from the detached leaf assay collected at 6 h, 1, 3 and 7 d post ACP exposure, and from leaf samples collected at 2, 4, 6, 8 and 10 wk after ACP inoculation to the citrus plants (Shi et al., 2017). The data analysis was done by JMP Genomics Fit Y by X function

Differential effector expression in citrus species with various HLB tolerance

To determine if effector expression patterns are different between HLB tolerant, susceptible, and resistant hosts at the same infection stage, the detached leaf study used in two highly susceptible citrus genotypes (Duncan and Washington navel orange), two tolerant genotypes (Cleopatra and citron), and two resistant genotypes (Pomeroy trifoliate and Carrizo citrange). Leaves of the six citrus genotypes were exposed to a 7-day infestation by the same ACP population for bacterial inoculation and only leaves which were Ca. L. asiaticus positive were used for effector analysis. The relative quantification (RQ) of effectors was compared within each citrus type to identify the higher expressed candidates that were likely to be critical for virulence. In the susceptible Duncan, only ‘4425’ had high expression while the level of other detected effectors were similar (Fig. 4a). Effectors ‘0460’ and ‘3695’ were the most expressed effectors in Washington Navel (Fig. 4b). In citron ‘3695’ showed the highest expression among all the effectors, and ‘0460’ and ‘0420’ were also highly ranked (Fig. 4c). Similarly ‘3695’ was the highest expressed effector in Cleopatra, and ‘4580’ and ‘5320’ had higher expression than most of the others (Fig. 4d). The trifoliate orange Pomeroy showed ‘0420’ and ‘0460’ as the highest expressing group, followed by ‘0525’ and ‘5315’ as the next group (Fig. 4e). In Carrizo, ‘0460’ and ‘4580’ showed significant higher levels than some other effectors, while the rest had no difference from each other (Fig. 4f). Therefore, we noted genotype-specific high expression of ‘4425’ in the susceptible Duncan, and genotype non-specific expression of ‘3695’ and ‘0460’ in multiple citrus. The two genotypes in each HLB-response group also show effector-expression profiles that seem to be different from each other.

Fig. 4
figure 4

The expression of Ca. L. asiaticus effectors in Duncan (a), Washington navel (b), citron (c), Cleopatra (d), Pomeroy (e), and Carrizo (f) after 7-day ACP infestation in a detached leaf assay. Three randomly selected Ca. L. asiaticus positive leaves were used for RNA isolation and expression analysis. The Ct values of each effector generated by RT-qPCR were normalized with citrus endogenous control UPL7 and transformed into relative quantification (RQ) values (2-ΔCt). The effectors with expression detections were subject to pair-wise comparison of standard least-square means (LS means) with Student’s t-test (p < 0.05). Expression levels indicated with the same letter are not significantly different. Bars are means ± standard error (n = 3)

Differential effector expression between leaf and root tissues in Ca. L. asiaticus-infected citrus

It is well known that infection of Ca. L. asiaticus occurs in citrus roots, which causes damage that directly affects tree health and also serves as source of inoculum to new foliar flushes through vascular movement of the bacterium [37]. Comparative analysis of the effector expression in leaf and root was conducted to study if effector expression has a tissue-specific pattern suggesting the bacterium employs different pathogenic strategies to colonize above- and below-ground plant organs. Ca. L. asiaticus-infected own-rooted citron, Duncan, and Cleopatra were used that had inoculated by ACP infestation in 2016 [9]. All plants were tested Ca. L. asiaticus positive in 2017 and after one year tolerant citron and Cleopatra plants were symptomatic but still maintained similar growth to their uninfected controls (Fig. 5a and c), whereas Duncan plants showed stunted growth and branch dieback (Fig. 5b).

Fig. 5
figure 5

Comparison of effector expression in the leaf and root of citron (a), Duncan (b) and Cleopatra (c) infected by Ca. L. asiaticus. Photographs display the health of Ca. L. asiaticus-infected plants (a-c) over a year after inoculation (right), in comparison to their uninfected controls (left). The Ca. L. asiaticus was quantified by qPCR amplifying 16 s rDNA from 100 ng of DNA template. A standard curve method with the function Log (copy number) = − 0.289*Ct + 11.66 was used to calculate 16 s copy numbers (d). The RNA was isolated from the same group of samples for effector analysis by RT-qPCR, from which the number of effectors detected was compared between leaf and root tissues in citron, Duncan and Cleopatra (E). Analysis was based on three biological replicates and was analyzed by Student’s t-test (p < 0.05) and significant difference between leaf/root pairs is indicated by an asterisk

Bacterial titers determined by 16 s quantification were lower in the roots of Duncan and Cleopatra than leaves, although this difference was insignificant in citron (Fig. 5d). Correspondingly, RT-qPCR analysis of effectors indicated there were lower numbers of effectors detected in roots than in leaves in all of the citrus types (Fig. 5e), which is in agreement with the previous observation on titer dependent effector mRNA detection (Fig. 3). For the effectors with detected mRNA, in citron ‘5640’ was the highest expressed effectors in both leaf and root tissues, and the expression of ‘3875’ and ‘4800’ were also among the highest in leaves (Fig. 6a-b). In Duncan, ‘3875’ and ‘4800’ ranked as the highest expressed effectors, although the difference was not significant for ‘4800’ in roots compared with other effectors detected (Fig. 6c-d). The expression of Ca. L. asiaticus effector candidates were generally higher in the leaves of Cleopatra than in citron and Duncan, with top expressed candidates including ‘4800’, ‘3875’ and ‘5640’ in leaf and they were the only effectors with mRNA detected in the roots of Cleopatra (Fig. 6e-f).

Fig. 6
figure 6

Comparison of effector expression in the leaf and root tissues of citron (a and b), Duncan (c and d) and Cleopatra (e and f). Leaf and root samples were collected from three plants of each citrus genotype that were inoculated 14–17 months prior via ACP infestation. In the expression analysis, the Ct values of each effector generated by RT-qPCR were normalized with citrus endogenous control UPL7 and transformed into relative quantification (RQ) values (2-delatCt). The effectors with expression detected were subject to pair-wise comparison of standard least-square means (LS means) with Student’s t-test (p < 0.05). Expression levels indicated with the same letter are not significantly different. Bars are means ± standard error (n = 3)


The secretome of a pathogenic bacterium represents an array of molecules that play offensive roles during colonization, among which effectors are an important class of proteins capable of suppressing defense and/or manipulating host physiology. Whole genome sequencing of Ca. L. asiaticus indicated that this intracellular fastidious bacterium lacks Type III and IV secretion systems typical for extracellular pathogens, but maintains genes encoding the general secretory pathway/Sec-translocon [3], which was proposed to be the major pathway for Ca. L. asiaticus effector secretion [10]. Due to common features such as presence of secretion SP, effectors can be predicted computationally [38] and this accelerates the discovery of key virulence factors and targeted host resistance/susceptibility genes [39]. In this study, a total of 28 effector candidates were identified, by filtering Ca. L. asiaticus genome for presence of SP, absence of transmembrane domain and relatively small protein size (Additional file 1: Table S1). This list of candidates contains the effectors identified during our previous screening with different parameters [11] and was also discovered in an analysis for Sec-translocon-dependent proteins [10], confirming the commonality of bacterial effector features and reliability of the bioinformatics methods utilized in this study.

One important aspect of pathogen effectors study is to analyze the expression levels during the initial microbe-plant interactions. Effector candidates with high levels of transcript suggests an active protein utilization by the pathogen during the infection processes, which may provide guidance on the selection of effectors for further functional studies. To determine Ca. L. asiaticus expression by RT-qPCR, we used SSPs for cDNA synthesis and different sets of SSPs for qPCR (Additional file 2: Figure S1), which resulted in higher specificity in target amplifications and better performance for low abundance targets than using random primers for cDNA, or the same reverse SSPs for cDNA synthesis and qPCR (data not shown). As a result, a total of 20 primer sets (Additional file 3: Table S2) provided good Ca. L. asiaticus-specific amplification in multiple citrus genotypes including the relatively distant species Poncirus trifoliata. The expression of Ca. L. asiaticus effectors have been reported in several studies. For example, Ca. L. asiaticus effector expression profile was compared between infected citrus and ACP, and revealed interesting candidates differentially expressed in the two hosts [40]. In another study, semi-quantitative RT-qPCR detected the expression of CLIBASIA_05315, CLIBASIA_00460, CLIBASIA_03230 and CLIBASIA_05640 in several citrus types, with CLIBASIA_05315 being promising as marker gene for disease early detection as it expressed in asymptomatic tissues [41]. However, because inoculation of bacteria was performed through grafting, plants materials in these studies were from relatively late infection stages which likely were at least several months following initial disease exposure [40]. Molecular events such as defense suppression by effectors and ETI-associated hypersensitive responses occur at early stages of contact, from hours to days after initial pathogen exposure [42]. Therefore we hypothesized that highly expressed Ca. L. asiaticus effectors during early host interactions are especially important for bacterial virulence, and their identification may lead to discovery of citrus resistance/susceptibility genes.

To test this hypothesis, we employed detached leaf inoculation to study effector expression within 7 days of initial bacterial contact (Fig. 2). Effective Ca. L. asiaticus inoculation was observed as soon as 6 h after infestation, resulting in a wide range of bacterial titers in multiple citrus genotypes (Fig. 2b-d). RT-qPCR analysis showed mRNA detection of only several effectors at various time points, which made it difficult to carry out mean comparisons for most of the genes. However, the expression of ‘4025’ in citron, and ‘5315’ in Duncan and Cleopatra were consistent at consecutive time points (Additional file 4: Table S3), suggesting possible roles during early host interactions, especially ‘5315’ which manipulates host cells for starch accumulation [12] and suppresses plant defenses [13]. Across all samples effector mRNA detection was titer-dependent, with higher number of effectors amplified in samples with higher bacterial quantification. This was further confirmed by analysis with greater sample sizes in citron, Duncan and Cleopatra (Fig. 3). This suggests that effector mRNA increases with bacterial cell number and thus effectors below the qPCR detection limit may become evident in citrus samples with high bacterial titer. It is also worth noting that the time points used in this study may not exactly correspond to ‘hours/days post inoculation (hpi/dpi)’ but rather to ‘inoculation access period (IAP)’ [36] within which inoculation occurs. This ongoing inoculation for seven days during our studies may explain ambiguous expression patterns of effectors in time point comparisons. Hence for simplicity we selected only 7-day IAP as representative of early bacterium-host interactions for the following studies.

To discover early and highly expressed effectors in citrus hosts with various HLB tolerance levels, six genotypes including citron, Duncan, Cleopatra, Pomeroy trifoliate, Washington navel, and Carrizo were subject to the detached leaf inoculation (7-day) by a common ACP population. Pair-wise comparisons ranked the effector expression levels to identify higher expressed candidates in each genotype (Fig. 4). Several effectors showed relative high expression in multiple genotypes, including ‘3695’ in citron, Cleopatra and Washington navel, ‘0460’ in citron, Washington navel, Pomeroy and Carrizo, ‘0420’ in citron and Pomeroy, and ‘4580’ in Cleopatra and Carrizo (Fig. 4a-f), suggesting these virulence factors are broadly active at the early infection stage. Among these effectors, the mRNA of ‘3695’ and ‘0460’ were detected at relatively high levels in both HLB-tolerant and susceptible citrus, suggesting they provide core virulence functions for bacterial colonization. Nevertheless, some effectors demonstrated tolerance/susceptibility-associated or host-specific high expression, including ‘0420’ in citron and Pomeroy, ‘4580’ in Cleopatra and Carrizo, ‘5320’ in Cleopatra, ‘4425’ in Duncan, ‘0525’ and ‘5315’ in Pomeroy, indicating host genetics may influence pathogen virulence factor expression. If confirmed, these tolerance- or resistance-associated effectors may be used as biomarkers to screen for HLB tolerance/resistance from citrus breeding materials. Further, the highly expressed effectors are good candidates for biochemical analysis to identify host binding targets, which may reveal important virulence mechanisms and lead to creation of resistant citrus.

A study was conducted using own-rooted citron, Duncan and Cleopatra plants that had been infected for more than a year, to evaluate effector expression in the pathogen after being fully established in tolerant and susceptible citrus hosts, assessing both leaf and root tissues (Fig. 5). Across all three citrus genotypes the expression patterns of effectors were similar between leaves but not roots (Fig. 6a-f). The number of effectors detected by RT-qPCR was less in roots than in leaves (Fig. 5d-e) which may be due to lower bacterial titers. Leaf expression of candidates ‘3875’, ‘4800’ and ‘5640’ were consistently among the highest regardless of citrus genotype/HLB tolerance levels (Fig. 6a, c and e). The expression profiles were markedly different from those of during the early infection for citron, Duncan and Cleopatra (Fig. 4a-c) which may suggest that Ca. L. asiaticus deploys different effectors over the time-course of infection. This is consistent with reports that this is a common strategy utilized by pathogens for successful colonization [15]. It is possible that our detached tissue-based study on early host interaction may be influenced by plant physiological changes associated with being detached from the mother plant. However, numerous researchers have employed the detached leaf inoculations to study outcomes of plant-pathogen interactions including accessing microbe pathogenicity and screening for resistant host and proved to produce consistent results to tests from whole plants [43,44,45,46,47]. In addition, at each sample time we also assessed mRNA for the citrus endogenous gene UPL7 and the expression was similar at all of time points and across all genotypes, suggesting metabolism was not compromised throughout the study period. In roots, high mRNA abundance of ‘5640’ in citron and ‘3875’ in Duncan and Cleopatra, even at low bacterial titers, may be evidence that their functions may contribute to HLB-related root damage such as root starch depletion and dieback reported to be causal rather than the consequence of HLB disease [37, 48]. Taken together, this part of the study identified effector candidates showing higher expression following established colonization in both HLB tolerant and susceptible citrus. The functional analysis of these candidates may enhance understanding of bacterial virulence and host interactions during chronic Ca. L. asiaticus infection.


A group of 28 candidate effectors were identified from the Ca. L. asiaticus genome via bioinformatics. Their transcriptional levels were studied in the infected citrus leaves exposed to ACP within 7 days from susceptible Duncan grapefruit and Washington navel orange, tolerant citron and Cleopatra mandarin, and resistant Pomeroy trifoliate and Carrizo citrange. Candidate effectors CLIBASIA_03695, CLIBASIA_00460, CLIBASIA_00420, CLIBASIA_04580, CLIBASIA_05320, CLIBASIA_04425, CLIBASIA_00525 and CLIBASIA_05315 showed relatively high expression in host-specific or -nonspecific manners. In citrus with HLB exposure for over one year, the expression of effectors was compared between leaves and roots tissues and indicated relatively high expression of CLIBASIA_03875, CLIBASIA_04800 and CLIBASIA_05640 in all leaf and some root tissues of citron, Duncan and Cleopatra genotypes. The identified Ca. L. asiaticus effectors candidates with high temporal and spatial expression levels may have roles in early bacterial colonization, host tolerance suppression and long-term survival that will be subject to further functional studies.


Plant materials

All of the plant materials were clonally produced and maintained at the US Horticultural Research Lab greenhouse facilities. Healthy seedlings of Duncan grapefruit, Cleopatra mandarin, citron, Pomeroy trifoliate, Washington naval orange and Carrizo citrange were grown in pots in the greenhouse that were watered daily and fertilized weekly. Horticultural mineral oil was applied as needed to control pests but no other insecticides were applied. For detached leaf inoculation, leaves for each citrus type were collected at similar maturity and size, rinsed and used immediately for described experiments. The leaf and root analysis used citron, Duncan and Cleopatra plants derived from seed and inoculated more than a year earlier, verified as Ca. L. asiaticus-infected, and maintained in the quarantine greenhouse.

Detached leaf inoculation and sample collection

The leaf samples were collected between 10 am to noon, gently rinsed and dried in the lab, with the petioles placed in water in 1.5 mL centrifuge tubes sealed with Parafilm. Each detached leaf, with its petiole in water, was placed in a 50 mL plastic centrifuge tube with a mesh top and was exposed to 10 psyllids collected from colonies maintained on Ca. L. asiaticus infected plants. Tubes containing leaves and psyllids were kept on the lab bench with supplemental light of 16 h per day. The light set up is previously described with light intensity at 9687 ± 3460 lm/m2 [49]. At the time of destructive sampling, the ACP were removed by vacuum and the leaves were wiped clean before sampling. A group of 10–15 leaves were collected each at 6 h, 1d, 3d and 7d post ACP infestation, or for each citrus genotype at the 7d time point. The midrib tissues were excised for PCR quantification of Ca. L. asiaticus 16 s rDNA, and the remaining leaf tissues were placed initially in liquid nitrogen and then stored at − 80 °C for RNA isolation and expression analysis.

Whole plant inoculation and sample collection

Ca. L. asiaticus-infected citron, Duncan and Cleopatra plants, all own-rooted from seed, were inoculated by ACP infestation during a previous study in 2016 and maintained in a greenhouse along with uninfected controls [9]. Briefly, six citrus plants with young flushes were infested by 400 Ca. L. asiaticus positive ACP or negative ACP (as mock inoculation) in dome cages for 2 weeks. Then insects were removed by aspiration, and plants were treated with carbaryl insecticide and greenhouse growth condition were resumed. A pool of 4 to 6 leaves were collected from each plant and used for leaf analysis of effector expression. Fibrous roots were sampled peripherally from the root mass of each plant, rinsed and stored at − 80 °C for analysis.

Determination of Ca. L. asiaticus bacterial titer

To quantify Ca. L. asiaticus in citrus leaves, DNA was isolated from midrib of leaves or fibrous roots using DNeasy plant mini kit (Qiagen, Gaithersburg, MD, USA), and the quantity and quality were determined by a NanoDrop Spectrophotometer (Thermo Scientific, Wilmington, DE, USA). For each sample, 100 ng of DNA was used for qPCR of Ca. L. asiaticus 16 s rDNA using Las Long primers and SYBR chemistry (Thermo Fisher Scientific, Waltham, MA, USA) and an ABI7500 thermal cycler (Applied Biosystems, Foster City, CA, USA). The threshold cycle (Ct) values were used to calculate bacterial titer using the standard curve method [9].

Expression analysis of effector candidates by RT-qPCR

The total RNA was isolated from Ca. L. asiaticus positive leaf or root samples using TriZol reagent (Invitrogen, Carlsbad, CA, USA) following the manufacturer’s instructions. On-column DNase treatment and RNA purification were performed using RNeasy Plant Mini Kit (Qiagen). The quantity and quality of RNA were determined using a NanoDrop Spectrophotometer (Thermo Scientific). Subsequently, 1 μg of RNA was used for cDNA synthesis by QuantiTect Reverse Transcription Kit (Qiagen). A mix of SSPs (Additional file 3: Table S2) targeting Ca. L. asiaticus effectors, 16 s and a citrus endogenous gene UPL7 was used for cDNA synthesis at a final concentration of 0.5 μM each. A 20 μL of reverse transcription reaction for each sample was performed according to the kit protocol. Subsequently, the cDNA was diluted to 5 ng/μL and 2 μL was used in each qPCR reaction together with forward and reverse SSPs (Additional file 3: Table S2) to amplify effector candidates.



Asian citrus psyllid

Ca. L. asiaticus:

Candidatus Liberibacter asiaticus

Ct :

Cycle threshold


Effector-triggered immunity


Effector-triggered susceptibility




Hours/days post inoculation


Inoculation access period


National Center for Biotechnology Information


Pathogen-associated molecular patterns


Polymerase chain reaction


Reverse transcription quantitative PCR


Signal peptide


Sequence specific primers


  1. Sechler A, et al. Cultivation of 'Candidatus Liberibacter asiaticus', 'Ca. L. africanus', and 'Ca. L. americanus' associated with Huanglongbing. Phytopathology. 2009;99(5):480–6.

    Article  CAS  PubMed  Google Scholar 

  2. Parker JK, et al. Viability of 'Candidatus Liberibacter asiaticus' prolonged by addition of Citrus juice to culture medium. Phytopathology. 2014;104(1):15–26.

    Article  CAS  PubMed  Google Scholar 

  3. Duan YP, et al. Complete genome sequence of Citrus Huanglongbing bacterium, 'Candidatus Liberibacter asiaticus' obtained through metagenomics. Mol Plant-Microbe Interact. 2009;22(8):1011–20.

    Article  CAS  PubMed  Google Scholar 

  4. Zheng Z, Deng X, Chen J. Whole-genome sequence of "Candidatus Liberibacter asiaticus" from Guangdong, China. Genome Announc. 2014;2(2).

  5. Zheng Z, et al. Whole-genome sequence of "Candidatus Liberibacter asiaticus" from a Huanglongbing-affected Citrus tree in Central Florida. Genome Announc. 2015;3(2).

  6. Kunta M, et al. Draft whole-genome sequence of "Candidatus Liberibacter asiaticus" strain TX2351 isolated from Asian Citrus psyllids in Texas, USA. Genome Announc. 2017;5(15).

  7. Jones JD, Dangl JL. The plant immune system. Nature. 2006;444(7117):323–9.

    Article  CAS  PubMed  Google Scholar 

  8. Zou HS, et al. The destructive Citrus pathogen, 'Candidatus Liberibacter asiaticus' encodes a functional Flagellin characteristic of a pathogen-associated molecular pattern. PLoS One. 2012;7(9).

  9. Shi QC, et al. Huanglongbing tolerance-associated genes are identified by comparative transcriptomics using bacterial flagellin 22 as a proxy to challenge citrus. Journal of Citrus Pathology. 2017:38.

  10. Prasad S, et al. SEC-Translocon dependent Extracytoplasmic proteins of Candidatus Liberibacter asiaticus. Front Microbiol. 2016;7:1989.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Pitino M, et al. Transient expression of Candidatus Liberibacter Asiaticus effector induces cell death in Nicotiana benthamiana. Front Plant Sci. 2016;7.

  12. Pitino M, Allen V, Duan Y. LasDelta5315 effector induces extreme starch accumulation and chlorosis as ca. Liberibacter asiaticus infection in Nicotiana benthamiana. Front Plant Sci. 2018;9:113.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Clark K, et al. An effector from the Huanglongbing-associated pathogen targets citrus proteases. Nat Commun. 2018;9(1):1718.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Hogenhout SA, et al. Emerging concepts in effector biology of plant-associated organisms. Mol Plant-Microbe Interact. 2009;22(2):115–22.

    Article  CAS  PubMed  Google Scholar 

  15. Win J, et al. Effector biology of plant-associated organisms: concepts and perspectives. Cold Spring Harb Symp Quant Biol. 2012;77:235–47.

    Article  CAS  PubMed  Google Scholar 

  16. Kourelis J, van der Hoorn RAL. Defended to the nines: 25 years of resistance gene cloning identifies nine mechanisms for R protein function. Plant Cell. 2018;30(2):285–99.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Hein I, et al. The zig-zag-zig in oomycete-plant interactions. Mol Plant Pathol. 2009;10(4):547–62.

    Article  CAS  PubMed  Google Scholar 

  18. Vleeshouwers VG, et al. Understanding and exploiting late blight resistance in the age of effectors. Annu Rev Phytopathol. 2011;49:507–31.

    Article  CAS  PubMed  Google Scholar 

  19. Vleeshouwers VGAA, Oliver RP. Effectors as tools in disease resistance breeding against biotrophic, Hemibiotrophic, and necrotrophic plant pathogens. Mol Plant-Microbe Interact. 2014;27(3):196–206.

    Article  CAS  PubMed  Google Scholar 

  20. Tai TH, et al. Expression of the Bs2 pepper gene confers resistance to bacterial spot disease in tomato. Proc Natl Acad Sci U S A. 1999;96(24):14153–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Zhao B, et al. A maize resistance gene functions against bacterial streak disease in rice. Proc Natl Acad Sci U S A. 2005;102(43):15383–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. van der Vossen EA, et al. The Rpi-blb2 gene from Solanum bulbocastanum is an Mi-1 gene homolog conferring broad-spectrum late blight resistance in potato. Plant J. 2005;44(2):208–22.

    Article  PubMed  Google Scholar 

  23. Yang S, et al. Alfalfa benefits from Medicago truncatula: the RCT1 gene from M. Truncatula confers broad-spectrum resistance to anthracnose in alfalfa. Proc Natl Acad Sci U S A. 2008;105(34):12164–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Folimonova SY, et al. Examination of the responses of different genotypes of Citrus to Huanglongbing (Citrus greening) under different conditions. Phytopathology. 2009;99(12):1346–54.

    Article  PubMed  Google Scholar 

  25. Cheema SS, Kapur SP, Chohan JS. Evaluation of rough lemon strains and other rootstocks against greening-disease of Citrus. Sci Hortic. 1982;18(1):71–5.

    Article  Google Scholar 

  26. Stover E, et al. Conventional Citrus of some Scion/rootstock combinations show field tolerance under high Huanglongbing disease pressure. Hortscience. 2016;51(2):127–32.

    Article  Google Scholar 

  27. Ramadugu C, et al. Long-term field evaluation reveals Huanglongbing resistance in Citrus relatives. Plant Dis. 2016;100(9):1858–69.

    Article  PubMed  Google Scholar 

  28. Miles GP, et al. Apparent tolerance to Huanglongbing in Citrus and Citrus-related germplasm. Hortscience. 2017;52(1):31–9.

    Article  Google Scholar 

  29. Albrecht U, Bowman KD. Tolerance of the trifoliate Citrus hybrid US-897 (Citrus reticulata Blanco × Poncirus trifoliata L. Raf.) to Huanglongbing. HortScience. 2011;46(1):16–22.

    Article  CAS  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

  31. Killiny N, et al. Metabolically speaking: possible reasons behind the tolerance of 'Sugar Belle' mandarin hybrid to huanglongbing. Plant Physiol Biochem. 2017;116:36–47.

    Article  CAS  PubMed  Google Scholar 

  32. Nielsen H, et al. Identification of prokaryotic and eukaryotic signal peptides and prediction of their cleavage sites. Protein Eng Des Sel. 1997;10(1):1–6.

    Article  CAS  Google Scholar 

  33. Sonnhammer EL, von Heijne G, Krogh A. A hidden Markov model for predicting transmembrane helices in protein sequences. Proc Int Conf Intell Syst Mol Biol. 1998;6:175–82.

    CAS  PubMed  Google Scholar 

  34. Krogh A, et al. Predicting transmembrane protein topology with a hidden Markov model: application to complete genomes. J Mol Biol. 2001;305(3):567–80.

    Article  CAS  Google Scholar 

  35. Petersen TN, et al. SignalP 4.0: discriminating signal peptides from transmembrane regions. Nat Methods. 2011;8(10):785–6.

    Article  CAS  PubMed  Google Scholar 

  36. Ammar el D, Walter AJ, Hall DG. New excised-leaf assay method to test inoculativity of Asian citrus psyllid (Hemiptera: Psyllidae) with Candidatus Liberibacter asiaticus associated with citrus huanglongbing disease. J Econ Entomol. 2013;106(1):25–35.

    Article  Google Scholar 

  37. Johnson EG, et al. 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. 2014;63(2):290–8.

    Article  Google Scholar 

  38. McDermott JE, et al. Computational prediction of type III and IV secreted effectors in gram-negative bacteria. Infect Immun. 2011;79(1):23–32.

    Article  CAS  PubMed  Google Scholar 

  39. Varden FA, et al. Taking the stage: effectors in the spotlight. Curr Opin Plant Biol. 2017;38:25–33.

    Article  CAS  PubMed  Google Scholar 

  40. Yan Q, et al. Global gene expression changes in Candidatus Liberibacter asiaticus during the transmission in distinct hosts between plant and insect. Mol Plant Pathol. 2013;14(4):391–404.

    Article  CAS  PubMed  Google Scholar 

  41. Pagliaccia D, et al. A pathogen secreted protein as a detection marker for Citrus Huanglongbing. Front Microbiol. 2017;8:2041.

    Article  PubMed  PubMed Central  Google Scholar 

  42. Hann DR, Rathjen JP. Early events in the pathogenicity of Pseudomonas syringae on Nicotiana benthamiana. Plant J. 2007;49(4):607–18.

    Article  CAS  PubMed  Google Scholar 

  43. Browne RA, Cooke BM. Development and evaluation of an in vitro detached leaf assay forc pre-screening resistance to Fusarium head blight in wheat. Eur J Plant Pathol. 2004;110(1):91–102.

    Article  Google Scholar 

  44. Pettitt TR, et al. A simple detached leaf assay provides rapid and inexpensive determination of pathogenicity of Pythium isolates to ‘all year round’ (AYR) chrysanthemum roots. Plant Pathol. 2011;60(5):946–56.

    Article  Google Scholar 

  45. Vleeshouwers VGAA, et al. A laboratory assay for Phytophthora infestans resistance in various Solanum species reflects the field situation. Eur J Plant Pathol. 1999;105(3):241–50.

    Article  Google Scholar 

  46. Huang S, et al. An accurate in vitro assay for high-throughput disease testing of Phytophthora infestans in potato. Plant Dis. 2005;89(12):1263–7.

    Article  PubMed  Google Scholar 

  47. Arraiano LS, Brading PA, Brown JKM. A detached seedling leaf technique to study resistance to Mycosphaerella graminicola (anamorph Septoria tritici) in wheat. Plant Pathol. 2001;50(3):339–46.

    Article  Google Scholar 

  48. Etxeberria E, et al. Anatomical distribution of abnormally high levels of starch in HLB-affected Valencia orange trees. Physiol Mol Plant Pathol. 2009;74(1):76–83.

    Article  CAS  Google Scholar 

  49. Hall DG, Moulton KM. Transmission rates of 'Candidatus Liberibacter asiaticus' to greenhouse seedlings by laboratory colonies of Asian Citrus psyllid (Hemiptera: Liviidae). J Econ Entomol. 2018;111(6):2546–52.

    PubMed  Google Scholar 

Download references


We thank Ellen Cochrane for support on use of various citrus plant materials, and Kathy Moulton for maintaining and providing ACP populations for Ca. L. asiaticus inoculation. Our thanks also go to Chelsea Veith, Spencer Marshall and Shelby Durden for lab technical assistance, Yolanda Avila for general lab supports and Shuo Zhang for assistance in data statistical analysis.


Funding was provided by the Citrus Research and Development Foundation (CRDF), Inc., Lake Alfred, FL, U.S.A. The funding body guided the design of the study and conduction of experiments through the funded project objectives.

Availability of data and materials

All data generated or analyzed during this study are included in this published article and its supplementary information files.

Author information

Authors and Affiliations



QS, MP, LC and ES conceived the project idea. QS and MP designed the experiments. QS, MP and LC conducted the experiments. QS, MP, LC, SZ, and JK analyzed the data. RS and DG contributed to materials and analysis tools. QS, MP and ES wrote the paper. All authors have read and approved the manuscript.

Corresponding author

Correspondence to Ed Stover.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interest

The authors declare that they have no conflicts of interest. USDA is an equal opportunity provider and employer. Mention of trade names or commercial products in this article is solely for the purpose of providing specific information and does not imply recommendation or endorsement by the United States Department of Agriculture.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

Table S1. The list of putative effectors identified by filtering Ca. L. asiaticus genome using bioinformatics. Supplemental Table 1-A contains the list of Ca. L. asiaticus genes with a signal peptide predicted by SignalP. Supplemental Table 1-B contains the list of Ca. L. asiaticus genes having signal peptide but not transmembrane domain, with the protein size smaller than 200 amino acids. (XLSX 22 kb)

Additional file 2:

Figure S1. The design of primers for reverse transcription (RT) and qPCR amplifying Ca. L. asiaticus effector candidates. The illustration uses CLIBASIA_05315 sequence as an example, where the RT primer (yellow) was selected towards the 3′ end of the sense strand (reverse complement to the selected), and the forward (red) and reverse (blue, reverse complement to the selected) primers for qPCR were located downstream along the cDNA synthesis direction. (DOCX 19 kb)

Additional file 3:

Table S2. Primer sequences for reverse transcription and qPCR analysis of effector expression. The list of primers used in this study for reverse transcription and qPCR (XLSX 10 kb)

Additional file 4:

Table S3. A time course study on the expression of Ca. L. asiaticus effectors in citron, Duncan and Cleopatra using the detached leaf assay. Three randomly selected Ca. L. asiaticus-infected samples were collected for each time point (Fig. 2) and were analyzed for expression by RT-qPCR. The data was presented as relative expression levels, based on normalization with citrus UPL7 (delta Ct), and transformed into 2-ΔCt. The table was formatted with color scales using Excel where values increase from green to red colors. The two highlighted columns were the two samples with high bacterial titer and larger number of detected candidate effectors. (XLSX 15 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

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Shi, Q., Pitino, M., Zhang, S. et al. Temporal and spatial detection of Candidatus Liberibacter asiaticus putative effector transcripts during interaction with Huanglongbing-susceptible, −tolerant, and -resistant citrus hosts. BMC Plant Biol 19, 122 (2019).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: