Search

With the on-going efforts in digitising museum collections, increased participation of citizen scientists, and greater accessibility to research data, accurately determining global patterns of diversity has become more achievable. Here, we used occurrence records from the Global Biodiversity Information Facility, with annotation of authoritative taxonomy, to evaluate the taxonomic richness of marine gastropods and to identify global species hotspots for this group. We also reviewed the availability of genetic resources within hotspots to detect potentially important regions where reference sequences for identifying these organisms are wanting. We find 33,268 unique and valid species under 3291 genera belonging to 380 gastropod families that have been recorded from 1662 to 2023. Globally, only 12.1% of reported species are linked to a COI barcode, whilst 27.1% of the families are represented with complete mitogenomes. Georeferenced records show the wide albeit disparate distribution of observations and species counts. The compiled dataset, published by organisations that are restricted to 55 countries, reveal a seeming lack of local submissions, but nonetheless demonstrates the growing contribution of citizen science platforms. We present 28 marine provinces across the globe as putative gastropod species hotspots and call for further work and stronger involvement, particularly within COI barcoding cold spots, to address the observed genetic reference inequity. Lastly, we highlight the important and relevant role of open and inclusive science to biodiversity monitoring and research.

Over the past decades the sampling of environmental DNA (eDNA) — encompassing the DNA of all organisms present in an environmental sample1 — has emerged as a technique for biodiversity monitoring and discovery in a diversity of environments. Avoiding the physical collection and identification of biota, this approach is praised for its independence of taxonomic expertise and has changed the way biologists study biodiversity. However, a common result in eDNA studies is the finding of unexpected taxa which are often removed by conservative bioinformatic filters or disregarded, since the authors are uncertain about the result and rarely have the interest, time, skills, and/or resources to return to the field and confirm with actual specimens2. Here, we report a case in which an eDNA discovery led to the physical localization of a member of the Micrognathozoa (Figure 1B) — a rare group of limnic micrometazoans, and the animal phylum to be discovered last3, which is the sister group to rotifers4,5. To this day, Micrognathozoa still comprises only a single named species from Greenland and a few additional disparate places

Over the past decades the sampling of environmental DNA (eDNA) — encompassing the DNA of all organisms present in an environmental sample1 — has emerged as a technique for biodiversity monitoring and discovery in a diversity of environments. Avoiding the physical collection and identification of biota, this approach is praised for its independence of taxonomic expertise and has changed the way biologists study biodiversity. However, a common result in eDNA studies is the finding of unexpected taxa which are often removed by conservative bioinformatic filters or disregarded, since the authors are uncertain about the result and rarely have the interest, time, skills, and/or resources to return to the field and confirm with actual specimens2. Here, we report a case in which an eDNA discovery led to the physical localization of a member of the Micrognathozoa (Figure 1B) — a rare group of limnic micrometazoans, and the animal phylum to be discovered last3, which is the sister group to rotifers4,5. To this day, Micrognathozoa still comprises only a single named species from Greenland and a few additional disparate places.

Information about the dietary composition of a species is crucial to understanding their position and role in the food web. Increasingly, molecular approaches such as DNA metabarcoding are used in studying trophic relationships, not least because they may alleviate problems such as low taxonomic resolution or underestimation of digestible taxa in the diet. Here, we used DNA metabarcoding with universal primers for cytochrome c oxidase I (COI) to study the diet composition of the northern shrimp (Pandalus borealis), an Arctic keystone species with large socio-economic importance. Across locations, jellyfish and chaetognaths were the most important components in the diet of P. borealis, jointly accounting for 40%–60% of the total read abundance. This dietary importance of gelatinous zooplankton contrasts sharply with published results based on stomach content analysis. At the same time, diet composition differed between fjord and shelf locations, pointing to different food webs supporting P. borealis in these two systems. Our study underlines the potential of molecular approaches to provide new insights into the diet of marine invertebrates that are difficult to obtain with traditional methods, and calls for a revision of the role of gelatinous zooplankton in the diet of the key Arctic species P. borealis, and in extension, Arctic food webs.

Ecosystem biomonitoring is crucial for proper management of natural communities during the Anthropocene era. With the advent of new sequencing technologies, DNA metabarcoding has been proposed as a game-changing tool for biomonitoring. In this Thesis we plead for the use of metabarcoding of a highly variable marker to infer not only the interspecies but also the intraspecies variability to assess both biogeographic, at the species level, and metaphylogeographic patterns, at the haplotype level. We focused on highly complex hard-substratum benthic littoral communities. The term "Metaphylogeography", coined in this Thesis, refers to the study of phylogeographic patterns of many species at the same time using metabarcoding data. However, as of the start of this Thesis, only a few studies had tested the metabarcoding method to directly characterize the whole eukaryotic community in highly diverse benthic ecosystems. This required to set up and calibrate methods for these communities as a prior step. We first evaluated both the sampling methods and the bioinformatic pipelines. We assessed the viability of detecting the environmental DNA released from the benthic community into the adjacent water layer using metabarcoding of COI with highly degenerated primers targeting the whole eukaryotic community. We sampled water from 0 to 20m from shallow rocky benthic communities and compared the DNA signal with the results obtained from metabarcoding directly the benthic communities by traditional quadrat sampling. We also designed a pipeline combining clustering and denoising methods to treat metabarcoding data of COI. We considered the entropy of each codon position of this coding fragment both to improve the detection of spurious sequences and to calibrate the best performing parameters of the software used. In addition, we created our own denoising program, DnoisE, to incorporate information on the codon position. This new code and parameter calibration were required as the co

Information about the dietary composition of a species is crucial to understanding their position and role in the food web. Increasingly, molecular approaches such as DNA metabarcoding are used in studying trophic relationships, not least because they may alleviate problems such as low taxonomic resolution or underestimation of digestible taxa in the diet. Here, we used DNA metabarcoding with universal primers for cytochrome c oxidase I (COI) to study the diet composition of the northern shrimp (Pandalus borealis), an Arctic keystone species with large socio-economic importance. Across locations, jellyfish and chaetognaths were the most important components in the diet of P. borealis, jointly accounting for 40%–60% of the total read abundance. This dietary importance of gelatinous zooplankton contrasts sharply with published results based on stomach content analysis. At the same time, diet composition differed between fjord and shelf locations, pointing to different food webs supporting P. borealis in these two systems. Our study underlines the potential of molecular approaches to provide new insights into the diet of marine invertebrates that are difficult to obtain with traditional methods, and calls for a revision of the role of gelatinous zooplankton in the diet of the key Arctic species P. borealis, and in extension, Arctic food webs.

Information about the dietary composition of a species is crucial to understanding their position and role in the food web. Increasingly, molecular approaches such as DNA metabarcoding are used in studying trophic relationships, not least because they may alleviate problems such as low taxonomic resolution or underestimation of digestible taxa in the diet. Here, we used DNA metabarcoding with universal primers for cytochrome c oxidase I (COI) to study the diet composition of the northern shrimp (Pandalus borealis), an Arctic keystone species with large socio-economic importance. Across locations, jellyfish and chaetognaths were the most important components in the diet of P. borealis, jointly accounting for 40%–60% of the total read abundance. This dietary importance of gelatinous zooplankton contrasts sharply with published results based on stomach content analysis. At the same time, diet composition differed between fjord and shelf locations, pointing to different food webs supporting P. borealis in these two systems. Our study underlines the potential of molecular approaches to provide new insights into the diet of marine invertebrates that are difficult to obtain with traditional methods, and calls for a revision of the role of gelatinous zooplankton in the diet of the key Arctic species P. borealis, and in extension, Arctic food webs.

DNA metabarcoding is broadly used in biodiversity studies encompassing a wide range of organisms. Erroneous amplicons, generated during amplification and sequencing procedures, constitute one of the major sources of concern for the interpretation of metabarcoding results. Several denoising programs have been implemented to detect and eliminate these errors. However, almost all denoising software currently available has been designed to process non-coding ribosomal sequences, most notably prokaryotic 16S rDNA. The growing number of metabarcoding studies using coding markers such as COI or RuBisCO demands a re-assessment and calibration of denoising algorithms. Here we present DnoisE, the first denoising program designed to detect erroneous reads and merge them with the correct ones using information from the natural variability (entropy) associated to each codon position in coding barcodes. We have developed an open-source software using a modified version of the UNOISE algorithm. DnoisE implements different merging procedures as options, and can incorporate codon entropy information either retrieved from the data or supplied by the user. In addition, the algorithm of DnoisE is parallelizable, greatly reducing runtimes on computer clusters. Our program also allows different input file formats, so it can be readily incorporated into existing metabarcoding pipelines.

Information about the dietary composition of a species is crucial to understanding their position and role in the food web. Increasingly, molecular approaches such as DNA metabarcoding are used in studying trophic relationships, not least because they may alleviate problems such as low taxonomic resolution or underestimation of digestible taxa in the diet. Here, we used DNA metabarcoding with universal primers for cytochrome c oxidase I (COI) to study the diet composition of the northern shrimp (Pandalus borealis), an Arctic keystone species with large socio-economic importance. Across locations, jellyfish and chaetognaths were the most important components in the diet of P. borealis, jointly accounting for 40%–60% of the total read abundance. This dietary importance of gelatinous zooplankton contrasts sharply with published results based on stomach content analysis. At the same time, diet composition differed between fjord and shelf locations, pointing to different food webs supporting P. borealis in these two systems. Our study underlines the potential of molecular approaches to provide new insights into the diet of marine invertebrates that are difficult to obtain with traditional methods, and calls for a revision of the role of gelatinous zooplankton in the diet of the key Arctic species P. borealis, and in extension, Arctic food webs.

Traditional morphological methods for species identification are highly time consuming, especially for small organisms, such as Foraminifera, a group of shell-building microbial eukaryotes. To analyze large amounts of samples more efficiently, species identification methods have extended to molecular tools in the last few decades. Although a wide range of phyla have good markers available, for Foraminifera only one hypervariable marker from the ribosomal region (18S) is widely used. Recently a new mitochondrial marker cytochrome oxidase subunit 1 (COI) has been sequenced. Here we investigate whether this marker has a higher potential for species identification compared to the ribosomal marker. We explore the genetic variability of both the 18S and COI markers in 22 benthic foraminiferal morphospecies (orders Miliolida and Rotaliida). Using single-cell DNA, the genetic variability within specimens (intra) and between specimens (inter) of each species was assessed using next-generation sequencing. Amplification success rate was twice as high for COI (151/200 specimens) than for 18S (73/200 specimens). The COI marker showed greatly decreased intra- and inter-specimen variability compared to 18S in six out of seven selected species. The 18S phylogenetic reconstruction fails to adequately cluster multiple species together in contrast to COI. Additionally, the COI marker helped recognize misclassified specimens difficult to morphologically identify to the species level. Integrative taxonomy, combining morphological and molecular characteristics, provides a robust picture of the foraminiferal species diversity. Finally, we suggest the use of a set of sequences (two or more) to describe species showing intra-genomic variability additionally to using multiple markers. Our findings highlight the potential of the newly discovered mitochondrial marker for molecular species identification and metabarcoding purposes.

In a recent paper, “Environmental DNA: What's behind the term? Clarifying the terminology and recommendations for its future use in biomonitoring,” Pawlowski et al. argue that the term eDNA should be used to refer to the pool of DNA isolated from environmental samples, as opposed to only extra-organismal DNA from macro-organisms. We agree with this view. However, we are concerned that their proposed two-level terminology specifying sampling environment and targeted taxa is overly simplistic and might hinder rather than improve clear communication about environmental DNA and its use in biomonitoring. This terminology is based on categories that are often difficult to assign and uninformative, and it overlooks a fundamental distinction within eDNA: the type of DNA (organismal or extra-organismal) from which ecological interpretations are derived.

Marine biodiversity underpins ecosystem health and societal well-being. Preservation of biodiversity hotspots is a global challenge. Molecular tools, like DNA barcoding and metabarcoding, hold great potential for biodiversity monitoring, possibly outperforming more traditional taxonomic methods. However, metabarcoding-based biodiversity assessments are limited by the availability of sequences in barcoding reference databases; a lack thereof results in high percentages of unassigned sequences. In this study we (i) present the current status of known vs. barcoded marine species at a global scale based on online taxonomic and genetic databases; and (ii) compare the current status with data from ten years ago. Then we analyzed occurrence data of marine animal species from five Large Marine Ecosystems (LMEs) classified as biodiversity hotspots, to identify any consistent disparities in COI barcoding coverage between geographic regions and at phylum level. Barcoding coverage varied among LMEs (from 36.8% to 62.4% COI-barcoded species) and phyla (from 4.8% to 74.7% COI-barcoded species), with Porifera, Bryozoa and Platyhelminthes being highly underrepresented, compared to Chordata, Arthropoda and Mollusca. We demonstrate that although barcoded marine species increased from 9.5% to 14.2% since the last assessment in 2011, about 15,000 (corresponding to 7.8% increase) new species were described from 2011 to 2021. The next ten years will thus be crucial to enroll concrete collaborative measures and long term initiatives (e.g., Horizon 2030, Ocean Decade) to populate barcoding libraries for the marine realm.

Background: The recent blooming of metabarcoding applications to biodiversity studies comes with some relevant methodological debates. One such issue concerns the treatment of reads by denoising or by clustering methods, which have been wrongly presented as alternatives. It has also been suggested that denoised sequence variants should replace clusters as the basic unit of metabarcoding analyses, missing the fact that sequence clusters are a proxy for species-level entities, the basic unit in biodiversity studies. We argue here that methods developed and tested for ribosomal markers have been uncritically applied to highly variable markers such as cytochrome oxidase I (COI) without conceptual or operational (e.g., parameter setting) adjustment. COI has a naturally high intraspecies variability that should be assessed and reported, as it is a source of highly valuable information. We contend that denoising and clustering are not alternatives. Rather, they are complementary and both should be used together in COI metabarcoding pipelines. Results: Using a COI dataset from benthic marine communities, we compared two denoising procedures (based on the UNOISE3 and the DADA2 algorithms), set suitable parameters for denoising and clustering, and applied these steps in different orders. Our results indicated that the UNOISE3 algorithm preserved a higher intra-cluster variability. We introduce the program DnoisE to implement the UNOISE3 algorithm taking into account the natural variability (measured as entropy) of each codon position in protein-coding genes. This correction increased the number of sequences retained by 88%. The order of the steps (denoising and clustering) had little influence on the final outcome. Conclusions: We highlight the need for combining denoising and clustering, with adequate choice of stringency parameters, in COI metabarcoding. We present a program that uses the coding properties of this marker to improve the denoising step. We recommend researche

Global environmental changes may have a profound impact on ecosystems. In this context, it is crucial to gather biological and ecological information of the main species in marine communities to predict and mitigate potential effects of shifts in their distribution, abundance, and interactions. Using genotyping by sequencing (GBS), we assessed the genetic structure of a keystone species in the Mediterranean shallow littoral ecosystems, the black sea urchin Arbacia lixula. This bioengineer species can shape their communities due to its grazing activity and it is experiencing an ongoing expansion with increasing temperatures. The population genomic analyses on 5,241 loci sequenced in 240 individuals from 11 Mediterranean sampled populations revealed that all populations were diverse and showed significant departure from equilibrium. Albeit genetic differentiation was in general shallow, a significant break separated the western and eastern Mediterranean populations, a break not detected in previous studies with less resolutive markers. Notably, no clear effect of the Almería-Oran front, an important break in the Atlanto-Mediterranean transition, could be detected among the western basin populations, where only a slight differentiation of the two northernmost populations was found. Despite the generally low levels of genetic differentiation found, we identified candidate regions for local adaptation by combining different genomic analysis with environmental data. Salinity, rather than temperature, seemed to be an important driver of genetic structure in A. lixula. Overall, from a population genomics standpoint, there is ample scope for A. lixula to continue thriving and adapting in the warming Mediterranean.

In a recent paper, “Environmental DNA: What's behind the term? Clarifying the terminology and recommendations for its future use in biomonitoring,” Pawlowski et al. argue that the term eDNA should be used to refer to the pool of DNA isolated from environmental samples, as opposed to only extra-organismal DNA from macro-organisms. We agree with this view. However, we are concerned that their proposed two-level terminology specifying sampling environment and targeted taxa is overly simplistic and might hinder rather than improve clear communication about environmental DNA and its use in biomonitoring. This terminology is based on categories that are often difficult to assign and uninformative, and it overlooks a fundamental distinction within eDNA: the type of DNA (organismal or extra-organismal) from which ecological interpretations are derived.

Reproductive rates of copepods are temperature-dependent, but poorly known for small copepods at low temperatures, hindering the predictions of population dynamics and secondary production in high-latitude ecosystems. We investigated egg hatching rates, hatching success and egg production of the small copepods Oithona similis and Microsetella norvegica (sac spawners) and Microcalanus pusillus (broadcast spawner) between March and August. Incubations were performed at ecologically relevant temperatures between 1.3 and 13.2°C, and egg production rates were calculated. All egg hatching rates were positively correlated to temperature, although with large species-specific differences. At the lowest temperatures, M. pusillus eggs hatched within 4 days, whereas the eggs from sac spawners took 3–8 weeks to hatch. The egg hatching success was ≤25% for M. pusillus, >75% for O. similis and variable for M. norvegica. The maximum weight-specific egg production rate (μg C μg−1 C d−1) of M. pusillus was higher (0.22) than O. similis (0.12) and M. norvegica (0.06). M. norvegica reproduction peaked at 6–8°C, the prevailing in situ temperatures during its reproductive period. The difference in reproductive rates indicates species-specific thermal plasticity for the three copepods, which could have implications for present and future population dynamics of the species in arctic fjords.

In the marine realm, biomonitoring using eDNA of benthic communities requires destructive direct sampling or the setting-up of settlement structures. Comparatively much less effort is required to sample the water column, which can be accessed remotely. In this study we assess the feasibility of obtaining information from the eukaryotic benthic communities by sampling the adjacent water layer. We studied two different rocky-substrate benthic communities with a technique based on quadrat sampling. We also took replicate water samples at distances from a few centimetres to 20 m from the benthic habitat. Using as marker a fragment of the Cytochrome c oxidase subunit I gene with universal primers, we obtained a total of 3,543 molecular operational taxonomic units (MOTUs) from the samples. The structure obtained in the two environments was markedly different, with Metazoa, Archaeplastida Rhodophyta and Stramenopiles being the most diverse group in benthic samples, and HacrobiaAlveolata, Metazoa and Alveolata Rhizaria in the water. Only 265 MOTUs (7.5%) were shared between benthos and water samples, and of these 180 MOTUs (5.1%) were identified as benthic MOTUs that left their DNA in the water. Most of them were found immediately adjacent to the benthos, and their number decreased and The distribution of these benthic shared MOTUs showed a decrease both in number of MOTUs and in number of reads as we moved apart from the benthic habitat. It was concluded that water eDNA, even in the close vicinity of the benthos, was a poor proxy for the analysis of benthic structure, and that direct sampling methods are required for monitoring these complex benthic communities via metabarcoding.

Sponges have recently been proposed as ideal candidates to act as natural samplers for environmental DNA due to their efficiency in filtering water. However, validation of the usefulness of DNA recovered from sponges to reveal vertebrate biodiversity patterns in Marine Protected Areas is still needed. Additionally, nothing is known about how different sponge species and morphologies influence the capture of environmental DNA and whether biodiversity patterns obtained from sponges are best described by quantitative or qualitative measures. In this study, we amplified and sequenced a vertebrate specific 12S barcode with a set of universal PCR primers (MiFish) for metabarcoding environmental DNA from fishes, to unveil fine-scale patterns of fish communities from natural-sampler DNA retrieved from 64 sponges (16 species) located in eutrophic and well-preserved coral reefs in Nha Trang Bay (central Vietnam). Ninety tropical fish species were identified from the sponges, corresponding to one third of the total local ichthyofauna reported from previous extensive conventional surveys. Significant differentiation in fish communities between eutrophic and well-preserved environments was observed, albeit eutrophication only explained a modest proportion of the variation between fish communities. Differences in efficiency of capturing fish environmental DNA among sponge species or morphologies were not observed. Overall, the majority of detected fish species corresponded to reef-associated small-sized species, as expected in coral reefs environments. Remarkably, pelagic, migratory, and deep-sea fish species were also recovered from sponge tissues, pointing out the ability of sponge natural sampled DNA to detect fishes that were not permanently associated to the biomes where the sponges were sampled. These results highlight the suitability of natural samplers as a cost-effective way to assess vertebrate diversity in hyper-diverse environments.

In the marine realm, biomonitoring using environmental DNA (eDNA) of benthic communities requires destructive direct sampling or the setting-up of settlement structures. Comparatively much less effort is required to sample the water column, which can be accessed remotely. In this study we assess the feasibility of obtaining information from the eukaryotic benthic communities by sampling the adjacent water layer. We studied two different rocky-substrate benthic communities with a technique based on quadrat sampling. We also took replicate water samples at four distances (0, 0.5, 1.5, and 20 m) from the benthic habitat. Using broad range primers to amplify a ca. 313 bp fragment of the cytochrome oxidase subunit I gene, we obtained a total of 3,543 molecular operational taxonomic units (MOTUs). The structure obtained in the two environments was markedly different, with Metazoa, Archaeplastida and Stramenopiles being the most diverse groups in benthic samples, and Hacrobia, Metazoa and Alveolata in the water. Only 265 MOTUs (7.5%) were shared between benthos and water samples and, of these, 180 (5.1%) were identified as benthic taxa that left their DNA in the water. Most of them were found immediately adjacent to the benthos, and their number decreased as we moved apart from the benthic habitat. It was concluded that water eDNA, even in the close vicinity of the benthos, was a poor proxy for the analysis of benthic structure, and that direct sampling methods are required for monitoring these complex communities via metabarcoding.

Information about the genomic processes underlying responses to temperature changes is still limited in non-model marine invertebrates. In this sense, transcriptomic analyses can help to identify genes potentially related to thermal responses. We here investigated, via RNA-seq, whole-transcriptomic responses to increased and decreased temperatures in a thermophilous keystone sea urchin, Arbacia lixula, whose populations are increasing in the Mediterranean. This species is a key driver of benthic communities’ structure due to its grazing activity. We found a strong response to experimentally induced cold temperature (7°C), with 1,181 differentially expressed transcripts relative to the control condition (13°C), compared to only 179 in the warm (22°C) treatment. A total of 84 (cold treatment) and three (warm treatment) gene ontology terms were linked to the differentially expressed transcripts. At 7°C the expression of genes encoding different heat shock proteins (HSPs) was upregulated, together with apoptotic suppressor genes (e.g., Bcl2), genes involved in the infection response and/or pathogen-recognition (e.g., echinoidin) and ATP-associated genes, while protein biosynthesis and DNA replication pathways were downregulated. At 22°C neither HSPs induction nor activation of the previously mentioned pathways were detected, with the exception of some apoptotic-related activities that were upregulated. Our results suggest a strong transcriptional response associated with low temperatures, and support the idea of low water temperature being a major limitation for A. lixula expansion across deep Mediterranean and northern Atlantic waters.

Among the complex ecosystems and habitats that form the deep sea, submarine canyons and open slope systems are regarded as potential hotspots of biodiversity. We assessed the spatial and temporal patterns of biodiversity in sediment communities of a NW Mediterranean Canyon and its adjacent open slope (Blanes Canyon) with DNA metabarcoding. We sampled three layers of sediment and four different depths (900–1750 m) at two seasons, and used a fragment of the mitochondrial gene cytochrome c oxidase subunit I (COI) as a metabarcoding marker. The final dataset contained a total of 15,318 molecular operational taxonomic units (MOTUs). Metazoa, Stramenopiles and Archaeplastida were the dominant taxa and, within metazoans, Arthropoda, Nematoda and Cnidaria were the most diverse. There was a trend towards decreasing MOTU richness and diversity in the first few cm (1 to 5) of the sediment, with only 26.3% of the MOTUs shared across sediment layers. Our results show the presence of heterogeneous communities in the studied area, which was significantly different between zones, depths and seasons. We compared our results with the ones presented in a previous study, obtained using the v7 region of the 18S rRNA gene in the same samples. There were remarkable differences in the total number of MOTUs and in the most diverse taxa. COI recovered a higher number of MOTUs, but more remained unassigned taxonomically. However, the broad spatio-temporal patterns elucidated from both datasets coincided, with both markers retrieving the same ecological information. Our results showed that COI can be used to accurately characterize the studied communities and constitute a high-resolution method to detect ecological shifts. We also highlight that COI reference databases for deep-sea organisms have important gaps, and their completeness is essential in order to successfully apply metabarcoding techniques.

Among the complex ecosystems and habitats that form the deep sea, submarine canyons and open slope systems are regarded as potential hotspots of biodiversity. We assessed the spatial and temporal patterns of biodiversity in sediment communities of a NW Mediterranean Canyon and its adjacent open slope (Blanes Canyon) with DNA metabarcoding. We sampled three layers of sediment and four different depths (900–1750 m) at two seasons, and used a fragment of the mitochondrial gene cytochrome c oxidase subunit I (COI) as a metabarcoding marker. The final dataset contained a total of 15,318 molecular operational taxonomic units (MOTUs). Metazoa, Stramenopiles and Archaeplastida were the dominant taxa and, within metazoans, Arthropoda, Nematoda and Cnidaria were the most diverse. There was a trend towards decreasing MOTU richness and diversity in the first few cm (1 to 5) of the sediment, with only 26.3% of the MOTUs shared across sediment layers. Our results show the presence of heterogeneous communities in the studied area, which was significantly different between zones, depths and seasons. We compared our results with the ones presented in a previous study, obtained using the v7 region of the 18S rRNA gene in the same samples. There were remarkable differences in the total number of MOTUs and in the most diverse taxa. COI recovered a higher number of MOTUs, but more remained unassigned taxonomically. However, the broad spatio-temporal patterns elucidated from both datasets coincided, with both markers retrieving the same ecological information. Our results showed that COI can be used to accurately characterize the studied communities and constitute a high-resolution method to detect ecological shifts. We also highlight that COI reference databases for deep-sea organisms have important gaps, and their completeness is essential in order to successfully apply metabarcoding techniques.

Aim: The Atlanto-Mediterranean edible purple sea urchin, Paracentrotus lividus, is a commercially exploited keystone species in benthic communities. Its browsing activity can deeply modify the littoral landscape, and changes in its abundance are of major conservation concern. This species is facing nowadays contrasting anthropogenic pressures linked to predator release, exploitation and sea warming. Management of this key species requires knowledge of its genetic structure, connectivity and local adaptation. Our goal was to assess the current global status of the species under a genomic perspective. Location: Atlanto-Mediterranean shores from Morocco and France to Turkey. Methods: We used genotyping by sequencing (GBS) of 241 individuals belonging to 11 populations spanning the known range of distribution of the species. We obtained 3,348 loci for population genomics and outlier analyses. Results: We identified significant genetic structure and a gradient matching the longitudinal position of the localities. A hierarchical analysis revealed two main clusters (Atlantic and Mediterranean) and subtler patterns of differentiation within them. Candidate markers for selection identified between and within these two main clusters were mostly different, likely indicating different selective pressures. Adaptation to maximum salinity and maximum temperature appeared as an important driver of the transition between Atlantic and Mediterranean basins. Other stressors, such as minimum temperature or range of temperature, seem to define the structuring within the Mediterranean. Main conclusions: Our study shows the potential of hierarchical analyses on population genomics to detect fine-scale genetic structure and adaptation signatures in marine species with long dispersal capabilities. Although genetic interchange occurs widely in Paracentrotus lividus, the species is sensitive to dispersal barriers, displays isolation by distance and faces local selective pressures associated to

Peat bogs located in high mountains are suitable places to study local environmental responses to climate variability. These ecosystems host a large number of eukaryotes with diverse taxonomic and functional diversity. We carried out a metabarcoding study using universal 18S and COI markers to explore the composition of past and present eukaryotic communities of a Pyrenean peat bog ecosystem. We assessed the molecular biodiversity of four different moss micro-habitats along a flood gradient in the lentic Bassa Nera system (Central Pyrenees). Five samples collected from different sediment depths at the same study site were also analysed, to test the suitability of these universal markers for studying paleoecological communities recovered from ancient DNA and to compare the detected DNA sequences to those obtained from the modern community. We also compared the information provided by the sedimentary DNA to the reconstruction from environmental proxies such as pollen and macro-remains from the same record. We successfully amplified ancient DNA with both universal markers from all sediment samples, including the deepest one (* 10,000 years old). Most of the metabarcoding reads obtained from sediment samples, however, were assigned to living edaphic organisms and only a small fraction of those reads was considered to be derived from paleoecological communities. Inferences from ancient sedimentary DNA were complementary to the reconstruction based on pollen and macro-remains, and the combined records reveal more detailed information. This molecular study yielded promising findings regarding the diversity of modern eukaryotic peat bog communities. Nevertheless, even though information about past communities could be retrieved from sediment samples, preferential amplification of DNA from living communities is a caveat for the use of universal metabarcoding markers in paleoecology.

The genetic structure of 13 populations of the amphiatlantic sea urchin Arbacia lixula, as well as temporal genetic changes in three of these localities, were assessed using ten hypervariable microsatellite loci. This thermophilous sea urchin is an important engineer species triggering the formation of barren grounds through its grazing activity. Its abundance seems to be increasing in most parts of the Mediterranean, probably favoured by warming conditions. Significant genetic differentiation was found both spatially and temporally. The main break corresponded to the separation of western Atlantic populations from those in eastern Atlantic and the Mediterranean Sea. A less marked, but significant differentiation was also found between Macaronesia (eastern Atlantic) and the Mediterranean. In the latter area, a signal of differentiation between the transitional area (Alboran Sea) and the rest of the Mediterranean was detected. However, no genetic structure is found within the Mediterranean (excluding Alboran) across the Siculo-Tunisian Strait, resulting from either enough gene flow to homogenize distance areas or/and a recent evolutionary history marked by demographic expansion in this basin. Genetic temporal variation at the Alboran Sea is as important as spatial variation, suggesting that temporal changes in hydrological features can affect the genetic composition of the populations. A picture of genetic homogeneity in the Mediterranean emerges, implying that the potential expansion of this keystone species will not be limited by intraspecific genetic features and/or potential impact of postulated barriers to gene flow in the region.

In the era of "Anthropocene defaunation," large species are often no longer detected in habitats where they formerly occurred. However, it is unclear whether this apparent missing, or "dark," diversity of megafauna results from local species extirpations or from failure to detect elusive remaining individuals. We find that despite two orders of magnitude less sampling effort, environmental DNA (eDNA) detects 44% more shark species than traditional underwater visual censuses and baited videos across the New Caledonian archipelago (south-western Pacific). Furthermore, eDNA analysis reveals the presence of previously unobserved shark species in human-impacted areas. Overall, our results highlight a greater prevalence of sharks than described by traditional survey methods in both impacted and wilderness areas. This indicates an urgent need for large-scale eDNA assessments to improve monitoring of threatened and elusive megafauna. Finally, our findings emphasize the need for conservation efforts specifically geared toward the protection of elusive, residual populations.

In the era of "Anthropocene defaunation," large species are often no longer detected in habitats where they formerly occurred. However, it is unclear whether this apparent missing, or "dark," diversity of megafauna results from local species extirpations or from failure to detect elusive remaining individuals. We find that despite two orders of magnitude less sampling effort, environmental DNA (eDNA) detects 44% more shark species than traditional underwater visual censuses and baited videos across the New Caledonian archipelago (south-western Pacific). Furthermore, eDNA analysis reveals the presence of previously unobserved shark species in human-impacted areas. Overall, our results highlight a greater prevalence of sharks than described by traditional survey methods in both impacted and wilderness areas. This indicates an urgent need for large-scale eDNA assessments to improve monitoring of threatened and elusive megafauna. Finally, our findings emphasize the need for conservation efforts specifically geared toward the protection of elusive, residual populations.

In the era of "Anthropocene defaunation," large species are often no longer detected in habitats where they formerly occurred. However, it is unclear whether this apparent missing, or "dark," diversity of megafauna results from local species extirpations or from failure to detect elusive remaining individuals. We find that despite two orders of magnitude less sampling effort, environmental DNA (eDNA) detects 44% more shark species than traditional underwater visual censuses and baited videos across the New Caledonian archipelago (south-western Pacific). Furthermore, eDNA analysis reveals the presence of previously unobserved shark species in human-impacted areas. Overall, our results highlight a greater prevalence of sharks than described by traditional survey methods in both impacted and wilderness areas. This indicates an urgent need for large-scale eDNA assessments to improve monitoring of threatened and elusive megafauna. Finally, our findings emphasize the need for conservation efforts specifically geared toward the protection of elusive, residual populations.

In the era of "Anthropocene defaunation," large species are often no longer detected in habitats where they formerly occurred. However, it is unclear whether this apparent missing, or "dark," diversity of megafauna results from local species extirpations or from failure to detect elusive remaining individuals. We find that despite two orders of magnitude less sampling effort, environmental DNA (eDNA) detects 44% more shark species than traditional underwater visual censuses and baited videos across the New Caledonian archipelago (south-western Pacific). Furthermore, eDNA analysis reveals the presence of previously unobserved shark species in human-impacted areas. Overall, our results highlight a greater prevalence of sharks than described by traditional survey methods in both impacted and wilderness areas. This indicates an urgent need for large-scale eDNA assessments to improve monitoring of threatened and elusive megafauna. Finally, our findings emphasize the need for conservation efforts specifically geared toward the protection of elusive, residual populations.

Biodiversity assessment of marine hard-bottom communities is hindered by the high diversity and size-ranges of the organisms present.Wedeveloped aDNAmetabarcoding protocol for biodiversity characterization of structurally complex natural marine hardbottom communities. We used two molecular markers: the ``Leray fragment'' of mitochondrial cytochrome c oxidase (COI), for which a novel primer set was developed, and the V7 region of the nuclear small subunit ribosomal RNA (18S). Eight different shallow marine littoral communities from two National Parks in Spain (one in the Atlantic Ocean and another in the Mediterranean Sea) were studied. Samples were sieved into three size fractions from where DNA was extracted separately. Bayesian clustering was used for delimiting molecular operational taxonomic units (MOTUs) and custom reference databases were constructed for taxonomic assignment. Despite applying stringent filters, we found high values for MOTU richness (2,510 and 9,679 MOTUs with 18S and COI, respectively), suggesting that these communities host a large amount of yet undescribed eukaryotic biodiversity. Significant gaps are still found in sequence reference databases, which currently prevent the complete taxonomic assignment of the detected sequences. In our dataset, 85% of 18S MOTUs and 64% of COI MOTUs could be identified to phylum or lower taxonomic level. Nevertheless, those unassigned were mostly rare MOTUs with low numbers of reads, and assigned MOTUs comprised over 90% of the total sequence reads. The identification rate might be significantly improved in the future, as reference databases are further completed. Our results show that marine metabarcoding, currently applied mostly to plankton or sediments, can be adapted to structurally complex hard bottom samples. Thus, eukaryotic metabarcoding emerges as a robust, fast, objective and affordable method to comprehensively characterize the diversity of marine benthic communities dominated by macroscopic seaw

We analysed with multigene (18S and COI) metabarcoding the effects of the proliferation of invasive seaweeds on rocky littoral communities in two Spanish Marine Protected Areas. The invasive algae studied were Caulerpa cylindracea, Lophocladia lallemandii and Asparagopsis armata. They are canopy-forming, landscape-dominant seaweeds, and we were interested in their effects on the underlying communities of meiobenthos and macrobenthos, separated in two size fractions through sieving. A new semiquantitative treatment of metabarcoding data is introduced. The results for both markers showed that the presence of the invasive seaweed had a significant effect on the understory communities for Lophocladia lallemandii and Asparagopsis armata but not for Caulerpa cylindracea. Likewise, changes in MOTU richness and diversity with invasion status varied in magnitude and direction depending on the alga considered. Our results showed that metabarcoding allows monitoring of the less conspicuous, but not least important, effects of the presence of dominant invasive seaweeds.

Sharks are charismatic predators that play a key role in most marine food webs. Their demonstrated vulnerability to exploitation has recently turned them into flagship species in ocean conservation. Yet, the assessment and monitoring of the distribution and abundance of such mobile species in marine environments remain challenging, often invasive and resource-intensive. Here we pilot a novel, rapid and non-invasive environmental DNA (eDNA) metabarcoding approach specifically targeted to infer shark presence, diversity and eDNA read abundance in tropical habitats. We identified at least 21 shark species, from both Caribbean and Pacific Coral Sea water samples, whose geographical patterns of diversity and read abundance coincide with geographical differences in levels of anthropogenic pressure and conservation effort. We demonstrate that eDNA metabarcoding can be effectively employed to study shark diversity. Further developments in this field have the potential to drastically enhance our ability to assess and monitor elusive oceanic predators, and lead to improved conservation strategies.

Sharks are charismatic predators that play a key role in most marine food webs. Their demonstrated vulnerability to exploitation has recently turned them into flagship species in ocean conservation. Yet, the assessment and monitoring of the distribution and abundance of such mobile species in marine environments remain challenging, often invasive and resource-intensive. Here we pilot a novel, rapid and non-invasive environmental DNA (eDNA) metabarcoding approach specifically targeted to infer shark presence, diversity and eDNA read abundance in tropical habitats. We identified at least 21 shark species, from both Caribbean and Pacific Coral Sea water samples, whose geographical patterns of diversity and read abundance coincide with geographical differences in levels of anthropogenic pressure and conservation effort. We demonstrate that eDNA metabarcoding can be effectively employed to study shark diversity. Further developments in this field have the potential to drastically enhance our ability to assess and monitor elusive oceanic predators, and lead to improved conservation strategies.

Sharks are charismatic predators that play a key role in most marine food webs. Their demonstrated vulnerability to exploitation has recently turned them into flagship species in ocean conservation. Yet, the assessment and monitoring of the distribution and abundance of such mobile species in marine environments remain challenging, often invasive and resource-intensive. Here we pilot a novel, rapid and non-invasive environmental DNA (eDNA) metabarcoding approach specifically targeted to infer shark presence, diversity and eDNA read abundance in tropical habitats. We identified at least 21 shark species, from both Caribbean and Pacific Coral Sea water samples, whose geographical patterns of diversity and read abundance coincide with geographical differences in levels of anthropogenic pressure and conservation effort. We demonstrate that eDNA metabarcoding can be effectively employed to study shark diversity. Further developments in this field have the potential to drastically enhance our ability to assess and monitor elusive oceanic predators, and lead to improved conservation strategies.

Sharks are charismatic predators that play a key role in most marine food webs. Their demonstrated vulnerability to exploitation has recently turned them into flagship species in ocean conservation. Yet, the assessment and monitoring of the distribution and abundance of such mobile species in marine environments remain challenging, often invasive and resource-intensive. Here we pilot a novel, rapid and non-invasive environmental DNA (eDNA) metabarcoding approach specifically targeted to infer shark presence, diversity and eDNA read abundance in tropical habitats. We identified at least 21 shark species, from both Caribbean and Pacific Coral Sea water samples, whose geographical patterns of diversity and read abundance coincide with geographical differences in levels of anthropogenic pressure and conservation effort. We demonstrate that eDNA metabarcoding can be effectively employed to study shark diversity. Further developments in this field have the potential to drastically enhance our ability to assess and monitor elusive oceanic predators, and lead to improved conservation strategies.

We assessed spatio-temporal patterns of diversity in deep-sea sediment communities using metabarcoding. We chose a recently developed eukaryotic marker based on the v7 region of the 18S rRNA gene. Our study was performed in a submarine canyon and its adjacent slope in the Northwestern Mediterranean Sea, sampled along a depth gradient at two different seasons. We found a total of 5,569 molecular operational taxonomic units (MOTUs), dominated by Metazoa, Alveolata and Rhizaria. Among metazoans, Nematoda, Arthropoda and Annelida were the most diverse. We found a marked heterogeneity at all scales, with important differences between layers of sediment and significant changes in community composition with zone (canyon vs slope), depth, and season. We compared the information obtained from metabarcoding DNA and RNA and found more total MOTUs and more MOTUs per sample with DNA (ca. 20% and 40% increase, respectively). Both datasets showed overall similar spatial trends, but most groups had higher MOTU richness with the DNA template, while others, such as nematodes, were more diverse in the RNA dataset. We provide metabarcoding protocols and guidelines for biomonitoring of these key communities in order to generate information applicable to management efforts.

Competition, predation, and facilitation mechanisms are the major drivers of biodiversity and community structure in marine benthic ecosystems. Habitat complexity is a determining factor of faunal richness and biodiversity in these communities. The structure of marine animal forests is originated by living threedimensional aggregations of modular animals. The persistence of these systems through time relies on the growth of existing individuals and the recruitment of new ones. Therefore, the present and future health of these valuable ecosystems may depend on the reproductive success of a few vulnerable species which might often be accomplished only under strict or very narrow conditions. Reproductive patterns of ecosystem engineers play a crucial role in determining the structure, function, and distribution of all kinds of marine animal forests at different scales. The reproductive strategies of these habitat-forming species may vary considerably. Though most ecosystem-engineering species are, to some extent, able to reproduce asexually, their sexual reproductive strategies are diverse. Dispersal ability strategies are selected as an adaptation to spatial heterogeneity and habitat stability and are important factors for the resilience of the ecosystems. Dispersal traits are essential for both small-scale population structuring and connectivity among distant populations. Disturbed habitats may be promptly recolonized by species with long-distance dispersal capability, but recolonization by species with low dispersal capacity might limit the full restoration of a disturbed ecosystem, especially in fragmented habitats with reduced connectivity between patches. Due to the lack of knowledge on the reproductive cycles of many marine invertebrates, the response of animal forest ecosystems to global change is, in general, unpredictable.

The “marine animal forests” are among the most diverse ecosystems in the Biosphere. However, exhaustive biodiversity assessment of these communities has been so far elusive. The real extent of biodiversity and its temporal and spatial variability patterns remain unknown for most animal forests, mainly due to the inability of traditional taxonomy methods to cope with such degree of diversity and structural complexity. The development of metabarcoding techniques has revolutionized biomonitoring. Using this approach, thousands of species present in any environmental sample can be detected by high-throughput DNA sequencing and identified using public databases. Though initially limited to homogeneous substrates such as plankton or sediments, the applications of metabarcoding have been recently extended to communities on heterogeneous complex hard bottom substrates. Here we present novel metabarcoding protocols, based on the use of short fragments of 18S rRNA or cytochrome c oxidase I genes as genetic markers. We aim to develop methods for robust, reproducible eukaryotic biodiversity assessment of structurally complex communities such as marine animal forests, allowing characterization of communities living on hard-bottom substrates or other marine benthic ecosystems. We propose some guidelines focusing on sampling techniques, sample preprocessing, DNA extraction, selection of genetic markers, and bioinformatic pipelines, including steps such as sequence filtering (removal of low quality reads), clustering algorithms for delimiting molecular operational taxonomic units, and automated taxonomic assignment using reference databases.We expect these recommendations will help marine ecologists to become familiar with the paradigm shift that metabarcoding represents in the way marine ecosystems will be monitored and managed in the next future.

We assessed spatio-temporal patterns of diversity in deep-sea sediment communities using metabarcoding. We chose a recently developed eukaryotic marker based on the v7 region of the 18S rRNA gene. Our study was performed in a submarine canyon and its adjacent slope in the Northwestern Mediterranean Sea, sampled along a depth gradient at two different seasons. We found a total of 5,569 molecular operational taxonomic units (MOTUs), dominated by Metazoa, Alveolata and Rhizaria. Among metazoans, Nematoda, Arthropoda and Annelida were the most diverse. We found a marked heterogeneity at all scales, with important differences between layers of sediment and significant changes in community composition with zone (canyon vs slope), depth, and season. We compared the information obtained from metabarcoding DNA and RNA and found more total MOTUs and more MOTUs per sample with DNA (ca. 20% and 40% increase, respectively). Both datasets showed overall similar spatial trends, but most groups had higher MOTU richness with the DNA template, while others, such as nematodes, were more diverse in the RNA dataset. We provide metabarcoding protocols and guidelines for biomonitoring of these key communities in order to generate information applicable to management efforts.

Several tagging methods have been tested to assess important aspects on ecology, such as population dynamic, movement pattern, and behavioral studies. The use of external and internal tags has been extensively used for individual identification in sea urchins. Different types of tags have been used in Echinometra lucunter, considered the most common sea urchin in the Caribbean. We evaluated whether PIT-tags and S-tags may affect retention, survival and growth rates under laboratory conditions, using the sea urchin E. lucunter as a model. The type of tag was critical in terms of retention rates, with significant differences between PIT-tag and S-tags groups (Kruskal–Wallis, H = 5.33, p = 0.021), and mortality exhibited similar rates among treatments. No significant association was detected between initial maximum test diameter and the number of weeks the individuals retained neither the PIT-tag (r = −0.09, p = 0.620), nor the S-tag (r = −0.175, p = 0.413), and no significant differences in growth rate were detected among treatments (F = 1.66, df = 2, p = 0.401). We do not recommend the use of S-tag due to low retention; but we support the use of PIT-tags in E. lucunter for their positive performance displayed in terms of retention and survival but with caution, particularly in those studies requiring 80% or higher retention.

Among echinoderms, asexual reproduction by fission occurs in few species. This strategy is considered a temporary response to stressful conditions and usually alternates with sexual reproduction events; thus, monoclonal populations are extremely rare. The occurrence of a single-clone population of the starfish Coscinasterias tenuispina at Llanc¸a` (NW Mediterranean) allowed us to study intra-clonal variation of the reproductive cycle during a two-year study. The few developed gonads (all male) were found in winter months, coinciding with the minimum photoperiod (q = -0.82; P\0.001) and lowest temperatures (q = -0.75; P\0.001), only in bestfed individuals, indicating that food availability influences individual ability for gonad development. Fissiparity happened throughout all the sampled period, but its rate increased with warm temperatures (q = 0.68; P\0.0001). In contrast to what has been reported in other species, no correlation between fission rates and population density was found. The population was maintained over time by asexual reproduction and remained monoclonal. Although sexual reproduction has probably not occurred in this all-male population for a long time, the ability to yearly produce mature gonads is retained by some individuals, indicating that potential to reproduce sexually may be preserved, even in the case of strictly asexual populations.

We studied the reproductive cycle of the sea urchin Arbacia lixula in a subtidal population from northeast Spain over four years using a gonadosomatic index (GSI) and gonad histology. Our results show that the GSI of A. lixula follows a seasonal cycle which peaks in May– July and attains its lowest values in October–November every year. The time course of the GSI matched closely the photoperiod cycle. We also found a remarkable interannual variability in the maximum value of GSI, which correlated with mean water temperature during the gonad growth period (winter and spring). Gonad histology was also in agreement with a single gametogenic cycle per year in this species. We explored the application of circular statistics to present and analyse gonadal development data, which allowed us to adequately handle the high intraindividual variability detected, with several developmental stages commonly found within the same gonad. The picture that emerged is one of a gametogenic timing driven by photoperiod, while the amount of reproductive output is determined by temperature. This is coherent with the tropical origin of the species and lends support to recent warnings about an increase in the abundance of this species in the Mediterranean as a result of global warming, with associated increased impact potential in sublittoral communities.

We studied the effects of temperature and pH on larval development, settlement and juvenile survival of a Mediterranean population of the sea urchin Arbacia lixula. Three temperatures (16, 17.5 and 19 °C) were tested at present pH conditions (pHT 8.1). At 19 °C, two pH levels were compared to reflect present average (pHT 8.1) and near-future average conditions (pHT 7.7, expected by 2100). Larvae were reared for 52-days to achieve the full larval development and complete the metamorphosis to the settler stage. We analyzed larval survival, growth, morphology and settlement success. We also tested the carry-over effect of acidification on juvenile survival after 3 days. Our results showed that larval survival and size significantly increased with temperature. Acidification resulted in higher survival rates and developmental delay. Larval morphology was significantly altered by low temperatures, which led to narrower larvae with relatively shorter skeletal rods, but larval morphology was only marginally affected by acidification. No carry-over effects between larvae and juveniles were detected in early settler survival, though settlers from larvae reared at pH 7.7 were significantly smaller than their counterparts developed at pH 8.1. These results suggest an overall positive effect of environmental parameters related to global change on the reproduction of A. lixula, and reinforce the concerns about the increasing negative impact on shallow Mediterranean ecosystems of this post-glacial colonizer.

We present the global phylogeography of the black sea urchin Arbacia lixula, an amphi-Atlantic echinoid with potential to strongly impact shallow rocky ecosystems. Sequences of the mitochondrial cytochrome c oxidase gene of 604 specimens from 24 localities were obtained, covering most of the distribution area of the species, including the Mediterranean and both shores of the Atlantic. Genetic diversity measures, phylogeographic patterns, demographic parameters and population differentiation were analysed. We found high haplotype diversity but relatively low nucleotide diversity, with 176 haplotypes grouped within three haplogroups: one is shared between Eastern Atlantic (including Mediterranean) and Brazilian populations, the second is found in Eastern Atlantic and the Mediterranean and the third is exclusively from Brazil. Significant genetic differentiation was found between Brazilian, Eastern Atlantic and Mediterranean regions, but no differentiation was found among Mediterranean sub-basins or among Eastern Atlantic sub-regions. The star-shaped topology of the haplotype network and the unimodal mismatch distributions of Mediterranean and Eastern Atlantic samples suggest that these populations have suffered very recent demographic expansions. These expansions could be dated 94–205 kya in the Mediterranean, and 31–67 kya in the Eastern Atlantic. In contrast, Brazilian populations did not show any signature of population expansion. Our results indicate that all populations of A. lixula constitute a single species. The Brazilian populations probably diverged from an Eastern Atlantic stock. The present-day genetic structure of the species in Eastern Atlantic and the Mediterranean is shaped by very recent demographic processes. Our results support the view (backed by the lack of fossil record) that A. lixula is a recent thermophilous colonizer which spread throughout the Mediterranean during a warm period of the Pleistocene, probably during the last interglacial. Impli

Asymmetry between males and females in the energy they invest initially in reproduction has resulted in the evolution of differing reproductive strategies (caring females vs. competitive males). However, parental care in many birds is shared by both sexes suggesting that male energy expenditure in agonistic behaviors and courtship feeding might compensate female costs of clutch production. Here, we tested the hypothesis that initial investment in reproduction by both sexes in the Yellow-legged Gull (Larus michahellis), a species with biparental care, is similar from a physiological perspective. In this income breeder, female and male reproductive investment during early breeding can be ultimately related to muscular activity (local foraging effort required for clutch production in females and courtship feeding and agonistic behaviors in the case of males). Thus, we evaluated sex-specific patterns of creatine kinase (CK, IU/L) levels in plasma, an indicator of physical effort associated with muscular activity dependent behaviors, through incubation as a reflection of the physiological response of both sexes to the reproductive investment they made up to clutch completion. Raw levels of CK were related to plasma levels of total proteins (TP, g/dL) to account for the differential physiological state of individuals when sampled (i.e. differential dehydratation). Female costs of clutch production were associated with post-laying levels of CK/TP. We grouped females according to their relative investment in clutch production: <15.8%, 15.8 to 17.3% and >17.3% of field metabolic rate; which showed increasing values of CK/TP (24.6, 53.1 and 66.0. IU/g, respectively). Moreover, we found similar CK/TP trends throughout incubation for both sexes (CK/TP = 50.2? [3.3 × days from laying]) suggesting similar physiological responses to reproductive effort and, therefore, analogous sex-specific initial investment. Thus, male investment in agonistic behaviors and courtship feeding appa