# Novel "candidatus liberibacter' species identified in the australian eggplant psyllid, acizzia solanicola

Jacqueline Morris, Jason Shiller, Rachel Mann, Grant Smith, Alan Yen,
Brendan Rodoni

► To cite this version:

Jacqueline Morris, Jason Shiller, Rachel Mann, Grant Smith, Alan Yen, et al.. Novel "candidatus liberibacter' species identified in the australian eggplant psyllid, acizzia solanicola. Microbial Biotechnology, 2017, 10 (4), pp.833-844. 10.1111/1751-7915.12707. hal-02624445

HAL Id: hal-02624445

https://hal.inrae.fr/hal-02624445

Submitted on 26 May 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Novel 'Candidatus Liberibacter' species identified in the Australian eggplant psyllid, *Acizzia solanicola*

Jacqueline Morris,<sup>1,2,3</sup> Jason Shiller,<sup>1,3,4</sup>

Rachel Mann,<sup>1,3</sup> Grant Smith,<sup>1,5,6</sup> Alan Yen<sup>1,2,3</sup> and
Brendan Rodoni<sup>1,2,3,*</sup>

## Summary

A novel candidate species of the liberibacter genus, 'Candidatus Liberibacter brunswickensis' (CLbr), was identified in the Australian eggplant psyllid, *Acizzia solanicola*. This is the first discovery of a species belonging to the liberibacter genus in Australia and the first report of a liberibacter species in the psyllid genus *Acizzia*. This new candidate liberibacter species has not been associated with plant disease, unlike other psyllid-vectored species in the genus including 'Candidatus Liberibacter asiaticus' (CLas), 'Candidatus Liberibacter africanus' (CLaf) and 'Ca. Liberibacter solanacearum' (CLso). This study describes novel generic liberibacter genus primers, used to screen Australian psyllids for the presence of microflora that may confound diagnosis of exotic pathogens. CLbr forms a unique clade in the liberibacter genus based on phylogenetic analysis of the 16S ribosomal ribonucleic acid (rRNA) region and multilocus sequence analysis (MLSA) of seven highly conserved genes, *dnaG*, *gyrB*, *mutS*, *nusG*, *rplA*, *rpoB*

and *tufB*. The MLSA mapping approach described in this article was able to discriminate between two 'Ca. Liberibacter' species within a metagenomic data set and represents a novel approach to detecting and differentiating unculturable species of liberibacter. Further, CLbr can confound the Li *et al.* (2006) quantitative PCR (qPCR) diagnostic tests for CLas and CLaf.

## Introduction

The liberibacter genus includes species of Gram-negative alpha Proteobacteria that are known plant pathogens and endophytic species. There are currently seven species which are found in most regions of the world, excluding mainland Australia: *Liberibacter crescens* (Lcr), 'Candidatus Liberibacter africanus' (CLaf), 'Ca. Liberibacter americanus' (CLam), 'Ca. Liberibacter asiaticus' (CLas), 'Ca. Liberibacter caribbeanus' (CLca), 'Ca. Liberibacter solanacearum' (CLso, synonymous to 'Ca. Liberibacter psyllorous', CLps) and 'Ca. Liberibacter europaeus' (CLeu) (Jagoueix *et al.*, 1994; Duan *et al.*, 2009; Liefting *et al.*, 2009a,b; Munyaneza *et al.*, 2010, 2014; Raddadi *et al.*, 2011; Leonard *et al.*, 2012; Phahladira *et al.*, 2012; Tahzima *et al.*, 2014; Keremane *et al.*, 2015). With the exception of Lcr, species in the liberibacter genus are associated with psyllids (Hemiptera: Psylloidea). In addition, Lcr is the only liberibacter that can be grown as an axenic culture (Leonard *et al.*, 2012). Due to the fastidious nature of the predominant liberibacter species, next-generation sequencing (NGS) approaches have played a large role in the investigation of this genus. NGS will aid the expansion of this genus and will continue to inform information on current species, their insect and plant host range.

CLso, CLas, CLam and CLaf are all exotic plant pathogens (EPPs) in Australia. In the event of an incursion, these phytopathogenic *Ca. Liberibacter* species would pose a severe threat to the Australian potato, tomato, carrot (CLso) and citrus industries (CLas, CLaf or CLam). None of the psyllid species that are known to host and vector the phytopathogenic 'Ca. Liberibacter' species are currently present in mainland Australia, yet Australia is a centre of psyllid diversity with an estimated 446 species present (Yen, 2002). Little is known about the microflora of native Australian psyllids, if the native microflora could confound diagnostic tests for

Received 5 December, 2016; accepted 27 February, 2017. *For correspondence. E-mail brendan.rodoni@ecodev.vic.gov.au; Tel. +61 3 9032 7319; Fax +61 3 9032 7604.

*Microbial Biotechnology* (2017) **10**(4), 833–844
doi:10.1111/1751-7915.12707

## Funding information

We would like to acknowledge the support of the Australian Government's Cooperative Research Centres Program.

<footnote n="1"><sup>1</sup>Plant Biosecurity Cooperative Research Centre, LPO Box 5012, Bruce, Australian Capital Territory 2617, Australia.</footnote>

<footnote n="2"><sup>2</sup>La Trobe University, AgriBio, 5 Ring Road, Bundoora, Victoria 3083, Australia.</footnote>

<footnote n="3"><sup>3</sup>Agriculture Victoria, AgriBio, 5 Ring Road, Bundoora, Victoria 3083, Australia.</footnote>

<footnote n="4"><sup>4</sup>INRA/Université d'Angers - IRHS Batiment C, 42 rue Georges Morel, Beaucouzé 49071, France.</footnote>

<footnote n="5"><sup>5</sup>Plant & Food Research Lincoln, Gerald St, Lincoln 7608, New Zealand.</footnote>

<footnote n="6"><sup>6</sup>Better Border Biosecurity, Lincoln 7608, New Zealand.</footnote>
phytopathogenic ‘Ca. Liberibacter’ species and whether native Australian psyllids could vector or host species of liberibacter. All diagnostic tests for species of liberibacter have been developed outside of Australia (Bové, 2006; Hocquellet et al., 1999; Jagoueix et al., 1994; Li et al., 2009; Morgan et al., 2012; Thompson et al., 2015; Wen et al., 2009, 2013). Therefore, Australian psyllid microflora that is not associated with plant diseases could yield false positives with the diagnostic tests, resulting in significant, unnecessary responses and economic and trade consequences. Thus, understanding the microflora of Australian psyllids is important for both biosecurity preparedness and response management.

*Acizzia solanicola*, the eggplant psyllid, is a native Australian psyllid that has broadened its host range from the native solanaceous host plant the rock nightshade, *Solanum pterophilum*, to eggplants, *Solanum melongena* (and other introduced ornamental and weed solanaceous plant species) (Taylor and Kent, 2013). Additional hosts for *A. solanicola* include wild tobacco bush, *S. mauritianum*, cape gooseberry, *Physalis peruviana*, and an undetermined species of angel’s trumpet *Brugmansia* (Kent and Taylor, 2010; Taylor and Kent, 2013). *Bactericera cockerelli*, the tomato potato psyllid, is able to vector CLso, feeding on solanaceous plants including the eggplant (Liefting et al., 2009a,b). Due to the overlap of host plant species with a known plant pathogen vector, *A. solanicola* was chosen as a case study to investigate whether native Australian psyllids may harbour microflora that are closely related to phytopathogenic species of liberibacter. This study is the first to use generic primers to screen native Australian psyllids for liberibacter species, followed by multilocus sequence analysis (MLSA) to detect a new candidate species of liberibacter that is not associated with disease symptoms.

## Results

### Development of generic Liberibacter genus primers

The generic liberibacter genus primer pair, LG774F-LG1463R, developed in this study successfully amplified the expected 684 base pair (bp) amplicons from the 16S ribosomal ribonucleic acid (rRNA) region of CLas, CLso and CLeu. When applied to DNA extracts from *A. solanicola* from Brunswick, Victoria, the LG774F-LG1463R primer pair amplified a 684-bp product in 17 of 37 individual psyllid DNA extracts tested. The universal 16S bacterial primer set, FD2-RP1 (Weisburg et al., 1991), successfully amplified a product of approximately 1500 bp for all 37 *A. solanicola* DNA extracts.

*Acizzia solanicola* DNA extracts from six additional sites tested from Victoria (VIC) and New South Wales (NSW) in south-eastern Australia were tested using the liberibacter generic primers, LG774F-LG1463R and universal 16S

bacterial primer set, FD2-RP1. All 31 extracts successfully amplified the 1500-bp product for FD2-RP1. None of the DNA extracts from the six additional sites amplified the 684-bp product for LG774F-LG1463R.

Four of the 17 LG774F-LG1463R positive PCR products from Brunswick *A. solanicola* were sequenced in both directions, and all four consensus sequences were identical to each other. Nucleotide Basic Local Alignment Search Tool (BLASTn) analysis of the longest consensus sequence (607 bp in length) determined that the LG774F-LG1463R amplicon from *A. solanicola*, had a per cent identity of 99.0% to CLso CLso-ZC1 (NR_074494.1 and CP002371.1), CLso Garden-City-KS-5 (FJ914619.1), CLso Garden-City-KS-1 (EU921626.1), CLso NZ083338 (EU935004.1), CLso TXZC018 (EU980389.1), CLps PRR1 (EU812559.1), CLps Tom100 (EU812558.1), CLso Tx15 (EU812556.1) and CLso NZ082226 (EU834130.1).

### Cloning and phylogenetic analyses of the 16S rRNA region of bacteria within *A. solanicola*

To obtain a longer region of 16S rRNA for phylogenetic analysis, the 16S rRNA FD2-RP1 amplicons (of all bacteria amplified) from two LG774F-LG1463R positive *A. solanicola* DNA extracts from Brunswick were cloned. Cloned colonies were screened for the presence of sequence belonging to the liberibacter genus using the LG774F-LG1463R primer pair. Positive clones were selected for Sanger sequencing of the full-length FD2-RP1 amplicon. BLASTn was used to compare the longest 16S rRNA consensus sequence generated (1464 bp) (KY077741) to the GenBank non-redundant (nr) database. The highest sequence identity was to CLaf PTSAPSY (NZ_CP004021.1) with 98.3% (100% query cover and E-value of 0.0). This was closely followed by CLas A4 (NZ_CP010804.1) with 98.1%, CLas Ishi-1 (NZ_AP014595.1) with 98.1%, CLas psy62 (NC_012985.3) with 98.1%, CLas gxpsy (CP04005.1) with 98.0% and CLas GuangXi-GL-1 (DQ778016.1) with 97.7% (all with 100% query cover and E-value of 0.0).

Phylogenetic analysis of the cloned 16S rRNA region revealed that the ‘*Ca. Liberibacter*’ species detected in *A. solanicola* is a relatively early branching, unique lineage sequence, differentiating prior to the CLas and CLaf clade (Fig. 1).

### Assembly of ‘Ca. Liberibacter’ genes for multilocus sequence analyses

The complete sequence of a further seven bacterial genes, dnaG (KY172970), gyrB (KY172971), mutS (KY172972), nusG (KY172973), rplA (KY172974), rpoB (KY172975) and tufB (KY172976), was assembled from the adapter and quality trimmed *A. solanicola*
<figure><img src="image_2.png" /><figcaption>**Fig. 1.** Phylogenetic analysis of the cloned 16S rRNA region of the *A. solanicola*-associated *Ca. Liberibacter* species. Unrooted maximum likelihood tree for the 16S rRNA region (1403 bp) was generated in MEGA6 [based on 1000 bootstraps and the Kimura 2-parameter model with Gamma distribution and Invariable sites (K2 + G+)] The tree includes the cloned 16S rRNA region of the *A. solanicola*-associated *Ca. Liberibacter* species, all known species of the liberibacter genus and close relatives as specified in Table 2.</figcaption></figure>

metagenome NextSeq data set (1 182 585 144 reads) (Fig. 2A). Of the assembled contigs for each gene, a single contig with full gene coverage and the highest nucleotide identity with corresponding genes from known 'Ca. Liberibacter' species was selected. BLASTn results revealed the highest nucleotide identity for six of the seven genes was to CLaf PTSAPSY (NZ_CP004021.1); dnaG (84.7%), gyrB (87.1%), mutS (84.7%), nusG (85.1%), rpoB (87.7%) and tufB (85.1%). The highest nucleotide identity for rplA (87.1%) was to CLas A4 (NZ_CP010804.1).

To validate the assembly of the seven selected genes from the *A. solanicola* metagenome NextSeq data set, simulated reads of CLaf PTSAPSY (NZ_CP004021.1) genome were run through the bioinformatics workflow (Fig. 2B). All seven genes were successfully reassembled for the simulated reads of CLaf PTSAPSY (NZ_CP004021.1) with a nucleotide identity of 100.0% to CLaf PTSAPSY (NZ_CP004021.1). When the simulated reads were spiked into the *A. solanicola* metagenome NextSeq data set (Fig. 2C), multiple contigs were assembled for each of the seven genes. The BLASTn searches confirmed the presence of assembled contigs

for each gene for both the *A. solanicola*-associated *Ca. Liberibacter* species and for CLaf PTSAPSY (NZ_CP004021.1) and clearly differentiated between the two species, with identical nucleotide sequence identity to the individual analyses.

*Multilocus sequence analysis and phylogenetic analyses for the *A. solanicola*-associated *Ca. Liberibacter* species*

No recombination events were recognized in Recombination Detection Program 4.56 (RPD4) (Martin *et al*., 2015) for the DNA sequences of the seven genes incorporated into the MLSA for the liberibacter genus analysed.

A distance matrix of the concatenated contigs representing all seven genes covering 12 379 bp determined that the *A. solanicola*-associated *Ca. Liberibacter* species has 86.1%, 84.1%, 82.2%, 75.7% and 69.5% sequence identity to CLaf, CLas, CLso, CLam and Lcr respectively (Table 1). The inclusion of four CLas genomes in the distance matrix showed high sequence identity within a species, with a sequence identity of 100%
<figure><img src="image_3.png" /><figcaption>**Fig. 2.** Phylogenetic analysis of the multilocus sequence of the *A. solanicola*-associated 'Ca. Liberibacter' species. Maximum likelihood tree of seven concatenated loci, *dnaG*, *gyrB*, *mutS*, *nusG*, *rplA*, *rpoB* and *tufB* (a total of 12 617 bp), was generated in MEGA6 [based on 1000 bootstraps, general time reversible model with gamma distribution and invariable sites (GTR+G+I)]. The tree includes the *A. solanicola*-associated 'Ca. Liberibacter' species, species the liberibacter genus and close relatives as specified in Table 2, this tree was rooted to *E. coli* str. O157:H7 Sakai.</figcaption></figure>

for the seven concatenated genes between the CLas isolates (Table 1, highlighted in bright blue).

Phylogenetic analysis of the concatenated genes differentiated the *A. solanicola*-associated 'Ca. Liberibacter' species into a separate clade within the liberibacter genus and was most closely related to CLaf PTSAPSY (NZ_CP004021.1) (Fig. 2). The concatenated phylogenetic tree presented the same trend in topology as the phylogenetic trees for individual genes, with the *A. solanicola*-associated 'Ca. Liberibacter' species consistently forming a single clade between CLaf and CLas (Fig. S1A–G). All genes had bootstrap support (≥ 97%) for the separation of the *A. solanicola*-associated 'Ca. Liberibacter' species with the exception of *nusG* and *rplA* which produced poor bootstrap support for the separation of CLaf and the *A. solanicola*-associated 'Ca. Liberibacter' species (59% for *nusG* and 32% for *rplA*). However, *nusG* also showed poor bootstrap support for the separation of CLas and CLaf (69%) (Fig. S1D). The bootstrap support significantly increased for the concatenated phylogeny (Fig. 2), showing very strong bootstrap support (100%) for the separation of all clades of bacteria that were analysed.

### Diagnostic quantitative PCR (qPCR) analyses of *A. solanicola*-associated 'Ca. Liberibacter' species

A positive qPCR result was generated from the *A. solanicola*-associated 'Ca. Liberibacter' DNA extracts using

the HLBas and HLBaf detection protocol (Li *et al.*, 2006). A negative qPCR result was generated from *A. solanicola*-associated 'Ca. Liberibacter' DNA using the HLBam (Li *et al.*, 2006) and LsoF (Li *et al.*, 2009) detection protocols. Negative results were generated for all *A. solanicola* DNA extracts that were negative for the *A. solanicola*-associated 'Ca. Liberibacter' species.

The nucleotide sequence alignment of the pathogenic diagnostic forward primers and the 16S rRNA target region on the *A. solanicola*-associated 'Ca. Liberibacter' species reveal six and nine identical base pairs at the 3' end binding site for HLBas and HLBaf respectively (Fig. S2). The reverse primer and probe for these detection protocols are identical.

### Discussion

A new candidate liberibacter species was detected in the native Australian eggplant psyllid, *A. solanicola*, based on the sequence and phylogenetic analyses of the 16S rRNA region and confirmed with MLSA of seven concatenated genes. The proposed name for this *A. solanicola*-associated species is 'Candidatus Liberibacter brunswickensis' (CLbr), after Brunswick, a suburb in Melbourne, Victoria, from where the psyllid host was first collected. CLbr was not associated with plant disease, instead the discovery arose during the screening of Australian native psyllids with newly developed generic liberibacter genus primers.
The generic liberibacter genus primers amplified CLas, CLeu and CLso *in vitro* and 14 sequences representing all known non-pathogenic and phytopathogenic species of liberibacter with sequence data available *in silico*, while excluding closely related bacterial genera (*in silico*). When used with DNA extracted from *A. solanicola*, a single amplicon of the expected size was generated from DNA extracts from 17 of 37 individual psyllids. Direct sequencing of these LG774F-LG1463R amplicons determined that the sequence was most closely related to CLso sequences (99% sequence identity), but was not an identical match. The high sequence identity of this amplicon to CLso highlights the risk of falsely identifying a liberibacter species on regions of the 16S rRNA, as is done in amplicon metagenomics. No additional *A. solanicola* DNA extracts from six sites across south-east Australia produced an amplicon using the liberibacter generic primer pair. These sites were selected based on collection sites previously reported (Kent and Taylor, 2010; Taylor and Kent, 2013). Kent and Taylor (2010) report damage to eggplants caused by *A. solanicola* feeding; however, the microflora of these psyllids was not tested and this may be due to the volume of psyllids feeding.

This study demonstrates that the generic liberibacter primers are a powerful tool to screen hosts for species of liberibacter. To confirm the phylogenetic position of the *A. solanicola*-associated ‘Ca. Liberibacter’ within the liberibacter genus, a larger region of the 16S rRNA was sequenced and the MLSA was conducted. Phylogenetic analysis of the 16S rRNA region separated CLbr (KY077741) into a unique clade, which branched relatively early and prior to the CLas and CLaf clades (Fig. 1). Significant branch bootstrap support (100%) was obtained with phylogenetic analysis of the concatenated genes (Fig. 2) confirming that CLbr separated into a unique clade and that this new species is most closely related to, but distinct from, CLaf PTSAPSY (NZ_CP004021.1). The phylogeny for each individual gene displayed the same topology as the concatenated tree for CLbr (Fig. S1A–G), however bootstrap values varied, supporting the need for the concatenated MLSA mapping approach.

Previous MLSA for unculturable bacteria, including the majority of species within the liberibacter genus, has been limited to cloning and chromosome walking techniques (Planet et al., 1995; Teixeira et al., 2008; Lin et al., 2011). An increase in the availability of next-generation sequencing techniques has allowed a small number of full or near complete (< 10 contigs) genomes of CLas (Duan et al., 2009; Lin et al., 2013; Katoh et al., 2014; Zheng et al., 2014; Wu et al., 2015), CLaf (Lin et al., 2015), CLam (Wulff et al., 2014) and CLso (Lin et al., 2011; Thompson et al., 2015) to be assembled. To date for the candidate liberibacter species, there are

seven complete liberibacter genomes and an extra five near complete genomes publically available. Although this is a small number, these genomes have allowed comparisons to Lcr which can be grown in axenic culture and have provided insights to the fastidious nature of the liberibacter genus (Leonard et al., 2012; Fagen et al., 2014; Wulff et al., 2014). However, it is common for genome sequences in the liberibacter genus to have multiple long repetitive regions and multiple variable prophage inserts, making assembly of a complete genome challenging, even with a reference genome present (Thompson et al., 2015; Wu et al., 2015).

Large complex data sets, such as the *A. solanicola* metagenome sequenced in this study, also present significant challenges to genome assembly. The MLSA mapping approach using known liberibacter reference genes to extract reads for assembly, enabled MLSA to be completed on an undescribed organism within a metagenomic data set. This approach was further validated by spiking the *A. solanicola* metagenomic data set with simulated reads from the complete genome of CLaf PTSAPSY (NZ_CP004021.1) followed by reassembly of the seven genes used for the MLSA analysis for both ‘Ca. Liberibacter’ species. The CLaf PTSAPSY genome was selected for validation as it is phylogenetically the closest species to CLbr and therefore the most likely to generate chimeric contigs during assembly. This method was able to differentiate a novel species of bacteria in the liberibacter genus and also provide assurance that multi-species sequence chimeras had not formed. Individual data sets for both CLbr (metagenomics data from *A. solanicola*, Fig. 3A) and CLaf PTSAPSY (NZ_CP004021.1) (simulated reads, Fig. 3B) assembled into identical contigs to those assembled from the spiked data set (Fig. 3C). This PCR-independent approach has the potential to identify multiple liberibacter species present in the one host, allows species identification without the need to obtain an axenic culture of bacteria and is a useful analytical tool for whole genome shotgun sequence data sets that contain obligate non-culturable organisms. The seven genes selected in this study are highly conserved within the liberibacter genus and were able to differentiate the novel ‘Ca. Liberibacter’ from known species. However, the seven genes are not suitable to separate haplotypes within a liberibacter species (Table 1 and Fig. 2). To apply the MLSA mapping approach to other unculturable bacterial genera, it is recommended to identify genes that are highly conserved within the genus.

The phylogenetic relationship of CLbr to CLaf and CLas may support the proposed theory of the Gondwana supercontinent origins of CLaf and CLas, as Australia was once a part of Gondwana before dislocation approximately 160 million years ago (Teixeira et al., 2008; Nelson et al., 2013; Bové, 2014). In addition, the *Acizzia*
<figure><img src="image_4.png" /><figcaption>**Fig. 3.** Graphical representation of the MLSA approach of raw reads being mapped to seven reference genes and the assembly of these reads for the NextSeq *A. solanicola* metagenome data set (A), CLaf only (B) and CLaf spiked NextSeq *A. solanicola* metagenome data set (C). The *A. solanicola*-associated 'Ca. Liberibacter' species and CLaf species are highlighted in blue and green respectively. All other reads in the metagenome data set are shaded black, and the reference genome CLas is shaded grey.</figcaption></figure>

genus is distributed from Australia, New Zealand, the Old World tropics and extending to North Africa, the Middle East and the Mediterranean (references within Taylor and Kent, 2013). This may infer that liberibacter species have been present within, and evolved with psyllids over time, although as a secondary symbiont this is likely to be a more recent event than that of the evolution of the primary symbiont, 'Ca. Carsonella' (Brumin et al., 2012; Sloan and Moran, 2012). Further research on the biology of the bacterium within its psyllid host is needed to determine the multiple co-speciation events. To date, *A. solanicola* has been detected in Australia, and most recently New Zealand (Taylor and Kent, 2013). Based on sequence analyses, species within the liberibacter genus do not group according to plant species in which they are disease causing or endophytic. This is shown with CLam, pathogenic to citrus, grouping closer to CLeu and Lcr endophytes in pear and papaya (Raddadi et al., 2011; Leonard et al., 2012), supporting the theory that these bacteria may have evolved with their insect hosts.

This is the first time a liberibacter species has been reported in mainland Australia. Without plant disease symptoms present or reported, this discovery arose due to the need to conduct surveillance for liberibacter species of biosecurity significance to Australian agriculture. CLbr gave a false-positive result for the diagnostic qPCR methods reported in (Li et al., 2006) for CLas and CLaf but not for CLso and CLam (Li et al., 2006, 2009). This shows that there is microflora in native Australian psyllids that is able to confound diagnostic tests developed

outside of Australia. Australia is a centre for psyllid diversity, and *A. solanicola* was the first psyllid species to be screened with the generic liberibacter genus primers. The detection of CLbr is an indication that the number of species that have been identified in the liberibacter genus is a huge underrepresentation and the genus will increase beyond the current eight recorded species. Further work into the biology, host range and distribution of CLbr is underway.

## Experimental procedures

### *Insect collection and morphological identification*

A colony of *A. solanicola* was collected from an eggplant, *Solanum melongena* Black Beauty in Brunswick (a metropolitan suburb of Melbourne), Victoria (VIC), Australia. No disease symptoms were seen on the eggplant, nor on the cape gooseberry, *Physalis peruviana* (cultivar unknown) on which an additional colony of *A. solanicola* was also established at the collection site. The eggplant and psyllid colony was transferred to the laboratory and kept in insect cages (BugDorm, Taiwan) at 20°C (± 2°C) and 60% (± 5%) relative humidity. Individuals from the psyllid colony were morphologically confirmed as *A. solanicola* (AgriBio, 5 Ring Road, Victoria 3083, Australia).

An additional 31 *A. solanicola* individuals were collected from *S. melongena* or *S. mauritianum* across six sites in south-eastern Australia. These sites include Bunchoora (VIC, nine psyllids), Kew (VIC, two psyllids), Bellingen (NSW, nine psyllids), Clybucca (NSW, five
psyllids), Dorrigo West (NSW, three psyllids) and Karangi (NSW, three psyllids).

## DNA extraction

DNA was extracted using the Qiagen DNeasy Blood and Tissue Kit (Qiagen, Germany) with the following amendments. Individual psyllids were placed into a 1.5-mL tube (Eppendorf, Germany) containing two 3 mm glass beads and 20 µl of proteinase K (Qiagen, Germany). Insects were homogenized using a Bead-Mill TissueLyser (Qiagen, Germany) at 30 megahertz (Mhz) for 1 min and were centrifuged at 17 000× gravitational force (g) for 20 s (s). The slurry was re-homogenized at 30 Mhz for 1 min, centrifuged a second time for 20 s, 180 µl of ATL buffer (Qiagen, Germany) added and the protocol completed according to manufacturer's instructions (Qiagen, Germany) with a final elution volume of 100 µl.

To allow an increased concentration of DNA to be extracted from an individual psyllid for Illumina sequencing, the cetyl trimethyl ammonium bromide (CTAB) method (Murray and Thompson, 1980) was used with the following amendments. Individual psyllids were placed into a 1.5-mL tube (Eppendorf, Germany) containing two 3 mm glass beads, 50 µl of CTAB buffer (2% CTAB, 1.4 M NaCl, 1% PVP-40, 0.02 M EDTA and 0.1 M Tris-HCl, pH 8.0) and insects were homogenized as described above, using a Bead-Mill TissueLyser (Qiagen, Germany). After precipitation, DNA was resuspended in 15 µl of UltraPure distilled water (Invitrogen, Netherlands).

DNA quantity and quality was estimated from each extraction using the Nanodrop ND-1000 Spectophotometer (Thermo Scientific, Germany) and Qubit 2.0 fluorometer (Invitrogen, Germany) and then stored at -20°C.

## Development of generic Liberibacter genus primers

To identify microflora in Australian psyllids that is phylogenetically close to the liberibacter genus, generic liberibacter genus-specific primers were designed to amplify sequences from known liberibacter species encompassing both non-pathogenic and pathogenic species, while excluding phylogenetically related bacterial genera. CLca was not able to be included in this analysis due to insufficient sequence data available (currently 1125 bp of the 16S rRNA region for both GenBank additions: KP012550.1 and KP012551.1). Full-length nucleotide sequences of the 16S rRNA region of the remaining six liberibacter species and species of phylogenetically related alpha Proteobacteria (Table 2 – Sequence alignment of the 16S rRNA region) were aligned using MUSCLE (Edgar, 2004) in MEGA6 (Tamura et al., 2013). PCR primers were designed using Primer3plus (http://primer3plus.com/) to target regions identified in the alignment that were conserved between liberibacter species, while excluding phylogenetically related alpha Proteobacteria.

Forward and reverse primers, LG774F (5'-GTAAAC GATGAGTGCTAGCTGTTGGG-3') and LG1463R (5'-CTGACCRTACCGTGGCCGG-3'), were selected to

<table><caption>**Table 2.** Bacteria used for sequence alignments and phylogenetic analyses.</caption><thead><tr><th rowspan="2">Bacterium species and strain/isolate</th><th rowspan="2">Accession number</th><th colspan="2">Sequence alignment and phylogenetic analysis</th><th rowspan="2">Class: Order</th></tr><tr><th>16S rRNA region</th><th>MLSA</th></tr></thead><tbody><tr><td>CLaf PTSAPSY</td><td>NZ_CP004021.1</td><td>●</td><td>●</td><td>Alpha Proteobacteria: Rhizobiales</td></tr><tr><td>CLam Sao Paulo</td><td>NC_022793.1</td><td>●</td><td>●</td><td>Alpha Proteobacteria: Rhizobiales</td></tr><tr><td>CLas gxpsy</td><td>NC_020549.1</td><td>●</td><td>●</td><td>Alpha Proteobacteria: Rhizobiales</td></tr><tr><td>CLas Ishi-1</td><td>NZ_AP014595.1</td><td>●</td><td>●</td><td>Alpha Proteobacteria: Rhizobiales</td></tr><tr><td>CLas A4</td><td>NZ_CP010804.1</td><td>●</td><td>●</td><td>Alpha Proteobacteria: Rhizobiales</td></tr><tr><td>CLas psy62</td><td>NC_012985.3</td><td>●</td><td>●</td><td>Alpha Proteobacteria: Rhizobiales</td></tr><tr><td>CLso CLso-ZC1</td><td>NC_014774.1</td><td>●</td><td>●</td><td>Alpha Proteobacteria: Rhizobiales</td></tr><tr><td>Lcr str. BT-1</td><td>NC_019907.1</td><td>●</td><td>●</td><td>Alpha Proteobacteria: Rhizobiales</td></tr><tr><td>Bartonella henselae str. Hourstron-1</td><td>NC_005956.1</td><td>●</td><td>●</td><td>Alpha Proteobacteria: Rhizobiales</td></tr><tr><td>Bartonella australis str. Aust/NH1</td><td>NC_020300.1</td><td>●</td><td>●</td><td>Alpha Proteobacteria: Rhizobiales</td></tr><tr><td>Agrobacterium tumefaciens str. Arch5</td><td>NZ_CP011246.1</td><td>●</td><td>●</td><td>Alpha Proteobacteria: Rhizobiales</td></tr><tr><td>Agrobacterium radiobacter str. K84</td><td>NC_011985.1</td><td>●</td><td>●</td><td>Alpha Proteobacteria: Rhizobiales</td></tr><tr><td>Sinorhizobium meliloti str. 1021</td><td>NC_003047.1</td><td>●</td><td>●</td><td>Alpha Proteobacteria: Rhizobiales</td></tr><tr><td>Sinorhizobium medicae str. WSM419</td><td>NC_009636.1</td><td>●</td><td>●</td><td>Alpha Proteobacteria: Rhizobiales</td></tr><tr><td>Rhizobium etli str. CFN42</td><td>NC_007761.1</td><td>●</td><td>●</td><td>Alpha Proteobacteria: Rhizobiales</td></tr><tr><td>CLaf Mpumalanga-UPCRI-06-0071</td><td>EU754741.1</td><td>●</td><td>●</td><td>Alpha Proteobacteria: Rhizobiales</td></tr><tr><td>CLaf str. LEI16SRRNB</td><td>L22533.1</td><td>●</td><td>●</td><td>Alpha Proteobacteria: Rhizobiales</td></tr><tr><td>CLso str. NZ083338</td><td>EU935004.1</td><td>●</td><td>●</td><td>Alpha Proteobacteria: Rhizobiales</td></tr><tr><td>CLso str. NZ082226</td><td>EU834130.1</td><td>●</td><td>●</td><td>Alpha Proteobacteria: Rhizobiales</td></tr><tr><td>CLps PRR1</td><td>EU812559.1</td><td>●</td><td>●</td><td>Alpha Proteobacteria: Rhizobiales</td></tr><tr><td>CLeu 94B</td><td>JX244260.1</td><td>●</td><td>●</td><td>Alpha Proteobacteria: Rhizobiales</td></tr><tr><td>Escherichia coli O157:H7 str. Sakai</td><td>NC_002695.1</td><td>●</td><td>●</td><td>Gamma Proteobacteria: Enterobacteriales</td></tr></tbody></table>
produce an amplicon of 684 bp and tested against DNA extracts of CLas, CLeu and CLso. DNA extractions of 37 individual *A. solanicola* psyllids isolated from the colony at different times were tested using these liberibacter generic primers. Reactions were performed in 25 µl containing 1 µl of DNA template, 10× PCR buffer (Invitrogen, Netherlands), 1.5 mM MgCl₂, 0.2 mM dNTPs, 0.25 µM of each primer and 1 U of Platinum Taq polymerase (Invitrogen, Netherlands). The PCR conditions were denaturation at 95°C for 3 mins, then 35 cycles of 95°C for 30 s, 62°C for 30 s and 72°C for 1 min, plus a final step at 72°C for 5 mins. A second PCR using the primer pair FD2 (5'-CCGAATTCGTCGACAACAGAGTTTGATCATGGCT CAG-3') and RP1 (5'-CCCGGGATCCAAGCTTACGGT TACCTTGTTACGACTT-3') was performed (Weisburg et al., 1991) to amplify the 16S rRNA region of all bacterial species present within the *A. solanicola* samples to ensure that DNA extraction had been successful. Amplification was performed in an Eppendorf Mastercycler pro PCR System (Fisher Scientific, United States of America), and PCR products were visualized by agarose gel electrophoresis (1.5% agarose gel) with SYBR Safe DNA gel stain (Invitrogen, Netherlands).

Amplicons were purified using the Wizard SV gel and PCR Clean-Up System (Promega, United States of America) according to manufacturer's instructions. Sequence of the purified PCR products was generated via capillary Sanger sequencing in both directions, performed by the Australian Genome Research Facility (AGRF), Melbourne, Australia. Consensus sequences were generated using Geneious 7.0 (Kearse et al., 2012) and used as queries against the GenBank nr database using BLASTn version 2.4.0 (Altschul et al., 1990) for homologous sequences, available at the National Centre for Biotechnology Information (NCBI).

## Cloning and sequencing of the 16S rRNA region of liberibacter-like bacteria within *A. solanicola*

The FD2-RP1 product of the 16S rRNA region (approximately 1500 bp) was cloned from *A. solanicola* DNA extracts that tested positive using the generic liberibacter genus primers LG774F-LG1463R. The amplicon generated using the generic bacterial primers for 16S rRNA region amplification, FD2-RP1, was cloned into the pGEM-T Easy vector (Promega, United States of America) following manufacturer's instructions. To obtain colonies with inserts closely related to the 16S rRNA region of liberibacter species, DNA from isolated colonies was screened by PCR using the liberibacter generic primers LG774F-LG1463R. Plasmid DNA from LG774F-LG1463R positive FD2-RP1 colonies was purified by centrifugation using the PureLink Quick Plasmid Miniprep Kit (Invitrogen, Netherlands) according to

manufacturer's instructions and sequenced via capillary Sanger sequencing (AGRF, Melbourne, Vic., Australia). Consensus sequences were analysed as described above.

## Library preparation and Illumina sequencing of *A. solanicola*

Illumina sequencing was performed on a single DNA extract of *A. solanicola* that tested positive using the liberibacter generic, LG774F-LG1463R, primer pair. The quality and quantity of the *A. solanicola* DNA was estimated using the Nanodrop ND-1000 Spectophotometer (Thermo Scientific) and double-stranded DNA (dsDNA) HS assay for the Qubit 2.0 fluorometer (Invitrogen). DNA was diluted to 0.2 ng µl⁻¹, and a library was prepared using the Nextera XT DNA Library Preparation kit (Illumina) following the manufacturer's instructions. Library quality was checked using the Qubit 2.0 fluorometer (Invitrogen, Netherlands), KAPA Library Quantification Kit (KAPA Biosystems, United States of America), and the 2200 TapeStation system (Agilent, United States of America) according to the High Sensitivity D1000 protocol. Initial paired-end sequencing (2 × 250 bp) was performed on an Illumina MiSeq (Reagent Kit v3) to establish library quality, optimal library concentration and estimate the required sequencing depth. A final paired-end sequencing run (2 x 150 bp) was performed on an Illumina NextSeq (Reagent Kit v2).

## In silico multilocus sequence analysis (MLSA)

Seven highly conserved genes were selected for MLSA of the unidentified 'Ca. Liberibacter'. These included DNA primase (dnaG), DNA mismatch repair protein (mutS) (Glynn et al., 2012), DNA gyrase subunit B (gyrB) (Marrero et al., 2013; Mee et al., 2015), Transcription termination/antitermination protein NusG (nusG), 50S ribosomal protein L1 (rpL1), DNA-directed RNA polymerase subunit beta (rpoB) and a single copy of Elongation factor Tu (tufB) (Teixeira et al., 2008).

Sequence reads from the NextSeq *A. solanicola* meta-genome data set were adapter and quality trimmed using TrimGalore (Krueger, 2012), and mapped to the seven genes listed above from the complete genome of CLas gxpsy (NC_020549.1) using BWA-MEM (Li, 2013) with relaxed parameters for a mismatch penalty (reduced from 4 to 2). Reads that mapped were assembled using SPAdes de novo version 3.6.0 (Bankevich et al., 2012) with kmer lengths of 107, 87, 71 and 31 (Fig. 3A). Resulting contigs were compared to known sequences in GenBank using the BLASTn algorithm (Altschul et al., 1990) to identify the 'Ca. Liberibacter' sequences and to confirm the gene identity.
To establish that one species of bacteria can be iden-
tified from a large metagenomic sequence data set using
the above MLSA mapping approach, and that the result-
ing contigs are not merely a chimera of bacterial reads
within the sample, one million simulated sequence reads
(150 bp in length) were generated using wgsim version
0.3.0 (Li, 2011) from the most closely related species to
the *A. solanicola*-associated ‘Ca. Liberibacter’ (based on
16S rRNA phylogenetic analysis), CLaf PTSAPSY (NZ_
CP004021.1). The simulated reads were mapped directly
to the seven reference genes of CLas gxpsy (NC_
020549.1) to ensure that the genes re-assembled from
the simulated reads (Fig. 3B). The simulated reads were
then spiked into the NextSeq *A. solanicola* metagenome
data set, and all reads were mapped to the seven refer-
ence loci listed above (Fig. 3C). For both mapping exer-
cises, the MLSA extraction, assembly and identification
were repeated as described above.

Phylogenetic analysis

Cloned 16S rRNA region sequence, the assembled con-
tigs of dnaG, gyrB, mutS, nusG, rplA, rpoB and tufB of
the new ‘Ca. Liberibacter’ from *A. solanicola* and the
concatenated multilocus sequence were aligned with
sequences of the liberibacter genus and sequences of
closely related bacteria using MUSCLE (Edgar, 2004)
according to Table 1. The concatenated liberibacter
genus sequences were screened for presence of recom-
bination with Recombination Detection Program 4.56
(RDP4) (Martin *et al.*, 2015) using the default parame-
ters for RDP, GENECONV, Chimaera, MaxChi, BOOT-
SCAN and SISCAN methods.

Phylogenetic analyses were performed for the 16S
rRNA region, the seven individual genes (Fig. S1) and
the concatenated MLSA using the program MEGA6 (Tamura *et al.*, 2013) employing the maximum likeli-
hood method. A bootstrap consensus tree was inferred
using 1000 samplings of the data, and bootstrap val-
ues below 70% were not shown. For each alignment,
24 different substitution models were compared using
the Find Best DNA Model Test in MEGA6. The best
substitution model was selected by the lowest BIC
score (Bayesian Information Criterion) and applied to
each phylogenetic analyses accordingly (Kimura, 1980;
Lanave *et al.*, 1984; Nei and Kumar, 2000; Zwickl and
Holder, 2004). Nearest-Neighbor-Interchange (NNI)
algorithm (Nei and Kumar, 2000), Neighbor Joining/Bio
Neighbor Joining (Saitou and Nei, 1987) and very
strong branch swap filter settings were used to obtain
the initial tree. *Escherichia coli* O157:H7 str. Sakai was
used as an out-group and a root placed on this branch
for the MLSA phylogenetic analyses (including individ-
ual gene trees, supporting information). In addition,

a distance matrix of the concatenated genes was
generated.

Diagnostic qPCR analyses of *A. solanicola*-associated
*Ca. Liberibacter* species

Two *A. solanicola* DNA extracts that were positive for the *A. solanicola*-associated ‘*Ca. Liberibacter*’ species and two *A. solanicola* DNA extracts that were negative for the *A. solanicola*-associated ‘*Ca. Liberibacter*’ species were tested using diagnostic qPCR protocols for CLas, CLaf, CLam and CLso (Li *et al.*, 2006, 2009) according to authors instructions using the QuantStudio 3 Real-Time PCR System and the QuantiTect Multiplex PCR NoROX Kit (Qiagen, Germany). To assess primer binding regions, the phytopathogenic ‘*Ca. Liberibacter*’ species-specific primers and probe were aligned to the target 16S rRNA region of *A. solanicola*-associated ‘*Ca. Liberibacter*’ using Mega6.

Acknowledgements

We would like to acknowledge the support of the Aus-
tralian Government's Cooperative Research Centres
Program; Annie Raser-Rowland for providing the initial
colony of *Acizzia solanicola* and eggplant for subse-
quent colony maintenance; Dr. Mark Blacket, Dr. Linda
Semeraro and Dr. Lia Liefting for advice on psyllid
extraction; Dr. Rebekah Frampton for providing 'Ca.
Liberibacter europaeus' and 'Ca. Liberibacter solana-
cearum' samples; Dr Baishi Hu for providing 'Ca.
Liberibacter asiaticus' samples; Dr Deborah Kent and
Dr Isabel Valenzuela-Gonzalez for guidance and sup-
port while collecting *Acizzia solanicola* in New South
Wales.

Conflict of Interest

None declared.

References

Altschul, S.F., Gish, W., Miller, W., Myers, E.W., and Lipman, D.J. (1990) Basic local alignment search tool. *J Mol Biol* **215**: 403–410.

Bankevich, A., Nurk, S., Antipov, D., Gurevich, A.A., Dvorkin, M., Kulikov, A.S., et al. (2012) SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. *J Comput Biol* **19**: 455–477.

Bové, J.M. (2006) Invited review Huanglongbing: a destruc-
tive, newly-emerging, century-old disease of citrus. *J
Plant Pathol* **88**: 7–37.

Bové, J.M. (2014) Huanglongbing or yellow shoot, a disease of Gondwanan origin: will it destroy citrus worldwide?
*Phytoparasitica* **42**: 579–583.
Brumin, M., Levy, M., and Ghanim, M. (2012) Transovarial transmission of *Rickettsia* spp. and organ-specific infection of the whitefly *Bemisia tabaci*. *Appl Environ Microbiol* **78**: 5565–5574.

Duan, Y., Zhou, L., Hall, D.G., Li, W., Doddapaneni, H., Lin, H., et al. (2009) Complete genome sequence of citrus huanglongbing bacterium, "*Candidatus* Liberibacter asiaticus" obtained through metagenomics. *Mol Plant-Microbe Interact* **22**: 1011–1020.

Edgar, R.C. (2004) MUSCLE: multiple sequence alignment with high accuracy and high throughput. *Nucleic Acids Research* **32**: 1792–1797.

Fagen, J.R., Leonard, M.T., McCullough, C.M., Edirisinghe, J.N., Henry, C.S., Davis, M.J., and Triplett, E.W. (2014) Comparative genomics of cultured and uncultured strains suggests genes essential for free-living growth of Liberibacter. *PLoS ONE* **9**: 1–11.

Glynn, J.M., Islam, M.S., Bai, Y., Lan, S., Wen, A., Gudmes-tad, N.C., et al. (2012) Multilocus sequence typing of "Candidatus Liberibacter solanacearum" isolates from North America and New Zealand. *J Plant Pathol* **94**: 223–228.

Hocquellet, A., Toorawa, P., Bové, J.-M., and Garnier, M. (1999) Detection and identification of the two *Candidatus* Liberobacter species associated with citrus huanglongbing by PCR amplification of ribosomal protein genes of the β operon. *Mol Cell Probes* **13**: 373–379.

Jagoueix, S., Bove, J.M., and Garnier, M. (1994) The phloem-limited bacterium of greening disease of citrus is a member of the alpha subdivision of the Proteobacteria. *Int J Syst Bacteriol* **44**: 379–386.

Katoh, H., Miyata, S., Inoue, H. and Iwanami, T. (2014) Unique features of a Japanese "Candidatus Liberibacter asiaticus" Strain revealed by whole genome sequencing. *PLoS ONE* **9**: e106109.

Kearse, M., Moir, R., Wilson, A., Stones-Havas, S., Cheung, M., Sturrock, S., et al. (2012) Geneious Basic: an integrated and extendable desktop software platform for the organization and analysis of sequence data. *Bioinformatics* (Oxford, England) **28**: 1647–1649.

Kent, D., and Taylor, G. (2010) Two new species of *Acizzia* Crawford (Hemiptera: Psyllidae) from the Solanaceae with a potential new economic pest of eggplant, *Solanum melongena*. *Aust J Entomol* **49**: 73–81.

Keremane, M.L., Ramadugu, C., Castaneda, A.Diaz., J. E., P., Chen, E.A., Duan, Y. P., Halbert, S. E. and Lee, R. F. (2015) Report of *Candidatus* Liberibacter caribbeanus, a new citrus- and psyllid-associated *Liberibacter* from Colombia, South America. In *American Phytopathological Society Annual Meeting*. URL http://www.apsnet.org/meetings/Documents/2015_meeting_abSTRACTS/aps2015ab0253.htm

Kimura, M. (1980) A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. *J Mol Evol* **16**: 111–120.

Krueger, F. (2012) Trim Galore! URL https://github.com/FelixKrueger/TrimGalore/blob/master/Trim_Galore_User_Guide.pdf

Lanave, C., Preparata, G., Sacone, C., and Serio, G. (1984) A new method for calculating evolutionary substitution rates. *J Mol Evol* **20**: 86–93.

Leonard, M.T., Fagen, J.R., Davis-Richardson, A.G., Davis, M.J., and Triplett, E.W. (2012) Complete genome

sequence of *Liberibacter crescens* BT-1. *Stand Genomic Sci* **7**: 271–283.

Li, H. (2011) Wgsim: Reads simulator. URL https://github.com/lh3/wgsim

Li, H. (2013) Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM. URL http://arxiv.org/abs/1303.3997

Li, W., Hartung, J.S., and Levy, L. (2006) Quantitative real-time PCR for detection and identification of *Candidatus* Liberibacter species associated with citrus huanglongbing. *J Microbiol Methods* **66**: 104–115.

Li, W., Abad, J.a., French-Monar, R.D., Rascoe, J., Wen, A. and Gudmestad, et al. (2009) Multiplex real-time PCR for detection, identification and quantification of "*Candidatus* Liberibacter solanacearum" in potato plants with zebra chip. *J Microbiol Methods* **78**: 59–65.

Liefting, L.W., Sutherland, P.W., Ward, L.I., Paice, K.L., Weir, B.S., and Clover, G.R.G. (2009a) A new "Candidatus Liberibacter" species associated with diseases of Solanaceous crops. *Plant Dis* **93**: 208–214.

Liefting, L.W., Weir, B.S., Pennycook, S.R., and Clover, G.R.G. (2009b) "Candidatus Liberibacter solanacearum", associated with plants in the family Solanaceae. *Int J Syst Evol Microbiol* **59**: 2274–2276.

Lin, H., Lou, B., Glynn, J.M., Doddapaneni, H., Civerolo, E.L., Chen, C., et al. (2011) The complete genome sequence of "Candidatus Liberibacter solanacearum", the bacterium associated with potato Zebra Chip disease. *PLoS ONE* **6**: 1–13.

Lin, H., Han, C.S., Liu, B., Lou, B., Bai, X., Deng, C., et al. (2013) Complete genome sequence of a chinese strain of "Candidatus Liberibacter asiaticus". *Genome Announc* **12**: 1–13.

Lin, H., Pietersen, G., Han, C., Read, D.A., Lou, B., Gupta, G. and Civerolo, E.L. (2015) Complete genome sequence of "Candidatus Liberibacter africanus," a bacterium associated with citrus Huanglongbing. *Genome Announc* **3**: 1–2.

Marrero, G., Schneider, K.L., Jenkins, D.M., and Alvarez, A.M. (2013) Phylogeny and classification of *Dickeya* based on multilocus sequence analysis. *Int J Syst Evol Microbiol* **63**: 3524–3539.

Martin, D.P., Murrell, B., Golden, M., Khoosal, A., and Muhire, B. (2015) RDP4: detection and analysis of recombi nation patterns in virus genomes. *Virus Evol* **1**: 1–5.

Mee, P.T., Weeks, A.R., Walker, P.J., Hoffmann, A.A., and Duchemin, J.-B. (2015) Detection of low-level *Cardinium* and *Wolbachia* infections in *Culicoides*. *Appl Environ Microbiol* **81**: 6177–6188.

Morgan, J.K., Zhou, L., Li, W., Shatters, R.G., Keremane, M., and Duan, Y.-P. (2012) 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* **26**: 90–98.

Munyaneza, J.E., Fisher, T.W., Sengoda, V.G., Garczynski, S.F., Nissinen, A., and Lemmetty, A. (2010) (Hemiptera: Triozidae) in Europe. *J Econ Entomol* **103**: 1060–1070.

Munyaneza, J.E., Sengoda, V.G., Sundheim, L., and Meadow, R. (2014) Survey of "Candidatus Liberibacter solanacearum" in carrot crops affected by the psyllid *Triozoa apicalis* (Hemiptera: Triozidae) in Norway. *J Plant Pathol* **96**: 397–402.
Murray, M.G., and Thompson, W.F. (1980) Rapid isolation of high molecular weight plant DNA. *Nucleic Acids Res* **8**: 4321–4325.

Nei, M., and Kumar, S. (2000) *Molecular evolution and phylogenetics*. Oxford University Press.

Nelson, W.R., Munyaneza, J.E., McCue, K.F., and Bové, J.M. (2013) The Pangaean origin of "Candidatus Liberibacter" species. *J Plant Pathol* **95**: 455–461.

Phahladira, M.N.B., Viljoen, R., and Pietersen, G. (2012) Widespread occurrence of "Candidatus Liberibacter africanus" subspecies capensis" in *Calodendrum capense* in South Africa. *Eur J Plant Pathol* **134**: 39–47.

Planet, P., Jagoueix, S., Bové, J.M., and Garnier, M. (1995) Detection and characterization of the African citrus greening Liberobacter by amplification, cloning, and sequencing of the *rplKAJL-rpoBC* operon. *Curr Microbiol* **30**: 137–141.

Raddadi, N., Gonella, E., Camerota, C., Pizzinat, A., Tedeschi, R., Crotti, E., et al. (2011) "Candidatus Liberibacter europaeus" sp. nov. that is associated with and transmitted by the psyllid *Cacopsylla pyri* apparently behaves as an endophyte rather than a pathogen. *Envi-ron Microbiol* **13**: 414–426.

Saitou, N., and Nei, M. (1987) The neighbor-joining method: a new method for reconstructing phylogenetic trees. *Mol Biol Evol* **4**: 406–425.

Sloan, D.B. and Moran, N.a. (2012) Genome reduction and co-evolution between the primary and secondary bacterial symbionts of psyllids. *Mol Biol Evol* **29**: 3781–3792.

Tahzima, R., Maes, M., Achbani, E.H., Swisher, K.D., Munyaneza, J.E., and De Jonghe, K. (2014) First report of "Candidatus Liberibacter solanacearum" on carrot in Africa. *Plant Dis* **98**: 1426.

Tamura, K., Stecher, G., Peterson, D., Filipski, A., and Kumar, S. (2013) MEGA6: Molecular Evolutionary Genetics Analysis version 6.0. *Mol Biol Evol* **30**: 2725–2729.

Taylor, G.S. and Kent, D.S. (2013) Potential economic pests of solanaceous crops: a new species of Solanum-feeding psyllid from Australia and first record from New Zealand of *Acizzia solanicola* (Hemiptera: Psyllidae). *Zootaxa* **3613**: 257–273.

Teixeira, D.C., Eveillard, S., Sirand-Pugnet, P., Wulff, A., Saillard, C., Ayres, A.J., and Bové, J.M. (2008) The *tufB-secE-nusG-rplKAJL-rpoB* gene cluster of the liberibacteria: sequence comparisons, phylogeny and speciation. *Int J Syst Evol Microbiol* **58**: 1414–1421.

Thompson, S.M., Johnson, C.P., Lu, A.Y., Frampton, R.A., Sullivan, K.L., Fiers, M.W.E.J., et al. (2015) Genomes of "Candidatus Liberibacter solanacearum" haplotype A from

New Zealand and the United States suggest significant genome plasticity in the species. *Phytopathology* **105**: 863–871.

Weisburg, W.G., Barns, S.M., Pelletie, D.A. and Lane, D.J. (1991) 16S ribosomal DNA amplification for phylogenetic study. *J Bacteriol* **173**: 697–703.

Wen, A., Mallik, I. and Alvarado, et al. (2009) Detection, distribution, and genetic variability of "Candidatus Liberibacter" species associated with Zebra Complex disease of potato in North America. *Plant Dis* **93**: 1102–1115.

Wen, A., Johnson, C., and Gudmestad, N.C. (2013) Development of a PCR assay for the rapid detection and differentiation of "Candidatus Liberibacter solanacearum" haplotypes and their spatiotemporal distribution in the United States. *Am J Potato Res* **90**: 229–236.

Wu, F., Zheng, Z., Deng, X., Cen, Y., Liang, G. and Chen, J. (2015) Draft genome sequence of "Candidatus Liberibacter asiaticus" from *Diaphorina citri* in Guangdong, China. *Genome Announc* **3**: 1–2.

Wulff, N.A., Zhang, S., Setubal, J.C., Almeida, N.F., Martins, E.C., Harakava, R., et al. (2014) The complete genome sequence of "Candidatus Liberibacter americanus", associated with Citrus huanglongbing. *Mol Plant-Microbe Interact* **27**: 163–176.

Yen, A.L. (2002) Short-range endemism and Australian Psylloidea (Insecta:Hemiptera) in the genera Glycaspis and Acizzia (Psyllidae). *Invertebr Syst* **16**: 631–636.

Zheng, Z., Deng, X. and Chen, J. (2014) Whole-genome sequence of "Candidatus Liberibacter asiaticus" from Guangdong, China. *Genome Announc* **2**: 1–2.

Zwickl, D., and Holder, M. (2004) Model parameterization, prior distributions, and the general time-reversible model in Bayesian phylogenetics. *Syst Biol* **53**: 877–888.

**Supporting information**

Additional Supporting Information may be found online in the supporting information tab for this article:

**Fig. S1.** Phylogenetic analysis for each gene in the MLSA.

**Fig. S2.** Sequence alignment for the 16S RNA region for the primer binding regions of the phytopathogenic 'Ca. Liberibacter' species-specific forward primers (HLBas, HLBaf, HLBam and LsoF) to the *A. solanicola*-associated 'Ca. Liberibacter' species.