RESEARCH ARTICLE

# A diagnostic real-time PCR assay for the rapid identification of the tomato-potato psyllid, *Bactericera cockerelli* (Šulc, 1909) and development of a psyllid barcoding database

J. C. Sumner-Kalkun<sup>1*</sup>, M. J. Sjölund<sup>1</sup>, Y. M. Arnsdorf<sup>1</sup>, M. Carnegie<sup>1</sup>, F. Highet<sup>1</sup>, D. Ouvrard<sup>2,3</sup>, A. F. C. Greenslade<sup>4</sup>, J. R. Bell<sup>4</sup>, R. Sigvald<sup>5</sup>, D. M. Kenyon<sup>1</sup>

1. SASA, Edinburgh, United Kingdom, 2. Department of Life Sciences, Natural History Museum, London, United Kingdom, 3. Entomology and invasive plants Unit, Plant Health Laboratory, ANSES, Montferrier-sur-Lez Cedex, France, 4. Rothamsted Insect Survey, Rothamsted Research, Harpenden, Hertfordshire, United Kingdom, 5. Department of Ecology, Swedish University of Agricultural Sciences, Uppsala, Sweden

* jason.sumner-kalkun@sasa.gov.scot

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

## Abstract

The accurate and rapid identification of insect pests is an important step in the prevention and control of outbreaks in areas that are otherwise pest free. The potato-tomato psyllid *Bactericera cockerelli* (Šulc, 1909) is the main vector of ‘Candidatus Liberibacter solanacearum’ on potato and tomato crops in North America and New Zealand; and is considered a threat for introduction in Europe and other pest-free regions. This study describes the design and validation of the first species-specific TaqMan probe-based real-time PCR assay, targeting the ITS2 gene region of *B. cockerelli*. The assay detected *B. cockerelli* genomic DNA from adults, immatures, and eggs, with 100% accuracy. This assay also detected DNA from cloned plasmids containing the ITS2 region of *B. cockerelli* with 100% accuracy. The assay showed 0% false positives when tested on genomic and cloned DNA from 73 other psyllid species collected from across Europe, New Zealand, Mexico and the USA. This included 8 other species in the *Bactericera* genus and the main vectors of ‘Candidatus Liberibacter solanacearum’ worldwide. The limit of detection for this assay at optimum conditions was 0.000001 ng DNA (~200 copies) of ITS2 DNA which equates to around a 1:10000 dilution of DNA from one single adult specimen. This assay is the first real-time PCR based method for accurate, robust, sensitive and specific identification of *B. cockerelli* from all life stages. It can be used as a surveillance and monitoring tool to further study this important crop pest and to aid the prevention of outbreaks, or to prevent their spread after establishment in new areas.

## 1. Introduction

The psyllid *Bactericera cockerelli* (Šulc, 1909), (commonly known as “Potato Psyllids” or “Tomato-Potato Psyllid”), is a major pest of cultivated Solanaceous crops including potato and tomato [1]. Feeding by this psyllid causes severe damage to potato plants including: deformed tubers; production of numerous small, poor quality tubers; curling of leaves and petioles; and

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

OPEN ACCESS

**Citation:** Sumner-Kalkun JC, Sjölund MJ, Arnsdorf YM, Carnegie M, Highet F, Ouvrard D, et al. (2020) A diagnostic real-time PCR assay for the rapid identification of the tomato-potato psyllid, *Bactericera cockerelli* (Šulc, 1909) and development of a psyllid barcoding database. PLoS ONE 15(3): e0230741. https://doi.org/10.1371/journal.pone.0230741

**Editor:** Sean Michael Prager, University of Saskatchewan College of Agriculture and Bioresources, CANADA

**Received:** December 18, 2019

**Accepted:** March 6, 2020

**Published:** March 26, 2020

**Peer Review History:** PLOS recognizes the benefits of transparency in the peer review process; therefore, we enable the publication of all of the content of peer review and author responses alongside final, published articles. The editorial history of this article is available here: https://doi.org/10.1371/journal.pone.0230741

**Copyright:** © 2020 Sumner-Kalkun et al. This is an open access article distributed under the terms of the [Creative Commons Attribution License](https://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
**Data Availability Statement:** All Co1 and ITS2 sequences area available from GenBank accession numbers (MT021761-MT021824; MT027551-MT027599; MT038907-MT038996; MT040955-MT040966). These will be made accessible on request.

**Funding:** DK: This work was supported by the EU Horizon2020 Programme under grant agreement No. 635646, POnTE (Pest Organisms Threatening Europe) https://ec.europa.eu/programmes/horizon2020/en and the Scottish Government [RRL/001/14]. JB: The Rothamsted Insect Survey, a National Capability, is funded by the Biotechnology and Biological Sciences Research Council under the Core Capability Grant BBS/E/C/ 000J0200. https://bbsrc.ukri.org/funding/The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

**Competing interests:** The authors have declared that no competing interests exist.

yellowing or purpling of leaves. This leads to stunted growth and loss of yield [2]. *Bactericera cockerelli* is also the main vector of ‘*Candidatus* Liberibacter solanacearum’ (Lso) which is associated with Zebra Chip in Central and North America and New Zealand [3–8]. *Bactericera cockerelli* is thought to originate from South-Western USA and Mexico [2,9] and from here has spread via natural and human-mediated dispersal to extend its range [10]. Outside America it is now established in New Zealand [11] and more recently Western Australia [12].

While *B. cockerelli* prefers to complete its life cycle on Solanaceous plants it can also complete development on species of Convolvulaceae (Bindweeds and Morning Glories) [13]. In addition, adult *B. cockerelli* have been found on over 40 species belonging to 20 families, however most of these are either casual, food or shelter plants on which the psyllid is unable to complete a full life cycle [2,9,14–19]. Four biotypes of *B. cockerelli* have been described according to polymorphisms in the mitochondrial cytochrome *c oxidase subunit I* (COI) gene and represent geographically distinct populations; central, western, north-western, and south-western [20,21]. Evidence suggest that these genetic types may differ in their ability to spread Lso [21,22].

The phloem-limited bacterium ‘*Candidatus* Liberibacter solanacearum’ (Lso) is a pathogen associated with Zebra Chip disease of potatoes [3,23–25] and disease in other Solanaceous crops such as cultivated tomato [1,3,26,27], pepper [28], eggplant [29], tobacco [30,31] and tomatillo [26]. Currently, *B. cockerelli* is the main vector of Lso in field and glasshouse-grown Solanaceous plants in the United States, Mexico, areas of Central America [27–30], Canada [32], New Zealand [5,6,25] and recently Ecuador [33]. Ten Lso haplotypes have been described, only three of which are associated with disease in Solanaceous plants. Haplotypes A, B, and F are associated with Zebra chip disease in America [3,34,35], whereas only haplotype A has been found in New Zealand [5,36]. Haplotype B has also been found in *Bactericera maculipennis* (Crawford) [37]. The remaining haplotypes are not vectored by *B. cockerelli* but by closely related species in the Triazidae family.

The impact of *B. cockerelli* and associated Lso transmission on agriculture is significant. Since its arrival in New Zealand circa 2005 via human-mediated dispersal it has caused millions of dollars of economic losses [6,21]. Similarly, management of *B. cockerelli* in the US is reported to have cost millions of dollars per year in major potato growing areas such as Texas [38] and the Pacific Northwest [39]. The introduction of *B. cockerelli* into potato growing regions in Europe or Asia would be devastating to the agricultural industry of those regions. If *B. cockerelli*, or a sufficient vector of Solanaceous Lso haplotypes, were to invade Europe it is estimated that the effects of Lso damage on potato and tomato would cost € 222 million per year and the negative impact of social welfare could cost an additional estimated € 114 million [40].

Currently, *B. cockerelli* is considered an A1 quarantine pest in the EPPO region [4]. Consignments of aubergine and *Capsicum* from Mexico infested with immature and adult stages of *B. cockerelli* were intercepted four times during UK border inspections between 2017–2018; indicating that there is a real threat of this pest making an incursion into the EPPO region if not properly monitored [41]. Monitoring and prevention of the spread of *B. cockerelli* is essential to prevent the risk of an outbreak of Lso on potato, tomato and other Solanaceous crops in areas where it is not currently found [42]. There is therefore an evident need for a rapid and accurate diagnostic test to identify *B. cockerelli* at all life stages not only as a tool to support import inspections, but also to assist monitoring, eradication and control strategies.

We designed a species-specific real-time PCR diagnostic assay to detect all life-stages of *B. cockerelli*, eggs, immatures and adults. The assay provides a rapid diagnostic test to quickly determine the presence of *B. cockerelli*, allowing for the early detection of invasions/introductions and aiding in the prevention of spread of this psyllid.
**2. Materials and methods**

**2.1. Specimen collection**

The assay was tested on 28 target adults *B. cockerelli* specimens and 73 non-target species cons-
sisting of 110 specimens see results section 3.1 for more info on samples. The classification fol-
lows Burckhardt & Ouvrard [43], and a complete taxonomic account of each species is given
in Ouvrard [20]. Psyllid identifications were confirmed against reference type specimens in
the NHM London collections. To account for intraspecific genetic variation, we obtained *B.
cockerelli* specimens from Mexico (Universidad Autónoma Agraria Antonio Narro) and USA
(USDA, Agricultural Research Services) from colony collections of each of the four recognised
biotypes of *B. cockerelli* in Central America, the Central, Western, Northwestern, and South-
western biotypes [19]. Specimens of *B. cockerelli* were also obtained from New Zealand lab-
reared colonies (Plant Research, New Zealand). Non-target specimens were mainly obtained
from 12.2 m suction-traps in the United Kingdom that form part of the Rothamsted Insect
Survey network described here [44]. Specimens were also obtained from suction-traps in Fin-
land, Germany, Spain and Sweden; as well as from field collections from Finland, Israel,
Mexico, Serbia, Spain, UK and USA. Non-target specimens from different regions of the USA
were used to test assay specificity on species that are commonly found in the same region as *B.
cockerelli*. As immatures and eggs are the most likely life stages that inspectors might find on
imported plant material, we also tested the assay on DNA extracted from immatures and eggs
from Mexico and the USA for validation.

**2.2. DNA extraction, PCR, and DNA sequencing for identification of
psyllids**

DNA for sequencing and assay validation was extracted from psyllids using a non-destructive
method first described in [45] and adapted from [46]. Psyllid specimens were preserved in
95% Ethanol: 5% Glycerol solution. Using a 15mm long, 0.15mm diameter stainless steel ento-
mological head-less pin (A3 size, Watkins and Doncaster) mounted in a holder, specimens
were initially pierced fully through the abdomen and half-way through the thorax from the
dorsal side while attempting to minimise damage to head, legs, wings, terminalia and other
body parts that are used for taxonomic identification. Pierced specimens were placed in a
microcentrifuge tube containing 180 µl of ATL buffer and 20 µl of proteinase-k as outlined in
the DNeasy Blood and Tissue Kit from Animal Tissues (Qiagen). Samples were placed in a
shaking incubator over-night (~8-10 hrs) at 56°C at 300 rpm. The subsequent steps from the
above mentioned protocol were followed and the psyllid integument voucher specimen was
stored in 95% Ethanol: 5% Glycerol for morphological identification. Psyllids were DNA bar-
coded using one or two gene regions. The internal transcribed spacer 2 (ITS2) and cytochrome
c oxidase subunit 1 (CO1) were amplified and sequenced for identification of different psyllid
species. For amplification of ITS2 primers CA55p8sFcm-F and CA28sB1d-R [47] were used;
and for amplification of CO1 gene regions arthropod barcoding Primers LCO1490 and
HCO2198 [48] were used. All reactions were performed in 20 µl consisting of: 10 µl 2x Type-It
Microsatellite PCR Kit Master Mix (Qiagen); 0.2 µM each forward and reverse primer; 7.2 µl
molecular grade water (Sigma-Aldrich) and 2 µl of psyllid template DNA. Reactions were run
on a Veriti 96-well thermal cycler (Applied Biosystems) using the following programs. ITS2:
95°C for 5 mins; 25 x cycles of (95°C for 30 s, 56°C for 90 s, 72°C for 30 s); and a final exten-
sion at 72°C for 10 mins. CO1: 94°C for 5 mins; 5 x cycles of (94°C for 30s, 45°C for 30s, 72°C
for 1 min); 25 x cycles of (94°C for 30s, 51°C for 1 min, 72°C for 1 min); and a final extension
of 72°C for 10 mins. PCR amplified gene regions were cleaned-up using EXO-SAP and
Ethanol precipitation, then sequenced using the BigDye Terminator Cycle Sequencing Kit (Applied Biosystems), forward and reverse complimentary DNA strands were sequenced separately for each sample and analysed using a 3500xL Genetic Analyser (Applied Biosystems).

**2.3. Bioinformatics and real-time PCR assay design**

Sequence editing, assembly and alignment were performed on “.AB1” trace files uploaded to Geneious R11 v 11.1.5 (Biomatters Ltd.). Contigs were assembled after trimming sections of low-quality sequence and aligning the complimentary strands using CLUSTAL-W multiple sequence alignment method [49]. Final contigs for each species and each gene region were aligned to identify variable areas suitable as targets for *B. cockerelli* specific primer and probe sets. Primers and probes were designed using manual selection of target-specific regions analysed using the “Basic Local Alignment Search Tool” (BLAST) [50] against the NCBI GenBank database [51] and processing of selected regions for suitability/ specificity in “Primer3” [52] and “Primer-BLAST” software [53]. Primer annealing temperature, hairpin formation, self-complementarity, GC content and were assessed using “Primer3” [52]. Potential amplification of non-specific insect species was checked using Primer BLAST which includes all psyllid species present in the GenBank database. Primer and probe sets were selected/rejected based on the following parameters: primer annealing temperature 59–62°C; primer annealing temperature + 8–10°C for probe annealing temperature; no more than 2°C difference in annealing temperature between primers, max probe length 30bp, no more than 3 Gs in a row in probe, amplicon length max 300bp and specificity to *B. cockerelli*.

**2.4. Real-time PCR set-up and standards**

To calculate standard curves DNA standards of *B. cockerelli* were prepared using dilution series of linearized cloned plasmid DNA. DNA was extracted as above using the non-destructive method, amplified and cloned into competent *Escherichia coli* cells using the TOPO TA cloning kit (Thermo-Fisher). DNA from successfully transformed colonies was extracted using “PureYield Plasmid Miniprep System” (Promega). For assay validation ITS2 DNA was cloned from other psyllid species (see results section 3.1). Stock DNA 10 ng/µl was linearised from cloned plasmid DNA using EcoRI restrictions enzyme (New England Biolabs), 0.5 µl of enzyme was added to 100 µl of stock DNA, this solution was incubated in a heat block (Thermomixer C, Eppendorf) at 37°C for 15 mins. The enzyme was then deactivated at 65°C for 20mins. Real-time PCRs were performed in 15 µl volumes including: 6.75 µl Jumpstart Taq Ready Mix (Sigma); 1.2 µl MgCl₂ (25mM); 0.45 µl of each primer; 0.15 µl probe; 4 µl of molecular grade water (Sigma); and 2 µl of template DNA. The standard real-time PCR cycle program was as follows. Hold stage: 50°C for 2 mins then; 95°C for 10 mins. PCR stage: 40 cycles of (95°C for 15 secs; X°C for 1 min), with primer annealing temperature X being 58, 60, 62, 64, or 68; depending on the experiment. Primer concentration, MgCl₂ concentration and temperature was adjusted for validation and optimization of the assay as described below. Reactions were performed on a “QuantStudio 6 Flex” (Applied Biosystems) real-time PCR machine and analysis was done on the “QuantStudio Real-Time PCR Software” (Applied Biosystems).

**2.5. Assay validation**

**2.5.1. Specificity.** The final primer and probe set was tested on genomic DNA from 47 *B. cockerelli* specimens from different life stages. These included the 4 US biotypes [17,54] and specimens from New Zealand to determine false negatives. The assay was tested for specificity against genomic DNA of 73 non-target psyllid species collected as mentioned above, to detect false positives. This included a total of 8 other closely related *Bactericera* spp. and the major
vectors of Lso on Apiaceous crops (*B. nigricornis*, *B. trigonica* and *Trioza apicalis*). Information regarding samples tested is in results section 3.1. The assay was also checked for cross-reaction against potato genomic DNA (*Solanum tuberosum*), 3 samples of *S. tuberosum* 'Maris Piper' were tested in replicates of 8. All reactions with non-target DNA were run in conjunction with a TaqMan Exogenous Internal Positive Control Reagent Kit (Applied Biosystems) to rule out the possibility that false positives were not obtained due to inhibition within the reaction. DNA from all non-target psyllids was sequenced in either ITS2, CO1 or both to ensure psyllid DNA was present in all reactions to rule out false negatives due to inefficient DNA extraction. Reactions were performed in duplicate at least, with a higher number of replicates for species closely related to *B. cockerelli*. False positives were defined as reactions with non-target DNA that showed fluorescence above the cycle threshold during 40 cycles; and false negatives were defined as reactions with *B. cockerelli* DNA that did not give a C~t~ after 40 cycles.

**2.5.2. Sensitivity.** Experiments were performed to determine the limit of detection of the assays. DNA standards were produced using *B. cockerelli* linearized cloned DNA from the ITS2 region. A nine point 10-fold dilution series starting with 10 ng/µl DNA up to 10^-8^ ng/µl of linearised plasmid DNA and genomic DNA was used to determine the limit of detection. 100ng/µl stock DNA concentration was initially checked using QuBit 4 Fluorometer (Invitrogen) and 5 µl was added to 45 µl of molecular grade water (Sigma-Aldrich) to dilute 1:10; eight subsequent dilutions were made. Linearised and non-linearised DNA was compared along with genomic DNA. The ability of the assay to detect immatures and eggs was also tested. DNA from various instars of immatures was extracted using the non-destructive protocol described above. Batches of 1 egg, 5 eggs and 10 eggs were extracted using the DNeasy Blood & Tissue kit (Qiagen) and initially broken with a pestle.

**2.5.3. Repeatability and reproducibility.** Variation in the performance of the assay between runs and within runs was assessed at a 0.2 µM primer concentration, with 1.5mM MgCl₂, and 60°C annealing temperature. Linearised plasmid DNA from *Escherichia coli* transformed with *B. cockerelli* ITS2 DNA was used. A six point 1:10 dilution series starting at 10ng/µl was used with each dilution being performed in triplicate. The same experiment was repeated 3x simultaneously. Runs and variations between the three experiments were recorded and analysed using QuantStudio 6 Real-Time PCR Software. An identical plate following the same plate set-up and reaction mix was run simultaneously on another QuantStudio 6 real-time PCR machine to compare inter-run variation.

**2.5.4. Robustness/optimization.** Amplification of target DNA, specificity and sensitivity at different MgCl₂ concentration, primer concentrations and annealing temperatures were performed to assess robustness. The assay was tested with 1.5, 3.5, 5.5, 7.5 and 9.5mM MgCl₂ concentration. For primers, 0.1, 0.2, 0.3, 0.5 and 1.0 µM concentrations were tested. The assay was also tested at different annealing temperatures 58, 60, 62, 64, 68°C across. For each tested parameter, optimization was performed across a nine point 1:10 dilution series starting at 10ng/µl DNA. All samples were tested in triplicates. Closely related *Bactericera* species were included in these assays to assess specificity under different assay conditions. After optimization of the assay a multifactorial robustness test was performed across two different real-time PCR machines to test the combined effects of small changes/errors in the PCR set-up. The assays were run on a "QuantStudio 6 Flex" (Applied Biosystems) and "CFX96 Real-Time System" (BioRad); results were analysed using "QuantStudio 6 Real-Time PCR Software" (Applied Biosystems) and "CFX Manager 3.1" (BioRad). The methodology used followed the European Network of GMO Laboratories (ENGL) recommendations [55].
## 3. Results

### 3.1. DNA extraction, PCR, and DNA sequencing for identification of psyllids

DNA from 110 psyllid specimens comprising 73 different species were extracted, amplified and sequenced successfully from either CO1 or ITS2 gene regions, or both (Table 1).

### 3.2. Bioinformatics and real-time PCR assay design

While differentiation within both the ITS2 and CO1 gene regions was sufficient to discriminate between psyllid species, the ITS2 gene region was more suitable for TaqMan assay design for *B. cockerelli*. Similarities between CO1 gene sequences between members of the *Bactericera* genus and *B. cockerelli* were higher than in the ITS2 region (average % similarity = 82.51 ± 0.68 for CO1 and 77.80 ± 4.79 for ITS2) (Table 2). The ITS2 region showed larger sections of variability along the gene on which to design primers and probes. Several primer and probe sets passed the selection criteria, but most were unsuitable due to high rate of false positives from closely related *Bactericera* species. The final primer and probe set Bcoc_JSK2 (Table 3) targets a 187bp region of the ITS2 gene (Fig 1).

### 3.3. Specificity and sensitivity

This assay did not amplify DNA from any of the 73 non-target psyllid species or *Solanum tuberosum* DNA when tested at 60°C with primer concentration 0.2 µM. Samples included nine closely related *Bactericera* species with similar ITS2 and CO1 sequences (Table 2). Under optimal conditions, false negatives = 0% for all non-target species tested with pure genomic DNA, giving a diagnostic specificity of 100%. Some suboptimal reaction conditions showed 33% false positives against high concentrations (10 ng / 1 ng) of *Bactericera albiventris* cloned DNA (see below). All *B. cockerelli* genomic DNA samples gave positive results (Table 4) giving 0% false negatives across 54 biological replicates and 147 technical replicates; resulting in a diagnostic sensitivity of 100%. These included *B. cockerelli* specimens from each of the four US biotypes as well as specimens from New Zealand. These specimens included adults, immature stages and eggs. The assay can amplify *B. cockerelli* DNA from both cloned and genomic samples. Under optimal conditions for PCR efficiency and specificity (60°C, 0.2 µM primer, 1.5 mM MgCl₂) the limit of detection was 0.000001 ng DNA across a range of different reaction parameters this equates to 200 copy numbers of ITS2 calculated using the following equation: Number of Copies = (ng DNA x 6.022x10²³) ÷ (length of plasmid (4656) + cloned fragment (700)bp) * 1x10⁹ * 660). The copy number calculator available at http://scienceprimer.com/copy-number-calculator-for-realtime-pcr was used. Diagnostic sensitivity was 100% on all DNA extracted from *B. cockerelli* immatures. False negatives from DNA from egg extractions were 0% for single eggs and 0% for batches of 3 and 10 eggs.

### 3.4. Repeatability and reproducibility

No significant differences were found between C~t~ means across the different replicates at different concentrations as tested by two-way ANOVA (F~5, 25~ = 0.54, *p* = 0.955). The assay also performed consistently across different machines and there was no significant difference between runs across the two machines as tested by two-way ANOVA (F~1, 5~ = 1.28, *p* = 0.279).
<table><caption>Table 1. Information on non-target psyllid species and plant specimens tested using the *B. cockerelli* real-time PCR assay Bcoc_JSK2 showing number of technical replicates and false positives.</caption><thead><tr><th>Family</th><th>Genus</th><th>Species</th><th>Voucher ID</th><th>Collection Location</th><th>Collection method</th><th>CO1 Ac#</th><th>ITS2 Ac#</th><th>Tech Reps</th><th>Voucher Location</th></tr></thead><tbody><tr><td rowspan="10">Aphalaridae</td><td rowspan="2">Aphalara</td><td>avicularis</td><td>160718.A. avi.23</td><td>Wellesbourne, UK</td><td>suction trap</td><td>MT021761</td><td>/</td><td>2</td><td>1</td></tr><tr><td>polygoni</td><td>160718.A. pol.22</td><td>Wellesbourne, UK</td><td>suction trap</td><td>/</td><td>MT038907</td><td>2</td><td>1</td></tr><tr><td rowspan="4">Blastopsylla</td><td>occidentalis</td><td>180312.Bl. occ.24</td><td>Salamanca, Spain</td><td>suction trap</td><td>MN272146</td><td>MN316692</td><td>2</td><td>3</td></tr><tr><td rowspan="3">Craspedolepta</td><td>gutierreziae</td><td>160825.5 US</td><td>Nevada, USA</td><td>field collection</td><td>MT021786</td><td>MT038962</td><td>2</td><td>1</td></tr><tr><td rowspan="2">minutissima</td><td>160825.1 US</td><td>Nevada, USA</td><td>field collection</td><td>MT021787</td><td>MT038963</td><td>2</td><td>1</td></tr><tr><td>160825.10 US</td><td>Oregon, USA</td><td>field collection</td><td>MT021789</td><td>/</td><td>2</td><td>1</td></tr><tr><td>nervosa</td><td>160728.Cra. ner.2</td><td>Gogarbank, UK</td><td>suction trap</td><td>MT021790</td><td>MT038964</td><td>2</td><td>1</td></tr><tr><td rowspan="3">Liviidae</td><td>pinicola</td><td>160825.2 US</td><td>Nevada, USA</td><td>field collection</td><td>/</td><td>MT038965</td><td>2</td><td>1</td></tr><tr><td>subpunctata</td><td>160421.C.sub.5</td><td>Gogarbank, UK</td><td>suction trap</td><td>MT021791</td><td>MT038966</td><td>2</td><td>1</td></tr><tr><td>Rhinocola</td><td>aceris</td><td>151014.R. ace.14</td><td>Wellesbourne, UK</td><td>suction trap</td><td>MT021810</td><td>MT038979</td><td>2</td><td>2</td></tr><tr><td rowspan="10">Psyllidae</td><td>Diaphorina</td><td>citri</td><td>160309.D.cit.6</td><td>Lab Colony, Vietnam</td><td>Lab Reared</td><td>MT021794</td><td>MT038969</td><td>2</td><td>1</td></tr><tr><td rowspan="3">Euphyllura</td><td>olivina</td><td>180125.Eup. oli.3</td><td>imports from Italy</td><td>imported Olea europeae</td><td>MT021797</td><td>MT038970</td><td>2</td><td>3</td></tr><tr><td rowspan="2">Livia</td><td>crefeldensis</td><td>180312.L.cre.5</td><td>Salamanca, Spain</td><td>suction trap</td><td>MN316678</td><td>MN272127</td><td>2</td><td>3</td></tr><tr><td>junci</td><td>160404.L.jun.1</td><td>Broom's Barn, UK</td><td>suction trap</td><td>MT021801</td><td>/</td><td>2</td><td>2</td></tr><tr><td rowspan="2">Arytaina</td><td>genistae</td><td>160825.6 US</td><td>Nevada, USA</td><td>field collection</td><td>MT021802</td><td>MT038973</td><td>2</td><td>1</td></tr><tr><td>gredi</td><td>151203.A. gen.2J</td><td>Ayr, UK</td><td>suction trap</td><td>/</td><td>MT038909</td><td>2</td><td>1</td></tr><tr><td rowspan="2">Arytainilla</td><td>gredi</td><td>180312.A.gre.1</td><td>Salamanca, Spain</td><td>suction trap</td><td>MN272123</td><td>MN316677</td><td>2</td><td>3</td></tr><tr><td>spartiophila</td><td>180716.A. spa.29</td><td>Edinburgh, UK</td><td>suction trap</td><td>MT021762</td><td>MT038908</td><td>2</td><td>3</td></tr><tr><td rowspan="2">Baeopelma</td><td>foersteri</td><td>151203.B. foe.1J</td><td>Ayr, UK</td><td>suction trap</td><td>/</td><td>MT038944</td><td>2</td><td>1</td></tr><tr><td>foersteri</td><td>160928.B.foe.2</td><td>SASA, UK</td><td>suction trap</td><td>MT021776</td><td>/</td><td>2</td><td>1</td></tr><tr><td rowspan="3">Cacopsylla</td><td>affinis</td><td>151203.C.aff.1</td><td>Wye, UK</td><td>suction trap</td><td>MT021777</td><td>MT038945</td><td>2</td><td>2</td></tr><tr><td rowspan="2">ambigua</td><td>160404.C. amb.4</td><td>Wye, UK</td><td>suction trap</td><td>/</td><td>MT038946</td><td>2</td><td>2</td></tr><tr><td>160104.C. amb.3</td><td>Preston, UK</td><td>suction trap</td><td>/</td><td>MT038947</td><td>2</td><td>1</td></tr><tr><td rowspan="4">americana</td><td rowspan="3">brunneipennis</td><td>160825.3 US</td><td>Nevada, USA</td><td>field collection</td><td>MT021778</td><td>MT038948</td><td>2</td><td>1</td></tr><tr><td>160309.C.bru.8</td><td>Wye, UK</td><td>suction trap</td><td>/</td><td>MT038949</td><td>2</td><td>2</td></tr><tr><td>160404.C.cra.3</td><td>Broom's Barn, UK</td><td>suction trap</td><td>MT021779</td><td>MT038950</td><td>2</td><td>2</td></tr><tr><td>mali</td><td>180910.C. mal.30</td><td>Elcho, UK</td><td>field collection</td><td>/</td><td>MT038951</td><td>2</td><td>3</td></tr><tr><td rowspan="2">melanoneura</td><td rowspan="2">moscovita</td><td>160718.C. mel.6</td><td>Kirton, UK</td><td>suction trap</td><td>/</td><td>MT038952</td><td>2</td><td>3</td></tr><tr><td>190109.C. mos.1</td><td>Germany</td><td>suction trap</td><td>/</td><td>/</td><td>2</td><td>3</td></tr><tr><td>peregrina</td><td>161024.C. per.11</td><td>Silwood Park, UK</td><td>suction trap</td><td>MT021780</td><td>MT038953</td><td>2</td><td>1</td></tr><tr><td>pruni</td><td>160203.C. pru.18</td><td>Wellesbourne, UK</td><td>suction trap</td><td>/</td><td>MT038954</td><td>2</td><td>2</td></tr></tbody></table>

(Continued)
<table><caption>Table 1. (Continued)</caption><thead><tr><th>Family</th><th>Genus</th><th>Species</th><th>Voucher ID</th><th>Collection Location</th><th>Collection method</th><th>CO1 Ac#</th><th>ITS2 Ac#</th><th>Tech Reps</th><th>Voucher Location</th></tr></thead><tbody><tr><td></td><td></td><td>pulchra</td><td>160718.C. pul.15</td><td>Elgin, UK</td><td>suction trap</td><td>/</td><td>MT038955</td><td>2</td><td>1</td></tr><tr><td></td><td></td><td>pyricola</td><td>160203.C. pco.2</td><td>Wye, UK</td><td>suction trap</td><td>MT021781</td><td>MT038956</td><td>2</td><td>2</td></tr><tr><td></td><td></td><td>saliceti</td><td>161024.C.sal.7</td><td>York, UK</td><td>suction trap</td><td>/</td><td>MT038958</td><td>2</td><td>1</td></tr><tr><td></td><td></td><td>sorbi</td><td>161024.C.sor.8</td><td>Preston, UK</td><td>suction trap</td><td>MT021782</td><td>MT038959</td><td>2</td><td>1</td></tr><tr><td></td><td></td><td>rhamnicola</td><td>151014.C.rha.8</td><td>Wellesbourne, UK</td><td>suction trap</td><td>/</td><td>MT038957</td><td>2</td><td>2</td></tr><tr><td></td><td></td><td>ulmi</td><td>171011.C. ulm.13</td><td>Germany</td><td>suction trap</td><td>MT021783</td><td>MT038960</td><td>2</td><td>3</td></tr><tr><td></td><td>Ceanothia</td><td>ceanothi</td><td>160825.9 US</td><td>Oregon, USA</td><td>field collection</td><td>MT021784</td><td>/</td><td>2</td><td>1</td></tr><tr><td></td><td>Chamaepsylla</td><td>hartigii</td><td>160728.Ch. har.1</td><td>Gogarbank, UK</td><td>suction trap</td><td>MT021785</td><td>MT038961</td><td>2</td><td>1</td></tr><tr><td></td><td>Euglyptoneura</td><td>fuscipennis</td><td>160825.7 US</td><td>Oregon, USA</td><td>field collection</td><td>MT021795</td><td>/</td><td>2</td><td>1</td></tr><tr><td></td><td></td><td>robusta</td><td>160825.8 US</td><td>Oregon, USA</td><td>field collection</td><td>MT021796</td><td>/</td><td>2</td><td>1</td></tr><tr><td></td><td>Heteropsylla</td><td>texana</td><td>160825.11 US</td><td>Texas, USA</td><td>field collection</td><td>MT021798</td><td>/</td><td>2</td><td>1</td></tr><tr><td></td><td>Psylla</td><td>alni</td><td>161019.P.aln.1</td><td>Sweden</td><td>suction trap</td><td>MT021804</td><td>/</td><td>2</td><td>1</td></tr><tr><td></td><td></td><td>buxi</td><td>180622.P. bux.22</td><td>Scotland, UK</td><td>suction trap</td><td>MT021806</td><td>MT038976</td><td>2</td><td>3</td></tr><tr><td></td><td></td><td>betulae</td><td>161123.P. bet.20</td><td>Jokioinen, Finland</td><td>suction trap</td><td>MT021805</td><td>MT038975</td><td>2</td><td>3</td></tr><tr><td></td><td>Psyllopsis</td><td>discrepans</td><td>151002.P.dis.8</td><td>Sweden</td><td>suction trap</td><td>MT021807</td><td>/</td><td>2</td><td>1</td></tr><tr><td></td><td></td><td>fraxini</td><td>180716.P.fri.33</td><td>Edinburgh, UK</td><td>suction trap</td><td>MT021808</td><td>MT038977</td><td>2</td><td>3</td></tr><tr><td></td><td></td><td>fraxinicola</td><td>160203.P.fra.6</td><td>Wellesbourne, UK</td><td>suction trap</td><td>MT021809</td><td>MT038978</td><td>2</td><td>2</td></tr><tr><td></td><td>Spanioneura</td><td>fonscolombii</td><td>180802.S. fon.29</td><td>Edinburgh, UK</td><td>field collection</td><td>/</td><td>MT038980</td><td>2</td><td>3</td></tr><tr><td rowspan="2">Spondyliaspidae</td><td rowspan="2">Ctenarytaina</td><td>spatulata</td><td>160404.Ct. spa.6</td><td>Wye, UK</td><td>suction trap</td><td>MT021792</td><td>MT038967</td><td>2</td><td>2</td></tr><tr><td>spatulata</td><td>161024.Ct. spa.5</td><td>Wye, UK</td><td>suction trap</td><td>MT021793</td><td>MT038968</td><td>2</td><td>1</td></tr><tr><td rowspan="8">Triozidae</td><td rowspan="8">Bactericera</td><td>albiventris</td><td>171214.B. alb.11</td><td>Jokioinen, Finland</td><td>suction trap</td><td>/</td><td>MT038910</td><td>5</td><td>3</td></tr><tr><td>curvatinervis</td><td>161123.B. cur.42</td><td>Jokioinen, Finland</td><td>suction trap</td><td>/</td><td>MT038911</td><td>5</td><td>3</td></tr><tr><td>dorsalis</td><td>160803.B.dor.2</td><td>Florida, USA</td><td>lab colony</td><td>MT021763</td><td>MT038912</td><td>5</td><td>3</td></tr><tr><td>maculipennis</td><td>190604.B. mac.1</td><td>Lab Colony, USA</td><td>Lab Reared</td><td>/</td><td>MT038913</td><td>2</td><td>3</td></tr><tr><td></td><td>190604.B. mac.2</td><td>Lab Colony, USA</td><td>Lab Reared</td><td>/</td><td>MT038914</td><td>2</td><td>3</td></tr><tr><td></td><td>190604.B. mac.3</td><td>Lab Colony, USA</td><td>Lab Reared</td><td>/</td><td>MT038915</td><td>2</td><td>3</td></tr><tr><td></td><td>190604.B. mac.4</td><td>Lab Colony, USA</td><td>Lab Reared</td><td>/</td><td>MT038916</td><td>2</td><td>3</td></tr><tr><td></td><td>190604.B. mac.5</td><td>Lab Colony, USA</td><td>Lab Reared</td><td>/</td><td>MT038917</td><td>2</td><td>3</td></tr><tr><td></td><td></td><td>190604.B. mac.6</td><td>Lab Colony, USA</td><td>Lab Reared</td><td>/</td><td>MT038918</td><td>2</td><td>3</td></tr><tr><td></td><td></td><td>190604.B. mac.7</td><td>Lab Colony, USA</td><td>Lab Reared</td><td>/</td><td>MT038919</td><td>2</td><td>3</td></tr><tr><td></td><td></td><td>nigricornis</td><td>170324.B. nig.18</td><td>Spain</td><td>field collection</td><td>MT021764</td><td>MT038920</td><td>5</td><td>3</td></tr></tbody></table>

(Continued)
<table><caption>Table 1. (Continued)</caption><thead><tr><th>Family</th><th>Genus</th><th>Species</th><th>Voucher ID</th><th>Collection Location</th><th>Collection method</th><th>CO1 Ac#</th><th>ITS2 Ac#</th><th>Tech Reps</th><th>Voucher Location</th></tr></thead><tbody><tr><td></td><td></td><td></td><td>170324.B. nig.22</td><td>Spain</td><td>field collection</td><td>MT021765</td><td>MT038921</td><td>5</td><td>3</td></tr><tr><td></td><td>salicivora</td><td></td><td>190116.B.sal.1</td><td>Elgin, UK</td><td>suction trap</td><td>/</td><td>/</td><td>6</td><td>3</td></tr><tr><td></td><td>striola</td><td></td><td>161123.B.str.9</td><td>Jokioinen, Finland</td><td>suction trap</td><td>/</td><td>MT038922</td><td></td><td></td></tr><tr><td></td><td>tremblayi</td><td></td><td>170731.B.tre.5</td><td>Belgrade, Serbia</td><td>field collection</td><td>/</td><td>MT038923</td><td>5</td><td>3</td></tr><tr><td></td><td></td><td></td><td>190604.B. tre.17</td><td>Spain</td><td>Lab Colony</td><td>/</td><td>MT038924</td><td>2</td><td>3</td></tr><tr><td></td><td></td><td></td><td>190604.B. tre.18</td><td>Spain</td><td>Lab Colony</td><td>/</td><td>MT038925</td><td>2</td><td>3</td></tr><tr><td></td><td></td><td></td><td>190604.B. tre.19</td><td>Spain</td><td>Lab Colony</td><td>/</td><td>MT038926</td><td>2</td><td>3</td></tr><tr><td></td><td></td><td></td><td>190604.B. tre.20</td><td>Spain</td><td>Lab Colony</td><td>/</td><td>MT038927</td><td>2</td><td>3</td></tr><tr><td></td><td></td><td></td><td>190604.B. tre.21</td><td>Spain</td><td>Lab Colony</td><td>/</td><td>MT038928</td><td>2</td><td>3</td></tr><tr><td></td><td>trigonica</td><td></td><td>170629.B.tri.16</td><td>Tunisia</td><td>field collection</td><td>MT021766</td><td>MT038929</td><td>3</td><td>3</td></tr><tr><td></td><td></td><td></td><td>170629.B.tri.17</td><td>Tunisia</td><td>field collection</td><td>/</td><td>MT038930</td><td>3</td><td>3</td></tr><tr><td></td><td></td><td></td><td>170629.B.tri.18</td><td>Tunisia</td><td>field collection</td><td>MT021767</td><td>MT038931</td><td>3</td><td>3</td></tr><tr><td></td><td></td><td></td><td>181010.B.tri.17</td><td>Spain</td><td>Lab Colony</td><td>MT021768</td><td>MT038932</td><td>2</td><td>3</td></tr><tr><td></td><td></td><td></td><td>181010.B.tri.18</td><td>Spain</td><td>Lab Colony</td><td>MT021769</td><td>MT038933</td><td>2</td><td>3</td></tr><tr><td></td><td></td><td></td><td>181010.B.tri.19</td><td>Spain</td><td>Lab Colony</td><td>/</td><td>MT038934</td><td>2</td><td>3</td></tr><tr><td></td><td></td><td></td><td>181010.B.tri.20</td><td>Spain</td><td>Lab Colony</td><td>MT021770</td><td>MT038935</td><td>2</td><td>3</td></tr><tr><td></td><td></td><td></td><td>181010.B.tri.21</td><td>Spain</td><td>Lab Colony</td><td>/</td><td>MT038936</td><td>2</td><td>3</td></tr><tr><td></td><td></td><td></td><td>190604.B.tri.23</td><td>Spain</td><td>Lab Colony</td><td>MT021771</td><td>MT038937</td><td>2</td><td>3</td></tr><tr><td></td><td></td><td></td><td>190604.B.tri.24</td><td>Spain</td><td>Lab Colony</td><td>/</td><td>MT038938</td><td>2</td><td>3</td></tr><tr><td></td><td></td><td></td><td>190604.B.tri.25</td><td>Spain</td><td>Lab Colony</td><td>MT021772</td><td>MT038939</td><td>2</td><td>3</td></tr><tr><td></td><td></td><td></td><td>190604.B.tri.26</td><td>Spain</td><td>Lab Colony</td><td>MT021773</td><td>MT038940</td><td>2</td><td>3</td></tr><tr><td></td><td></td><td></td><td>190604.B.tri.27</td><td>Spain</td><td>Lab Colony</td><td>MT021774</td><td>MT038941</td><td>2</td><td>3</td></tr><tr><td></td><td></td><td></td><td>190604.B.tri.28</td><td>Spain</td><td>Lab Colony</td><td>/</td><td>MT038942</td><td>2</td><td>3</td></tr><tr><td></td><td></td><td></td><td>190604.B.tri.29</td><td>Spain</td><td>Lab Colony</td><td>MT021775</td><td>MT038943</td><td>2</td><td>3</td></tr><tr><td>Heterotrioza</td><td>chenopodii</td><td></td><td>160203.H. che.11</td><td>Kirton, UK</td><td>suction trap</td><td>/</td><td>MT038971</td><td>2</td><td>2</td></tr><tr><td></td><td></td><td></td><td>160825.12 US</td><td>Washington, USA</td><td>field collection</td><td>MT021799</td><td>/</td><td>2</td><td>1</td></tr><tr><td>Lauritrioza</td><td>alacris</td><td></td><td>160816.L.ala.2</td><td>Spain</td><td>suction trap</td><td>MT021800</td><td>MT038972</td><td>2</td><td>1</td></tr><tr><td>Powellia</td><td>vitreoradiata</td><td></td><td>161024.P.vit.10</td><td>Kirton, UK</td><td>suction trap</td><td>MT021803</td><td>MT038974</td><td>2</td><td>1</td></tr><tr><td>Trioza</td><td>albifrons</td><td></td><td>160825.18.US</td><td>Nevada, USA</td><td>field collection</td><td>MT021811</td><td>MT038981</td><td>2</td><td>1</td></tr><tr><td></td><td>anthrisci</td><td></td><td>150708.T. ant.11</td><td>Jokioinen, Finland</td><td>field collection</td><td>MT021812</td><td>/</td><td>2</td><td>3</td></tr><tr><td></td><td>apicalis</td><td></td><td>161019.T{oapi.5</td><td>Sweden</td><td>field collection</td><td>MT021813</td><td>/</td><td>2</td><td>3</td></tr><tr><td></td><td>buxtoni</td><td></td><td>170324.T. bux.11</td><td>Israel</td><td>field collection</td><td>MT021814</td><td>MT038982</td><td>2</td><td>3</td></tr><tr><td></td><td>centranthi</td><td></td><td>161024.T.cen.9</td><td>Wye, UK</td><td>suction trap</td><td>MT021815</td><td>/</td><td>2</td><td>1</td></tr><tr><td></td><td>cerastii</td><td></td><td>171214.T. cer.32</td><td>Vikki, Finland</td><td>suction trap</td><td>MT021816</td><td>MT038983</td><td>2</td><td>3</td></tr><tr><td></td><td>dispar</td><td></td><td>160718.T. dis.26</td><td>Hellfreda, Sweden</td><td>suction trap</td><td>MT021817</td><td>/</td><td>2</td><td>1</td></tr><tr><td></td><td>erytreae</td><td></td><td>160808.ICA.19</td><td>Spain</td><td>Lab Colony</td><td>/</td><td>MT038984</td><td>2</td><td>1</td></tr><tr><td></td><td>flavipennis</td><td></td><td>160421.T{fla.3</td><td>Sweden</td><td>suction trap</td><td>MT021818</td><td>MT038985</td><td>2</td><td>1</td></tr></tbody></table>

(Continued)
<table><caption>**Table 1.** (Continued)</caption><thead><tr><th>Family</th><th>Genus</th><th>Species</th><th>Voucher ID</th><th>Collection Location</th><th>Collection method</th><th>CO1 Ac#</th><th>ITS2 Ac#</th><th>Tech Reps</th><th>Voucher Location</th></tr></thead><tbody><tr><td></td><td></td><td>galii</td><td>160203.T. gal.23</td><td>Wellesbourne, UK</td><td>suction trap</td><td>/</td><td>MT038986</td><td>2</td><td>2</td></tr><tr><td></td><td></td><td>remota</td><td>160718.T. rem.8</td><td>Sweden</td><td>suction trap</td><td>/</td><td>MT038987</td><td>2</td><td>1</td></tr><tr><td></td><td></td><td></td><td>180424.T. rem.1</td><td>Dundee, UK</td><td>Suction trap</td><td>MT021819</td><td>MT038988</td><td>3</td><td>3</td></tr><tr><td></td><td></td><td></td><td>180424.T. rem.6</td><td>Dundee, UK</td><td>Suction trap</td><td>MT021820</td><td>MT038989</td><td>3</td><td>3</td></tr><tr><td></td><td></td><td></td><td>180424.T. rem.16</td><td>Dundee, UK</td><td>Suction trap</td><td>MT021821</td><td>MT038990</td><td>3</td><td>3</td></tr><tr><td></td><td></td><td></td><td>180424.T. rem.18</td><td>Dundee, UK</td><td>Suction trap</td><td>MT021822</td><td>MT038991</td><td>3</td><td>3</td></tr><tr><td></td><td></td><td></td><td>180424.T. rem.19</td><td>Dundee, UK</td><td>Suction trap</td><td>/</td><td>MT038992</td><td>3</td><td>3</td></tr><tr><td></td><td></td><td></td><td>190116.T. rem.7</td><td>UK</td><td>Suction trap</td><td>MT021823</td><td>MT038993</td><td>3</td><td>3</td></tr><tr><td></td><td></td><td>rhamni</td><td>151002.T. rha.13</td><td>Sweden</td><td>suction trap</td><td>MT021824</td><td>MT038994</td><td>2</td><td>1</td></tr><tr><td></td><td></td><td>tatrensis</td><td>160718.T. tat.27</td><td>Sweden</td><td>suction trap</td><td>/</td><td>MT038995</td><td>2</td><td>1</td></tr><tr><td></td><td></td><td>urticae</td><td>160816.T. urt.17</td><td>Spain</td><td>field collection</td><td>/</td><td>MT038996</td><td>2</td><td>1</td></tr></tbody></table>

All non-target species gave 0% false positives. GenBank Accession numbers are included for ITS2 and CO1 regions if sequencing was successful. Voucher Location: 1 = 1; 2 = 2 Research Insect Survey; 3 = SASA Hemipteran DNA Database. All DNA samples are stored in the SASA Hemipteran DNA database. "/" = no sequence obtained.

https://doi.org/10.1371/journal.pone.0230741.t001

<table><caption>**Table 2.** Closely related *Bactericera* species tested with Bcoc_JSK2 assay.</caption><thead><tr><th rowspan="2">Species</th><th colspan="3">ITS2</th><th colspan="3">CO1</th></tr><tr><th>% Similarity</th><th>bp</th><th>GC content %</th><th>% Similarity</th><th>bp</th><th>GC content %</th></tr></thead><tbody><tr><td>B. trigonica</td><td>78.96</td><td>662</td><td>59.3</td><td>82.88</td><td>509</td><td>35.4</td></tr><tr><td>B. tremblayi</td><td>79.16</td><td>665</td><td>59.1</td><td>82.97</td><td>682</td><td>33</td></tr><tr><td>B. curvatinervis</td><td>80.30</td><td>655</td><td>58</td><td>82.23</td><td>678</td><td>34.7</td></tr><tr><td>B. nigricornis</td><td>81.16</td><td>668</td><td>59.3</td><td>81.28</td><td>521</td><td>36.7</td></tr><tr><td>B. albiventris</td><td>76.67</td><td>667</td><td>59.2</td><td>83.41</td><td>663</td><td>32.9</td></tr><tr><td>B. dorsalis</td><td>65.59</td><td>560</td><td>61.3</td><td>82.31</td><td>685</td><td>32.6</td></tr><tr><td>B. maculipennis</td><td>80.67</td><td>674</td><td>61.6</td><td>nd</td><td>nd</td><td>nd</td></tr><tr><td>B. salicivora</td><td>nd</td><td>nd</td><td>nd</td><td>nd</td><td>nd</td><td>nd</td></tr><tr><td>B. striola</td><td>79.91</td><td>663</td><td>59.1</td><td>nd</td><td>nd</td><td>nd</td></tr><tr><td>B.cockerelli</td><td>N/A</td><td>569</td><td>61.0</td><td>N/A</td><td>595</td><td>32.6</td></tr></tbody></table>

ITS similarity = % identity to DNA sample 150727.B.coc.02. CO1 similarity = % identity to a consensus sequences of all *B. cockerelli* sequences obtained during this study. CO1 genes showed higher similarity and fewer variable regions compared to ITS2 regions. Highest % similarity to *B. cockerelli* in the ITS2 region was found in *B. nigricornis* (81.16) and to *B. albiventris* in the CO1 region (83.41). The Bcoc_JSK2 assay does not give false positives with any of the species listed here. (nd = not determined due to sequencing failing).

https://doi.org/10.1371/journal.pone.0230741.t002
<table><caption>Table 3. Final oligonucleotide sequences for the Bcoc_JSK2 TaqMan real-time PCR assay to identify *B. cockerelli*. The assay targets a 187 bp region of the ITS2 gene region.</caption><thead><tr><th>Oligo Name</th><th>Function</th><th>Sequence 5'-3'</th><th>Tm</th><th>Length (bp)</th></tr></thead><tbody><tr><td>Bcoc_JSK2-f</td><td>forward primer</td><td>GAGGTCTCCTCATCGTGCGT</td><td>61</td><td>25</td></tr><tr><td>Bcoc_JSK2-r</td><td>reverse primer</td><td>GGACGAGCATTGCTGCTGC</td><td>62.2</td><td>23</td></tr><tr><td>Bcoc_JSK2-p</td><td>probe (FAM-BHQ)</td><td>GCAAACGCGGCACAAAGTACCGCGC</td><td>70.9</td><td>25</td></tr></tbody></table>

https://doi.org/10.1371/journal.pone.0230741.t003

## 3.5. Robustness/optimization

The assays amplified *B. cockerelli* DNA at all primer concentrations, MgCl₂ concentrations and annealing temperatures with varying levels of efficiency, precision, and sensitivity (S1–S3 Tables). At primer concentration 0.5 µM, the assay was less sensitive only amplifying down to 0.0001 ng DNA. At higher primer concentrations (1.0 µM), the assay showed higher sensitivity, but efficiency was outside the range for acceptable use. The assay performed optimally at 0.2 µM primer concentration showing good efficiency and high sensitivity (0.000001 ng DNA) (S1 Table). Generally, standard deviation of the C~t~ was lower at higher DNA concentrations and some of the primer concentrations showed SD slightly above the accepted level for

<figure><img src="image_3.png" /><figcaption>Fig 1. CLUSTAL-W alignment of ITS2 regions from closely related *Bactericera* species showing variable regions and the gene target for the Bcoc_JSK2 primer and probe set. Bases shades with black show differences to *B. cockerelli* sequence. Colour highlights locations of forward primer (blue highlight); reverse primer (green highlight) and probe (yellow highlight). The probe and reverse primer are reverse compliments of the highlighted regions here.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0230741.g001
<table><caption>**Table 4.** Information on *Bactericera cockerelli* samples tested with Bcoc_JSK2 assay including genomic DNA from adults, immatures, single eggs and egg batches.
Location of samples collection is also included. All samples gave 100% positives. Accession numbers for CO1 and ITS2 (MT027551-MT027599) regions are included. "/" = no sequence obtained.</caption><thead><tr><th>Sample name</th><th>Life Stage</th><th>Origin</th><th>C<sub>t</sub> ave</th><th>Tech reps</th><th>CO1 Ac#</th><th>ITS2 Ac#</th><th>DNA Source</th></tr></thead><tbody><tr><td>181119.B.coc.06</td><td>1 egg</td><td>Mexico</td><td>29.80</td><td>2</td><td>/</td><td>MT027568</td><td>Genomic</td></tr><tr><td>191003.B.coc.01</td><td>1 egg</td><td>Mexico</td><td>33.41</td><td>3</td><td>/</td><td>MT027592</td><td>Genomic</td></tr><tr><td>191003.B.coc.02</td><td>1 egg</td><td>Mexico</td><td>24.95</td><td>3</td><td>/</td><td>MT027593</td><td>Genomic</td></tr><tr><td>191003.B.coc.03</td><td>1 egg</td><td>Mexico</td><td>33.79</td><td>3</td><td>/</td><td>MT027594</td><td>Genomic</td></tr><tr><td>191003.B.coc.04</td><td>1 egg</td><td>Mexico</td><td>22.43</td><td>6</td><td>/</td><td>MT027595</td><td>Genomic</td></tr><tr><td>181119.B.coc.07</td><td>5 eggs</td><td>Mexico</td><td>24.42</td><td>2</td><td>/</td><td>MT027569</td><td>Genomic</td></tr><tr><td>181119.B.coc.21</td><td>5 eggs</td><td>Mexico</td><td>28.32</td><td>2</td><td>/</td><td>MT027582</td><td>Genomic</td></tr><tr><td>181119.B.coc.08</td><td>10 eggs</td><td>Mexico</td><td>29.61</td><td>2</td><td>/</td><td>MT027570</td><td>Genomic</td></tr><tr><td>181119.B.coc.22</td><td>10 eggs</td><td>Mexico</td><td>26.43</td><td>2</td><td>/</td><td>MT027583</td><td>Genomic</td></tr><tr><td>181119.B.coc.03</td><td>immature</td><td>Mexico</td><td>22.56</td><td>2</td><td>/</td><td>MT027565</td><td>Genomic</td></tr><tr><td>181119.B.coc.04</td><td>immature</td><td>Mexico</td><td>22.33</td><td>2</td><td>/</td><td>MT027566</td><td>Genomic</td></tr><tr><td>181119.B.coc.05</td><td>immature</td><td>Mexico</td><td>21.46</td><td>2</td><td>/</td><td>MT027567</td><td>Genomic</td></tr><tr><td>181119.B.coc.11</td><td>immature</td><td>Mexico</td><td>23.16</td><td>2</td><td>/</td><td>MT027573</td><td>Genomic</td></tr><tr><td>181119.B.coc.12</td><td>immature</td><td>Mexico</td><td>24.15</td><td>2</td><td>/</td><td>MT027574</td><td>Genomic</td></tr><tr><td>181119.B.coc.13</td><td>immature</td><td>Mexico</td><td>23.94</td><td>2</td><td>/</td><td>MT027575</td><td>Genomic</td></tr><tr><td>181119.B.coc.14</td><td>immature</td><td>Mexico</td><td>25.75</td><td>2</td><td>/</td><td>MT027576</td><td>Genomic</td></tr><tr><td>181119.B.coc.16</td><td>immature</td><td>Mexico</td><td>23.49</td><td>2</td><td>/</td><td>MT027578</td><td>Genomic</td></tr><tr><td>181119.B.coc.18</td><td>immature</td><td>Mexico</td><td>22.45</td><td>2</td><td>/</td><td>MT027580</td><td>Genomic</td></tr><tr><td>181119.B.coc.19</td><td>immature</td><td>Mexico</td><td>23.50</td><td>2</td><td>/</td><td>MT027581</td><td>Genomic</td></tr><tr><td>190604.B.coc.13</td><td>immature</td><td>Mexico</td><td>24.96</td><td>2</td><td>/</td><td>MT027588</td><td>Genomic</td></tr><tr><td>190604.B.coc.14</td><td>immature</td><td>Mexico</td><td>25.09</td><td>2</td><td>/</td><td>MT027589</td><td>Genomic</td></tr><tr><td>190604.B.coc.15</td><td>immature</td><td>Mexico</td><td>28.37</td><td>2</td><td>/</td><td>MT027590</td><td>Genomic</td></tr><tr><td>150727.B.coc.02</td><td>Adult</td><td>South Western, USA</td><td>22.18</td><td>2</td><td>MT040955</td><td>MG719775</td><td>Genomic</td></tr><tr><td>150827.B.coc.02</td><td>Adult</td><td>South Western, USA</td><td>22.18</td><td>2</td><td>MT040956</td><td>MT027597</td><td>Genomic</td></tr><tr><td>150827.B.coc.03</td><td>Adult</td><td>Central USA</td><td>24.49</td><td>6</td><td>MT040957</td><td>MT027598</td><td>Genomic</td></tr><tr><td>150827.B.coc.04</td><td>Adult</td><td>North Western, USA</td><td>24.77</td><td>2</td><td>MT040958</td><td>MT027599</td><td>Genomic</td></tr><tr><td>150827.B.coc.06</td><td>Adult</td><td>North Western, USA</td><td>23.68</td><td>2</td><td>MT040960</td><td>MT027552</td><td>Genomic</td></tr><tr><td>150827.B.coc.12</td><td>Adult</td><td>Western, USA</td><td>20.39</td><td>2</td><td>MT040961</td><td>MT027596</td><td>Genomic</td></tr><tr><td>150827.B.coc.17</td><td>Adult</td><td>South Western, USA</td><td>19.65</td><td>2</td><td>MT040962</td><td>MT027553</td><td>Genomic</td></tr><tr><td>160725.B.coc.05</td><td>Adult</td><td>Central, USA</td><td>21.45</td><td>2</td><td>MT040963</td><td>/</td><td>Genomic</td></tr><tr><td>160726.B.coc.01</td><td>Adult</td><td>New Zealand</td><td>21.56</td><td>2</td><td>/</td><td>MT027557</td><td>Genomic</td></tr><tr><td>160726.B.coc.02</td><td>Adult</td><td>New Zealand</td><td>21.02</td><td>2</td><td>/</td><td>MT027558</td><td>Genomic</td></tr><tr><td>160726.B.coc.03</td><td>Adult</td><td>New Zealand</td><td>20.48</td><td>2</td><td>/</td><td>MT027559</td><td>Genomic</td></tr><tr><td>160726.B.coc.04</td><td>Adult</td><td>New Zealand</td><td>21.98</td><td>2</td><td>/</td><td>MT027560</td><td>Genomic</td></tr><tr><td>160726.B.coc.05</td><td>Adult</td><td>New Zealand</td><td>19.43</td><td>2</td><td>/</td><td>MT027561</td><td>Genomic</td></tr><tr><td>160726.B.coc.06</td><td>Adult</td><td>New Zealand</td><td>20.96</td><td>2</td><td>/</td><td>MT027562</td><td>Genomic</td></tr><tr><td>180731.B.coc.04</td><td>Adult</td><td>North Western, USA</td><td>24.42</td><td>6</td><td>MT040964</td><td>/</td><td>Genomic</td></tr><tr><td>180731.B.coc.05</td><td>Adult</td><td>Western, USA</td><td>22.91</td><td>6</td><td>MT040965</td><td>/</td><td>Genomic</td></tr><tr><td>180731.B.coc.06</td><td>Adult</td><td>Western, USA</td><td>27.14</td><td>6</td><td>MT040966</td><td>/</td><td>Genomic</td></tr><tr><td>181119.B.coc.01</td><td>Adult</td><td>Mexico</td><td>21.47</td><td>2</td><td>/</td><td>MT027563</td><td>Genomic</td></tr><tr><td>181119.B.coc.02</td><td>Adult</td><td>Mexico</td><td>19.98</td><td>2</td><td>/</td><td>MT027564</td><td>Genomic</td></tr><tr><td>181119.B.coc.09</td><td>Adult</td><td>Mexico</td><td>21.83</td><td>2</td><td>/</td><td>MT027571</td><td>Genomic</td></tr><tr><td>181119.B.coc.10</td><td>Adult</td><td>Mexico</td><td>19.48</td><td>2</td><td>/</td><td>MT027572</td><td>Genomic</td></tr><tr><td>181119.B.coc.15</td><td>Adult</td><td>Mexico</td><td>21.27</td><td>2</td><td>/</td><td>MT027577</td><td>Genomic</td></tr><tr><td>181119.B.coc.17</td><td>Adult</td><td>Mexico</td><td>23.74</td><td>2</td><td>/</td><td>MT027579</td><td>Genomic</td></tr><tr><td>190604.B.coc.09</td><td>Adult</td><td>USDA, Lab Colony</td><td>21.51</td><td>2</td><td>/</td><td>MT027584</td><td>Genomic</td></tr></tbody></table>

(Continued)
<table><caption>**Table 4.** (Continued)</caption><thead><tr><th>Sample name</th><th>Life Stage</th><th>Origin</th><th>C<sub>t</sub> ave</th><th>Tech reps</th><th>CO1 Ac#</th><th>ITS2 Ac#</th><th>DNA Source</th></tr></thead><tbody><tr><td>190604.B.coc.10</td><td>Adult</td><td>Mexico</td><td>20.33</td><td>2</td><td>/</td><td>MT027585</td><td>Genomic</td></tr><tr><td>190604.B.coc.11</td><td>Adult</td><td>Mexico</td><td>22.67</td><td>2</td><td>/</td><td>MT027586</td><td>Genomic</td></tr><tr><td>190604.B.coc.12</td><td>Adult</td><td>Mexico</td><td>24.37</td><td>2</td><td>/</td><td>MT027587</td><td>Genomic</td></tr><tr><td>190604.B.coc.16</td><td>Adult</td><td>Mexico</td><td>27.15</td><td>2</td><td>/</td><td>MT027591</td><td>Genomic</td></tr><tr><td>150827.B.coc.05.col.04</td><td>transformed <i>E. coli</i></td><td>Lab</td><td>11.23</td><td>6</td><td>MT040959</td><td>MT027551</td><td>Cloned, 10ng</td></tr><tr><td>160725.B.coc.01.col.06</td><td>transformed <i>E. coli</i></td><td>Lab</td><td>11.55</td><td>6</td><td>/</td><td>MT027554</td><td>Cloned, 10ng</td></tr><tr><td>160725.B.coc.06.col.04</td><td>transformed <i>E. coli</i></td><td>Lab</td><td>11.78</td><td>6</td><td>/</td><td>MT027555</td><td>Cloned, 10ng</td></tr><tr><td>160725.B.coc.07.col.08</td><td>transformed <i>E. coli</i></td><td>Lab</td><td>11.67</td><td>6</td><td>/</td><td>MT027556</td><td>Cloned, 10ng</td></tr></tbody></table>

https://doi.org/10.1371/journal.pone.0230741.t004

quantitative real-time PCR, however this module is intended for qualitative use. At high DNA concentrations all primer concentrations are suitable for use with Bcoc_JSK2 primer and probe set to detect *B. cockerelli* but 0.2 µM is recommended for best results. The assay did not amplify non-target DNA from the 8 other *Bactericera* species tested at the different primer concentrations (0.1, 0.2, 0.3, 0.5 and 1.0 µM).

The MgCl2 concentration of the assay made only small differences to the overall performance of the assay (S2 Table) and the assay was able to amplify *B. cockerelli* DNA at low concentrations (0.000001 ng) at each MgCl2 concentration. The precision of the assay was lower at higher MgCl2 concentrations 7.5mM and 9.5mM (S2 Table).

Sensitivity was slightly higher at 64°C giving 33.33% (n = 3) positives for only 20 copies of *B. cockerelli* DNA (0.0000001 ng), however at 64°C and 66°C 33.33% (n = 3) false positives were found with 10ng and 1 ng of *B. albiventris* cloned DNA (S3 Table). Reactions at 58°C were 10 to 100-fold less sensitive than reactions at 64°C. For best sensitivity and specificity, it is suggested that assays using the Bcoc_JSK2 primer and probe set should be performed at 60°C or 62°C. While higher temperatures appear to be more sensitive, they are not recommended on unknown samples due to the small likelihood of returning false positives with *B. albiventris* and possibly other un-tested *Bactericera* spp.

It is recommended that this assay be performed at 60°C– 62°C, with a MgCl2 concentration of 1.5mM and a primer concentration of 0.2 µM. To test the robustness of these conditions a multifactorial approach was taken [55]. The assay performed satisfactorily across the different treatments and was shown to be robust and unaffected by small changes in assay set-up (S4 Table). Each treatment gave 100% positives for amplification of *B. cockerelli* genomic DNA.

## 4. Discussion

The Tomato-Potato psyllid is an economically damaging pest of solanaceous plants that has spread by human mediated dispersal. It causes feeding damage to plants but also is the major vector of ‘*Candidatus* Liberibacter solanacearum’ (Lso), a phloem limited bacterium that is associated with disease in solanaceous and apiaceous plants. Management of this insect pest requires accurate identification of *B. cockerelli*, this is often difficult if eggs or immature life stages only are available for identification. Hitherto, identification of *B. cockerelli* required either considerable expertise in psyllid taxonomy or the lengthy process of DNA barcoding [54].

We have designed and validated the first species-specific, quantitative real-time PCR TaqMan assay for *B. cockerelli* by using the comparison of 73 non-target species to identify unique gene regions that were suitable for primer/probe design and species differentiation. The genus *Bactericera* currently contains 160 species [20] and <1% of these have been tested in the current study due to the difficulty in obtaining other specimens from the field or lab colonies.
However Europe is home to 26 different species of *Bactericera* [20], 30% of which have been tested for false positives using this assay. Psyllid species that were tested are most commonly found in potato and carrot fields in Europe and the wider EPPO region which should minimize the potential for false positives and ensure the assay is efficient at detecting outbreaks in European fields. The assay was also tested on nine closely related *Bactericera* species. The number of species used in our study is relatively high compared to other reported TaqMan assays for plant pests that report lower numbers of non-target species [56,57].

The assay is based on a 187 bp region of the ITS2 gene which was suitable as it contained high interspecific variation consisting of stretches of insertions and deletions (INDELs). The ITS2 region has been used to distinguish species phylogenetically and to identify cryptic species in the *Cacopsylla pruni* complex [47]. DNA sequences obtained from this study will improve psyllid representation on online DNA databases, reducing the chance of Type II errors (i.e. misidentification due to lack of conspecific references) [58]. The *B. cockerelli* sequences on which we tested this assay (and many of the non-target psyllid species) were from different geographic locations to account for intraspecific variation. *Bactericera cockerelli* specimens from the four USA biotypes and specimens from New Zealand all gave 100% true positives.

The success rates of eradications are dependent on the length of time between introduction, detection, and implementation of eradication measures as Lso displays a short transmission time from *B. cockerelli* to potatoes [4,25]. Feasibly, methodology described in this study could be used to extract DNA from a specimen and test for *B. cockerelli* positives within 6–12 hrs or quicker. This is faster than identification by DNA barcoding and could aid in eradications/ prevention of incursions. This time could be reduced further if the real-time assay is used in conjunction with faster DNA extraction protocols.

There are currently no methods described within the EPPO “agreed diagnostic protocol for identification of *B. cockerelli*” [4]. In addition, the current EPPO control system for *B. cockerelli* and Lso [4] highlights the importance of identifying psyllid eggs and immatures on various plant materials during inspections and monitoring but gives minimal guidelines for achieving this. Validation of this assay demonstrates that it would be a reliable and accurate tool for use in this area and it will therefore be prepared for consideration by the EPPO diagnostic panel. This assay is also useful for monitoring *B. cockerelli* occurrence at several spatial scales, from local border checks to regional surveys which use different trapping methods (water, sticky, suction, aerial balloon traps) where no host plant data is available. Given the sensitivity of this assay it should be possible to detect *B. cockerelli* DNA from insect fragments (e.g. legs, heads) if DNA extraction is adequate. However, further validation should be performed to ensure the assay performs adequately on samples obtained from different traps. This assay should be tested on additional congeneric species and other closely related Triozidae psyllids. Another limitation of this assay is that it cannot yet be taken out into the field, making it less portable than LAMP assays or other NGS sequencing techniques such as Nanopore technology.

In conclusion a rapid, specific, robust, repeatable and reliable real-time PCR assay has now been validated and can be used to detect the important pest *B. cockerelli*. This will be an important tool for providing much-needed support to prevent new outbreaks. The assay can be implemented by practitioners with molecular biology experience and does not require personnel to have classical taxonomic knowledge of insects or psyllids; making this tool more accessible than traditional methods. The assay can be used to complement field surveillance and may facilitate further ecological studies of *B. cockerelli* requiring the identification of immatures and eggs. The strength of this assay lies in the collaboration of molecular biologists and classical taxonomists working together to build a reliable database for DNA barcoding of psyllids.
**Supporting information**

**S1 Table.** Assay performance across a range of primer concentrations at 60°C and 1.5mM MgCl₂. Optimum primer concentration was 0.2 µM showing the best combination of r², slope, efficiency, and sensitivity.
(DOCX)

**S2 Table.** Performance of *B. cockerelli* real-time PCR assay at different magnesium chloride (MgCl₂) concentrations.
(DOCX)

**S3 Table.** Summary of standard curves from optimisation of temperature on Bcoc_JSK2 real-time PCR assay for identification of B. cockerelli. All DNA concentrations tested above the limit of detection (10ng, 1 ng, 0.1ng, 0.01ng, 0.001ng, 0.0001ng, 0.00001ng) gave 100% positives across 3 x replicates. LOD is given for each temperature. All non-target Bactericera species tested at different DNA concentration gave 0% false positives except for B. albiventris cloned DNA which cross reacted at 64 and 66°C. (*reactions at 64°C gave 33.33% positives at 20 copy numbers).
(DOCX)

**S4 Table.** Set-up and results of multifactorial robustness experiment testing the Bcoc_JSK2 assay on *B. cockerelli* genomic DNA. All treatments showed 100% positives despite small changes to the overall set-up.
(DOCX)

**Acknowledgments**

We thank A. Fereres & C. A. Antolínez Delgado (Institute of Agricultural Sciences, CSIC, Spain), A. Nissinen (Natural Resources Institute Finland), J. Munyaneza, R. Cooper, M. Heidt, K. Swisher Grimm (USDA Agricultural Research Services), S. Bulman (Plant and Food Research, New Zealand), A. Jensen, S. Halbert (Florida Department of Agriculture & Consumer Services, Dept. of Plant Industry) and Alberto Flores (Universidad Autónoma Agraria Antonio Narro) for specimens; and thank C. Jeffries, L. Webster, V. Mulholland, and A. Reid (SASA) for providing advice. We also thank SASA Potato Genotyping team for providing Potato DNA.

**Author Contributions**

**Conceptualization:** J. C. Sumner-Kalkun, M. J. Sjölund, F. Highet, J. R. Bell, D. M. Kenyon.

**Data curation:** J. C. Sumner-Kalkun, M. J. Sjölund, Y. M. Arnsdorf, M. Carnegie, F. Highet, D. Ouvrard, A. F. C. Greenslade, J. R. Bell, R. Sigvald.

**Formal analysis:** J. C. Sumner-Kalkun, M. J. Sjölund, Y. M. Arnsdorf, M. Carnegie, D. Ouvrard, A. F. C. Greenslade.

**Funding acquisition:** J. R. Bell, D. M. Kenyon.

**Investigation:** J. C. Sumner-Kalkun, M. J. Sjölund, Y. M. Arnsdorf, D. Ouvrard, A. F. C. Greenslade.

**Methodology:** J. C. Sumner-Kalkun, M. J. Sjölund, F. Highet, J. R. Bell, D. M. Kenyon.

**Project administration:** J. C. Sumner-Kalkun, M. J. Sjölund, F. Highet, J. R. Bell, D. M. Kenyon.
**Resources:** M. J. Sjölund, M. Carnegie, D. Ouvrard, A. F. C. Greenslade, J. R. Bell, R. Sigvald, D. M. Kenyon.

**Software:** J. C. Sumner-Kalkun.

**Supervision:** J. C. Sumner-Kalkun, M. J. Sjölund, F. Highet, J. R. Bell, D. M. Kenyon.

**Validation:** J. C. Sumner-Kalkun, M. J. Sjölund, Y. M. Arnsdorf.

**Visualization:** J. C. Sumner-Kalkun, M. J. Sjölund.

**Writing – original draft:** J. C. Sumner-Kalkun, M. J. Sjölund, F. Highet.

**Writing – review & editing:** J. C. Sumner-Kalkun, M. J. Sjölund, F. Highet, D. Ouvrard, A. F. C. Greenslade, J. R. Bell.

## References

1. Munyaneza JE, Sengoda VG, Crosslin JM, Garzon-Tiznado JA, Cardenas-Valenzuela OG. First report of "Candidatus Liberibacter solanacearum" in tomato plants in Mexico. Plant Dis. 2009; 93: 1076. https://doi.org/10.1094/PDIS-93-10-1076A PMID: 30754366

2. Wallis RL. Ecological studies on the Potato Psyllid as a pest of potatoes. USDA Tech Bull. 1955; 1107: 25.

3. Hansen AK, Trumble JT, Stouthamer R, Paine TD. A new huanglongbing species, "Candidatus Liberibacter" bacter psyllaurous," found to infect tomato and potato, is vectored by the psyllid *Bactericera cockerelli* (Sulc). Appl Environ Microbiol. 2008; 74: 5862–5865. https://doi.org/10.1128/AEM.01268-08 PMID: 18676707

4. EPPO. PM 9/25 (1) Bactericera cockerelli and 'Candidatus Liberibacter solanacearum.' EPPO Bull. 2017; 47: 513–523. https://doi.org/10.1111/epp.12442

5. Liefting LW, Perez-Egusquiza ZC, Clover GRG, Anderson JAD. A new 'Candidatus Liberibacter' species in Solanum tuberosum in New Zealand. Plant Dis. 2008; 92: 1474–1474. https://doi.org/10.1094/pdis-92-10-1474a PMID: 30769561

6. Liefting LW, Sutherland PW, Ward LI, Paice KL, Weir BS, Clover GRG. A new 'Candidatus Liberibacter' species associated with diseases of solanaceous crops. Plant Dis. 2009; 93: 208–214. https://doi.org/10.1094/PDIS-93-3-0208 PMID: 30764179

7. Munyaneza JE. Zebra Chip disease of potato: Biology, epidemiology, and management. Am J Potato Res. 2012; 89: 329–350. https://doi.org/10.1007/s12230-012-9262-3

8. Crosslin JM, Hamm PB, Eggers JE, Rondon SI, Sengoda VG, Munyaneza JE. First report of Zebra Chip disease and "Candidatus Liberibacter solanacearum" on potatoes in Oregon and Washington State. Plant Dis. 2012; 96: 452. https://doi.org/10.1094/PDIS-10-11-0894 PMID: 30727109

9. Pletsch D. The Potato Psyllid *Paratrioza cockerelli* (Sulc), its biology and control. Bull Mont Agric Exp Stn. 1947; 446: 1–95.

10. Ouvrard D. 2018 Psyl'list—The World Psylloidea Database. 2019. Available: http://www.hemiptera-databases.com/psyllist

11. Teulon DA, Workman PJ, Thomas KL, Nielsen MC. *Bactericera cockerelli*: Incursion, dispersal and current distribution on vegetable crops in New Zealand. New Zeal Plant Prot. 2009; 62: 136–144.

12. Australian Government Department of Agriculture and Water Resources. Detection of Bactericera cockerelli (Tomato-potato psyllid) in Western Australia. Canberra; 2017. Available: https://www.ippc.int/en/countries/australia/pestreports/2017/02/detection-of-bactericera-cockerelli-tomato-potato-psyllid-in-western-australia/

13. Kaur N, Cooper WR, Duringer JM, Badillo-Vargas IE, Esparza-Díaz G, Rashed A, et al. Survival and development of Potato Psyllid (Hemiptera: Triozidae) on Convolvulaceae: Effects of a plant-fungus symbiosis (Periglandula). PLoS One. 2018; 13: e1–19. https://doi.org/10.1371/journal.pone.0201506 PMID: 30204748

14. Knowlton G, Thomas W. Host plants of the Potato Psyllid. J Econ Entomol. 1934; 27: 547.

15. Wallis RL. Potato Psyllid selection of host plants. J Econ Entomol. 1951; 44: 815–817. https://doi.org/10.1093/jee/44.5.815

16. Butler CD, Trumble JT. The Potato Psyllid, *Bactericera cockerelli* (Sulc) (Hemiptera: Triozidae): life history, relationship to plant diseases, and management strategies. Terr Arthropod Rev. 2012; 5: 87–111. https://doi.org/10.1163/187498312X634266
17. Prager SM, Esquivel I, Trumble JT. Factors influencing host plant choice and larval performance in *Bactericera cockerelli*. PLoS One. 2014; 9. https://doi.org/10.1371/journal.pone.0094047 PMID: 24710468

18. Martin NA. Host plants of the Potato / Tomato Psyllid: A cautionary tale. Weta. 2008; 35: 12–16.

19. Cooper WR, Horton DR, Miliczky E, Wohleb CH, Waters TD. The weed link in Zebra Chip epidemiology: suitability of non-crop Solanaceae and Convolvulaceae to Potato Psyllid and "*Candidatus* Liberibacter solanacearum." Am J Potato Res. 2019; 262–271. https://doi.org/10.1007/s12230-019-09712-z

20. Swisher KD, Henne DC, Crosslin JM. Identification of a fourth haplotype of *Bactericera cockerelli* (Hemiptera: Triozidae) in the United States. J Insect Sci. 2014; 14: 1–7. https://doi.org/10.1093/jis/14.1.1

21. Swisher KD, Munyaneza JE, Crosslin JM. High resolution melting analysis of the Cytochrome Oxidase I gene identifies three haplotypes of the Potato Psyllid in the United States. Environ Entomol. 2012; 41: 1019–1028. https://doi.org/10.1603/EN12066

22. Swisher KD, Sengoda VG, Dixon J, Echegaray E, Murphy AF, Rondon SI, et al. Haplotypes of the potato psyllid, *Bactericera cockerelli*, on the wild host plant, Solanum dulcamara, in the pacific North-western United States. Am J Potato Res. 2013; 90: 570–577. https://doi.org/10.1007/s12230-013-9330-3

23. Munyaneza JE, Crosslin JM, Upton JE. Association of *Bactericera cockerelli* (Homoptera: Psyllidae) with "Zebra Chip" a New Potato Disease in Southwestern United States and Mexico. J Econ Entomol. 2007; 100: 656–663. https://doi.org/10.1603/0022-0493(2007)100[656:aobchp]2.0.co;2 PMID: 17598522

24. Munyaneza JE, Fisher TW, Sengoda VG, Garczynski SF, Nissinen A, Lemmetty A. Association of "Candidatus Liberibacter solanacearum" with the psyllid, *Triozoa apicalis* (Hemiptera: Triozidae) in Europe. J Econ Entomol. 2010; 103: 1060–1070. https://doi.org/10.1603/ec10027 PMID: 20857712

25. Munyaneza JE, Sengoda VG, Stegmark R, Arvidsson AK, Anderbrant O, Yuvaraj JK, et al. First report of "Candidatus Liberibacter solanacearum" associated with psyllid-affected carrots in Sweden. Plant Dis. 2012; 96: 453–453. https://doi.org/10.1094/PDIS-10-11-0871 PMID: 30727133

26. Liefting LW, Weir BS, Pennycook SR, Clover GRG. "Candidatus Liberibacter solanacearum", associated with plants in the family Solanaceae. Int J Syst Evol Microbiol. 2009; 59: 2274–2276. https://doi.org/10.1099/ijs.0.007377-0 PMID: 19620372

27. Bextine B, Aguilar E, Rueda A, Caceres O, Sengoda VG, McCue KF, et al. First report of "Candidatus Liberibacter solanacearum" on tomato in El Salvador. Plant Dis. 2013; 97: 1245. https://doi.org/10.1094/PDIS-03-13-0248-PDN PMID: 30722460

28. Munyaneza JE, Sengoda VG, Aguilar E, Bextine B, McCue KF. First report of "Candidatus Liberibacter solanacearum" on pepper in Honduras. Plant Dis. 2014; 98: 154. https://doi.org/10.1094/PDIS-06-13-0598-PDN PMID: 30708588

29. Munyaneza JE, Sengoda VG, Aguilar E, Bextine BR, McCue KF. First report of "Candidatus Liberibacter solanacearum" infecting eggplant in Honduras. Plant Dis. 2013; 97: 1654. https://doi.org/10.1094/PDIS-06-13-0641-PDN PMID: 30716844

30. Aguilar E, Sengoda VG, Bextine B, McCue KF, Munyaneza JE. First report of "Candidatus Liberibacter solanacearum" on tobacco in Honduras. Plant Dis. 2013; 97: 1376. https://doi.org/10.1094/PDIS-04-13-0453-PDN PMID: 30722168

31. Munyaneza JE, Sengoda VG, Aguilar E, Bextine B, McCue KF. First report of "Candidatus Liberibacter solanacearum" associated with psyllid-infested tobacco in Nicaragua. Plant Dis. 2013; 97: 1244. https://doi.org/10.1094/PDIS-03-13-0247-PDN PMID: 30722439

32. Johnson D. Canadian Potato Psyllid and Zebra Chip Monitoring Network Report. News Lett Can Potato Psyllid Zebra Chip Monit Netw. 2017; 1–19. https://doi.org/10.15713/ins.mmj.3

33. Castillo Carrillo C, Fu Z, Burckhardt D. First record of the Tomato Potato Psyllid *Bactericera cockerelli* from South America. Bull Insectology. 2019; 72: 85–91.

34. Swisher Grimm KD, Garczynski SF. Identification of a new haplotype of *Candidatus* Liberibacter solanacearum* in *Solanum tuberosum*. Plant Dis. 2019; 103: 468–474. https://doi.org/10.1094/PDIS-06-18-0937-RE PMID: 30632471

35. Nelson WR, Fisher TW, Munyaneza JE. Haplotypes of "Candidatus Liberibacter solanacearum" suggest long-standing separation. Eur J Plant Pathol. 2011; 130: 5–12. https://doi.org/10.1007/s10658-010-9737-3

36. Wen A, Mallik I, Alvarado VY, Pasche JS, Wang X, Li W, et al. Detection, distribution, and genetic variability of *Candidatus Liberibacter* species associated with zebra complex disease of potato in North America. Plant Dis. 2009; 93: 1102–1115. https://doi.org/10.1094/PDIS-93-11-1102 PMID: 30754588

37. Borges KM, Cooper WR, Garczynski SF, Thinakaran J, Jensen AS, Horton DR, et al. "Candidatus Liberibacter solanacearum" associated with the psyllid, *Bactericera maculipennis* (Hemiptera: Triozidae). Environ Entomol. 2017; 46: 210–216. https://doi.org/10.1093/ee/nvw174 PMID: 28108600
38. CNAS. Economic Impacts of Zebra Chip on the Texas Potato Industry. Texas; 2006.

39. Greenway G, Rondon SI. Economic impacts of Zebra Chip in Idaho, Oregon, and Washington. Am J Potato Res. 2018; 95: 362–367.

40. Soliman T. Economic impact assessment of invasive plant pests in the European Union. Wageningen University, Wageningen, Netherlands. 2012.

41. DEFRA. Pest Alert: Bactericera cockerelli. 2017. Available: https://planthealthportal.defra.gov.uk/ph-api/pests/27077/notices/6682/documents/4183/document

42. EPPO. Bactericera cockerelli. EPPO Bull. 2013; 43: 202–208. https://doi.org/10.1111/epp.12044

43. Burckhardt D, Ouvrard D. A revised classification of the jumping plant-lice (Hemiptera: Psylloidea). Zoon taxa. 2012; 34: 1–34.

44. Bell JR, Alderson L, Izera D, Kruger T, Parker S, Pickup J, et al. Long-term phenological trends, species accumulation rates, aphid traits and climate: Five decades of change in migrating aphids. J Anim Ecol. 2015; 84: 21–34. https://doi.org/10.1111/1365-2656.12282 PMID: 25123260

45. Sjolund MJ, Ouvrard D, Kenyon D, Hightet F. Developing an RT-PCR assay for the identification of psyllid species. Proc Crop Prot North Britain. 2016; 279–282.

46. Percy DM. Radiation, diversity, and host-plant interactions among island and continental legume-feeding psyllids. Evolution (N Y). 2003; 57: 2540–2556. https://doi.org/10.1111/j.0014-3820.2003.tb01498.x PMID: 14686530

47. Peccoud J, Labonne G, Sauvion N. Molecular test to assign individuals within the Cacopsylla pruni complex. PLoS One. 2013; 8: e1–8. https://doi.org/10.1371/journal.pone.0072454 PMID: 23977301

48. EPPO. PM 7/129 (1) DNA barcoding as an identification tool for a number of regulated pests. EPPO Bull. 2016; 46: 501–537. https://doi.org/10.1111/epp.12344

49. Thompson JD, Higgins DG, Gibson TJ. CLUSTAL W: Improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994; 22: 4673–4680. https://doi.org/10.1093/nar/22.22.4673 PMID: 7984417

50. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990; 215: 403–410. https://doi.org/10.1016/S0022-2836(05)80360-2 PMID: 2231712

51. Clark K, Karsch-Mizrachi I, Lipman DJ, Ostell J, Sayers EW. GenBank. Nucleic Acids Res. 2015/11/20. 2016; 44: D67–D72. https://doi.org/10.1093/nar/gkv1276 PMID: 26590407

52. Untergasser A, Cutcutache I, Koressaar T, Ye J, Faircloth B, Remm M, et al. Primer3—new capabilities and interfaces. Nucleic Acids Res. 2012; 40: e115. Available: http://www.ncbi.nlm.nih.gov/pmc/articles/PMC3424584/ doi: 10.1093/nar/gks596 PMID: 22730293

53. Ye J, Coulouris G, Zaretskaya I, Cutcutache I, Rozen S, Madden T. Primer-BLAST: A tool to design target-specific primers for polymerase chain reaction. BMC Bioinformatics. 2012; 13: 134.

54. Liu D, Trumble JT, Stouthamer R. Genetic differentiation between eastern populations and recent introductions of Potato Psyllid (Bactericera cockerelli) into western North America. Entomol Exp Appl. 2006; 118: 177–183. https://doi.org/10.1111/j.1570-7458.2006.00383.x

55. European Network of GMO Laboratories. Definition of minimum performance requirements for analytical methods of GMO testing. 2015. https://doi.org/10.2788/65827

56. Li D, Fan Q-H, Waite DW, Gunawardana D, George S, Kumarasinghe L. Development and validation of a Real-Time PCR assay for rapid detection of Two-Spotted Spider Mite, *Tetranychus urticae* (Acari: Tetranychidae). PLoS One. 2015; 10: e0131887. https://doi.org/10.1371/journal.pone.0131887 PMID: 26147599

57. Dhami MK, Dsouza M, Waite DW, Anderson D, Li D. Real-Time PCR assay for the identification of the Brown Marmorated Stink Bug (Halyomorpha halys). Front Mol Biosci. 2016; 3: 1–6. https://doi.org/10.3389/fmolb.2016.00001

58. Virgilio M, Jordaens K, Breman FC, Backeljau T, de Meyer M. Identifying insects with incomplete DNA barcode libraries, African fruit flies (Diptera: Tephritidae) as a test case. PLoS One. 2012; 7. https://doi.org/10.1371/journal.pone.0031581 PMID: 22359600