Skip to main content

Repertoire of novel sequence signatures for the detection of Candidatus Liberibacter asiaticus by quantitative real-time PCR



Huanglongbing (HLB) or citrus greening is a devastating disease of citrus. The gram-negative bacterium Candidatus Liberibacter asiaticus (Las) belonging to the α-proteobacteria is responsible for HLB in North America as well as in Asia. Currently, there is no cure for this disease. Early detection and quarantine of Las-infected trees are important management strategies used to prevent HLB from invading HLB-free citrus producing regions. Quantitative real-time PCR (qRT-PCR) based molecular diagnostic assays have been routinely used in the detection and diagnosis of Las. The oligonucleotide primer pairs based on conserved genes or regions, which include 16S rDNA and the β-operon, have been widely employed in the detection of Las by qRT-PCR. The availability of whole genome sequence of Las now allows the design of primers beyond the conserved regions for the detection of Las explicitly.


We took a complimentary approach by systematically screening the genes in a genome-wide fashion, to identify the unique signatures that are only present in Las by an exhaustive sequence based similarity search against the nucleotide sequence database. Our search resulted in 34 probable unique signatures. Furthermore, by designing the primer pair specific to the identified signatures, we showed that most of our primer sets are able to detect Las from the infected plant and psyllid materials collected from the USA and China by qRT-PCR. Overall, 18 primer pairs of the 34 are found to be highly specific to Las with no cross reactivity to the closely related species Ca. L. americanus (Lam) and Ca. L. africanus (Laf).


We have designed qRT-PCR primers based on Las specific genes. Among them, 18 are suitable for the detection of Las from Las-infected plant and psyllid samples. The repertoire of primers that we have developed and characterized in this study enhanced the qRT-PCR based molecular diagnosis of HLB.


Huanglongbing (HLB) or citrus greening is the most devastating disease of citrus, threatening the citrus industry worldwide, and leading to massive reduction in fruit production as well as death of infected trees [1]. The causal agents of HLB are three closely related gram-negative, phloem-limited α-proteobacteria Candidatus Liberibacter species [2, 3]. The heat tolerant strain Ca. L. asiaticus (Las) is the most widespread in Asia as well as in the USA whereas Ca. L. americanus (Lam) is mostly limited to South America [24]. Ca. L. africanus (Laf) is heat sensitive and localized to the African continent. All the three Liberibacter species are currently uncultured and are known to reside in the sieve tubes of the plant phloem [5] or in the gut of the phloem-feeding psyllids [6]. Psyllids are the natural vectors in transmitting the bacteria between plants [1, 6]. The Asian psyllid, Diaphorina citri Kuwayama (Homoptera: Psyllidae) is responsible for transmitting Las and Lam in Asia and America, while the African citrus psyllid, Trioza erytreae Del Guercio (Homoptera: Psyllidae), is the natural vector of Laf in Africa [7]. The characteristic symptoms of the infected plants include the yellow shoots, foliar blotchy mottles, along with poor flowering and stunting [1]. HLB also results in poorly colored, unpleasant tasting, reduced size fruit that shows staining of vascular columella and seed abortion [1]. Generally the fruit may remain partially green, for this reason HLB is also called citrus greening [1]. Chronically infected trees are sparsely foliated and display extensive twig or limb die-back and eventually die within three to five years [1]. Moreover, the disorders induced in diseased plants vary with cultivar, tree maturity, time of infection, stages of disease and other abiotic or biotic agents that affect the tree [1]. HLB symptoms also share certain similarities to nutrient deficiency [1], citrus stubborn disease caused by Spiroplasma citri[8] and a HLB-like disease caused by a phytoplasma [9, 10]. Early diagnosis and differentiation of Las infections from those defects and agents mentioned above, is thus critical to reducing the spread and devastation of this disease locally and via international trade, as well as minimizing the economic impact of potential false positive diagnoses.

Importantly, HLB and the Asian citrus psyllid (D. citri) are expanding to new citrus production areas. Currently, Asian citrus psyllid has been found in Florida, Texas, California, Arizona, Hawaii, Louisiana, Georgia, and Alabama in the USA, as well as in parts of South and Central America, Mexico, and the Caribbean. Meanwhile, HLB has not only been identified in Florida, Louisiana, South Carolina, Louisiana, Georgia, Texas and California of the USA; it has also been discovered in Cuba, Belize, Jamaica, Mexico, and other countries in the Caribbean [11]. While HLB and D. citri have been found in different producing areas, the number of infected trees and the psyllid vector population vary dramatically among different regions. Thus, different strategies of management of HLB are recommended for different regions, according to the corresponding severity of HLB and occurrence of psyllid vectors.

Currently, no efficient management strategy is available to control HLB. For the recently Las-infected citrus producing areas such as California, prevention and eradication of HLB are the most efficient and cost-effective approaches. Additionally, Las infected trees are most often found to be asymptomatic during the early stage of infection. Thus, accurate early detection of Las in citrus plants and psyllids is critical for enacting containment measures in non-endemic citrus producing areas. For the citrus producing areas without HLB, such as the Mediterranean region, accurate detection is critical for the success of quarantine measures against Ca. Liberibacter.

Methods such as biological indexing using graft, dodder transmission [12], isothermal loop amplification (LAMP) [13], electron microscopy [1], DNA probes [14], enzyme-linked immunosorbent assays (ELISA) [15], conventional PCR [1622] and quantitative real-time PCR (qRT-PCR) [2226] are used for the diagnosis and confirmation of HLB. Although diagnostic tools like conventional PCR and LAMP showed good sensitivity, they were not consistent in detection of Las bacterium from infected plant and psyllid materials [6, 13, 25]. The current HLB diagnostic detection mainly employs qRT-PCR based methods due to their sensitive and quantitative nature. The initial qRT-PCR oligonucleotide primer sets for the detection of Las, targeted rplKAJL-rpoBC operon (β-operon: CQULA04f/r) [26], 16S ribosomal RNA gene (rDNA) (HLBasf/r) [23], EUB338f/EUB518r [27], ALF518f/ EUB518r [27] or species specific variable regions. EUB338f/EUB518r primers are universal to Eubacteria [27], while ALF518f/EUB518r primers identify α-proteobacteria universally [27] including Las, therefore not specific. Furthermore, the primers based on the conserved 16S and β-operon regions are popular but nevertheless have been shown to pose a potential specificity issue, as both false negatives and false positives have been reported [28]. Therefore, efforts have been directed towards developing effective qRT-PCR primers that target other non-conserved sequences. Recent studies made use of intragenic repeat regions of the prophage sequence for the detection of Las by qRT-PCR [25]. However, the intragenic repeat regions of the prophage sequence were also identified in Lam. Therefore, these primer pairs, hyvi/hyvii did not distinguish between Las and Lam, posing a specificity issue [25]. Consequently, primer pairs that specifically detect Las and make clear distinction among other phylogenetically closely related bacteria are essential.

Here we took a complimentary approach to identify the genes that are unique to Las by a bioinformatic analysis with the goal of expanding the arsenal of tools for Las detection. The advancement in the genome sequencing of Las [29] provides an opportunity to design primers based on species specific sequences for the detection of Las. We designed the oligonucleotide primer pairs specific to the identified unique genic signatures. We further validated their specificities and selectivity against closely related strains that demonstrated the application to Las-infected tissues and insect vectors by a qRT-PCR.

Results and discussion

Recently, the whole genome sequences of Las [29, 30] have been sequenced. This allows for systematic screening of unique Las genes in a genome-wide fashion. The availability of the genome sequences of the closely related species Lam [31], L. crescens (Lcr) [32] and Ca. L. solanacearum (Lso) [33], further effectively helps in identification of unique regions, by minimizing the cross-species reactions, thereby enhancing the diagnostic identification of Las in a more distinct manner.

Bioinformatic analysis

Several high-throughput applications have been developed recently to design diagnostic primers using the whole genome sequence information including KPATH, Insignia, TOFI, and TOPSI [3440]. Among them, KPATH, Insignia, and TOPSI have the potential to be used for design of real-time PCR primers for qRT-PCR based assays for Las, whereas TOFI is used to design signatures for microarray-based assays. These methods mentioned above can be basically categorized into alignment-free and alignment-based approaches. The alignment-free approach uses both coding and non-coding regions of the genome and is useful for the genomes with less accurate sequence information, but generally result in high false positive rates as it does not involve pre-screening of the selected genomic loci for their discriminatory ability [37]. The alignment-based approach involves pre-screening of the selected genomic loci for their discriminatory ability [34]. This approach does not consider the genome annotation of genic and non-genic information, but rather aligns bigger regions of the genome, hence prone to lose shorter discriminatory sequence regions. Additionally, discriminatory ability of the selected regions are screened bioinformatically only on limited number of closely related species, which provide more opportunities for false positives. We therefore took a complementary bioinformatics approach by pre-screening shorter genic regions against the nucleotide sequence database (nt) at NCBI, to identify all the possible unique genic regions from the Las genome. The natural selection acts more strongly on genic region, hence use of discriminatory sequences in this region results in less false positives as the organisms are under selection pressure [41]. Additionally, pre-screening against the nt is more effective as it contains the largest pool of well-annotated nucleotide sequences from different organisms. We envisioned that these two steps would result in more specific detection of target organism with less false positives, hence are included in our bioinformatics approach.

There are ~1100 genes assigned to the Las genome. Therefore, manual searching of each of these sequences against the nt database using BLAST program [42, 43] is a laborious and time consuming procedure. Hence, we automated this sequence similarity search step by developing a standalone PERL script (Additional file 1). This script performed the similarity searches for each of the Las gene against the specified database with hard-coded parameters for the BLAST program. Further, manual analysis of the resulting BLAST search output files is also laborious and time consuming; we therefore, automated this step by developing a second PERL script (Additional file 2). This script automatically parsed all the BLAST output files and returned the Las sequences for which, no hits were found in other organisms. We refer to these sequences as probable unique sequences, because there are nearly no identical sequences found in other organisms (Figure 1).

Figure 1

Pictorial representation of the bioinformatics strategy employed to churn out the unique genic regions from Las genome. The input and output of each step are shown in oval or square boxes. Actions taken are noted to the left side of the arrow mark, while the information used is indicated to the right side of the arrow.

We performed the sequence similarity searches first by using stringent E-value of ≤ 1 × 10-3 against nt database (Figure 1). This search resulted in ~200 sequences that are unique to Las. This set of sequences is relatively high to validate experimentally; therefore, to further reduce the number of unique sequences, we performed the second sequence similarity search with a relaxed E-value of ≤ 1. This search resulted in 38 unique sequences. The E-value of ≤ 1 excludes the sequences with even little similarity to other organisms. Therefore, the resulting 38 unique sequences are considered unique to Las and constitute the promising candidates for qRT-PCR based detection (Figure 1).

We further searched the 38 unique sequences of Las against the phylogenetically closely related Lso, Lam, and Lcr. Because these organisms are closely related, we used the stringent E-value threshold of ≤ 1 × 10-3 for this similarity search. In order to achieve this E-value, the sequences need to be highly similar between the Las, Lso, Lam, and Lcr. Therefore, this close species filter procedure potentially eliminates all the Las sequence targets that could lead to false positive results in qRT-PCR based molecular diagnostic assays. Consequently, we further eliminated four conserved sequences from the list of 38 unique sequences, resulting in a total of 34 potential sequence signatures. We could not apply this close species filter step against Laf genome as its genome is yet to be sequenced.

Five (~15%) of the 34 unique gene sequences namely CLIBASIA_05545, CLIBASIA_05555, CLIBASIA_05560, CLIBASIA_05575 and CLIBASIA_05605 are in the prophage region of the Las genome. All these five unique sequences are located upstream of the genomic locus CLIBASIA_05610 encoding a phage terminase. There are possibly 30 genes that represent the complete prophage genome within the Las genome [25, 44], of which 16 open reading frames (ORFs) are upstream of the phage terminase, while the remaining 13 ORFs are downstream. The prophage genes CLIBASIA_05610 (primer pair 766 F and 766R) and CLIBASIA_05538 (primer pair LJ900F and LJ900R) have been targeted in previous studies by both conventional as well as qRT-PCR based assays [25, 44].

We further analyzed the genomic orientation of the 34 unique genes. This analysis revealed that 15 (~44%) of them are oriented on the sense strand, while the remaining 19 (~56%) were present on the anti-sense strand (Additional file 3: Figure S1). The sequence length of these unique genes ranged from 93 to 2595 base pairs (bp) (Additional file 4: Table S1).

Designing of Las specific primers and experimental validation of the specificity and sensitivity of qRT-PCR assay to detect Las

Based on the genome sequence of Las strain psy62, we designed 34 qRT-PCR primer pairs that specifically target the 34 unique sequences identified in our bioinformatic analyses (Additional file 4: Table S1). We designed the melting temperature (Tm) of each of these primers to range from 59°C to 65°C with an optimum of 62°C. The GC content of the primers ranged from 35% to 65% with an optimum of 50%. The PCR amplicon sizes for each primer set are between 84 to 185 bp (Additional file 4: Table S1).

In addition to the novel primers designed in this work, we also used a set of control primers that have been previously used in a qRT-PCR based detection of Las. These known primers include 16S rDNA pairs specific to the three different Candidatus Liberibacter species (HLBasf/r: Las, HLBamf/r: Lam and HLBaf/r: Laf) [23], β-operon (CQULA04f/r: β-operon) [26], intragenic repeats regions of the prophage sequence (LJ900f/r: Prophage) [25], and the primer pair specific to the plant cytochrome oxidase (COXf/r: COX) gene [23] as a positive endogenous control.

We performed qRT-PCR assays to test the specificity of the designed primers using total DNA extracted from Las-infected citrus plants as a template. To further validate the specificity of these primers, we also included total DNA from the phylogenetically closely related species Lam and Laf in our test. Additionally, DNA extracted from healthy citrus plant was used as a negative control, whereas water served as a no template control. The results of qRT-PCR assays are listed in Table 1.

Table 1 Specificity and sensitivity of the novel primers in the detection of Las as shown by qRT-PCR assay

Most of our novel custom designed primer pairs targeting the unique gene sequences were indeed found to be highly specific to Las, as assessed by qRT-PCR assays (Table 1). Among the 34 primer pairs, 29 produced amplicons only when Las-infected citrus plant DNA was used as a template, with an average CT value ranged from 19.48 to 27.47. Two primer pairs, P13 and P15, didn’t produce any amplicons under the standard conditions tested. The other three primer pairs, P19, P27 and P28, produced amplicons when Las or Laf infected plant DNA was used as a template, indicating P19, P27 and P28 could be used to detect both Las and Laf. We were unable to filter for cross-reactivity of P19, P27 and P28 in the bioinformatic analysis, because the Laf genome sequence is currently unavailable. With the exception of these three primer sets that showed amplicons with Laf template, none of the other primer sets produced any amplicons with DNA of Lam, Laf, and healthy citrus or water as template, which further confirms the specificity of these primers to the Las.

We further evaluated the specificity of these primer sets using DNA templates from various citrus associated fungal and bacterial pathogens including Colletotrichum acutatum KLA-207, Elsinoe fawcettii, Xanthomonas axonopodis pv. citrumelo 1381, X. citri subsp. citri strains 306, Aw, and A*. Only two primers sets, P20 and P21 showed unspecific amplification against template DNA extracted from fungal pathogen C. acutatum KLA-207 (Table 1). C. acutatum causes citrus blossom blight, post-bloom fruit drop and anthracnose symptoms that are phenotypically distinguishable from citrus HLB. The P20 and P21 were not filtered by the bioinformatic analysis since C. acutatum genome sequence was unavailable in the database. Because of the complexity of the natural microbial community and the limited number of sequences available in the current nucleotide sequence database, it is impossible to completely filter out all the potential false positives bioinformatically. However, false positives could be identified experimentally by combining the different sets of primer pairs by a consensus approach [37]. We eliminated these two primer sets from further evaluation in this study.

The melting temperature analysis of the amplicons produced from our novel primer set with Las as a template indicated that amplicons were of a single species. This suggests that there is no off target amplification for our primer pairs on the Las genome. Overall, the experimental validation of the 34 novel primer sets specific to unique targets revealed that 27 (~80%) of these targets are indeed specific to the Las genome (Table 1). This demonstrates the significance of the bioinformatics strategy employed here for identifying the suitable target regions for the detection of the bacteria by qRT-PCR based methods. These 27 novel primer pairs were selected for further characterization.

To test the sensitivity of our designed novel primers, serial dilutions of Las-infected psyllid DNA was used as a template in the qRT-PCR assay. This serial dilution qRT-PCR assay indicated that most of our novel primer pairs were able to detect Las up to 104 dilutions from the initial template DNA concentration, which is comparable to that of the primer set targeting Las 16S rDNA (Table 1). However, lower sensitivity was observed in the case of primer pairs P9, P12, P14 and P22, which were eliminated from further study. The remaining 23 primer pairs were able to detect Las up to 104 dilutions, with a correlation co-efficient (R2 >0.94) between the CT values and dilutions (Table 1). This demonstrates the high sensitivity of these 23 primers in the detection of Las.

qRT-PCR detection of Las from plant and psyllid DNA samples isolated from diverse locations in USA and China

In order to further demonstrate the degree of applicability of the 23 primer pairs in the detection of Las from infected biological material, we performed qRT-PCR on the various Las-infected plant and psyllid DNA samples. Considering the potential variation in nucleotide sequences of Las isolates in different geographic locations that might affect our detection due to the potential nucleotides changes of the target unique genes, we collected Las-infected plant DNA samples as tabulated in Table 2, from not only USA, but also from China, where Las was reported more than 100 years ago [1]. We tested the 23 primer pairs on 17 Las-infected plant DNA samples. Of these 17, 12 were collected from different locations in Florida, USA (Figure 2, Table 2), and the remaining five were collected from different locations in China (Table 2). Additionally, Las-infected psyllid DNA samples collected from five different locations in Florida, USA, were also included in the qRT-PCR assays (Table 3, Figure 2).

Table 2 qRT-PCR detection of Las from plant samples that were collected from different locations in USA and China
Figure 2

Schematic representation of the plant and the psyllid samples collected from Florida. Las-infected plant DNA samples were collected from 12 different locations and psyllids from 5 different locations in Florida, USA. The color shaded symbols for representative plant and psyllid samples are based on their average infection level across all the primer pairs tested based on CT values.

Table 3 qRT-PCR detection of Las from psyllid DNA samples that were collected from different locations in Florida, USA

All the 23 primer pairs detected Las from all 12 Florida HLB diseased plant samples (Table 2) and 5 psyllid DNA samples (Table 3) in a qRT-PCR assay, which further validated the detection applicability of our novel primers (Figure 2). However, 4 of the 23 primer pairs (P1, P7, P8 and P10) failed to produce amplicons with the infected plant DNA sample from Jiangxi and Guangdong Province, China (Table 2). Primer pair P3 produced no amplicon with Jiangxi sample, and produced unspecific amplicon with the Guangdong sample (with an altered PCR product size, data not shown). Interestingly, all these 5 primer pairs target the genes located in prophage region of the Las genome (Additional file 3). These primers (P1, P3, P7, P8 and P10) based on prophage genes could detect Las from Florida, but not from Jiangxi and Guangdong province, China. This is consistent with previous report [44], that prophage was detected in only 15.8% of the 120 HLB diseased citrus samples acquired in Guangdong Province, China, but was detected in 97.4% of the 39 Las positive citrus samples acquired in Yunnan Province, China. This suggests that those prophage genes are not universally present in all strains of Las. Alternately, the prophage sequences were found to be highly variable among the strains tested.


We have successfully designed 18 novel primer pairs, which are specific to Las. These primers will provide an additional arsenal to qRT-PCR based detection of Las-infected plants and psyllids. Compared to the commonly used primers based on 16S rDNA and β-operon, the 18 primers developed in this study have comparable sensitivity. Moreover, these primers could successfully differentiate Las from Lam, Laf and other common microbes associated with citrus.



The nucleotide sequences of Las with accession number NC_012985 [29, 45], Lso with accession number NC_014774 [33], Lcr with accession number NC_019907 and comprehensive nucleotide (nt) database (26th July 2012) were downloaded from the NCBI ftp server ( The stand-alone BLAST [42, 43] was used to search the Las genes against nt, Lso and Lcr databases using a custom-made PERL script 1 (Additional file 1) by varying the E-value with all other parameters kept to a default value. The output files of the BLAST searches were further parsed using a second custom-made PERL script 2 (Additional file 2).

Plant and psyllid materials and extraction of DNA

Las infected citrus leaf samples with typical visible symptoms were collected from 2 years old infected sweet orange (Citrus sinensis) plants maintained at the Citrus Research and Education Center (CREC), Lake Alfred, Florida, USA. As a negative control, the leaves from healthy citrus plants were collected from pathogen-free seedlings grown in the healthy plant greenhouse maintained at CREC, Lake Alfred, Florida, USA. The Laf and Lam infected samples were obtained from South Africa and Brazil respectively. The total DNA from the leaves of citrus was extracted using the protocol mentioned elsewhere [46]. Briefly, the leaves were washed under tap water and surface sterilized in 35% bleach (2% active Chlorine) and 70% (v/v) ethanol for 2 min each. The sterilized leaves were further rinsed three times in sterile water. The midribs from the leaf samples were separated and cut into small pieces. Approximately 100 mg of midrib pieces were used from each sample to extract the DNA using the Wizard® genomics DNA purification kit (Promega, Madison, WI, USA). The extracted DNA was suspended in 100 μl H2O.

Las infected psyllids (Diaphorina citri) were maintained on confirmed Las-infected sweet orange plants at the CREC, Lake Alfred, FL, USA. In this work, 16 psyllids (around 20 mg) were pooled and the total DNA was extracted using a DNeasy Blood & Tissue Kit (Qiagen, Valencia, CA). The extracted DNA was suspended in 100 μl H2O. The quality and quantity of the extracted DNA was determined using a NanoDrop™ 1000 spectrophotometer (NanoDrop Technologies, Inc., Wilmington, DE).

Quantitative real-time polymerase chain reaction (qRT-PCR)

Gene specific primers were designed using PrimerQuestSM from Integrated DNA technologies (IDT), Coralville, Iowa (Additional file 4: Table S1). qRT-PCR experiments were performed using ABI PRISM 7500 FAST Real-time PCR System (Applied Biosystems, Foster City, CA, US) in a 96-well plate by using an absolute quantification protocol. The reaction mixture in each well contained 12.5 μL 2x FAST SYBR® Green PCR Master Mix reagent (Applied Biosystems), 2 μL DNA template (~30 ng), 0.625 μL of 10 μM of each gene-specific primer pair in a final volume of 25 μL. The standard thermal profile for all amplifications was followed, which involved 95°C for 20 min followed by 40 cycles of 95 °C for 3 sec, and 50°C for 30 sec. All assays were performed in triplicates.

Melting curve analysis was performed using ABI PRISM 7500 FAST Real-time PCR System Software version SDS v1.4 21 CFR Part 11 Module (Applied Biosystems®) to characterize the amplicons produced in a PCR reaction.


  1. 1.

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

    Google Scholar 

  2. 2.

    do Carmo Teixeira D, Luc Danet J, Eveillard S, Cristina Martins E, de Jesus Junior WC, Takao Yamamoto P, Aparecido Lopes S, Beozzo Bassanezi R, Juliano Ayres A, Saillard C, Bové JM: Citrus huanglongbing in Sao Paulo State, Brazil: PCR detection of the 'Candidatus’ Liberibacter species associated with the disease. Mol Cell Probes. 2005, 19 (3): 173-179. 10.1016/j.mcp.2004.11.002.

    Article  PubMed  Google Scholar 

  3. 3.

    Jagoueix S, Bové JM, Garnier M: Comparison of the 16S/23S ribosomal intergenic regions of “Candidatus Liberobacter asiaticum” and “Candidatus Liberobacter africanum”, the two species associated with citrus huanglongbing (greening) disease. Int J Syst Bacteriol. 1997, 47 (1): 224-227. 10.1099/00207713-47-1-224.

    Article  CAS  PubMed  Google Scholar 

  4. 4.

    Lopes SA, Frare GF, Bertolini E, Cambra M, Fernandes NG, Ayres AJ, Marin DR, Bové JM: Liberibacters associated with citrus Huanglongbing in Brazil: 'Candidatus Liberibacter asiaticus’ is heat tolerant, 'Ca. L. americanus’ is heat sensitive. Plant Dis. 2009, 93 (3): 257-262. 10.1094/PDIS-93-3-0257.

    Article  Google Scholar 

  5. 5.

    Tatineni S, Sagaram US, Gowda S, Robertson CJ, Dawson WO, Iwanami T, Wang N: In planta distribution of 'Candidatus Liberibacter asiaticus’ as revealed by polymerase chain reaction (PCR) and real-time PCR. Phytopathology. 2008, 98 (5): 592-599. 10.1094/PHYTO-98-5-0592.

    Article  CAS  PubMed  Google Scholar 

  6. 6.

    Manjunath KL, Halbert SE, Ramadugu C, Webb S, Lee RF: Detection of 'Candidatus Liberibacter asiaticus’ in Diaphorina citri and its importance in the management of citrus huanglongbing in Florida. Phytopathology. 2008, 98 (4): 387-396. 10.1094/PHYTO-98-4-0387.

    Article  CAS  PubMed  Google Scholar 

  7. 7.

    McClean APD, Oberholzer PCJ: Citrus psylla, a vector of the greening disease of sweet orange. South African J of Agricultural Sci. 1965, 8: 297-298.

    Google Scholar 

  8. 8.

    Shi J, Pagliaccia D, Morgan R, Qiao Y, Pan S, Vidalakis G, Ma W: Novel diagnosis for Citrus Stubborn Disease by detection of a Spiroplasma citri-secreted protein. Phytopathology. 2014, 104 (2): 188-195. 10.1094/PHYTO-06-13-0176-R.

    Article  CAS  PubMed  Google Scholar 

  9. 9.

    Chen J, Pu X, Deng X, Liu S, Li H, Civerolo E: A phytoplasma related to 'Candidatus phytoplasma asteri’ detected in citrus showing Huanglongbing (yellow shoot disease) symptoms in Guangdong, P. R. China. Phytopathology. 2009, 99 (3): 236-242. 10.1094/PHYTO-99-3-0236.

    Article  CAS  PubMed  Google Scholar 

  10. 10.

    Teixeira DC, Wulff NA, Martins EC, Kitajima EW, Bassanezi R, Ayres AJ, Eveillard S, Saillard C, Bové JM: A phytoplasma closely related to the pigeon pea witches’-broom phytoplasma (16Sr IX) is associated with citrus huanglongbing symptoms in the state of Sao Paulo, Brazil. Phytopathology. 2008, 98 (9): 977-984. 10.1094/PHYTO-98-9-0977.

    Article  CAS  PubMed  Google Scholar 

  11. 11.

    Wang N, Trivedi P: Citrus huanglongbing: a newly relevant disease presents unprecedented challenges. Phytopathology. 2013, 103 (7): 652-665. 10.1094/PHYTO-12-12-0331-RVW.

    Article  PubMed  Google Scholar 

  12. 12.

    Gottwald TR, da Graca JV, Bassanezi RB: Citrus huanglongbing: the pathogen and its impact. Plant Health Progress. 2007, doi:10.1094/PHP-2007-0906-1001-RV

    Google Scholar 

  13. 13.

    Okuda M, Mitsuhito M, Tanaka Y, Subandiyah S, Iwanami T: Characterization of the tufB-secE-nusG-rplKSJL-ropB gene cluster of the citrus greening organism and detection by loop-mediated isothermal amplification. Plant Dis. 2005, 89 (7): 705-711. 10.1094/PD-89-0705.

    Article  CAS  Google Scholar 

  14. 14.

    Villechanoux S, Garnier M, Renaudin J, Bové J: Detection of several strains of the bacterium-like organism of citrus greening disease by DNA probes. Curr Microbiol. 1992, 24 (2): 89-95. 10.1007/BF01570903.

    Article  CAS  Google Scholar 

  15. 15.

    Garnier M, Martin-Gros G, Bové JM: Monoclonal antibodies against the bacterial-like organism associated with citrus greening disease. Ann Inst Pasteur Microbiol. 1987, 138 (6): 639-650. 10.1016/0769-2609(87)90142-6.

    Article  CAS  PubMed  Google Scholar 

  16. 16.

    Gurtler V, Stanisich VA: New approaches to typing and identification of bacteria using the 16S-23S rDNA spacer region. Microbiology. 1996, 142 (1): 3-16. 10.1099/13500872-142-1-3.

    Article  PubMed  Google Scholar 

  17. 17.

    Hocquellet A, Toorawa P, Bové JM, Garnier M: Detection and identification of the two Candidatus Liberobacter species associated with citrus huanglongbing by PCR amplification of ribosomal protein genes of the beta operon. Mol Cell Probes. 1999, 13 (5): 373-379. 10.1006/mcpr.1999.0263.

    Article  CAS  PubMed  Google Scholar 

  18. 18.

    Hung TH, Wu ML, Su HJ: Development of a rapid method for the diagnosis of Citrus Greening Disease using the Polymerase Chain Reaction. J Phytopathol. 1999, 147 (10): 599-604. 10.1046/j.1439-0434.1999.00435.x.

    Article  Google Scholar 

  19. 19.

    Jagoueix S, Bové JM, Garnier M: PCR detection of the two 'Candidatus’ Liberobacter species associated with greening disease of citrus. Mol Cell Probes. 1996, 10 (1): 43-50. 10.1006/mcpr.1996.0006.

    Article  CAS  PubMed  Google Scholar 

  20. 20.

    Shang S, Fu J, Dong G, Hong W, Du L, Yu X: Establishment and analysis of specific DNA patterns in 16S-23S rRNA gene spacer regions for differentiating different bacteria. Chin Med J (Engl). 2003, 116 (1): 129-133.

    CAS  Google Scholar 

  21. 21.

    Fujikawa T, Iwanami T: Sensitive and robust detection of citrus greening (huanglongbing) bacterium “Candidatus Liberibacter asiaticus” by DNA amplification with new 16S rDNA-specific primers. Mol Cell Probes. 2012, 26 (5): 194-197. 10.1016/j.mcp.2012.06.001.

    Article  CAS  PubMed  Google Scholar 

  22. 22.

    Teixeira DC, Saillard C, Couture C, Martins EC, Wulff NA, Eveillard-Jagoueix S, Yamamoto PT, Ayres AJ, Bové JM: Distribution and quantification of Candidatus Liberibacter americanus, agent of huanglongbing disease of citrus in Sao Paulo State, Brasil, in leaves of an affected sweet orange tree as determined by PCR. Mol Cell Probes. 2008, 22 (3): 139-150. 10.1016/j.mcp.2007.12.006.

    Article  CAS  PubMed  Google Scholar 

  23. 23.

    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-115. 10.1016/j.mimet.2005.10.018.

    Article  CAS  PubMed  Google Scholar 

  24. 24.

    Li W, Li D, Twieg E, Hartung JS, Levy L: Optimized quantification of unculturable Candidatus Liberibacter spp. in host plants using real-time PCR. Plant Dis. 2008, 92 (6): 854-861. 10.1094/PDIS-92-6-0854.

    Article  CAS  Google Scholar 

  25. 25.

    Morgan JK, Zhou L, Li W, Shatters RG, Keremane M, Duan YP: Improved real-time PCR detection of 'Candidatus Liberibacter asiaticus’ from citrus and psyllid hosts by targeting the intragenic tandem-repeats of its prophage genes. Mol Cell Probes. 2012, 26 (2): 90-98. 10.1016/j.mcp.2011.12.001.

    Article  CAS  PubMed  Google Scholar 

  26. 26.

    Wang Z, Yin Y, Hu H, Yuan Q, Peng G, Xia Y: Development and application of molecular-based diagnosis for 'Candidatus Liberibacter asiaticus’, the causal pathogen of citrus huanglongbing. Plant Pathol. 2006, 55 (5): 630-638. 10.1111/j.1365-3059.2006.01438.x.

    Article  CAS  Google Scholar 

  27. 27.

    Lane D: 16 s/23s rRNA sequencing. Nucleic acid techniques in bacterial systematics. Edited by: Stackebrandt E, Goodfellow M. 1991, West Sussex, United Kingdom: John Wiley & Sons, 115-175.

    Google Scholar 

  28. 28.

    Nageswara-Rao M, Irey M, Garnsey SM, Gowda S: Candidate gene markers for Candidatus Liberibacter asiaticus for detecting citrus greening disease. J Biosci. 2013, 38 (2): 229-237. 10.1007/s12038-013-9315-x.

    Article  CAS  PubMed  Google Scholar 

  29. 29.

    Duan Y, Zhou L, Hall DG, Li W, Doddapaneni H, Lin H, Liu L, Vahling CM, Gabriel DW, Williams KP, Dickerman A, Sun Y, Gottwald T: Complete genome sequence of citrus huanglongbing bacterium, 'Candidatus Liberibacter asiaticus’ obtained through metagenomics. MPMI. 2009, 22 (8): 1011-1020. 10.1094/MPMI-22-8-1011.

    Article  CAS  PubMed  Google Scholar 

  30. 30.

    Lin H, Han CS, Liu B, Lou B, Bai X, Deng C, Civerolo EL, Gupta G: Complete genome sequence of a Chinese strain of “Candidatus Liberibacter asiaticus”. Genome Announc. 2013, 1 (2): doi:10.1128/genomeA.00184-13.

    Google Scholar 

  31. 31.

    Lin H, Coletta-Filho HD, Han CS, Lou B, Civerolo EL, Machado MA, Gupta G: Draft genome sequence of “Candidatus Liberibacter americanus” bacterium associated with Citrus Huanglongbing in Brazil. Genome Announc. 2013, 1 (3): doi:10.1128/genomeA.00275-13

    Google Scholar 

  32. 32.

    Leonard MT, Fagen JR, Davis-Richardson AG, Davis MJ, Triplett EW: Complete genome sequence of Liberibacter crescens BT-1. Stand Genomic Sci. 2012, 7 (2): 271-283. 10.4056/sigs.3326772.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  33. 33.

    Lin H, Lou B, Glynn JM, Doddapaneni H, Civerolo EL, Chen C, Duan Y, Zhou L, Vahling CM: The complete genome sequence of 'Candidatus Liberibacter solanacearum’, the bacterium associated with potato zebra chip disease. PLoS One. 2011, 6 (4): e19135-10.1371/journal.pone.0019135.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  34. 34.

    Ho CC, Yuen KY, Lau SK, Woo PC: Rapid identification and validation of specific molecular targets for detection of Escherichia coli O104:H4 outbreak strain by use of high-throughput sequencing data from nine genomes. J Clin Microbiol. 2011, 49 (10): 3714-3716. 10.1128/JCM.05062-11.

    PubMed Central  Article  PubMed  Google Scholar 

  35. 35.

    Phillippy AM, Ayanbule K, Edwards NJ, Salzberg SL: Insignia: a DNA signature search web server for diagnostic assay development. Nucleic Acids Res. 2009, 37 (suppl 2): W229-W234.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  36. 36.

    Phillippy AM, Mason JA, Ayanbule K, Sommer DD, Taviani E, Huq A, Colwell RR, Knight IT, Salzberg SL: Comprehensive DNA signature discovery and validation. PLoS Comput Biol. 2007, 3 (5): e98-10.1371/journal.pcbi.0030098.

    PubMed Central  Article  PubMed  Google Scholar 

  37. 37.

    Pritchard L, Holden NJ, Bielaszewska M, Karch H, Toth IK: Alignment-free design of highly discriminatory diagnostic primer sets for Escherichia coli O104:H4 outbreak strains. PLoS One. 2012, 7 (4): e34498-10.1371/journal.pone.0034498.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  38. 38.

    Slezak T, Kuczmarski T, Ott L, Torres C, Medeiros D, Smith J, Truitt B, Mulakken N, Lam M, Vitalis E, Zemla A, Zhou CE, Gardner S: Comparative genomics tools applied to bioterrorism defence. Brief Bioinform. 2003, 4 (2): 133-149. 10.1093/bib/4.2.133.

    Article  CAS  PubMed  Google Scholar 

  39. 39.

    Vijaya Satya R, Kumar K, Zavaljevski N, Reifman J: A high-throughput pipeline for the design of real-time PCR signatures. BMC Bioinforma. 2010, 11: 340-10.1186/1471-2105-11-340.

    Article  Google Scholar 

  40. 40.

    Vijaya Satya R, Zavaljevski N, Kumar K, Bode E, Padilla S, Wasieloski L, Geyer J, Reifman J: In silico microarray probe design for diagnosis of multiple pathogens. BMC Genomics. 2008, 9: 496-10.1186/1471-2164-9-496.

    PubMed Central  Article  PubMed  Google Scholar 

  41. 41.

    Nielsen R: Molecular signatures of natural selection. Annu Rev Genet. 2005, 39: 197-218. 10.1146/annurev.genet.39.073003.112420.

    Article  CAS  PubMed  Google Scholar 

  42. 42.

    Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215 (3): 403-410.

    Article  CAS  PubMed  Google Scholar 

  43. 43.

    Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25 (17): 3389-3402. 10.1093/nar/25.17.3389.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  44. 44.

    Liu R, Zhang P, Pu X, Xing X, Chen J, Deng X: Analysis of a prophage gene frequency revealed population variation of 'Candidatus Liberibacter asiaticus’ from two citrus-growing provinces in China. Plant Dis. 2010, 95 (4): 431-435.

    Article  Google Scholar 

  45. 45.

    Tyler HL, Roesch LF, Gowda S, Dawson WO, Triplett EW: Confirmation of the sequence of 'Candidatus Liberibacter asiaticus’ and assessment of microbial diversity in Huanglongbing-infected citrus phloem using a metagenomic approach. MPMI. 2009, 22 (12): 1624-1634. 10.1094/MPMI-22-12-1624.

    Article  CAS  PubMed  Google Scholar 

  46. 46.

    Kim JS, Wang N: Characterization of copy numbers of 16S rDNA and 16S rRNA of Candidatus Liberibacter asiaticus and the implication in detection in planta using quantitative PCR. BMC Research Notes. 2009, 2: 37-10.1186/1756-0500-2-37.

    PubMed Central  Article  PubMed  Google Scholar 

Download references


We thank Dr. Nelson A. Wulff, Fundecitrus – Fundo de Defesa da Citricultura, Sao Paulo, Brazil, for kindly providing the Lam DNA. DNA samples of fungal pathogens Colletotrichum acutatum KLA-207, Elsinoe fawcettii were kindly provided by Dr. Kuang-Ren Chung. We also thank Vladimir Kolbasov for the technical assistance in DNA isolation. This work was supported by Citrus Research and Development Foundation.

Author information



Corresponding author

Correspondence to Nian Wang.

Additional information

Competing interests

We declare no competing interests.

Authors’ contributions

NW conceived and coordinated the work and wrote the manuscript. SK designed, performed bioinformatic analysis and wrote the manuscript. SK, QY and NR performed qRT-PCR experiments. SK, QY, XD, CR, TE, MR, MI, GP, and CR participated in experimental design, manuscript writing and provided reagents. All authors read and approved the final manuscript.

Electronic supplementary material

Additional file 1:PERL script 1 facilitates the similarity search in an automated fashion. This script performs similarity searches against the specified nucleotide sequence database using a stand-alone BLAST program for each of the input gene sequences from the Las genome. (TXT 4 KB)

Additional file 2:PERL script 2 facilitates the identification of unique genes to Las. This script facilitates the identification of unique genes by automatically parsing all the BLAST output files generated from the Additional file 1 PERL script 1and returns the unique gene sequences with no similarity to the DNA sequences of other organisms. (TXT 3 KB)

Additional file 3: Figure S1: Snapshot of the unique genes identified by bioinformatics is shown in the context of the whole genome of Las. The absolute positions of the regions are shown. The novel unique regions of Las identified in this study are shown in bluish green, while the currently known targets are colored in green. (PDF 1 MB)

Additional file 4: Table S1: Custom designed primer pairs specific to the unique sequences of Las identified by bioinformatic analysis. The forward and reverse primer pair for each of the unique genic regions is given. The product size for each of the primers is shown along with the %GC content. (DOC 62 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Rights and permissions

Reprints and Permissions

About this article

Cite this article

Kogenaru, S., Yan, Q., Riera, N. et al. Repertoire of novel sequence signatures for the detection of Candidatus Liberibacter asiaticus by quantitative real-time PCR. BMC Microbiol 14, 39 (2014).

Download citation


  • Detection system
  • Diagnostic
  • Candidatus Liberibacter asiaticus
  • Greening
  • Huanglongbing
  • Bacteria
  • Psyllid
  • Citrus