Paula Arribas

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.

Reduced salinities compromise the thermal tolerance of hypersaline specialist diving beetles.

Inland saline waters are globally threatened habitats that harbour unique assemblages of specialist invertebrates. In many Mediterranean regions, irrigation associated with intensive agriculture is lowering the salinity of these habitats, resulting in the loss of their specialist biota, although the mechanisms by which reductions in salinity lead to species loss are poorly understood. In the present study, the effects of reduced salinity on the temperature tolerance and thermal acclimatory abilities of two related species of hypersaline water beetles, Nebrioporus baeticus (Schaum) and Nebrioporus ceresyi (Aubé), are explored. Both upper (UTL) and lower thermal limits (LTL) are assessed, and both salinity and temperature are found to influence the thermal biology of Nebrioporus. Mean UTLs are greater in individuals of both species acclimated at high salinities, with salinity appearing to be more important than acclimation temperature in determining UTL. In both taxa, the lowest mean LTLs are recorded in individuals acclimated at the highest salinities and lowest temperatures; temperature‐dependent acclimation is only reported after exposure to relatively high salinities. The data show that salinity influences the thermal tolerance and acclimatory ability of these hypersaline beetles, and that lowered salinity compromises the ability of adult Nebrioporus to cope with both heat and cold. Such an effect may partly explain why specialist species are lost from hypersaline habitats subject to salinity reductions, and suggests that ongoing reduction in salinity may compromise the ability of such specialist taxa to cope with rapid climate change.

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.

The comparative osmoregulatory ability of two water beetle genera whose species span the fresh-hypersaline gradient in inland waters (Coleoptera: Dytiscidae, Hydrophilidae).

A better knowledge of the physiological basis of salinity tolerance is essential to understanding the ecology and evolutionary history of organisms that have colonized inland saline waters. Coleoptera are amongst the most diverse macroinvertebrates in inland waters, including saline habitats; however, the osmoregulatory strategies they employ to deal with osmotic stress remain unexplored. Survival and haemolymph osmotic concentration at different salinities were examined in adults of eight aquatic beetle species which inhabit different parts of the fresh—hypersaline gradient. Studied species belong to two unrelated genera which have invaded saline waters independently from freshwater ancestors; Nebrioporus (Dytiscidae) and Enochrus (Hydrophilidae). Their osmoregulatory strategy (osmoconformity or osmoregulation) was identified and osmotic capacity (the osmotic gradient between the animal’s haemolymph and the external medium) was compared between species pairs co-habiting similar salinities in nature. We show that osmoregulatory capacity, rather than osmoconformity, has evolved independently in these different lineages. All species hyperegulated their haemolymph osmotic concentration in diluted waters; those living in fresh or low-salinity waters were unable to hyporegulate and survive in hyperosmotic media (> 340 mosmol kg-1). In contrast, the species which inhabit the hypo-hypersaline habitats were effective hyporegulators, maintaining their haemolymph osmolality within narrow limits (ca. 300 mosmol kg-1) across a wide range of external concentrations. The hypersaline species N. ceresyi and E. jesusarribasi tolerated conductivities up to 140 and 180 mS cm-1, respectively, and maintained osmotic gradients over 3500 mosmol kg-1, comparable to those of the most effective insect osmoregulators known to date. Syntopic species of both genera showed similar osmotic capacities and in general, osmotic responses correlated well with upper salinity levels occupied by individual species in nature. Therefore, osmoregulatory capacity may mediate habitat segregation amongst congeners across the salinity gradient.

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.

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.

Aquatic insects dealing with dehydration: do desiccation resistance traits differ in species with contrasting habitat preferences?

Background Desiccation resistance shapes the distribution of terrestrial insects at multiple spatial scales. However, responses to drying stress have been poorly studied in aquatic groups, despite their potential role in constraining their distribution and diversification, particularly in arid and semi-arid regions. Methods We examined desiccation resistance in adults of four congeneric water beetle species (Enochrus, family Hydrophilidae) with contrasting habitat specificity (lentic vs. lotic systems and different salinity optima from fresh- to hypersaline waters). We measured survival, recovery capacity and key traits related to desiccation resistance (fresh mass, % water content, % cuticle content and water loss rate) under controlled exposure to desiccation, and explored their variability within and between species. Results Meso- and hypersaline species were more resistant to desiccation than freshwater and hyposaline ones, showing significantly lower water loss rates and higher water content. No clear patterns in desiccation resistance traits were observed between lotic and lentic species. Intraspecifically, water loss rate was positively related to specimens’ initial % water content, but not to fresh mass or % cuticle content, suggesting that the dynamic mechanism controlling water loss is mainly regulated by the amount of body water available. Discussion Our results support previous hypotheses suggesting that the evolution of desiccation resistance is associated with the colonization of saline habitats by aquatic beetles. The interespecific patterns observed in Enochrus also suggest that freshwater species may be more vulnerable than saline ones to drought intensification expected under climate change in semi-arid regions such as the Mediterranean Basin.

The contribution of mitochondrial metagenomics to large-scale data mining and phylogenetic analysis of Coleoptera

A phylogenetic tree at the species level is still far off for highly diverse insect orders, including the Coleoptera, but the taxonomic breadth of public sequence databases is growing. In addition, new types of data may contribute to increasing taxon coverage, such as metagenomic shotgun sequencing for assembly of mitogenomes from bulk specimen samples. The current study explores the application of these techniques for large-scale efforts to build the tree of Coleoptera. We used shotgun data from 17 different ecological and taxonomic datasets (5 unpublished) to assemble a total of 1942 mitogenome contigs of >3000 bp. These sequences were combined into a single dataset together with all mitochondrial data available at GenBank, in addition to nuclear markers widely used in molecular phylogenetics. The resulting matrix of nearly 16,000 species with two or more loci produced trees (RAxML) showing overall congruence with the Linnaean taxonomy at hierarchical levels from suborders to genera. We tested the role of full-length mitogenomes in stabilizing the tree from GenBank data, as mitogenomes might link terminals with non-overlapping gene representation. However, the mitogenome data were only partly useful in this respect, presumably because of the purely automated approach to assembly and gene delimitation, but improvements in future may be possible by using multiple assemblers and manual curation. In conclusion, the combination of data mining and metagenomic sequencing of bulk samples provided the largest phylogenetic tree of Coleoptera to date, which represents a summary of existing phylogenetic knowledge and a defensible tree of great utility, in particular for studies at the intra-familial level, despite some shortcomings for resolving basal nodes.