Geographic isolation and larval dispersal shape seascape genetic patterns differently according to spatial scale

Genomic differentiation of three pico-phytoplankton species in the Mediterranean Sea

For more than a decade, high-throughput sequencing has transformed the study of marine planktonic communities and has highlighted the extent of protist diversity in these ecosystems. Nevertheless, little is known relative to their genomic diversity at the species-scale as well as their major speciation mechanisms. An increasing number of data obtained from global scale sampling campaigns is becoming publicly available, and we postulate that metagenomic data could contribute to deciphering the processes shaping protist genomic differentiation in the marine realm. As a proof of concept, we developed a findable, accessible, interoperable and reusable (FAIR) pipeline and focused on the Mediterranean Sea to study three a priori abundant protist species: Bathycoccus prasinos, Pelagomonas calceolata and Phaeocystis cordata. We compared the genomic differentiation of each species in light of geographic, environmental and oceanographic distances. We highlighted that isolation-by-environment shapes the genomic differentiation of B. prasinos, whereas P. cordata is impacted by geographic distance (i.e. isolation-by-distance). At present time, the use of metagenomics to accurately estimate the genomic differentiation of protists remains challenging since coverages are lower compared to traditional population surveys. However, our approach sheds light on ecological and evolutionary processes occurring within natural marine populations and paves the way for future protist population metagenomic studies.

Spatial coalescent connectivity through multi-generation dispersal modelling predicts gene flow across marine phyla

Gene flow governs the contemporary spatial structure and dynamic of populations as well as their long-term evolution. For species that disperse using atmospheric or oceanic flows, biophysical models allow predicting the migratory component of gene flow, which facilitates the interpretation of broad-scale spatial structure inferred from observed allele frequencies among populations. However, frequent mismatches between dispersal estimates and observed genetic diversity prevent an operational synthesis for eco-evolutionary projections. Here we use an extensive compilation of 58 population genetic studies of 47 phylogenetically divergent marine sedentary species over the Mediterranean basin to assess how genetic differentiation is predicted by Isolation-By-Distance, single-generation dispersal and multi-generation dispersal models. Unlike previous approaches, the latter unveil explicit parents-to-offspring links (filial connectivity) and implicit links among siblings from a common ancestor (coalescent connectivity). We find that almost 70 % of observed variance in genetic differentiation is explained by coalescent connectivity over multiple generations, significantly outperforming other models. Our results offer great promises to untangle the eco-evolutionary forces that shape sedentary population structure and to anticipate climate-driven redistributions, altogether improving spatial conservation planning.

Biophysical models of dispersal contribute to seascape genetic analyses

Dispersal is generally difficult to directly observe. Instead, dispersal is often inferred from genetic markers and biophysical modelling where a correspondence indicates that dispersal routes and barriers explain a significant part of population genetic differentiation. Biophysical models are used for wind-driven dispersal in terrestrial environments and for propagules drifting with ocean currents in the sea. In the ocean, such seascape genetic or seascape genomic studies provide promising tools in applied sciences, as actions within management and conservation rely on an understanding of population structure, genetic diversity and presence of local adaptations, all dependent on dispersal within the metapopulation. Here, we surveyed 87 studies that combine population genetics and biophysical models of dispersal. Our aim was to understand if biophysical dispersal models can generally explain genetic differentiation. Our analysis shows that genetic differentiation and lack of genetic differentiation can often be explained by dispersal, but the realism of the biophysical model, as well as local geomorphology and species biology also play a role. The review supports the use of a combination of both methods, and we discuss our findings in terms of recommendations for future studies and pinpoint areas where further development is necessary, particularly on how to compare both approaches. This article is part of the theme issue 'Species' ranges in the face of changing environments (part I)'.

Evidence of spatial genetic structure in a snow leopard population from Gansu, China

Understanding the spatial structure of genetic diversity provides insights into a populations' genetic status and enables assessment of its capacity to counteract the effects of genetic drift. Such knowledge is particularly scarce for the snow leopard, a conservation flagship species of Central Asia mountains. Focusing on a snow leopard population in the Qilian mountains of Gansu Province, China, we characterised the spatial genetic patterns by incorporating spatially explicit indices of diversity and multivariate analyses, based on different inertia levels of Principal Component Analysis (PCA). We compared two datasets differing in the number of loci and individuals. We found that genetic patterns were significantly spatially structured and were characterised by a broad geographical division coupled with a fine-scale cline of differentiation. Genetic admixture was detected in two adjoining core areas characterised by higher effective population size and allelic diversity, compared to peripheral localities. The power to detect significant spatial relationships depended primarily on the number of loci, and secondarily on the number of PCA axes. Spatial patterns and indices of diversity highlighted the cryptic structure of snow leopard genetic diversity, likely driven by its ability to disperse over large distances. In combination, the species' low allelic richness and large dispersal ability result in weak genetic differentiation related to major geographical features and isolation by distance. This study illustrates how cryptic genetic patterns can be investigated and analysed at a fine spatial scale, providing insights into the spatially variable isolation effects of both geographic distance and landscape resistance.

Climate differently influences the genomic patterns of two sympatric marine fish species

Climate influences population genetic variation in marine species. Capturing these impacts remains challenging for marine fishes which disperse over large geographical scales spanning steep environmental gradients. It requires the extensive spatial sampling of individuals or populations, representative of seascape heterogeneity, combined with a set of highly informative molecular markers capable of revealing climatic-associated genetic variations. We explored how space, dispersal and environment shape the genomic patterns of two sympatric fish species in the Mediterranean Sea, which ranks among the oceanic basins most affected by climate change and human pressure. We hypothesized that the population structure and climate-associated genomic signatures of selection would be stronger in the less mobile species, as restricted gene flow tends to facilitate the fixation of locally adapted alleles. To test our hypothesis, we genotyped two species with contrasting dispersal abilities: the white seabream Diplodus sargus and the striped red mullet Mullus surmuletus. We collected 823 individuals and used genotyping by sequencing (GBS) to detect 8,206 single nucleotide polymorphisms (SNPs) for the seabream and 2,794 for the mullet. For each species, we identified highly differentiated genomic regions (i.e. outliers) and disentangled the relative contribution of space, dispersal and environmental variables (climate, marine primary productivity) on the outliers' genetic structure to test the prevalence of gene flow and local adaptation. We observed contrasting patterns of gene flow and adaptive genetic variation between the two species. The seabream showed a distinct Alboran sea population and panmixia across the Mediterranean Sea. The mullet revealed additional differentiation within the Mediterranean Sea that was significantly correlated to summer and winter temperatures, as well as marine primary productivity. Functional annotation of the climate-associated outlier SNPs then identified candidate genes involved in heat tolerance that could be examined to further predict species' responses to climate change. Our results illustrate the key steps of a comparative seascape genomics study aiming to unravel the evolutionary processes at play in marine species, to better anticipate their response to climate change. Defining population adaptation capacities and environmental niches can then serve to incorporate evolutionary processes into species conservation planning.

Species ecology explains the spatial components of genetic diversity in tropical reef fishes

Generating genomic data for 19 tropical reef fish species of the Western Indian Ocean, we investigate how species ecology influences genetic diversity patterns from local to regional scales. We distinguish between the α, β and γ components of genetic diversity, which we subsequently link to six ecological traits. We find that the α and γ components of genetic diversity are strongly correlated so that species with a high total regional genetic diversity display systematically high local diversity. The α and γ diversity components are negatively associated with species abundance recorded using underwater visual surveys and positively associated with body size. Pelagic larval duration is found to be negatively related to genetic β diversity supporting its role as a dispersal trait in marine fishes. Deviation from the neutral theory of molecular evolution motivates further effort to understand the processes shaping genetic diversity and ultimately the diversification of the exceptional diversity of tropical reef fishes.

Comparative population genetics of swimming crab host (Portunus pelagicus) and common symbiotic barnacle (Octolasmis angulata) in Vietnam

Background

By comparing spatial geographical structures of host populations with that of their symbionts light can be shed on their biological interactions, and the degree of congruence between host and symbiont phylogeographies should reflect their life histories and especially dispersal mechanisms.

Methods

Here, we analyzed the genetic diversity and structure of a host, the blue swimming crab, Portunus pelagicus, and its symbiotic pedunculate barnacle Octolasmis angulata from six location sites representing three geographic regions (north, central and south) along the Vietnam coastline. High levels of congruence in their phylogeographic patterns were expected as they both undergo planktonic larval stages.

Results

Based on the COI mtDNA markers, O. angulata populations showed higher genetic diversity in comparison with their host P. pelagicus (number of haplotype/individuals, haplotype and nucleotide diversity are 119/192, 0.991 ± 0.002 and 0.02; and 89/160, 0.913 ± 0.02 and 0.015, respectively). Pairwise Fst and AMOVA analyses showed a more pronounced population structure in the symbiotic barnacle than in its crab host. The DAPC analyses identified three genetic clusters. However, both haplotype networks and scatter plots supported connectivity of the host and the symbiotic barnacle throughout their distribution range, except for low subdivision of southern population. Isolation by distance were detected only for the symbiont O. angulata (R² = 0.332, P = 0.05), while dbMEM supported spatial structure of both partners, but only at MEM-1 (Obs. 0.2686, P < 0.01 and Obs. 0.2096, P < 0.01, respectively).

Drivers for genetic structure at different geographic scales for Pacific red snapper (Lutjanus peru) and yellow snapper (Lutjanus argentiventris) in the tropical eastern Pacific

The tropical eastern Pacific (TEP) is a highly dynamic region and a model system to study how habitat discontinuities affect the distribution of shorefishes, particularly for species that display ontogenetic habitat shifts, including snappers (Lutjanidae). To evaluate the genetic structure of the Pacific red snapper (Lutjanus peru) and the yellow snapper (Lutjanus argentiventris) throughout their distribution range along the TEP, 13 and 11 microsatellite loci were analysed, respectively. The genetic diversity of L. peru (N = 446) and L. argentiventris (N = 170) was evaluated in 10 and 5 localities, respectively, showing slightly higher but non-significant values in the Gulf of California for both species. The genetic structure analysis identified the presence of significant genetic structure in both species, but the locations of the identified barriers for the gene flow differed between species. The principal driver for the genetic structure at large scales >2500 km was isolation by distance. At smaller scales (<250 km), the habitat discontinuity for juveniles and adults and the environmental differences throughout the distribution range represented potential barriers to gene flow between populations for both species.

Scale-dependent signatures of local adaptation in a foundation tree species

Understanding local adaptation is critical for conservation management under rapidly changing environmental conditions. Local adaptation inferred from genotype-environment associations may show different genomic patterns depending on the spatial scale of sampling, due to differences in the slope of environmental gradients and the level of gene flow. We compared signatures of local adaptation across the genome of mountain ash (Eucalyptus regnans) at two spatial scales: A species-wide data set and a topographically-complex subregional data set. We genotyped 367 individual trees at over 3700 single-nucleotide polymorphisms (SNPs), quantified patterns of spatial genetic structure among populations, and used two analytical methods to identify loci associated with at least one of three environmental variables at each spatial scale. Together, the analyses identified 549 potentially adaptive SNPs at the subregion scale, and 435 SNPs at the range-wide scale. A total of 39 genic or near-genic SNPs, associated with 28 genes, were identified at both spatial scales, although no SNP was identified by both methods at both scales. We observed that nongenic regions had significantly higher homozygote excess than genic regions, possibly due to selective elimination of inbred genotypes during stand development. Our results suggest that strong environmental selection occurs in mountain ash, and that the identification of putatively adaptive loci can differ substantially depending on the spatial scale of analyses. We also highlight the importance of multiple adaptive genetic architectures for understanding patterns of local adaptation across large heterogenous landscapes, with comparison of putatively adaptive loci among spatial scales providing crucial insights into the process of adaptation.

Population genomics and history of speciation reveal fishery management gaps in two related redfish species (Sebastes mentella and Sebastes fasciatus)

Understanding the processes shaping population structure and reproductive isolation of marine organisms can improve their management and conservation. Using genomic markers combined with estimation of individual ancestries, assignment tests, spatial ecology, and demographic modeling, we (i) characterized the contemporary population structure, (ii) assessed the influence of space, fishing depth, and sampling years on contemporary distribution, and (iii) reconstructed the speciation history of two cryptic redfish species, Sebastes mentella and S. fasciatus. We genotyped 860 individuals in the Northwest Atlantic Ocean using 24,603 filtered single nucleotide polymorphisms (SNPs). Our results confirmed the clear genetic distinctiveness of the two species and identified three ecotypes within S. mentella and five populations in S. fasciatus. Multivariate analyses highlighted the influence of spatial distribution and depth on the overall genomic variation, while demographic modeling revealed that secondary contact models best explained inter- and intragenomic divergence. These species, ecotypes, and populations can be considered as a rare and wide continuum of genomic divergence in the marine environment. This acquired knowledge pertaining to the evolutionary processes driving population divergence and reproductive isolation will help optimizing the assessment of demographic units and possibly to refine fishery management units.

New genomic resources for three exploited Mediterranean fishes

Extensive fishing has led to fish stock declines throughout the last decades. While clear stock identification is required for designing management schemes, stock delineation is problematic due to generally low levels of genetic structure in marine species. The development of genomic resources can help to solve this issue. Here, we present the first mitochondrial and nuclear draft genome assemblies of three economically important Mediterranean fishes, the white seabream, the striped red mullet, and the comber. The assemblies are between 613 and 785 Mbp long and contain between 27,222 and 32,375 predicted genes. They were used as references to map Restriction-site Associated DNA markers, which were developed with a single-digest approach. This approach provided between 15,710 and 21,101 Single Nucleotide Polymorphism markers per species. These genomic resources will allow uncovering subtle genetic structure, identifying stocks, assigning catches to populations and assessing connectivity. Furthermore, the annotated genomes will help to characterize adaptive divergence.

Unraveling hierarchical genetic structure in a marine metapopulation: A comparison of three high-throughput genotyping approaches

Marine metapopulations often exhibit subtle population structure that can be difficult to detect. Given recent advances in high-throughput sequencing, an emerging question is whether various genotyping approaches, in concert with improved sampling designs, will substantially improve our understanding of genetic structure in the sea. To address this question, we explored hierarchical patterns of structure in the coral reef fish Elacatinus lori using a high-resolution approach with respect to both genetic and geographic sampling. Previously, we identified three putative E. lori populations within Belize using traditional genetic markers and sparse geographic sampling: barrier reef and Turneffe Atoll; Glover's Atoll; and Lighthouse Atoll. Here, we systematically sampled individuals at ~10 km intervals throughout these reefs (1,129 individuals from 35 sites) and sequenced all individuals at three sets of markers: 2,418 SNPs; 89 microsatellites; and 57 nonrepetitive nuclear loci. At broad spatial scales, the markers were consistent with each other and with previous findings. At finer spatial scales, there was new evidence of genetic substructure, but our three marker sets differed slightly in their ability to detect these patterns. Specifically, we found subtle structure between the barrier reef and Turneffe Atoll, with SNPs resolving this pattern most effectively. We also documented isolation by distance within the barrier reef. Sensitivity analyses revealed that the number of loci (and alleles) had a strong effect on the detection of structure for all three marker sets, particularly at small spatial scales. Taken together, these results illustrate empirically that high-throughput genotyping data can elucidate subtle genetic structure at previously-undetected scales in a dispersive marine fish.

Ocean currents and the population genetic signature of fish migrations

Animal migrations are a fascinating and global phenomenon, yet they are often difficult to study and sometimes poorly understood. Here, we build on classic ecological theory by hypothesizing that some enigmatic spawning migrations across coastal marine habitats can be inferred from the population genetic signature of larval dispersal by ocean currents. We test this assumption by integrating spatially realistic simulations of alternative spawning migration routes, associated patterns of larval dispersal, and associated variation in the population genetic structure of eastern Australian sea mullet (Mugil cephalus). We then use simulation results to assess the implications of alternative spawning destinations for larval replenishment, and we contrast simulated against measured population genetic variation. Both analyses suggest that the spawning migrations of M. cephalus in eastern Australia are likely to be localized (approximately 100 km along the shore), and that spawning is likely to occur in inshore waters. Our conclusions are supported by multiple lines of evidence available through independent studies, but they challenge the more traditional assumption of a single, long-distance migration event with subsequent offshore spawning in the East Australian Current. More generally, our study operationalizes classic theory on the relationship between fish migrations, ocean currents, and reproductive success. However, rather than confirming the traditionally assumed adaptation of migratory behavior to dominant ocean current flow, our findings support the concept of a genetically measurable link between fish migrations and local oceanographic conditions, specifically water temperature and coastal retention of larvae. We believe that future studies using similar approaches for high resolution and spatially realistic ecological-genetic scenario testing can help rapidly advance our understanding of key ecological processes in many other marine species.