# Guiding classical biological control of an invasive mealybug using integrative taxonomy

Aleixandre Beltra, Pia Addison, Juan Antonio Avalos, Didier Crochard, Ferran Garcia Mari, Emilio Guerrieri, Jan H. Giliomee, Thibaut Malausa, Cristina Navarro Campos, Ferran Palero, et al.

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

## To cite this version:

Aleixandre Beltra, Pia Addison, Juan Antonio Avalos, Didier Crochard, Ferran Garcia Mari, et al. Guiding classical biological control of an invasive mealybug using integrative taxonomy. PLoS ONE, 2015, 10 (6), 10.1371/journal.pone.0128685. [hal-02631542](https://hal.inrae.fr/hal-02631542)

HAL Id: hal-02631542

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

Submitted on 27 May 2020

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published 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_2.png" /></figure>
RESEARCH ARTICLE

# Guiding Classical Biological Control of an Invasive Mealybug Using Integrative Taxonomy

Aleixandre Beltrà¹, Pia Addison², Juan Antonio Ávalos¹, Didier Crochard³, Ferran Garcia-Mari¹, Emilio Guerrieri⁴, Jan H. Giliomee⁵, Thibaut Malausa³, Cristina Navarro-Campos¹*, Ferran Palero³, Antonia Soto¹*

1 Institut Agroforestal Mediterrani, Universitat Politècnica de València, València, Spain, 2 Department of Conservation Ecology and Entomology, Stellenbosch University, Stellenbosch, South Africa, 3 INRA, Univ. Nice Sophia Antipolis, CNRS, UMR 1355-7254 Institut Sophia Agrobiotech, 06900 Sophia Antipolis, France,

4 Istituto per la Protezione Sostenibile delle Piante, Consiglio Nazionale delle Ricerche, Portici, Italy,

5 Centre for Invasion Biology, Department of Botany & Zoology, Stellenbosch University, Stellenbosch, South Africa

* Current address: R&D Department, Biobest Belgium N.V., Westerlo, Belgium

* asoto@eaf.upv.es

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

CrossMark
click for updates

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

## OPEN ACCESS

**Citation:** Beltrà A, Addison P, Ávalos JA, Crochard D, Garcia-Mari F, Guerrieri E, et al. (2015) Guiding Classical Biological Control of an Invasive Mealybug Using Integrative Taxonomy. PLoS ONE 10(6): e0128685. doi:10.1371/journal.pone.0128685

**Academic Editor:** Heike Lutermann, University of Pretoria, SOUTH AFRICA

**Received:** February 24, 2015

**Accepted:** April 29, 2015

**Published:** June 5, 2015

**Copyright:** © 2015 Beltrà et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** The analysed sequences were deposited in Genbank under accession numbers KP771926-KP771972. Mealybug slides are available at the Polytechnic University of Valencia (Valencia, Spain). Parasitoid slides are deposited at Università degli Studi di Napoli Federico II (Portici, Italy).

**Funding:** This work was supported by the European Union Seventh Framework Programme FP7-IRSES #269196 "IPRABIO" and FP7-IAPP #324475 "COLBICS" and FP7-IRSES #612566 "BIOMODICS". The funders had no role in the study design, data

## Abstract

*Delottococcus aberiae* De Lotto (Hemiptera: Pseudococcidae) is a mealybug of Southern African origin that has recently been introduced into Eastern Spain. It causes severe distortions on young citrus fruits and represents a growing threat to Mediterranean citrus production. So far, biological control has proven unsatisfactory due to the absence of efficient natural enemies in Spain. Hence, the management of this pest currently relies only on chemical control. The introduction of natural enemies of *D. aberiae* from the native area of the pest represents a sustainable and economically viable alternative to reduce the risks linked to pesticide applications. Since biological control of mealybugs has been traditionally challenged by taxonomic misidentification, an intensive survey of *Delottococcus* spp. and their associated parasitoids in South Africa was required as a first step towards a classical biological control programme. Combining morphological and molecular characterization (integrative taxonomy) a total of nine mealybug species were identified in this study, including three species of *Delottococcus*. Different populations of *D. aberiae* were found on wild olive trees, in citrus orchards and on plants of *Chrysanthemoides monilifera*, showing intra-specific divergences according to their host plants. Interestingly, the invasive mealybug populations from Spanish orchards clustered together with the population on citrus from Limpopo Province (South Africa), sharing COI haplotypes. This result pointed to an optimum location to collect natural enemies against the invasive mealybug. A total of 14 parasitoid species were recovered from *Delottococcus* spp. and identified to genus and species level, by integrating morphological and molecular data. A parasitoid belonging to the genus *Anagyrus*, collected from *D. aberiae* in citrus orchards in Limpopo, is proposed here as a good biological control agent to be introduced into Spain.
collection and analysis, decision to publish, or
preparation of the manuscript.

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

**Introduction**

Many mealybug species (Hemiptera: Pseudococcidae) are major pests that cause significant losses in crops and ornamental plants [1]. Due to their small size and cryptic behaviour, they are often unnoticed and become invasive species that spread through the international trade of fruits and ornamentals [2,3]. Population outbreaks are frequent when mealybugs are introduced into new areas without their specific natural enemies and therefore classical biological control programmes have been widely used for their management [4]. These programmes generally rely on the importation of encyrtid parasitoids from their native area, given that their high specificity enables optimum results to be achieved with low risk of parasitizing non-target hosts [4-6]. Encyrtid parasitoids have allowed for the effective control of important mealybug outbreaks, such as *Maconelicoccus hirsutus* (Green) [7] and *Paracoccus marginatus* Williams & Granara de Willink in the Caribbean [8,9] and *Rastrococcus invadens* Williams in West Africa [10,11]. Besides these successes, there are well-documented examples of biological control programmes against mealybugs that have failed [4,5,12]. According to Moore [4], the most common causes of these failures are host misidentification, hyperparasitism and low acclimation capacity. Taxonomy has proven to be crucial in mealybug biological control and the difficulties associated with morphological identification have delayed the implementation of several control programmes. For example, the misidentification of *Planococcus kenyae* (Le Pель-ley) as *P. lilacinus* (Cockerell) led to the unsuccessful introduction of several parasitoids from Southeast Asia into Kenya [5]. Similarly, the misidentification of *Phenacoccus manihoti* Matile-Ferrero led to the ineffective introduction of *Phenacoccus herreni* Cox & Williams parasitoids into West Africa [13]. Both cases were later amended through the correct identification of the target mealybugs and the introduction of host-specific parasitoids [5,14]. Taxonomic expertise is also required to identify the candidate species of parasitoids for biological control [15]. Indeed, some biological and behavioural characteristics relevant for biological control such as host preference may differ in closely-related natural enemies [16].

Morphological identification of mealybugs and encyrtids share similar difficulties: high number of undescribed species, reduced number of experienced taxonomists and presence of cryptic species. These difficulties can be addressed by applying integrative taxonomy, which combines multiple disciplines such as phylogeography, comparative anatomy, population genetics, ecology and behavioural biology to solve taxonomic problems and delimit species [17,18]. In recent years, the integration of molecular techniques for the characterization of mealybugs and encyrtids has provided a new approach for correct identification at species level [19-25]. Among these techniques, DNA barcoding has been shown to be particularly useful because it allows for fast and accurate identification of previously sequenced species, in addition to the flagging of cryptic species and providing important insights into population genetics and molecular phylogenetics [26,27]. From an applied point of view, DNA barcoding can be a key tool for assessing the specific area of origin of the target pest and selecting its coevolved natural enemies [28-30].

*Delottococcus aberiae* (De Lotto) (Hemiptera: Pseudococcidae) is an invasive mealybug from Southern Africa that was detected in Eastern Spain in 2009 causing serious damage in citrus orchards [31]. This polyphagous species feeds on tropical and subtropical crops such as coffee, guava, citrus, persimmon, and pear [32-34]. Like other species of mealybugs damaging citrus, *D. aberiae* reduces plant vigour and excretes honeydew which promotes the growth of sooty-mould fungi. However, when *D. aberiae* develops on young fruits it causes severe distortions leading to major crop losses. Since its establishment, surveys in citrus orchards revealed the absence of parasitoids and the inadequacy of generalist predators for controlling outbreaks of *D. aberiae* in spring and summer [35]. Thus, the management of *D. aberiae* still relies on the
application of broad-spectrum insecticides such as chlorpyrifos. The economic and environ-
mental impacts of chemical control, and its potential interference with the biological control of
other citrus pests, compelled us to develop additional management strategies. Among them,
classical biological control, which has been successfully used against invasive scale insects in
Spanish citrus [36,37], appeared feasible and affordable.

The success of an effective biological control program for *D. aberiae* could be challenged by the misidentification of the mealybug and/or its natural enemies. Indeed, Miller and Giliomee [34] suggested the presence of cryptic species within morphospecies of *D. aberiae* and only one parasitoid has been tentatively associated with *D. aberiae* [38]. Therefore, this study was a first step towards classical biological control of *D. aberiae*, in which we surveyed mealybug populations in Spanish and South African citrus orchards and natural ecosystems to characterize *Delottococcus* spp. and their parasitoids. Specifically, we used integrative taxonomy to: i) discriminate *D. aberiae* and closely related species; ii) estimate the intraspecific genetic distances among populations of *D. aberiae*; and iii) identify candidate parasitoids for biological control of *Delottococcus* species.

**Materials and Methods**

Mealybug and parasitoid survey

A total of 25 sites were surveyed across Eastern Spain and the South African provinces of the Western Cape, Mpumalanga and Limpopo, between 2012 and 2014 (Table 1). Sampling sites comprised natural ecosystems, citrus orchards and botanical gardens. Some of these sites were selected following previous records of *Delottococcus* spp. [34]. Mealybugs were collected and placed in small plastic vials with 70% ethanol and preserved at -20°C for molecular identification. When populations of *Delottococcus* were tentatively identified in the field, mealybug infested twigs and leaves were collected for two hours. The material was placed into sampling bags and examined in the laboratory with a dissecting microscope. Mummified mealybugs were isolated to 3 x 0.8 cm glass vials covered with a cotton plug and kept in the laboratory at room temperature (20 ± 5°C) and natural photoperiod. Vials were checked daily for parasitoid emergence. Upon emergence, 70% ethanol was added into the tube to kill adult parasitoids and vials were stored at -20°C.

All samplings were carried out on private land and non-protected areas, except those con-
ducted in Jonkershoek Nature Reserve which were under permit number 0056-AAA041 00028
from Cape-Nature (Table 1). All private properties were surveyed under permission of their
owners. No specific permission was required for sampling insects in other areas. The samplings
did not involve endangered or protected species.

**Morphological and molecular characterization of insects**

The characterization of mealybugs and parasitoids was carried out in the following steps: i) morphological identification of all the mealybug populations surveyed examining five individuals of each population; ii) molecular analysis of *Delottococcus* and closely-related genera; iii) molecular analysis of parasitoid specimens emerged from *Delottococcus* populations; iv) morphological identification of all the mealybug and parasitoid specimens whose DNA was successfully sequenced.

Mealybug morphological identification was performed according to the procedures de-
scribed by Williams & Granara de Willink [39] with modifications. A small ventral incision
was cut behind the hind leg of each specimen and it was heated to 80°C in KOH for 25 minutes
and washed in distilled water for 15 minutes. Once the body contents were removed, the speci-
men was stained for one hour in a saturated solution of fuchsin in a 1:1:1 mixture of water,
<table><caption>Table 1. Collection localities, mealybugs and parasitoids surveyed from South Africa and Spain.</caption><thead><tr><th rowspan="2">Province</th><th rowspan="2">City</th><th colspan="5">Sampling site</th><th colspan="2">Mealybugs</th><th colspan="2">Parasitoids</th></tr><tr><th>Host plant</th><th>GPS coordinates</th><th>Protection status</th><th>Collection date</th><th>Population</th><th>Species</th><th>Species</th><th>Species</th></tr></thead><tbody><tr><td rowspan="12">Western Cape</td><td>Stellenbosch</td><td>Olea europaea subsp. africana</td><td>-33.945104,18.842711</td><td>Non protected area</td><td>25/01/2012</td><td>1</td><td>Delottococcus aberiae</td><td>Anagyrus aurantifrons (2)</td><td></td></tr><tr><td>Stellenbosch</td><td>Olea europaea subsp. africana</td><td>-33.942719,18.859448</td><td>Non protected area</td><td>28/01/2012</td><td>2</td><td>Delottococcus aberiae</td><td>Lamennaisia sp. (5)</td><td></td></tr><tr><td>Stellenbosch</td><td>Olea europaea subsp. africana</td><td>-33.933266,18.886614</td><td>Non protected area</td><td>2/02/2012</td><td>3</td><td>Delottococcus aberiae</td><td>Pachyneuron sp. (1)</td><td></td></tr><tr><td>Stellenbosch</td><td>Olea europaea subsp. africana</td><td>-33.940886, 18.858011</td><td>Non protected area</td><td>8/02/2012</td><td>4</td><td>Delottococcus aberiae (6)</td><td>Aenasius comperei (3)</td><td></td></tr><tr><td>Stellenbosch</td><td>Olea europaea subsp. africana</td><td>-33.93834,18.879361</td><td>Non protected area</td><td>9/02/2012</td><td>5</td><td>Delottococcus aberiae</td><td>Anagyrus aurantifrons (2), Cynipoidea (1), Lamennaisia sp. (1), Proctotrupoidea (2)</td><td></td></tr><tr><td>Stellenbosch</td><td>Olea europaea subsp. africana</td><td>-33.93729,18.875188</td><td>Non protected area</td><td>26/02/2012</td><td>6</td><td>Delottococcus aberiae</td><td>Anagyrus aurantifrons (5)</td><td></td></tr><tr><td>Stellenbosch</td><td>Olea europaea subsp. africana</td><td>-33.929616, 18.851437</td><td>Non protected area</td><td>12/03/2012</td><td>7</td><td>Delottococcus aberiae (6)</td><td>-</td><td></td></tr><tr><td>Kirstenbosch</td><td>Olea europaea subsp. africana</td><td>-33.986828,18.435936</td><td>Non protected area</td><td>3/02/2012</td><td>8</td><td>Delottococcus aberiae</td><td>-</td><td></td></tr><tr><td>Kirstenbosch</td><td>Olea europaea subsp. africana</td><td>-33.986828,18.435936</td><td>Non protected area</td><td>28/02/2012</td><td>9</td><td>Delottococcus aberiae</td><td>-</td><td></td></tr><tr><td>Jonkershoek</td><td>Olea europaea subsp. africana</td><td>-33.968122,18.933896</td><td>Nature reserve</td><td>31/01/2012</td><td>10</td><td>Delottococcus aberiae</td><td>-</td><td></td></tr><tr><td>Paarl</td><td>Olea europaea subsp. africana</td><td>-33.762416,18.933198</td><td>Non protected area</td><td>16/02/2012</td><td>11</td><td>Delottococcus aberiae</td><td>-</td><td></td></tr><tr><td>Wellington</td><td>Olea europaea subsp. africana</td><td>-33.7775, 18.951111</td><td>Non protected area</td><td>28/02/2012</td><td>12</td><td>Delottococcus aberiae</td><td>-</td><td></td></tr><tr><td>Citrusdal</td><td>Citrus sinensis</td><td>-32.61393,18.709717</td><td>Private land</td><td>26/01/2012</td><td>13</td><td>Planococcus citri</td><td>-</td><td></td></tr><tr><td>Citrusdal</td><td>Citrus sinensis</td><td>-32.41127,18.790741</td><td>Private land</td><td>26/01/2012</td><td>14</td><td>Planococcus citri</td><td>-</td><td></td></tr><tr><td>Stellenbosch</td><td>Citrus sinensis</td><td>-33.944624,18.870885</td><td>Non protected area</td><td>12/02/2012</td><td>15</td><td>Pseudococcus longispinus; Planococcus citri</td><td>-</td><td></td></tr><tr><td>Vermont</td><td>Chrysanthemoides monilifera</td><td>-34.415478,19.177537</td><td>Private land</td><td>1/01/2011</td><td>16</td><td>Delottococcus aberiae (7)</td><td>-</td><td></td></tr><tr><td>Vermont</td><td>Chrysanthemoides monilifera</td><td>-34.415478,19.177537</td><td>Private land</td><td>28/02/2012</td><td>17</td><td>Vryburgia transvaalensis (4)</td><td>Coccophagus sp. (1)</td><td></td></tr><tr><td>Jonkershoek</td><td>Phylica pubescens</td><td>-33.968122,18.933896</td><td>Nature reserve</td><td>5/02/2012</td><td>18</td><td>Delottococcus phylicus (6)</td><td>Anagyrus sp. 2 (1), Chartocerus sp. 1 (8), Rhopus notuis (4), Anagyrus sp. 1 (1)</td><td></td></tr><tr><td>Kirstenbosch</td><td>Phylica pubescens</td><td>-33.982513,18.453941</td><td>Non protected area</td><td>28/02/2012</td><td>19</td><td>Delottococcus phylicus</td><td>Anagyrus sp. 1 (2), Proctotrupoidea (1)</td><td></td></tr><tr><td>Kirstenbosch</td><td>Leucadendron argenteum</td><td>-33.982513,18.453941</td><td>Non protected area</td><td>28/02/2012</td><td>20</td><td>Delottococcus confusus (6)</td><td>Dendrocerus sp. (1), Prochyloneurus sp. (1)</td><td></td></tr><tr><td>Porterville</td><td>Protea magnifica</td><td>-32.931433,19.040717</td><td>Private land</td><td>6/03/2012</td><td>21</td><td>Delottococcus confusus (6)</td><td>Chartocerus sp. 1 (1), Chartocerus sp. 2 (1)</td><td></td></tr></tbody></table>

(Continued)
<table><caption>**Table 1.** (Continued)</caption><thead><tr><th rowspan="2">Province</th><th rowspan="2">City</th><th colspan="5">Sampling site</th><th colspan="2">Mealybugs</th><th>Parasitoids</th></tr><tr><th>Host plant</th><th>GPS coordinates</th><th>Protection status</th><th>Collection date</th><th>Population</th><th>Species</th><th>Species</th></tr></thead><tbody><tr><td rowspan="6">Mpumalanga</td><td rowspan="2">Nelspruit</td><td>Citrus sinensis</td><td>-25.4759,31.003375</td><td>Private land</td><td>1/03/2012</td><td>22</td><td>Paracoccus burnerae (6)</td><td>-</td></tr><tr><td>Citrus spp.</td><td>-25.435485,30.970631</td><td>Private land</td><td>1/03/2012</td><td>23</td><td>Paracoccus burnerae</td><td>-</td></tr><tr><td rowspan="3">Nelspruit</td><td>Citrus spp.</td><td>-25.435485,30.970631</td><td>Private land</td><td>2/03/2012</td><td>24</td><td>Ferrisia virgata</td><td>-</td></tr><tr><td>Citrus spp.</td><td>-25.435485,30.970631</td><td>Private land</td><td>2/03/2012</td><td>25</td><td>Planococcus citri</td><td>-</td></tr><tr><td>Citrus spp.</td><td>-25.435485,30.970631</td><td>Private land</td><td>2/03/2012</td><td>26</td><td>Paracoccus burnerae</td><td>-</td></tr><tr><td>Nelspruit</td><td>Olea europaea</td><td>-25.462495,30.94677</td><td>Non protected area</td><td>2/03/2012</td><td>27</td><td>Nairobia bifrons</td><td>-</td></tr><tr><td rowspan="4">Limpopo</td><td rowspan="2">Letsitele</td><td>Citrus x paradisi</td><td>-23.853205,30.388875</td><td>Private land</td><td>29/01/2014</td><td>28</td><td>Delottococcus aberiae (2)</td><td>-</td></tr><tr><td>Citrus x paradisi</td><td>-23.848194,30.401205</td><td>Private land</td><td>29/01/2014</td><td>29</td><td>Delottococcus aberiae</td><td>-</td></tr><tr><td rowspan="2">Letsitele</td><td>Citrus x paradisi</td><td>-23.798969,30.436491</td><td>Private land</td><td>29/01/2014</td><td>30</td><td>Delottococcus aberiae (3)</td><td>Anagyrus sp. 1 (3)</td></tr><tr><td>Citrus x paradisi</td><td>-23.839636,30.453972</td><td>Private land</td><td>29/01/2014</td><td>31</td><td>Delottococcus aberiae (3)</td><td>-</td></tr><tr><td>Comunitat Valenciana</td><td>Quart de les Valls</td><td>Citrus reticulata x Citrus sinensis</td><td>39.745544, -0.296638</td><td>Private land</td><td>16/07/2012</td><td>32</td><td>Delottococcus aberiae (5)</td><td>-</td></tr></tbody></table>

Number of individuals sequenced (n).

doi:10.1371/journal.pone.0128685.t001

lactic acid and glycerol. Following this, the specimen was transferred into acetic acid for one hour to fix the dye and then moved into clove oil for one hour. Insects were slide-mounted in Canada balsam. Mealybugs were identified to species level using available taxonomic keys [32,34,40–43]. The slides are available for examination at the Polytechnic University of Valencia (Valencia, Spain) for any interested researcher who visited this location.

Parasitoid identification was carried out as follows: Every specimen sequenced was placed for 24h in a xylene and absolute ethanol 1:1 solution, transferred into amyl acetate for 24 h, dried up in amyl acetate until evaporation and mounted on card with water-soluble glue. Selected card-mounted specimens were slide-mounted following Noyes [44]. In brief, wings were removed from card-mounted specimens and mounted in Canada balsam with no further passage. The remaining insect was transferred for 10 minutes to a KOH 10% solution at 100°C, moved to acetic acid for five minutes at room temperature, rinsed with distilled water and dehydrated in a progressive series of ethanol (from 70% to 100%). Once in absolute ethanol a drop of clove oil was added waiting the complete evaporation of ethanol. The insect was then transferred to the slide in a Canada balsam drop, dissected and heated at 100° overnight. All slides are deposited at Università degli Studi di Napoli Federico II (Portici, Italy) and are available for examination for any interested researcher who visited this location. Parasitoids were identified to species or genus level with the aid of available keys and by comparing them with type material or material authoritatively identified and preserved at the Natural History Museum of London (UK).

Ten mealybug populations and 44 parasitoid specimens were selected for further molecular analyses (Table 1). DNA was extracted without crushing the specimen body using the DNeasy Tissue Kit (QIAGEN, Hilden, Germany) [20] for mealybugs and the prepGEM Insect Kit (ZyGEM, Lane Hamilton, New Zealand) for parasitoids. DNA was amplified from three different genes: the cytochrome oxidase subunit 1 (mtDNA), the 28S ribosomal gene (nuclear
<table><caption>Table 2. PCR primers used in this study to amplify mealybug and parasitoid DNA.</caption><thead><tr><th>Target group</th><th>Locus</th><th>Primer names</th><th>Primer sequences</th><th>Annealing temperature</th><th>PCR product length (bp)</th><th>Reference</th></tr></thead><tbody><tr><td rowspan="3">Pseudococcidae</td><td>COI</td><td>PcoF1-LepR1</td><td>CCTTCAACTAATCATAAAAATATYAG / TAAACTTCTGGATGTCCAAAAAATCA</td><td>54°C</td><td>~700bp</td><td>[61,62]</td></tr><tr><td>28S (D10)</td><td>S3690-A4394</td><td>GAGAGTTMAASAGTACGTGAAAC / TCGGARGGAACCAGCTACTA</td><td>58°C</td><td>~800bp</td><td>[63]</td></tr><tr><td>rpS15-16ST</td><td>leuA-U16S</td><td>GTATCTAGAGGNATHCAYCARGAYGGNG / GCCGTMCGACTWGCATGTG</td><td>60°C</td><td>~1000bp</td><td>[64]</td></tr><tr><td rowspan="2">Chalcidoidea</td><td>COI</td><td>LCO1490-HCO2198</td><td>GGTCAACAAATCATAAGAGATTGG / TAAACTTCAGGGTGACCAAAAAAATCA</td><td>50°C</td><td>~700bp</td><td>[65]</td></tr><tr><td>28S (D2)</td><td>28S-D2 (F)- 28S-D2 (R)</td><td>CGTGTTGCTTGATAGTCAGC / TCAAGACGGGTCCTGAAAGT</td><td>58°C</td><td>~600bp</td><td>[66]</td></tr></tbody></table>

doi:10.1371/journal.pone.0128685.t002

genome), and the 16s-leuA region (from the genome of the mealybug endosymbiont *Tremblaya princeps*) (Table 2). PCR was performed with a 23µl reaction mixture and 2µl of diluted DNA (1–20 ng). The reagent concentrations were 12.5µl of 1X QIAGEN Multiplex PCR buffer and 0.2µM of each primer (primers in Table 2). PCR was carried out as follows: initial denaturation at 95°C for 15 minutes, followed by 35 cycles of denaturation at 94°C for 30s, annealing for 90s at a temperature of 50°C–60°C, depending on the primers (Table 2), elongation at 72°C for 60s, followed by a final extension at 72°C for 10 minutes. The final products were separated by electrophoresis with the QIAxcel Advanced System (QIAGEN, Hilden, Germany) for quality checking. PCR products were then sequenced in both directions using an ABI 3130XL automatic sequencer (Applied Biosystems, Foster City, CA, USA) at Genoscreen (Lille, France) or Beckman Genomics (Takeley, United Kingdom). Consensus sequences were generated and analysed with Seqscape v2.5 software (Applied Biosystems), and alignments were manually edited with Bioedit [45]. When a sequence of a specimen displayed genetic variation at one or more site(s), it was considered as a new haplotype. The analysed sequences were deposited in Genbank under accession numbers as shown in S1 and S2 Tables.

## Mealybug phylogenetic analysis and intraspecific variability

In order to carry out the phylogenetic analysis, DNA sequences from representative species of other Pseudococcidae were obtained in our laboratory by sequencing or from public databases (S1 Table). Alignments of the sequence data-sets were conducted using the program Muscle v3.6 (Edgar, 2004) with default parameters. To avoid alignment ambiguity, gaps and hyper-variable regions were excluded using GBlocks [46] with the following parameters: minimum number of sequences for a conserved or flanking position: 32, maximum number of contiguous non-conserved positions: 8, minimum length of a block: 10, and allowed gap positions: with half. Single-gene alignments were then concatenated and the best-fit model of DNA evolution was selected in MEGA6 [47]. Models with the lowest BIC (Bayesian Information Criterion) scores are considered to better describe the DNA substitution pattern of our alignment. Non-uniformity of evolutionary rates among sites was modelled by using a discrete Gamma distribution (+G) with 5 rate categories and by assuming that a certain fraction of sites are evolutionarily invariable (+I). After selecting for the best-fit DNA substitution model, Bayesian inference was applied using the BEAST software [48] to infer phylogenetic relationships among samples. Two independent runs starting from a random tree, using estimated base frequencies under the best-fit model and a Yule tree prior were used. Markov chains were run for 10,000,000 generations, sampling every 1,000th tree. All MCMC runs were assessed using Tracer v1.5, the graphical tool for visualization and diagnostics of MCMC output, and with a
10% burn-in. Finally, the sample of trees obtained from the MCMC runs after discarding the burn-in was summarized onto a single consensus tree using TreeAnnotator [48].

The cytochrome oxidase subunit 1 (COI) is commonly used in DNA barcoding studies to distinguish between species and it has been shown to be particularly useful for mealybugs [20,21,49]. Therefore, the COI-gene intraspecific divergence among mealybug populations collected on different hosts was also estimated through Maximum Composite Likelihood [50]. Estimates of intraspecific divergence and the corresponding standard errors were obtained using MEGA6 [47].

# Results

## Mealybug species identified

A total of nine different mealybug species were identified from the 24 sites surveyed in South Africa (Table 1). Five mealybug species were collected from citrus orchards: *D. aberiae*, *Plano-coccus citri* (Risso), *Pseudococcus longispinus* (Targioni Tozzetti), *Paracoccus burnerae* (Brain) and *Ferrisia virgata* (Cockerell). Two species were recovered on olive trees: *D. aberiae* and *Nairobi bifrons* De Lotto, even though these species did not co-occur in the same location. *Delottococcus aberiae* was present on wild olive trees (*Olea europaea* subsp. *africana*) in Western Cape natural ecosystems, while *N. bifrons* was occasionally found on *Olea europaea* in Nelsonspruit botanical gardens. Moreover, other mealybug species were collected from different host plants: *Vryburgia transvaalensis* (Brain) and *D. aberiae* on *Chrysanthemoides monilifera*, *Delot-tococcus phylicus* (De Lotto) on *Phylica pubescens*; and *Delottococcus confusus* (De Lotto) on *Leucadendron argenteum* and *Protea magnifica*. New DNA sequences are provided for the species *D. aberiae*, *D. confusus*, *D. phylicus*, *P. burnerae* and *V. transvaalensis*.

## Mealybug phylogenetic analysis and intraspecific variability

The dataset for the concatenated 28S, 16S and cytochrome oxidase alignments had an initial length of 2531 bp. A total of 2353 positions were kept for further analyses after running GBlocks (92% of the original 2531 positions). The nucleotide substitution model selected in MEGA6 was TrN93+G (BIC = 11720), with the estimated alpha parameter for the gamma distribution (α = 0.058) indicating a significant heterogeneity on the DNA substitution among sites. The effective sample size for each parameter under Bayesian inference was always >200, indicating a convergence of the MCMC runs. The consensus phylogenetic tree showed a highly significant clustering of all specimens of *D. aberiae*, with a significant support for the monophyly of this clade (Fig 1). Reciprocal monophyly was found among populations of *D. aberiae* obtained from different host plants. The Spanish populations of *D. aberiae* collected from citrus clustered significantly with the South African populations from citrus in the Limpopo Province. The specimens of *D. confusus* collected in our study also present strongly supported clades that correspond to mealybugs collected from different host plants. However, and contrary to the case of *D. aberiae*, host plant and geography cannot be disentangled within our dataset for *D. confusus*. Finally, *D. phylicus* was found to cluster with species belonging to the genus *Vryburgia*.

The Maximum Composite Likelihood genetic divergence among different species of *Delot-tococcus* ranged between 5.3% and 6.4% (Table 3). Genetic divergence was estimated among populations of *D. aberiae* to further characterize intraspecific patterns according to plant hosts. Specifically, citrus populations diverged 1.1% from those on wild olive trees and 1.8% from those on *C. monilifera*. These values are much larger than those found among populations of *D. confusus* from *L. argenteum* and *P. magnifica* (divergence = 0.2%).
<figure><img src="image_5.png" /><figcaption>Fig 1. Phylogenetic relationships among the mealybug taxa surveyed in study, as revealed by a Bayesian consensus tree obtained using the BEAST package. For each sample of the dataset, the morphological identification, the host plant, the sampling site and the specimen code are given.</figcaption></figure>

doi:10.1371/journal.pone.0128685.g001

## Parasitoid species identified

Fourteen parasitoid species emerged from *Delottococcus* spp. (Table 1). Most of these parasitoids were recovered from populations of *D. aberiae*, namely: four species of Encyrtidae, *Anagyrus aurantifrons* Compere (new host record), *Anagyrus* sp. 1 (new host record), *Aenasius comperei*

<table><caption>Table 3. Estimates of evolutionary divergence over sequence pairs between groups.</caption><thead><tr><th rowspan="2">Num.</th><th rowspan="2">Mealybug</th><th rowspan="2">Host</th><th colspan="8">Num.</th></tr><tr><th>1</th><th>2</th><th>3</th><th>4</th><th>5</th><th>6</th><th>7</th><th>8</th></tr></thead><tbody><tr><td>1</td><td>Delottococcus aberiae</td><td>Olea europaea</td><td></td><td><strong>0.003</strong></td><td><strong>0.003</strong></td><td>0.005</td><td>0.006</td><td>0.009</td><td>0.014</td><td>0.010</td></tr><tr><td>2</td><td>Delottococcus aberiae</td><td>Citrus x paradisi</td><td><strong>1.1%</strong></td><td></td><td><strong>0.005</strong></td><td>0.009</td><td>0.009</td><td>0.012</td><td>0.014</td><td>0.012</td></tr><tr><td>3</td><td>Delottococcus aberiae</td><td>Chrysanthemoides monilifera</td><td><strong>1.2%</strong></td><td><strong>1.8%</strong></td><td></td><td>0.009</td><td>0.009</td><td>0.011</td><td>0.015</td><td>0.012</td></tr><tr><td>4</td><td>Delottococcus confusus</td><td>Protea magnifica</td><td>3.0%</td><td>4.8%</td><td>4.9%</td><td></td><td>0.001</td><td>0.009</td><td>0.013</td><td>0.011</td></tr><tr><td>5</td><td>Delottococcus confusus</td><td>Leucadendron argenteum</td><td>3.2%</td><td>4.7%</td><td>4.6%</td><td>0.2%</td><td></td><td>0.009</td><td>0.013</td><td>0.011</td></tr><tr><td>6</td><td>Delottococcus phylicus</td><td>Phylica pubescens</td><td>5.3%</td><td>6.4%</td><td>6.2%</td><td>5.8%</td><td>5.7%</td><td></td><td>0.010</td><td>0.012</td></tr><tr><td>7</td><td>Vryburgia transvaalensis</td><td>Chrysanthemoides monilifera</td><td>5.9%</td><td>6.1%</td><td>6.4%</td><td>5.7%</td><td>5.8%</td><td>5.0%</td><td></td><td>0.017</td></tr><tr><td>8</td><td>Paracoccus burnerae</td><td>Citrus sinensis</td><td>7.3%</td><td>7.4%</td><td>7.6%</td><td>8.1%</td><td>7.7%</td><td>8.2%</td><td>8.5%</td><td></td></tr></tbody></table>

Standard error estimate(s) are shown above the diagonal. Numbers in bold denote the estimates between populations of *Delottococcus aberiae*.

doi:10.1371/journal.pone.0128685.t003
(Kerrich), *Lamennaisia* sp. (new host record); one species of Pteromalidae, *Pachyneuron* sp.; and two species of the superfamily Proctotrupoidea and the family Cynipidae. Two species of *Anagyrus* (*Anagyrus* sp. 1 and *Anagyrus* sp. 2) were found parasitizing *D. phylicus* together with *Rhopus notius* Prinsloo (Encyrtidae) (new host records), *Chartocerus* sp. (Signiphoridae) and one species of Proctotrupoidea. Finally, *Prochiloneurus* sp. (Encyrtidae) (new host record), two species of *Chartocerus* (Signiphoridae) and *Dendrocerus* sp. (Ceraphronoidea) were collected from *D. confusus* (Table 1). All these species were processed with molecular analysis and we obtained 59 sequences characterizing these specimens (accession numbers in S2 Table).

Discussion

The success of a biocontrol programme against *D. aberiae* in Spanish citrus orchards could be impaired by the misidentification of either mealybug and/or its natural enemies. The main aim of the current research was to characterize the diversity of *D. aberiae* and closely-related species with that of their natural enemies, collected from the native area of the mealybug. A total of nine mealybug species were identified in this survey, three of them belonging to the genus *Delottococcus* (i.e. *D. aberiae*, *D. phylicus* and *D. confusus*). In the Western Cape area (SW within South Africa), *D. aberiae* was mainly found on wild olive trees (homogeneously distributed at low densities) and on the roots of the flowering shrub *Chrysanthemoides monilifera* as reported by Miller and Giliomee [34], but it could not be found in citrus orchards. However, populations of *D. aberiae* were successfully recovered from citrus orchards (*Citrus x paradisi*) in the Limpopo area (NE within South Africa) where some outbreaks had been previously detected by Citrus Research International (Table 1). The irregular distribution of *D. aberiae* in South African citrus orchards was expected, considering that in this country it is a secondary pest of citrus that can go unnoticed for years [33,34].

Molecular data on *Delottococcus* and *Vryburgia* are scarce, so our results represent an important contribution to characterize the diversity of South African mealybugs. Hardy et al. [51] proposed for the first time the existence of a South African clade composed by the genera *Diversicrus* De Lotto, *Vryburgia* De Lotto, *Lenania* De Lotto, and some species of the paraphyletic genera *Paracoccus* Ezzat & McConnell, *Paraputo* Laing, and *Erium* Cockerell. In a preliminary study, Beltrà et al. [31] also found that the introduced populations of *D. aberiae* from Spain were closely related to those of introduced *Vryburgia rimariae* Tranfaglia from France, which reinforced the idea of a South African clade. The existence of this clade is supported by our study, following an intensive survey of South African mealybugs, which included populations of several species of *Delottococcus*. The Bayesian phylogenetic tree confirmed that *Delottococcus* and *Vryburgia* are paraphyletic genera, which is in agreement with the fact that none of the characters used to define *Delottococcus* are consistently present in all the species of this genus [34,52]. These sequences will enable inexperienced taxonomists to perform precise identifications using DNA comparison and therefore contribute to the characterization of some South African mealybugs as a complement to the initial works of Pieterse et al. [53] and Sethusa et al. [54]. The new sequences could also be used in further studies to develop a multiplex PCR protocol for fast identification of citrus mealybugs in quarantine controls including the South African species *D. aberiae* and *P. burnerae*.

The intraspecific variation found among populations of *D. aberiae* ranged from 1.1% to 1.8% in the COI locus, which might not be high enough to state conclusively that the populations collected on citrus, wild olives and *C. monilifera* constitute different species. Although Park et al. [21] found an average intraspecific genetic divergence of 0.97% within scale insect species, one quarter of their species showed divergences larger than 2.0%. In another study, Rung et al. [49] reported intraspecific genetic divergences from 1.90% to 1.98% among cryptic
species of genus *Planococcus*. Nevertheless, the variation found among populations of *D. aberiae* from South Africa should not be ignored, because it is a key aspect for the collection of specific and efficient natural enemies. Parasitoids are usually adapted to local host populations and can be more effective in parasitizing local genotypes [55]. Many encyrtids show specific interactions with mealybugs and their coevolution plays an important role on their ability to overcome defensive strategies of their hosts [6,56–58]. Indeed, specific strains of the encyrtids *Anagyrus* sp. near *pseudococci* (Girault) have shown to be more effective in parasitizing specific populations of *Planococcus ficus* (Signoret) [59]. Our results showed that populations of *D. aberiae* from Spain were closest to those found in Limpopo citrus orchards, sharing identical COI haplotypes. Therefore, within the framework of a classical biological control programme, this geographic area should be considered as a first choice for collecting parasitoids to be introduced into Europe against *D. aberiae*. Furthermore, the genetic analyses also provided some insights into the possible introduction pathway of *D. aberiae*, suggesting that fruit trade could have been involved in the mealybug invasion. This is in agreement with historical records showing that citrus fruit importation is one of the most frequent pathways of introduction of scale insects into Europe and that the first record of *D. aberiae* in Spain was located close to a citrus warehouse [3,35].

Parasitoid identification in the current study greatly benefited from combining molecular and morphological data analysis. This technique was particularly useful for matching males and females of different, though closely related, parasitoid species emerging from the same host. The parasitoid complex of *Delottococcus* spp. collected in this survey consisted of 14 parasitoid species. One species, namely *A. comperei*, has already been reported from *Delottococcus* spp.

<table><caption>Table 4. Parasitoids of *Delottococcus* spp. recorded in previous works.</caption><thead><tr><td>Mealybug host</td><td>Parasitoid species</td></tr></thead><tbody><tr><td>Delottococcus sp.</td><td>Aenasius comperei [13]</td></tr><tr><td>Delottococcus sp.</td><td>Gyranusoidea klugei [67]</td></tr><tr><td>Delottococcus sp.</td><td>Gyranusoidea litura[68,69]</td></tr><tr><td>Delottococcus sp.</td><td>Pseudococcobius dolus [70]</td></tr><tr><td>Delottococcus sp.</td><td>Pseudococcobius vibex [70]</td></tr><tr><td>Delottococcus sp.</td><td>Aenasius sp. [69]</td></tr><tr><td>Delottococcus sp.</td><td>Anagyrus sp. [69]</td></tr><tr><td>Delottococcus sp.</td><td>Anagyrus nigrescens [71]</td></tr><tr><td>Delottococcus sp.</td><td>Aphycus sp. [69]</td></tr><tr><td>Delottococcus aberiae</td><td>Aenasius comperei [13]</td></tr><tr><td>Delottococcus proteae</td><td>Leptomastix dactylopii [72]</td></tr><tr><td>Delottococcus quaesitus</td><td>Aenasius comperei [13]</td></tr><tr><td>Delottococcus quaesitus</td><td>Cheiloneurus carinatus [69]</td></tr><tr><td>Delottococcus quaesitus</td><td>Coccidoxenoides perminutus [69]</td></tr><tr><td>Delottococcus quaesitus</td><td>Gyranusoidea citrina [68,69]</td></tr><tr><td>Delottococcus quaesitus</td><td>Leptomastidea usta [73]</td></tr><tr><td>Delottococcus quaesitus</td><td>Leptomastix dactylopii [72,74]</td></tr><tr><td>Delottococcus quaesitus</td><td>Aphycus sp. [69]</td></tr><tr><td>Delottococcus quaesitus</td><td>Anagyrus sp. [69]</td></tr><tr><td>Delottococcus trichiliae</td><td>Aenasius sp. [69]</td></tr><tr><td>Delottococcus trichiliae</td><td>Anagyrus sp. [69]</td></tr><tr><td>Delottococcus trichiliae</td><td>Clauselina sp. [69]</td></tr><tr><td>Delottococcus taigae</td><td>Leptomastidea bifasciata [75]</td></tr></tbody></table>

doi:10.1371/journal.pone.0128685.t004
[38,60] whilst *Anagyrus*, *Lamennaisia*, *Rhopus*, and *Prochiloneurus* spp, are new records for *Delottococcus* spp. Our data integrate previous parasitoid records already available for *Delotto-coccus* spp. (Table 4). Among the parasitoids recovered, the species of *Anagyrus*, *Aenasius*, and *Rhopus* might be of special interest because these genera have been widely used in mealybug bio-logical control. *Anagyrus* sp. 1 should be considered as the most promising biological control candidate for introduction into Eastern Spain because it parasitized the haplotypes of *D. aberiae* found in South African citrus orchards. Further research includes the detailed taxonomical characterization of this species and the completion of laboratory bioassays to assess its host specificity and performance on parasitizing Spanish haplotypes of *D. aberiae*.

**Supporting Information**

S1 Table. Complete list of mealybug samples with corresponding haplotypes: codes of voucher slide mounted specimens, species, population code (see Table 1) and Genbank accession numbers for haplotypes.

(DOCX)

S2 Table. Complete list of parasitoid samples with corresponding haplotypes: codes of voucher slide-mounted specimens, species, population code (see Table 1) and Genbank accession numbers for haplotypes.

(DOCX)

**Acknowledgments**

We want to thank Ian Millar, Douglass Miller and Gerhard Prinsloo for assistance with the morphological identification of some mealybug and parasitoid species.

**Author Contributions**

Conceived and designed the experiments: AB PA FG JG TM CN AS. Performed the experiments: AB JA DC EG JG CN FP. Analyzed the data: AB TM FP AS. Contributed reagents/materials/analysis tools: PA EG JG TM FP. Wrote the paper: AB FP EG. Revised the text: PA FG JG TM CN AS.

**References**

1. McKenzie HL. Mealybugs of California: With Taxonomy, Biology, and Control of North American Species (Homoptera, Coccoidea, Pseudococcidae). 1st ed. Berkeley: University of California Press; 1967.

2. Miller DR, Miller GL, Watson GW. Invasive species of mealybugs (Hemiptera: Pseudococcidae) and their threat to US agriculture. Proc Entomol Soc Wash. 2002; 104: 825-836.

3. Pellizzari G, Germain JF. Scales (Hemiptera, Superfamily Coccoidea). BioRisk. 2010; 4: 475-510.

4. Moore D. Agents used for biological control of mealybugs (Pseudococcidae). Biocontrol News Inf. 1988; 9: 209-225.

5. Bartlett BR. Pseudococcidae. In: Clausen CO, editor. Introduced parasites and predators of arthropod pests and weeds: a world review. Washington: Agricultural Research Service USDA; 1978. pp. 137-170.

6. Charles JG. Using parasitoids to infer a native range for the obscure mealybug, *Pseudococcus viburni*, in South America. BioControl. 2011; 56: 155-161.

7. Kairo MTK, Pollard GV, Peterkin DD, Lopez VF. Biological control of the Hibiscus Mealybug, *Maconelli-coccus hirsutus* Green (Hemiptera: Pseudococcidae) in the Caribbean. Integr Pest Manag Rev. 2000; 5: 241-254.

8. Muniappan R, Meyerdirk DE, Sengebau FM, Berringer D, Reddy GVP. Classical biological control of the papaya mealybug, *Paracoccus marginatus* (Hemiptera: Pseudococcidae) in the Republic of Palau. Florida Entomol. 2006; 89: 212-217.
9. Amarasekare KG, Mannion CM, Epsky ND. Efficiency and establishment of three introduced parasitoids of the mealybug *Paracoccus marginatus* (Hemiptera: Pseudococcidae). *Biol Control*. 2009; 51: 91–95.

10. Agricola U, Agounké D, Fischer H, Moore D. The control of *Rastrococcus invadens* Williams (Hemiptera: Pseudococcidae) in Togo by the introduction of *Gyranusoidea tebygi* Noyes (Hymenoptera: Encyrtidae). *Bull Entomol Res*. 1989; 79: 671–678.

11. Neuenschwander P, Boavida C, Bokonon-Ganta A, Gado A, Herren HR. Establishment and spread of *Gyranusoidea tebygi* Noyes and *Anagyrus mangicola* Noyes (hymenoptera: encyritidae), 2 biological-control agents released against the mango mealybug *Rastrococcus invadens* Williams (homoptera: pseudococcidae) in Africa. *Biocontrol Sci Technol*. 1994; 4: 61–69.

12. Stiling P. Why do natural enemies fail in classical biological control programs? *Am Entomol*. 1993; 39: 31–37.

13. Noyes JS, Hayat M. Oriental mealybug parasitoids of the Anagyrini (Hymenoptera: Encyrtidae). 1st ed. Wallingford: CAB International; 1994.

14. Neuenschwander P. Biological Control of the Cassava Mealybug in Africa: A Review. *Biol Control*. 2001; 21: 214–229. PMID: 11156855

15. Danks HV. Systematics in Support of Entomology. Annu Rev Entomol. 1988; 33: 271–294.

16. Rosen D, DeBach P. Systematics, morphology and biological control. Entomophaga. 1973; 18: 215–222.

17. Schlick-Steiner BC, Steiner FM, Seifert B, Stauffer C, Christian E, Crozier RH. Integrative taxonomy: a multisource approach to exploring biodiversity. *Annu Rev Entomol*. 2010; 55:421–438. doi: 10.1146/annurev-ento-112408-085432 PMID: 19737081

18. Dayrat B. Towards integrative taxonomy. *Biol J Linn Soc*. 2005; 85:407–415.

19. Triapitsyn SV, González D, Vickerman DB, Noyes JS, White EB. Morphological, biological, and molecular comparisons among the different geographical populations of *Anagyrus pseudococci* (Hymenoptera: Encyrtidae), parasitoids of *Planococcus* spp. (Hemiptera: Pseudococcidae), with notes on *Anagyrus dactylopii*. *Biol Control*. 2007; 41: 14–24.

20. Malausa T, Fenis A, Warot S, Germain J, Ris N, Prado E, et al. DNA markers to disentangle complexes of cryptic taxa in mealybugs (Hemiptera: Pseudococcidae). *J Appl Entomol*. 2011; 135: 142–155. doi: 10.1152/japplphysiol.01408.2010 PMID: 21527663

21. Park D-S, Suh S-J, Hebert PDN, Oh H-W, Hong K-J. DNA barcodes for two scale insect families, mealybugs (Hemiptera: Pseudococcidae) and armored scales (Hemiptera: Diaspididae). *Bull Entomol Res*. 2011; 101: 429–434. doi: 10.1017/S0007485310000714 PMID: 21272395

22. Zhang YZ, Si SL, Zheng JT, Li HL, Fang Y, Zhu CD, et al. DNA barcoding of endoparasitoid wasps in the genus *Anicetus* reveals high levels of host specificity (Hymenoptera: Encyrtidae). *Biol Control*. 2011; 58: 182–191.

23. Correa MCG, Germain JF, Malausa T, Zaviezo T. Molecular and morphological characterization of mealybugs (Hemiptera: Pseudococcidae) from Chilean vineyards. *Bull Entomol Res*. 2012; 102: 524–530. doi: 10.1017/S0007485312000053 PMID: 22361038

24. Pacheco da Silva VC, Bertin A, Blin A, Germain JF, Bernardi D, Rignol G, et al. Molecular and morphological identification of mealybug species (Hemiptera: Pseudococcidae) in Brazilian vineyards. *PLoS One*. 2014; 9: e103267. doi: 10.1371/journal.pone.0103267 PMID: 25062012

25. Rugman-Jones PF, Forster LD, Guerrieri E, Luck RF, Morse JG, Monti MM, et al. Taxon-specific multiplex-PCR for quick, easy, and accurate identification of encyrtid and aphelinid parasitoid species attacking soft scale insects in California citrus groves. *BioControl*. 2011; 56: 265–275.

26. Hebert PDN, Cywinska A, Ball SL, deWaard JR. Biological identifications through DNA barcodes. *Proc Biol Sci*. 2003; 270: 313–21. PMID: 12614582

27. Hajibabaei M, Singer GAC, Hebert PDN, Hickey DA. DNA barcoding: how it complements taxonomy, molecular phylogenetics and population genetics. *Trends Genet*. 2007; 23: 167–172. PMID: 17316886

28. Unruh TR, Wolley JB. Molecular methods in classical biological control. In: Bellows TS, Fisher TW, editors. Handbook of Biological Control. San Diego: Academic Press; 1999. pp. 57–85.

29. Goolsby JA, Ciomperlik MA, Legaspi BC, Legaspi JC, Wendel LE. Laboratory and field evaluation of exotic parasitoids of *Bemisia tabaci* (Gennadius) (Biotype "B") (Homoptera: Aleyrodidae) in the Lower Rio Grande Valley of Texas. *Biol Control*. 1998; 12: 127–135. PMID: 9673712

30. Kirk AA, Lacey LA, Brown JK, Ciomperlik MA, Goolsby JA, Vacek DC, et al. Variation in the Bemisia tabaci s. l. species complex (Hemiptera: Aleyrodidae) and its natural enemies leading to successful biological control of Bemisia biotype B in the USA. *Bull Entomol Res*. 2000; 90: 317–327. PMID: 11020790
31. Beltrà A, Soto A, Malausa T. Molecular and morphological characterisation of Pseudococcidae surveyed on crops and ornamental plants in Spain. Bull Entomol Res. 2012; 102: 165-172. doi: 10.1017/S0007485311000514 PMID: 22008190

32. De Lotto G. New Pseudococcidae (Homoptera: Coccoidea) from Africa. Bull Br Museum Nat Hist Entomol. British Museum (Natural History); 1961; 10: 211-238.

33. Hattingh V, Cilliers C, Bedford E. Citrus mealybugs. In: Bedford ECG, Van den Berg MA, De Villiers EA, editors. Citrus Pests in the Republic of South Africa. Nelspruit: Agricultural Research Council; 1998. pp. 112-120.

34. Miller DR, Giliomee JH. Systematic revision of the mealybug genus *Delottococcus* Cox & Ben-Dov (Hemiptera: Pseudococcidae). Afr Entomol. 2011; 19: 614-640.

35. Beltrà A, Garcia-Marí A, Soto A. El cotonet de Les Valls, *Delottococcus aberiae*, nueva plaga de los cítricos. Levante Agrícola. 2013; 419: 348-352.

36. Jacas J, Urbaneja A. Biological control in citrus in Spain: from classical to conservation biological control. In: Ciancio A, Mukerji KG, editors. Integrated Management of Arthropod Pests and Insect Borne Diseases. Dordrecht: Springer; 2010. pp. 61-72.

37. Garcia-Marí F. Plagas de los cítricos: gestión integrada en países de clima mediterráneo. 1st ed. Valencia: Phytoma; 2012.

38. Prinsloo G. A review of the encyrtid wasp tribe Aenasiini, with descriptions of new Afrotropical taxa (Hemiptera: Chalcidoidea). J Nat Hist. 1988; 22: 1465-1482.

39. Williams DJ, Granara de Willink MC. Mealybugs of Central and South America. London: CAB International; 1992.

40. Brain C. The Coccidae of South Africa. Trans R Soc South Africa. 1915; 5: 65-194.

41. De Lotto G. Observations on African mealybugs (Hemiptera: Coccoidea). Bull Br Museum Nat Hist Entomol. London; 1964; 14: 343-397.

42. De Lotto G. On some African mealybugs (Homoptera: Coccoidea: Pseudococcidae). J Entomol Soc South Afr. 1977; 40: 13-36.

43. Wakgari WM, Giliomee JH. Description of adult and immature females of six mealybug species (Hemiptera: Pseudococcidae) found on citrus in South Africa. Afr Entomol. 2005; 13: 281-332.

44. Noyes JS. Collecting and preserving chalcid wasps (Hymenoptera: Chalcidoidea). J Nat Hist. London 1982; 16: 315-334.

45. Hall T. BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp Ser. 1999; 41: 95-98.

46. Castresana J. Selection of conserved blocks from multiple alignments for their use in phylogenetic analysis. Mol Biol Evol. 2000; 17: 540-552. PMID: 10742046

47. Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: Molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013; 30: 2725-2729. doi: 10.1093/molbev/mst197 PMID: 24132122

48. Drummond AJ, Rambaut A. BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol Biol. 2007; 7: 214. PMID: 17996036

49. Rung A, Scheffer SJ, Evans G, Miller D. Molecular identification of two closely related species of mealybugs of the genus *Planococcus* (Homoptera omPseudococcidae). Ann Entomol Soc Am. 2008; 101:525-532.

50. Tamura K, Nei M, Kumar S. Prospects for inferring very large phylogenies by using the neighbor-joining method. Proc Natl Acad Sci USA. 2004; 101: 11030-11035. PMID: 15258291

51. Hardy NB, Gullan PJ, Hodgson CJ. A subfamily-level classification of mealybugs (Hemiptera: Pseudococcidae) based on integrated molecular and morphological data. Syst Entomol. 2008; 33: 51-71.

52. Cox JM, Ben-Dov Y. Planococcine mealybugs of economic importance from the Mediterranean Basin and their distinction from a new African genus (Hemiptera: Pseudococcidae). Bull Entomol Res. 1986; 76: 481-489.

53. Pieterse W, Muller DL, van Vuuren BJ. A molecular identification approach for five species of mealybug (Hemiptera: Pseudococcidae) on citrus fruit exported from South Africa. Afr Entomol. 2010; 18: 23-28.

54. Sethusa MT, Millar IM, Yessoufou K, Jacobs A, van der Bank M, van der Bank H. DNA Barcode efficacy for the identification of economically important scale insects (Hemiptera: Coccoidea) in South Africa. African Entomol. 2014; 22: 257-266.

55. Hufbauer RA, Roderick GK. Microevolution in biological control: Mechanisms, patterns, and processes. Biol Control. 2005; 35: 227-239.

56. Foldi I. Defense strategies in scale insects: Phylogenetic inference and evolutionary scenarios (Hemiptera, Coccoidea). Mem Museum Natl Hist Nat. 1997; 173: 203-230.
57. Gross P. Insect behavioral and morphological defenses against parasitoids. Annu Rev Entomol. 1993; 38: 251-273.

58. Blumberg D. Parasitoid encapsulation as a defense mechanism in the Coccoidea (Homoptera) and its importance in biological control. Biol Control. Orlando, Fla.: Academic Press; 1997; 8:225-236.

59. Sforza R, Quaglietti B, Daane KM. Survey and efficacy of the Western Europe *Anagyrus* near *pseudococci* over *Planococcus ficus*, a pest of the California vineyards. In: Mason PG, Gillespie DR, Vincent C, editors. Proceedings of the 4th International Symposium on Biological Control of arthropods; 2013 4-8 March 2013; Pucón: Chile. Pucón: Agriculture and Agri-Food Canada; 2013. pp. 46-48.

60. Kerrich G. On the classification of the Anagyrine Encyrtidae, with a revision of some of the genera (Hymenoptera: Chalcidoidea). Bull Br Mus. 1967; 20: 143-250.

61. Park DS, Leem YJ, Hahn KW, Suh SJ, Hong KJ, Oh HW. Molecular identification of mealybugs (Hemiptera: Pseudococcidae) found on Korean pears. J Econ Entomol. 2010; 103: 25-33. PMID: 20214364

62. Hebert PDN, Penton EH, Burns JM, Janzen DH, Hallwachs W. Ten species in one: DNA barcoding reveals cryptic species in the Neotropical skipper butterfly *Astraptes fulgerator*. Proc Natl Acad Sci USA. 2004; 101: 14812-14817. PMID: 15465915

63. Sequeira AS, Normark BB, Farrell BD. Evolutionary assembly of the conifer fauna: distinguishing ancient from recent associations in bark beetles. Proc Biol Sci. 2000; 267: 2359-2366. PMID: 11133024

64. Baumann L, Thao MLL, Hess JM, Johnson MW, Baumann P. The genetic properties of the primary endosymbionts of mealybugs differ from those of other endosymbionts of plant sap-sucking insects. Appl Environ Microbiol. 2002; 68: 3198-3205. PMID: 12088995

65. Folmer O, Black M, Hoeh W, Lutz R, Vrijenhoek R. DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates. Mol Mar Biol Biotechnol. 1994; 3: 294-299. PMID: 7881515

66. Heraty J, Hawks D, Kostecki JS, Carmichael A. Phylogeny and behaviour of the Gollumiellinae, a new subfamily of the ant-parasitic Eucharitidae (Hymenoptera: Chalcidoidea). Syst Entomol. 2004; 29: 544-559.

67. Prinsloo GL, Annecke DP. New Encyrtidae (Hymenoptera) from South West Africa. J Entomol Soc South Afr. 1976; 39: 190-192.

68. Prinsloo GL. The southern African species of *Gyranusoidea* Compere (Hymenoptera: Encyrtidae). J Entomol Soc South Afr. 1983; 46: 103-113.

69. Prinsloo GL. A parasitoid-host index of Afrotropical Encyrtidae (Hymenoptera: Chalcidoidea). 1st ed. Pretoria: Department of Agriculture Republic of South Africa; 1983.

70. Prinsloo GL. Revision of the mealybug parasitoids of the genus *Pseudococcobius* Timberlake (Hymenoptera: Encyrtidae) from South Africa. African Entomol. 2003; 11: 77-89.

71. Prinsloo GL. Poorly known and newly recorded species of mealybug parasitoids of the genus *Anagyrus* Howard (Hymenoptera: Encyrtidae) from South Africa. African Plant Prot. 1998; 4: 81-90.

72. Japoshvili G, Celik H. Fauna of Encyrtidae, parasitoids of coccids in Golcuk Natural Park. Entomol Hell. 2010; 19: 132-136.

73. Prinsloo GL. The Afrotropical species of *Leptomastidea* Mercet (Hymenoptera: Encyrtidae) parasitoids of mealybugs. J Hymenopt Res. 2001; 10: 145-162.

74. Anga J, Noyes J. A Revision of the Asian and Malagasy species of the genus *Leptomastix* (Hymenoptera, Encyrtidae), parasitoids of mealybugs (Homoptera: Pseudococcidae). Bull Nat Hist Museum London. 1999; 68: 93-128.

75. Trjapitzin V. Parasitic Hymenoptera of the fam. Encyrtidae of Palaearitics. 1st ed. Leningrad: Zoologicheskim Institutom Akademii Nauk SSR; 1989.