# Preliminary Investigation of Species Diversity of Rice Hopper Parasitoids in Southeast Asia

Christina Sann, Franziska Wemheuer, Alexis Beaurepaire, Rolf Daniel, Silvio Erler, Stefan Vidal

► To cite this version:

Christina Sann, Franziska Wemheuer, Alexis Beaurepaire, Rolf Daniel, Silvio Erler, et al.. Preliminary Investigation of Species Diversity of Rice Hopper Parasitoids in Southeast Asia. Insects, 2018, 9 (1), 10.3390/insects9010019. [https://doi.org/10.3390/insects9010019](https://doi.org/10.3390/insects9010019)

HAL Id: hal-02621044

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

Submitted on 26 May 2020

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

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

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

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Article

# Preliminary Investigation of Species Diversity of Rice Hopper Parasitoids in Southeast Asia

Christina Sann <fnref n="1" />,* ID, Franziska Wemheuer <fnref n="1" />,<fnref n="2" /> ID, Alexis Beaurepaire <fnref n="3" />,<fnref n="4" />, Rolf Daniel <fnref n="5" /> ID,
Silvio Erler <fnref n="3" /> ID and Stefan Vidal <fnref n="1" />,*

Received: 13 November 2017; Accepted: 26 January 2018; Published: 9 February 2018

**Abstract:** Ongoing intensification of rice production systems in Southeast Asia is causing devastating yield losses each year due to rice hoppers. Their continuing development of immunity to resistant rice varieties and pesticide applications further complicates this problem. Hence, there is a high demand for biological control agents of rice hoppers. Egg parasitoid wasps are among the most important natural enemies of rice hoppers, such as *Nilaparvata lugens* and *Nephotettix* spp. However, our knowledge of their diversity is still very limited, due to their small size and the lack of available morphological information. Classifying these parasitoids is the first step to properly understanding their role in the rice agroecosystem. We used traditional morphological identification, as well as DNA sequencing of the 28S rRNA and the COI genes, to investigate the diversity of four important hopper egg parasitoid genera in the Philippines. Parasitoids of the genera *Anagrus*, *Oligosita*, *Gonatocerus*, and *Paracentrobia* were collected in eight study landscapes located in Luzon. Our findings illustrate that characterization of species diversity using morphological and molecular analyses were concordant only for the genus *Paracentrobia*. The genera *Anagrus* and *Gonatocerus* exhibited more genetic diversity than estimated with the morphological analysis, while the opposite was observed for *Oligosita*. This is the first study investigating the molecular diversity of rice hopper parasitoids in the Philippines. More research combining morphological, behavioral, and molecular methods, as well as the establishment of a comprehensive DNA database, are urgently needed to assess the performance and suitability of these organisms as biocontrol agents.

**Keywords:** DNA barcoding; genetic diversity; hymenopteran parasitoids; *Nephotettix*; *Nilaparvata lugens*; rice

## 1. Introduction

Rice is the main food resource for more than half of the world's population [1,2]. This makes the rice production system of Asia one of the most important food production systems on Earth [3]. The brown planthopper (BPH; *Nilaparvata lugens*, Stål 1854) and the green leafhopper (GLH; *Nephottetix* spp.) are among the economically most important rice pests. These insects cause immense damage in the Asian rice paddies through xylem sap feeding, resulting in wilting, and subsequently, the death of the rice crops [4], as well as through the transfer of devastating viruses among fields [5,6].

<footnote n="1">¹ Agricultural Entomology, Department of Crop Sciences, University of Göttingen, Grisebachstr. 6, 37077 Göttingen, Germany; fwemheu@gwdg.de</footnote>

<footnote n="2">² Applied Marine and Estuarine Ecology, School of Biological, Earth and Environmental Sciences, The University of New South Wales (UNSW), Sydney, NSW 2052, Australia</footnote>

<footnote n="3">³ Molecular Ecology, Institute of Biology, Martin-Luther-University Halle-Wittenberg, Hoher Weg 4, 06099 Halle (Saale), Germany; alexis.beaurepaire@inra.fr (A.B.); silvio.erler@zoologie.uni-halle.de (S.E.)</footnote>

<footnote n="4">⁴ Abeilles et Environment, INRA PACA, 84914 Avignon CEDEX 9, France</footnote>

<footnote n="5">⁵ Department of Genomic and Applied Microbiology, Institute of Microbiology and Genetics, University of Göttingen, Grisebachstr. 8, 37077 Göttingen, Germany; rdaniel@gwdg.de</footnote>

<footnote>* Correspondence: Christina.Sann@agr.uni-goettingen.de (C.S.); svidal@gwdg.de (S.V.)</footnote>
So far, the introduction of rice varieties resistant to BPH and GLH has not been successful, as these pests rapidly adapt to the new varieties [7,8]. The misuse of pesticide applications further enhances the problem by disturbing the rice agroecosystem, which can increase the outbreak risk of hoppers [7,9,10].

As an alternative, integrated pest management (IPM) strategies can be used to maximize the pest's mortality by enhancing the role of their biological control agents. BPH and GLH are typically attacked by a wide range of natural enemies, such as spiders, predatory bugs, dragonflies, and egg parasitoid wasps from the families Mymaridae and Trichogrammatidae [11]. Among these enemies, parasitoid species are of particular interest. These organisms are very mobile, and can disperse over large distances [12]. The adults feed on pollen, nectar or the honeydew of their sap-feeding hosts, whereas their larvae develop in the hopper eggs and disrupt the hoppers' life cycles at the earliest possible stage [11]. In two previous studies on egg parasitoids of rice hoppers, egg parasitism levels of more than 60% were observed [13,14]. These findings are supported by Drechsler and Settele [15], who showed that parasitoids can play a major role in controlling hopper pests in rice agroecosystems.

Despite their importance, knowledge about species composition and diversity of parasitoid wasps in the tropics is still limited (see Nishida, Wongsiri and Wongsiri [16], as well as Gurr et al. [11]). Morphological species identification requires extensive taxonomic expertise, and is hampered by the small size (<1.5 mm) of the wasps as well as the limited amount of literature. To date, molecular information on rice hopper parasitoids from the Philippines is completely lacking. Molecular methods have become a promising tool to resolve species identities [17–19]. Different markers, such as the mitochondrial cytochrome c oxidase I (COI) and fragments of the small or large subunits of the ribosomal RNA genes (i.e., 18S or 28S rRNA genes), have been applied to construct hymenopteran and dipteran parasitoid phylogenies [20,21]. However, molecular analyses not only improve species identification, but may also reveal cryptic species [22,23]. Smith et al. [24] suggested combining barcoding with morphology and natural history. A similar conclusion was drawn by Padial and colleagues [25] in their review on an integrative taxonomy for the improvement of species discovery and description.

In a previous study, we investigated parasitoid wasps of *Nilaparvata lugens* (BPH) and *Nephotettix* spp. (GLH) in eight rice production landscapes located in Luzon, Philippines [26]. We found that BPH was parasitized by the Chalcid genera *Anagrus* and *Oligosita*, while GLH was parasitized by the chalcidoid genera *Gonatocerus* and *Paracentrobia*. In the present study, we analyzed the diversity of these genera with traditional morphological and molecular techniques. The present study provides the first molecular identification of parasitoid wasps in rice paddies in the Philippines, combined with a morphological identification based on dichotomous keys.

## 2. Materials and Methods

### 2.1. *Study System and Sampling*

This study was embedded in the project LEGATO, which focused on sustainable rice production [27]. Parasitoids were collected in eight study landscapes located in the Laguna province, Luzon, Philippines (Supplementary Figure S1). Sampling took place during the rice growing and fallow periods of the dry season from February to June 2013. In brief, rice plants of the variety Taichung Native (1) (TN1) were grown in a greenhouse for 6 weeks, trimmed to three tillers, and covered with small tubular insect cages (85 cm high, 15 cm diameter). Greenhouse cultures of BPH and GLH were reared on TN1, as previously described by Heinrichs et al. [28]. The BPH culture exclusively consisted of *Nilaparvata lugens*, while the GLH culture consisted of *Nephotettix virescens* (Distant) and *Nephotettix nigropictus* (Ståhl). Both hopper populations originated from wild individuals caught in the rice fields of the Laguna Province. Four gravid females of either BPH or GLH were released into each cage for 48 h to lay eggs. Subsequently, the plants were transferred to three different plots per study site to cover the naturally occurring landscape diversity in the Laguna province. Three plants infested with BPH eggs and three plants infested with GLH eggs were distributed randomly within each plot.
After 72 h, all plants were returned to the greenhouse and placed back inside separate insect cages. The parasitoids emerged 13–17 days after the field exposure of the eggs. Due to positive phototaxis behavior, parasitoids aggregated in small transparent glass vials, which were placed upside down on top of the sleeved cages. Adult parasitoids were collected daily within this period and immediately stored at –80 °C. The dead parasitoids were identified to genus under a stereo microscope using morphological keys [29–31]. In total, 19,455 parasitoids were collected and separated into the four genera *Paracentrobia*, *Gonatocerus*, *Oligosita*, and *Anagrus* [26].

## 2.2. Preliminary Morphological Identification

Fifty parasitoids from each of the four identified genera (*Paracentrobia*, *Gonatocerus*, *Oligosita*, and *Anagrus*) were randomly selected. The 200 parasitoids were permanently slide-mounted using Faure's medium, and identified to species level using a microscope. A specialist at the International Rice Research Institute (IRRI) helped with the identification using the key published in Heong and Hardy [30].

## 2.3. DNA Extraction and Amplification

Total genomic DNA was extracted from 267 newly selected whole single female parasitoid individuals (Supplementary Table S2) employing the Phire Animal Tissue Direct PCR Kit (Thermo Scientific, Waltham, MA, USA). Parasitoids were suspended in 20 µL TE buffer (100 mM Tris, 10 mM EDTA). To increase DNA extraction efficiency, 1 µL proteinase K (20 mg mL⁻¹) was added. Parasitoid samples were carefully homogenized with sterile micro pestles, and incubated at room temperature overnight.

Two independent gene fragments were amplified from the extracted DNA: one located on the COI subunit I (COI I, amplicon length around 670 bp) and the other one on the expansion regions D2–D3 of the 28S ribosomal subunit (28S-D2/28S-D3, amplicon length approximately 610 bp). The COI region was amplified using the primer pair HCO2198/LCO1490 [32]. The PCR reaction mixture (25 µL) contained 2.5 µL of 10-fold Ex Taq Buffer (Takara Biotechnology, co., LTD, Dalian, China), 25 mM MgCl₂, 2.5 mM of each of the four dNTPs (deoxynucleotide triphosphates), 10 µM of each primer, 1 U TaKaRa Ex Taq polymerase (Takara Biotechnology), and approximately 25 ng of parasitoid DNA. The following thermal cycling scheme was used: initial denaturation at 94 °C for 3 min, 5 cycles of denaturation at 94 °C for 45 s, annealing at 45 °C for 45 s, followed by elongation at 72 °C for 1 min, and 25 cycles of denaturation at 94 °C for 45 s, annealing at 50 °C for 45 s, followed by another elongation at 72 °C for 1 min. The final extension was carried out at 72 °C for 5 min.

The D2–D3 region of the 28S rRNA gene was amplified as described for the PCR above, using the primers D2-3549 [33] and D2-4068 [34]. The following thermal cycling scheme was used: initial denaturation at 94 °C for 3 min, 30 cycles of denaturation at 94 °C for 45 s, annealing at 58 °C for 45 s, followed by elongation at 72 °C for 1 min. The final extension was carried out at 72 °C for 6 min. Negative controls were performed by using the reaction mixture without template. PCR products were controlled for appropriate size and then purified using the peqGOLD Gel Extraction kit as recommended by the manufacturer (Peqlab, Erlangen, Germany; now VWR).

Sequencing was performed at the Göttingen Genomics Laboratory using an ABI 3730xl system and BigDye terminator chemistry version 3.1 (Thermo Fisher Scientific).

## 2.4. Data Analysis

Quality control and trimming of the forward and reverse DNA sequences were performed with Gap5 v 1.2.14-r [35]. DNA sequence alignment was performed using MEGA 6 [36]. Additional Chalcidoidea (the super family containing all four genera) sequences were used as reference material. They were retrieved by using the National Centre for Biotechnology Information (NCBI) BLAST tool [37] and the BOLD Identification System (http://www.boldsystems.org) (Supplementary Table S3). The best model to calculate parasitoid phylogeny was determined assuming partial deletion,
site coverage cut-off of 95%, and the branch swap filter set to "very strong". Maximum likelihood trees were generated under the assumption of the best fitting model with 1000 bootstraps. Finally, median-joining networks were constructed using the software NETWORK v. 4.6.1.2 [38].

The genetic diversity between and within species was estimated using the reference sequences from the superfamily Chalcidoidea (Supplementary Table S3). The appropriate models for the calculation of the pairwise molecular distances between individuals were determined with MEGA 6. The 28S rRNA distances were calculated based on the K2P + G model. COI distances were calculated based on the T92 + G model. The values were compared using a Kruskal-Wallis test, followed by a Dunn's post hoc test [39] in R version 3.2.3 [40].

To estimate, whether the genetic diversity of the parasitoids was covered by the sample size used in the analyses, we performed a rarefaction analysis (Supplementary Figures S4 and S5) in R version 3.2.3 using the packages *vegan* [41] and *drc* [42].

## 2.5. Sequence Data Deposition

All sequence data was deposited in the NCBI GenBank under accession number MG785407-MG785510, MG904875-MG904939, MG911990-MG912010.

# 3. Results

## 3.1. Morphological Analysis

In total, 200 parasitoids were morphologically determined. Only one species was found for the genera *Gonatocerus* and *Paracentrobia*, while three species were found for the genera *Anagrus* and *Oligosita*: *Anagrus flaveolus* (11 specimens), *Anagrus optabilis* (6 specimens), *Anagrus frequens* (33 specimens), *Oligosita aesopi* (41 specimens), *Oligosita naias* (7 specimens), *Oligosita shibuyae* (2 specimens), *Gonatocerus orientalis* (50 specimens), and *Paracentrobia andoi* (50 specimens) (Table 1, Supplementary Table S1).

<table><caption>**Table 1.** Species identified by the morphological analysis using 50 specimens per genera (GLH: green leafhopper, BPH: brown planthopper).</caption><thead><tr><td>Species</td><td>Host</td><td>Examined Individuals</td></tr></thead><tbody><tr><td>Paracentrobia andoi</td><td>GLH</td><td>50</td></tr><tr><td>Oligosita aesopi</td><td>BPH</td><td>41</td></tr><tr><td>Oligosita naias</td><td>BPH</td><td>7</td></tr><tr><td>Oligosita shibuyae</td><td>BPH</td><td>2</td></tr><tr><td>Anagrus frequens</td><td>BPH</td><td>33</td></tr><tr><td>Anagrus optabilis</td><td>BPH</td><td>6</td></tr><tr><td>Anagrus flaveolus</td><td>BPH</td><td>11</td></tr><tr><td>Gonatocerus orientalis</td><td>GLH</td><td>50</td></tr></tbody></table>

## 3.2. Molecular Analysis

Of the 267 specimens used for the molecular analysis, we successfully sequenced a total of 105 parasitoid samples (Supplementary Table S2). We failed to sequence 162 specimens, due to the poor quality of the extracted DNA, or a failure to extract DNA from the samples. A total of 86 (COI) and 105 (28S) sequences were used to create the maximum likelihood trees. The final dataset for comparison of the two gene fragments included the sequences from 74 parasitoid individuals.

Rarefaction analysis for the local hopper egg parasitoids genetic diversity revealed that the majority of the 28S rRNA gene haplotypes was recovered by the surveying effort (Supplementary Figure S2). By contrast, the rarefaction curve for the COI gene haplotypes was not saturated (Supplementary Figure S3).
### 3.3. *Comparison between the Morphological and Molecular Approach*

The molecular analyses based on COI and 28S rRNA genes revealed that the sequences clearly segregated according to the morphologically pre-assigned genera (Figures 1–4, Supplementary Figures S4 and S5). *Oligosita* exhibited the lowest genetic diversity, with three highly similar sequences in the COI gene analysis and only one haplotype in the 28S rRNA gene analysis. In contrast to that, three *Oligosita* species were identified using the morphological approach.

*Paracentrobia* sequences varied only slightly, with two highly similar haplotypes in the 28S rRNA gene analysis and three highly similar haplotypes in the COI gene analysis (Figures 1–4, Supplementary Figures S4 and S5). This result indicates that all samples of *Paracentrobia* that we analyzed morphologically and molecularly belonged to a single species.

The sequences for both genes from *Gonatocerus* exhibited more variability compared to the morphological data alone. Three clades were unambiguously identified, with one cluster occurring prevalently (88.2% for the 28S gene, 83.3% for the COI gene) compared to the other two clusters (Figures 1–4, Supplementary Figures S4 and S5). *Anagrus* exhibited more variability, and was the most genetically diverse genus: 7 and 12 different haplotypes were identified in the 28S rRNA gene analysis and in the COI gene analysis, respectively (Figures 1 and 2, Supplementary Figures S4 and S5, Supplementary Table S2).

<figure><img src="image_3.png" /><figcaption>**Figure 1.** Maximum likelihood tree for the 28S rRNA gene sequences (bootstraps = 1000, TN93 model). Specimens with the same haplotype were pooled. A full list of the specimens used can be found in Supplementary Table S2. Maximum likelihood bootstrap values are given for each node.</figcaption></figure>
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Maximum likelihood tree for the COI sequences (bootstraps = 1000, GTR model). Specimens with the same haplotype were pooled. A full list of the specimens used can be found in Supplementary Table S2. Maximum likelihood bootstrap values are given for each node.</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>**Figure 3.** Overall 28S rRNA haplotype distances. Network representing the number of substitutions between the different 28S rRNA gene sequences obtained from different parasitoid genera: *Paracentrobia* (blue); *Oligosita* (green); *Anagrus* (orange); and *Gonatocerus* (red). The sample size for each haplotype is written in the respective circle. Each diamond represents a substitution. The numbers next to the dotted lines are the number of substitutions not represented in detail in the figure.</figcaption></figure>
<figure><img src="image_6.png" /><figcaption>**Figure 4.** Overall cytochrome c oxidase I (COI) haplotype distance. Network representing the number of substitutions between the different COI sequences obtained from different parasitoid genera: *Paracentrobia* (blue); *Oligosita* (green); *Anagrus* (orange); and *Gonatocerus* (red). The sample size for each haplotype is written in the respective circle. Each diamond represents a substitution. The numbers next to the dotted lines are the number of substitutions not represented in detail in the figure.</figcaption></figure>

The degrees of genetic diversity found in the genera *Anagrus* and *Gonatocerus* sequences were particularly high. Within group pairwise distances were 0.075 ± 0.012 SE (28S rRNA gene) and 0.050 ± 0.009 SE (COI gene) for *Anagrus*, and 0.069 ± 0.012 SE (28S rRNA gene) and 0.100 ± 0.016 SE (COI gene) for *Gonatocerus*. By contrast, within group pairwise distances were 0.000 ± 0.000 (28S rRNA gene) and 0.005 ± 0.003 SE (COI gene) for *Oligosita*, and 0.001 ± 0.001 SE (28S rRNA gene) and 0.009 ± 0.004 SE (COI gene) for *Paracentrobia* (Figure 5). *Paracentrobia* and *Oligosita* were significantly different from the congeneric Chalcidoidea, but not from the conspecific Chalcidoidea (Figure 5, Table 2). The genetic differences within *Gonatocerus* and *Anagrus* were in the same order of magnitude as the congeneric Chalcidoidea data.

<figure><img src="image_7.png" /><figcaption>**Figure 5.** Pairwise genetic distances between individuals from different species of the same genus in the Chalcidoidea (congeneric), different members of the same species in the genus Chalcidoidea (conspecific), and individuals collected in this study. Pairwise distances were calculated by using the K2P+G (Kimura 2-parameter model with gamma distribution) model for the 28S rRNA gene sequences and the T92+G (Tamura 3-parameter model with gamma distribution) model for the COI sequences.</figcaption></figure>
<table><caption>**Table 2.** Results of the Dunn's test (*p*) with degrees of freedom (*df*) for pairwise comparison of the pairwise genetic distances calculated within the parasitoid genera examined in this study, and the pairwise genetic differences calculated for the congeneric and conspecific Chalcidoidea sequences (n.s. stands for not significant).</caption><thead><tr><th rowspan="2">Gene</th><th rowspan="2">Parasitoid Group</th><th colspan="2">Paracentrobia spp.</th><th colspan="2">Oligosita spp.</th><th colspan="2">Anagrus spp.</th><th colspan="2">Gonatocerus spp.</th></tr><tr><th>p</th><th>df</th><th>p</th><th>df</th><th>p</th><th>df</th><th>p</th><th>df</th></tr></thead><tbody><tr><td rowspan="2">28S gene</td><td>Chalcidoidea congeneric</td><td>0.003</td><td>2</td><td>0.003</td><td>2</td><td>0.021</td><td>2</td><td>n.s.</td><td>2</td></tr><tr><td>Chalcidoidea conspecific</td><td>n.s.</td><td>2</td><td>n.s.</td><td>2</td><td>&lt;0.001</td><td>2</td><td>0.05</td><td>2</td></tr><tr><td rowspan="2">COI gene</td><td>Chalcidoidea congeneric</td><td>0.017</td><td>2</td><td>0.002</td><td>2</td><td>&lt;0.001</td><td>2</td><td>n.s.</td><td>2</td></tr><tr><td>Chalcidoidea conspecific</td><td>n.s.</td><td>2</td><td>n.s.</td><td>2</td><td>&lt;0.001</td><td>2</td><td>0.006</td><td>2</td></tr></tbody></table>

## 4. Discussion

The molecular analyses of the diversity within the four Chalcidoidea genera studied showed discrepancies with the morphological analysis. This finding is in line with previous studies on hymenopteran parasitoids [23,43,44]. For example, Mottern and Heraty [45] found that one species of the parasitoid previously described as *Cales noacki* were actually ten different *Cales* species. Although previous studies also showed that DNA barcoding can be misleading, resulting in potentially more species than actually exist [46,47], the COI gene, as the classical barcoding gene, has shown its utility in many studies [43,48,49], and constitutes the basis for the BOLD database (www.barcodinglife.org). The 28S rRNA gene is valuable to distinguish among closely related species of Chalcidoidea [50–52]. In addition, the 28S rRNA gene is generally more conserved than the COI gene for this group of insects [44,45,53]. This observation is supported by our results, as we found that sampling saturation was only reached for the 28S rRNA gene but not for the COI gene.

Once a reliable and comprehensive DNA database is established, the identification of hopper parasitoids using DNA sequencing will be a vital step towards assessing the performance and suitability of these organisms as biocontrol agents [54,55]. Our study shows that an accurate identification relies on both molecular and morphological techniques. Similar conclusions were drawn previously [21,23,56].

For instance, to date, *Anagrus* species are described as generalist parasitoids that can switch between alternative hopper species in rice agroecosystems [57,58]. Yet, the strong diversity found in the present study by using molecular tools suggests that the *Anagrus* genus is a complex of cryptic species that may be host specialists. This distinction is of high importance to control the hoppers effectively. In addition, the discovery of unknown cryptic species could expand the list of potential biological control agents [45].

Although we were not able to validate the existence of new species of *Anagrus* and *Gonatocerus* by mating tests and/or further morphometric analyses, our data strongly indicates that the two genera include unknown species. The level of within-genus genetic difference, for both 28S rRNA and COI gene sequences among samples from the *Anagrus* and *Gonatocerus* genera, exceeded the threshold of 0.17–2% within-sequence variation, which is generally accepted to delineate individuals of the same species [17,49,59]. According to the literature [17,49,59], genera from the Chalcidoidea superfamily diverge by 5.8–11.25%. Thus, the level of genetic differences we observed for *Anagrus* and *Gonatocerus* suggests that these two genera were composed of individuals from different species. When comparing the genetic analyses to the morphological results, we found strong evidence that there was at least one species of *Anagrus* and *Gonatocerus* that was not detected by the morphological analysis. Interestingly, we observed the opposite for *Oligosita* samples. In the latter genus, it is possible that one or more species were not detected due to the low success rate of DNA amplification. It is also possible that we observed a discrepancy between morphological and molecular analyses, because they were performed on a different set of specimens. As one *Oligosita* species was much more prevalent than the other two in the morphological dataset, we cannot exclude that we analyzed the most common species only, in the molecular analysis.
**5. Conclusions**

In conclusion, our results clearly demonstrate that molecular identification should be used in
combination with morphological methods and mating tests for assessing the diversity of rice hopper
parasitoids. However, further studies using an integrative approach are needed to cover the whole
diversity of parasitoids, as well as to find sustainable solutions to problems caused by BPH and GLH.
To validate the potential application of parasitoid wasps as biocontrol agents, it is of crucial importance
to have a comprehensive knowledge on their ecology and diversity. We hope that this study will
encourage further research by providing the first barcodes for egg parasitoid species from rice paddies
in Southeast Asia.

**Supplementary Materials:** The following are available online at http://www.mdpi.com/2075-4450/9/1/19/s1, Figure S1: Distribution of the eight study sites throughout the Laguna province, Figure S2: Rarefaction curve for the 28S rRNA gene, Figure S3: Rarefaction curve for the COI gene, Figure S4: Maximum likelihood tree for the 28S rRNA sequences, Figure S5: Maximum likelihood tree for the COI sequences, Table S1: Morphometric results for parasitoids throughout the Laguna province, Philippines, Table S2: Sequences obtained from parasitoids of the rice fields in the Philippines, Table S3: Sequences retrieved from the NCBI GenBank.

**Acknowledgments:** We thank Bianca Tappe, Daniel Kretzschmar and Heike Rollwage for their technical assistance. We also want to thank Jo Catinding from the International Rice Research Institute (IRRI) for the help to identify the insects and David Beck for proof reading this manuscript. The work has been funded by the LEGATO project [27], which is part of the BMBF (German Federal Ministry of Education and Research) Framework Program Research for Sustainable Development (FONA, FKZ: 01LL0917A & 01LL0917N).

**Author Contributions:** S.V. conceived and designed the experiments; C.S. performed the experiment; C.S. and F.W. did the laboratory work, C.S., A.B., F.W. and S.E. analyzed the data; R.D. contributed materials and analysis tools; C.S. wrote the first draft of the manuscript. All authors contributed to critical manuscript revision, read and approved the submitted version.

**Conflicts of Interest:** The authors declare no conflict of interest.

**References**

1. FAO. Available online: http://www.fao.org/faostat/en/#home (accessed on 30 March 2017).

2. IRRI—Trends in Global Rice Consumption. Available online: http://irri.org/rice-today/trends-in-global-rice-consumption (accessed on 1 June 2017).

3. Gurr, G.M.; Lu, Z.; Zheng, X.; Xu, H.; Zhu, P.; Chen, G.; Yao, X.; Cheng, J.; Zhu, Z.; Catindig, J.L.; et al. Multi-country evidence that crop diversification promotes ecological intensification of agriculture. *Nat. Plants* **2016**, *2*, 16014. [CrossRef] [PubMed]

4. Heong, K.L.; Aquino, G.-B.; Barrion, A.T. Population dynamics of plant- and leafhoppers and their natural enemies in rice ecosystems in the Philippines. *Crop Prot.* **1992**, *11*, 371–379. [CrossRef]

5. Ling, K.C. *Rice Virus Diseases*; International Rice Research Institute: Los Banos, Philippines, 1972.

6. Cabauatan, P.Q.; Cabunagan, R.C.; Choi, I.-R. Rice viruses transmitted by the brown planthopper *Nilaparvata lugens* Stal. In *Planthoppers: New Threats to the Sustainability of Intensive Rice Production Systems in Asia*; Heong, K.L., Hardy, B., Eds.; International Rice Research Institute: Los Banos, Philippines, 2009; pp. 357–368.

7. Catindig, J.L.; Arida, G.S.; Baehaki, S.E.; Bentur, J.S.; Cuong, L.Q.; Norowi, M.; Rattanakarn, W.; Sriratanasak, W.; Xia, J.; Lu, Z. Situations of planthoppers in Asia. In *Planthoppers: New Threats to the Sustainability of Intensive Rice Production Systems in Asia*; Heong, K.L., Hardy, B., Eds.; International Rice Research Institute: Los Banos, Philippines, 2009; pp. 191–220.

8. Seo, B.Y.; Jung, J.K.; Choi, B.-R.; Park, H.M.; Lee, B.H. Resistance-breaking ability and feeding behavior of the brown planthopper *Nilaparvata lugens*, recently collected in Korea. In *Planthoppers: New Threats to the Sustainability of Intensive Rice Production Systems in Asia*; International Rice Research Institute: Los Banos, Philippines, 2009; pp. 303–314.

9. Heong, K.L.; Schoenly, K.G. Impact of insecticides on herbivore-natural enemy communities in tropical rice ecosystems. In *Ecotoxicology*; Springer: Boston, MA, USA, 1998; pp. 381–403, ISBN 978-1-4613-7653-8.
10. Heong, K.L. Are planthopper problems caused by a breakdown in ecosystem services? In *Planthoppers: New Threats to the Sustainability of Intensive Rice Production Systems in Asia*; Heong, K.L., Hardy, B., Eds.; International Rice Research Institute: Los Banos, Philippines, 2009; pp. 221–232, ISBN 978-971-22-0251-3.

11. Gurr, G.M.; Liu, J.; Read, D.M.Y.; Catindig, J.L.A.; Cheng, J.A.; Lan, L.P.; Heong, K.L. Parasitoids of Asian rice planthopper (Hemiptera: Delphacidae) pests and prospects for enhancing biological control by ecological engineering. *Ann. Appl. Biol.* **2011**, *158*, 149–176. [CrossRef]

12. Antolin, M.F.; Strong, D.R. Long-distance dispersal by a parasitoid (*Anagrus delicatus*, Mymaridae) and its host. *Oecologia* **1987**, *73*, 288–292. [CrossRef] [PubMed]

13. Fowler, S.V.; Claridge, M.F.; Morgan, J.C.; Peries, I.D.R.; Nugaliyadde, L. Egg mortality of the brown planthopper, *Nilaparvata lugens* (Homoptera: Delphacidae) and green leafhoppers, *Nephotettix* spp. (Homoptera: Cicadellidae), on rice in Sri Lanka. *Bull. Entomol. Res.* **1991**, *81*, 161–167. [CrossRef]

14. Watanabe, T.; Wada, T.; Salleh, N.M.N. Parasitic Activities of Egg Parasitoids on the rice planthoppers, *Nilaparvata lugens* (STAL) and *Sogatella furcifera* (HORVATH) (Homoptera: Delphacidae), in the Muda Area, Peninsular Malaysia. *Appl. Entomol. Zool.* **1992**, *27*, 205–211. [CrossRef]

15. Drechsler, M.; Settele, J. Predator–prey interactions in rice ecosystems: Effects of guild composition, trophic relationships, and land use changes—A model study exemplified for Philippine rice terraces. *Ecol. Model.* **2001**, *137*, 135–159. [CrossRef]

16. Nishida, T.; Wongsiri, T.; Wongsiri, N. Species composition, population trends and egg parasitism of planthopper and leafhopper rice pests of Thailand. *Plant Prot. Bull. FAO* **1976**, *24*, 22–26.

17. Chesters, D.; Wang, Y.; Yu, F.; Bai, M.; Zhang, T.-X.; Hu, H.-Y.; Zhu, C.-D.; Li, C.-D.; Zhang, Y.-Z. The integrative taxonomic approach reveals host specific species in an encyrtid parasitoid species Complex. *PLoS ONE* **2012**, *7*, e37655. [CrossRef] [PubMed]

18. Hebert, P.D.N.; Penton, E.H.; Burns, J.M.; Janzen, D.H.; 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. [CrossRef] [PubMed]

19. Smith, M.A.; Rodriguez, J.J.; Whitfield, J.B.; Deans, A.R.; Janzen, D.H.; Hallwachs, W.; Hebert, P.D.N. Extreme diversity of tropical parasitoid wasps exposed by iterative integration of natural history, DNA barcoding, morphology, and collections. *Proc. Natl. Acad. Sci. USA* **2008**, *105*, 12359–12364. [CrossRef] [PubMed]

20. Schwarzfeld, M.D.; Broad, G.R.; Sperling, F.A.H. Molecular phylogeny of the diverse parasitoid wasp genus *Ophion Fabricius* (Hymenoptera: Ichneumonidae: Ophioninae): Phylogeny of *Ophion. Syst. Entomol.* **2016**, *41*, 191–206. [CrossRef]

21. Munro, J.B.; Heraty, J.M.; Burks, R.A.; Hawks, D.; Mottern, J.; Cruaud, A.; Rasplus, J.-Y.; Jansta, P. A molecular phylogeny of the Chalcidoidea (Hymenoptera). *PLoS ONE* **2011**, *6*, e27023. [CrossRef] [PubMed]

22. Smith, M.A.; Woodley, N.E.; Janzen, D.H.; Hallwachs, W.; Hebert, P.D.N. DNA barcodes reveal cryptic host-specificity within the presumed polyphagous members of a genus of parasitoid flies (Diptera: Tachinidae). *Proc. Natl. Acad. Sci. USA* **2006**, *103*, 3657–3662. [CrossRef] [PubMed]

23. Zhou, Q.-S.; Xi, Y.-Q.; Yu, F.; Zhang, X.; Li, X.-J.; Liu, C.-L.; Niu, Z.-Q.; Zhu, C.-D.; Qiao, G.-X.; Zhang, Y.-Z. Application of DNA barcoding to the identification of Hymenoptera parasitoids from the soybean aphid (*Aphis glycines*) in China: DNA barcoding to identify parasitoids from *Aphis glycines* in China. *Insect Sci.* **2014**, *21*, 363–373. [CrossRef] [PubMed]

24. Smith, M.A.; Wood, D.M.; Janzen, D.H.; Hallwachs, W.; Hebert, P.D.N. DNA barcodes affirm that 16 species of apparently generalist tropical parasitoid flies (Diptera, Tachinidae) are not all generalists. *Proc. Natl. Acad. Sci. USA* **2007**, *104*, 4967–4972. [CrossRef] [PubMed]

25. Padial, J.M.; Miralles, A.; De la Riva, I.; Vences, M. The integrative future of taxonomy. *Front. Zool.* **2010**, *7*, 16. [CrossRef] [PubMed]

26. Sann, C.; Theodorou, P.; Heong, K.L.; Villareal, S.; Settele, J.; Vidal, S.; Westphal, C. Hopper parasitoids do not significantly benefit from non-crop habitats in rice production landscapes. *Agric. Ecosyst. Environ.* **2018**, *254*, 224–232. [CrossRef]

27. Settele, J.; Spangenberg, J.H.; Heong, K.L.; Burkhard, B.; Bustamante, J.V.; Cabbigat, J.; Van Chien, H.; Escalada, M.; Grescho, V.; Hai, L.H.; et al. Agricultural landscapes and ecosystem services in South-East Asia—The LEGATO-Project. *Basic Appl. Ecol.* **2015**, *16*, 661–664. [CrossRef]

28. Heinrichs, E.A.; Medrano, F.G.; Rapusas, H.R. *Genetic Evaluation for Insect in Rice*; International Rice Research Institute: Los Banos, Philippines, 1985; ISBN 971-10-4110-3.
29. Shepard, B.M.; Barrion, A.T.; Listinger, J.A. *Friends of the Rice Farmer. Helpful Insects, Spiders and Pathogens*; International Rice Research Institute: Los Banos, Philippines, 1987; ISBN 971-10-4162-6.

30. Dupo, A.L.; Barrion, A.T. Taxonomy and general biology of delphacid planthoppers in rice agroecosystems. In *Planthoppers: New Threats to the Sustainability of Intensive Rice Production Systems in Asia*; Heong, K.L., Hardy, B., Eds.; International Rice Research Institute: Los Banos, Philippines, 2009; pp. 3–156, ISBN 978-971-22-0251-3.

31. Heinrichs, E.A. *Biology and Management of Rice Insects*; Wiley Eastern Limited: New Delhi, India, 1994; ISBN 81-224-0581-2.

32. 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. [PubMed]

33. Campbell, B.C.; Steffen-Campbell, J.D.; Werren, J.H. Phylogeny of the Nasonia species complex (Hymenoptera: Pteromalidae) inferred from an internal transcribed spacer (ITS2) and 28S rDNA sequences. *Insect Mol. Biol.* **1993**, *2*, 225–237. [CrossRef] [PubMed]

34. Campbell, B.C.; Heraty, J.M.; Rasplus, J.Y.; Chan, K.; Steffen-Campbell, J.D.; Babcock, C. Molecular systematic of the Chalcidoidea using 28S-D2 rDNA. In *Hymenoptera: Evolution, Biodiversity and Biological Control*; Heong, K.L., Hardy, B., Eds.; CSIRO Publishing: Collingwood, Australia, 2000; pp. 59–73, ISBN 0-643-06610-1.

35. Bonfield, J.K.; Whitwham, A. Gap5—Editing the billion fragment sequence assembly. *Bioinformatics* **2010**, *26*, 1699–1703. [CrossRef] [PubMed]

36. Tamura, K.; Dudley, J.; Nei, M.; Kumar, S. MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. *Mol. Biol. Evol.* **2007**, *24*, 1596–1599. [CrossRef] [PubMed]

37. Altschul, S.F.; Gish, W.; Miller, W.; Myers, E.W.; Lipman, D.J. Basic local alignment search tool. *J. Mol. Biol.* **1990**, *215*, 403–410. [CrossRef]

38. Bandelt, H.J.; Forster, P.; Rohl, A. Median-joining networks for inferring intraspecific phylogenies. *Mol. Biol. Evol.* **1999**, *16*, 37–48. [CrossRef] [PubMed]

39. Dunn, O.J. Multiple comparisons using rank sums. *Technometrics* **1964**, *6*, 241–252. [CrossRef]

40. R Development Core Team. *The R Foundation for Statistical Computing, Ver. 3.0.2*; Vienna University of Technology: Vienna, Austria, 2013. Available online: http://www.r-project.org/ (accessed on 3 May 2017).

41. Dixon, P. VEGAN, a package of R functions for community ecology. *J. Veg. Sci.* **2003**, *14*, 927–930. [CrossRef]

42. Ritz, C.; Baty, F.; Streibig, J.C.; Gerhard, D. Dose-response analysis using R. *PLoS ONE* **2015**, *10*, e0146021. [CrossRef] [PubMed]

43. Smith, M.; Fernández-Triana, J.L.; Eveleigh, E.; Gómez, J.; Guclu, C.; Hallwachs, W.; Hebert, P.D.N.; Hrcek, J.; Huber, J.T.; Janzen, D.; et al. DNA barcoding and the taxonomy of Microgastrinae wasps (Hymenoptera, Braconidae): Impacts after 8 years and nearly 20000 sequences. *Mol. Ecol. Resour.* **2013**, *13*, 168–176. [CrossRef] [PubMed]

44. Kenyon, S.G.; Buerki, S.; Hansson, C.; Alvarez, N.; Benrey, B. Uncovering cryptic parasitoid diversity in Horismenus (Chalcidoidea, Eulophidae). *PLoS ONE* **2015**, *10*, e0136063. [CrossRef] [PubMed]

45. Mottern, J.L.; Heraty, J.M. Revision of the *Cales noacki* species complex (Hymenoptera, Chalcidoidea, Aphelinidae). *Syst. Entomol.* **2014**, *39*, 354–379. [CrossRef]

46. Brower, A.V.Z. Problems with DNA barcodes for species delimitation: "Ten species" of *Astraptes fulgerator* reassessed (Lepidoptera: Hesperiidae). *Syst. Biodivers.* **2006**, *4*, 127–132. [CrossRef]

47. Collins, R.A.; Cruickshank, R.H. The seven deadly sins of DNA barcoding. *Mol. Ecol. Resour.* **2013**, *13*, 969–975. [CrossRef] [PubMed]

48. Gutiérrez-Arellano, D.; Gutiérrez-Arellano, C.; Zaldívar-Riverón, A. DNA Barcoding of the parasitoid wasp subfamily Doryctinae (Hymenoptera: Braconidae) from Chamela, Mexico. *Biodivers. Data J.* **2015**, *3*, e5109. [CrossRef]

49. Hebert, P.D.N.; Cywinska, A.; Ball, S.L.; deWaard, J.R. Biological identifications through DNA barcodes. *Proc. R. Soc. B Biol. Sci.* **2003**, *270*, 313–321. [CrossRef] [PubMed]

50. Babcock, C.S.; Heraty, J.M.; De Barro, P.J.; Driver, F.; Schmidt, S. Preliminary phylogeny of *Encarsia Förster* (Hymenoptera: Aphelinidae) based on morphology and 28S rDNA. *Mol. Phylogenet. Evol.* **2001**, *18*, 306–323. [CrossRef] [PubMed]
51. Manzari, S.; Polaszek, A.; Belshaw, R.; Quicke, D.L.J. Morphometric and molecular analysis of the *Encarsia inaron* species-group (Hymenoptera: Aphelinidae), parasitoids of whiteflies (Hemiptera: Aleyrodidae). *Bull. Entomol. Res.* **2002**, *92*, 165–176. [CrossRef] [PubMed]

52. Triapitsyn, S.V.; Vickerman, D.B.; Heraty, J.M.; Logarzo, G.A. A new species of *Gonatocerus* (Hymenoptera: Mymaridae) parasitic On *Proconiine Sharpshooters* (Hemiptera: Cicadellidae) in the New World. *Zootaxa* **2006**, *1158*, 55–67.

53. Gillespie, J.J. A secondary structural model of the 28S rRNA expansion segments D2 and D3 for chalcidoid wasps (Hymenoptera: Chalcidoidea). *Mol. Biol. Evol.* **2005**, *22*, 1593–1608. [CrossRef] [PubMed]

54. Greenstone, M.H. Molecular methods for assessing insect parasitism. *Bull. Entomol. Res.* **2006**, *96*, 1–13. [CrossRef] [PubMed]

55. Gariepy, T.D.; Kuhlmann, U.; Gillott, C.; Erlandson, M. Parasitoids, predators and PCR: The use of diagnostic molecular markers in biological control of Arthropods. *J. Appl. Entomol.* **2007**, *131*, 225–240. [CrossRef]

56. Desneux, N.; Starý, P.; Delebecque, C.J.; Gariepy, T.D.; Barta, R.J.; Hoelmer, K.A.; Heimpel, G.E. Cryptic species of parasitoids attacking the soybean aphid (Hemiptera: Aphididae) in Asia: Binodoxys communis and binodoxys koreanus (Hymenoptera: Braconidae: Aphidiinae). *Ann. Entomol. Soc. Am.* **2009**, *102*, 925–936. [CrossRef]

57. Zheng, X.; Yu, X.; Lu, Z.; Chen, J.; Xu, H.; Ju, R. Parasitization adaptability of *Anagrus optabilis* on *Nilaparvata lugens*. *Chin. J. Biol. Control* **2003**, *19*, 136–138.

58. Zheng, X.; Lu, Y.; Zhu, P.; Zhang, F.; Tian, J.; Xu, H.; Chen, G.; Nansen, C.; Lu, Z. Use of banker plant system for sustainable management of the most important insect pest in rice fields in China. *Sci. Rep.* **2017**, *7*, 45581. [CrossRef] [PubMed]

59. Meier, R.; Zhang, G.; Ali, F. The use of mean instead of smallest interspecific distances exaggerates the size of the “Barcoding Gap” and leads to misidentification. *Syst. Biol.* **2008**, *57*, 809–813. [CrossRef] [PubMed]

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

© 2018 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).