PHFinder: assisted detection of point heteroplasmy in Sanger sequencing chromatograms

Heteroplasmy is the presence of two or more organellar genomes (mitochondrial or plastid DNA) in an organism, tissue, cell or organelle. Heteroplasmy can be detected by visual inspection of Sanger sequencing chromatograms, where it appears as multiple peaks of fluorescence at a single nucleotide position. Visual inspection of chromatograms is both consuming and highly subjective, as heteroplasmy is difficult to differentiate from background noise. Few software solutions are available to automate the detection of point heteroplasmies, and those that are available are typically proprietary, lack customization or are unsuitable for automated heteroplasmy assessment in large datasets. Here, we present PHFinder, a Python-based, open-source tool to assist in the detection of point heteroplasmies in large numbers of Sanger chromatograms. PHFinder automatically identifies point heteroplasmies directly from the chromatogram trace data. The program was tested with Sanger sequencing data from 100 humpback whales (Megaptera novaeangliae) tissue samples with known heteroplasmies. PHFinder detected most (90%) of the known heteroplasmies thereby greatly reducing the amount of visual inspection required. PHFinder is flexible and enables explicit specification of key parameters to infer double peaks (i.e., heteroplasmies).

DNA methylation markers of age(ing) in non-model animals

Inferring the chronological and biological age of individuals is fundamental to population ecology and our understanding of ageing itself, its evolution, and the biological processes that affect or even cause ageing. Epigenetic clocks based on DNA methylation (DNAm) at specific CpG sites show a strong correlation with chronological age in humans, and discrepancies between inferred and actual chronological age predict morbidity and mortality. Recently, a growing number of epigenetic clocks have been developed in non-model animals and we here review these studies. We also conduct a meta-analysis to assess the effects of different aspects of experimental protocol on the performance of epigenetic clocks for non-model animals. Two measures of performance are usually reported, the R2 of the association between the predicted and chronological age, and the mean/median absolute deviation (MAD) of estimated age from chronological age, and we argue that only the MAD reflects accuracy. R2 for epigenetic clocks based on the HorvathMammalMethylChip4 was higher and the MAD scaled to age range lower, compared with other DNAm quantification approaches. Scaled MAD tended to be lower among individuals in captive populations, and decreased with an increasing number of CpG sites. We conclude that epigenetic clocks can predict chronological age with relatively high accuracy, suggesting great potential in ecological epigenetics. We discuss general aspects of epigenetic clocks in the hope of stimulating further DNAm-based research on ageing, and perhaps more importantly, other key traits.

Wild pedigrees inform mutation rates and historic abundance in baleen whales

Phylogeny-based estimates suggesting a low germline mutation rate (μ) in baleen whales have influenced research ranging from assessments of whaling impacts to evolutionary cancer biology. We estimated μ directly from pedigrees in four baleen whale species for both the mitochondrial control region and nuclear genome. The results suggest values higher than those obtained through phylogeny-based estimates and similar to pedigree-based values for primates and toothed whales. Applying our estimate of μ reduces previous genetic-based estimates of preexploitation whale abundance by 86% and suggests that μ cannot explain low cancer rates in gigantic mammals. Our study shows that it is feasible to estimate μ directly from pedigrees in natural populations, with wide-ranging implications for ecological and evolutionary research.

Threatened North African seagrass meadows have supported green turtle populations for millennia

"Protect and restore ecosystems and biodiversity" is the second official aim of the current UN Ocean Decade (2021 to 2030) calling for the identification and protection of critical marine habitats. However, data to inform policy are often lacking altogether or confined to recent times, preventing the establishment of long-term baselines. The unique insights gained from combining bioarchaeology (palaeoproteomics, stable isotope analysis) with contemporary data (from satellite tracking) identified habitats which sea turtles have been using in the Eastern Mediterranean over five millennia. Specifically, our analysis of archaeological green turtle (Chelonia mydas) bones revealed that they likely foraged on the same North African seagrass meadows as their modern-day counterparts. Here, millennia-long foraging habitat fidelity has been directly demonstrated, highlighting the significance (and long-term dividends) of protecting these critical coastal habitats that are especially vulnerable to global warming. We highlight the potential for historical ecology to inform policy in safeguarding critical marine habitats.

Long-term stability in the circumpolar foraging range of a Southern Ocean predator between the eras of whaling and rapid climate change

Assessing environmental changes in Southern Ocean ecosystems is difficult due to its remoteness and data sparsity. Monitoring marine predators that respond rapidly to environmental variation may enable us to track anthropogenic effects on ecosystems. Yet, many long-term datasets of marine predators are incomplete because they are spatially constrained and/or track ecosystems already modified by industrial fishing and whaling in the latter half of the 20th century. Here, we assess the contemporary offshore distribution of a wide-ranging marine predator, the southern right whale (SRW, Eubalaena australis), that forages on copepods and krill from ~30°S to the Antarctic ice edge (>60°S). We analyzed carbon and nitrogen isotope values of 1,002 skin samples from six genetically distinct SRW populations using a customized assignment approach that accounts for temporal and spatial variation in the Southern Ocean phytoplankton isoscape. Over the past three decades, SRWs increased their use of mid-latitude foraging grounds in the south Atlantic and southwest (SW) Indian oceans in the late austral summer and autumn and slightly increased their use of high-latitude (>60°S) foraging grounds in the SW Pacific, coincident with observed changes in prey distribution and abundance on a circumpolar scale. Comparing foraging assignments with whaling records since the 18th century showed remarkable stability in use of mid-latitude foraging areas. We attribute this consistency across four centuries to the physical stability of ocean fronts and resulting productivity in mid-latitude ecosystems of the Southern Ocean compared with polar regions that may be more influenced by recent climate change.

How genomics can help biodiversity conservation

The availability of public genomic resources can greatly assist biodiversity assessment, conservation, and restoration efforts by providing evidence for scientifically informed management decisions. Here we survey the main approaches and applications in biodiversity and conservation genomics, considering practical factors, such as cost, time, prerequisite skills, and current shortcomings of applications. Most approaches perform best in combination with reference genomes from the target species or closely related species. We review case studies to illustrate how reference genomes can facilitate biodiversity research and conservation across the tree of life. We conclude that the time is ripe to view reference genomes as fundamental resources and to integrate their use as a best practice in conservation genomics.

Seagrass ecosystem multifunctionality under the rise of a flagship marine megaherbivore

Large grazers (megaherbivores) have a profound impact on ecosystem functioning. However, how ecosystem multifunctionality is affected by changes in megaherbivore populations remains poorly understood. Understanding the total impact on ecosystem multifunctionality requires an integrative ecosystem approach, which is especially challenging to obtain in marine systems. We assessed the effects of experimentally simulated grazing intensity scenarios on ecosystem functions and multifunctionality in a tropical Caribbean seagrass ecosystem. As a model, we selected a key marine megaherbivore, the green turtle, whose ecological role is rapidly unfolding in numerous foraging areas where populations are recovering through conservation after centuries of decline, with an increase in recorded overgrazing episodes. To quantify the effects, we employed a novel integrated index of seagrass ecosystem multifunctionality based upon multiple, well-recognized measures of seagrass ecosystem functions that reflect ecosystem services. Experiments revealed that intermediate turtle grazing resulted in the highest rates of nutrient cycling and carbon storage, while sediment stabilization, decomposition rates, epifauna richness, and fish biomass are highest in the absence of turtle grazing. In contrast, intense grazing resulted in disproportionally large effects on ecosystem functions and a collapse of multifunctionality. These results imply that (i) the return of a megaherbivore can exert strong effects on coastal ecosystem functions and multifunctionality, (ii) conservation efforts that are skewed toward megaherbivores, but ignore their key drivers like predators or habitat, will likely result in overgrazing-induced loss of multifunctionality, and (iii) the multifunctionality index shows great potential as a quantitative tool to assess ecosystem performance. Considerable and rapid alterations in megaherbivore abundance (both through extinction and conservation) cause an imbalance in ecosystem functioning and substantially alter or even compromise ecosystem services that help to negate global change effects. An integrative ecosystem approach in environmental management is urgently required to protect and enhance ecosystem multifunctionality.

Strong and lasting impacts of past global warming on baleen whales and their prey

Global warming is affecting the population dynamics and trophic interactions across a wide range of ecosystems and habitats. Translating these real-time effects into their long-term consequences remains a challenge. The rapid and extreme warming period that occurred after the Last Glacial Maximum (LGM) during the Pleistocene-Holocene transition (7-12 thousand years ago) provides an opportunity to gain insights into the long-term responses of natural populations to periods with global warming. The effects of this post-LGM warming period have been assessed in many terrestrial taxa, whereas insights into the impacts of rapid global warming on marine taxa remain limited, especially for megafauna. In order to understand how large-scale climate fluctuations during the post-LGM affected baleen whales and their prey, we conducted an extensive, large-scale analysis of the long-term effects of the post-LGM warming on abundance and inter-ocean connectivity in eight baleen whale and seven prey (fish and invertebrates) species across the Southern and the North Atlantic Ocean; two ocean basins that differ in key oceanographic features. The analysis was based upon 7032 mitochondrial DNA sequences as well as genome-wide DNA sequence variation in 100 individuals. The estimated temporal changes in genetic diversity during the last 30,000 years indicated that most baleen whale populations underwent post-LGM expansions in both ocean basins. The increase in baleen whale abundance during the Holocene was associated with simultaneous changes in their prey and climate. Highly correlated, synchronized and exponential increases in abundance in both baleen whales and their prey in the Southern Ocean were indicative of a dramatic increase in ocean productivity. In contrast, the demographic fluctuations observed in baleen whales and their prey in the North Atlantic Ocean were subtle, varying across taxa and time. Perhaps most important was the observation that the ocean-wide expansions and decreases in abundance that were initiated by the post-LGM global warming, continued for millennia after global temperatures stabilized, reflecting persistent, long-lasting impacts of global warming on marine fauna.

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.

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.

Decadal shift in foraging strategy of a migratory southern ocean predator

Rapid anthropogenic environmental change is expected to impact a host of ecological parameters in Southern Ocean ecosystems. Of critical concern are the consequences of these changes on the range of species that show fidelity to migratory destinations, as philopatry is hypothesized to help or hinder adaptation to climate change depending on the circumstances. Many baleen whales show philopatry to feeding grounds and are also capital breeders that meet migratory and reproductive costs through seasonal energy intake. Southern right whales (Eubalaena australis, SRWs) are capital breeders that have a strong relationship between reproductive output and foraging success. The population dynamics of South Africa's population of SRWs are characterized by two distinct periods: the 1990s, a period of high calving rates; and the late 2010s, a period associated with lowered calving rates. Here we use analyses of stable carbon (δ13 C) and nitrogen (δ15 N) isotope values from SRW biopsy samples (n = 122) collected during these two distinct periods to investigate foraging ecology of the South African population of SRWs over a time period coincident with the demographic shift. We show that South African SRWs underwent a dramatic northward shift, and diversification, in foraging strategy from 1990s to 2010s. Bayesian mixing model results suggest that during the 1990s, South African SRWs foraged on prey isotopically similar to South Georgia/Islas Georgias del Sur krill. In contrast, in the 2010s, South African SRWs foraged on prey isotopically consistent with the waters of the Subtropical Convergence, Polar Front and Marion Island. We hypothesize that this shift represents a response to changes in preferred habitat or prey, for example, the decrease in abundance and southward range contraction of Antarctic krill. By linking reproductive decline to changing foraging strategies for the first time in SRWs, we show that altering foraging strategies may not be sufficient to adapt to a changing ocean.

Population recovery changes population composition at a major southern Caribbean juvenile developmental habitat for the green turtle, Chelonia mydas

Understanding the population composition and dynamics of migratory megafauna at key developmental habitats is critical for conservation and management. The present study investigated whether differential recovery of Caribbean green turtle (Chelonia mydas) rookeries influenced population composition at a major juvenile feeding ground in the southern Caribbean (Lac Bay, Bonaire, Caribbean Netherlands) using genetic and demographic analyses. Genetic divergence indicated a strong temporal shift in population composition between 2006-2007 and 2015-2016 (ϕ_ST = 0.101, P < 0.001). Juvenile recruitment (<75.0 cm straight carapace length; SCL) from the north-western Caribbean increased from 12% to 38% while recruitment from the eastern Caribbean region decreased from 46% to 20% between 2006-2007 and 2015-2016. Furthermore, the product of the population growth rate and adult female abundance was a significant predictor for population composition in 2015-2016. Our results may reflect early warning signals of declining reproductive output at eastern Caribbean rookeries, potential displacement effects of smaller rookeries by larger rookeries, and advocate for genetic monitoring as a useful method for monitoring trends in juvenile megafauna. Furthermore, these findings underline the need for adequate conservation of juvenile developmental habitats and a deeper understanding of the interactions between megafaunal population dynamics in different habitats.

Long-term isolation at a low effective population size greatly reduced genetic diversity in Gulf of California fin whales

The Gulf of California, Mexico is home to many cetacean species, including a presumed resident population of fin whales, Balaenoptera physalus. Past studies reported very low levels of genetic diversity among Gulf of California fin whales and a significant level of genetic differentiation from con-specifics in the eastern North Pacific. The aim of the present study was to assess the degree and timing of the isolation of Gulf of California fin whales in a population genetic analysis of 18 nuclear microsatellite genotypes from 402 samples and 565 mitochondrial control region DNA sequences (including mitochondrial sequences retrieved from NCBI). The analyses revealed that the Gulf of California fin whale population was founded ~2.3 thousand years ago and has since remained at a low effective population size (~360) and isolated from the eastern North Pacific (N_em between 0.89-1.4). The low effective population size and high degree of isolation implied that Gulf of California fin whales are vulnerable to the negative effects of genetic drift, human-caused mortality and habitat change.

Return to the Sea, Get Huge, Beat Cancer: An Analysis of Cetacean Genomes Including an Assembly for the Humpback Whale (Megaptera novaeangliae)

Cetaceans are a clade of highly specialized aquatic mammals that include the largest animals that have ever lived. The largest whales can have ∼1,000× more cells than a human, with long lifespans, leaving them theoretically susceptible to cancer. However, large-bodied and long-lived animals do not suffer higher risks of cancer mortality than humans-an observation known as Peto's Paradox. To investigate the genomic bases of gigantism and other cetacean adaptations, we generated a de novo genome assembly for the humpback whale (Megaptera novaeangliae) and incorporated the genomes of ten cetacean species in a comparative analysis. We found further evidence that rorquals (family Balaenopteridae) radiated during the Miocene or earlier, and inferred that perturbations in abundance and/or the interocean connectivity of North Atlantic humpback whale populations likely occurred throughout the Pleistocene. Our comparative genomic results suggest that the evolution of cetacean gigantism was accompanied by strong selection on pathways that are directly linked to cancer. Large segmental duplications in whale genomes contained genes controlling the apoptotic pathway, and genes inferred to be under accelerated evolution and positive selection in cetaceans were enriched for biological processes such as cell cycle checkpoint, cell signaling, and proliferation. We also inferred positive selection on genes controlling the mammalian appendicular and cranial skeletal elements in the cetacean lineage, which are relevant to extensive anatomical changes during cetacean evolution. Genomic analyses shed light on the molecular mechanisms underlying cetacean traits, including gigantism, and will contribute to the development of future targets for human cancer therapies.

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

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.

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.

Shift of grey seal subspecies boundaries in response to climate, culling and conservation

Identifying the processes that drive changes in the abundance and distribution of natural populations is a central theme in ecology and evolution. Many species of marine mammals have experienced dramatic changes in abundance and distribution due to climatic fluctuations and anthropogenic impacts. However, thanks to conservation efforts, some of these species have shown remarkable population recovery and are now recolonizing their former ranges. Here, we use zooarchaeological, demographic and genetic data to examine processes of colonization, local extinction and recolonization of the two northern European grey seal subspecies inhabiting the Baltic Sea and North Sea. The zooarchaeological and genetic data suggest that the two subspecies diverged shortly after the formation of the Baltic Sea approximately 4200 years bp, probably through a gradual shift to different breeding habitats and phenologies. By comparing genetic data from 19th century pre-extinction material with that from seals currently recolonizing their past range, we observed a marked spatiotemporal shift in subspecies boundaries, with increasing encroachment of North Sea seals on areas previously occupied by the Baltic Sea subspecies. Further, both demographic and genetic data indicate that the two subspecies have begun to overlap geographically and are hybridizing in a narrow contact zone. Our findings provide new insights into the processes of colonization, extinction and recolonization and have important implications for the management of grey seals across northern Europe.

Adapting to a warmer ocean--seasonal shift of baleen whale movements over three decades

Global warming poses particular challenges to migratory species, which face changes to the multiple environments occupied during migration. For many species, the timing of migration between summer and winter grounds and also within-season movements are crucial to maximise exploitation of temporarily abundant prey resources in feeding areas, themselves adapting to the warming planet. We investigated the temporal variation in the occurrence of fin (Balaenoptera physalus) and humpback whales (Megaptera novaeangliae) in a North Atlantic summer feeding ground, the Gulf of St. Lawrence (Canada), from 1984 to 2010 using a long-term study of individually identifiable animals. These two sympatric species both shifted their date of arrival at a previously undocumented rate of more than 1 day per year earlier over the study period thus maintaining the approximate 2-week difference in arrival of the two species and enabling the maintenance of temporal niche separation. However, the departure date of both species also shifted earlier but at different rates resulting in increasing temporal overlap over the study period indicating that this separation may be starting to erode. Our analysis revealed that the trend in arrival was strongly related to earlier ice break-up and rising sea surface temperature, likely triggering earlier primary production. The observed changes in phenology in response to ocean warming are a remarkable example of phenotypic plasticity and may partly explain how baleen whales were able to survive a number of changes in climate over the last several million years. However, it is questionable whether the observed rate of change in timing can be maintained. Substantial modification to the distribution or annual life cycle of these species might be required to keep up with the ongoing warming of the oceans.

Utility of telomere length measurements for age determination of humpback whales

This study examines the applicability of telomere length measurements by quantitative PCR as a tool for minimally invasive age determination of free-ranging cetaceans. We analysed telomere length in skin samples from 28 North Atlantic humpback whales (Megaptera novaeangliae), ranging from 0 to 26 years of age. The results suggested a significant correlation between telomere length and age in humpback whales. However, telomere length was highly variable among individuals of similar age, suggesting that telomere length measured by quantitative PCR is an imprecise determinant of age in humpback whales. The observed variation in individual telomere length was found to be a function of both experimental and biological variability, with the latter perhaps reflecting patterns of inheritance, resource allocation trade-offs, and stochasticity of the marine environment.

Fin whale MDH-1 and MPI allozyme variation is not reflected in the corresponding DNA sequences

The appeal of genetic inference methods to assess population genetic structure and guide management efforts is grounded in the correlation between the genetic similarity and gene flow among populations. Effects of such gene flow are typically genomewide; however, some loci may appear as outliers, displaying above or below average genetic divergence relative to the genomewide level. Above average population, genetic divergence may be due to divergent selection as a result of local adaptation. Consequently, substantial efforts have been directed toward such outlying loci in order to identify traits subject to local adaptation. Here, we report the results of an investigation into the molecular basis of the substantial degree of genetic divergence previously reported at allozyme loci among North Atlantic fin whale (Balaenoptera physalus) populations. We sequenced the exons encoding for the two most divergent allozyme loci (MDH-1 and MPI) and failed to detect any nonsynonymous substitutions. Following extensive error checking and analysis of additional bioinformatic and morphological data, we hypothesize that the observed allozyme polymorphisms may reflect phenotypic plasticity at the cellular level, perhaps as a response to nutritional stress. While such plasticity is intriguing in itself, and of fundamental evolutionary interest, our key finding is that the observed allozyme variation does not appear to be a result of genetic drift, migration, or selection on the MDH-1 and MPI exons themselves, stressing the importance of interpreting allozyme data with caution. As for North Atlantic fin whale population structure, our findings support the low levels of differentiation found in previous analyses of DNA nucleotide loci.

Inferring recent historic abundance from current genetic diversity

Recent historic abundance is an elusive parameter of great importance for conserving endangered species and understanding the pre-anthropogenic state of the biosphere. The number of studies that have used population genetic theory to estimate recent historic abundance from contemporary levels of genetic diversity has grown rapidly over the last two decades. Such assessments often yield unexpectedly large estimates of historic abundance. We review the underlying theory and common practices of estimating recent historic abundance from contemporary genetic diversity, and critically evaluate the potential issues at various estimation steps. A general issue of mismatched spatio-temporal scales between the estimation itself and the objective of the estimation emerged from our assessment; genetic diversity-based estimates of recent historic abundance represent long-term averages, whereas the objective typically is an estimate of recent abundance for a specific population. Currently, the most promising approach to estimate the difference between recent historic and contemporary abundance requires that genetic data be collected from samples of similar spatial and temporal duration. Novel genome-enabled inference methods may be able to utilize additional information of dense genome-wide distributions of markers, such as of identity-by-descent tracts, to infer recent historic abundance from contemporary samples only.

Empirical evaluation of humpback whale telomere length estimates; quality control and factors causing variability in the singleplex and multiplex qPCR methods

BACKGROUND

Telomeres, the protective cap of chromosomes, have emerged as powerful markers of biological age and life history in model and non-model species. The qPCR method for telomere length estimation is one of the most common methods for telomere length estimation, but has received recent critique for being too error-prone and yielding unreliable results. This critique coincides with an increasing awareness of the potentials and limitations of the qPCR technique in general and the proposal of a general set of guidelines (MIQE) for standardization of experimental, analytical, and reporting steps of qPCR. In order to evaluate the utility of the qPCR method for telomere length estimation in non-model species, we carried out four different qPCR assays directed at humpback whale telomeres, and subsequently performed a rigorous quality control to evaluate the performance of each assay.

RESULTS

Performance differed substantially among assays and only one assay was found useful for telomere length estimation in humpback whales. The most notable factors causing these inter-assay differences were primer design and choice of using singleplex or multiplex assays. Inferred amplification efficiencies differed by up to 40% depending on assay and quantification method, however this variation only affected telomere length estimates in the worst performing assays.

CONCLUSION

Our results suggest that seemingly well performing qPCR assays may contain biases that will only be detected by extensive quality control. Moreover, we show that the qPCR method for telomere length estimation can be highly precise and accurate, and thus suitable for telomere measurement in non-model species, if effort is devoted to optimization at all experimental and analytical steps. We conclude by highlighting a set of quality controls which may serve for further standardization of the qPCR method for telomere length estimation, and discuss some of the factors that may cause variation in qPCR experiments.

Reliability of genetic bottleneck tests for detecting recent population declines

The identification of population bottlenecks is critical in conservation because populations that have experienced significant reductions in abundance are subject to a variety of genetic and demographic processes that can hasten extinction. Genetic bottleneck tests constitute an appealing and popular approach for determining if a population decline has occurred because they only require sampling at a single point in time, yet reflect demographic history over multiple generations. However, a review of the published literature indicates that, as typically applied, microsatellite-based bottleneck tests often do not detect bottlenecks in vertebrate populations known to have experienced declines. This observation was supported by simulations that revealed that bottleneck tests can have limited statistical power to detect bottlenecks largely as a result of limited sample sizes typically used in published studies. Moreover, commonly assumed values for mutation model parameters do not appear to encompass variation in microsatellite evolution observed in vertebrates and, on average, the proportion of multi-step mutations is underestimated by a factor of approximately two. As a result, bottleneck tests can have a higher probability of 'detecting' bottlenecks in stable populations than expected based on the nominal significance level. We provide recommendations that could add rigor to inferences drawn from future bottleneck tests and highlight new directions for the characterization of demographic history.

How many genetic markers to tag an individual? An empirical assessment of false matching rates among close relatives

Genetic identification of individuals is now commonplace, enabling the application of tagging methods to elusive species or species that cannot be tagged by traditional methods. A key aspect is determining the number of loci required to ensure that different individuals have non-matching multi-locus genotypes. Closely related individuals are of particular concern because of elevated matching probabilities caused by their recent co-ancestry. This issue may be addressed by increasing the number of loci to a level where full siblings (the relatedness category with the highest matching probability) are expected to have non-matching multi-locus genotypes. However, increasing the number of loci to meet this "full-sib criterion" greatly increases the laboratory effort, which in turn may increase the genotyping error rate resulting in an upward-biased mark-recapture estimate of abundance as recaptures are missed due to genotyping errors. We assessed the contribution of false matches from close relatives among 425 maternally related humpback whales, each genotyped at 20 microsatellite loci. We observed a very low (0.5-4%) contribution to falsely matching samples from pairs of first-order relatives (i.e., parent and offspring or full siblings). The main contribution to falsely matching individuals from close relatives originated from second-order relatives (e.g., half siblings), which was estimated at 9%. In our study, the total number of observed matches agreed well with expectations based upon the matching probability estimated for unrelated individuals, suggesting that the full-sib criterion is overly conservative, and would have required a 280% relative increase in effort. We suggest that, under most circumstances, the overall contribution to falsely matching samples from close relatives is likely to be low, and hence applying the full-sib criterion is unnecessary. In those cases where close relatives may present a significant issue, such as unrepresentative sampling, we propose three different genotyping strategies requiring only a modest increase in effort, which will greatly reduce the number of false matches due to the presence of related individuals.

Detecting populations in the 'ambiguous' zone: kinship-based estimation of population structure at low genetic divergence

Identifying population structure is one of the most common and important objectives of spatial analyses using population genetic data. Population structure is detected either by rejecting the null hypothesis of a homogenous distribution of genetic variation, or by estimating low migration rates. Issues arise with most current population genetic inference methods when the genetic divergence is low among putative populations. Low levels of genetic divergence may be as a result of either high ongoing migration or historic high migration but no current, ongoing migration. We direct attention to recent developments in the use of the tempo-spatial distribution of closely related individuals to detect population structure or estimate current migration rates. These 'kinship-based' approaches complement more traditional population-based genetic inference methods by providing a means to detect population structure and estimate current migration rates when genetic divergence is low. However, for kinship-based methods to become widely adopted, formal estimation procedures applicable to a range of species life histories are needed.

Genetic analyses of historic and modern marbled murrelets suggest decoupling of migration and gene flow after habitat fragmentation

The dispersal of individuals among fragmented populations is generally thought to prevent genetic and demographic isolation, and ultimately reduce extinction risk. In this study, we show that a century of reduction in coastal old-growth forests, as well as a number of other environmental factors, has probably resulted in the genetic divergence of marbled murrelets (Brachyramphus marmoratus) in central California, despite the fact that 7 per cent of modern-sampled murrelets in this population were classified as migrants using genetic assignment tests. Genetic differentiation appears to persist because individuals dispersing from northern populations contributed relatively few young to the central California population, as indicated by the fact that migrants were much less likely to be members of parent-offspring pairs than residents (10.5% versus 45.4%). Moreover, a recent 1.4 per cent annual increase in the proportion of migrants in central California, without appreciable reproduction, may have masked an underlying decline in the resident population without resulting in demographic rescue. Our results emphasize the need to understand the behaviour of migrants and the extent to which they contribute offspring in order to determine whether dispersal results in gene flow and prevents declines in resident populations.

Characterizing dispersal patterns in a threatened seabird with limited genetic structure

Genetic assignment methods provide an appealing approach for characterizing dispersal patterns on ecological time scales, but require sufficient genetic differentiation to accurately identify migrants and a large enough sample size of migrants to, for example, compare dispersal between sexes or age classes. We demonstrate that assignment methods can be rigorously used to characterize dispersal patterns in a marbled murrelet (Brachyramphus marmoratus) population from central California that numbers approximately 600 individuals and is only moderately differentiated (F(ST) approximately 0.03) from larger populations to the north. We used coalescent simulations to select a significance level that resulted in a low and approximately equal expected number of type I and II errors and then used this significance level to identify a population of origin for 589 individuals genotyped at 13 microsatellite loci. The proportion of migrants in central California was greatest during winter when 83% of individuals were classified as migrants compared to lower proportions during the breeding (6%) and post-breeding (8%) seasons. Dispersal was also biased toward young and female individuals, as is typical in birds. Migrants were rarely members of parent-offspring pairs, suggesting that they contributed few young to the central California population. A greater number of migrants than expected under equilibrium conditions, a lack of individuals with mixed ancestry, and a small number of potential source populations (two), likely allowed us to use assignment methods to rigorously characterize dispersal patterns for a population that was larger and less differentiated than typically thought required for the identification of migrants.

Characterizing source-sink dynamics with genetic parentage assignments

Source-sink dynamics have been suggested to characterize the population structure of many species, but the prevalence of source-sink systems in nature is uncertain because of inherent challenges in estimating migration rates among populations. Migration rates are often difficult to estimate directly with demographic methods, and indirect genetic methods are subject to a variety of assumptions that are difficult to meet or to apply to evolutionary timescales. Furthermore, such methods cannot be rigorously applied to high-gene-flow species. Here, we employ genetic parentage assignments in conjunction with demographic simulations to infer the level of immigration into a putative sink population. We use individual-based demographic models to estimate expected distributions of parent-offspring dyads under competing sink and closed-population models. By comparing the actual number of parent-offspring dyads (identified from multilocus genetic profiles) in a random sample of individuals taken from a population to expectations under these two contrasting demographic models, it is possible to estimate the rate of immigration and test hypotheses related to the role of immigration on population processes on an ecological timescale. The difference in the expected number of parent-offspring dyads between the two population models was greatest when immigration into the sink population was high, indicating that unlike traditional population genetic inference models, the highest degree of statistical power is achieved for the approach presented here when migration rates are high. We used the proposed genetic parentage approach to demonstrate that a threatened population of Marbled Murrelets (Braclhyrarmphus marmotus) appears to be supplemented by a low level of immigration (approximately 2-6% annually) from other populations.

Radiation and speciation of pelagic organisms during periods of global warming: the case of the common minke whale, Balaenoptera acutorostrata

How do populations of highly mobile species inhabiting open environments become reproductively isolated and evolve into new species? We test the hypothesis that elevated ocean-surface temperatures can facilitate allopatry among pelagic populations and thus promote speciation. Oceanographic modelling has shown that increasing surface temperatures cause localization and reduction of upwelling, leading to fragmentation of feeding areas critical to pelagic species. We test our hypothesis by genetic analyses of populations of two closely related baleen whales, the Antarctic minke whale (Balaenoptera bonaerensis) and common minke whale (Balaenoptera acutorostrata) whose current distributions and migration patterns extent are largely determined by areas of consistent upwelling with high primary production. Phylogeographic and population genetic analyses of mitochondrial DNA control-region nucleotide sequences collected from 467 whales sampled in four different ocean basins were employed to infer the evolutionary relationship among populations of B. acutorostrata by rooting an intraspecific phylogeny with a population of B. bonaerensis. Our findings suggest that the two species diverged in the Southern Hemisphere less than 5 million years ago (Ma). This estimate places the speciation event during a period of extended global warming in the Pliocene. We propose that elevated ocean temperatures in the period facilitated allopatric speciation by disrupting the continuous belt of upwelling maintained by the Antarctic Circumpolar Current. Our analyses revealed that the current populations of B. acutorostrata likely diverged after the Pliocene some 1.5 Ma when global temperatures had decreased and presumably coinciding with the re-establishment of the polar-equatorial temperature gradient that ultimately drives upwelling. In most population samples, we detected genetic signatures of exponential population expansions, consistent with the notion of increasing carrying capacity after the Pliocene. Our hypothesis that prolonged periods of global warming facilitate speciation in pelagic marine species that depend on upwelling should be tested by comparative analyses in other pelagic species.

Elevation and connectivity define genetic refugia for mountain sheep as climate warms

Global warming is predicted to affect the evolutionary potential of natural populations. We assessed genetic diversity of 25 populations of desert bighorn sheep (Ovis canadensis nelsoni) in southeastern California, where temperatures have increased and precipitation has decreased during the 20th century. Populations in low-elevation habitats had lower genetic diversity, presumably reflecting more fluctuations in population sizes and founder effects. Higher-elevation habitats acted as reservoirs of genetic diversity. However, genetic diversity was also affected by population connectivity, which has been disrupted by human development. Restoring population connectivity may be necessary to buffer the effects of climate change on this desert-adapted ungulate.

Identification of management units using population genetic data

The identification of management units (MUs) is central to the management of natural populations and is crucial for monitoring the effects of human activity upon species abundance. Here, we propose that the identification of MUs from population genetic data should be based upon the amount of genetic divergence at which populations become demographically independent instead of the current criterion that focuses on rejecting panmixia. MU status should only be assigned when the observed estimate of genetic divergence is significantly greater than a predefined threshold value. We emphasize the need for a demographic interpretation of estimates of genetic divergence given that it is often the dispersal rate of individuals that is the parameter of immediate interest to conservationists rather than the historical amount of gene flow.

DNA registers of legally obtained wildlife and derived products as means to identify illegal takes

The exploitation and sale of wildlife species that are endangered in only part of their range present regulators with the critical challenge of separating legal from illegal takes. Wildlife DNA registers created from tissue samples of legally obtained individual wildlife specimens can address this problem by allowing managers to identify unregistered (presumably illegally obtained) specimens. We tested the effectiveness of the only current, fully operational wildlife DNA register of individual genetic profiles collected from legally caught minke whales (Balaenoptera acutorostrata). Twenty minke whale tissue samples collected at markets in Norway and 2 additional samples collected from beached minke whales in Denmark were genotyped at 12 loci used by the Norwegian minke whale DNA register Genetic profiles of these samples then were compared against the 2676 individual profiles deposited in the Norwegian register The high number of genetic markers used to identify individuals in our study allowed consistent matching of sample and reference profiles despite an overall error rate (due to experimental and interlaboratory data standardization) estimated at 0.015 per locus. Of the 22 test samples only the 2 Danish samples failed to match an existing profile in the Norwegian minke whale DNA register Our results show that the basic principle of wildlife DNA registers can work in a real-life situation. The strength of wildlife DNA registers lies in their ability to unambiguously identify unregistered specimens with the aid of sensitive genetic methods that enable analysis of highly processed or degraded tissue samples. Our study also highlights a number of methodological problems such as laboratory errors and interlaboratory data standardization, which need be addressed to ensure a successful implementation of wildlife DNA registers.

Discerning between recurrent gene flow and recent divergence under a finite-site mutation model applied to North Atlantic and Mediterranean Sea fin whale (Balaenoptera physalus) populations

Genetic divergence among conspecific subpopulations can be due to either low recurrent gene flow or recent divergence and no gene flow. Here we present a modification of an earlier method developed by Nielsen and Wakeley (2001), which accommodates a finite-site mutation model, to assess which of the two models of divergence is most likely given the observed data. We apply the method to nucleotide sequence data collected from the variable part of the mitochondrial control region in fin whales (Balaenoptera physalus) from the Atlantic coast off Spain and the Mediterranean Sea. Our estimations strongly favor a model of recurrent gene flow over a model of recent divergence and zero gene flow. We estimated the migration rate at two females per generation. While the estimated rate is high by evolutionary standards, exchange rates of this order of magnitude is low from an ecological and conservation perspective and entirely consistent with the current paucity of fin whale sightings in the Strait of Gibraltar today. Intensive commercial shore-based whaling during the 1920s removed substantial numbers of fin whales in the Strait of Gibraltar and this local population has seemingly since failed to recover.