Linking genetic kinship and demographic analyses to characterize dispersal:Methods and application to Blanding's turtle

Characterizing how frequently, and at what life stages and spatial scales, dispersal occurs can be difficult, especially for species with cryptic juvenile periods and long reproductive life spans. Using a combination of mark-recapture information, microsatellite genetic data, and demographic simulations, we characterize natal and breeding dispersal patterns in the long-lived, slow-maturing, and endangered Blanding's turtle (Emydoidea blandingii), focusing on nesting females. We captured and genotyped 310 individual Blanding's turtles (including 220 nesting females) in a central Wisconsin population from 2010 to 2013, with additional information on movements among 3 focal nesting areas within this population available from carapace-marking conducted from 2001 to 2009. Mark-recapture analyses indicated that dispersal among the 3 focal nesting areas was infrequent (<0.03 annual probability). Dyads of females with inferred first-order relationships were more likely to be found within the same nesting area than split between areas, and the proportion of related dyads declined with increasing distance among nesting areas. The observed distribution of related dyads for nesting females was consistent with a probability of natal dispersal at first breeding between nearby nesting areas of approximately 0.1 based on demographic simulations. Our simulation-based estimates of infrequent female dispersal were corroborated by significant spatial genetic autocorrelation among nesting females at scales of <500 m. Nevertheless, a lack of spatial genetic autocorrelation among non-nesting turtles (males and females) suggested extensive local connectivity, possibly mediated by male movements or long-distance movements made by females between terrestrial nesting areas and aquatic habitats. We show here that coupling genetic and demographic information with simulations of individual-based population models can be an effective approach for untangling the contributions of natal and breeding dispersal to spatial ecology

Mind the gut:Genomic insights to population divergence and gut microbial composition of two marine keystone species

BACKGROUND: Deciphering the mechanisms governing population genetic divergence and local adaptation across heterogeneous environments is a central theme in marine ecology and conservation. While population divergence and ecological adaptive potential are classically viewed at the genetic level, it has recently been argued that their microbiomes may also contribute to population genetic divergence. We explored whether this might be plausible along the well-described environmental gradient of the Baltic Sea in two species of sand lance (Ammodytes tobianus and Hyperoplus lanceolatus). Specifically, we assessed both their population genetic and gut microbial composition variation and investigated not only which environmental parameters correlate with the observed variation, but whether host genome also correlates with microbiome variation. RESULTS: We found a clear genetic structure separating the high-salinity North Sea from the low-salinity Baltic Sea sand lances. The observed genetic divergence was not simply a function of isolation by distance, but correlated with environmental parameters, such as salinity, sea surface temperature, and, in the case of A. tobianus, possibly water microbiota. Furthermore, we detected two distinct genetic groups in Baltic A. tobianus that might represent sympatric spawning types. Investigation of possible drivers of gut microbiome composition variation revealed that host species identity was significantly correlated with the microbial community composition of the gut. A potential influence of host genetic factors on gut microbiome composition was further confirmed by the results of a constrained analysis of principal coordinates. The host genetic component was among the parameters that best explain observed variation in gut microbiome composition. CONCLUSIONS: Our findings have relevance for the population structure of two commercial species but also provide insights into potentially relevant genomic and microbial factors with regards to sand lance adaptation across the North Sea-Baltic Sea environmental gradient. Furthermore, our findings support the hypothesis that host genetics may play a role in regulating the gut microbiome at both the interspecific and intraspecific levels. As sequencing costs continue to drop, we anticipate that future studies that include full genome and microbiome sequencing will be able to explore the full relationship and its potential adaptive implications for these species

Demographic changes in Pleistocene sea turtles were driven by past sea level fluctuations affecting feeding habitat availability

Pleistocene environmental changes are generally assumed to have dramatically affected species’ demography via changes in habitat availability, but this is challenging to investigate due to our limited knowledge of how Pleistocene ecosystems changed through time. Here, we tracked changes in shallow marine habitat availability resulting from Pleistocene sea level fluctuations throughout the last glacial cycle (120–14 thousand years ago; kya) and assessed correlations with past changes in genetic diversity inferred from genome-wide SNPs, obtained via ddRAD sequencing, in Caribbean hawksbill turtles, which feed in coral reefs commonly found in shallow tropical waters. We found sea level regression resulted in an average 75% reduction in shallow marine habitat availability during the last glacial cycle. Changes in shallow marine habitat availability correlated strongly with past changes in hawksbill turtle genetic diversity, which gradually declined to ~1/4th of present-day levels during the Last Glacial Maximum (LGM; 26–19 kya). Shallow marine habitat availability and genetic diversity rapidly increased after the LGM, signifying a population expansion in response to warming environmental conditions. Our results suggest a positive correlation between Pleistocene environmental changes, habitat availability and species’ demography, and that demographic changes in hawksbill turtles were potentially driven by feeding habitat availability. However, we also identified challenges associated with disentangling the potential environmental drivers of past demographic changes, which highlights the need for integrative approaches. Our conclusions underline the role of habitat availability on species’ demography and biodiversity, and that the consequences of ongoing habitat loss should not be underestimated

The era of reference genomes in conservation genomics

Progress in genome sequencing now enables the large-scale generation of reference genomes. Various international initiatives aim to generate reference genomes representing global biodiversity. These genomes provide unique insights into genomic diversity and architecture, thereby enabling comprehensive analyses of population and functional genomics, and are expected to revolutionize conservation genomics

Population Structure of Humpback Whales from Their Breeding Grounds in the South Atlantic and Indian Oceans

Although humpback whales are among the best-studied of the large whales, population boundaries in the Southern Hemisphere (SH) have remained largely untested. We assess population structure of SH humpback whales using 1,527 samples collected from whales at fourteen sampling sites within the Southwestern and Southeastern Atlantic, the Southwestern Indian Ocean, and Northern Indian Ocean (Breeding Stocks A, B, C and X, respectively). Evaluation of mtDNA population structure and migration rates was carried out under different statistical frameworks. Using all genetic evidence, the results suggest significant degrees of population structure between all ocean basins, with the Southwestern and Northern Indian Ocean most differentiated from each other. Effective migration rates were highest between the Southeastern Atlantic and the Southwestern Indian Ocean, followed by rates within the Southeastern Atlantic, and the lowest between the Southwestern and Northern Indian Ocean. At finer scales, very low gene flow was detected between the two neighbouring sub-regions in the Southeastern Atlantic, compared to high gene flow for whales within the Southwestern Indian Ocean. Our genetic results support the current management designations proposed by the International Whaling Commission of Breeding Stocks A, B, C, and X as four strongly structured populations. The population structure patterns found in this study are likely to have been influenced by a combination of long-term maternally directed fidelity of migratory destinations, along with other ecological and oceanographic features in the region

Mitochondrial Control Region and microsatellite analyses on harbour porpoise (Phocoena phocoena) unravel population differentiation in the Baltic Sea and adjacent waters

The population status of the harbour porpoise (Phocoena phocoena) in the Baltic area has been a continuous matter of debate. Here we present the by far most comprehensive genetic population structure assessment to date for this region, both with regard to geographic coverage and sample size: 497 porpoise samples from North Sea, Skagerrak, Kattegat, Belt Sea, and Inner Baltic Sea were sequenced at the mitochondrial Control Region and 305 of these specimens were typed at 15 polymorphic microsatellite loci. Samples were stratified according to sample type (stranding vs. by-caught), sex, and season (breeding vs. non-breeding season). Our data provide ample evidence for a population split between the Skagerrak and the Belt Sea, with a transition zone in the Kattegat area. Among other measures, this was particularly visible in significant frequency shifts of the most abundant mitochondrial haplotypes. A particular haplotype almost absent in the North Sea was the most abundant in Belt Sea and Inner Baltic Sea. Microsatellites yielded a similar pattern (i.e., turnover in occurrence of clusters identified by STRUCTURE). Moreover, a highly significant association between microsatellite assignment and unlinked mitochondrial haplotypes further indicates a split between North Sea and Baltic porpoises. For the Inner Baltic Sea, we consistently recovered a small, but significant separation from the Belt Sea population. Despite recent arguments that separation should exceed a predefined threshold before populations shall be managed separately, we argue in favour of precautionary acknowledging the Inner Baltic porpoises as a separate management unit, which should receive particular attention, as it is threatened by various factors, in particular local fishery measures. © Springer Science+Business Media B.V. 2009

Decline in abundance and apparent survival rates of fin whales (Balaenoptera physalus) in the northern Gulf of St. Lawrence

A. Schleimer was supported by the Luxembourg National Research Fund (FNR; AFR/11256673).Estimates of abundance and survivorship provide quantifiable measures to monitor populations and to define and understand their conservation status. This study investigated changes in abundance and survival rates of fin whales (Balaenoptera physalus) in the northern Gulf of St. Lawrence in the context of anthropogenic pressures and changing environmental conditions. A long‐term data set, consisting of 35 years of photo‐identification surveys and comprising more than 5,000 identifications of 507 individuals, formed the basis of this mark–recapture study. Based on model selection using corrected Akaike Information Criterion, the most parsimonious Cormack–Jolly–Seber model included a linear temporal trend in noncalf apparent survival rates with a sharp decline in the last 5 years of the study and a median survival rate of 0.946 (95% confidence interval (CI) 0.910–0.967). To account for capture heterogeneity due to divergent patterns of site fidelity, agglomerative hierarchical cluster analysis was employed to categorize individuals based on their annual and survey site fidelity indices. However, the negative trend in survivorship remained and was corroborated by a significant decline in the estimated super‐population size from 335 (95% CI 321–348) individuals in 2004–2010 to 291 (95% CI 270–312) individuals in 2010–2016. Concurrently, a negative trend was estimated in recruitment to the population, supported by a sharp decrease in the number of observed calves. Ship strikes and changes in prey availability are potential drivers of the observed decline in fin whale abundance. The combination of clustering methods with mark–recapture represents a flexible way to investigate the effects of site fidelity on demographic variables and is broadly applicable to other individual‐based studies.Publisher PDFPeer reviewe

Fin whale (Balaenoptera physalus) mitogenomics: A cautionary tale of defining sub-species from mitochondrial sequence monophyly

The advent of massive parallel sequencing technologies has resulted in an increase of studies based upon complete mitochondrial genome DNA sequences that revisit the taxonomic status within and among species. Spatially distinct monophyly in such mitogenomic genealogies, i.e., the sharing of a recent common ancestor among con-specific samples collected in the same region has been viewed as evidence for subspecies. Several recent studies in cetaceans have employed this criterion to suggest subsequent intraspecific taxonomic revisions. We reason that employing intra-specific, spatially distinct monophyly at non-recombining, clonally inherited genomes is an unsatisfactory criterion for defining subspecies based upon theoretical (genetic drift) and practical (sampling effort) arguments. This point was illustrated by a re-analysis of a global mitogenomic assessment of fin whales, Balaenoptera physalus spp., published by Archer et al. (2013), which proposed to further subdivide the Northern Hemisphere fin whale subspecies, B. p. physalus. The proposed revision was based upon the detection of spatially distinct monophyly among North Atlantic and North Pacific fin whales in a genealogy based upon complete mitochondrial genome DNA sequences. The extended analysis conducted in this study (1676 mitochondrial control region, 162 complete mitochondrial genome DNA sequences and 20 microsatellite loci genotyped in 380 samples) revealed that the apparent monophyly among North Atlantic fin whales reported by Archer et al. (2013) to be due to low sample sizes. In conclusion, defining sub-species from monophyly (i.e., the absence of para- or polyphyly) can lead to erroneous conclusions due to relatively 'trivial' aspects, such as sampling. Basic population genetic processes (i.e., genetic drift and migration) also affect the time to the most recent common ancestor and hence the probability that individuals in a sample are monophyletic

Investigating Population Genetic Structure in a Highly Mobile Marine Organism: The Minke Whale Balaenoptera acutorostrata acutorostrata in the North East Atlantic

Inferring the number of genetically distinct populations and their levels of connectivity is of key importance for the sustainable management and conservation of wildlife. This represents an extra challenge in the marine environment where there are few physical barriers to gene-flow, and populations may overlap in time and space. Several studies have investigated the population genetic structure within the North Atlantic minke whale with contrasting results. In order to address this issue, we analyzed ten microsatellite loci and 331 bp of the mitochondrial D-loop on 2990 whales sampled in the North East Atlantic in the period 2004 and 2007–2011. The primary findings were: (1) No spatial or temporal genetic differentiations were observed for either class of genetic marker. (2) mtDNA identified three distinct mitochondrial lineages without any underlying geographical pattern. (3) Nuclear markers showed evidence of a single panmictic population in the NE Atlantic according STRUCTURE's highest average likelihood found at K = 1. (4) When K = 2 was accepted, based on the Evanno's test, whales were divided into two more or less equally sized groups that showed significant genetic differentiation between them but without any sign of underlying geographic pattern. However, mtDNA for these individuals did not corroborate the differentiation. (5) In order to further evaluate the potential for cryptic structuring, a set of 100 in silico generated panmictic populations was examined using the same procedures as above showing genetic differentiation between two artificially divided groups, similar to the aforementioned observations. This demonstrates that clustering methods may spuriously reveal cryptic genetic structure. Based upon these data, we find no evidence to support the existence of spatial or cryptic population genetic structure of minke whales within the NE Atlantic. However, in order to conclusively evaluate population structure within this highly mobile species, more markers will be required

Accurate Inference of Subtle Population Structure (and Other Genetic Discontinuities) Using Principal Coordinates

Accurate inference of genetic discontinuities between populations is an essential component of intraspecific biodiversity and evolution studies, as well as associative genetics. The most widely-used methods to infer population structure are model-based, Bayesian MCMC procedures that minimize Hardy-Weinberg and linkage disequilibrium within subpopulations. These methods are useful, but suffer from large computational requirements and a dependence on modeling assumptions that may not be met in real data sets. Here we describe the development of a new approach, PCO-MC, which couples principal coordinate analysis to a clustering procedure for the inference of population structure from multilocus genotype data.PCO-MC uses data from all principal coordinate axes simultaneously to calculate a multidimensional "density landscape", from which the number of subpopulations, and the membership within subpopulations, is determined using a valley-seeking algorithm. Using extensive simulations, we show that this approach outperforms a Bayesian MCMC procedure when many loci (e.g. 100) are sampled, but that the Bayesian procedure is marginally superior with few loci (e.g. 10). When presented with sufficient data, PCO-MC accurately delineated subpopulations with population F(st) values as low as 0.03 (G'(st)>0.2), whereas the limit of resolution of the Bayesian approach was F(st) = 0.05 (G'(st)>0.35).We draw a distinction between population structure inference for describing biodiversity as opposed to Type I error control in associative genetics. We suggest that discrete assignments, like those produced by PCO-MC, are appropriate for circumscribing units of biodiversity whereas expression of population structure as a continuous variable is more useful for case-control correction in structured association studies