Next Article in Journal
Biological Adaptations Associated with Dehydration in Mosquitoes
Next Article in Special Issue
Life on the Edge: Ecological Genetics of a High Arctic Insect Species and Its Circumpolar Counterpart
Previous Article in Journal
Insecticidal and Antifeedant Activities of Malagasy Medicinal Plant (Cinnamosma sp.) Extracts and Drimane-Type Sesquiterpenes against Aedes aegypti Mosquitoes
Previous Article in Special Issue
Comparative Analysis of Intra- and Inter-Specific Genomic Variability in the Peach Potato Aphid, Myzus persicae
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Insects
Volume 10
Issue 11
10.3390/insects10110374
insects-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Unraveling the Genetic Structure of the Coconut Scale Insect Pest (Aspidiotus rigidus Reyne) Outbreak Populations in the Philippines

by
Joeselle M. Serrana
1,2,
Naoto Ishitani
1,2,
Thaddeus M. Carvajal
1,2,
Billy Joel M. Almarinez
2,3,
Alberto T. Barrion
2,3,
Divina M. Amalin
2,3 and
Kozo Watanabe
1,2,*
1
Department of Civil and Environmental Engineering, Ehime University, Bunkyo-cho 3, Matsuyama 790-8577, Japan
2
Biological Control Research Unit, Center for Natural Sciences and Environmental Research, De La Salle University, 2401 Taft Avenue, Manila 1004, Philippines
3
Biology Department, College of Science, De La Salle University, 2401 Taft Avenue, Manila 1004, Philippines
*
Author to whom correspondence should be addressed.
Insects 2019, 10(11), 374; https://0-doi-org.brum.beds.ac.uk/10.3390/insects10110374
Submission received: 12 September 2019 / Revised: 14 October 2019 / Accepted: 24 October 2019 / Published: 26 October 2019
(This article belongs to the Special Issue Cutting-Edge Evolution of Insects: Factors Governing Ecological Divergence of Natural Populations Assessed Using Molecular and Genetic Markers)
Browse Figures
Versions Notes

Abstract

:
Despite the fact that massive outbreaks of the coconut scale insect pest, Aspidiotus rigidus Reyne (Hemiptera: Diaspididae) are inflicting significant economic losses to the Philippines’ coconut industry, little is known about the population and dispersal history of this invasive pest in the country. Here, we examined the genetic diversity, structure and demographic history of A. rigidus sampled from localities with reported outbreaks from 2014 to 2017. We analyzed the genetic structure of the outbreak populations using mitochondrial COI and nuclear EF-1α markers. Both markers and all methods of population genetic structure analyses indicate clear differentiation among the A. rigidus populations separating the north from the southern regions of the Philippines. Very low or no genetic differentiation was observed within and amongst the populations per geographic region indicating two unrelated outbreak events of the pest originating from two genetically uniform populations isolated in each respective region. Historical data supports the resurgence of an established A. rigidus population in the south which could have been driven by sudden climatic changes or human-induced habitat imbalance. Our study provides valuable information on the genetic differentiation of the two A. rigidus groups that would be useful for developing and implementing biological control strategies against this pest in the Philippines.

1. Introduction

Insect pest outbreaks are characterized by an explosive increase in an insect population occurring over a relatively short period [1]. Large and rapid alterations in the environment or changes in the intrinsic genetic or physiological properties of individual organisms within a population can result to the resurgence of insect pests to outbreak-level status [2,3]. Likewise, insect outbreaks may occur when non-native species have no natural enemies or a few inefficient natural enemies, and if the local beneficial species are unable to suppress them in the area of introduction [4,5]. The invasive success of pest species may be determined by both the biology and environmental factors promoting its spread in a suitable area [6,7]. A better understanding of the source population, route and the mechanism of spread could provide valuable insights for designing and implementing quarantine strategies to understand the invasion success and decline of outbreak populations [4,8].
In 2009, the Philippines suffered a devastating coconut scale insect outbreak damaging the coconut palms on the provinces of Luzon (northern region of the Philippines) and currently on some areas in Mindanao (southern region) inflicting significant economic losses to the country’s coconut industry. The diaspidid scale insect Aspidiotus rigidus Reyne (Hemiptera: Diaspididae) resides on the underside of the leaf, blocks the stomata and sucks plant sap strongly reducing the plant’s photosynthetic activity leading to a characteristic yellowing and drying of the leaves. Severely infested coconut palms dry up and die within six months or less [9]. Prior to the renewed interest on A. rigidus due to the outbreak in the Philippines, historical and observational data on the spread of the invasive coconut pest has been scarce. Other than the most recent observation published by Watson et al. [10] the last known study on the biology of this invasive species was conducted by Reyne [9,11] with full documentation of the outbreak on the island of Sangi (North Celebes) in Indonesia from mid-1925 to 1928. The outbreak recorded by Reyne [9] came to an end naturally after two years due to reduced female fecundity and high mortality of the immature stages. The decrease in the A. rigidus population may have been associated with natural enemies that, over time, regulated the pest population. However, the outbreaks from the localities infested with A. rigidus in the Philippines took longer to recover [10], e.g., six years for the northern province of Batangas, and are still on-going for the southern areas, i.e., Basilan, Zamboanga Peninsula, and the Caraga Region.
The introduction of A. rigidus to the Philippines, and its spread was believed to be either by wind or accidental transportation of infested plants, coconut planting materials and products [10]. Infestation in the northern provinces of the Philippines spread rapidly from its initial local report in Tanauan, Batangas in the Calabarzon Region (Luzon) in 2009, reaching nearby coconut planted areas throughout the region. These outbreaks lasted for at least three years [10] and were reported to be manageable by 2015 [12]. The more recent outbreak in the southern region, specifically in Basilan, started in early 2013 [10], implying its direct connection with the northern outbreak. However, given the means of spread via wind wherein crawlers are dispersed from one area to another [10], it is highly improbable that the infestation from the northern region reached the infested southern islands meanwhile moving past other provinces planted with coconut palms en route. In addition, transport of infested plants from the northern provinces was highly unlikely given the national attention focused on quarantine and management strategies against the spread of the coconut scale insect during the outbreak [12,13].
It is also likely that A. rigidus has been in the country as a minor pest and regulated by natural enemies. Based on historical reports, Lever [14] reported sightings of A. rigidus in the Philippines as long ago as 1969, and Velasquez [15] recounted that the pest was probably highly confined in the southern part of the country. It was more likely that the source of the observed A. rigidus came from the island of Sangi in Indonesia given its relative closeness to Mindanao. In time, the immigrant A. rigidus could have established a resident population complemented by natural enemies limiting its colonization outside the area of introduction. Changes in anthropogenic, biotic interactions or climatic factors can influence a population’s rise to an outbreak level [3,16]. The recent outbreak observed in the southern part of the Philippines may have been caused by a sudden rise in the abundance of the supposedly established A. rigidus population due to factors such as human-induced habitat imbalance e.g., excessive use of pesticide affecting the natural enemies controlling the pest population, or climate change such as a prolonged dry spell, which may induce changes in the local biotic community.
Inference of the source population, route and the mechanism of spread of A. rigidus in the Philippines needs further assessment and confirmation. Tracing the history of an invasion or identifying the geographic origin of a pest population can be done by characterizing population-level genetic variation using molecular markers e.g., [17,18,19,20]. Sequencing selected gene fragments, e.g., mitochondrial COI is a traditional population genetic tool providing insights on dispersal pathways and population structure. The mitochondrial cytochrome oxidase (mtCOI) gene and the nuclear protein-encoding gene elongation factor 1α (EF-1α) have been commonly used in studies investigating the origin [21,22], or inference of phylogenetic relationships [23,24] of various invasive diaspidid species.
Here, we aimed to assess the population genetic structure and demography of the outbreak populations of the coconut scale insect, A. rigidus in the Philippines. Given the historical documentation of the pest in the southern region and the relatively extensive and rapid spread, but faster recovery, of the infestations in the northern region compared to the southern outbreaks i.e., Basilan and Zamboanga Peninsula, we hypothesize the presence of two distinct genetic groups for the outbreak events isolated within each geographic region. A population genetics approach is a useful tool to examine whether the northern and southern coconut scale insect outbreaks originated from immigrant or resident populations. To test this hypothesis, we utilized sequences of the mtCOI gene and the nuclear protein-encoding gene - EF-1α to investigate the genetic structure and diversity of A. rigidus populations from localities with documented outbreak-level infestations in the Philippines from 2014 to 2017. Furthermore, we employed a coalescent genealogy approach to provide additional evidence on the demographic relationship of the outbreak A. rigidus populations between the northern and southern geographic regions in the Philippines.

2. Materials and Methods

2.1. Sample Collection

Aspidiotus rigidus populations were sampled at seven localities with reported coconut scale insect outbreak across the Philippines from 2014 to 2017 (Figure 1). The northern localities sampled were Orani, Bataan (BT; N14.769786, E120.454510), Nagcarlan (NG; N14.158930, E121.413670) and San Pablo (SP; N14.056420, E121.333300), Laguna, Tanauan (TN; N14.098870, E121.091330) and Talisay (TL; N14.093340, E121.010730) Batangas. The southern localities are Basilan (BS; N06.707853, E121.983358) and Zamboanga (ZB; N06.993166, E121.927963). Detailed information regarding location and sample collection is collected in Table 1.
Mature female scale insects were identified as A. rigidus based on the characteristic distribution of egg skins, which for this species occurs along the posterior or pygidial half of the insect body (Figure 2). Non-parasitized adult females were carefully selected from infested leaves and preserved in 95% ethanol before molecular analysis. To further confirm identification as well as the purity of samples, A. rigidus collected from Orani, Bataan were reared at the rearing facility of the Biological Control Research Unit (BCRU) located at De La Salle University (DLSU), Science and Technology Complex (Binan City, Laguna) on mangosteen (Garcinia mangostana L.), a differential host of A. rigidus observed not to support multiple generations of A. destructor. A phylogenetic analysis was employed (see the molecular analysis section below) to confirm the identifications of the field-collected samples by comparing them to A. destructor and mangosteen-reared A. rigidus sequences.

2.2. DNA Extraction, PCR Amplification, and Sequencing

Genomic DNA was extracted individually using the DNeasy Blood & Tissue Kit (Qiagen, Hilden, Germany) following the manufacturer’s guidelines. Extraction was performed by crushing the insect body of each individual in separate microcentrifuge tubes using a micropestle. DNA concentration and quality were assessed by spectrophotometry (NanoDrop 2000 spectrophotometer, ThermoScientific, Wilmington, DE, USA). The mitochondrial COI gene was amplified using the forward primer PcoF1 designed for scale insects by Park et al. [25] and the standard reverse primer LepR1. The nuclear gene EF-1α was amplified using the forward primer EF-1α of Morse and Normark [26] paired with the EF2 reverse primer [27]. PCR reactions were performed in a 25 μL reaction mixture containing 10× Buffer, 2.5 mM dNTP mixture, 25 mM MgCl2, 10 pmol of each primer, 1 U of Taq DNA polymerase (TaKaRa Bio Inc., Japan), and 2 to 50 ng of template DNA. PCR thermocycling was performed in a T100™ Thermal Cycler (Bio-Rad, CA, USA). Following the conditions from Park et al. [28], the mtCOI gene was amplified with an initial denaturation step at 95 °C for 5 min, followed by 5 cycles of 94 °C for 40 s, annealing at 45 °C for 40 s, extension at 72 °C for 1 min and 10 s, and another 35 cycles of denaturation at 94 °C for 40 s, annealing at 51 °C for 40 s, extension at 72 °C for 1 min and 10 s, and a 5 min final extension at 72 °C after the last cycle. Furthermore, after an initial denaturation at 95 °C, with a denaturation at 95 °C for 30 s and extension at 72 °C for 2 min every cycle, a touch-down procedure was performed for the amplification of the EF-1α gene following the protocol of Morse and Normark [26] in which the initial annealing temperature of 58 °C was decreased by 2 °C every three cycles until a final temperature of 42 °C was reached, then held for 18 cycles followed by a 5 min final extension at 72 °C. PCR products were visualized in 1.5% agarose gels stained with Midori Green Direct (Nippon Genetics Co. Ltd., city, Japan), and cleaned using the QIAquick PCR Purification Kit (Qiagen, Tokyo, Japan). Samples were sent to Eurofins Genomics (Eurofins Genomics Co., Ltd., city, Japan) for Sanger sequencing to produce both forward and reverse fragments.

2.3. Genetic Diversity and Population Structure

Sequences were assembled using CodonCode Aligner v. 5.1.5 (CodonCode Corporation, Dedham, MA, USA). Concatenated sequences were deposited in GenBank with the following accession numbers: MK862684-MK862988 for mtCOI and MN603810-MN603884 for EF-1α before the subsequent molecular analysis, the sequences were aligned via MAFFT v. 7.409 [29], and the ambiguously aligned regions were excluded using GBlocks 0.91b [30]. Sequence polymorphisms for both mtCOI and EF-1α gene were assessed. The number of variable sites (S) and haplotypes (h), average number of nucleotide difference (k), haplotype diversity (Hd), and nucleotide diversity (Pi) of the two marker genes were calculated in DnaSP v. 6.10.04 [31]. The hierarchical analysis of molecular variance (AMOVA) was implemented in Arlequin v. 3.5.2.2 [32]. Two geographic groups were defined: The northern (Luzon) and the southern (Mindanao) groups prior to the analysis. Population pairwise FST was computed in Arlequin v. 3.5.2.2 using 1000 permutations.

2.4. Demographic Inference

Tajima’s D and Fu’s Fs statistic tests were estimated to infer demographic history and dynamics in each population, for the two geographic groups, and for all populations grouped together in Arlequin v. 3.5.2.2 for both mtCOI and EF-1α datasets. Also, Fu and Li’s D* and F* test statistics were computed in DnaSP v. 5.0 to determine departures from the mutation-drift equilibrium [33]. Parameters of demographic expansion included moment estimators of time to the expansion Tau, effective population size before expansion (Theta0, θ0), and effective population size after expansion (Theta1, θ1) between the observed and expected mismatches. The adjustment to a model of population expansion was estimated from the sum of squared deviation (SSD) and the raggedness index (r) in Arlequin v. 3.5.2.2.

2.5. Gene Flow Analysis and Median-joining Networks of Haplotypes

To test migration history between the two geographic groups, we calculated Bayes factors from the marginal likelihoods estimated in MIGRATE v. 4.4.0 [34,35] based on both mtCOI and EF-1α datasets. MIGRATE-N utilizes marginal likelihoods to compare and order structured population models [35]. The program provides estimates of historic gene-flow with the assumption that populations have reached mutation-migration-drift equilibrium. We tested eight possible models of migration history. Model 1 allows migration between the two groups, with the populations assumed to exist for a very long time. Model 2 presents a migration from northern to southern group, while model 3 presents the opposite situation and with the populations assumed to exist for a very long time. Model 4 assumed one panmictic population encompassing the northern and southern groups. Model 5 allows divergence among populations within southern group splitting from the northern group, and migration from north to south, with the northern group existing for a long time and the southern group recently splitting off. Model 6 is a mirror image of model 5. Model 7 is similar to model 5 except that no interaction occurred between the two groups after the split. Model 8 represented the opposite situation of model 7. Similar parameters were used to run all models.
A Bayesian search strategy was performed with the following parameters: one long chain (10,000 trees) with a burn-in of 5000 iterations. A static heating scheme with 4 chains was applied using temperature parameters set by default with a swapping interval of one. Bayes factors were calculated via “BF” implemented in carlopacioni/mtraceR, a package for analyzing migrate-n outputs in R v. 3.5.2. Log Bayes factors (LBF) of all models were calculated by comparing against the model that has the highest log-likelihood. The models are ranked based on LBF and calculated model probability.
Median-joining (MJ) networks [36] for the two markers were constructed to estimate the genealogical relationship in A. rigidus haplotypes via PopART v. 1.7 [37].

2.6. Phylogenetic Analysis

Since no site variation was observed between the sequences in each population (except for the EF-1α sequences from Basilan), six representative samples per population, a total of 48 mtCOI sequences along with 42 plus all 16 Basilan EF-1α sequences were chosen for the phylogenetic analysis. Sequences of A. destructor collected from coconut palms, identified based on the circular distribution of egg skins were selected as an outgroup (accession number: MK862989 for mtCOI and MN603885 for EF-1α). The Akaike information criterion corrected for sample size (AICc) was implemented to find the best fitting evolutionary model for phylogenetic reconstruction via jModelTest v. 2.1.10 [38]. The evolutionary model for the mtCOI sequences was TIM2+G, while TrNef was the model for the EF-1α sequences. Maximum likelihood (ML) tree inference was performed in RAxML-NG v. 0.5.1 [39] with 1000 bootstrap replicates.

3. Results

3.1. Genetic Diversity and Population Structure

All samples identified based on the characteristic distribution of egg skins were confirmed as A. rigidus by comparing the sequences with the mangosteen-reared A. rigidus, and A. destructor sequences. DNA sequence analysis of all the concatenated 647-bp mtCOI sequences of 305 individuals from seven A. rigidus outbreak populations collected from 2014 to 2017 in the Philippines, with the mangosteen-reared samples revealed only two distinct haplotypes (h), separated by 31 polymorphic sites (s) (Figure 3). Haplotype diversity (Hd) was calculated to be 0.50 +/− 0.005 SD. Average number of nucleotide difference (k) was 15.5 and nucleotide diversity (Pi) was 0.024 +/− 0.00025 SD. No variation was found in the sequences of samples collected per populations. Reared A. rigidus and samples collected from the five populations of northern group, i.e., Orani, Bataan (BT), Nagcarlan (NG) and San Pablo (SP), Laguna, Tanauan (TN) and Talisay (TL), Batangas are grouped into one haplotype. Samples from the two populations of southern group, i.e., Basilan (BS) and Zamboanga (ZB) were grouped together in the second mtCOI haplotype (Table 2). For the nuclear EF-1α gene, 75 concatenated sequences of 1007-bp length generated 14 polymorphic sites (s), with four haplotypes (h). Like the mtCOI sequences, all samples from the northern group clustered into one haplotype. For the southern group, samples from ZB grouped into one haplotype, while BS were separated into two haplotypes. Haplotype diversity (Hd) was calculated to be 0.48 +/− 0.064 SD. Average number of nucleotide difference (k) was 5.7 and nucleotide diversity (Pi) was 0.0057 +/− 0.00061 SD. Except for BS, all other localities have no sequence variation per locality.
In addition, genetic diversity parameters were calculated per geographic group. The mtCOI sequences for both groups, and the EF-1α sequence of the northern group showed no sequence variation. The EF-1α sequences from the southern group had three haplotypes with an estimated Hd of 0.688 +/− 0.039 SD, with a k value of 0.087 and a low Pi value of 0.00086 +/− 0.00010 (Table 2). Both the median-joining haplotype network and ML inferred trees present distinct two and four haplotypes for the mtCOI and EF-1α dataset, respectively (Figure 3 and Figure S1).
AMOVA analysis indicated a highly structured genetic variability of 100% and 97.63% variations among the groups for mtCOI and EF-1α dataset, respectively. There were zero, or a relatively small percentage of variation among populations within groups and within populations. Except for the source of variation among populations within groups in the mtCOI data, AMOVA showed that significant genetic structure occurred in A. rigidus at various hierarchical levels (Table 3). Pairwise FST values varied from zero to 1.00 for mtCOI, and zero to 0.98 for the EF-1α dataset. The differentiation between populations was only significant when the comparison was between a northern and a southern population (Table S1). Moreover, pairwise FST values between the two groups revealed a high value of 1.00 and 0.97 for the mtCOI and EF-1α dataset, respectively. As shown in the ML trees (Figure S1), the phylogenetic analyses of both markers were consistent with the results of the analyses above. Samples clustered according to their geographic group, with the BS and ZB EF-1α sequences in three separate nodes.

3.2. Demographic History and Gene Flow

For both datasets, neutrality test computation for all samples showed positive values for Tajima’s D, Fu and Li’s D*, Fu’s Fs, and Fu and Li’s F*, and were significant for the first two parameters (Table 2). Estimations per population for these parameters were mostly zero or positive but not significant, suggesting neither population expansion nor purifying selection in these populations. Estimations of the SSD and r parameters both returned zero values, except for the EF-1α sequences from Basilan, with a significant SSD of 0.0283 (p < 0.001) and a non-significant r of 0.2871. Other demographic parameters such as Tau, θ0 and θ1 index, are presented in Table S2.
Results of the analysis in migrate-n were presented in Table S3. We found contrasting results for the two markers employed using Bayes factors to compare the eight models of dispersal. For the mtCOI dataset, model 7 was ranked best with a probability of 0.996. For the EF-1α dataset, model 3 was ranked best with a probability of 1.000.

4. Discussion

The main aim of the present study was to elucidate the genetic structure and demography of the coconut scale insect pest A. rigidus from selected localities in the Philippines with reported heavy infestations collected from 2014 to 2017. Both the mtCOI and EF-1α markers and all methods of population structure analyses revealed strong differentiation among the A. rigidus populations separating the northern (Luzon) outbreak from the southern (Mindanao) region. The separation of the populations by geographic region and the observed lack of genetic variability within populations were represented graphically in the median-joining network and phylogenetic analysis employed in the study (Table 1 and Figure 3).

4.1. Genetic Structure of A. rigidus: Evidence of Coconut Scale Insect “Superclones” in the Philippines

Our results indicate the existence of two mitochondrial, and four nuclear haplotypes (one northern and three southern). Genetic population clusters result from multiple source populations contributing to an insect pest outbreak [8]. However, we only observed two clusters separating the outbreak populations into their respective geographic regions. Also, genetic variation was either absent or very low within and amongst the populations of the northern and the southern region, implying that populations from each region consisted of a single genotype. Hence, the presence of two distinct A. rigidus single genotype populations or “superclones” [40] in the Philippines which supports our hypothesis on the occurrence of two genetically unrelated outbreak events in the country.
Armoured diaspid scale insects were observed to have unusual genetic systems such as diverse types of paternal genome elimination and parthenogenesis [22]. Specifically, several aspidiotine insects have obligate [24,41] or cryptic parthenogenetic populations as observed in Aspidiotus nerii Bouche [21,22]. Accordingly, A. rigidus was observed to reproduce parthenogenetically [9,10]. Yellow winged adult males are seen in outbreak populations but the sex ratio varies widely with males thought to be non-functional [9]. Parthenogenetic reproduction has been thought to be the leading driver to the dominance of “superclones” across space and time in other insects such as aphids (Hemiptera: Aphididae) [40,42]. As exemplified with our findings, some invasive insect pests have been found to depend on clonal population structures to successfully invade and multiply in a broad range of niches e.g., [43]. A highly specialized clonal genotype of a strictly asexual population of the pea aphid, Acyrthosiphon pisum Harris was the main reason influencing the demographic success of the pest in central Chile [44]. Cifuentes et al. [45] found no genetic variation within populations of the tomato leaf miner, Tuta absoluta Meyrick (Lepidoptera: Gelechiidae) spreading through South America and ultimately reaching the Mediterranean Basin, with only a single genetic type. Likewise, a well-established invasive population of the oleander aphid, Aphis nerii Fonscolombe were reported as having extremely low genetic diversity in the southern United States, with a “superclone” population supposedly obligatorily asexual [46]. Caron et al. [47] reported two widespread, invasive and strictly parthenogenetic “superclones” of the sawfly, Nematus oligospilus Forster (Hymenoptera; Tenthredinidae) dominating willows in three countries in the southern hemisphere i.e., South Africa, New Zealand, and Australia.

4.2. Demographic History: Resurgence of Resident Population or Recent Introduction?

Based on historical reports, Lever [14] claimed that the more invasive coconut scale insect A. destructor rigidus (now A. rigidus) reported by Reyne [9] in Indonesia was also present in the Philippines. However, Velasquez [15] did not observe its occurrence across the archipelago and reported a highly probable confinement of the pest in the southern region of the country. Given this historical evidence, we assume that the southern populations have existed for a long time. This scenario might have supported our dispersal model for the mitochondrial sequences, except that the first reported sighting of A. rigidus in Tanauan, Batangas, Luzon was in 2009 [10] with no historical evidence of resident populations in the past. This suggests that the northern populations were most probably a recent introduction event from a different source.
A Bayesian search strategy was performed to assess the migration history between the northern and southern populations. However, our results from the mitochondrial and nuclear datasets are difficult to reconcile. We reiterate that in this historic gene-flow analysis, the two groups were assumed to have reached mutation-migration-drift equilibrium. Despite the contrasting results, both models indicate that sequences from the two groups do not belong to one panmictic population. Even so, given the difference in the divergence or migration pattern of the models for each marker, we were unable to unequivocally infer the source of each outbreak population. Methodological assumptions [48] in the program Migrate-n, just like other coalescent-based approaches, did not take into account another source of migrants, or that ancestral variation may come from populations that were not considered in the analysis. Hence, the inference of the possible source of the northern outbreak population needs further exploration.
On the other hand, lower genetic variation is expected for younger populations due to founder effects and genetic bottlenecks during colonization and establishment [49]. Invasive or recently introduced species have been reported to exhibit reduced genetic variation e.g., [50,51]. Introduced populations are usually small so decreased genetic diversity is expected, and are often less variable than the source population which contributes to the invasive success of the species [45]. The nuclear marker revealed the existence of three southern haplotypes, with samples from Basilan having two distinct haplotypes. Genetic variation amongst the populations was very low and the Zamboanga EF-1α sequences are differentiated compared with that from Basilan with few nucleotide substitutions (Figure 3B).
Genetic variation was already low amongst the southern populations for the nuclear DNA, but in comparison to the uniformly genetic northern population (Table 2 and Table 3), our findings infer that the southern A. rigidus population was relatively older than the northern one. Alongside previous historical reports, the level of genetic variation between the geographic regions supports our hypothesis of an existing resident A. rigidus population in the southern part of the Philippines. Local insect populations have the potential to outbreak due to anthropogenic and environmental changes [1,3]. Similar observations on insect pests have been reported in the literature. A notable example by Kobayashi et al. [8] reveals that the multiple nationwide outbreaks of the native populations of the mirid bug, Stenotus rubrovittatus Matsumura (Heteroptera: Miridae) in Japan were induced by changes in the agro-ecosystem without invasion of populations from other areas. Populations of the pest were also genetically isolated by distance and separated into genetic clusters occupying spatially segregated regions. Additionally, temporal fluctuations of pest insects in agroecosystems could be driven by various factors [2]. Pesticide application may induce the resurgence of native pest insect populations by reducing the abundance of natural enemies or by the removal of competitive species in the area e.g., [52,53]. Weather conditions can also trigger insect outbreaks due to the dramatic changes in pest abundance. Ward and Aukema [54] reported that the cyclic outbreaks of the native tree-killing bark beetle, Dendroctonus simplex LeConte (Coleoptera: Curculionidae) on Tamarack (Larix laricina (Du Roi) K. Koch) in Minnesota, USA are climate-driven, being specifically associated with warmer and dryer years, more likely in areas with prior defoliation. Schwartzberg et al. [55] simulated climate warming and observed warming-induced phenological shifts in the forest tent caterpillar, Malacosoma disstria Hübner (Lepidoptera: Lasiocampidae) in relation to the phenology of its host trees. These findings illustrate the mechanisms by which anthropogenic and climatic changes induce outbreaks from native insect pests.

4.3. Future Prospects: Applications for Biological Control

The information regarding the amount of genetic variation in pest species could be useful because insects with reduced genetic variation have been reported to be more easily controlled than those with more genetic diversity [56]. For instance, populations with multiple resistant phenotypes are more difficult to control via pathogenic biological control agents compared to populations lacking diversity [57]. In 2014, Almarinez et al. [58] discovered a new encyrtid species, Comperiella calauanica Barrion, Almarinez, Amalin (Hymenoptera: Encyrtidae) parasitizing A. rigidus infesting the coconut palms in a forest-based resort in Calauan, Laguna (Luzon, Philippines) with no history of chemical insecticide applications. The parasitoid has been observed to be host-specific, with up to 80% parasitisation rate. The use of C. calauanica for biological control implementations through conservation and augmentation in the field have been conducted in the Philippines [58,59]. Utilizing the natural variation of the pest insect, or lack thereof, would be vital in the development and implementation of biological control strategies for integrated pest management programmes against this coconut scale insect pest.

5. Conclusions

The current opinion about the origin of the coconut scale insect outbreak in the Philippines was a recent introduction of A. rigidus from other countries of native range and spread via wind dispersal or importation of infested planting material from the northern region to the south given the timeline of the outbreak reports. However, our results strongly suggest the separation of two distinct groups, the northern and southern A. rigidus from the outbreak populations collected from 2014 to 2017 in the Philippines. Very low or no genetic differentiation was observed within and amongst the populations per geographic region, indicating two unrelated outbreak events of the pest species originating from two genetically uniform or “superclone” populations currently isolated in each respective region. Historical data supports our assumption on the current resurgence of an established A. rigidus population in the south. Given the lack of historical information supporting the existence of an established A. rigidus populations in the northern region, we disregard the possible resurgence of a persisting native population and suggest that the outbreak possibly resulted from a recent introduction of a non-native population. Assessment of the possible source population of the northern outbreaks needs further exploration.
The use of mtCOI and the nuclear EF-1α markers showed no or very low genetic differentiation for all A. rigidus populations. Other robust and more informative genetic markers such as microsatellites could provide further genetic information during study of the invasive coconut scale population. Further studies should also include more expansive sampling, taking into consideration other possible sources of A. rigidus such as Indonesia [10] and Vietnam [24]. This would provide a more robust and stringent population and gene flow estimation of A. rigidus in the Philippines. Nevertheless, our findings provided an initial important genetic basis and information for designing and implementing biological control strategies against the invasive coconut scale insect pest A. rigidus in the Philippines.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2075-4450/10/11/374/s1, Figure S1: Phylogeny of the outbreak Aspidiotus rigidus populations in the Philippines based on mtCOI and EF-1α Table S1: Summary of the computations of the parameters for demographic expansion for the EF-1α dataset, Table S2: Population pairwise FST values between the outbreak A. rigidus populations for mtCOI and EF-1α data, Table S3: Bayes factors and log marginal likelihoods.

Author Contributions

Conceptualization, K.W., D.M.A., T.M.C., B.J.M.A. and J.M.S.; fieldwork, B.J.M.A. and D.M.A.; laboratory work, N.I.; formal analysis, N.I. and J.M.S.; data curation, T.M.C., writing—original draft preparation, J.M.S.; writing—review and editing, A.T.B., D.M.A. and K.W.; project administration and supervision, D.M.A. and K.W.; funding acquisition, K.W.

Funding

This research was funded by the Sumitomo Electric Industries Group Corporate Social Responsibility Foundation, Ehime University’s Research Unit Program, and the Japan Society for the Promotion of Science (JSPS) Core-to-Core program (Asia-Africa Science Platforms).

Acknowledgments

We are grateful to Leslie Ann Ormenita and Reynaldo Majaducon of the Biological Control Research Unit of the Center for Natural Sciences and Environmental Research of De La Salle University for their technical support during the field work and preparation of samples. Many thanks also to Maribet Gamboa for her valuable suggestions on the bioinformatics analyses conducted in this study.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Berryman, A.A. The theory and classification of outbreaks. In Insect Outbreaks; Academic Press, Inc.: San Diego, CA, USA, 1987; pp. 3–30. [Google Scholar]
  2. Risch, S.J. Agricultural ecology and insect outbreaks. In Insect Outbreaks; Academic Press, Inc.: San Diego, CA, USA, 1987; pp. 217–238. [Google Scholar]
  3. Ziska, L.H.; Blumenthal, D.M.; Runion, G.B.; Hunt, E.R.; Diaz-Soltero, H. Invasive species and climate change: An agronomic perspective. Clim. Chang. 2011, 105, 13–42. [Google Scholar] [CrossRef]
  4. Handley, L.J.L.; Estoup, A.; Evans, D.M.; Thomas, C.E.; Lombaert, E.; Facon, B.; Aebi, A.; Roy, H.E. Ecological genetics of invasive alien species. BioControl 2011, 56, 409–428. [Google Scholar] [CrossRef] [Green Version]
  5. Strayer, D.L.; D’Antonio, C.M.; Essl, F.; Fowler, M.S.; Geist, J.; Hilt, S.; Jarić, I.; Jöhnk, K.; Jones, C.G.; Latzka, A.W.; et al. Boom-bust dynamics in biological invasions: Towards an improved application of the concept. Ecol. Lett. 2017, 20, 1337–1350. [Google Scholar] [CrossRef] [PubMed]
  6. Prentis, P.J.; Wilson, J.R.; Dormontt, E.E.; Richardson, D.M.; Lowe, A.J. Adaptive evolution in invasive species. Trends. Plant. Sci. 2008, 13, 288–294. [Google Scholar] [CrossRef]
  7. Renault, D.; Laparie, M.; McCauley, S.J.; Bonte, D. Environmental adaptations, ecological filtering, and dispersal central to insect invasions. Annu. Rev. Ent. 2018, 63, 345–368. [Google Scholar] [CrossRef]
  8. Kobayashi, T.; Sakurai, T.; Sakakibara, M.; Watanabe, T. Multiple origins of outbreak populations of a native insect pest in an agro-ecosystem. Bull. Entomol. Res. 2011, 101, 313–324. [Google Scholar] [CrossRef]
  9. Reyne, A. Studies on a serious outbreak of Aspidiotus destructor rigidus in the coconut-palms of Sangi (North Celebes). Tijdschr. Entomol. 1948, 89, 83–123. [Google Scholar]
  10. Watson, G.W.; Adalla, C.B.; Shepard, B.M.; Carner, G.R. Aspidiotus rigidus Reyne (Hemiptera: Diaspididae): A devastating pest of coconut in the Philippines. Agric. For. Entomol. 2015, 17, 1–8. [Google Scholar] [CrossRef]
  11. Reyne, A. Notes on the biology of Comperiella unifasciata Ishii and its host Aspidiotus destructor rigidus nov. subspec. Tijdschr. Entomol. 1947, 88, 294–302. [Google Scholar]
  12. Manohar, E.C. Infestation Status, IPM Initiatives and Recovery Risk Assessment Protocol. In Proceedings of the NSTW Agri-Aqua Forum: Coping up with Cocolisap, SMX Convention Center, Mall of Asia, Pasay City, Philippines, 26 July 2015. [Google Scholar]
  13. Javier, E. Science to the Rescue of Coconut Industry. Philippine Daily Inquirer. 2014. Available online: http://opinion.inquirer.net/76268/science-to-the-rescue-of-coconut-industry (accessed on 8 February 2017).
  14. Lever, R.J.A.W. Pests of Coconut Palm, 2nd ed.; United Nations Food & Agriculture Organization: Rome, Italy, 1979. [Google Scholar]
  15. Velasquez, F.J. Some Philippine armored scale insects of the tribe Aspidiotini (Diaspididae, Homoptera). Philipp. Entomol. 1971, 2, 89–154. [Google Scholar]
  16. Wilby, A.; Thomas, M.B. Natural enemy diversity and pest control: Patterns of pest emergence with agricultural intensification. Ecol. Lett. 2002, 5, 353–360. [Google Scholar] [CrossRef]
  17. Rugman-Jones, P.F.; Hoddle, M.S.; Phillips, P.A.; Jeong, G.; Stouthamer, R. Strong genetic structure among populations of the invasive avocado pest Pseudacysta perseae (Heidemann) (Hemiptera: Tingidae) reveals the source of introduced populations. Biol. Invasions 2002, 14, 1079–1100. [Google Scholar] [CrossRef]
  18. Kébé, K.; Alvarez, N.; Espíndola, A.; Justy, F.; Olivieri, I.; Sembène, M. Insights into the genetic structure of the cowpea pest Callosobruchus maculatus in Africa. J. Pest. Sci. 2016, 89, 449–458. [Google Scholar] [CrossRef]
  19. Yang, W.Y.; Tang, X.T.; Ju, R.T.; Zhang, Y.; Du, Y.Z. The population genetic structure of Corythucha ciliata (Say) (Hemiptera: Tingidae) provides insights into its distribution and invasiveness. Sci. Rep. 2017, 7, 635. [Google Scholar] [CrossRef] [PubMed]
  20. Zhang, L.; Cai, W.; Luo, J.; Zhang, S.; Li, W.; Wang, C.; Lv, L.; Cui, J. Population genetic structure and expansion patterns of the cotton pest Adelphocoris fasciaticollis. J. Pest. Sci. 2018, 91, 539–550. [Google Scholar] [CrossRef]
  21. Provencher, L.M.; Morse, G.E.; Weeks, A.R.; Normark, B.B. Parthenogenesis in the Aspidiotus nerii complex (Hemiptera: Diaspididae): A single origin of a worldwide, polyphagous lineage associated with Cardinium bacteria. Ann. Entomol. Soc. Am. 2005, 98, 629–635. [Google Scholar] [CrossRef]
  22. Andersen, J.C. A Phylogenetic Analysis of Armored Scale Insects, based upon Nuclear, Mitochondrial, and Endosymbiont Gene Sequences. Master’s Thesis, University of Massachusetts Amherst, Amherst, MA, USA, 2014. [Google Scholar]
  23. Andersen, J.C.; Gruwell, M.E.; Morse, G.E.; Normark, B.B. Cryptic diversity in the Aspidiotus nerii complex in Australia. Ann. Entomol. Soc. Am. 2010, 103, 844–854. [Google Scholar] [CrossRef]
  24. Schneider, S.A.; Okusu, A.; Normark, B.B. Molecular phylogenetics of Aspidiotini armored scale insects (Hemiptera: Diaspididae) reveals rampant paraphyly, curious species radiations, and multiple origins of association with Melissotarsus ants (Hymenoptera: Formicidae). Mol. Phylogenet. Evol. 2018, 129, 291–303. [Google Scholar] [CrossRef]
  25. Park, D.S.; Leem, Y.J.; Hahn, K.W.; Suh, S.J.; Hong, K.J.; Oh, H.W. Molecular identification of mealybugs (Hemiptera: Pseudococcidae) found on Korean pears. J. Econ. Entomol. 2010, 103, 25–33. [Google Scholar] [CrossRef]
  26. Morse, G.E.; Normark, B.B. A molecular phylogenetic study of armoured scale insects (Hemiptera: Diaspididae). Syst. Entomol. 2006, 31, 338–349. [Google Scholar] [CrossRef]
  27. Palumbi, S.R. Nucleic acids II: The polymerase chain reaction. In Molecular Systematics; Oxford University Press: Oxford, England, 1996; pp. 205–247. [Google Scholar]
  28. Park, D.S.; Suh, S.J.; Hebert, P.D.; 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. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Katoh, K.; Standley, D.M. MAFFT multiple sequence alignment software version 7: Improvements in performance and usability. Mol. Biol. Evol. 2013, 30, 772–780. [Google Scholar] [CrossRef] [PubMed]
  30. Talavera, G.; Castresana, J. Improvement of phylogenies after removing divergent and ambiguously aligned blocks from protein sequence alignments. Syst. Biol. 2007, 56, 564–577. [Google Scholar] [CrossRef] [PubMed]
  31. Rozas, J.; Ferrer-Mata, A.; Sánchez-DelBarrio, J.C.; Guirao-Rico, S.; Librado, P.; Ramos-Onsins, S.E.; Sánchez-Gracia, A. DnaSP 6: DNA sequence polymorphism analysis of large data sets. Mol. Biol. Evol. 2017, 34, 3299–3302. [Google Scholar] [CrossRef] [PubMed]
  32. Excoffier, L.; Lischer, H.E. Arlequin suite ver 3.5: A new series of programs to perform population genetics analyses under Linux and Windows. Mol. Ecol. Resour. 2010, 10, 564–567. [Google Scholar] [CrossRef] [PubMed]
  33. Fu, Y.X.; Li, W.H. Statistical tests of neutrality of mutations. Genetics 1993, 133, 693–709. [Google Scholar]
  34. Beerli, P. Comparison of Bayesian and maximum-likelihood inference of population genetic parameters. Bioinformatics 2005, 22, 341–345. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Beerli, P.; Palczewski, M. Unified framework to evaluate panmixia and migration direction among multiple sampling locations. Genetics 2010, 185, 313–326. [Google Scholar] [CrossRef]
  36. Bandelt, H.J.; Forster, P.; Röhl, A. Median-joining networks for inferring intraspecific phylogenies. Mol. Biol. Evol. 1999, 16, 37–48. [Google Scholar] [CrossRef]
  37. Leigh, J.W.; Bryant, D. POPART: Full-feature software for haplotype network construction. Methods Ecol. Evol. 2015, 6, 1110–1116. [Google Scholar] [CrossRef]
  38. Darriba, D.; Taboada, G.L.; Doallo, R.; Posada, D. jModelTest 2: More models, new heuristics and parallel computing. Nat. Methods 2012, 9, 772. [Google Scholar] [CrossRef] [PubMed]
  39. Kozlov, A.; Darriba, D.; Flouri, T.; Morel, B.; Stamatakis, A. RAxML-NG: A fast, scalable, and user-friendly tool for maximum likelihood phylogenetic inference. Bioinformatics 2019. [Google Scholar] [CrossRef] [PubMed]
  40. Abbot, P. A closer look at the spatial architecture of aphid clones. Mol. Ecol. 2011, 20, 4587–4589. [Google Scholar] [CrossRef] [PubMed]
  41. Normark, B.B.; Johnson, N.A. Niche explosion. Genetica 2011, 139, 551–564. [Google Scholar] [CrossRef] [PubMed]
  42. Loxdale, H.D.; Edwards, O.; Tagu, D.; Vorburger, C. Population genetic issues: New insights using conventional molecular markers and genomics tools. In Aphids as Crop Pests, 2nd ed.; van Emden, H.F., Harrington, R., Eds.; CABI: Oxford, UK, 2017; pp. 50–80. [Google Scholar]
  43. Figueroa, C.C.; Simon, J.C.; Le Gallic, J.F.; Prunier-Leterme, N.; Briones, L.M.; Dedryver, C.A.; Niemeyer, H.M. Genetic structure and clonal diversity of an introduced pest in Chile, the cereal aphid Sitobion avenae. Heredity 2005, 95, 24. [Google Scholar] [CrossRef]
  44. Peccoud, J.; Figueroa, C.C.; Silva, A.X.; Ramirez, C.C.; Mieuzet, L.; Bonhomme, J.; Stoeckel, S.; Plantegenest, M.; Simon, J.C. Host range expansion of an introduced insect pest through multiple colonizations of specialized clones. Mol. Ecol. 2008, 17, 4608–4618. [Google Scholar] [CrossRef]
  45. Cifuentes, D.; Chynoweth, R.; Bielza, P. Genetic study of Mediterranean and South American populations of tomato leafminer Tuta absoluta (Povolny, 1994) (Lepidoptera: Gelechiidae) using ribosomal and mitochondrial markers. Pest Manag. Sci. 2011, 67, 1155–1162. [Google Scholar] [CrossRef]
  46. Harrison, J.S.; Mondor, E.B. Evidence for an invasive aphid “superclone”: Extremely low genetic diversity in oleander aphid (Aphis nerii) populations in the southern United States. PLoS ONE 2011, 6, e17524. [Google Scholar] [CrossRef]
  47. Caron, V.; Ede, F.J.; Sunnucks, P. Unravelling the paradox of loss of genetic variation during invasion: Superclones may explain the success of a clonal invader. PLoS ONE 2014, 9, e97744. [Google Scholar] [CrossRef]
  48. Knowles, L.L.; Carstens, B.C.; Keat, M.L. Coupling genetic and ecological-niche models to examine how past population distributions contribute to divergence. Curr. Biol. 2007, 17, 940–946. [Google Scholar] [CrossRef]
  49. Hewitt, G.M. Genetic consequences of climatic oscillations in the Quaternary. Philos. Trans. R Soc. Lond. B Biol. Sci. 2004, 359, 183–195. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Tsutsui, N.D.; Suarez, A.V.; Holway, D.A.; Case, T.J. Reduced genetic variation and the success of an invasive species. Proc. Nat. Acad. Sci. USA 2000, 97, 5948–5953. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Navia, D.; De Moraes, G.J.; Roderick, G.; Navajas, M. The invasive coconut mite Aceria guerreronis (Acari: Eriophyidae): Origin and invasion sources inferred from mitochondrial (16S) and nuclear (ITS) sequences. Bull. Entom. Res. 2005, 95, 505–516. [Google Scholar] [CrossRef]
  52. Lu, Y.; Wu, K.; Jiang, Y.; Xia, B.; Li, P.; Feng, H.; Wyckhuys, K.A.; Guo, Y. Mirid bug outbreaks in multiple crops correlated with wide-scale adoption of Bt cotton in China. Science 2010, 328, 1151–1154. [Google Scholar] [CrossRef]
  53. Bommarco, R.; Miranda, F.; Bylund, H.; Björkman, C. Insecticides suppress natural enemies and increase pest damage in cabbage. J. Econ. Entomol. 2011, 104, 782–791. [Google Scholar] [CrossRef]
  54. Ward, S.F.; Aukema, B.H. Anomalous outbreaks of an invasive defoliator and native bark beetle facilitated by warm temperatures, changes in precipitation, and interspecific interactions. Ecography 2019. [Google Scholar] [CrossRef]
  55. Schwartzberg, E.G.; Jamieson, M.A.; Raffa, K.F.; Reich, P.B.; Montgomery, R.A.; Lindroth, R.L. Simulated climate warming alters phenological synchrony between an outbreak insect herbivore and host trees. Oecologia 2014, 175, 1041–1049. [Google Scholar] [CrossRef]
  56. Ahern, R.G.; Hawthorne, D.J.; Raupp, M.J. Founder effects and phenotypic variation in Adelges cooleyi, an insect pest introduced to the eastern United States. Biol. Invasions 2009, 11, 959. [Google Scholar] [CrossRef]
  57. Burdon, J.J.; Thrall, P.H. Genetic structure of natural plant and pathogen populations. In Genetics, Evolution, and Biological Control; Ehler, L.E., Sforza, R., Mateille, T., Eds.; CABI Publishing: Wallingford, UK, 2003; pp. 1–17. [Google Scholar]
  58. Almarinez, B.J.M.; Amalin, D.M.; Carandang, J.S.R.; Navasero, M.V.; Navasero, M.M. First Philippine record of the parasitoid, Comperiella sp. (Hymenoptera: Encyrtidae): A potential biocontrol agent against Aspidiotus rigidus (Hemiptera: Diaspididae). J. Appl. Entomol. 2015, 139, 237–240. [Google Scholar] [CrossRef]
  59. Navasero, M.M.; Corpuz-Raros, L.A.; Navasero, M.V.; Cayabyab, B.F. Survey of mites (Acari) associated with outbreak populations of Aspidiotus rigidus Reyne (Hemiptera: Diaspididae) on coconut in Calabarzon. Philipp. Entomol. 2015, 29, 104–113. [Google Scholar]
Insects 10 00374 g001 550
Figure 1. Map of the seven localities with reported A. rigidus outbreaks in the Philippines from 2014 to 2017. The insect rearing facility of the Biological Control Research Unit of De La Salle University labeled “AR”. Dots indicate sampling locations. Northern localities: Orani, Bataan (BT), Nagcarlan (NG) and San Pablo (SP), Laguna, Tanauan (TN) and Talisay (TL), Batangas; Southern localities: Basilan (BS) and Zamboanga (ZB). See Table 1 for the more detailed information regarding location and sample collection information.
Figure 1. Map of the seven localities with reported A. rigidus outbreaks in the Philippines from 2014 to 2017. The insect rearing facility of the Biological Control Research Unit of De La Salle University labeled “AR”. Dots indicate sampling locations. Northern localities: Orani, Bataan (BT), Nagcarlan (NG) and San Pablo (SP), Laguna, Tanauan (TN) and Talisay (TL), Batangas; Southern localities: Basilan (BS) and Zamboanga (ZB). See Table 1 for the more detailed information regarding location and sample collection information.
Insects 10 00374 g001
Insects 10 00374 g002 550
Figure 2. Representative adult female A. rigidus from different outbreak areas: (A) Southern Tagalog Region (Laguna, Cavite, and Batangas); (B) Orani, Bataan; (C) Basilan; and (D) Zamboanga City. The arrows point to the egg skins, which for this species is characteristically distributed along the posterior or pygidial half of the insect body.
Figure 2. Representative adult female A. rigidus from different outbreak areas: (A) Southern Tagalog Region (Laguna, Cavite, and Batangas); (B) Orani, Bataan; (C) Basilan; and (D) Zamboanga City. The arrows point to the egg skins, which for this species is characteristically distributed along the posterior or pygidial half of the insect body.
Insects 10 00374 g002
Insects 10 00374 g003 550
Figure 3. Median-joining network of the A. rigidus populations from 305 individuals for the mtCOI gene (A), and 75 individuals for the protein-coding EF-1α gene (B), showing location and frequency of haplotypes. Each circle represents an observed haplotype; circle size indicates the number of individuals observed; the colors correspond to sampling localities. The total number of mutations, Eta presented as hatch marks.
Figure 3. Median-joining network of the A. rigidus populations from 305 individuals for the mtCOI gene (A), and 75 individuals for the protein-coding EF-1α gene (B), showing location and frequency of haplotypes. Each circle represents an observed haplotype; circle size indicates the number of individuals observed; the colors correspond to sampling localities. The total number of mutations, Eta presented as hatch marks.
Insects 10 00374 g003
Table
Table 1. Sampling localities of the outbreak A. rigidus populations. N, number of individuals with mtCOI and EF-1α sequences; H, haplotypes indicated in Figure 3.
Table 1. Sampling localities of the outbreak A. rigidus populations. N, number of individuals with mtCOI and EF-1α sequences; H, haplotypes indicated in Figure 3.
LocalityCodeCollection DatemtCOIEF-1α
NHNH
Northern RegionDLSU-STC, Laguna aARJuly 201735H18H1
Orani, BataanBTSeptember 201521H10—
Nagcarlan, LagunaNGJanuary 201524H111H1
San Pablo, LagunaSPDecember 201421H18H1
Talisay, BatangasTLDecember 201435H118H1
Tanauan, BatangasTNDecember 201429H18H1
Southern RegionIsabela, BasilanBSNovember 201644H216H2; H3
Zamboanga City, ZamboangaZBApril 201796H26H4
aA. rigidus reared on mangosteen at the DLSU-STC BCRU rearing facility from samples collected on the outbreak population in Orani, Bataan.
Table
Table 2. Parameters of genetic diversity and demographic analysis of the two population groups.
Table 2. Parameters of genetic diversity and demographic analysis of the two population groups.
GeneGroupNShHaplotype akHd (SD)Pi (SD)DbFsD *bF *
mtCOINorthern165——H1———————
Southern140——H2———————
All305312—15.450.4980 (0.005)0.0239 (0.00025)5.85 ***59.592.04 **4.43
EF-1αNorthern53——H1———————
Southern2223H2; H3; H40.870.6880 (0.039)0.0009 (0.00010)1.330.990.851.13
All75144—5.710.4770 (0.064)0.0057 (0.00061)2.83 **12.811.55 *2.38
a Haplotype data by DnaSP v. 6.10.04; b Parameters with statistical test: * indicates p < 0.05; ** indicates p < 0.02; *** indicates p < 0.01.
Table
Table 3. Partitioning of genetic variation at different hierarchical levels.
Table 3. Partitioning of genetic variation at different hierarchical levels.
GeneSource of Variationdf.Sum of
Squares		Variance
Components		Percentage
of Variation		Fixation Indices
mtCOIAmong groups12347.86915.50000 Va100FCT = 1.00000 *
Among populations within groups600.00000 Vb0FSC = 0.00000
Within populations29700.00000 Vc0FST = 1.00000 **
EF-1αAmong groups1202.0106.45405 Va97.63FCT = 0.97630 *
Among populations within groups55.1990.09876 Vb1.49FSC = 0.63040 **
Within populations683.9380.05790 Vc0.88FST = 0.99124 **
* indicates p < 0.05; ** indicates p < 0.01.

Share and Cite

MDPI and ACS Style

Serrana, J.M.; Ishitani, N.; Carvajal, T.M.; Almarinez, B.J.M.; Barrion, A.T.; Amalin, D.M.; Watanabe, K. Unraveling the Genetic Structure of the Coconut Scale Insect Pest (Aspidiotus rigidus Reyne) Outbreak Populations in the Philippines. Insects 2019, 10, 374. https://0-doi-org.brum.beds.ac.uk/10.3390/insects10110374

AMA Style

Serrana JM, Ishitani N, Carvajal TM, Almarinez BJM, Barrion AT, Amalin DM, Watanabe K. Unraveling the Genetic Structure of the Coconut Scale Insect Pest (Aspidiotus rigidus Reyne) Outbreak Populations in the Philippines. Insects. 2019; 10(11):374. https://0-doi-org.brum.beds.ac.uk/10.3390/insects10110374

Chicago/Turabian Style

Serrana, Joeselle M., Naoto Ishitani, Thaddeus M. Carvajal, Billy Joel M. Almarinez, Alberto T. Barrion, Divina M. Amalin, and Kozo Watanabe. 2019. "Unraveling the Genetic Structure of the Coconut Scale Insect Pest (Aspidiotus rigidus Reyne) Outbreak Populations in the Philippines" Insects 10, no. 11: 374. https://0-doi-org.brum.beds.ac.uk/10.3390/insects10110374

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop