Genomic signatures of population bottleneck and recovery in Northwest Atlantic pinnipeds

Genetic diversity and structure of a recent fish invasion: Tench (Tinca tinca) in eastern North America

Introduced and geographically expanding populations experience similar eco-evolutionary challenges, including founder events, genetic bottlenecks, and novel environments. Theory predicts that reduced genetic diversity resulting from such phenomena limits the success of introduced populations. Using 1900 SNPs obtained from restriction-site-associated DNA sequencing, we evaluated hypotheses related to the invasion history and connectivity of an invasive population of Tench (Tinca tinca), a Eurasian freshwater fish that has been expanding geographically in eastern North America for three decades. Consistent with the reported history of a single introduction event, our findings suggest that multiple introductions from distinct genetic sources are unlikely as Tench had a small effective population size (~114 [95% CI = 106-123] individuals), no strong population subdivision across time and space, and evidence of a recent genetic bottleneck. The large genetic neighbourhood size (220 km) and weak within-population genetic substructure suggested high connectivity across the invaded range, despite the relatively large area occupied. There was some evidence for a small decay in genetic diversity as the species expanded northward, but not southward, into new habitats. As eradicating the species within a ~112 km radius would be necessary to prevent recolonization, eradicating Tench is likely not feasible at watershed-and possibly local-scales. Management should instead focus on reducing abundance in priority conservation areas to mitigate adverse impacts. Our study indicates that introduced populations can thrive and exhibit relatively high levels of genetic diversity despite severe bottlenecks (<1.5% of the ancestral effective population size) and suggests that landscape heterogeneity and population demographics can generate variability in spatial patterns of genetic diversity within a single range expansion.

Sequence diversity and differences at the highly duplicated MHC-I gene reflect viral susceptibility in sympatric pinniped species

Differences in disease susceptibility among species can result from rapid host-pathogen coevolution and differences in host species ecology that affect the strength and direction of natural selection. Among two sympatric pinniped species that differ in sociality and putative disease exposure, we investigate observed differences in susceptibility through an analysis of a highly variable, duplicated gene family involved in the vertebrate immune response. Using high-throughput amplicon sequencing, we characterize diversity at the two exons that encode the peptide binding region of the major histocompatibility complex class I (MHC-I) gene in harbor (N = 60) and gray (N = 90) seal populations from the Northwest Atlantic. Across species, we identified 106 full-length exon 2 and 103 exon 3 sequence variants and a minimum of 11 duplicated MHC-I loci. The sequence variants clustered in 15 supertypes defined by the physiochemical properties of the peptide binding region, including a putatively novel Northwest Atlantic MHC-I diversity sublineage. Trans-species polymorphisms, dN/dS ratios, and evidence of gene conversion among supertypes are consistent with balancing selection acting on this gene. High functional redundancy suggests particularly strong selection among gray seals at the novel Northwest Atlantic MHC-I diversity sublineage. At exon 2, harbor seals had a significantly greater number of variants per individual than gray seals, but fewer supertypes. Supertype richness and private supertypes are hypothesized to contribute to observed differences in disease resistance between species, as consistently, across the North Atlantic and many disease outbreaks, gray seals appear to be more resistant to respiratory viruses than harbor seals.

The reconstruction of invasion histories with genomic data in light of differing levels of anthropogenic transport

Unravelling the history of range shifts is key for understanding past, current and future species distributions. Anthropogenic transport of species alters natural dispersal patterns and directly affects population connectivity. Studies have suggested that high levels of anthropogenic transport homogenize patterns of genetic differentiation and blur colonization pathways. However, empirical evidence of these effects remains elusive. We compared two range-shifting species (Microcosmus squamiger and Ciona robusta) to examine how anthropogenic transport affects our ability to reconstruct colonization pathways using genomic data. We first investigated shipping networks from the 18th century onwards, cross-referencing these with regions where the species have records to infer how each species has potentially been affected by different levels of anthropogenic transport. We then genotyped thousands of single-nucleotide polymorphisms from 280 M. squamiger and 190 C. robusta individuals collected across their extensive species' ranges and reconstructed colonization pathways. Differing levels of anthropogenic transport did not preclude the elucidation of population structure, though specific inferences of colonization pathways were difficult to discern in some of the considered scenario sets. We conclude that genomic data in combination with information of underlying introduction drivers provide key insights into the historic spread of range-shifting species.

This article is part of the theme issue ‘Species’ ranges in the face of changing environments (part I)’.

Origin and expansion of the world's most widespread pinniped: range-wide population genomics of the harbour seal (Phoca vitulina)

The harbour seal (Phoca vitulina) is the most widely distributed pinniped, occupying a wide variety of habitats and climatic zones across the Northern Hemisphere. Intriguingly, the harbour seal is also one of the most philopatric seals, raising questions as to how it colonised virtually the whole of the Northern Hemisphere. To shed light on the origin, remarkable range expansion, population structure and genetic diversity of this species, we used genotyping-by-sequencing to analyse ~13,500 biallelic SNPs from 286 individuals sampled from 22 localities across the species' range. Our results point to a Northeast Pacific origin, colonisation of the North Atlantic via the Canadian Arctic, and subsequent stepping-stone range expansions across the North Atlantic from North America to Europe, accompanied by a successive loss of genetic diversity. Our analyses further revealed a deep divergence between modern North Pacific and North Atlantic harbour seals, with finer-scale genetic structure at regional and local scales consistent with strong philopatry. The study provides new insights into the harbour seal's remarkable ability to colonise and adapt to a wide range of habitats. Furthermore, it has implications for current harbour seal subspecies delineations and highlights the need for international and national red lists and management plans to ensure the protection of genetically and demographically isolated populations.

Informing conservation strategies with museum genomics: Long-term effects of past anthropogenic persecution on the elusive European wildcat

Like many carnivore species, European wildcats (Felis silvestris) have suffered severe anthropogenic population declines in the past, resulting in a strong population bottleneck at the beginning of the 20th century. In Germany, the species has managed to survive its near extinction in small isolated areas and is currently recolonizing former habitats owing to legal protection and concerted conservation efforts. Here, we SNP-genotyped and mtDNA-sequenced 56 historical and 650 contemporary samples to assess the impact of massive persecution on genetic diversity, population structure, and hybridization dynamics of wildcats. Spatiotemporal analyses suggest that the presumed postglacial differentiation between two genetically distinct metapopulations in Germany is in fact the result of the anthropogenic bottleneck followed by re-expansion from few secluded refugia. We found that, despite the bottleneck, populations experienced no severe genetic erosion, nor suffered from elevated inbreeding or showed signs of increased hybridization with domestic cats. Our findings have significant implications for current wildcat conservation strategies, as the data analyses show that the two presently recognized wildcat population clusters should be treated as a single conservation unit. Although current populations appear under no imminent threat from genetic factors, fostering connectivity through the implementation of forest corridors will facilitate the preservation of genetic diversity and promote long-term viability. The present study documents how museum collections can be used as essential resource for assessing long-term anthropogenic effects on natural populations, for example, regarding population structure and the delineation of appropriate conservation units, potentially informing todays' species conservation.

Emergence and radiation of distemper viruses in terrestrial and marine mammals

Canine distemper virus (CDV) and phocine distemper virus (PDV) are major pathogens to terrestrial and marine mammals. Yet little is known about the timing and geographical origin of distemper viruses and to what extent it was influenced by environmental change and human activities. To address this, we (i) performed the first comprehensive time-calibrated phylogenetic analysis of the two distemper viruses, (ii) mapped distemper antibody and virus detection data from marine mammals collected between 1972 and 2018, and (iii) compiled historical reports on distemper dating back to the eighteenth century. We find that CDV and PDV diverged in the early seventeenth century. Modern CDV strains last shared a common ancestor in the nineteenth century with a marked radiation during the 1930s-1950s. Modern PDV strains are of more recent origin, diverging in the 1970s-1980s. Based on the compiled information on distemper distribution, the diverse host range of CDV and basal phylogenetic placement of terrestrial morbilliviruses, we hypothesize a terrestrial CDV-like ancestor giving rise to PDV in the North Atlantic. Moreover, given the estimated timing of distemper origin and radiation, we hypothesize a prominent role of environmental change such as the Little Ice Age, and human activities like globalization and war in distemper virus evolution.

Genome-wide SNP analysis of three moose subspecies at the southern range limit in the contiguous United States

Genome-wide evaluations of genetic diversity and population structure are important for informing management and conservation of trailing-edge populations. North American moose (Alces alces) are declining along portions of the southern edge of their range due to disease, species interactions, and marginal habitat, all of which may be exacerbated by climate change. We employed a genotyping by sequencing (GBS) approach in an effort to collect baseline information on the genetic variation of moose inhabiting the species’ southern range periphery in the contiguous United States. We identified 1920 single nucleotide polymorphisms (SNPs) from 155 moose representing three subspecies from five states: A. a. americana (New Hampshire), A. a. andersoni (Minnesota), and A. a. shirasi (Idaho, Montana, and Wyoming). Molecular analyses supported three geographically isolated clusters, congruent with currently recognized subspecies. Additionally, while moderately low genetic diversity was observed, there was little evidence of inbreeding. Results also indicated > 20% shared ancestry proportions between A. a. shirasi samples from northern Montana and A. a. andersoni samples from Minnesota, indicating a putative hybrid zone warranting further investigation. GBS has proven to be a simple and effective method for genome-wide SNP discovery in moose and provides robust data for informing herd management and conservation priorities. With increasing disease, predation, and climate related pressure on range edge moose populations in the United States, the use of SNP data to identify gene flow between subspecies may prove a powerful tool for moose management and recovery, particularly if hybrid moose are more able to adapt.

Genetic association with boldness and maternal performance in a free-ranging population of grey seals (Halichoerus grypus)

Individual variation in quantitative traits clearly influence many ecological and evolutionary processes. Moderate to high heritability estimates of personality and life-history traits suggest some level of genetic control over these traits. Yet, we know very little of the underlying genetic architecture of phenotypic variation in the wild. In this study, we used a candidate gene approach to investigate the association of genetic variants with repeated measures of boldness and maternal performance traits (weaning mass and lactation duration) collected over an 11- and 28-year period, respectively, in a free-ranging population of grey seals on Sable Island National Park Reserve, Canada. We isolated and re-sequenced five genes: dopamine receptor D4 (DRD4), serotonin transporter (SERT), oxytocin receptor (OXTR), and melanocortin receptors 1 (MC1R) and 5 (MC5R). We discovered single nucleotide polymorphisms (SNPs) in each gene; and, after accounting for loci in linkage disequilibrium and filtering due to missing data, we were able to test for genotype-phenotype relationships at seven loci in three genes (DRD4, SERT, and MC1R). We tested for association between these loci and traits of 180 females having extreme shy-bold phenotypes using mixed-effects models. One locus within SERT was significantly associated with boldness (effect size = 0.189) and a second locus within DRD4 with weaning mass (effect size = 0.232). Altogether, genotypes explained 6.52-13.66% of total trait variation. Our study substantiates SERT and DRD4 as important determinants of behaviour, and provides unique insight into the molecular mechanisms underlying maternal performance variation in a marine predator.

Persistent panmixia despite extreme habitat loss and population decline in the threatened tricolored blackbird (Agelaius tricolor)

Habitat loss and alteration has driven many species into decline, often to the point of requiring protection and intervention to avert extinction. Genomic data provide the opportunity to inform conservation and recovery efforts with details about vital evolutionary processes with a resolution far beyond that of traditional genetic approaches. The tricolored blackbird (Agelaius tricolor) has suffered severe losses during the previous century largely due to anthropogenic impacts on their habitat. Using a dataset composed of a whole genome paired with reduced representation libraries (RAD-Seq) from samples collected across the species' range, we find evidence for panmixia using multiple methods, including PCA (no geographic clustering), admixture analyses (ADMIXTURE and TESS conclude K = 1), and comparisons of genetic differentiation (average FST = 0.029). Demographic modeling approaches recovered an ancient decline that had a strong impact on genetic diversity but did not detect any effect from the known recent decline. We also did not detect any evidence for selection, and hence adaptive variation, at any site, either geographic or genomic. These results indicate that species continues to have high vagility across its range despite population decline and habitat loss and should be managed as a single unit.

Rates of increase in gray seal (Halichoerus grypus atlantica) pupping at recolonized sites in the United States, 1988-2019

Gray seals were historically distributed along the northeastern coast of the United States, but bounties and lack of protection reduced numbers and they were rarely observed for most of the 20th century. Once protections were enacted, the population started to rebound. Here, we describe the recolonization and recovery of gray seals in the United States, focusing on the re-establishment of pupping sites. We fit individual generalized linear models to various time series (1988–2019) to estimate rates of increase in observed pup counts at four of the more data-rich sites. Annual rate of increase at individual sites ranged from −0.2% (95% CI: −2.3–1.9%) to 26.3% (95% CI: 21.6–31.4%). The increase in sites and number of pups born in the United States is driven by population growth and immigration from Canadian colonies and is part of a larger recovery of the Northwest Atlantic population. Wildlife protection, a healthy source population, habitat availability, and species traits that allow for dispersal and high productivity were all important factors in this recovery.