Bayesian inference of a historical bottleneck in a heavily exploited marine mammal

Ancient DNA indicates a century of overhunting did not reduce genetic diversity in Pacific Walruses (Odobenus rosmarus divergens)

Pacific Walruses (Odobenus rosmarus divergens [Illiger 1815]) are gregarious marine mammals considered to be sentinels of the Arctic because of their dependence on sea ice for feeding, molting, and parturition. Like many other marine mammal species, their population sizes were decimated by historical overhunting in the nineteenth and twentieth centuries. Although they have since been protected from nearly all commercial hunting pressure, they now face rapidly accelerating habitat loss as global warming reduces the extent of summer sea ice in the Arctic. To investigate how genetic variation was impacted by overhunting, we obtained mitochondrial DNA sequences from historic Pacific Walrus samples in Alaska that predate the period of overhunting, as well as from extant populations. We found that genetic variation was unchanged over this period, suggesting Pacific Walruses are resilient to genetic attrition in response to reduced population size, and that this may be related to their high vagility and lack of population structure. Although Pacific Walruses will almost certainly continue to decline in number as the planet warms and summer sea ice is further reduced, they may be less susceptible to the ratcheting effects of inbreeding that typically accompany shrinking populations.

Genetic diversity and demographic history of the leopard seal: A Southern Ocean top predator

Leopard seals (Hydrurga leptonyx) are top predators that can exert substantial top-down control of their Antarctic prey species. However, population trends and genetic diversity of leopard seals remain understudied, limiting our understanding of their ecological role. We investigated the genetic diversity, effective population size and demographic history of leopard seals to provide fundamental data that contextualizes their predatory influence on Antarctic ecosystems. Ninety leopard seals were sampled from the northern Antarctic Peninsula during the austral summers of 2008-2019 and a 405bp segment of the mitochondrial control region was sequenced for each individual. We uncovered moderate levels of nucleotide (π = 0.013) and haplotype (Hd = 0.96) diversity, and the effective population size was estimated at around 24,000 individuals (NE = 24,376; 95% CI: 16,876-33,126). Consistent with findings from other ice-breeding pinnipeds, Bayesian skyline analysis also revealed evidence for population expansion during the last glacial maximum, suggesting that historical population growth may have been boosted by an increase in the abundance of sea ice. Although leopard seals can be found in warmer, sub-Antarctic locations, the species' core habitat is centered on the Antarctic, making it inherently vulnerable to the loss of sea ice habitat due to climate change. Therefore, detailed assessments of past and present leopard seal population trends are needed to inform policies for Antarctic ecosystems.

Phylogeography of the Atlantic Blue Crab Callinectes sapidus (Brachyura: Portunidae) in the Americas versus the Mediterranean Sea: Determining Origins and Genetic Connectivity of a Large-Scale Invasion

The American blue crab Callinectes sapidus is a particularly successful invader in estuarine ecosystems worldwide. Despite increasing awareness of its potential harm, the invasion history and underlying genetic diversity of this species within the Mediterranean Sea remain unknown. This study constitutes the first large-scale approach to study phylogeographic patterns of C. sapidus in Europe, facilitated by the first comparison of all currently available COI sequence data. For this investigation, 71 individuals of C. sapidus were newly analyzed and the entire COI gene was sequenced and used for a comparative phylogeographic analyses. For the first time, two separately used adjacent regions of this gene were combined in a single dataset. This allowed emphasizing the prevalence of three geographically defined lineages within the native range: (1) eastern North America, including the Gulf of Mexico, (2) the Caribbean, and (3) Brazil. New data from the Mediterranean reveal that non-native populations of C. sapidus are characterized by a conspicuously low genetic diversity (except for Turkey, where stocking took place), and that there is surprisingly low connectivity among established populations. The occurrence of strong genetic bottlenecks suggests few founder individuals. This confirms that, even under a scenario of restricted large-scale gene flow, a very limited number of invasive individuals is sufficient for a massive impact.

Intronic primers reveal unexpectedly high major histocompatibility complex diversity in Antarctic fur seals

The major histocompatibility complex (MHC) is a group of genes comprising one of the most important components of the vertebrate immune system. Consequently, there has been much interest in characterising MHC variation and its relationship with fitness in a variety of species. Due to the exceptional polymorphism of MHC genes, careful PCR primer design is crucial for capturing all of the allelic variation present in a given species. We therefore developed intronic primers to amplify the full-length 267 bp protein-coding sequence of the MHC class II DQB exon 2 in the Antarctic fur seal. We then characterised patterns of MHC variation among mother-offspring pairs from two breeding colonies and detected 19 alleles among 771 clone sequences from 56 individuals. The distribution of alleles within and among individuals was consistent with a single-copy, classical DQB locus showing Mendelian inheritance. Amino acid similarity at the MHC was significantly associated with genome-wide relatedness, but no relationship was found between MHC heterozygosity and genome-wide heterozygosity. Finally, allelic diversity was several times higher than reported by a previous study based on partial exon sequences. This difference appears to be related to allele-specific amplification bias, implying that primer design can strongly impact the inference of MHC diversity.

Demographic Reconstruction of Antarctic Fur Seals Supports the Krill Surplus Hypothesis

Much debate surrounds the importance of top-down and bottom-up effects in the Southern Ocean, where the harvesting of over two million whales in the mid twentieth century is thought to have produced a massive surplus of Antarctic krill. This excess of krill may have allowed populations of other predators, such as seals and penguins, to increase, a top-down hypothesis known as the 'krill surplus hypothesis'. However, a lack of pre-whaling population baselines has made it challenging to investigate historical changes in the abundance of the major krill predators in relation to whaling. Therefore, we used reduced representation sequencing and a coalescent-based maximum composite likelihood approach to reconstruct the recent demographic history of the Antarctic fur seal, a pinniped that was hunted to the brink of extinction by 18th and 19th century sealers. In line with the known history of this species, we found support for a demographic model that included a substantial reduction in population size around the time period of sealing. Furthermore, maximum likelihood estimates from this model suggest that the recovered, post-sealing population at South Georgia may have been around two times larger than the pre-sealing population. Our findings lend support to the krill surplus hypothesis and illustrate the potential of genomic approaches to shed light on long-standing questions in population biology.

Phylogeography of Otaria flavescens (Carnivora: Pinnipedia): unravelling genetic connectivity at the southernmost limit of its distribution

The Pleistocene glacial period shaped the current genetic structure of numerous species. The last glacial dynamics has been proposed to have split the South American sea lion, Otaria flavescens, into two Evolutionarily Significant Units (ESUs), one on each side of the continent. However, previous studies have not provided genetic information on colonies found along 3000 km of coastline of the southernmost limit of the species distribution, where gene flow could occur. We conducted an exhaustive phylogeographical analysis of O. flavescens using a mtDNA marker, including, for the first time, data from colonies living south of latitude 45° S, in the Argentinian provinces of Santa Cruz and Tierra del Fuego. Our results indicated the presence of five Conservation Units across the distribution range of O. flavescens and suggest that the Patagonian population must have expanded about 150 000 BP. We found evidence for gene flow across the entire species range, supporting a scenario of secondary contact in Tierra del Fuego where representatives of the oldest lineages coexist. The presence of gene flow between oceans leads us to reject the assumption of complete reciprocal monophyly for mtDNA between the presumed ESUs, suggesting that the species constitutes a single Evolutionarily Significant Unit.

Local Drivers of Change in Southern Ocean Ecosystems: Human Activities and Policy Implications

Local drivers are human activities or processes that occur in specific locations, and cause physical or ecological change at the local or regional scale. Here, we consider marine and land-derived pollution, non-indigenous species, tourism and other human visits, exploitation of marine resources, recovery of marine mammals, and coastal change as a result of ice loss, in terms of their historic and current extent, and their interactions with the Southern Ocean environment. We summarise projected increases or decreases in the influence of local drivers, and projected changes to their geographic range, concluding that the influence of non-indigenous species, fishing, and the recovery of marine mammals are predicted to increase in the future across the Southern Ocean. Local drivers can be managed regionally, and we identify existing governance frameworks as part of the Antarctic Treaty System and other instruments which may be employed to mitigate or limit their impacts on Southern Ocean ecosystems.

Spatial and temporal genetic variation in an exploited reef fish: The effects of exploitation on cohort genetic structure

Many coral reef fishes are fished, often resulting in detrimental genetic effects; however, reef fishes often show unpredictable patterns of genetic variation, which potentially mask the effects of fishing. Our goals were to characterize spatial and temporal genetic variation and determine the effects of fishing on an exploited reef fish, Plectropomus leopardus, Lacepède (the common coral trout). To determine population structure, we genotyped 417 Great Barrier Reef coral trout from four populations sampled in 2 years (1996 and 2004) at nine microsatellite loci. To test for exploitation effects, we additionally genotyped 869 individuals from a single cohort (ages 3-5) across eight different reefs, including fished and control populations. Genetic structure differed substantially in the two sampled years, with only 1 year exhibiting isolation by distance. Thus, genetic drift likely plays a role in shaping population genetic structure in this species. Although we found no loss of genetic diversity associated with exploitation, our relatedness patterns show that pulse fishing likely affects population genetics. Additionally, genetic structure in the cohort samples likely reflected spatial variation in recruitment contributing to genetic structure at the population level. Overall, we show that fishing does impact coral reef fishes, highlighting the importance of repeated widespread sampling to accurately characterize the genetic structure of reef fishes, as well as the power of analysing cohorts to avoid the impacts of recruitment-related genetic swamping. The high temporal and spatial variability in genetic structure, combined with possible selection effects, will make conservation/management of reef fish species complex.

Wintertime overlaps between female Antarctic fur seals (Arctocephalus gazella) and the krill fishery at South Georgia, South Atlantic

The diet of Antarctic fur seals (Arctocephalus gazella) at South Georgia is dominated by Antarctic krill (Euphausia superba). During the breeding season, foraging trips by lactating female fur seals are constrained by their need to return to land to provision their pups. Post-breeding, seals disperse in order to feed and recover condition; estimates indicate c.70% of females remain near to South Georgia, whilst others head west towards the Patagonian Shelf or south to the ice-edge. The krill fishery at South Georgia operates only during the winter, providing the potential for fur seal: fishery interaction during these months. Here we use available winter (May to September) tracking data from Platform Terminal Transmitter (PTT) tags deployed on female fur seals at Bird Island, South Georgia. We develop habitat models describing their distribution during the winters of 1999 and 2003 with the aim of visualising and quantifying the degree of spatial overlap between female fur seals and krill harvesting in South Georgia waters. We show that spatial distribution of fur seals around South Georgia is extensive, and that the krill fishery overlaps with small, highly localised areas of available fur seal habitat. From these findings we discuss the implications for management, and future work.

Have whales returned to a historical hotspot of industrial whaling? The pattern of southern right whale Eubalaena australis recovery at South Georgia

Around 176500 whales were killed in the sub-Antarctic waters off South Georgia (South Atlantic) between 1904 and 1965. In recent decades, whales have once again become summer visitors, with the southern right whale (SRW) the most commonly reported species until 2011. Here, we assess the distribution, temporal pattern, health status and likely prey of SRWs in these waters, combining observations from a summertime vessel-based expedition to South Georgia, stable isotope data collected from SRWs and putative prey and sightings reports collated by the South Georgia Museum. The expedition used directional acoustics and visual surveys to localise whales and collected skin biopsies and photo-IDs. During 76 h of visual observation effort over 19 expedition days, SRWs were encountered 15 times (~31 individuals). Photo-IDs, combined with publicly contributed images from commercial vessels, were reconciled and quality-controlled to form a catalogue of 6 fully (i.e. both sides) identified SRWs and 26 SRWs identified by either left or right sides. No photo-ID matches were found with lower-latitude calving grounds, but 3 whales had gull lesions supporting a direct link with Península Valdés, Argentina. The isotopic position of SRWs in the South Georgia food web suggests feeding on a combination of copepod and krill species. Opportunistic reports of SRW sightings and associated group sizes remain steady over time, while humpback whales provide a strong contrast, with increased sighting rates and group sizes seen since 2013. These data suggest a plateau in SRWs and an increasing humpback whale presence in South Georgia waters following the cessation of whaling.

An 85K SNP Array Uncovers Inbreeding and Cryptic Relatedness in an Antarctic Fur Seal Breeding Colony

High density single nucleotide polymorphism (SNP) arrays allow large numbers of individuals to be rapidly and cost-effectively genotyped at large numbers of genetic markers. However, despite being widely used in studies of humans and domesticated plants and animals, SNP arrays are lacking for most wild organisms. We developed a custom 85K Affymetrix Axiom array for an intensively studied pinniped, the Antarctic fur seal (Arctocephalus gazella). SNPs were discovered from a combination of genomic and transcriptomic resources and filtered according to strict criteria. Out of a total of 85,359 SNPs tiled on the array, 75,601 (88.6%) successfully converted and were polymorphic in 270 animals from a breeding colony at Bird Island in South Georgia. Evidence was found for inbreeding, with three genomic inbreeding coefficients being strongly intercorrelated and the proportion of the genome in runs of homozygosity being non-zero in all individuals. Furthermore, analysis of genomic relatedness coefficients identified previously unknown first-degree relatives and multiple second-degree relatives among a sample of ostensibly unrelated individuals. Such “cryptic relatedness” within fur seal breeding colonies may increase the likelihood of consanguineous matings and could therefore have implications for understanding fitness variation and mate choice. Finally, we demonstrate the cross-amplification potential of the array in three related pinniped species. Overall, our SNP array will facilitate future studies of Antarctic fur seals and has the potential to serve as a more general resource for the wider pinniped research community.

Antarctic fur seal (Arctocephalus gazella) annual migration and temporal patterns of on-shore occurrence of leucistic individuals on King George Island

Non-invasive tracking the on-shore occurrence of the atypically pigmented animals and determination of land residency duration of leucistic seals would help us find out more about the rotation of the migrating population. During seven austral summer seasons (2011–2018), by counting the animals every 10 days at the Cape Lions Rump shore, King George Island, South Shetlands, in the Antarctic Specially Protected Area No. 151 and the adjacent ice-free land (31.52 km2) we registered fourteen leucistic individuals per a total of 43,919 animals. Moreover, daily monitoring of local fauna resulted in further 33 leucistic animals (together 47, in all seasons). Whilst the results of 10-day censuses of the total population were similar inter-seasonally, a tendency for increased occurrence of leucistic individuals in successive seasons was revealed. Generally, the number of animals increased significantly as season progressed. Since leucistic individuals stayed on-shore for 1–2 days usually, it can be hypothesised that the observation of migrating Antarctic fur seals every 3 days does not involve the same individuals. Also, additional every 5-day censuses taken in one season in ASPA 151 resulted in a higher seasonal number of animals, which proves that more frequent counts help us estimate population abundance more efficiently. Thus, every 5-day counts are proposed as a feasible and justified method of population monitoring.

The genetic legacy of extreme exploitation in a polar vertebrate

Understanding the effects of human exploitation on the genetic composition of wild populations is important for predicting species persistence and adaptive potential. We therefore investigated the genetic legacy of large-scale commercial harvesting by reconstructing, on a global scale, the recent demographic history of the Antarctic fur seal (Arctocephalus gazella), a species that was hunted to the brink of extinction by 18th and 19th century sealers. Molecular genetic data from over 2,000 individuals sampled from all eight major breeding locations across the species' circumpolar geographic distribution, show that at least four relict populations around Antarctica survived commercial hunting. Coalescent simulations suggest that all of these populations experienced severe bottlenecks down to effective population sizes of around 150-200. Nevertheless, comparably high levels of neutral genetic variability were retained as these declines are unlikely to have been strong enough to deplete allelic richness by more than around 15%. These findings suggest that even dramatic short-term declines need not necessarily result in major losses of diversity, and explain the apparent contradiction between the high genetic diversity of this species and its extreme exploitation history.

A method to estimate pre-exploitation population size

Antarctic fur seals (Arctocephalus gazella) were commercially exploited on the subantarctic island of South Georgia for over 100 years and nearly driven to extinction. Since the cessation of harvesting, however, their populations have rebounded, and they are now often considered a nuisance species whose impact on the terrestrial landscape should be mitigated. Any evaluation of their current population requires the context provided by their historic, pre-exploitation abundance, lest ecologists fall prey to shifting baseline syndrome in which their perspective on current abundance is compared only with an altered state resulting from past anthropogenic disturbance. Estimating pre-exploitation abundance is critical to defining species recovery and setting recovery targets, both of which are needed for the International Union for the Conservation of Nature's recent efforts to develop a green list of recovering species. To address this issue, we reconstructed the South Georgia fur seal harvest from 1786 to 1908 from ship logbooks and other historical records and interpolated missing harvest data as necessary with a generalized linear model fit to the historical record. Using an approximate Bayesian computation framework, harvest data, and a stochastic age-structured population model, we estimated the pre-exploitation abundance of Antarctic fur seals on South Georgia was 2.5 million females (95% CI 1.5-3.5 million). This estimate is similar to recent abundance estimates, and suggests current populations, and the ecological consequences of so many fur seals on the island, may be similar to conditions prior to human harvest. Although the historic archive on the fur sealing era is unavoidably patchy, the use of archival records is essential for reconstructing the past and, correspondingly, to understanding the present. Article impact statement: Defining species recovery requires an understanding of baseline population state, which can be estimated through statistical methods.

Complete mitochondrial genome of the Antarctic fur seal (Arctocephalus gazella)

The Antarctic fur seal (Arctocephalus gazella) is an abundant Antarctic otariid. Here, we present the complete mitochondrial DNA sequence of this species, which includes 13 protein-coding genes, 22 transfer RNA genes, 2 ribosomal RNA genes, and the control region for a total length of 16,156 bp. A phylogenetic analysis including all 25 publically available pinniped mitogenomes nested the Antarctic fur seal within the Otariid clade, which was clearly resolved from the Phocidae and Odobenidae.

Genetic Characteristics of the Japanese Serow Capricornis crispus in the Kii Mountain Range, Central Japan

The Japanese serow, Capricornis crispus, is an indigenous bovid species exclusively inhabiting mountain regions in the main Japanese islands, excepting Hokkaido. It had decreased in abundance to its lowest level due to overhunting and deforestation, with its distribution severely fragmented from the middle of the 20th century, many populations of C. crispus currently facing the risk of extinction. The Kii Mountain Range (KM) on Honshu is one such location that has seen a drastic population decline of C. crispus. In this study, we examined genetic characteristics of C. crispus in KM and neighboring regions of the Chubu district, using mtDNA and microsatellite markers, in order to devise strategies for its conservation. Results for mtDNA were characterized by low nucleotide diversity with five endemic and two dominant haplotypes shared by individuals in neighboring regions. A Bayesian skyline plot indicated a gradual increase after the last glacial maximum. For microsatellites, the genetic diversity of C. crispus in KM was comparable to Shizuoka and higher than Shikoku. Recent genetic bottlenecks were strongly suggested in C. crispus in KM. Bayesian clustering showed a genetic cline between KM and neighboring regions, where multivariate analysis suggested three local populations. A Mantel test indicated male-biased dispersal. These results indicate that C. crispus in KM and neighboring regions constitute multiple local populations, connected through restricted gene flow. For the conservation of C. crispus, it is important to define small-scale conservation units, among which genetic connectivity should be facilitated to prevent further loss of genetic diversity.

Genetic population structure and demography of an apex predator, the tiger shark Galeocerdo cuvier

Population genetics has been increasingly applied to study large sharks over the last decade. Whilst large shark species are often difficult to study with direct methods, improved knowledge is needed for both population management and conservation, especially for species vulnerable to anthropogenic and climatic impacts. The tiger shark, Galeocerdo cuvier, is an apex predator known to play important direct and indirect roles in tropical and subtropical marine ecosystems. While the global and Indo-West Pacific population genetic structure of this species has recently been investigated, questions remain over population structure and demographic history within the western Indian (WIO) and within the western Pacific Oceans (WPO). To address the knowledge gap in tiger shark regional population structures, the genetic diversity of 286 individuals sampled in seven localities was investigated using 27 microsatellite loci and three mitochondrial genes (CR,COI, and cytb). A weak genetic differentiation was observed between the WIO and the WPO, suggesting high genetic connectivity. This result agrees with previous studies and highlights the importance of the pelagic behavior of this species to ensure gene flow. Using approximate Bayesian computation to couple information from both nuclear and mitochondrial markers, evidence of a recent bottleneck in the Holocene (2,000-3,000 years ago) was found, which is the most probable cause for the low genetic diversity observed. A contemporary effective population size as low as 111 [43,369] was estimated during the bottleneck. Together, these results indicate low genetic diversity that may reflect a vulnerable population sensitive to regional pressures. Conservation measures are thus needed to protect a species that is classified as Near Threatened.

Genetic signatures of small effective population sizes and demographic declines in an endangered rattlesnake, Sistrurus catenatus

Endangered species that exist in small isolated populations are at elevated risk of losing adaptive variation due to genetic drift. Analyses that estimate short-term effective population sizes, characterize historical demographic processes, and project the trajectory of genetic variation into the future are useful for predicting how levels of genetic diversity may change. Here, we use data from two independent types of genetic markers (single nucleotide polymorphisms [SNPs] and microsatellites) to evaluate genetic diversity in 17 populations spanning the geographic range of the endangered eastern massasauga rattlesnake (Sistrurus catenatus). First, we use SNP data to confirm previous reports that these populations exhibit high levels of genetic structure (overall Fst = 0.25). Second, we show that most populations have contemporary Ne estimates <50. Heterozygosity-fitness correlations in these populations provided no evidence for a genetic cost to living in small populations, though these tests may lack power. Third, model-based demographic analyses of individual populations indicate that all have experienced declines, with the onset of many of these declines occurring over timescales consistent with anthropogenic impacts (<200 years). Finally, forward simulations of the expected loss of variation in relatively large (Ne = 50) and small (Ne = 10) populations indicate they will lose a substantial amount of their current standing neutral variation (63% and 99%, respectively) over the next 100 years. Our results argue that drift has a significant and increasing impact on levels of genetic variation in isolated populations of this snake, and efforts to assess and mitigate associated impacts on adaptive variation should be components of the management of this endangered reptile.

Demographic histories and genetic diversity across pinnipeds are shaped by human exploitation, ecology and life-history

A central paradigm in conservation biology is that population bottlenecks reduce genetic diversity and population viability. In an era of biodiversity loss and climate change, understanding the determinants and consequences of bottlenecks is therefore an important challenge. However, as most studies focus on single species, the multitude of potential drivers and the consequences of bottlenecks remain elusive. Here, we combined genetic data from over 11,000 individuals of 30 pinniped species with demographic, ecological and life history data to evaluate the consequences of commercial exploitation by 18th and 19th century sealers. We show that around one third of these species exhibit strong signatures of recent population declines. Bottleneck strength is associated with breeding habitat and mating system variation, and together with global abundance explains much of the variation in genetic diversity across species. Overall, bottleneck intensity is unrelated to IUCN status, although the three most heavily bottlenecked species are endangered. Our study reveals an unforeseen interplay between human exploitation, animal biology, demographic declines and genetic diversity.

Historical hunting has caused documented declines in pinnipeds, but the extent to which hunting caused genetic bottlenecks among species was unknown. Here, the authors show evidence of severe bottlenecks in several pinniped species, particularly those that breed on land.

A global cline in a colour polymorphism suggests a limited contribution of gene flow towards the recovery of a heavily exploited marine mammal

Evaluating how populations are connected by migration is important for understanding species resilience because gene flow can facilitate recovery from demographic declines. We therefore investigated the extent to which migration may have contributed to the global recovery of the Antarctic fur seal (Arctocephalus gazella), a circumpolar distributed marine mammal that was brought to the brink of extinction by the sealing industry in the eighteenth and nineteenth centuries. It is widely believed that animals emigrating from South Georgia, where a relict population escaped sealing, contributed to the re-establishment of formerly occupied breeding colonies across the geographical range of the species. To investigate this, we interrogated a genetic polymorphism (S291F) in the melanocortin 1 receptor gene, which is responsible for a cream-coloured phenotype that is relatively abundant at South Georgia and which appears to have recently spread to localities as far afield as Marion Island in the sub-Antarctic Indian Ocean. By sequencing a short region of this gene in 1492 pups from eight breeding colonies, we showed that S291F frequency rapidly declines with increasing geographical distance from South Georgia, consistent with locally restricted gene flow from South Georgia mainly to the South Shetland Islands and Bouvetøya. The S291F allele was not detected farther afield, suggesting that although emigrants from South Georgia may have been locally important, they are unlikely to have played a major role in the recovery of geographically more distant populations.

Reduced representation sequencing detects only subtle regional structure in a heavily exploited and rapidly recolonizing marine mammal species

Next-generation reduced representation sequencing (RRS) approaches show great potential for resolving the structure of wild populations. However, the population structure of species that have shown rapid demographic recovery following severe population bottlenecks may still prove difficult to resolve due to high gene flow between subpopulations. Here, we tested the effectiveness of the RRS method Genotyping-By-Sequencing (GBS) for describing the population structure of the New Zealand fur seal (NZFS, Arctocephalus forsteri), a species that was heavily exploited by the 19th century commercial sealing industry and has since rapidly recolonized most of its former range from a few isolated colonies. Using 26,026 neutral single nucleotide polymorphisms (SNPs), we assessed genetic variation within and between NZFS colonies. We identified low levels of population differentiation across the species range (<1% of variation explained by regional differences) suggesting a state of near panmixia. Nonetheless, we observed subtle population substructure between West Coast and Southern East Coast colonies and a weak, but significant (p = 0.01), isolation-by-distance pattern among the eight colonies studied. Furthermore, our demographic reconstructions supported severe bottlenecks with potential 10-fold and 250-fold declines in response to Polynesian and European hunting, respectively. Finally, we were able to assign individuals treated as unknowns to their regions of origin with high confidence (96%) using our SNP data. Our results indicate that while it may be difficult to detect population structure in species that have experienced rapid recovery, next-generation markers and methods are powerful tools for resolving fine-scale structure and informing conservation and management efforts.

RAD Sequencing and a Hybrid Antarctic Fur Seal Genome Assembly Reveal Rapidly Decaying Linkage Disequilibrium, Global Population Structure and Evidence for Inbreeding

Recent advances in high throughput sequencing have transformed the study of wild organisms by facilitating the generation of high quality genome assemblies and dense genetic marker datasets. These resources have the potential to significantly advance our understanding of diverse phenomena at the level of species, populations and individuals, ranging from patterns of synteny through rates of linkage disequilibrium (LD) decay and population structure to individual inbreeding. Consequently, we used PacBio sequencing to refine an existing Antarctic fur seal (Arctocephalus gazella) genome assembly and genotyped 83 individuals from six populations using restriction site associated DNA (RAD) sequencing. The resulting hybrid genome comprised 6,169 scaffolds with an N50 of 6.21 Mb and provided clear evidence for the conservation of large chromosomal segments between the fur seal and dog (Canis lupus familiaris). Focusing on the most extensively sampled population of South Georgia, we found that LD decayed rapidly, reaching the background level by around 400 kb, consistent with other vertebrates but at odds with the notion that fur seals experienced a strong historical bottleneck. We also found evidence for population structuring, with four main Antarctic island groups being resolved. Finally, appreciable variance in individual inbreeding could be detected, reflecting the strong polygyny and site fidelity of the species. Overall, our study contributes important resources for future genomic studies of fur seals and other pinnipeds while also providing a clear example of how high throughput sequencing can generate diverse biological insights at multiple levels of organization.

Effective population size and the genetic consequences of commercial whaling on the humpback whales (Megaptera novaeangliae) from Southwestern Atlantic Ocean

Genotypes of 10 microsatellite loci of 420 humpback whales from the Southwestern Atlantic Ocean population were used to estimate for the first time its contemporary effective (N_e) and census (N_c) population sizes and to test the genetic effect of commercial whaling. The results are in agreement with our previous studies that found high genetic diversity for this breeding population. Using an approximate Bayesian computation approach, the scenario of constant N_e was significantly supported over scenarios with moderate to strong size changes during the commercial whaling period. The previous generation N_c (N_e multiplied by 3.6), which should corresponds to the years between around 1980 and 1990, was estimated between ~2,600 and 6,800 whales (point estimate ~4,000), and is broadly compatible with the recent abundance surveys extrapolated to the past using a growth rate of 7.4% per annum. The long-term N_c in the constant scenario (point estimate ~15,000) was broadly compatible (considering the confidence interval) with pre-whaling catch records estimates (point estimate ~25,000). Overall, our results shown that the Southwestern Atlantic Ocean humpback whale population is genetically very diverse and resisted well to the strong population reduction during commercial whaling.

Ecological and genetic impact of the 2011 Tohoku Earthquake Tsunami on intertidal mud snails

Natural disturbances often destroy local populations and can considerably affect the genetic properties of these populations. The 2011 Tohoku Earthquake Tsunami greatly damaged local populations of various coastal organisms, including the mud snail Batillaria attramentaria, which was an abundant macroinvertebrate on the tidal flats in the Tohoku region. To evaluate the impact of the tsunami on the ecology and population genetic properties of these snails, we monitored the density, shell size, and microsatellite DNA variation of B. attramentaria for more than ten years (2005–2015) throughout the disturbance event. We found that the density of snails declined immediately after the tsunami. Bayesian inference of the genetically effective population size (N_e) demonstrated that the N_e declined by 60–99% at the study sites exposed to the tsunami. However, we found that their genetic diversity was not significantly reduced after the tsunami. The maintenance of genetic diversity is essential for long-term survival of local populations, and thus, the observed genetic robustness could play a key role in the persistence of snail populations in this region which has been devastated by similar tsunamis every 500–800 years. Our findings have significant implications for understanding the sustainability of populations damaged by natural disturbances.

Low Spatial Genetic Differentiation Associated with Rapid Recolonization in the New Zealand Fur Seal Arctocephalus forsteri

Population declines resulting from anthropogenic activities are of major consequence for the long-term survival of species because the resulting loss of genetic diversity can lead to extinction via the effects of inbreeding depression, fixation of deleterious mutations, and loss of adaptive potential. Otariid pinnipeds have been exploited commercially to near extinction with some species showing higher demographic resilience and recolonization potential than others. The New Zealand fur seal (NZFS) was heavily impacted by commercial sealing between the late 18th and early 19th centuries, but has recolonized its former range in southern Australia. The species has also recolonized its former range in New Zealand, yet little is known about the pattern of recolonization. Here, we first used 11 microsatellite markers (n = 383) to investigate the contemporary population structure and dispersal patterns in the NZFS (Arctocephalus forsteri). Secondly, we model postsealing recolonization with 1 additional mtDNA cytochrome b (n = 261) marker. Our data identified 3 genetic clusters: an Australian, a subantarctic, and a New Zealand one, with a weak and probably transient subdivision within the latter cluster. Demographic history scenarios supported a recolonization of the New Zealand coastline from remote west coast colonies, which is consistent with contemporary gene flow and with the species' high resilience. The present data suggest the management of distinct genetic units in the North and South of New Zealand along a genetic gradient. Assignment of individuals to their colony of origin was limited (32%) with the present data indicating the current microsatellite markers are unlikely sufficient to assign fisheries bycatch of NZFSs to colonies.

Born blonde: a recessive loss-of-function mutation in the melanocortin 1 receptor is associated with cream coat coloration in Antarctic fur seals

Although the genetic basis of color variation has been extensively studied in humans and domestic animals, the genetic polymorphisms responsible for different color morphs remain to be elucidated in many wild vertebrate species. For example, hypopigmentation has been observed in numerous marine mammal species but the underlying mutations have not been identified. A particularly compelling candidate gene for explaining color polymorphism is the melanocortin 1 receptor (MC1R), which plays a key role in the regulation of pigment production. We therefore used Antarctic fur seals (Arctocephalus gazella) as a highly tractable marine mammal system with which to test for an association between nucleotide variation at the MC1R and melanin‐based coat color phenotypes. By sequencing 70 wild‐type individuals with dark‐colored coats and 26 hypopigmented individuals with cream‐colored coats, we identified a nonsynonymous mutation that results in the substitution of serine with phenylalanine at an evolutionarily highly conserved structural domain. All of the hypopigmented individuals were homozygous for the allele coding for phenylalanine, consistent with a recessive loss‐of‐function allele. In order to test for cryptic population structure, which can generate artefactual associations, and to evaluate whether homozygosity at the MC1R could be indicative of low genome‐wide heterozygosity, we also genotyped all of the individuals at 50 polymorphic microsatellite loci. We were unable to detect any population structure and also found that wild‐type and hypopigmented individuals did not differ significantly in their standardized multilocus heterozygosity. Such a lack of association implies that hypopigmented individuals are unlikely to suffer disproportionately from inbreeding depression, and hence, we have no reason to believe that they are at a selective disadvantage in the wider population.

Population structure and historical demography of South American sea lions provide insights into the catastrophic decline of a marine mammal population

Understanding the causes of population decline is crucial for conservation management. We therefore used genetic analysis both to provide baseline data on population structure and to evaluate hypotheses for the catastrophic decline of the South American sea lion (Otaria flavescens) at the Falkland Islands (Malvinas) in the South Atlantic. We genotyped 259 animals from 23 colonies across the Falklands at 281 bp of the mitochondrial hypervariable region and 22 microsatellites. A weak signature of population structure was detected, genetic diversity was moderately high in comparison with other pinniped species, and no evidence was found for the decline being associated with a strong demographic bottleneck. By combining our mitochondrial data with published sequences from Argentina, Brazil, Chile and Peru, we also uncovered strong maternally directed population structure across the geographical range of the species. In particular, very few shared haplotypes were found between the Falklands and South America, and this was reflected in correspondingly low migration rate estimates. These findings do not support the prominent hypothesis that the decline was caused by migration to Argentina, where large-scale commercial harvesting operations claimed over half a million animals. Thus, our study not only provides baseline data for conservation management but also reveals the potential for genetic studies to shed light upon long-standing questions pertaining to the history and fate of natural populations.

High genetic diversity in the endangered and narrowly distributed amphibian species Leptobrachium leishanense

Threatened species typically have a small or declining population size, which make them highly susceptible to loss of genetic diversity through genetic drift and inbreeding. Genetic diversity determines the evolutionary potential of a species; therefore, maintaining the genetic diversity of threatened species is essential for their conservation. In this study, we assessed the genetic diversity of the adaptive major histocompatibility complex (MHC) genes in an endangered and narrowly distributed amphibian species, Leptobrachium leishanense in Southwest China. We compared the genetic variation of MHC class I genes with that observed in neutral markers (5 microsatellite loci and cytochrome b gene) to elucidate the relative roles of genetic drift and natural selection in shaping the current MHC polymorphism in this species. We found a high level of genetic diversity in this population at both MHC and neutral markers compared with other threatened amphibian species. Historical positive selection was evident in the MHC class I genes. The higher allelic richness in MHC markers compared with that of microsatellite loci suggests that selection rather than genetic drift plays a prominent role in shaping the MHC variation pattern, as drift can affect all the genome in a similar way but selection directly targets MHC genes. Although demographic analysis revealed no recent bottleneck events in L. leishanense, additional population decline will accelerate the dangerous status for this species. We suggest that the conservation management of L. leishanense should concentrate on maximizing the retention of genetic diversity through preventing their continuous population decline. Protecting their living habitats and forbidding illegal hunting are the most important measures for conservation of L. leishanense.

Ancient and contemporary DNA reveal a pre-human decline but no population bottleneck associated with recent human persecution in the kea (Nestor notabilis)

The impact of population bottlenecks is an important factor to consider when assessing species survival. Population declines can considerably limit the evolutionary potential of species and make them more susceptible to stochastic events. New Zealand has a well documented history of decline of endemic avifauna related to human colonization. Here, we investigate the genetic effects of a recent population decline in the endangered kea (Nestor notabilis). Kea have undergone a long-lasting persecution between the late 1800s to 1970s where an estimated 150,000 kea were culled under a governmental bounty scheme. Kea now number 1,000–5,000 individuals in the wild and it is likely that the recent population decline may have reduced the genetic diversity of the species. Comparison of contemporary (n = 410), historical (n = 15) and fossil samples (n = 4) showed a loss of mitochondrial diversity since the end of the last glaciation (Otiran Glacial) but no loss of overall genetic diversity associated with the cull. Microsatellite data indicated a recent bottleneck for only one population and a range-wide decline in N_e dating back some 300 – 6,000 years ago, a period predating European arrival in NZ. These results suggest that despite a recent human persecution, kea might have experienced a large population decline before stabilizing in numbers prior to human settlement of New Zealand in response to Holocene changes in habitat distribution. Our study therefore highlights the need to understand the respective effects of climate change and human activities on endangered species dynamics when proposing conservation guidelines.

Naturally rare versus newly rare: demographic inferences on two timescales inform conservation of Galápagos giant tortoises

Long-term population history can influence the genetic effects of recent bottlenecks. Therefore, for threatened or endangered species, an understanding of the past is relevant when formulating conservation strategies. Levels of variation at neutral markers have been useful for estimating local effective population sizes (N e ) and inferring whether population sizes increased or decreased over time. Furthermore, analyses of genotypic, allelic frequency, and phylogenetic information can potentially be used to separate historical from recent demographic changes. For 15 populations of Galápagos giant tortoises (Chelonoidis sp.), we used 12 microsatellite loci and DNA sequences from the mitochondrial control region and a nuclear intron, to reconstruct demographic history on shallow (past ∽100 generations, ∽2500 years) and deep (pre-Holocene, >10 thousand years ago) timescales. At the deep timescale, three populations showed strong signals of growth, but with different magnitudes and timing, indicating different underlying causes. Furthermore, estimated historical N e of populations across the archipelago showed no correlation with island age or size, underscoring the complexity of predicting demographic history a priori. At the shallow timescale, all populations carried some signature of a genetic bottleneck, and for 12 populations, point estimates of contemporary N e were very small (i.e., < 50). On the basis of the comparison of these genetic estimates with published census size data, N e generally represented ∽0.16 of the census size. However, the variance in this ratio across populations was considerable. Overall, our data suggest that idiosyncratic and geographically localized forces shaped the demographic history of tortoise populations. Furthermore, from a conservation perspective, the separation of demographic events occurring on shallow versus deep timescales permits the identification of naturally rare versus newly rare populations; this distinction should facilitate prioritization of management action.

Contrasting genetic diversity and population structure among three sympatric Madagascan shorebirds: parallels with rarity, endemism, and dispersal

Understanding the relative contributions of intrinsic and extrinsic factors to population structure and genetic diversity is a central goal of conservation and evolutionary genetics. One way to achieve this is through comparative population genetic analysis of sympatric sister taxa, which allows evaluation of intrinsic factors such as population demography and life history while controlling for phylogenetic relatedness and geography. We used ten conserved microsatellites to explore the population structure and genetic diversity of three sympatric and closely related plover species in southwestern Madagascar: Kittlitz's plover (Charadrius pecuarius), white-fronted plover (C. marginatus), and Madagascar plover (C. thoracicus). Bayesian clustering revealed strong population structure in the rare and endemic Madagascar plover, intermediate population structure in the white-fronted plover, and no detectable population structure in the geographically widespread Kittlitz's plover. In contrast, allelic richness and heterozygosity were highest for the Kittlitz's plover, intermediate for the white-fronted plover and lowest for the Madagascar plover. No evidence was found in support of the “watershed mechanism” proposed to facilitate vicariant divergence of Madagascan lemurs and reptiles, which we attribute to the vagility of birds. However, we found a significant pattern of genetic isolation by distance among populations of the Madagascar plover, but not for the other two species. These findings suggest that interspecific variation in rarity, endemism, and dispersal propensity may influence genetic structure and diversity, even in highly vagile species.

Demographic response of cutlassfish (Trichiurus japonicus and T. nanhaiensis) to fluctuating palaeo-climate and regional oceanographic conditions in the China seas

Glacial cycles of the Quaternary have heavily influenced the demographic history of various species. To test the evolutionary impact of palaeo-geologic and climatic events on the demographic history of marine taxa from the coastal Western Pacific, we investigated the population structure and demographic history of two economically important fish (Trichiurus japonicus and T. nanhaiensis) that inhabit the continental shelves of the East China and northern South China Seas using the mitochondrial cytochrome b sequences and Bayesian Skyline Plot analyses. A molecular rate of 2.03% per million years, calibrated to the earliest flooding of the East China Sea shelf (70-140 kya), revealed a strong correlation between population sizes and primary production. Furthermore, comparison of the demographic history of T. japonicus populations from the East China and South China Seas provided evidence of the postglacial development of the Changjiang (Yangtze River) Delta. In the South China Sea, interspecific comparisons between T. japonicus and T. nanhaiensis indicated possible evolutionary responses to changes in palaeo-productivity that were influenced by East Asian winter monsoons. This study not only provides insight into the demographic history of cutlassfish but also reveals potential clues regarding the historic productivity and regional oceanographic conditions of the Western Pacific marginal seas.

Demographic histories and genetic diversities of Fennoscandian marine and landlocked ringed seal subspecies

Island populations are on average smaller, genetically less diverse, and at a higher risk to go extinct than mainland populations. Low genetic diversity may elevate extinction probability, but the genetic component of the risk can be affected by the mode of diversity loss, which, in turn, is connected to the demographic history of the population. Here, we examined the history of genetic erosion in three Fennoscandian ringed seal subspecies, of which one inhabits the Baltic Sea 'mainland' and two the 'aquatic islands' composed of Lake Saimaa in Finland and Lake Ladoga in Russia. Both lakes were colonized by marine seals after their formation c. 9500 years ago, but Lake Ladoga is larger and more contiguous than Lake Saimaa. All three populations suffered dramatic declines during the 20th century, but the bottleneck was particularly severe in Lake Saimaa. Data from 17 microsatellite loci and mitochondrial control-region sequences show that Saimaa ringed seals have lost most of the genetic diversity present in their Baltic ancestors, while the Ladoga population has experienced only minor reductions. Using Approximate Bayesian computing analyses, we show that the genetic uniformity of the Saimaa subspecies derives from an extended founder event and subsequent slow erosion, rather than from the recent bottleneck. This suggests that the population has persisted for nearly 10,000 years despite having low genetic variation. The relatively high diversity of the Ladoga population appears to result from a high number of initial colonizers and a high post-colonization population size, but possibly also by a shorter isolation period and/or occasional gene flow from the Baltic Sea.

Museum DNA reveals the demographic history of the endangered Seychelles warbler

The importance of evolutionary conservation – how understanding evolutionary forces can help guide conservation decisions – is widely recognized. However, the historical demography of many endangered species is unknown, despite the fact that this can have important implications for contemporary ecological processes and for extinction risk. Here, we reconstruct the population history of the Seychelles warbler (Acrocephalus sechellensis) – an ecological model species. By the 1960s, this species was on the brink of extinction, but its previous history is unknown. We used DNA samples from contemporary and museum specimens spanning 140 years to reconstruct bottleneck history. We found a 25% reduction in genetic diversity between museum and contemporary populations, and strong genetic structure. Simulations indicate that the Seychelles warbler was bottlenecked from a large population, with an ancestral N_e of several thousands falling to <50 within the last century. Such a rapid decline, due to anthropogenic factors, has important implications for extinction risk in the Seychelles warbler, and our results will inform conservation practices. Reconstructing the population history of this species also allows us to better understand patterns of genetic diversity, inbreeding and promiscuity in the contemporary populations. Our approaches can be applied across species to test ecological hypotheses and inform conservation.

Global population structure and demographic history of the grey seal

Although the grey seal Halichoerus grypus is one of the most familiar and intensively studied of all pinniped species, its global population structure remains to be elucidated. Little is also known about how the species as a whole may have historically responded to climate-driven changes in habitat availability and anthropogenic exploitation. We therefore analysed samples from over 1500 individuals collected from 22 colonies spanning the Western and Eastern Atlantic and the Baltic Sea regions, represented by 350 bp of the mitochondrial hypervariable region and up to nine microsatellites. Strong population structure was observed at both types of marker, and highly asymmetrical patterns of gene flow were also inferred, with the Orkney Islands being identified as a source of emigrants to other areas in the Eastern Atlantic. The Baltic and Eastern Atlantic regions were estimated to have diverged a little over 10 000 years ago, consistent with the last proposed isolation of the Baltic Sea. Approximate Bayesian computation also identified genetic signals consistent with postglacial population expansion across much of the species range, suggesting that grey seals are highly responsive to changes in habitat availability.

Mixed fortunes: ancient expansion and recent decline in population size of a subtropical montane primate, the Arunachal macaque Macaca munzala

Quaternary glacial oscillations are known to have caused population size fluctuations in many temperate species. Species from subtropical and tropical regions are, however, considerably less studied, despite representing most of the biodiversity hotspots in the world including many highly threatened by anthropogenic activities such as hunting. These regions, consequently, pose a significant knowledge gap in terms of how their fauna have typically responded to past climatic changes. We studied an endangered primate, the Arunachal macaque Macaca munzala, from the subtropical southern edge of the Tibetan plateau, a part of the Eastern Himalaya biodiversity hotspot, also known to be highly threatened due to rampant hunting. We employed a 534 bp-long mitochondrial DNA sequence and 22 autosomal microsatellite loci to investigate the factors that have potentially shaped the demographic history of the species. Analysing the genetic data with traditional statistical methods and advance Bayesian inferential approaches, we demonstrate a limited effect of past glacial fluctuations on the demographic history of the species before the last glacial maximum, approximately 20,000 years ago. This was, however, immediately followed by a significant population expansion possibly due to warmer climatic conditions, approximately 15,000 years ago. These changes may thus represent an apparent balance between that displayed by the relatively climatically stable tropics and those of the more severe, temperate environments of the past. This study also draws attention to the possibility that a cold-tolerant species like the Arunachal macaque, which could withstand historical climate fluctuations and grow once the climate became conducive, may actually be extremely vulnerable to anthropogenic exploitation, as is perhaps indicated by its Holocene ca. 30-fold population decline, approximately 3,500 years ago. Our study thus provides a quantitative appraisal of these demographically important events, emphasising the ability to potentially infer the occurrence of two separate historical events from contemporary genetic data.

Detecting past changes of effective population size

Understanding and predicting population abundance is a major challenge confronting scientists. Several genetic models have been developed using microsatellite markers to estimate the present and ancestral effective population sizes. However, to get an overview on the evolution of population requires that past fluctuation of population size be traceable. To address the question, we developed a new model estimating the past changes of effective population size from microsatellite by resolving coalescence theory and using approximate likelihoods in a Monte Carlo Markov Chain approach. The efficiency of the model and its sensitivity to gene flow and to assumptions on the mutational process were checked using simulated data and analysis. The model was found especially useful to provide evidence of transient changes of population size in the past. The times at which some past demographic events cannot be detected because they are too ancient and the risk that gene flow may suggest the false detection of a bottleneck are discussed considering the distribution of coalescence times. The method was applied on real data sets from several Atlantic salmon populations. The method called VarEff (Variation of Effective size) was implemented in the R package VarEff and is made available at https://qgsp.jouy.inra.fr and at http://cran.r-project.org/web/packages/VarEff.

High genetic diversity in a small population: the case of Chilean blue whales

It is generally assumed that species with low population sizes have lower genetic diversities than larger populations and vice versa. However, this would not be the case for long-lived species with long generation times, and which populations have declined due to anthropogenic effects, such as the blue whale (Balaenoptera musculus). This species was intensively decimated globally to near extinction during the 20th century. Along the Chilean coast, it is estimated that at least 4288 blue whales were hunted from an apparently pre-exploitation population size (k) of a maximum of 6200 individuals (Southeastern Pacific). Thus, here, we describe the mtDNA (control region) and nDNA (microsatellites) diversities of the Chilean blue whale aggregation site in order to verify the expectation of low genetic diversity in small populations. We then compare our findings with other blue whale aggregations in the Southern Hemisphere. Interestingly, although the estimated population size is small compared with the pre-whaling era, there is still considerable genetic diversity, even after the population crash, both in mitochondrial (N = 46) and nuclear (N = 52) markers (Hd = 0.890 and Ho = 0.692, respectively). Our results suggest that this diversity could be a consequence of the long generation times and the relatively short period of time elapsed since the end of whaling, which has been observed in other heavily-exploited whale populations. The genetic variability of blue whales on their southern Chile feeding grounds was similar to that found in other Southern Hemisphere blue whale feeding grounds. Our phylogenetic analysis of mtDNA haplotypes does not show extensive differentiation of populations among Southern Hemisphere blue whale feeding grounds. The present study suggests that although levels of genetic diversity are frequently used as estimators of population health, these parameters depend on the biology of the species and should be taken into account in a monitoring framework study to obtain a more complete picture of the conservation status of a population.

Demographic history of a recent invasion of house mice on the isolated Island of Gough

Island populations provide natural laboratories for studying key contributors to evolutionary change, including natural selection, population size and the colonization of new environments. The demographic histories of island populations can be reconstructed from patterns of genetic diversity. House mice (Mus musculus) inhabit islands throughout the globe, making them an attractive system for studying island colonization from a genetic perspective. Gough Island, in the central South Atlantic Ocean, is one of the remotest islands in the world. House mice were introduced to Gough Island by sealers during the 19th century and display unusual phenotypes, including exceptionally large body size and carnivorous feeding behaviour. We describe genetic variation in Gough Island mice using mitochondrial sequences, nuclear sequences and microsatellites. Phylogenetic analysis of mitochondrial sequences suggested that Gough Island mice belong to Mus musculus domesticus, with the maternal lineage possibly originating in England or France. Cluster analyses of microsatellites revealed genetic membership for Gough Island mice in multiple coastal populations in Western Europe, suggesting admixed ancestry. Gough Island mice showed substantial reductions in mitochondrial and nuclear sequence variation and weak reductions in microsatellite diversity compared with Western European populations, consistent with a population bottleneck. Approximate Bayesian computation (ABC) estimated that mice recently colonized Gough Island (~100 years ago) and experienced a 98% reduction in population size followed by a rapid expansion. Our results indicate that the unusual phenotypes of Gough Island mice evolved rapidly, positioning these mice as useful models for understanding rapid phenotypic evolution.

Greece: a Balkan subrefuge for a remnant red deer (cervus elaphus) population

A number of phylogeographic studies have revealed the existence of multiple ice age refugia within the Balkan Peninsula, marking it as a biodiversity hotspot. Greece has been reported to harbor genetically differentiated lineages from the rest of Balkans for a number of mammal species. We therefore searched for distinct red deer lineages in Greece, by analyzing 78 samples originating from its last population in Parnitha Mountain (Central Greece). Additionally, we tested the impact of human-induced practices on this population. The presence of 2 discrete mtDNA lineages was inferred: 1) an abundant one not previously sampled in the Balkans and 2) a more restricted one shared with other Balkan populations, possibly the result of successful translocations of Eastern European individuals. Microsatellite-based analyses of 14 loci strongly support the existence of 2 subpopulations with relative frequencies similar to mitochondrial analyses. This study stresses the biogeographic importance of Central Greece as a separate Last Glacial Maximum period refugium within the Balkans. It also delineates the possible effects that recent translocations of red deer populations had on the genetic structuring within Parnitha. We suggest that the Greek red deer population of Parnitha is genetically distinct, and restocking programs should take this genetic evidence into consideration.

Laboratory colonisation and genetic bottlenecks in the tsetse fly Glossina pallidipes

Background

The IAEA colony is the only one available for mass rearing of Glossina pallidipes, a vector of human and animal African trypanosomiasis in eastern Africa. This colony is the source for Sterile Insect Technique (SIT) programs in East Africa. The source population of this colony is unclear and its genetic diversity has not previously been evaluated and compared to field populations.

Methodology/Principal Findings

We examined the genetic variation within and between the IAEA colony and its potential source populations in north Zimbabwe and the Kenya/Uganda border at 9 microsatellites loci to retrace the demographic history of the IAEA colony. We performed classical population genetics analyses and also combined historical and genetic data in a quantitative analysis using Approximate Bayesian Computation (ABC). There is no evidence of introgression from the north Zimbabwean population into the IAEA colony. Moreover, the ABC analyses revealed that the foundation and establishment of the colony was associated with a genetic bottleneck that has resulted in a loss of 35.7% of alleles and 54% of expected heterozygosity compared to its source population. Also, we show that tsetse control carried out in the 1990's is likely reduced the effective population size of the Kenya/Uganda border population.

Conclusions/Significance

All the analyses indicate that the area of origin of the IAEA colony is the Kenya/Uganda border and that a genetic bottleneck was associated with the foundation and establishment of the colony. Genetic diversity associated with traits that are important for SIT may potentially have been lost during this genetic bottleneck which could lead to a suboptimal competitiveness of the colony males in the field. The genetic diversity of the colony is lower than that of field populations and so, studies using colony flies should be interpreted with caution when drawing general conclusions about G. pallidipes biology.

There is only one mass reared laboratory colony of Glossina pallidipes, a vector of human African trypanosomiasis and arguably the main vector of animal African trypanosomiasis in eastern Africa. This colony is the main one used for basic research on this species and is intended to be used for Sterile Insect Technique (SIT) programs for control of field populations. The origins of this colony are not clear and the extent to which it is genetically representative of the species is unknown. Using population genetics analyses to compare the current colony with two potential source populations we have shown that the colony is from the Kenya/Uganda border and that its foundation and establishment in the laboratory were associated with a genetic bottleneck, i.e. reduction of genetic variation due to increased genetic drift in a population of reduced size. As a consequence, the genetic diversity of the colony is lower than that of G. pallidipes field populations.