Research Article |
Corresponding author: Ulrike Aspöck ( ulrike.aspoeck@nhm-wien.ac.at ) Academic editor: Dominique Zimmermann
© 2017 Horst Aspöck, Ulrike Aspöck, Axel Gruppe, Marcia Sittenthaler, Elisabeth Haring.
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:
Aspöck H, Aspöck U, Gruppe A, Sittenthaler M, Haring E (2017) Anthropogenic dispersal of a snakefly (Insecta, Neuropterida) – a singular phenomenon or a model case in Raphidioptera? Deutsche Entomologische Zeitschrift 64(2): 123-131. https://doi.org/10.3897/dez.64.19859
|
The Mediterranean snakefly Raphidia mediterranea H. Aspöck, U. Aspöck & Rausch, 1977 – known from many parts of the Balkan Peninsula, several Aegean islands, southern parts of Italy, northwest of Anatolia and a few localities in Eastern Europe, yet not recorded in Central Europe – was surprisingly found with an astoundingly high population density on bushes in the yard of an old farmhouse at a comparatively high altitude (800 m) in Upper Austria, north of the Danube River, in 2013. This spectacular phenomenon was again observed in the following years (2014, 2015, 2016, 2017), and in 2016 the suspicion that the larvae develop in the straw of the thatched roof of the farmhouse could be confirmed by findings of larvae, pupae, and exuviae.It is most likely, that the occurrence of this Raphidia species in Austria is to be traced back to a human-caused introduction at some point in time. It remains, however, unknown when and specifically how this event might have occurred.
Morphologically no substantial differences were found between specimens from Greece, Italy and Upper Austria. In addition, the genetic uniformity (using 3 genes: cox1, cox3, and 28S) of the populations was verified. This supports an earlier hypothesis that the occurrence of the species, as well as that of R. mediterranea in Italy, Anatolia and perhaps elsewhere, may be related to importation of goods involving wood or soil.
A molecular genetic analysis of several Raphidia species confirmed the present morphology-based concept of their systematic position.
The means of dispersal of Raphidioptera are largely unknown. We do not know of any other similar cases of anthropogenic dispersal of a snakefly, but it cannot be excluded that human activities may have played a greater role in the dispersal of Raphidioptera than previously assumed. Phylogenomic studies would therefore be promising to solve some of these questions.
Raphidia mediterranea , man-made dispersal, neozoon, Raphidia , phylogeny
By the 1960s, it was revealed that the Balkan Peninsula was a distribution (and evolution) centre of the order Raphidioptera with an incredibly high number of species (
Numerous adults of both sexes of Raphidia mediterranea from many localities in Greece, Italy, Anatolia, as well as specimens of R. ophiopsis from Upper Austria and other parts of Central Europe were compared with imagines from Pelmberg (Upper Austria) based on the well-known morphological characters of male and female genitalia (
The distribution map was provided with ArcGis/ArcMap ver. 10.3.1.4959 based on the distribution records provided in the Suppl. material
For DNA analysis samples of four individuals of Raphidia mediterranea were selected, which had been collected in Pelmberg (Upper Austria), Gargano (Italy) and Zachlorou (Peloponnesus, Greece). Moreover, five representatives of the genus were included: Raphidia ophiopsis Linnaeus, 1758, Raphidia alcoholica H. Aspöck & U. Aspöck, 1969, Raphidia ulrikae H. Aspöck, 1964, Raphidia ariadne H. Aspöck & U. Aspöck, 1964, and Raphidia ligurica Albarda, 1891. A list of specimens analysed with exact localities is given in Table
Specimens analysed genetically in the present study and GenBank accession numbers for Cox1, Cox3 and 28S sequences.
Species | Sampling locality | Labcode | GenBank | ||
Cox1 | Cox3 | 28S | |||
Raphidia mediterranea | Austria, Upper Austria, Pelmberg near Hellmonsödt | Ramed-1 | MF975675 | MF975666 | MF975657 |
Raphidia mediterranea | Austria, Upper Austria, Pelmberg near Hellmonsödt | Ramed-2 | MF975676 | MF97567 | MF975658 |
Raphidia mediterranea | Greece, Peloponnesus, Kato Zachlorou | Ramed-3 | MF975677 | MF975668 | MF975659 |
Raphidia mediterranea | Italy, Puglia, Gargano, Punta Lunga | Ramed-4 | MF975678 | MF975669 | MF975660 |
Raphidia ariadne | Greece, Crete, W Omalos | Raari-1 | MF975679 | MF975672 | MF975661 |
Raphidia alcoholica | Greece, Phokis, S Pendayi | Raalc-1 | MF975680 | MF975671 | MF975662 |
Raphidia ophiopsis | Germany, Bayerischer Wald, Ruckwiesberg | Raoph-1 | MF975681 | MF975670 | MF975663 |
Raphidia ulrikae | Austria, Styria, Gulsen near Kraubath | Raulr-1 | MF975682 | MF975673 | MF975664 |
Raphidia ligurica | Italy, Sila Grande, near Viváio | Ralig-1 | MF975683 | MF975674 | MF975665 |
Two mitochondrial marker sequences were amplified using primers listed in Table
DNA extraction was performed using the DNeasy-Blood and Tissue Kit (QIAGEN, Hilden, Germany) according to the manufacturer’s instructions. The final volume of elution buffer was 40 µl. DNA solutions were stored in aliquots to avoid too frequent thawing. Control extractions with pure extraction buffer (without tissue) were prepared. PCR was carried out in an Eppendorf Thermocycler in a volume of 25 µl, containing Taq Polymerase (1.25 units/reaction; QIAGEN, Hilden, Germany) 1 µM of each primer, and 0.2 mM of each dNTP, 1.5 mM MgCl2, 5 µl Q-Solution, 2.5 µl 10x PCR buffer and 1 µl of template DNA. The PCR protocols were as follows: cox3: initial denaturation 94° C (3 min); 35 cycles: 94° C (60 s) / 50° C (30 s) / 72° C (60 sec); final extension at 72° C (10 min). cox1: initial denaturation 94° C (3 min); 35 cycles: 94° C (60 s) / 50° C (30 s) / 72° C (60 sec); final extension at 72° C (10 min). 28S: initial denaturation 94° C (3 min); 35 cycles: 94° C (60 s) / 55° C (30 s) / 72° C (60 sec); final extension at 72° C (10 min). Negative PCR controls were carried out to screen for contaminated reagents: (1) control extractions without tissue and subsequent PCR reactions (i.e. instead of template DNA); (2) PCR reactions with distilled water instead of template. PCR products were purified with the QIAquick PCR Purification Kit (QIAGEN, Hilden, Germany) prior to sequencing. Sequencing (both directions) was performed at Microsynth (Vienna, Austria) using the PCR primes as well as various internal primers (Table
Gene | Primer | Sequence (5´-3´) | Reference |
---|---|---|---|
Cox3 external primers | |||
Arth-cox3-fwd | 5’-TAGTTGATTATAGACCATGACC-3’ | a | |
Arth- cox3-rev | 5’-ACATCAACAAAATGTCAATATCA-3’ | a | |
Cox3 internal primers | |||
Raph- cox3-1+ | 5’-ACAAATTCCTTTATTAAATAC-3’ | present study | |
Raph- cox3-2- | 5’-CAWGTAATTGTTAATCCTGA-3’ | present study | |
Cox1 external primers | |||
Tyr-myr-1+ | 5’-CCCATAAATAAATTTACAGTTTA-3’ | present study | |
Leu-Myr-1- | 5’-GCACTATTCTGCCATATTAG-3’ | present study | |
Cox1 internal primers | |||
Raph- cox1-int1+ | 5’-TAGCAGGAGCTATCACTATACT-3’ | present study | |
Raph- cox1-int2- | 5’-ATATAAACTTCTGGATGTCC-3’ | present study | |
Raph- cox1-int5+ | 5’-CGAATACCTTTATTTGTATGATC-3’ | present study | |
Raph- cox1-int4- | 5’-AGAATAGGGTCTCCTCCTCC-3’ | present study | |
28S internal primers | |||
Raph-28S1+ | 5’-CAGGGGTAAACCTGAGAAA-3’ | b | |
Raph-28S-4- | 5’-AGCGCCAGTTCTGCTTACC-3’ | b | |
28S external primers | |||
Raph28S-3+ | 5’-AGCTTTGGGTACTTTCAGGA-3’ | b | |
Raph28S-2- | 5’-ACATGCTAGACTCCTTGGT-3’ | b |
Raw sequences were manually aligned in BioEdit v.7.1.3 (
Bayesian Inference (BI) was used for calculating phylogenetic trees. For BI the best fitting substitution model was determined for each of the three genes as well as codon positions of the protein coding genes by jModelTest v.2.1.5 (
In former studies (
The DNA sequence analysis revealed that the four specimens of R. mediterranea are identical in cox1 and 28S, while in cox3 a single substitution differentiating Ramed-4 from the other (identical) sequences was found. In general, the variation within 28S was extremely low. Except R. ligurica, which shows distances to the other ingroup taxa of 1.08 and 1.35%, respectively, sequences of all other ingroup species differ with p distances below 1% or are even identical. Between Raphidia and Agulla 28S distances ranged from 6.0 to 6.5%. Concerning the mitochondrial marker sequences, p distances between R. mediterranea and R. ophiopsis (the closest relative) were 5.08% in cox1and 5.62% in cox3, while the other species of Raphidia differed between 8.63–14.83% (cox1) and 7.46–16.94% (cox3) from R. mediterranea. Distances between species in the various gene sequences are illustrated by the NJ trees in Suppl. material
To assess the systematic position of R. mediterranea not only on the basis of morphological characters, we performed a phylogenetic analysis based on three genes (cox1, cox3, 28S). The two mt sequences resulted in trees in which the sister group of R. mediterranea is R. ophiopsis. In most analyses, R. alcoholica is the sister group of those two lineages, followed by R. ariadne; only in the BI tree of cox1 the relationships were unresolved (Suppl. material
The discovery of an isolated and unusually large population of Raphidia mediterranea – a Mediterranean species which has never been found elsewhere in Central Europe – in a farmhouse in a comparatively climatically unfavourable part of Upper Austria had raised the question concerning the origin of this population. It was assumed that morphological and/or genetic differences would be found, if the species had reached Upper Austria long ago by natural means of expansion of the distribution range. Therefore, specimens of the population from Upper Austria were compared with specimens from Greece and Italy. In the present study, the morphology-based results were clearly confirmed by molecular genetic analyses: The four specimens of R. mediterranea (two from Pelmberg (Austria), one from Greece, one from Italy) had almost identical sequences. It is legitimate to conclude that these populations were not separated long ago. The other species of Raphidia are clearly separated (see Fig.
Substantial differences could not be found in morphological characters, particularly in male and female genitalia, or in the sequences of three genes (cox1, cox3, 28S). This implies that all presently known and examined populations of R. mediterranea originated from a single glacial refugium. This refugium can reasonably be assumed to be in the south of the Balkan Peninsula as a part of the large balkanopontomediterranean refugium (
In Greece, the species occurs in many regions, in various habitats at altitudes of 10 – 1200 m and often in high population densities. Thus, it might have been occasionally transported to new habitats by ships carrying wood or soil. This could have occurred already in antiquity and throughout the past centuries.
The discovery of the isolated population of R. mediterranea in Upper Austria and the failure to detect any morphological or genetic differences between these vastly distant populations supports our previous hypothesis (
Concerning the systematic position of R. mediterranea, the phylogenetic analysis based on three genes confirmed our view of the systematics of R. ophiopsis, R. mediterranea, R. alcoholica, R. ariadne, R. ulrikae and R. ligurica (
We know little about the formation of distribution patterns of Raphidioptera. Snakeflies are generally characterised by low, in many cases extremely low, expansivity, and many species have hardly enlarged their distribution beyond their glacial refugial areas. In Central Europe, 16 species of Raphidioptera (13 Raphidiidae and 3 Inocelliidae) occur, 10 of these are of Mediterranean origin and presumably have reached Central Europe after the last glacial period, i.e. within the past 10,000 years (
Cordial thanks to Mag. Harald Bruckner (NHMW) for the photographs and for preparing the distribution map. The technical assistance of cand. med. Alexandra Szewczyk, BA, (Vienna) is gratefully acknowledged. We thank Prof. Dr. Sven Bacher (Fribourg) and Prof. Dr. Peter Duelli (Zürich) for their constructive reviews. We are very thankful to DDr. John Plant (Guilford, Connecticut) for polishing the English. We gratefully acknowledge the Museum für Naturkunde Berlin for waiving the author’s fees.
BI and NJ trees
Data type: molecular data
Explanation note: BI and NJ trees calculated with cox1, cox3 and 28S sequences, respectively. Posterior probability values (BI trees) and bootstrap values (in %, NJ trees) are indicated at the nodes.
Records of Raphidia mediterranea
Data type: Microsoft Excel Worksheet (.xlsx)
Explanation note: Records of Raphidia mediterranea documented in the distribution map (Fig.