Ignacio Ribera

A Comprehensive Phylogeny of Beetles Reveals the Evolutionary Origins of a Superradiation

Beetles represent almost one-fourth of all described species, and knowledge about their relationships and evolution adds to our understanding of biodiversity. We performed a comprehensive phylogenetic analysis of Coleoptera inferred from three genes and nearly 1900 species, representing more than 80% of the worlds recognized beetle families. We defined basal relationships in the Polyphaga supergroup, which contains over 300,000 species, and established five families as the earliest branching lineages. By dating the phylogeny, we found that the success of beetles is explained neither by exceptional net diversification rates nor by a predominant role of herbivory and the Cretaceous rise of angiosperms. Instead, the pre-Cretaceous origin of more than 100 present-day lineages suggests that beetle species richness is due to high survival of lineages and sustained diversification in a variety of niches.

The Effect of Geographical Scale of Sampling on DNA Barcoding

Abstract Eight years after DNA barcoding was formally proposed on a large scale, CO1 sequences are rapidly accumulating from around the world. While studies to date have mostly targeted local or regional species assemblages, the recent launch of the global iBOL project (International Barcode of Life), highlights the need to understand the effects of geographical scale on Barcodings goals. Sampling has been central in the debate on DNA Barcoding, but the effect of the geographical scale of sampling has not yet been thoroughly and explicitly tested with empirical data. Here, we present a CO1 data set of aquatic predaceous diving beetles of the tribe Agabini, sampled throughout Europe, and use it to investigate how the geographic scale of sampling affects 1) the estimated intraspecific variation of species, 2) the genetic distance to the most closely related heterospecific, 3) the ratio of intraspecific and interspecific variation, 4) the frequency of taxonomically recognized species found to be monophyletic, and 5) query identification performance based on 6 different species assignment methods. Intraspecific variation was significantly correlated with the geographical scale of sampling (R-square = 0.7), and more than half of the species with 10 or more sampled individuals (N = 29) showed higher intraspecific variation than 1% sequence divergence. In contrast, the distance to the closest heterospecific showed a significant decrease with increasing geographical scale of sampling. The average genetic distance dropped from > 7% for samples within 1 km, to < 3.5% for samples up to > 6000 km apart. Over a third of the species were not monophyletic, and the proportion increased through locally, nationally, regionally, and continentally restricted subsets of the data. The success of identifying queries decreased with increasing spatial scale of sampling; liberal methods declined from 100% to around 90%, whereas strict methods dropped to below 50% at continental scales. The proportion of query identifications considered uncertain (more than one species < 1% distance from query) escalated from zero at local, to 50% at continental scale. Finally, by resampling the most widely sampled species we show that even if samples are collected to maximize the geographical coverage, up to 70 individuals are required to sample 95% of intraspecific variation. The results show that the geographical scale of sampling has a critical impact on the global application of DNA barcoding. Scale-effects result from the relative importance of different processes determining the composition of regional species assemblages (dispersal and ecological assembly) and global clades (demography, speciation, and extinction). The incorporation of geographical information, where available, will be required to obtain identification rates at global scales equivalent to those in regional barcoding studies. Our result hence provides an impetus for both smarter barcoding tools and sprouting national barcoding initiatives—smaller geographical scales deliver higher accuracy.

Speciation of Iberian diving beetles in Pleistocene refugia (Coleoptera, Dytiscidae)

The Mediterranean basin is an area of high diversity and endemicity, but the age and origin of its fauna are still largely unknown. Here we use species‐level phylogenies based on ≈ 1300 base pairs of the genes 16S rRNA and cytochrome oxidase I to establish the relationships of 27 of the 34 endemic Iberian species of diving beetles in the family Dytiscidae, and to investigate their level of divergence. Using a molecular clock approach, 18–19 of these species were estimated to be of Pleistocene origin, with four to six of them from the Late Pleistocene (≈ 100 000 years). A second, lower speciation frequency peak was assigned to Late Miocene or Early Pliocene. Analysis of the distributional ranges showed that endemic species placed in the tip nodes of the trees are significantly more likely to be allopatric with their sisters than endemic species at lower node levels. Allopatric sister species are also significantly younger than sympatric clades, in agreement with an allopatric mode of speciation and limited subsequent range movement. These results strongly suggest that for some taxa Iberian populations were isolated during the Pleistocene long enough to speciate, and apparently did not expand their ranges to recolonize areas north of the Pyrenees. This is in contradiction to observations from fossil beetles in areas further north, which document large range movements associated with the Pleistocene glacial cycles hypothesized to suppress population isolation and allopatric speciation.

Biogeography and conservation of Iberian water beetles

An analysis is made of the main biogeographical trends and the conservation status of the aquatic Coleoptera from the Iberian peninsula. More than 600 species in 100 genera and 20 families are recognised, although knowledge of some of them is still scarce. Overall, the level of endemism is about 20%, with a highly uneven taxonomic distribution. Five main geographical regions are recognised based on geological and topographical criteria: Pyrenean, Cantabrian mountains, Hercynian, South East, and South West. These regions have a significantly distinct fauna when all species, species of the family Hydraenidae, all Iberian endemics, and Hydraenidae endemics are considered. For other groups this geographical division is not relevant, as the species follow a random distribution. The main biogeographical division is between the northern and the southern regions, with closer similarity between the Pyrenees and France, and between the South East and Morocco, than with other Iberian regions. The conservation status of the species included in the IUCN red list is reviewed, suggesting the exclusion of three species previously considered to be endangered or vulnerable, and the inclusion of two new species as vulnerable. The status of the rare species is discussed, with the distinction between local endemics, possible relict species, and species with discontinuous distributions due to habitat loss or habitat fragmentation. The habitats with the highest conservation value are medium altitude freshwater streams in the Hercynian mountains of central and western Iberia, the pre-Pyrenees and the Baetic cordilleras, together with saline streams in the South East and the central Ebro Valley. Some general guidelines are suggested for their preservation.

Mitochondrial Cox1 Sequence Data Reliably Uncover Patterns of Insect Diversity But Suffer from High Lineage-Idiosyncratic Error Rates

Background The demand for scientific biodiversity data is increasing, but taxonomic expertise is often limited or not available. DNA sequencing is a potential remedy to overcome this taxonomic impediment. Mitochondrial DNA is most commonly used, e.g., for species identification (“DNA barcoding”). Here, we present the first study in arthropods based on a near-complete species sampling of a family-level taxon from the entire Australian region. We aimed to assess how reliably mtDNA data can capture species diversity when many sister species pairs are included. Then, we contrasted phylogenetic subsampling with the hitherto more commonly applied geographical subsampling, where sister species are not necessarily captured. Methodology/Principal Findings We sequenced 800 bp cox1 for 1,439 individuals including 260 Australian species (78% species coverage). We used clustering with thresholds of 1 to 10% and general mixed Yule Coalescent (GMYC) analysis for the estimation of species richness. The performance metrics used were taxonomic accuracy and agreement between the morphological and molecular species richness estimation. Clustering (at the 3% level) and GMYC reliably estimated species diversity for single or multiple geographic regions, with an error for larger clades of lower than 10%, thus outperforming parataxonomy. However, the rates of error were higher for some individual genera, with values of up to 45% when very recent species formed nonmonophyletic clusters. Taxonomic accuracy was always lower, with error rates above 20% and a larger variation at the genus level (0 to 70%). Sørensen similarity indices calculated for morphospecies, 3% clusters and GMYC entities for different pairs of localities was consistent among methods and showed expected decrease over distance. Conclusion/Significance Cox1 sequence data are a powerful tool for large-scale species richness estimation, with a great potential for use in ecology and β-diversity studies and for setting conservation priorities. However, error rates can be high in individual lineages.

Discovery of Aspidytidae, a new family of aquatic Coleoptera.

The six extant aquatic families of Hydradephaga (Coleoptera) known so far represent a diverse group of beetles morphologically highly modified for life in the water. We report the discovery of a new genus with two species from South Africa and China, which differ greatly from all extant families, but resemble the Jurassic–Cretaceous †Liadytidae (the dagger symbol indicates that the taxa are known only as fossils). Based on a combined phylogenetic analysis of molecular and morphological data we erect a new family, Aspidytidae, which is the sister group of Dytiscidae plus Hygrobiidae. We propose a new scenario for the evolution of swimming behaviour in adephagan beetles, in which the transition into the aquatic environment is followed by complex and repeated changes in lifestyles, including the secondary complete loss of swimming ability in Aspidytidae.

Molecular phylogeny of Pacific Island Colymbetinae: radiation of New Caledonian and Fijian species (Coleoptera, Dytiscidae)

We present a molecular phylogeny and taxonomic review of the Pacific island colymbetine diving beetles, focusing on the Fijian and New Caledonian faunas. Four new species are described: Rhantus monteithi and R. poellerbauerae from New Caledonia, and R. kini and R. bula from Fiji. We also describe the 3rd instar larvae of R. monteithi and R. poellerbauerae spp. nov., assigned to adults using mtDNA sequence data and discuss larval characters in the light of phylogeny. The phylogenetic hypotheses derived from both parsimony and Bayesian inference based on 3508 aligned nucleotides from a combination of mitochondrial (cox1, cob and rrnL‐tRNALeu‐Nad1) and nuclear genes (18S rRNA and H3) reveal a clade comprising R. novaecaledoniae, R. alutaceus, R. pseudopacificus, R. monteithi sp. nov. and R. poellerbauerae sp. nov., which agrees with the R. pacificus group sensu Balke (1993). Carabdytes upin was included within this clade, possibly indicating paraphyly of the genus Rhantus. Rhantus annectens, R. bacchusi, R. supranubicus, R. suturalis, R. simulans, and the Palearctic R. exsoletus, R. latitans and R. bistriatus formed a clade corresponding to the R. suturalis group sensu Balke (2001). Rhantus vitiensis, previously assigned to the R. pacificus group, was included in the R. suturalis clade. We find some support for a scenario where the Pacific was colonized out of the Northern hemisphere only during the past c. 12 million years, rejecting a Gondwanan origin of the morphologically isolated endemics. The new species are all characterized by mtDNA haplotype clusters, the degree of divergence between sister species pairs ranging from 1.3 to 7%, while R. novaecaledoniae individuals from all over New Caledonia apparently form one morphospecies, with moderate genetic diversity (up to 2.3% mtDNA divergence between populations). The sisters R. pollerbauerae sp. nov. + R. monteithi sp. nov. occur sympatrically on Mont Panié but appear ecologically segregated, while the sisters R. vitiensis + R. bula sp. nov. were encountered syntopically on Viti Levu. Comparing genetic and morphological data of Fijian Rhantus and Copelatus diving beetles, we here show that even in island radiations it is not per se possible to know if mitochondrial DNA barcoding would perform well (Rhantus: YES, Copelatus: NO). At the same time we show that fixed cutoff values, as sometimes used to discriminate between barcodes, thus species, might be meaningless. We underpin the importance of morphology for sustainable exploration of global diversity.

Mitochondrial DNA phylogeography and population history of Meladema diving beetles on the Atlantic Islands and in the Mediterranean basin (Coleoptera, Dytiscidae)

The phylogeny and population history of Meladema diving beetles (Coleoptera, Dytiscidae) were examined using mitochondrial DNA sequence from 16S ribosomal RNA and cytochrome oxidase I genes in 51 individuals from 22 populations of the three extant species (M. imbricata endemic to the western Canary Islands, M. lanio endemic to Madeira and M. coriacea widespread in the Western Mediterranean and on the western Canaries), using a combination of phylogenetic and nested clade analyses. Four main lineages are observed within Meladema, representing the three recognized species plus Corsican populations of M. coriacea. Phylogenetic analyses demonstrate the sister relationship of the two Atlantic Island taxa, and suggest the possible paraphyly of M. coriacea. A molecular clock approach reveals that speciation within the genus occurred in the Early Pleistocene, indicating that the Atlantic Island endemics are not Tertiary relict taxa as had been proposed previously. Our results point to past population bottlenecks in all four lineages, with recent (Late‐Middle Pleistocene) range expansion in non‐Corsican M. coriacea and M. imbricata. Within the Canary Islands, M. imbricata seems to have independently colonized La Gomera and La Palma from Tenerife (although a colonization of La Palma from La Gomera cannot be discarded), and M. coriacea has independently colonized Tenerife and Gran Canaria from separate mainland lineages. In the Mediterranean basin this species apparently colonized Corsica on a single occasion, relatively early in its evolutionary history (Early Pleistocene), and has colonized Mallorca recently on multiple occasions. On the only island where M. coriacea and M. imbricata are broadly sympatric (Tenerife), we report evidence of bidirectional hybridization between the two species.

Phylogeny and historical biogeography of Agabinae diving beetles (Coleoptera) inferred from mitochondrial DNA sequences

The Agabinae, with more than 350 species, is one of the most diverse lineages of diving beetles (Dytiscidae). Using the mitochondrial genes 16S rRNA and cytochrome oxidase I we present a phylogenetic analysis based on 107 species drawn mostly from the four main Holarctic genera. Two of these genera (Ilybius and Ilybiosoma) are consistently recovered as monophyletic with strong support, Platambus is never recovered as monophyletic, and Agabus is found paraphyletic with respect to several of the species groups of Platambus. Basal relationships among the main lineages are poorly defined, although within each of them relationships are in general robust and very consistent across the parameter space, and in agreement with previous morphological analyses. In each of the two most diverse lineages (Ilybius and Agabus including part of Platambus) there is a basal split between Palearctic and Nearctic clades, estimated to have occurred in the late Eocene. The Palearctic clade in turn splits into a Western Palearctic clade and a clade containing mostly Eastern Palearctic species, and assumed to be ancestrally Eastern Palearctic but with numerous transitions to a Holarctic or Nearctic distribution. These results suggest an asymmetry in the colonization routes, as there are very few cases of transcontinental range expansions originating from the Nearctic or the Western Palearctic. According to standard clock estimates, we do not find any transcontinental shift during the Pliocene, but numerous speciation events within each of the continental or subcontinental regions.

Agabus alexandrae sp. n. from Morocco, with a molecular phylogeny of the Western Mediterranean species of the A. guttatus group (Coleoptera: Dytiscidae)

A new species of the Agabus guttatus group, A. alexandrae sp. n., is described from the Moyen Atlas in Morocco. A phylogeny of most of the Western Mediterranean species of the group based on Cytochrome Oxydase I and 16S rRNA sequences (approximately 1,250 bp) shows that A. alexandrae sp. n. is sister to the remaining species of the group. The evolutionary history of the species of the group is tentatively reconstructed, based on a Maximum Likelihood estimate of divergence time of approximately 2% per MY for the combined 16S and COI genes, similar to other estimates of evolutionary rates of the mitochondrial DNA of Coleoptera.