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

RAD-sequencing improves the genetic characterization of a threatened tree peony (Paeonia ludlowii) endemic to China: Implications for conservation

Compared with traditional genetic markers, genomic approaches have proved valuable to the conservation of endangered species. Paeonia ludlowii having rarely and pure yellow flowers, is one of the world's most famous tree peonies. However, only several wild populations remain in the Yarlung Zangbo Valley (Nyingchi and Shannan regions, Xizang) in China due to increasing anthropogenic impact on the natural habitats. We used genome-wide single nucleotide polymorphisms to elucidate the spatial pattern of genetic variation, population structure and demographic history of P. ludlowii from the fragmented region comprising the entire range of this species, aiming to provide a basis for conserving the genetic resources of this species. Unlike genetic uniformity among populations revealed in previous studies, we found low but varied levels of intra-population genetic diversity, in which lower genetic diversity was detected in the population in Shannan region compared to those in Nyingzhi region. These spatial patterns may be likely associated with different population sizes caused by micro-environment differences in these two regions. Additionally, low genetic differentiation among populations (Fst = 0.0037) were detected at the species level. This line of evidence, combined with the result of significant genetic differentiation between the two closest populations and lack of isolation by distance, suggested that shared ancestry among now remnant populations rather than contemporary genetic connectivity resulted in subtle population structure. Demographic inference suggested that P. ludlowii probably experienced a temporal history of sharp population decline during the period of Last Glacial Maximum, and a subsequent bottleneck event resulting from prehistoric human activities on the Qinghai-Tibet Plateau. All these events, together with current habitat fragment and excavation might contribute to the endangered status of P. ludlowii. Our study improved the genetic characterization of the endangered tree peony (P. ludlowii) in China, and these genetic inferences should be considered when making different in situ and ex situ conservation actions for P. ludlowii in this evolutionary hotspot region.

Sponge diversification in marine lakes: Implications for phylogeography and population genomic studies on sponges

The relative influence of geography, currents, and environment on gene flow within sessile marine species remains an open question. Detecting subtle genetic differentiation at small scales is challenging in benthic populations due to large effective population sizes, general lack of resolution in genetic markers, and because barriers to dispersal often remain elusive. Marine lakes can circumvent confounding factors by providing discrete and replicated ecosystems. Using high-resolution double digest restriction-site-associated DNA sequencing (4826 Single Nucleotide Polymorphisms, SNPs), we genotyped populations of the sponge Suberites diversicolor (n = 125) to test the relative importance of spatial scales (1-1400 km), local environmental conditions, and permeability of seascape barriers in shaping population genomic structure. With the SNP dataset, we show strong intralineage population structure, even at scales <10 km (average F ST = 0.63), which was not detected previously using single markers. Most variation was explained by differentiation between populations (AMOVA: 48.8%) with signatures of population size declines and bottlenecks per lake. Although the populations were strongly structured, we did not detect significant effects of geographic distance, local environments, or degree of connection to the sea on population structure, suggesting mechanisms such as founder events with subsequent priority effects may be at play. We show that the inclusion of morphologically cryptic lineages that can be detected with the COI marker can reduce the obtained SNP set by around 90%. Future work on sponge genomics should confirm that only one lineage is included. Our results call for a reassessment of poorly dispersing benthic organisms that were previously assumed to be highly connected based on low-resolution markers.

Molecular ecology meets systematic conservation planning

Integrative and proactive conservation approaches are critical to the long-term persistence of biodiversity. Molecular data can provide important information on evolutionary processes necessary for conserving multiple levels of biodiversity (genes, populations, species, and ecosystems). However, molecular data are rarely used to guide spatial conservation decision-making. Here, we bridge the fields of molecular ecology (ME) and systematic conservation planning (SCP) (the 'why') to build a foundation for the inclusion of molecular data into spatial conservation planning tools (the 'how'), and provide a practical guide for implementing this integrative approach for both conservation planners and molecular ecologists. The proposed framework enhances interdisciplinary capacity, which is crucial to achieving the ambitious global conservation goals envisioned for the next decade.

Landscape genomics of the American lobster (Homarus americanus)

In marine species experiencing intense fishing pressures, knowledge of genetic structure and local adaptation represent a critical information to assist sustainable management. In this study, we performed a landscape genomics analysis in the American lobster to investigate the issues pertaining to the consequences of making use of putative adaptive loci to reliably infer population structure and thus more rigorously delineating biological management units in marine exploited species. Toward this end, we genotyped 14,893 single nucleotide polymorphism (SNPs) in 4190 lobsters sampled across 96 sampling sites distributed along 1000 km in the northwest Atlantic in both Canada and the USA. As typical for most marine species, we observed a weak, albeit highly significant genetic structure. We also found that adaptive genetic variation allows detecting fine-scale population structure not resolved by neutral genetic variation alone. Using the recent genome assembly of the American lobster, we were able to map and annotate several SNPs located in functional genes potentially implicated in adaptive processes such as thermal stress response, salinity tolerance and growth metabolism pathways. Taken together, our study indicates that weak population structure in high gene flow systems can be resolved at various spatial scales, and that putatively adaptive genetic variation can substantially enhance the delineation of biological management units of marine exploited species.

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.

Population assignment tests uncover rare long-distance marine larval dispersal events

Long-distance dispersal (LDD) is consequential to metapopulation ecology and evolution. In systems where dispersal is undertaken by small propagules, such as larvae in the ocean, documenting LDD is especially challenging. Genetic parentage analysis has gained traction as a method for measuring larval dispersal, but such studies are generally spatially limited, leaving LDD understudied in marine species. We addressed this knowledge gap by uncovering LDD with population assignment tests in the coral reef fish Elacatinus lori, a species whose short-distance dispersal has been well-characterized by parentage analysis. When adults (n = 931) collected throughout the species' range were categorized into three source populations, assignment accuracy exceeded 99%, demonstrating low rates of gene flow between populations in the adult generation. After establishing high assignment confidence, we assigned settlers (n = 3,828) to source populations. Within the settler cohort, <0.1% of individuals were identified as long-distance dispersers from other populations. These results demonstrate an exceptionally low level of connectivity between E. lori populations, despite the potential for ocean currents to facilitate LDD. More broadly, these findings illustrate the value of combining genetic parentage analysis and population assignment tests to uncover short- and long-distance dispersal, respectively.

Seascape Genomics of the Sugar Kelp Saccharina latissima along the North Eastern Atlantic Latitudinal Gradient

Temperature is one of the most important range-limiting factors for many seaweeds. Driven by the recent climatic changes, rapid northward shifts of species' distribution ranges can potentially modify the phylogeographic signature of Last Glacial Maximum. We explored this question in detail in the cold-tolerant kelp species Saccharina latissima, using microsatellites and double digest restriction site-associated DNA sequencing ( ddRAD-seq) derived single nucleotide polymorphisms (SNPs) to analyze the genetic diversity and structure in 11 sites spanning the entire European Atlantic latitudinal range of this species. In addition, we checked for statistical correlation between genetic marker allele frequencies and three environmental proxies (sea surface temperature, salinity, and water turbidity). Our findings revealed that genetic diversity was significantly higher for the northernmost locality (Spitsbergen) compared to the southern ones (Northern Iberia), which we discuss in light of the current state of knowledge on phylogeography of S. latissima and the potential influence of the recent climatic changes on the population structure of this species. Seven SNPs and 12 microsatellite alleles were found to be significantly associated with at least one of the three environmental variables. We speculate on the putative adaptive functions of the genes associated with the outlier markers and the importance of these markers for successful conservation and aquaculture strategies for S. latissima in this age of rapid global change.

Adaptive divergence, neutral panmixia, and algal symbiont population structure in the temperate coral Astrangia poculata along the Mid-Atlantic United States

Astrangia poculata is a temperate scleractinian coral that exists in facultative symbiosis with the dinoflagellate alga Breviolum psygmophilum across a range spanning the Gulf of Mexico to Cape Cod, Massachusetts. Our previous work on metabolic thermal performance of Virginia (VA) and Rhode Island (RI) populations of A. poculata revealed physiological signatures of cold (RI) and warm (VA) adaptation of these populations to their respective local thermal environments. Here, we used whole-transcriptome sequencing (mRNA-Seq) to evaluate genetic differences and identify potential loci involved in the adaptive signature of VA and RI populations. Sequencing data from 40 A. poculata individuals, including 10 colonies from each population and symbiotic state (VA-white, VA-brown, RI-white, and RI-brown), yielded a total of 1,808 host-associated and 59 algal symbiont-associated single nucleotide polymorphisms (SNPs) post filtration. Fst outlier analysis identified 66 putative high outlier SNPs in the coral host and 4 in the algal symbiont. Differentiation of VA and RI populations in the coral host was driven by putatively adaptive loci, not neutral divergence (Fst = 0.16, p = 0.001 and Fst = 0.002, p = 0.269 for outlier and neutral SNPs respectively). In contrast, we found evidence of neutral population differentiation in B. psygmophilum (Fst = 0.093, p = 0.001). Several putatively adaptive host loci occur on genes previously associated with the coral stress response. In the symbiont, three of four putatively adaptive loci are associated with photosystem proteins. The opposing pattern of neutral differentiation in B. psygmophilum, but not the A. poculata host, reflects the contrasting dynamics of coral host and algal symbiont population connectivity, dispersal, and gene by environment interactions.

Complementarity and discriminatory power of genotype and otolith shape in describing the fine-scale population structure of an exploited fish, the common sole of the Eastern English Channel

Marine organisms show population structure at a relatively fine spatial scale, even in open habitats. The tools commonly used to assess subtle patterns of connectivity have diverse levels of resolution and can complement each other to inform on population structure. We assessed and compared the discriminatory power of genetic markers and otolith shape to reveal the population structure on evolutionary and ecological time scales of the common sole (Solea solea), living in the Eastern English Channel (EEC) stock off France and the UK. First, we genotyped fish with Single Nucleotide Polymorphisms to assess population structure at an evolutionary scale. Then, we tested for spatial segregation of the subunits using otolith shape as an integrative tracer of life history. Finally, a supervised machine learning framework was applied to genotypes and otolith phenotypes to probabilistically assign adults to subunits and assess the discriminatory power of each approach. Low but significant genetic differentiation was found among subunits. Moreover, otolith shape appeared to vary spatially, suggesting spatial population structure at fine spatial scale. However, results of the supervised discriminant analyses failed to discriminate among subunits, especially for otolith shape. We suggest that the degree of population segregation may not be strong enough to allow for robust fish assignments. Finally, this study revealed a weak yet existing metapopulation structure of common sole at the fine spatial scale of the EEC based on genotypes and otolith shape, with one subunit being more isolated. Our study argues for the use of complementary tracers to investigate marine population structure.

Natal philopatry increases relatedness within groups of coral reef cardinalfish

A central issue in evolutionary ecology is how patterns of dispersal influence patterns of relatedness in populations. In terrestrial organisms, limited dispersal of offspring leads to groups of related individuals. By contrast, for most marine organisms, larval dispersal in open waters is thought to minimize kin associations within populations. However, recent molecular evidence and theoretical approaches have shown that limited dispersal, sibling cohesion and/or differential reproductive success can lead to kin association and elevated relatedness. Here, we tested the hypothesis that limited dispersal explains small-scale patterns of relatedness in the pajama cardinalfish Sphaeramia nematoptera. We used 19 microsatellite markers to assess parentage of 233 juveniles and pairwise relatedness among 527 individuals from 41 groups in Kimbe Bay, Papua New Guinea. Our findings support three predictions of the limited dispersal hypothesis: (i) elevated relatedness within groups, compared with among groups and elevated relatedness within reefs compared with among reefs; (ii) a weak negative correlation of relatedness with distance; (iii) more juveniles than would be expected by chance in the same group and the same reef as their parents. We provide the first example for natal philopatry at the group level causing small-scale patterns of genetic relatedness in a marine fish.

A fair fight between molecular marker types in a seascape genetics setting

From its inception, population genetics has been nearly as concerned with the genetic data type-to which analyses are brought to bear-as it is with the analysis methods themselves. The field has traversed allozymes, microsatellites, segregating sites in multilocus alignments and, currently, single nucleotide polymorphisms (SNPs) generated by high-throughput genomic sequencing methods, primarily whole genome sequencing and reduced representation library (RRL) sequencing. As each emerging data type has gained traction, it has been compared to existing methods, based on its relative ability to discern population structural complexity at increasing levels of resolution. However, this is usually done by comparing the gold standard in one data type to the gold standard in the new data type. These gold standards frequently differ in power and in sampling density, both across a genome and throughout a spatial range. In this issue of Molecular Ecology, D'Aloia et al. apply the high-throughput approach as fully as possible to microsatellites, nuclear loci and SNPs genotyped through an RRL method; this is coupled with a spatially dense sampling scheme. Completing a battery of population genetics analyses across data types (including a series of down-sampled data sets), the authors find that SNP data are slightly more sensitive to fine-scale genetic structure, and the results are more resilient to down-sampling than microsatellites and nonrepetitive nuclear loci. However, their results are far from an unqualified victory for RRL SNP data over all previous data types: the authors note that modest additions to the microsatellites and nuclear loci data sets may provide the necessary analytical power to delineate the fine-scale genetic structuring identified by SNPs. As always, as the field begins to fully embrace the newest thing, good science reminds us that traditional data types are far from useless, especially when combined with a well-designed sampling scheme.