Genetic and phenotypic differentiation between the critically endangered Balearic shearwater and neighboring colonies of its sibling species

Species delimitation using genomic data to resolve taxonomic uncertainties in a speciation continuum of pelagic seabirds

Speciation is a continuous and complex process shaped by the interaction of numerous evolutionary forces. Despite the continuous nature of the speciation process, the implementation of conservation policies relies on the delimitation of species and evolutionary significant units (ESUs). Puffinus shearwaters are globally distributed and threatened pelagic seabirds. Due to remarkable morphological status the group has been under intense taxonomic debate for the past three decades. Here, we use double digest Restriction-Site Associated DNA sequencing (ddRAD-Seq) to genotype species and subspecies of North Atlantic and Mediterranean Puffinus shearwaters across their entire geographical range. We assess the phylogenetic relationships and population structure among and within the group, evaluate species boundaries, and characterise the genomic landscape of divergence. We find that current taxonomies are not supported by genomic data and propose a more accurate taxonomy by integrating genomic information with other sources of evidence. Our results show that several taxon pairs are at different stages of a speciation continuum. Our study emphasises the potential of genomic data to resolve taxonomic uncertainties, which can help to focus management actions on relevant taxa, even if they do not necessarily coincide with the taxonomic rank of species.

The First Genome of the Balearic Shearwater (Puffinus mauretanicus) Provides a Valuable Resource for Conservation Genomics and Sheds Light on Adaptation to a Pelagic lifestyle

The Balearic shearwater (Puffinus mauretanicus) is the most threatened seabird in Europe and a member of the most speciose group of pelagic seabirds, the order Procellariiformes, which exhibit extreme adaptations to a pelagic lifestyle. The fossil record suggests that human colonisation of the Balearic Islands resulted in a sharp decrease of the Balearic shearwater population size. Currently, populations of the species continue to be decimated mainly due to predation by introduced mammals and bycatch in longline fisheries, with some studies predicting its extinction by 2070. Here, using a combination of short and long reads, we generate the first high-quality reference genome for the Balearic shearwater, with a completeness amongst the highest across available avian species. We used this reference genome to study critical aspects relevant to the conservation status of the species and to gain insights into the adaptation to a pelagic lifestyle of the order Procellariiformes. We detected relatively high levels of genome-wide heterozygosity in the Balearic shearwater despite its reduced population size. However, the reconstruction of its historical demography uncovered an abrupt population decline potentially linked to a reduction of the neritic zone during the Penultimate Glacial Period (∼194-135 ka). Comparative genomics analyses uncover a set of candidate genes that may have played an important role into the adaptation to a pelagic lifestyle of Procellariiformes, including those for the enhancement of fishing capabilities, night vision, and the development of natriuresis. The reference genome obtained will be the crucial in the future development of genetic tools in conservation efforts for this Critically Endangered species.

Sea surface temperature, rather than land mass or geographic distance, may drive genetic differentiation in a species complex of highly dispersive seabirds

Seabirds, particularly Procellariiformes, are highly mobile organisms with a great capacity for long dispersal, though simultaneously showing high philopatry, two conflicting life-history traits that may lead to contrasted patterns of genetic population structure. Landmasses were suggested to explain differentiation patterns observed in seabirds, but philopatry, isolation by distance, segregation between breeding and nonbreeding zones, and oceanographic conditions (sea surface temperatures) may also contribute to differentiation patterns. To our knowledge, no study has simultaneously contrasted the multiple factors contributing to the diversification of seabird species, especially in the gray zone of speciation. We conducted a multilocus phylogeographic study on a widespread seabird species complex, the little shearwater complex, showing highly homogeneous morphology, which led to considerable taxonomic debate. We sequenced three mitochondrial and six nuclear markers on all extant populations from the Atlantic (lherminieri) and Indian Oceans (bailloni), that is, five nominal lineages from 13 populations, along with one population from the eastern Pacific Ocean (representing the dichrous lineage). We found sharp differentiation among populations separated by the African continent with both mitochondrial and nuclear markers, while only mitochondrial markers allowed characterizing the five nominal lineages. No differentiation could be detected within these five lineages, questioning the strong level of philopatry showed by these shearwaters. Finally, we propose that Atlantic populations likely originated from the Indian Ocean. Within the Atlantic, a stepping-stone process accounts for the current distribution. Based on our divergence time estimates, we suggest that the observed pattern of differentiation mostly resulted from historical and current variation in sea surface temperatures.

Integrating Sequence Capture and Restriction Site-Associated DNA Sequencing to Resolve Recent Radiations of Pelagic Seabirds

The diversification of modern birds has been shaped by a number of radiations. Rapid diversification events make reconstructing the evolutionary relationships among taxa challenging due to the convoluted effects of incomplete lineage sorting (ILS) and introgression. Phylogenomic data sets have the potential to detect patterns of phylogenetic incongruence, and to address their causes. However, the footprints of ILS and introgression on sequence data can vary between different phylogenomic markers at different phylogenetic scales depending on factors such as their evolutionary rates or their selection pressures. We show that combining phylogenomic markers that evolve at different rates, such as paired-end double-digest restriction site-associated DNA (PE-ddRAD) and ultraconserved elements (UCEs), allows a comprehensive exploration of the causes of phylogenetic discordance associated with short internodes at different timescales. We used thousands of UCE and PE-ddRAD markers to produce the first well-resolved phylogeny of shearwaters, a group of medium-sized pelagic seabirds that are among the most phylogenetically controversial and endangered bird groups. We found that phylogenomic conflict was mainly derived from high levels of ILS due to rapid speciation events. We also documented a case of introgression, despite the high philopatry of shearwaters to their breeding sites, which typically limits gene flow. We integrated state-of-the-art concatenated and coalescent-based approaches to expand on previous comparisons of UCE and RAD-Seq data sets for phylogenetics, divergence time estimation, and inference of introgression, and we propose a strategy to optimize RAD-Seq data for phylogenetic analyses. Our results highlight the usefulness of combining phylogenomic markers evolving at different rates to understand the causes of phylogenetic discordance at different timescales. [Aves; incomplete lineage sorting; introgression; PE-ddRAD-Seq; phylogenomics; radiations; shearwaters; UCEs.].

Identifying mechanisms of genetic differentiation among populations in vagile species: historical factors dominate genetic differentiation in seabirds

Elucidating the factors underlying the origin and maintenance of genetic variation among populations is crucial for our understanding of their ecology and evolution, and also to help identify conservation priorities. While intrinsic movement has been hypothesized as the major determinant of population genetic structuring in abundant vagile species, growing evidence indicates that vagility does not always predict genetic differentiation. However, identifying the determinants of genetic structuring can be challenging, and these are largely unknown for most vagile species. Although, in principle, levels of gene flow can be inferred from neutral allele frequency divergence among populations, underlying assumptions may be unrealistic. Moreover, molecular studies have suggested that contemporary gene flow has often not overridden historical influences on population genetic structure, which indicates potential inadequacies of any interpretations that fail to consider the influence of history in shaping that structure. This exhaustive review of the theoretical and empirical literature investigates the determinants of population genetic differentiation using seabirds as a model system for vagile taxa. Seabirds provide a tractable group within which to identify the determinants of genetic differentiation, given their widespread distribution in marine habitats and an abundance of ecological and genetic studies conducted on this group. Herein we evaluate mitochondrial DNA (mtDNA) variation in 73 seabird species. Lack of mutation-drift equilibrium observed in 19% of species coincided with lower estimates of genetic differentiation, suggesting that dynamic demographic histories can often lead to erroneous interpretations of contemporary gene flow, even in vagile species. Presence of land across the species sampling range, or sampling of breeding colonies representing ice-free Pleistocene refuge zones, appear to be associated with genetic differentiation in Tropical and Southern Temperate species, respectively, indicating that long-term barriers and persistence of populations are important for their genetic structuring. Conversely, biotic factors commonly considered to influence population genetic structure, such as spatial segregation during foraging, were inconsistently associated with population genetic differentiation. In light of these results, we recommend that genetic studies should consider potential historical events when identifying determinants of genetic differentiation among populations to avoid overestimating the role of contemporary factors, even for highly vagile taxa.

Patterns of at-sea behaviour at a hybrid zone between two threatened seabirds

Patterns of behavioural variation and migratory connectivity are important characteristics of populations, particularly at the edges of species distributions, where processes involved in influencing evolutionary trajectories, such as divergence, mutual persistence, and natural hybridization, can occur. Here, we focused on two closely related seabird species that breed in the Mediterranean: Balearic shearwaters (Puffinus mauretanicus) and Yelkouan shearwaters (Puffinus yelkouan). Genetic and phenotypic evidence of hybridization between the two species on Menorca (the eastern and westernmost island in the breeding ranges of the two shearwaters, respectively) has provided important insights into relationships between these recently diverged species. Nevertheless, levels of behavioural and ecological differentiation amongst these populations remain largely unknown. Using geolocation and stable isotopes, we compared the at-sea movement behaviour of birds from the Menorcan ‘hybrid’ population with the nearest neighbouring populations of Balearic and Yelkouan shearwaters. The Menorcan population displayed a suite of behavioural features intermediate to those seen in the two species (including migration strategies, breeding season movements and limited data on phenology). Our findings provide new evidence to support suggestions that the Menorcan population is admixed, and indicate a role of non-breeding behaviours in the evolutionary trajectories of Puffinus shearwaters in the Mediterranean.

Comparing multiple criteria for species identification in two recently diverged seabirds

Correct species identification is a crucial issue in systematics with key implications for prioritising conservation effort. However, it can be particularly challenging in recently diverged species due to their strong similarity and relatedness. In such cases, species identification requires multiple and integrative approaches. In this study we used multiple criteria, namely plumage colouration, biometric measurements, geometric morphometrics, stable isotopes analysis (SIA) and genetics (mtDNA), to identify the species of 107 bycatch birds from two closely related seabird species, the Balearic (Puffinus mauretanicus) and Yelkouan (P. yelkouan) shearwaters. Biometric measurements, stable isotopes and genetic data produced two stable clusters of bycatch birds matching the two study species, as indicated by reference birds of known origin. Geometric morphometrics was excluded as a species identification criterion since the two clusters were not stable. The combination of plumage colouration, linear biometrics, stable isotope and genetic criteria was crucial to infer the species of 103 of the bycatch specimens. In the present study, particularly SIA emerged as a powerful criterion for species identification, but temporal stability of the isotopic values is critical for this purpose. Indeed, we found some variability in stable isotope values over the years within each species, but species differences explained most of the variance in the isotopic data. Yet this result pinpoints the importance of examining sources of variability in the isotopic data in a case-by-case basis prior to the cross-application of the SIA approach to other species. Our findings illustrate how the integration of several methodological approaches can help to correctly identify individuals from recently diverged species, as each criterion measures different biological phenomena and species divergence is not expressed simultaneously in all biological traits.

Population structure and dispersal patterns within and between Atlantic and Mediterranean populations of a large-range pelagic seabird

Dispersal is critically linked to the demographic and evolutionary trajectories of populations, but in most seabird species it may be difficult to estimate. Using molecular tools, we explored population structure and the spatial dispersal pattern of a highly pelagic but philopatric seabird, the Cory's shearwater Calonectris diomedea. Microsatellite fragments were analysed from samples collected across almost the entire breeding range of the species. To help disentangle the taxonomic status of the two subspecies described, the Atlantic form C. d. borealis and the Mediterranean form C. d. diomedea, we analysed genetic divergence between subspecies and quantified both historical and recent migration rates between the Mediterranean and Atlantic basins. We also searched for evidence of isolation by distance (IBD) and addressed spatial patterns of gene flow. We found a low genetic structure in the Mediterranean basin. Conversely, strong genetic differentiation appeared in the Atlantic basin. Even if the species was mostly philopatric (97%), results suggest recent dispersal between basins, especially from the Atlantic to the Mediterranean (aprox. 10% of migrants/generation across the last two generations). Long-term gene flow analyses also suggested an historical exchange between basins (about 70 breeders/generation). Spatial analysis of genetic variation indicates that distance is not the main factor in shaping genetic structure in this species. Given our results we recommend gathering more data before concluded whether these taxa should be treated as two species or subspecies.