Integrative analysis of DNA phylogeography and morphology of the European rose chafer (Cetonia aurata) to infer species taxonomy and patterns of postglacial colonisation in Europe

Species Diagnosis and DNA Taxonomy

The use of DNA has helped to improve and speed up species identification and delimitation. However, it also provides new challenges to taxonomists. Incongruence of outcome from various markers and delimitation methods, bias from sampling and skewed species distribution, implemented models, and the choice of methods/priors may mislead results and also may, in conclusion, increase elements of subjectivity in species taxonomy. The lack of direct diagnostic outcome from most contemporary molecular delimitation approaches and the need for a reference to existing and best sampled trait reference systems reveal the need for refining the criteria of species diagnosis and diagnosability in the current framework of nomenclature codes and good practices to avoid nomenclatorial instability, parallel taxonomies, and consequently more and new taxonomic impediment.

From DNA barcodes to ecology: Meta-analysis of central European beetles reveal link with species ecology but also to data pattern and gaps

DNA barcoding has been used worldwide to identify biological specimens and to delimit species. It represents a cost-effective, fast, and efficient way to assess biodiversity with help of the public Barcode of Life Database (BOLD) accounting for more than 236,000 animal species and more than 10 million barcode sequences. Here, we performed a meta-analysis of available barcode data of central European Coleoptera to detect intraspecific genetic patterns among ecological groups in relation to geographic distance with the aim to investigate a possible link between infraspecific variation and species ecology. We collected information regarding feeding style, body size, as well as habitat and biotope preferences. Mantel tests and two variants of Procrustes analysis, both involving the Principal Coordinates Neighborhood Matrices (PCNM) approach, were applied on genetic and geographic distance matrices. However, significance levels were too low to further use the outcome for further trait investigation: these were in mean for all ecological guilds only 7.5, 9.4, or 15.6% for PCNM + PCA, NMDS + PCA, and Mantel test, respectively, or at best 28% for a single guild. Our study confirmed that certain ecological traits were associated with higher species diversity and foster stronger genetic differentiation. Results suggest that increased numbers of species, sampling localities, and specimens for a chosen area of interest may give new insights to explore barcode data and species ecology for the scope of conservation on a larger scale. We performed a meta-analysis of available barcode data of central European beetles to detect intraspecific genetic patterns among ecological groups in relation to geographic distance, regarding feeding style, body size, as well as habitat and biotope preferences. Our study confirmed that certain ecological traits were associated with higher species diversity and foster stronger genetic differentiation. However, significance levels were too low to further use the outcome for further trait investigation.

Multiple species delimitation approaches with COI barcodes poorly fit each other and morphospecies – An integrative taxonomy case of Sri Lankan Sericini chafers (Coleoptera: Scarabaeidae)

DNA taxonomy including barcoding and metabarcoding is widely used to explore the diversity in biodiversity hotspots. In most of these hotspot areas, chafers are represented by a multitude of species, which are well defined by the complex shape of male genitalia. Here, we explore how well COI barcode data reflect morphological species entities and thus their usability for accelerated species inventorization. We conducted dedicated field surveys in Sri Lanka to collect the species‐rich and highly endemic Sericini chafers (Coleoptera: Scarabaeidae). Congruence among results of a series of protocols for de novo species delimitation and with morphology‐based species identifications was investigated. Different delimitation methods, such as the Poisson tree processes (PTP) model, Statistical Parsimony Analysis (TCS), Automatic Barcode Gap Discovery (ABGD), Assemble Species by Automatic Partitioning (ASAP), and Barcode Index Number (BIN) assignments, resulted in different numbers of molecular operational taxonomic units (MOTUs). All methods showed both over‐splitting and lumping of morphologically identified species. Only 18 of the observed 45 morphospecies perfectly matched MOTUs from all methods. The congruence of delimitation between MOTUs and morphospecies expressed by the match ratio was low, ranging from 0.57 to 0.67. TCS and multirate PTP (mPTP) showed the highest match ratio, while (BIN) assignment resulted in the lowest match ratio and most splitting events. mPTP lumped more species than any other method. Principal coordinate analysis (PCoA) on a match ratio‐based distance matrix revealed incongruent outcomes of multiple DNA delimitation methods, although applied to the same data. Our results confirm that COI barcode data alone are unlikely to correctly delimit all species, in particular, when using only a single delimitation approach. We encourage the integration of various approaches and data, particularly morphology, to validate species boundaries.

Excluding spatial sampling bias does not eliminate oversplitting in DNA-based species delimitation analyses

DNA barcoding and DNA-based species delimitation are major tools in DNA taxonomy. Sampling has been a central debate in this context, because the geographical composition of samples affects the accuracy and performance of DNA barcoding. Performance of complex DNA-based species delimitation is to be tested under simpler conditions in absence of geographic sampling bias. Here, we present an empirical dataset sampled from a single locality in a Southeast-Asian biodiversity hotspot (Laos: Phou Pan mountain). We investigate the performance of various species delimitation approaches on a megadiverse assemblage of herbivorous chafer beetles (Coleoptera: Scarabaeidae) to infer whether species delimitation suffers in the same way from exaggerate infraspecific variation despite the lack of geographic genetic variation that led to inconsistencies between entities from DNA-based and morphology-based species inference in previous studies. For this purpose, a 658 bp fragment of the mitochondrial cytochrome c oxidase subunit 1 (cox1) was analyzed for a total of 186 individuals of 56 morphospecies. Tree-based and distance-based species delimitation methods were used. All approaches showed a rather limited match ratio (max. 77%) with morphospecies. Poisson tree process (PTP) and statistical parsimony network analysis (TCS) prevailingly over-splitted morphospecies, while 3% clustering and Automatic Barcode Gap Discovery (ABGD) also lumped several species into one entity. ABGD revealed the highest congruence between molecular operational taxonomic units (MOTUs) and morphospecies. Disagreements between morphospecies and MOTUs have to be explained by historically acquired geographic genetic differentiation, incomplete lineage sorting, and hybridization. The study once again highlights how important morphology still is in order to correctly interpret the results of molecular species delimitation.

DNA Barcoding and geographical scale effect: The problems of undersampling genetic diversity hotspots

DNA barcoding identification needs a good characterization of intraspecific genetic divergence to establish the limits between species. Yet, the number of barcodes per species is many times low and geographically restricted. A poor coverage of the species distribution range may hamper identification, especially when undersampled areas host genetically distinct lineages. If so, the genetic distance between some query sequences and reference barcodes may exceed the maximum intraspecific threshold for unequivocal species assignation. Taking a group of Quercus herbivores (moths) in Europe as model system, we found that the number of DNA barcodes from southern Europe is proportionally very low in the Barcoding of Life Data Systems. This geographical bias complicates the identification of southern query sequences, due to their high intraspecific genetic distance with respect to barcodes from higher latitudes. Pairwise intraspecific genetic divergence increased along with spatial distance, but was higher when at least one of the sampling sites was in southern Europe. Accordingly, GMYC (General Mixed Yule Coalescent) single-threshold model retrieved clusters constituted exclusively by Iberian haplotypes, some of which could correspond to cryptic species. The number of putative species retrieved was more reliable than that of multiple-threshold GMYC but very similar to results from ABGD and jMOTU. Our results support GMYC as a key resource for species delimitation within poorly inventoried biogeographic regions in Europe, where historical factors (e.g., glaciations) have promoted genetic diversity and singularity. Future European DNA barcoding initiatives should be preferentially performed along latitudinal gradients, with special focus on southern peninsulas.

Hiding in plain sight: DNA barcoding suggests cryptic species in all ‘well-known’ Australian flower beetles (Scarabaeidae: Cetoniinae)

DNA barcode data is presented for Australian cetoniine flower beetles to aid with species discovery and guide revisionary taxonomy. Sequences of the COI gene’s DNA barcode region were acquired from 284 cetoniine specimens, covering 68 described species and 33 genera. This equates to 48% of the known species and 83% of the genera which occur in Australia. Results suggest up to 27 putative undescribed species in our sample, only 11 of which were suspected to be undescribed before this study, leaving 16 unexpected (“cryptic”) species. The Australian cetoniine fauna may hence be increased by up to 19%. An unanticipated result of the work is that each of the five most visible and commonly collected Australian cetoniine species, Eupoecila australasiae (Donovan, 1805), Neorrhina punctatum (Donovan, 1805), Glycyphana (Glycyphaniola) stolata (Fabricius, 1781), Chondropyga dorsalis (Donovan, 1805) and Bisallardiana gymnopleura (Fischer, 1823), have unexpectedly high diversity in DNA barcode sequences and were consequently split into multiple clusters, possibly indicating the presence of cryptic species.

Roots of the taxonomic impediment: Is the "integrativeness" a remedy?

The use of the "integrative approach" for classification of organisms since its formal establishment in 2005 has become a recurrent theme of zoosystematics. A bibliometric survey of the publications on integrative taxonomy of animals, which is aimed at exploring the most popular areas of research and characterizing the practical systematists' attitudes to this new approach, is presented. An analysis of 582 papers, which appeared between 2005 and 2017 in journals indexed by Scopus and the Web of Science Core Collection, has illustrated the gradual growth of the popularity of integrative taxonomy as well as some biases in the representation of higher taxa in "integrated" studies. It has been shown that the "integrative" papers have more chance of appearing in a top-ranking journal and gain relatively more citations as compared with non-integrative papers. The obtained results are discussed in the context of the "taxonomic impediment" problem thought to be a consequence of the institutional crisis of traditional taxonomy, which has been vividly debated over the past decades.

Genetic and morphological divergence among three closely related Phrynocephalus species (Agamidae)

Background

The Qinghai-Tibetan Plateau (QTP) is the world’s highest and largest plateau, but the role of its uplift in the evolution of species or biotas still remains poorly known. Toad-headed lizards of the reproductively bimodal genus Phrynocephalus are a clade of agamids, with all viviparous species restricted to the QTP and adjacent regions. The eastern part of the range of the viviparous taxa is occupied by three closely related but taxonomically controversial species, P. guinanensis, P. putjatia and P. vlangalii. Here, we combined genetic (mitochondrial ND4 gene and nine microsatellite loci), morphological (11 mensural and 11 meristic variables), and ecological (nine climatic variables) data to explore possible scenarios that may explain the discordance between genetic and morphological patterns, and to test whether morphological divergence is associated with local adaptation.

Results

We found weak genetic differentiation but pronounced morphological divergence, especially between P. guinanensis and P. vlangalii. Genetically, the species boundary was not so clear between any species pair. Morphologically, the species boundary was clear between P. guinanensis and P. vlangalii but not between other two species pairs. Body size and scale characters accounted best for morphological divergence between species. Morphological divergence was related to habitat types that differ climatically.

Conclusions

Our study provides evidence for genetic and morphological divergence among the three closely related viviparous species of Phrynocephalus lizards, and supports the idea that natural selection in spatially heterogeneous environments can lead to population divergence even in the presence of gene flow. Our study supports the hypothesis that the evolutionary divergence between viviparous Phrynocephalus species was a consequence of environmental change after the uplift of the QTP.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1443-y) contains supplementary material, which is available to authorized users.

Phylogeography and DNA-based species delimitation provide insight into the taxonomy of the polymorphic rose chafer Protaetia (Potosia) cuprea species complex (Coleoptera: Scarabaeidae: Cetoniinae) in the Western Palearctic

The development of modern methods of species delimitation, unified under the "integrated taxonomy" approach, allows a critical examination and re-evaluation of complex taxonomic groups. The rose chafer Protaetia (Potosia) cuprea is a highly polymorphic species group with a large distribution range. Despite its overall commonness, its taxonomy is unclear and subject to conflicting hypotheses, most of which largely fail to account for its evolutionary history. Based on the sequences of two mitochondrial markers from 65 individuals collected across the species range, and a detailed analysis of morphological characters including a geometric morphometry approach, we infer the evolutionary history and phylogeography of the P. cuprea species complex. Our results demonstrate the existence of three separate lineages in the Western Palearctic region, presumably with a species status. However, these lineages are in conflict with current taxonomic concepts. None of the 29 analyzed morphological characters commonly used in the taxonomy of this group proved to be unambiguously species- or subspecies- specific. The geometric morphometry analysis reveals a large overlap in the shape of the analyzed structures (pronotum, meso-metaventral projection, elytra and aedeagus), failing to identify either the genetically detected clades or the classical species entities. Our results question the monophyly of P. cuprea in regard to P. cuprina, as well as the species status of P. metallica. On the other hand, we found support for the species status of the Sicilian P. hypocrita. Collectively, our findings provide a new and original insight into the taxonomy and phylogeny of the P. cuprea species complex. At the same time, the results represent the first attempt to elucidate the phylogeography of these polymorphic beetles.

Local Hotspots of Endemism or Artifacts of Incorrect Taxonomy? The Status of Microendemic Pill Millipede Species of the Genus Glomeris in Northern Italy (Diplopoda, Glomerida)

Local endemic species with their unique evolutionary history always stirred the interest of scientists. One such area especially rich in endemics is northern Italy. In case of pill millipedes of the genus Glomeris Latreille, 1803, only a single species is found in northern Europe, while 22 country-endemics alone are known from Italy. Many of these endemics, however, have not been studied in several decades; therefore we aimed to determine whether this diversity is the result of overlooked synonymies or natural processes. A focus was placed on the local endemics that are in some aspects morphologically similar to the widespread and variable G. klugii Brandt, 1833. The local endemics Glomeris larii Verhoeff, 1921, G. primordialis Verhoeff, 1930, G. oblongoguttata Verhoeff, 1894, G. oropensis Verhoeff, 1936, G. transalpina Koch, 1836, G. romana Verhoeff, 1900, G. ligurica Latzel, 1884 and G. apuana Verhoeff, 1911 were included in a molecular analysis incorporating ribosomal nuclear (28S) and mitochondrial (COI) genes. Individuals were sequenced and compared to 31 specimens from 18 localities of G. klugii. The final dataset included 657 base pairs for 56 terminals in the COI, and 14 terminals with 1068 base pairs in the combined 28S and COI analysis. Our analysis shows intraspecific distances of up to 5% in the COI gene in G. klugii that are not strictly correlated to geography or color pattern. G. larii is discovered to be genetically and morphologically identical to G. klugii and is synonymised with the latter. Interspecific distances in our dataset vary between 6.7 to 15.9%, with the lowest (6.7-9.0%) between G. primordialis and G. klugii. Our analysis confirms the species status of the local endemics G. primordialis, G. oblongoguttata, G. oropensis, G. transalpina, G. ligurica and G. apuana. We also confirm the synonymy of G. undulata Koch, 1844 under G. klugii. G. genuensis Latzel, 1886 is indistinguishable from G. ligurica.

Novel members of the adipokinetic hormone family in beetles of the superfamily Scarabaeoidea

Eight beetle species of the superfamily Scarabaeoidea were investigated with respect to peptides belonging to the adipokinetic hormone (AKH) family in their neurohemal organs, the corpora cardiaca (CC). The following beetle families are represented: Scarabaeidae, Lucanidae, and Geotrupidae. AKH peptides were identified through a heterospecific trehalose-mobilizing bioassay and by sequence analyses, using liquid chromatography coupled to positive electrospray mass spectrometry (LC-ESI-MS) and analysis of the tandem MS² spectra obtained by collision-induced dissociation. All the beetle species have octapeptide AKHs; some have two AKHs, while others have only one. Novel AKH members were found in Euoniticellus intermedius and Circellium bacchus (family Scarabaeidae), as well as in Dorcus parallelipipedus (family Lucanidae). Two species of the family Geotrupidae and two species of the Scarabaeidae subfamily Cetoniinae contain one known AKH peptide, Melme-CC, while E. intermedius produces a novel peptide code named Euoin-AKH: pEINFTTGWamide. Two AKH peptides were each identified in CC of C. bacchus and D. parallelipipedus: the novel Cirba-AKH: pEFNFSAGWamide and the known peptide, Scade-CC-I in the former, and the novel Dorpa-AKH: pEVNYSPVW amide and the known peptide, Melme-CC in the latter. Kheper bonelli (subfamily Scarabaeinae) also has two AKHs, the known Scade-CC-I and Scade-CC-II. All the novel peptides were synthesized and the amino acid sequence assignments were unequivocally confirmed by co-elution of the synthetic peptides with their natural equivalent, and identical MS parameters of the two forms. The novel synthetic peptides are all active in inducing hypertrehalosemia in cockroaches.

Molecular phylogeny of the genus Dicronocephalus (Coleoptera, Scarabaeidae, Cetoniinae) based on mtCOI and 16S rRNA genes

The seven species belonging to the genus Dicronocephalus are a very interesting group with a unique appearance and distinct sexual dimorphism. Only one species among them, Dicronocephalusadamsi, has been known in the Korean fauna. This species is recognized as having a wide distribution from Tibet to Korean Peninsula and is currently represented by two subspecies that have separated geographical ranges. The phylogenetic relationships of Dicronocephalusadamsi were still unclear. The phylogeny of Dicronocephalus is reconstructed with a phylogenetic study of five species including four subspecies based on a molecular approach using mitochondrial COI and 16S rRNA genes. Our results are compared with the results obtained by previous authors based on morphological characters. They show that the tested taxa are divided into two major clades. Clade A consists of two species (Dicronocephalusadamsi + Dicranocephalusyui) and Clade B includes the others (Dicronocephalusdabryi + Dicranocephalusuenoi + Dicranocephaluswallichii). This result generally supports Kurosawa's proposal except that Dicronocephalusdabryi and Dicranocephalusuenoi are newly recognized as members of a monophyletic group. We propose that Dicronocephalusadamsidrumonti is a junior subjective synonym of Dicronocephalusadamsiadamsi. These results show that three members of the Dicranocephaluswallichii group should be treated as species rather than subspecies. However, further research including analyses of different genetic markers is needed to reconfirm our results.

Streamlining DNA barcoding protocols: automated DNA extraction and a new cox1 primer in arachnid systematics

Background

DNA barcoding is a popular tool in taxonomic and phylogenetic studies, but for most animal lineages protocols for obtaining the barcoding sequences—mitochondrial cytochrome C oxidase subunit I (cox1 AKA CO1)—are not standardized. Our aim was to explore an optimal strategy for arachnids, focusing on the species-richest lineage, spiders by (1) improving an automated DNA extraction protocol, (2) testing the performance of commonly used primer combinations, and (3) developing a new cox1 primer suitable for more efficient alignment and phylogenetic analyses.

Methodology

We used exemplars of 15 species from all major spider clades, processed a range of spider tissues of varying size and quality, optimized genomic DNA extraction using the MagMAX Express magnetic particle processor—an automated high throughput DNA extraction system—and tested cox1 amplification protocols emphasizing the standard barcoding region using ten routinely employed primer pairs.

Results

The best results were obtained with the commonly used Folmer primers (LCO1490/HCO2198) that capture the standard barcode region, and with the C1-J-2183/C1-N-2776 primer pair that amplifies its extension. However, C1-J-2183 is designed too close to HCO2198 for well-interpreted, continuous sequence data, and in practice the resulting sequences from the two primer pairs rarely overlap. We therefore designed a new forward primer C1-J-2123 60 base pairs upstream of the C1-J-2183 binding site. The success rate of this new primer (93%) matched that of C1-J-2183.

Conclusions

The use of C1-J-2123 allows full, indel-free overlap of sequences obtained with the standard Folmer primers and with C1-J-2123 primer pair. Our preliminary tests suggest that in addition to spiders, C1-J-2123 will also perform in other arachnids and several other invertebrates. We provide optimal PCR protocols for these primer sets, and recommend using them for systematic efforts beyond DNA barcoding.