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

Neutral markers reveal complex population structure across the range of a widespread songbird

Understanding how both contemporary and historical physical barriers influence gene flow is key to reconstructing evolutionary histories and can allow us to predict species' resilience to changing environmental conditions. During the last glacial maximum (LGM), many high latitude North American bird species were forced into glacial refugia, including mountain bluebirds (Silia currucoides). Within their current breeding range, mountain bluebirds still experience a wide variety of environmental conditions and barriers that may disrupt gene flow and isolate populations. Using single nucleotide polymorphisms (SNPs) obtained through restriction site-associated DNA sequencing, we detected at least four genetically distinct mountain bluebird populations. Based on this structure, we determined that isolation-by-distance, the northern Rocky Mountains, and discontinuous habitat are responsible for the low connectivity and the overall history of each population going back to the last glacial maximum. Finally, we identified five candidate genes under balancing selection and three loci under diversifying selection. This study provides the first look at connectivity and gene flow across the range of these high-altitude and high latitude songbirds.

Using conservation genetics to prioritise management options for an endangered songbird

Genetic data can be highly informative for answering questions relevant to practical conservation efforts, but remain one of the most neglected aspects of species recovery plans. Framing genetic questions with reference to practical and tractable conservation objectives can help bypass this limitation of the application of genetics in conservation. Using a single-nucleotide polymorphism dataset from reduced-representation sequencing (DArTSeq), we conducted a genetic assessment of remnant populations of the endangered forty-spotted pardalote (Pardalotus quadragintus), a songbird endemic to Tasmania, Australia. Our objectives were to inform strategies for the conservation of genetic diversity in the species and estimate effective population sizes and patterns of inter-population movement to identify management units relevant to population conservation and habitat restoration. We show population genetic structure and identify two small populations on mainland Tasmania as 'satellites' of larger Bruny Island populations connected by migration. Our data identify management units for conservation objectives relating to genetic diversity and habitat restoration. Although our results do not indicate the immediate need to genetically manage populations, the small effective population sizes we estimated for some populations indicate that they are vulnerable to genetic drift, highlighting the urgent need to implement habitat restoration to increase population size and to conduct genetic monitoring. We discuss how our genetic assessment can be used to inform management interventions for the forty-spotted pardalote and show that by assessing contemporary genetic aspects, valuable information for conservation planning and decision-making can be produced to guide actions that account for genetic diversity and increase chances of recovery in species of conservation concern.

Museum Genomics Provide Evidence for Persistent Genetic Differentiation in a Threatened Seabird Species in the Western Atlantic

Connectivity among wildlife populations facilitates exchange of genetic material between groups. Changes to historical connectivity patterns resulting from anthropogenic activities can therefore have negative consequences for genetic diversity, particularly for small or isolated populations. DNA obtained from museum specimens can enable direct comparison of temporal changes in connectivity among populations, which can aid in conservation planning and contribute to the understanding of population declines. However, museum DNA can be degraded and only available in low quantities, rendering it challenging for use in population genomic analyses. Applications of genomic methodologies such as targeted sequencing address this issue by enabling capture of shared variable sites, increasing quantity and quality of recovered genomic information. We used targeted sequencing of ultra-conserved Elements (UCEs) to evaluate potential changes in connectivity and genetic diversity of roseate terns (Sterna dougallii) with a breeding distribution in the northwestern Atlantic and the Caribbean. Both populations experienced range contractions and population declines due to anthropogenic activity in the 20th century, which has the potential to alter historical connectivity regimes. Instead, we found that the two populations were differentiated historically as well as contemporaneously, with little evidence of migration between them for either time period. We also found no evidence for temporal changes in genetic diversity, although these interpretations may have been limited due to sequencing artifacts caused by the degraded nature of the museum samples. Population structuring in migratory seabirds is typically reflective of low rates of divergence and high connectivity among geographically segregated subpopulations. Our contrasting results suggest the potential presence of ecological mechanisms driving population differentiation, and highlight the value of targeted sequencing on DNA derived from museum specimens to uncover long-term patterns of genetic differentiation in wildlife populations.

Gene Pool Boundaries for the Yosemite Toad (Anaxyrus canorus) Reveal Asymmetrical Migration Within Meadow Neighborhoods

The Yosemite toad (Anaxyrus [Bufo] canorus) is a federally threatened species of meadow-specializing amphibian endemic to the high-elevation Sierra Nevada Mountains of California. The species is one of the first amphibians to undergo a large demographic collapse that was well-documented, and is reputed to remain in low abundance throughout its range. Recent phylogeographic work has demonstrated that Pleistocene toad lineages diverged and then admixed to differing extents across an elevational gradient. Although lineage divisions may have significant effects on evolutionary trajectories over large spatial and temporal scales, present-day population dynamics must be delineated in order to manage and conserve the species effectively. In this study, we used a double-digest RADseq dataset to address three primary questions: (1) Are single meadows or neighborhoods of nearby meadows most correlated with population boundaries? (2) Does asymmetrical migration occur among neighborhoods of nearby meadows? (3) What topographic or hydrological variables predict such asymmetrical migration in these meadow neighborhoods? Hierarchical STRUCTURE and AMOVA analyses suggested that populations are typically circumscribed by a single meadow, although 84% of meadows exist in neighborhoods of at least two meadows connected by low levels of migration, and over half (53%) of neighborhoods examined display strong asymmetrical migration. Meadow neighborhoods often contain one or more large and flat “hub” meadows that experience net immigration, surrounded by smaller and topographically rugged “satellite” meadows with net emigration. Hubs tend to contain more genetic diversity and could be prioritized for conservation and habitat management and as potential sources for reestablishment efforts.

Effects of sample age on data quality from targeted sequencing of museum specimens: what are we capturing in time?

Background

Next generation sequencing (NGS) can recover DNA data from valuable extant and extinct museum specimens. However, archived or preserved DNA is difficult to sequence because of its fragmented, damaged nature, such that the most successful NGS methods for preserved specimens remain sub-optimal. Improving wet-lab protocols and comprehensively determining the effects of sample age on NGS library quality are therefore of vital importance. Here, I examine the relationship between sample age and several indicators of library quality following targeted NGS sequencing of ~ 1300 loci using 271 samples of pinned moth specimens (Helicoverpa armigera) ranging in age from 5 to 117 years.

Results

I find that older samples have lower DNA concentrations following extraction and thus require a higher number of indexing PCR cycles during library preparation. When sequenced reads are aligned to a reference genome or to only the targeted region, older samples have a lower number of sequenced and mapped reads, lower mean coverage, and lower estimated library sizes, while the percentage of adapters in sequenced reads increases significantly as samples become older. Older samples also show the poorest capture success, with lower enrichment and a higher improved coverage anticipated from further sequencing.

Conclusions

Sample age has significant, measurable impacts on the quality of NGS data following targeted enrichment. However, incorporating a uracil-removing enzyme into the blunt end-repair step during library preparation could help to repair DNA damage, and using a method that prevents adapter-dimer formation may result in improved data yields.

Projecting Suitability and Climate Vulnerability of Bhutanitis thaidina (Blanchard) (Lepidoptera: Papilionidae) with Conservation Implications

Bhutanitis thaidina is an endemic, rare, and protected swallowtail in China. Deforestation, habitat fragmentation, illegal commercialised capture, and exploitation of larval food plants are believed to be the four major causes of population decline of B. thaidina in the recent decade. However, little attention was paid to the impact of climate change. This study used ecological niche factor analysis and species distribution model to analyse the current suitable areas for B. thaidina with BioClim variables as well as its future suitable areas under four future climate scenarios (represented by four Representative Concentration Pathways: RCP2.6, RCP4.5, RCP6.0, and RCP8.5). Statistical analysis was carried out to compare the possible area and altitude changes to the distribution of B. thaidina under changing climate. Our analyses showed that the suitable areas for B. thaidina are fragmented under the current climate, with four suitable centres in northwestern Yunnan, northeastern Yunnan and northwestern Guizhou, the western margin of Sichuan Basin, and Qinling mountains. Apart from further habitat fragmentation under climate change, slight range expansion (average 6.0–8.9%) was detected under the RCP2.6 and RCP4.5 scenarios, while more range contraction (average 1.3–26.9%) was detected under the RCP6.0 and RCP8.5 scenarios, with the two southern suitable centres suffering most. Also, a tendency of contraction (2,500–3,500 m) and upslope shift (~600 m) in suitable altitude range were detected. The findings of this study supported the climate-vulnerable hypothesis of B. thaidina, especially under future climate like the RCP6.0 and RCP8.5 scenarios, in terms of contraction in suitable areas and altitude ranges. Conservation priority should be given to northwestern Yunnan, northeastern Yunnan, and northwestern Guizhou to alleviate the stress of massive habitat loss and extinction. Refugial areas should be established in all four suitable centres to maintain genetic diversity of B. thaidina in China.

Historical DNA as a tool to address key questions in avian biology and evolution: A review of methods, challenges, applications, and future directions

Museum specimens play a crucial role in addressing key questions in systematics, evolution, ecology, and conservation. With the advent of high-throughput sequencing technologies, specimens that have long been the foundation of important biological discoveries can inform new perspectives as sources of genomic data. Despite the many possibilities associated with analyzing DNA from historical specimens, several challenges persist. Using avian systems as a model, we review DNA extraction protocols, sequencing technologies, and capture methods that are helping researchers overcome some of these difficulties. We highlight empirical examples in which researchers have used these technologies to address fundamental questions related to avian conservation and evolution. Increasing accessibility to new sequencing technologies will provide researchers with tools to tap into the wealth of information contained within our valuable natural history collections.

Density-dependent fitness, not dispersal movements, drives temporal variation in spatial genetic structure in dark-eyed juncos (Junco hyemalis)

Some studies have found that dispersal rates and distances increase with density, indicating that density-dependent dispersal likely affects spatial genetic structure. In an 11-year mark-recapture study on a passerine, the dark-eyed junco, we tested whether density affected dispersal distance and/or fine-scale spatial genetic structure. Contrary to expectations, we found no effect of predispersal density on dispersal distance or the proportion of locally produced juveniles returning to the population from which they hatched. However, even though density did not affect dispersal distance or natal return rates, we found that density still did affect spatial genetic structure. We found significant positive spatial genetic structure at low densities of (postdispersal) adults but not at high densities. In years with high postdispersal (adult) densities that also had high predispersal (juvenile) densities in the previous year, we found negative spatial genetic structure, indicating high levels of dispersal. We found that density also affected fitness of recruits, and fitness of immigrants, potentially linking these population parameters with the spatial genetic structure detected. Immigrants and recruits rarely nested in low postdispersal density years. In contrast, in years with high postdispersal density, recruits were common and immigrants had equal success to local birds, so novel genotypes diluted the gene pool and effectively eliminated positive spatial genetic structure. In relation to fine-scale spatial genetic structure, fitness of immigrants and new recruits is poorly understood compared to dispersal movements, but we conclude that it can have implications for the spatial distribution of genotypes in populations.

Museum specimens of terrestrial vertebrates are sensitive indicators of environmental change in the Anthropocene

Natural history museums and the specimen collections they curate are vital scientific infrastructure, a fact as true today as it was when biologists began collecting and preserving specimens over 200 years ago. The importance of museum specimens in studies of taxonomy, systematics, ecology and evolutionary biology is evidenced by a rich and abundant literature, yet creative and novel uses of specimens are constantly broadening the impact of natural history collections on biodiversity science and global sustainability. Excellent examples of the critical importance of specimens come from their use in documenting the consequences of environmental change, which is particularly relevant considering the alarming rate at which we now modify our planet in the Anthropocene. In this review, we highlight the important role of bird, mammal and amphibian specimens in documenting the Anthropocene and provide examples that underscore the need for continued collection of museum specimens.This article is part of the theme issue 'Biological collections for understanding biodiversity in the Anthropocene'.

Habitat Fragmentation Reduces Genetic Diversity and Connectivity of the Mexican Spotted Owl: A Simulation Study Using Empirical Resistance Models

We evaluated how differences between two empirical resistance models for the same geographic area affected predictions of gene flow processes and genetic diversity for the Mexican spotted owl (Strix occidentalis lucida). The two resistance models represented the landscape under low- and high-fragmentation parameters. Under low fragmentation, the landscape had larger but highly concentrated habitat patches, whereas under high fragmentation, the landscape had smaller habitat patches that scattered across a broader area. Overall habitat amount differed little between resistance models. We tested eight scenarios reflecting a factorial design of three factors: resistance model (low vs. high fragmentation), isolation hypothesis (isolation-by-distance, IBD, vs. isolation-by-resistance, IBR), and dispersal limit of species (200 km vs. 300 km). Higher dispersal limit generally had a positive but small influence on genetic diversity. Genetic distance increased with both geographic distance and landscape resistance, but landscape resistance displayed a stronger influence. Connectivity was positively related to genetic diversity under IBR but was less important under IBD. Fragmentation had a strong negative influence on the spatial patterns of genetic diversity and effective population size (Ns). Despite habitats being more concentrated and less widely distributed, the low-fragmentation landscape had greater genetic diversity than the high-fragmentation landscape, suggesting that highly concentrated but larger habitat patches may provide a genetic refuge for the Mexican spotted owl.

Isolating the roles of movement and reproduction on effective connectivity alters conservation priorities for an endangered bird

Movement is important for ecological and evolutionary theory as well as connectivity conservation, which is increasingly critical for species responding to environmental change. Key ecological and evolutionary outcomes of movement, such as population growth and gene flow, require effective dispersal: movement that is followed by successful reproduction. However, the relative roles of movement and postmovement reproduction for effective dispersal and connectivity remain unclear. Here we isolate the contributions of movement and immigrant reproduction to effective dispersal and connectivity across the entire breeding range of an endangered raptor, the snail kite (Rostrhamus sociabilis plumbeus). To do so, we unite mark-resight data on movement and reproduction across 9 years and 27 breeding patches with an integrated model that decomposes effective dispersal into its hierarchical levels of movement, postmovement breeding attempt, and postmovement reproductive success. We found that immigrant reproduction limits effective dispersal more than movement for this endangered species, demonstrating that even highly mobile species may have limited effective connectivity due to reduced immigrant reproduction. We found different environmental limitations for the reproductive component of effective dispersal compared with movement, indicating that different conservation strategies may be needed when promoting effective dispersal rather than movement alone. We also demonstrate that considering immigrant reproduction, rather than movement alone, alters which patches are the most essential for connectivity, thereby changing conservation priorities. These results challenge the assumption that understanding movement alone is sufficient to infer connectivity and highlight that connectivity conservation may require not only fostering movement but also successful reproduction of immigrants.

Eco-evolutionary effects on population recovery following catastrophic disturbance

Fine-scale genetic diversity and contemporary evolution can theoretically influence ecological dynamics in the wild. Such eco-evolutionary effects might be particularly relevant to the persistence of populations facing acute or chronic environmental change. However, experimental data on wild populations is currently lacking to support this notion. One way that ongoing evolution might influence the dynamics of threatened populations is through the role that selection plays in mediating the 'rescue effect', the ability of migrants to contribute to the recovery of populations facing local disturbance and decline. Here, we combine experiments with natural catastrophic events to show that ongoing evolution is a major determinant of migrant contributions to population recovery in Trinidadian guppies (Poecilia reticulata). These eco-evolutionary limits on migrant contributions appear to be mediated by the reinforcing effects of natural and sexual selection against migrants, despite the close geographic proximity of migrant sources. These findings show that ongoing adaptive evolution can be a double-edged sword for population persistence, maintaining local fitness at a cost to demographic risk. Our study further serves as a potent reminder that significant evolutionary and eco-evolutionary dynamics might be at play even where the phenotypic status quo is largely maintained generation to generation.

Unlocking the vault: next-generation museum population genomics

Natural history museum collections provide unique resources for understanding how species respond to environmental change, including the abrupt, anthropogenic climate change of the past century. Ideally, researchers would conduct genome-scale screening of museum specimens to explore the evolutionary consequences of environmental changes, but to date such analyses have been severely limited by the numerous challenges of working with the highly degraded DNA typical of historic samples. Here, we circumvent these challenges by using custom, multiplexed, exon capture to enrich and sequence ~11,000 exons (~4 Mb) from early 20th-century museum skins. We used this approach to test for changes in genomic diversity accompanying a climate-related range retraction in the alpine chipmunks (Tamias alpinus) in the high Sierra Nevada area of California, USA. We developed robust bioinformatic pipelines that rigorously detect and filter out base misincorporations in DNA derived from skins, most of which likely resulted from postmortem damage. Furthermore, to accommodate genotyping uncertainties associated with low-medium coverage data, we applied a recently developed probabilistic method to call single-nucleotide polymorphisms and estimate allele frequencies and the joint site frequency spectrum. Our results show increased genetic subdivision following range retraction, but no change in overall genetic diversity at either nonsynonymous or synonymous sites. This case study showcases the advantages of integrating emerging genomic and statistical tools in museum collection-based population genomic applications. Such technical advances greatly enhance the value of museum collections, even where a pre-existing reference is lacking and points to a broad range of potential applications in evolutionary and conservation biology.

Gene flow and demographic history of leopards (Panthera pardus) in the central Indian highlands

Gene flow is a critical ecological process that must be maintained in order to counteract the detrimental effects of genetic drift in subdivided populations, with conservation benefits ranging from promoting the persistence of small populations to spreading adaptive traits in changing environments. We evaluated historical and contemporary gene flow and effective population sizes of leopards in a landscape in central India using noninvasive sampling. Despite the dramatic changes in land-use patterns in this landscape through recent times, we did not detect any signs that the leopard populations have been through a genetic bottleneck, and they appear to have maintained migration-drift equilibrium. We found that historical levels of gene flow (mean m h = 0.07) were significantly higher than contemporary levels (mean m c = 0.03), and populations with large effective population sizes (Satpura and Kanha Tiger Reserves) are the larger exporters of migrants at both timescales. The greatest decline in historical versus contemporary gene flow is between pairs of reserves that are currently not connected by forest corridors (i.e., Melghat-Pench m h - m c = 0.063; and Kanha-Satpura m h - m c = 0.054). We attribute this reduction in gene flow to accelerated fragmentation and habitat alteration in the landscape over the past few centuries, and suggest protection of forest corridors to maintain gene flow in this landscape.

Museum genomics: low-cost and high-accuracy genetic data from historical specimens

Natural history collections are unparalleled repositories of geographical and temporal variation in faunal conditions. Molecular studies offer an opportunity to uncover much of this variation; however, genetic studies of historical museum specimens typically rely on extracting highly degraded and chemically modified DNA samples from skins, skulls or other dried samples. Despite this limitation, obtaining short fragments of DNA sequences using traditional PCR amplification of DNA has been the primary method for genetic study of historical specimens. Few laboratories have succeeded in obtaining genome-scale sequences from historical specimens and then only with considerable effort and cost. Here, we describe a low-cost approach using high-throughput next-generation sequencing to obtain reliable genome-scale sequence data from a traditionally preserved mammal skin and skull using a simple extraction protocol. We show that single-nucleotide polymorphisms (SNPs) from the genome sequences obtained independently from the skin and from the skull are highly repeatable compared to a reference genome.

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.

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.