Research Article |
Corresponding author: Saša S. Stanković ( sasasta@gmail.com ) Academic editor: Jovana M. Jasso-Martínez
© 2024 Vladimir Žikić, Milana Mitrović, Saša S. Stanković, José L. Fernández-Triana, Maja Lazarević, Kees van Achterberg, Dawid Marczak, Marijana Ilić Milošević, Mark R. Shaw.
This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Citation:
Žikić V, Mitrović M, Stanković SS, Fernández-Triana JL, Lazarević M, van Achterberg K, Marczak D, Milošević MI, Shaw MR (2024) An integrative taxonomic study of north temperate Cotesia Cameron (Hymenoptera, Braconidae, Microgastrinae) that form silken cocoon balls, with the description of a new species. Journal of Hymenoptera Research 97: 255-276. https://doi.org/10.3897/jhr.97.116378
|
Using CO1 sequence analysis, we investigated the relationships of Western Palearctic and Nearctic Cotesia that spin aggregated cocoons in the shape of a ball, and as adults are morphologically very similar. The analysis included the conceptual taxa C. tibialis, C. ofella, C. vanessae, C. ruficrus, C. xylina and C. yakutatensis, as well as the newly described species C. trivaliae sp. nov. The examined specimens of C. tibialis, C. ofella, C. vanessae, C. ruficrus and C. trivaliae sp. nov. were collected in several European countries, and C. xylina and C. yakutatensis in Canada and the USA. Molecular analyses showed that C. ruficrus is not closely related to the other studied taxa. Based on the genetic distances as well as biology and morphology, C. vanessae and C. ofella are confirmed as solid taxa. The species C. yakutatensis comprises two entities. Having 8 haplotypes, C. tibialis also emerges as a species complex, divided into two clusters. With 26 detected haplotypes, C. xylina shows the highest diversity, being composed of three segregates. The conceptual species C. tibialis, C. xylina and C. yakutatensis seem to be species complexes containing several candidates for recognition as distinct species. One from the European C. tibialis complex is here described as new, and the impediments to be overcome before the description of further species are outlined.
DNA barcoding, genetic distance, hosts, species aggregates
Cotesia Cameron, 1891 is a large genus of parasitoid wasps with about 340 species described worldwide (
Cocoon-spinning by Cotesia species usually takes place externally on, under, around, or near the dying caterpillar, which can live several days after the emergence of the parasitoid larvae (
There is a major problem in properly understanding the host repertoires of parasitoids based on the published literature that has been uncritically compiled in sources such as Taxapad (
Morphologically and by cocoon architecture, Cotesia ofella (Nixon, 1974) is very similar to C. tibialis. There are not many published data on the hosts for this parasitoid, but certainly its host repertoire includes noctuid species; in this case, Acronicta aceris (Linnaeus, 1758), A. rumicis (Linnaeus, 1758) (
Another species morphologically similar to C. tibialis is Cotesia berberis (Nixon, 1974). Nothing is known about the biology of this rarely found species, recorded from just three countries (
Unlike the previous three species, C. vanessae is predominantly recorded from some Nymphalidae (Nymphalini) as well as certain Noctuidae. Definite summer hosts are caterpillars of the nymphalids Aglais urticae (Linnaeus, 1758), Vanessa atalanta (Linnaeus, 1758) and V. cardui (Linnaeus, 1758) (
Cotesia ruficrus (Haliday, 1834), a taxon with cosmopolitan distribution recorded in all regions, spins a cocoon mass that more weakly conceals the individual cocoons but overlaps with the above species in the host repertoire. In the north temperate area, this parasitoid is frequently recorded from the pest noctuids Helicoverpa armigera (Hübner, 1808), Leucania loreyi (Duponchel, 1827), Mythimna separata Walker, 1865 and Spodoptera exigua (Hübner, 1808), but it undoubtedly has a much wider host repertoire (MRS, unpublished).
In the Nearctic region, there are at least two more conceptual species that make ball-like cocoon masses, C. xylina (Say, 1836) and C. yakutatensis (Ashmead, 1902). Hosts parasitized by C. xylina include the following species of Noctuinae: Mamestra configurata Walker, 1856, (
Considering the morphological resemblance of these seven nominal Cotesia species, similar architecture of the cocoons, as well as sometimes overlapping host repertoires, it is of interest, and also importance in view of the biocontrol potential of some taxa in this group, to examine their relationships and especially their integrity.
Samples were mostly collected from cocoon masses (Fig.
Photographs were taken using a Leica Flexacam C3 on a Leica M165C stereomicroscope with a magnification of 7.3×. For micrographs, we used a Leica DFC490 camera (Leica Microsystems, Wetzlar, Germany), adapted to a microscope Leica 2500 (Leica Microsystems, Wetzlar, Germany), at a total magnification of 5.0–20.0×. The equipment is in the Laboratory of Zoology at the Faculty of Sciences and Mathematics, Department of Biology and Ecology, University of Niš, Serbia.
Our material is deposited in the collections of the Faculty of Sciences and Mathematics, University of Niš, Serbia; the Naturalis Biodiversity Center, Leiden, The Netherlands; the National Museums of Scotland, Edinburgh, UK; and the Canadian National Collection of Insects, Arachnids and Nematodes, Agriculture and Agri-Food Canada, Ottawa, Ontario, Canada. Many of the specimens whose sequences were obtained from public databases are housed in depositories less accessible to us and, partly for that reason but also because many lack host data, a detailed morphological analysis of all segregates has been postponed until such a time that it can be done in conjunction with both DNA barcoding and host data.
Genomic DNA was extracted from whole specimens of parasitoids using Dneasy Blood and Tissue Kit (Qiagen Inc., Valencia, CA) according to the manufacturer’s instructions. For material reared from hosts, caterpillar identity was confirmed by DNA barcoding following the emergence of parasitoids. The genomic DNA used for analysis was extracted from the caterpillar head.
Taxonomic relationships among the relevant species within the genus Cotesia were investigated using sequence data of the barcoding region of the mitochondrial cytochrome oxidase subunit I (CO1). Standard primer pairs LCO1490/HCO2198 in combination with other primers (Table
The list of primers used to retrieve barcoding fragments of CO1 in Cotesia samples.
Primer name | 5’ primer sequence 3’ | Primer direction | Reference |
---|---|---|---|
LCO1490 | GGTCAACAAATCATAAAGATATTGG | Forward |
|
HCO2198 | TAAACTTCAGGCTGACCAAAAAATCA | Reverse | |
LCO1490puc | TTTCAACWAATCATAAAGATATTGG | Forward |
|
HCO2198puc | TAAACTTCWGGRTGWCCAAARAATCA | Reverse | |
LCO1490hem | TTTCAACTAAYCATAARGATATYGG | Forward |
|
HCO2198hem | TAAACYTCDGGATGBCCAAARAATCA | Reverse | |
Aph2Fd | ATAATTGGWGGATTTGGWAATTG | Forward |
|
Aph2Rd | GTWCTAATAAAATTAATWGCWCC | Reverse | |
Lys2Rd | GTWCTAATAAAATTAATTGCHCC | Reverse | |
Pr2Fd | ATAATTGGAGGRTTTGGWAATTG | Forward |
In addition, 63 sequences of CO1 barcoding fragments of Cotesia parasitoids were obtained from the public databases GenBank (https://www.ncbi.nlm.nih.gov/genbank/) and the Barcoding of Life Data Systems (BOLD, http://www.boldsystems.org/) and included in the molecular analysis (Suppl. material
Sequenced CO1 fragments were manually edited in FinchTV ver. 1.4.0 (https://digitalworldbiology.com/FinchTV) and aligned using the ClustalW program integrated into MEGA5 (
A condensed Bayesian tree inferred from the CO1 barcoding fragments of Cotesia specimens. Bayesian posterior probabilities are shown above branches; scale bar indicates substitutions per site (0.02). Potential scale reduction factors (PSRF) were all approximately equal to one. Description of Cotesia CO1 barcode sequences included in the analysis is given in Suppl. material
A Bayesian tree inferred from the Cotesia CO1 barcoding haplotypes. Bayesian posterior probabilities are shown above branches; scale bar indicates substitutions per site (0.02). Potential scale reduction factors (PSRF) were all approximately equal to one. Description of Cotesia haplotypes is given in Suppl. material
CO1 sequences were aligned, trimmed to the same size of 501 bp (a compromise to achieve uniformity for a large number of sequences), and analysed for haplotype diversity and evolutionary distances. The haplotype data file was generated using DnaSP ver. 5.10.01 (
Terminology of body characters follows van
Preliminary morphological sorting into possible taxa is reflected in Suppl. material
From the trimmed CO1 sequences, 56 haplotypes were determined with 108 variable sites detected (Suppl. material
All three Bayesian inference trees (Figs
Median-joining network designed for 56 CO1 haplotypes of different Cotesia segregates. Black dots are median vectors representing the missing unsampled intermediary haplotype(s). Mutational steps are marked with short black lines. Haplotypes H1, H2 (yellow circles) = ‘ruficrus’; H3, H4 (light green circles) = ‘vanessae’; H5 (grey circle) = ‘cf. tibialis white cocoons’; H6, H7 (orange circles) = ‘trivaliae sp. nov.’; H8–H12 (blue circles) = ‘tibialis 1’; H13–H15 (light blue circles) = ‘tibialis 2’; H16, H17 (purple circles) = ‘ofella’; H18–H20 (pink circles) = ‘xylina 1’; H21–H39 (turquois circles) = ‘xylina 2’; H40–H43 (white circles) = ‘xylina 3’; H44–H46 (green circles) = ‘yakutatensis 1’; H47–H56 (red circles) = ‘yakutatensis 2’.
Average between-groups genetic distances (%) calculated using the Tamura-Nei method.
Group | [1] | [2] | [3] | [4] | [5] | [6] | [7] | [8] | [9] | [10] | [11] |
---|---|---|---|---|---|---|---|---|---|---|---|
[1] ‘ruficrus’ | |||||||||||
[2] ‘vanessae’ | 11.5 | ||||||||||
[3] ‘cf. tibialis white cocoons’ | 9.3 | 6.8 | |||||||||
[4] trivaliae sp. nov. | 8.7 | 7.1 | 3.0 | ||||||||
[5] ‘tibialis 1’ | 8.6 | 6.5 | 1.4 | 2.3 | |||||||
[6] ‘tibialis 2’ | 7.8 | 6.4 | 1.8 | 2.3 | 1.1 | ||||||
[7] ‘ofella’ | 8.4 | 7.8 | 2.8 | 3.3 | 3.0 | 2.9 | |||||
[8] ‘xylina 1’ | 9.3 | 6.9 | 1.1 | 3.1 | 1.6 | 1.9 | 2.9 | ||||
[9] ‘xylina 2’ | 8.6 | 6.5 | 1.8 | 2.3 | 1.2 | 1.1 | 2.9 | 1.9 | |||
[10] xylina 3’ | 8.6 | 6.5 | 1.8 | 2.0 | 1.2 | 1.1 | 2.7 | 1.9 | 1.2 | ||
[11] ‘yakutatensis 1’ | 7.7 | 6.8 | 3.3 | 3.1 | 3.1 | 2.2 | 4.0 | 3.4 | 3.0 | 3.0 | |
[12] ‘yakutatensis 2’ | 8.5 | 6.9 | 4.3 | 4.4 | 4.1 | 3.2 | 5.1 | 4.4 | 4.0 | 4.0 | 3.4 |
Four sequences determined as C. vanessae originating from Spain (3) and Canada (1) clustered together and formed the ‘vanessae’ group. Two haplotypes, H3 and H4 were identified, with a genetic distance of 0.4% (Fig.
Cotesia cf. tibialis white cocoons’ consisted of sequences extracted from three independent samples originating from Poland. Since all three sequences were identical, they appear on the tree under the same haplotype (H5). The average genetic distance of this group from ‘vanessae’, ‘ruficrus’, ‘Cotesia sp.’ (= trivaliae sp. nov.), ‘ofella’, ‘yakutatensis 1’ and ‘yakutatensis 2’ is relatively large (Table
Sixteen CO1 sequences of specimens from Serbia (3), Austria (3), Slovenia (2), Finland (1), Germany (1) and Poland (6), initially determined as Cotesia cf. tibialis (3) or C. tibialis (12), were grouped within the same clade ‘tibialis 1’. Five haplotypes (H8, H9, H10, H11, H12) were identified with in-group average distance of 0.4%. This group separates from ‘ruficrus’, ‘vanessae’, ‘ofella’, ‘yakutatensis 1’ and ‘yakutatensis 2’ in a range from 3.0 to 8.6% (Table
Eight barcode CO1 sequences originating from Serbia, two from the Netherlands and one from Austria clustered within the ‘tibialis 2’ group (Suppl. material
With only one haplotype (H5), ‘cf. tibialis white cocoons’ is inserted into the ‘tibialis’/‘xylina’/‘ofella’ part of the tree. Comparing the average genetic distances between this group and all ‘tibialis’ and ‘xylina’ segregates, ‘cf. tibialis white cocoons’ is closest to ‘xylina 1’ with an average genetic distance of 1.1%, then to ‘tibialis 1’ with an average genetic distance of 1.4%, and equally distant from the groups ‘tibialis 2’, ‘xylina 2’ and ‘xylina 3’ with an average genetic distance of 1.8%.
Within the ‘ofella’ group there are two haplotypes (H16, H17) that differ among themselves with an average distance of 0.2%. This group of haplotypes is separated from other groups with a range of 2.9 to 8.4% (Suppl. material
Analysis of 32 specimens initially determined as C. xylina, originating from Canada (31) and the USA (1) revealed the separation of barcode sequences into three groups to which the following names were assigned: ‘xylina 1’, ‘xylina 2’ and ‘xylina 3’ (Fig.
Barcode sequences of 16 specimens determined as C. yakutatensis originating from Canada (15) and the USA (1) were analysed. In total, 13 haplotypes were determined (Suppl. material
Seven sequences clustered together revealing a new taxon Cotesia trivaliae sp. nov., described below. This includes three sequences from Poland (host unknown) and one from Slovenia, ex Orthosia sp. (Suppl. material
The new species shares some morphological similarities with C. tibialis, including the shape of the first metasomal tergite (T1) laterally, with a medial keel near the proximal part, and the apical truncation of the hypopygium (features used by
Holotype : Poland ♀; Kampinos National Park, Granica; 05.VI.2018; ex cocoon mass in grassland; V. Žikić leg.; dry mounted. Paratypes: Poland 41 ♀ 29 ♂; Kampinos National Park, Granica; 05.VI.2018; ex same brood (3 ♀ 3 ♂ dry mounted, 3 ♀ 3 ♂ microscopic slide mounted, 1 ♀ barcoded, the rest kept 34 ♀ 23 ♂ in 96% alcohol); Poland 40 ♂; Kampinos National Park, Granica, 05.VI.2018, (1 ♂ barcoded, 3 ♂ dry mounted, 1 ♂ microscopic slide mounted, 35 ♂ in 96% alcohol); “same data as for preceding” 15 ♀ 12 ♂ (1 ♀ barcoded, 3 ♀ 3 ♂ dry mounted, 3 ♀ 1 ♂ microscopic slide mounted, 8 ♀ 8 ♂ in 96% alcohol); “same data as for preceding” 54 ♀ 6 ♂ (1 ♀ barcoded, 3 ♀ 3 ♂ dry mounted, 3 ♀ 1 ♂ microscopic slide mounted, the rest kept in 96% alcohol); V. Žikić leg; Slovenia 1 ♀; Ljubljana; ex Orthosia sp.; 01.VI.2018; Š. Modic leg.; (1 ♀ barcoded, slide mounted); Great Britain 28 ♀ 16 ♂ Scotland, Berwickshire, Foulden ex Orthosia gracilis, collected 10.VI.2017, (17 ♀ 5 ♂ barcoded; 8 ♀ 11 ♂ barcoded; 3 ♀ from a further brood of low emergence; VII.2017; M. R. Shaw leg.; Great Britain 1 ♀ 1 ♂; Scotland, Fife, Fleecefauld, ex O. gracilis collected 14.VII.2012, emerged 04.VIII.2012 (1 ♀ 1 ♂ barcoded); M. R. Shaw leg.; Great Britain 3♀; England, Kent, Swanscombe, ex cocoon mass collected on Phragmites australis 10.IV.2017, emerged 02.V.2017 (3♀ barcoded) M. Jennings leg.
The holotype ♀ of C. trivaliae sp. nov., and paratypes have been deposited in the collection of the Faculty of Sciences and Mathematics, the University of Niš, Serbia. A single female from Slovenia has been deposited in the collection of the Faculty of Sciences and Mathematics, the University of Niš, Serbia. Additionally, 3 ♀ 3 ♂ from Poland, as well as all specimens from Great Britain have been deposited at the Department of Natural Sciences, National Museums of Scotland.
The new species is named in honour of the gothic rock band Trivalia.
The currently known distribution of the new species is Poland, Slovenia and the UK (England and Scotland).
Female: (Fig.
Head
: In frontal view (Fig.
Cotesia trivaliae sp. nov. female A habitus B head frontal view C mesosoma and T1 lateral view D head dorsal view E fore wing F hind wing G mesoscutum dorsal view H ovipositor lateral view I metasoma lateral view J metasoma dorsal view (T1-T3) K propodeum dorsal view L 5th tarsomere of front leg, arrow points spine M hind leg outer face N antenna. Scale bars: 500 µm (A–K, M, N); 100 µm (L).
Mesosoma
: Mesoscutum in dorsal view (Fig.
Legs
: Fore leg spine on 5th tarsomere present, extremely weak; hard to see (Fig.
Wings
: Fore wing length 2.5 mm (range 2.4–2.5 mm), and 2.5 times as long as wide at maximal distances (Fig.
Metasoma
: In dorsal aspect T1 heavily rugose all over, widened behind (Fig.
Colour : Head, mesosoma including tegula, all coxae and metasoma black. Trochanter brownish, rest of leg parts yellowish; hind femora apically with small dark spot. Antennae brown. Palpi yellow. Wing venation distinctly yellowish to brown, metacarp, pterostigma and radialis brown.
Male: (Fig.
Cocoon mass
: (Fig.
Analysis of the barcode sequences of selected parasitoid species from the genus Cotesia shows clear separation of some taxa, as well as close relatedness among a significant number of others. What is common to the presently analysed species of Cotesia that make ball-like cocoon masses is that all of them parasitize members of the family Noctuidae, probably exclusively except for C. vanessae (C. ruficrus is considered to be outside this strict grouping, see below). The Bayesian tree (Figs
Molecular analyses confirmed evident delineation as separate entities for the haplotypes of the specimens determined as C. ruficrus, C. vanessae, C. ofella, C. trivaliae sp. nov. (see below) and, for C. yakutatensis, two entities ‘yakutatensis 1’ and ‘yakutatensis 2’. The large genetic difference between C. ruficrus and the other Cotesia is reflected in as many as 30 mutations from the bifurcation spot, as shown in the haplotype network (Fig.
All the Cotesia species studied here exclusively use noctuid larvae to complete their annual life cycle except C. vanessae which (at least in Europe) parasitises a restricted group of hosts from the Nymphalidae (Nymphalinae) during summer. These, however, do not overwinter in the larval stage, and C. vanessae critically depends, like the others, on overwintering as larvae inside overwintering noctuid larvae that either diapause or feed only sporadically through the winter, especially low-feeding or sub-surface resting species in the subfamilies Noctuinae and Plusiinae. It seems possible that the summer hosts of C. vanessae reflect an extension from its ancestral host repertoire. Cotesia vanessae differs from the other 11 defined groups with a genetic distance in the range of 6.4–7.8%.
In the case of C. yakutatensis specimens of North American origin included in this analysis, there is a clear discrimination between ‘yakutatensis 1’ and ‘yakutatensis 2’ as two separate lineages. Their within-group genetic distances are significantly lower than the between-group distance (Table
DNA barcoding showed no clear discrimination between the specimens from groups ‘tibialis 1’, ‘tibialis 2’, ‘cf. tibialis white cocoons’, ‘xylina 1’, ‘xylina 2’ and ‘xylina 3’. Cotesia tibialis is already considered to be a complex of species with variable morphology of the body, wings and cocoon (e.g.,
There is a definite separation between the groups ‘tibialis 1’ and ‘tibialis 2’. Thus ‘tibialis 1’ is connected with ‘cf. tibialis white cocoons’ and ‘xylina 1’, whereas ‘tibialis 2’ is genetically close to the segregates ‘xylina 2’ and ‘xylina 3’. Based on the admittedly small number of identified caterpillars for the samples included in the analysis, it seems that this separation is reflected by the hosts. While ‘tibialis 1’ was reared from several hosts (the overwintering hosts Noctua interposita (Hübner, 1790) and Xestia xanthographa (Denis & Schiffermüller, 1775) both collected in Austria, and the summer hosts Cucullia chamomillae (Denis & Schiffermüller, 1775) from Finland and Anarta myrtilli (Linnaeus, 1761) from Germany), ‘tibialis 2’ was reared only from Mythimna conigera (Denis & Schiffermüller, 1775) collected in Serbia.
Finally, within Cotesia tibialis (sensu lato), we found a third group ‘cf. tibialis white cocoons’. It stands in a common lineage with ‘xylina 1’ but is distant from the other ‘xylina’ groups and also ‘tibialis 2’, though still not enough to be identified as a potentially new taxon (Table
Cotesia ofella could be discriminated as a separate lineage, being closest to C. trivaliae sp. nov. (Fig.
Seven specimens clustered in the separate group initially marked as ‘Cotesia sp.’ but described above as C. trivaliae sp. nov. are the only ones collected from caterpillars of the genus Orthosia (the single species O. gracilis, which is one of the few in the genus that is low-feeding rather than arboreal). Based on the general morphology of the adult body and the architecture and colour of the cocoon, they were originally identified as, or close to, C. tibialis. They are closely related to the tibialis–xylina groups, but genetically, morphologically, and biologically distinct enough to be described as a new species. It is crucial that we are able to describe this species as new because it differs morphologically from the concept of C. tibialis and synonyms as defined by their type specimens, and also because our reared material presents a consistent host repertoire that (although possibly incompletely known) is distinct. But the situation regarding the other probable new species involved in the C. tibialis, C. xylina and C. yakutatensis complexes is far more problematic. In the first place, we have no molecular data for the relevant type specimens that would be needed to fix the use of available names, and neither permissions nor resources to try to obtain that. Second, biological data for the specimens analysed genetically is almost entirely lacking, so the important consideration of host repertoire differences cannot be addressed at this time. We recognise that rearing fresh material is a necessity for further progress, with next-generation techniques to obtain DNA from the old specimens in collections deemed to have reliable host data a further possibility – and we hope that this preliminary study will prompt and help to enable that. Third, a thorough morphological study of the sequenced specimens (and the relevant types, and reared specimens beyond our own) is necessary; but it is premature to conduct that until we can incorporate more reared material that still needs to be obtained. Thus, our study is held at this preliminary phase, but it is hoped that our results will promote further research that is clearly warranted.
In this extensive work, the following helped us a lot: Aleksandra Petrović, Slobodan Stevčić, Korana Kocić, Jelisaveta Čkrkić, Andjeljko Petrović, Željko Tomanović, Špela Modic. The authors thank Nigel Spring for the photo (Fig.
Supplementary information
Data type: docx
Explanation note: table S1. The list of Cotesia specimens submitted to the molecular analyses. fig. S1. A Bayesian tree inferred from the CO1 barcoding fragments of Cotesia specimens.