ORIGINAL RESEARCH article

Front. For. Glob. Change, 09 June 2023

Sec. Pests, Pathogens and Invasions

Volume 6 - 2023 | https://doi.org/10.3389/ffgc.2023.1167541

Taxonomic relationship among four European Physokermes species (Hemiptera: Coccomorpha) based on nuclear and mitochondrial DNA

  • 1. Department of Biology, Lund University, Lund, Sweden

  • 2. Department of Physical Geography and Ecosystem Science, Lund University, Lund, Sweden

Abstract

Introduction:

Scale insects (Hemiptera: Coccomorpha) feed on a wide variety of agricultural crops and forest and ornamental trees worldwide. These pest insects damage plants not only by causing sap loss but also by reducing the plant's photosynthetic activity. This is because the honeydew they produce acts as a substrate for mold, which covers leaf surfaces. In the last decades, several outbreaks of Physokermes spp. (soft-scale insects) have occurred throughout Europe and have partly been attributed to unusual weather conditions or climate change, as some species seem to be expanding their distribution range. However, the small size of these insects and their large intraspecific morphological variation have hindered the identification of the species responsible for outbreaks.

Methods:

In this study, mitochondrial DNA (cytochrome c oxidase subunit I, COI), ribosomal RNA (28S), and nuclear (elongation factor 1α, EF1α) DNA markers were used to reconstruct the phylogenetic relationships of four Physokermes species sampled throughout Europe in 2013–2015.

Results and discussion:

The results allowed us to clearly distinguish P. hellenicus and P. inopinatus from each other, as these appeared in well-supported clades in the phylogenetic trees and from P. piceae and P. hemicryphys. However, P. hemicryphus appeared in a single clade in trees based on 28S and EF1α but among P. piceae in the COI tree. Further investigations are therefore required to determine the taxonomic status of P. piceae and P. hemicryphys, which seem to comprise a species complex.

1. Introduction

Scale insects (Hemiptera: Coccomorpha) are important plant pests feeding on a wide range of agricultural crops as well as on forest and ornamental trees worldwide (Park et al., ). In addition to directly affecting host plants through sap loss, the honeydew produced by scale insects covers leaf surfaces, allowing mold development and thus reducing photosynthetic activity (Kozár et al., ). Because they are small and frequently found living within the plant, scale insects are often unseen. This contributes to their invasive potential and presents a huge challenge for plant quarantine (Andersen et al., ) and a major phytosanitary problem. The genus Physokermes comprises an important group of soft-scale insects infesting coniferous plants in the Holarctic region. In total, 13 species are currently recognized within Physokermes (Papanastasiou et al., ), but only four are known to occur in Europe: the widespread P. hemicryphus (Dalman) and P. piceae (Schrank), the more easterly distributed P. inopinatus (Danzig & Kozár), and P. hellenicus, (Kozár & Gounari) described from Abies cephalonica (Pinaceae) in Greece (Kozár et al., ). Although commonly known as the Hungarian spruce scale, P. inopinatus caused a mass infestation in Sweden in 2010 (McCarthy and Skovsgaard, ), suggesting a considerable northward expansion of this species' distribution, possibly linked to climate change or the unrestricted transport of plants within Europe (Kozár et al., ). Indeed, the attack by P. inopinatus was triggered by a sequence of droughts that contributed to population expansions in the relatively dry summers of 2008–2009, leading to the severe mass infestation of 2010 (McCarthy and Skovsgaard, ; Olsson et al., ; Gertsson and Isacsson, ; Winde et al., ).

However, the identity of the species causing the 2010 outbreak in Sweden was not immediately determined. This was because morphological identification of scale insects requires a high level of expertise and is time consuming. It mostly relies on characters present only in adult females, which is difficult when the morphological variation observed between females of closely related species corresponds to natural intraspecific variation instead of interspecific variation (Malausa et al., ; Park et al., ). The limitations of this methodology have also raised concern about the possible non-identification of cryptic species complexes (Andersen et al., ), and thus, several DNA markers have been developed to distinguish between closely related species of scale insects, mainly mealybugs and armored scales (Kondo et al., ; Andersen et al., ; Gwiazdowski et al., ; Malausa et al., and references therein; Park et al., ; ; Abd-Rabou et al., ; Beltrà et al., ; Deng et al., ; Sethusa et al., ). Such markers include a set of gene regions providing reliable and rapid identification of scale insect species, even when only immature and male specimens are available. Although DNA barcoding based on cytochrome c oxidase subunit I (COI) sequence variation has become a popular tool for species delimitation in several insect groups (Hajibabaei et al., ), and therefore an ideal candidate for identifying scale insects (e.g., Kondo et al., ; Malausa et al., ; Park et al., ; Abd-Rabou et al., ; Beltrà et al., ; Deng et al., ; Sethusa et al., ), only 2,172 COI sequences of Coccidae are deposited in GenBank (http://www.ncbi.nlm.nih.gov; accessed November 8, 2022). Gene 28S has also been used to identify insect species (Campbell et al., ; Monaghan et al., ; Smith et al., ). Although lacking sufficient variation to delimitate some species (Park et al., ; Deng et al., ), gene 28S has been proposed as a complementary marker to COI in scale insects (Sethusa et al., ). This supports van Nieukerken et al. () suggestion of including at least one nuclear marker in addition to COI in studies concerning cryptic lineages. Thus, in this study, three genes from three independent loci were amplified: COI (mitochondrial DNA), 28S (D2 and D3 regions, ribosomal RNA gene), and elongation factor 1alpha (EF1α, nuclear DNA). These genes represent a range of different rates of evolution: COI is a relatively fast-evolving mitochondrial gene (Simon et al., ), 28S is relatively conserved and represented in the nuclear genome as numerous tandem repeats (Hoy, ), and EF1α is a low-copy number nuclear protein-coding gene with conserved exons and variable introns (Hoy, ).

The DNA sequence analysis performed in this study for several Physokermes spp. populations across Europe aimed to (i) scan for cryptic species within taxonomically recognized Physokermes species; (ii) study the phylogenetic relationships between the identified taxa; and (iii) obtain a series of molecular markers suitable for the reliable identification of closely related Physokermes spp. in Europe, especially those with invasive pest status, thus contributing to more efficient pest management and quarantine procedures. Finding cryptic species will also provide the necessary background and molecular evidence to conduct thorough morphological studies on Physokermes specimens representative of several species. This will allow laying the foundations for a contemporary taxonomic revision and for the recognition of morphological characters important to species differentiation, similar to what has been done for Chionaspis spp. by Vea et al. () following the work of Gwiazdowski et al. (). To date, only one molecular study of Physokermes spp. has been published, and it included COI sequence analysis of P. hellenicus from different Greek populations (Papanastasiou et al., ).

2. Materials and methods

2.1. Taxon sampling

In total, 203 P. inopinatus, P. piceae, P. hemicryphus, and P. hellenicus individuals were collected across their distribution ranges in Europe, i.e., from Turkey, in the southeast, to Sweden, in the northwest (Table 1). Sampling occurred in 2013–2015, and insects were immediately placed in 95% ethanol and kept at −20°C for DNA preservation. Taxonomic identity was determined by experienced scale insect specialists (see Acknowledgments) based on discriminating morphological characters published in reference studies (e.g., Kosztarab and Kozár, ; Kozár et al., ).

Table 1

SpeciesLocationCountryCoordinatesHostNSample code in
Phylogenetic trees		Sequences (number)
COI28SEF1α
P. inopinatusSödra SandbySwedenPicea glauca var. conica20Pino1.n; Pino2.n121317
P. piceaeMålillaSweden57.337920, 15.754570Picea abies30Ppic1Mn-Pic6Mn826072
TilburgThe Netherlands51.540504, 5.016668Picea abies2Ppic4.n
AncaraTurkeyNot availablePicea excelsa20Ppic5.n-Ppic6.n
Szabadság Square, BudapestHungary47.528272, 18.954156Picea abies4Ppic2.n
VillányHungaryNot availablePicea abies8Ppic3.n
KaszóHungary46.320881, 17.223469Picea abies4Ppic21.n
NagykanizsaHungary46.401989, 17.053061Picea abies6Ppic24.n
BosárkányHungary47.689524, 17.251699Picea abies6Ppic27.n
MosonmagyaróvárHungary47.884160, 17.25471Picea abies7Ppic28.n
VágáshutaHungary48.423353, 21.5394Picea abies7Ppic38.n
MikóházaHungary48.466467, 21.588708Picea abies6Ppic40.n
SzékelyvajaRomania46.440833, 24.64866Picea abies6Ppic42.n
LukailencfalvaRomania46.463667, 24.531167Picea abies5Ppic43.n
FarádHungary47.604667, 17.204583Picea abies6Ppic47.n
BerzenceHungary46.206400, 17.151934Picea abies6Ppic48.n
BalavásárRomania46.409167, 24.704167Picea abies7Ppic41.n
P. hemicryphusLundSweden55.714360, 13.196340Picea glauca var. conica24Phem1.n, Phem2.n, Phem3.n, Phem4.n, Phem2L, Phem5L, Phem6L14164
MålillaSweden57.337920, 15.754570Picea abies2Phem1M, Phem2M
Gellért Hill, BudapestHungaryNot availableNot available15Phem5.n, Phem6.n
P. hellenicusMount ParnithaGreece38.172806, 23.731111Abies cephalonica10Phel1.1n957
Total20311794100

Collection sites and the number of Physokermes spp. individuals examined in the present study, and the number of cytochrome oxidase I (COI), 28S, and elongation factor 1 α (EF1α) sequences obtained.

2.2. DNA extraction, amplification, and sequencing

Total genomic DNA was isolated from the 203 individuals using the QIAamp? DNA Mini Kit (QIAGEN, Hilden, Germany), following the manufacturer's protocol for the purification of total DNA from insects. The nucleotide sequences including the barcoding region of COI, 28S, and EF1α were amplified via PCR using the 5′-tailed primers and profiles described in Table 2. Amplification reaction mixtures (total volume = 20 μL) contained 10 μl DreamTaqTM PCR Master Mix (FermentasTM, Waltham, MA, USA), 1 μl each primer, 6 μl DNA-free water, and 2 μl template DNA. To check for contamination, negative controls (sterile water instead of template DNA) were included in each reaction array. The resulting PCR products were then visualized under UV light in 1% agarose gels stained with SYBRTM Safe DNA Gel Stain (ThermoFisher Scientific, Waltham, MA, USA). Those representing successful reactions were prepared for sequencing by digesting 5 μl of the PCR mixture with 0.5 μl Exonuclease I, and 1 μl FastAPTM Thermosensitive Alkaline Phosphatase (both from FermentasTM) during a thermal cycle consisting of 15 min at 37°C and 15 min at 85°C. PCR products were then sequenced in both directions, to ensure high-quality reads, using PCR primers (Table 2) and the BigDye Terminator Cycle Reaction Kit (Applied Biosystems, Foster City, CA, USA).

Table 2

RegionPrimerSequenceSourcePCR
COIC1-J-2753ywrGTAAACCTAACATTTTTYCCWCARCAGwiazdowski et al., 1′ @ 95°C; 30′' @ 95°C, 1′ @ 40-50°C, 2′ @72°C;
annealing T was decreased 2°C every 2 cycles from 50°C to the final T of 40°C;
10′ @72°C
COIC2-N-3662CCACAAATTTCTGAACATTGACCGwiazdowski et al.,
COI barcodePcoF1CCTTCAACTAATCATAAAAATATYAGPark et al., 2′ @ 95°C; 5x (40′' @ 94°C, 40′' @ 45°C, 1′ @72°C);
40x (40′' @ 94°C, 40′' @ 51°C, 1′ @ 72°C);
5′ @72°C
COI barcodeLepR1TAAACTTCTGGATGTCCAAAAAATCAPark et al.,
D2-D3 28Ss3660GAGAGTTMAASAGTACGTGAAACMorse and Normark, @ 95°C; 30′' @ 95°C, 1′ @58°C, 2′ @ 72°C;
annealing T was lowered 2°C every 3 cycles until final annealing of 42°C,
which was used for 18 Cycles; 10′ @ 72C
D2-D3 28S28bTCGGAAGGAACCAGCTACTAMorse and Normark,
EF1αEF1-a(a)GATGCTCCGGGACAYAGAMorse and Normark,
EF1αEF2ATGTGAGCGGTGTGGCAATCCAAMorse and Normark,

Primers and PCR profiles used to amplify cytochrome oxidase I (COI), 28S, and elongation factor 1 α (EF1α) gene fragments in Physokermes spp.

2.3. Phylogenetic analyses

Contiguous sequences were assembled in Geneious® 8.1.15 (Biomatters, Auckland, New Zealand), aligned using the default parameters and cost matrix of the Geneious® alignment algorithm as included in Geneious® 8.1.15 and edited by eye to maximize blocks of sequence identity. Sequences were checked for homology by comparison with those deposited in GenBank, using the Basic Local Alignment Search Tool (BLAST) available on the NCBI website (NCBI, http://www.ncbi.nlm.nih.gov/Genbank). Translation into amino acids was also performed to confirm data correctness.

Using PartitionFinder 2 (Lanfear et al., ), the best partitioning scheme and model of nucleotide substitution were selected for each dataset based on the Bayesian information criterion. All possible combinations of codon positions were considered as candidate partitions (i.e., from a single partition for the entire sequence to each codon position treated as a partition), but only the nucleotide substitution models supported by MrBayes 3.2.6 (Ronquist and Huelsenbeck, ) were included in the candidate model pool. This procedure was repeated five times to ensure correct estimation. The resulting models and partition schemes were used in phylogenetic reconstructions based on maximum likelihood (ML) and Bayesian Inference (BI), performed in Garli 2.1 (Zwickl, ) and MrBayes 3.2.6, respectively. Because considerable statistical inconsistency can be found when using the maximum parsimony method of phylogenetic reconstruction on this type of data (different rates of evolution), it was not performed in this study.

In Garli, clade stability was evaluated by non-parametric bootstrapping (five independent search replicates of 10,000 generations each). In MrBayes, and starting from a random tree, alternative tree topologies were evaluated through the estimation of posterior probabilities by Markov chain Monte Carlo re-sampling. These were performed in two runs of four chains each (with the default settings implemented in MrBayes) for 10 × 106 generations with a sampling frequency of 100 (average standard deviation of split frequencies <0.01). These samples were then used to generate 50% majority-rule consensus trees. Commands “sumt” and “burnin = 5,000” were set to discard trees that did not reach the likelihood score detected with the “sump” command. Only bootstrap and posterior probabilities values >80% were considered to support clades.

The species delimitation plugin implemented in Geneious® was used to evaluate different scenarios of lineage divergence within Physokermes spp. based on a neighbor-joining (NJ) tree constructed for COI sequences. The distance-based species divergence criterion relies on the postulate that intraspecific genetic variation is small relative to interspecific variation. In this study, two clusters containing a pair of nearest neighbors with >2% distance from each other were considered different species, as this threshold has been proven efficient to discriminate other insect groups (Foottit et al., ; Zhu et al., ). Only monophyletic clades showing bootstrap support >80% were tested.

3. Results and discussion

The BLAST of all sequences against the NCBI database revealed that some (most of them corresponding to 28S sequences) matched those of internal parasites, parasitoids, or microorganisms infecting Physokermes spp.; hence, these sequences were excluded from further analysis. This higher-than-expected ratio of non-specific amplification might be due to the binding of one or both PCR primers to a sequence other than the target one, which is frequently a problem in the amplification of conserved sequences, such as ribosomal DNA. The 117 COI sequences (total length = 712 bp) and the 94 28S sequences (total length = 842 bp) presented the highest sequence homology to Physokermes jezoensis Siraiwa isolate S2, whereas the 100 EF1α sequences (total length = 1,056 bp) presented the highest sequence homology to Paralecanium sp. (Coccomorpha) isolate D0764B. These BLAST hits are probably the result of the very low number of Physokermes species nucleotide sequences deposited in GenBank (n = 22, Available online at: http://www.ncbi.nlm.nih.gov; accessed November 8, 2022). In addition, the handling of specimens for morphological identification might have damaged the DNA quality leading to amplified fragments shorter than ideal for BLAST hits. Nevertheless, the most similar sequences were used as outgroups in the subsequent phylogenetic analyses. The sequences obtained and analyzed in the present study are deposited in GenBank under accession numbers OQ433782–898 (COI), OQ470427–526 (EF1α), and OQ657360–453 (28S).

According to PartitionFinder, the best nucleotide substitution model for Physokermes spp. COI evolution was a generalized time reversible (GTR) with a gamma-distributed rate variation (G) among sites, whereas 28S evolved according to the Hasegawa–Kishino–Yano (HKY) model with a proportion of invariable sites (I). A partitioned model of evolution was found for EF1α sequences, with the first two positions evolving according to Felsenstein's 1981 model (F81) with G and the third according to HKY. The ML and BI analyses performed for all datasets (Figure 1) revealed similar results, although BI usually presented higher support values, and therefore, these are the trees presented here. One clade comprising P. hellenicus and another comprising P. inopinatus were clearly evidenced and well-supported in all trees. However, only 28S and EF1α BI trees evidenced clades comprising P. hemicryphus alone (Figures 1B, C, respectively), with this species appearing mostly among P. piceae and P. inopinatus. In addition, P. hemicryphus samples from Hungary (Phem5.n and Phem6.n) appeared in a different clade from P. hemicryphus collected in Sweden (all other Phem codes). These populations do not seem to constitute host races, because they were both sampled from Picea abies, but might have been isolated by distance. Thus, the 28S nuclear gene is complementary to COI in Physokermes spp. identification, as proposed by Sethusa et al. ().

Figure 1

The NJ tree produced for the COI sequences (Figure 2) also evidenced well-supported clades containing P. hellenicus and P. inopinatus, and three other clades corresponding to P. piceae from Kaszó in Hungary (Ppic21.1 and Ppic21.2) and from Balavásár (Ppic41.2 and Ppic41.3) and Lukailencfalva (P.pic43.1, Ppic43.2, and Ppic43.3) in Romania. All these clades were monophyletic and their distance to the nearest neighbor varied from 3.5% (between P. piceae clades from Romania) to 6.3% (between P. inopinatus and P. hellenicus). These results support the distinction between P. inopinatus and P. hellenicus, as indicated in the ML and BI analyses. Moreover, the genetic distance found is similar to the congeneric 5.8% average mitochondrial DNA divergence (Pons et al., ) and within the 0.0–7.95% range (Hebert et al., ) reported for insects. The clade support and genetic distances found within P. piceae from Hungary and Romania (>2%) suggest the existence of cryptic species. Similar to other scale insects, Physokermes spp. females are flightless, males have a short lifespan, and both sexes have low dispersal capacity, all contributing to population differentiation in different areas and/or host species and promoting cryptic speciation, as observed in Aspidiotus nerii in Australia (Andersen et al., ) and in Chionaspis heterophyllae and C. pinifoliae in North America (Gwiazdowski et al., ). Because all the genetically differentiated P. piceae were sampled from P. abies, host species does not appear to be the factor contributing to such differentiation. Although explanations for the recent outbreaks remain elusive, particularly given the limited dispersal ability of Physokermes spp., changes in climate or weather conditions inducing rapid reproduction of endemic (and thus escaping detection) populations are more likely than sudden range expansions.

Figure 2

Overall, the results obtained here (i) clearly distinguished the closely related P. hellenicus, P. inopinatus, and (P. piceae + P. hemicryphus); (ii) revealed the existence of intraspecific variation and, possibly, of unknown complexes of species; and (iii) confirmed the value of using 28S in addition to COI for identifying scale insects. However, the differentiation of P. hemicryphus suggested in the nuclear data (particularly in 28S) analyses was not supported in the species delimitation analysis, as P. hemicryphus was not genetically distant from P. piceae. However, the sample size of P. hemicryphus and its geographic range was much smaller than that of P. piceae, and thus, further studies are needed to confirm the taxonomic status of P. hemicryphus.

Statements

Data availability statement

The gene sequences for this study can be found in GeneBank at https://www.ncbi.nlm.nih.gov/nuccore/?term=OQ433782:OQ433898[accn] for COI; https://www.ncbi.nlm.nih.gov/nuccore/?term=OQ4704273:OQ470526[accn] for EF1α; https://www.ncbi.nlm.nih.gov/nuccore/?term=OQ657360:OQ657453[accn] for 28S.

Author contributions

OA, AJ, and IW secured funding. IW ensured access to specimens. JM performed the genetic analyses and wrote the first draft of the manuscript. All authors conceived the study, contributed to the manuscript, and approved the submitted version.

Funding

The study is a contribution to the Strategic Research Area Biodiversity and Ecosystem Services in a Changing Climate (BECC), and it was supported by the Carl Trygger Foundation for Scientific Research (CTS 14:17), the Tranemåla Foundation, and the Swedish Research Council Formas (Svenska Forskningsrådet Formas) (2010-822).

Acknowledgments

We would like to thank Éva Szita, Selma Ülgentük, Maurica Jansen, Gyuri Csóka, and Iosif Papanastasiou for generously supplying the insect samples, making our investigation possible and Carl-Axel Gertsson for additional species identity verification and general support to our work.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

Summary

Keywords

Hungarian spruce scale, phylogeny, invasive species, forest pest, climate change

Citation

Marques JF, Winde I, Jönsson AM and Anderbrant O (2023) Taxonomic relationship among four European Physokermes species (Hemiptera: Coccomorpha) based on nuclear and mitochondrial DNA. Front. For. Glob. Change 6:1167541. doi: 10.3389/ffgc.2023.1167541

Received

16 February 2023

Accepted

11 May 2023

Published

09 June 2023

Volume

6 - 2023

Edited by

Quan Lu, Chinese Academy of Forestry, China

Reviewed by

San-an Wu, Beijing Forestry University, China; Zhuhong Wang, Fujian Agriculture and Forestry University, China

Updates

Copyright

© 2023 Marques, Winde, Jönsson and Anderbrant.

This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Olle Anderbrant

†Present address: Inis Winde, Department of Aquatic Sciences and Assessment, Uppsala, Sweden

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Share on WeChat

Scan with WeChat to share this article

Article metrics