Carmelo Andújar

Winding up the molecular clock in the genus Carabus (Coleoptera: Carabidae): assessment of methodological decisions on rate and node age estimation.

BackgroundRates of molecular evolution are known to vary across taxa and among genes, and this requires rate calibration for each specific dataset based on external information. Calibration is sensitive to evolutionary model parameters, partitioning schemes and clock model. However, the way in which these and other analytical aspects affect both the rates and the resulting clade ages from calibrated phylogenies are not yet well understood. To investigate these aspects we have conducted calibration analyses for the genus Carabus (Coleoptera, Carabidae) on five mitochondrial and four nuclear DNA fragments with 7888 nt total length, testing different clock models and partitioning schemes to select the most suitable using Bayes Factors comparisons.ResultsWe used these data to investigate the effect of ambiguous character and outgroup inclusion on both the rates of molecular evolution and the TMRCA of Carabus. We found considerable variation in rates of molecular evolution depending on the fragment studied (ranging from 5.02% in cob to 0.26% divergence/My in LSU-A), but also on analytical conditions. Alternative choices of clock model, partitioning scheme, treatment of ambiguous characters, and outgroup inclusion resulted in rate increments ranging from 28% (HUWE1) to 1000% (LSU-B and ITS2) and increments in the TMRCA of Carabus ranging from 8.4% (cox1-A) to 540% (ITS2). Results support an origin of the genus Carabus during the Oligocene in the Eurasian continent followed by a Miocene differentiation that originated all main extant lineages.ConclusionsThe combination of several genes is proposed as the best strategy to minimise both the idiosyncratic behaviors of individual markers and the effect of analytical aspects in rate and age estimations. Our results highlight the importance of estimating rates of molecular evolution for each specific dataset, selecting for optimal clock and partitioning models as well as other methodological issues potentially affecting rate estimation.

Phylogenetic community ecology of soil biodiversity using mitochondrial metagenomics

High‐throughput DNA methods hold great promise for the study of taxonomically intractable mesofauna of the soil. Here, we assess species diversity and community structure in a phylogenetic framework, by sequencing total DNA from bulk specimen samples and assembly of mitochondrial genomes. The combination of mitochondrial metagenomics and DNA barcode sequencing of 1494 specimens in 69 soil samples from three geographic regions in southern Iberia revealed >300 species of soil Coleoptera (beetles) from a broad spectrum of phylogenetic lineages. A set of 214 mitochondrial sequences longer than 3000 bp was generated and used to estimate a well‐supported phylogenetic tree of the order Coleoptera. Shorter sequences, including cox1 barcodes, were placed on this mitogenomic tree. Raw Illumina reads were mapped against all available sequences to test for species present in local samples. This approach simultaneously established the species richness, phylogenetic composition and community turnover at species and phylogenetic levels. We find a strong signature of vertical structuring in soil fauna that shows high local community differentiation between deep soil and superficial horizons at phylogenetic levels. Within the two vertical layers, turnover among regions was primarily at the tip (species) level and was stronger in the deep soil than leaf litter communities, pointing to layer‐mediated drivers determining species diversification, spatial structure and evolutionary assembly of soil communities. This integrated phylogenetic framework opens the application of phylogenetic community ecology to the mesofauna of the soil, among the most diverse and least well‐understood ecosystems, and will propel both theoretical and applied soil science.

Integrative taxonomy and conservation of cryptic beetles in the Mediterranean region (Hydrophilidae)

Arribas, P., Andújar, C., Sánchez‐Fernández, D., Abellán, P. & Millán, A. (2012). Integrative taxonomy and conservation of cryptic beetles in the Mediterranean region (Hydrophilidae). —Zoologica Scripta, 00, 000–000.

Tempo and mode of the multiple origins of salinity tolerance in a water beetle lineage

Salinity is one of the most important drivers of the distribution, abundance and diversity of organisms. Previous studies on the evolution of saline tolerance have been mainly centred on marine and terrestrial organisms, while lineages inhabiting inland waters remain largely unexplored. This is despite the fact that these systems include a much broader range of salinities, going from freshwater to more than six times the salinity of the sea (i.e. >200 g/L). Here, we study the pattern and timing of the evolution of the tolerance to salinity in an inland aquatic lineage of water beetles (Enochrus species of the subgenus Lumetus, family Hydrophilidae), with the general aim of understanding the mechanisms by which it was achieved. Using a time‐calibrated phylogeny built from five mitochondrial and two nuclear genes and information about the salinity tolerance and geographical distribution of the species, we found that salinity tolerance appeared multiple times associated with periods of global aridification. We found evidence of some accelerated transitions from freshwater directly to high salinities, as reconstructed with extant lineages. This, together with the strong positive correlation found between salinity tolerance and aridity of the habitats in which species are found, suggests that tolerance to salinity may be based on a co‐opted mechanism developed originally for drought resistance.

Integration of conflict into integrative taxonomy: fitting hybridization in species delimitation of Mesocarabus (Coleoptera: Carabidae)

In species differentiation, characters may not diverge synchronously, and there are also processes that shuffle character states in lineages descendant from a common ancestor. Species are thus expected to show some degree of incongruence among characters; therefore, taxonomic delimitation can benefit from integrative approaches and objective strategies that account for character conflict. We illustrate the potential of exploiting conflict for species delimitation in a study case of ground beetles of the subgenus Carabus (Mesocarabus), where traditional taxonomy does not accurately delimit species. The molecular phylogenies of four mitochondrial and three nuclear genes, cladistic analysis of the aedeagus, ecological niche divergence and morphometry of pronotal shape in more than 500 specimens of Mesocarabus show that these characters are not fully congruent. For these data, a three‐step operational strategy is proposed for species delimitation by (i) delineating candidate species based on the integration of incongruence among conclusive lines of evidence, (ii) corroborating candidate species with inconclusive lines of evidence and (iii) refining a final species proposal based on an integrated characterization of candidate species based on the evolutionary analysis of incongruence. This procedure provided a general understanding of the reticulate process of hybridization and introgression acting on Mesocarabus and generated the hypothesis of seven Mesocarabus species, including two putative hybrid lineages. Our work emphasizes the importance of incorporating critical analyses of character and phylogenetic conflict to infer both the evolutionary history and species boundaries through an integrative taxonomic approach.

Congruence test of molecular clock calibration hypotheses based on Bayes factor comparisons

This work was supported by the Spanish Ministry of Science and Innovation (CGL2006/06706 and CGL2009-10906 to C.A. and J.S. and CGL2008-00007 to J.G.-Z., the latter also with support of the European Regional Development Fund). CA received support from an FPU predoctoral studentship (Spanish Ministry of Education).

The mitochondrial genome of Iberobaenia (Coleoptera: Iberobaeniidae): first rearrangement of protein-coding genes in the beetles

Abstract The complete mitochondrial genome of the recently discovered beetle family Iberobaeniidae is described and compared with known coleopteran mitogenomes. The mitochondrial sequence was obtained by shotgun metagenomic sequencing using the Illumina Miseq technology and resulted in an average coverage of 130 × and a minimum coverage of 35×. The mitochondrial genome of Iberobaeniidae includes 13 protein-coding genes, 2 rRNAs, 22 tRNAs genes, and 1 putative control region, and showed a unique rearrangement of protein-coding genes. This is the first rearrangement affecting the relative position of protein-coding and ribosomal genes reported for the order Coleoptera.

Metabarcoding of freshwater invertebrates to detect the effects of a pesticide spill

Biomonitoring underpins the environmental assessment of freshwater ecosystems and guides management and conservation. Current methodology for surveys of (macro)invertebrates uses coarse taxonomic identification where species‐level resolution is difficult to obtain. Next‐generation sequencing of entire assemblages (metabarcoding) provides a new approach for species detection, but requires further validation. We used metabarcoding of invertebrate assemblages with two fragments of the cox1 “barcode” and partial nuclear ribosomal (SSU) genes, to assess the effects of a pesticide spill in the River Kennet (southern England). Operational taxonomic unit (OTU) recovery was tested under 72 parameters (read denoising, filtering, pair merging and clustering). Similar taxonomic profiles were obtained under a broad range of parameters. The SSU marker recovered Platyhelminthes and Nematoda, missed by cox1, while Rotifera were only amplified with cox1. A reference set was created from all available barcode entries for Arthropoda in the BOLD database and clustered into OTUs. The River Kennet metabarcoding produced matches to 207 of these reference OTUs, five times the number of species recognized with morphological monitoring. The increase was due to the following: greater taxonomic resolution (e.g., splitting a single morphotaxon “Chironomidae” into 55 named OTUs); splitting of Linnaean binomials into multiple molecular OTUs; and the use of a filtration‐flotation protocol for extraction of minute specimens (meiofauna). Community analyses revealed strong differences between “impacted” vs. “control” samples, detectable with each gene marker, for each major taxonomic group, and for meio‐ and macrofaunal samples separately. Thus, highly resolved taxonomic data can be extracted at a fraction of the time and cost of traditional nonmolecular methods, opening new avenues for freshwater invertebrate biodiversity monitoring and molecular ecology.

Gondwanian relicts and oceanic dispersal in a cosmopolitan radiation of euedaphic ground beetles

Anillini are a tribe of minute, euedaphic ground beetles (Carabidae) characterized by the loss of eyes, loss of wings and high levels of local endemism. Despite their presumed low dispersal, they have a nearly cosmopolitan distribution, including isolated islands such as New Zealand and New Caledonia. We used a time calibrated molecular phylogeny to test, first, if the tribe as currently understood is monophyletic and, second, whether the time of divergence is compatible with an early vicariant diversification after the breakup of Gondwana. We sequenced portions of 6 mitochondrial and 3 nuclear genes for 66 specimens in 17 genera of Anillini plus 39 outgroups. The resulting phylogenetic tree was used to estimate the time of diversification using two independent calibration schemes, by applying molecular rates for the related genus Carabus or by dating the tree with fossil and geological information. Rates of molecular evolution and lineage ages were mostly concordant between both calibration schemes. The monophyly of Anillini was well-supported, and its age was consistent with a Gondwanian origin of the main lineages and an initial diversification at ca. 100Ma representing the split between the eyed Nesamblyops (New Zealand) and the remaining Anillini. The subsequent diversification, including the split of the Nearctic Anillinus and the subsequent splits of Palaearctic lineages, was dated to between 80 and 100Ma and thus was also compatible with a tectonic vicariant origin. On the contrary, the estimated age of the New Caledonian blind Orthotyphlus at ca. 30±20Ma was incompatible with a vicariant origin, suggesting the possibility of trans-oceanic dispersal in these endogean beetles.

Lessons from genome skimming of arthropod-preserving ethanol

Field‐collected specimens of invertebrates are regularly killed and preserved in ethanol, prior to DNA extraction from the specimens, while the ethanol fraction is usually discarded. However, DNA may be released from the specimens into the ethanol, which can potentially be exploited to study species diversity in the sample without the need for DNA extraction from tissue. We used shallow shotgun sequencing of the total DNA to characterize the preservative ethanol from two pools of insects (from a freshwater habitat and terrestrial habitat) to evaluate the efficiency of DNA transfer from the specimens to the ethanol. In parallel, the specimens themselves were subjected to bulk DNA extraction and shotgun sequencing, followed by assembly of mitochondrial genomes for 39 of 40 species in the two pools. Shotgun sequencing from the ethanol fraction and read‐matching to the mitogenomes detected ~40% of the arthropod species in the ethanol, confirming the transfer of DNA whose quantity was correlated to the biomass of specimens. The comparison of diversity profiles of microbiota in specimen and ethanol samples showed that ‘closed association’ (internal tissue) bacterial species tend to be more abundant in DNA extracted from the specimens, while ‘open association’ symbionts were enriched in the preservative fluid. The vomiting reflex of many insects also ensures that gut content is released into the ethanol, which provides easy access to DNA from prey items. Shotgun sequencing of DNA from preservative ethanol provides novel opportunities for characterizing the functional or ecological components of an ecosystem and their trophic interactions.