Conservation genetics of population bottlenecks: the role of chance, selection, and history

Determining evolutionary origin and phylogenetic relationships of mallard-like ducks of Oceania, greater Indonesia, and the Philippines with ddRAD-seq data

AIM

We aim to determine the evolutionary origins and population genetics of mallard-like ducks of Oceania, greater Indonesia, and the Philippines.

LOCATION

Oceania, greater Indonesia, and the Philippines.

TAXON

Mallard (Anas platyrhynchos), Pacific black duck (A. superciliosa spp.), and Philippine duck (A. luzonica) METHODS: Thousands of nuclear ddRAD-seq loci and the mitochondrial DNA control region were assayed across individuals representative of each species' range. We assessed population structure and phylogenetic relationships, as well as estimated demographic histories to reconstruct the biogeographical history of each species.

RESULTS

Philippine and Pacific black ducks represent unique genetic lineages that diverged from the mallard 1-2 million years ago. We find no support for the Philippine duck representing a hybrid species as once posited; however, their low levels of genetic diversity requires further attention. We find a lack of substructure among Philippine ducks. However, we found pronounced differentiation between subspecies of Pacific black ducks, especially between A. s. superciliosa from New Zealand and A. s. rogersi from Australia, Papua New Guinea, and Timor-Leste, Indonesia. Anas superciliosa pelewensis gave mixed results; individuals from the Solomon Islands were differentiated from the other subspecies, but those from the island of Aunu'u, American Samoa, were genetically more similar to A. s. rogersi than A. s. pelewensis samples from the Solomon Islands. Finally, we find limited evidence of interspecific gene flow at evolutionary scales, and mallard introgression among contemporary samples.

MAIN CONCLUSIONS

Mallard-like ducks radiated across Oceania, greater Indonesia, and the Philippines within the last 2 million years. Only the Pacific black duck showed unique sub-structuring that largely followed known sub-species ranges, except for A. s. pelewensis. We posit that the high interrelatedness among Solomon Island samples suggests that their genetic distinctiveness may simply be the result of high levels of genetic drift. In contrast, we conclude that mainland Australian Pacific black ducks were the most likely source for the recent colonization of American Samoa. As a result, our findings suggest that either the A. s. pelewensis subspecies designations and/or its geographical range may require re-evaluation. Continued re-evaluation of evolutionary and taxonomic relationships is necessary when attempting to reconstruct and understand biogeographical histories, with important implications towards any attempts to implement conservation strategies.

Predicting species invasiveness with genomic data: Is genomic offset related to establishment probability?

Predicting the risk of establishment and spread of populations outside their native range represents a major challenge in evolutionary biology. Various methods have recently been developed to estimate population (mal)adaptation to a new environment with genomic data via so-called Genomic Offset (GO) statistics. These approaches are particularly promising for studying invasive species but have still rarely been used in this context. Here, we evaluated the relationship between GO and the establishment probability of a population in a new environment using both in silico and empirical data. First, we designed invasion simulations to evaluate the ability to predict establishment probability of two GO computation methods (Geometric GO and Gradient Forest) under several conditions. Additionally, we aimed to evaluate the interpretability of absolute Geometric GO values, which theoretically represent the adaptive genetic distance between populations from distinct environments. Second, utilizing public empirical data from the crop pest species Bactrocera tryoni, a fruit fly native from Northern Australia, we computed GO between "source" populations and a diverse range of locations within invaded areas. This practical application of GO within the context of a biological invasion underscores its potential in providing insights and guiding recommendations for future invasion risk assessment. Overall, our results suggest that GO statistics represent good predictors of the establishment probability and may thus inform invasion risk, although the influence of several factors on prediction performance (e.g., propagule pressure or admixture) will need further investigation.

Fire-driven animal evolution in the Pyrocene

Fire regimes are a major agent of evolution in terrestrial animals. Changing fire regimes and the capacity for rapid evolution in wild animal populations suggests the potential for rapid, fire-driven adaptive animal evolution in the Pyrocene. Fire drives multiple modes of evolutionary change, including stabilizing, directional, disruptive, and fluctuating selection, and can strongly influence gene flow and genetic drift. Ongoing and future research in fire-driven animal evolution will benefit from further development of generalizable hypotheses, studies conducted in highly responsive taxa, and linking fire-adapted phenotypes to their underlying genetic basis. A better understanding of evolutionary responses to fire has the potential to positively influence conservation strategies that embrace evolutionary resilience to fire in the Pyrocene.

Population demographic history and evolutionary rescue: Influence of a bottleneck event

Rapid environmental change presents a significant challenge to the persistence of natural populations. Rapid adaptation that increases population growth, enabling populations that declined following severe environmental change to grow and avoid extinction, is called evolutionary rescue. Numerous studies have shown that evolutionary rescue can indeed prevent extinction. Here, we extend those results by considering the demographic history of populations. To evaluate how demographic history influences evolutionary rescue, we created 80 populations of red flour beetle, Tribolium castaneum, with three classes of demographic history: diverse populations that did not experience a bottleneck, and populations that experienced either an intermediate or a strong bottleneck. We subjected these populations to a new and challenging environment for six discrete generations and tracked extinction and population size. Populations that did not experience a bottleneck in their demographic history avoided extinction entirely, while more than 20% of populations that experienced an intermediate or strong bottleneck went extinct. Similarly, among the extant populations at the end of the experiment, adaptation increased the growth rate in the novel environment the most for populations that had not experienced a bottleneck in their history. Taken together, these results highlight the importance of considering the demographic history of populations to make useful and effective conservation decisions and management strategies for populations experiencing environmental change that pushes them toward extinction.

The Relationship between Inbreeding and Fitness Is Different between Two Genetic Lines of European Bison

The pedigree data for European bison are published in the European Bison Pedigree Book, established one hundred years ago. The species passed a serious bottleneck and was kept in captivity for more than 30 years. After reintroduction, some individuals were captured and moved to enclosures, which caused gaps in pedigree data. To prevent the underestimation of inbreeding value, only animals with a fully known pedigree were used in the analysis. European bison were divided into two genetic lines, Lowland-Caucasian (LC) and Lowland (LB), with different numbers of founders (all 12 vs. 7 of them). The relation between inbreeding and survival up to one month of age, calculated using logistic regression, gave different results for every genetic line. In the LC line (N = 5441), the average inbreeding level was equal to 0.253 and the influence into survival was negative (exp(B) = 0.190), but in the LB line (N = 1227), the inbreeding level was much higher (0.410) but the influence into survival was positive (exp(B) = 6.596). It could be assumed that the difference between lines is a result of purging in the first period of species restitution.

Telemetry and Accelerometer Tracking of Green Toads in an Urban Habitat: Methodological Notes and Preliminary Findings

Advancements in tracking technologies provide an increasingly important tool in animal monitoring and conservation that can describe animal spatial behavior in native habitats and uncover migratory routes that otherwise may be difficult or impossible to map. In addition, high-resolution accelerometer sensors provide powerful insights into animal activity patterns and can help to identify specific behaviors from accelerometer profiles alone. Previously, such accelerometers were restricted to larger animals due to size and mass constraints. However, recent advances make it possible to use such devices on smaller animals such as the European green toad (Bufotes viridis), the focus of our current study. We deploy custom made tracking devices, that consist of very-high-frequency transmitters and tri-axial accelerometers, to track toads in their native urban environment in Vienna (Austria). A total of nine toads were tracked, ranging from three to nine tracking days per individual during the post-breeding season period. We demonstrate that our devices could reliably monitor toad movement and activity during the observation period. Hence, we confirmed the predominantly nocturnal activity patterns and recorded low overall movement at this urban site. Accelerometer data revealed that toads exhibited brief but intense activity bursts between 10 pm and midnight, resting periods during the night and intermittent activity during the day. Positional tracking alone would have missed the major activity events as they rarely resulted in large positional displacements. This underscores the importance of and value in integrating multiple tracking sensors for studies of movement ecology. Our approach could be adapted for other amphibians or other animals with mass constraints and may become standard monitoring equipment in the near future.

Conservation genomics of an endangered arboreal mammal following the 2019-2020 Australian megafire

The impacts of a changing climate threaten species, populations and ecosystems. Despite these significant and large-scale impacts on threatened species, many remain understudied and have little to no genetic information available. The greater glider, Petauroides volans, is an endangered species highly sensitive to the predicted changes in temperature under a changing climate and was recently severely impacted by a megafire natural disaster (85% estimated population loss). Baseline genetic data is essential for conservation management and for detecting detrimental changes in fire-effected populations. We collected genetic samples within 2 years post the 2019-2020 catastrophic Australian bushfires to examine adaptive potential, baseline genetic diversity and population structure, across their southern range in the state of New South Wales. Population genomic analyses were conducted using 8493 genome-wide SNPs for 86 greater glider individuals across 14 geographic locations. Substantial genetic structure was detected across locations, with low genetic diversity and effective population sizes observed in isolated areas. Additionally, we found signals of putative adaptation in response to temperature in greater gliders using a genotype-environment association analysis. These findings have important implications for the management of greater glider populations by identifying at-risk populations and identifying adaptive potential. We demonstrate the importance of baseline genetic information for endangered species as a practical approach to conservation. This is particularly important given the threat that changes in temperatures and megafire events, as predicted under a changing climate, poses for this species.

African wild dogs (Lycaon pictus) from the Kruger National Park, South Africa are currently not inbred but have low genomic diversity

African wild dogs (Lycaon pictus) have undergone severe population reductions and are listed as endangered on the International Union for Conservation of Nature Red List. Small, isolated populations have the potential to suffer from threats to their genetic diversity that may impact species viability and future survival. This study provides the first set of population-wide genomic data to address conservation concerns for this endangered species. Whole genome sequencing data were generated for 71 free-ranging African wild dogs from the Kruger National Park (KNP), South Africa, and used to estimate important population genomic parameters. Genomic diversity metrics revealed that variation levels were low; however, this African wild dog population showed low levels of inbreeding. Very few first- and second-order relationships were observed in this cohort, with most relationships falling into the third-order or distant category. Patterns of homozygosity could have resulted from historical inbreeding or a loss in genome variation due to a population bottleneck. Although the results suggest that this stronghold African wild dog population maintains low levels of inbreeding, likely due to their cooperative breeding system, it may lead to a continuous population decline when a reduced number of suitable mates are available. Consequently, the low genomic variation may influence species viability over time. This study highlights the importance of assessing population genomic parameters to set conservation priorities. Future studies should include the investigation of the potential of this endangered species to adapt to environmental changes considering the low genomic diversity in this population.

Characterization of polymorphic microsatellite markers for the Japanese endangered land snail Mandarina

Objective

Mandarina is an endangered land snail genus of the oceanic Ogasawara archipelago. On Chichijima Island, the largest inhabited island in Ogasawara, this genus is almost extinct in the wild due to predation by invasive species. Although ex situ conservation programs started in 2010, genetic diversity and population structure remain unclear due to a lack of genetic markers with sufficient genetic variation. In this study, we designed polymorphic microsatellite markers of Mandarina to enable genetic analysis and to develop appropriate conservation plans.

Results

Twenty-three polymorphic microsatellite markers were identified from the genomic DNA of wild samples of Mandarina mandarina. We assessed the genetic diversity of each marker. In 16 markers, neither linkage disequilibrium nor deviation from Hardy–Weinberg equilibrium was detected. These 16 markers were tested for multiplex PCR using low-density DNA extracted non-lethally from captive samples of M. mandarina, M. chichijimana and M. suenoae. Of the 16 markers, 15, 12 and 9 were usable for multiplex PCR, respectively. Genetic analysis using these microsatellite loci will be an important resource for the conservation of Mandarina.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13104-022-06147-4.

Food Restriction Reveals Individual Differences in Insulin-Like Growth Factor-1 Reaction Norms

Most organisms have to cope with unpredictable environmental challenges such as fluctuations in nutritional resources. Insulin-like growth factor-1 (IGF-1) is an evolutionarily conserved hormone that is highly sensitive to the individual nutritional status and regulates major life-history traits including lifespan and reproduction across vertebrates. We investigated the role of IGF-1 during periods of food shortages by altering between two feeding regimes (110 and 70% of daily food intake) after a period of ad libitum feeding in captive bearded reedlings (Panurus biarmicus). Each dietary treatment was repeated twice. Birds lost mass under food restriction, but the magnitude of mass change depended on the preceding dietary conditions. Moreover, bearded reedlings showed large, repeatable individual differences in their IGF-1 reaction norms with some individuals increasing IGF-1 levels in response to a restricted diet, whereas others showed no responses or decreased IGF-1 levels. This variation was explained by differences in average body mass: heavier individuals had higher IGF-1 levels during the control treatment and were more likely to decrease IGF-1 levels in response to the dietary restriction than did lighter ones. This result uncovers an individual by environment interaction (I × E) and may have important implications for the evolution of IGF-1 related hormonal phenotypes in this species.

Pinanga palms revisited 20 years on: what can changes in Pinanga species populations tell us about rainforest understory palm persistence?

Borneo is the centre of diversity of the palm genus Pinanga. At least 13 understory species have been recorded in the Ulu Temburong National Park in Brunei, but little is known of their persistence. Changes in populations of Pinanga understory palms may be indicative of more widespread changes due to climate change, such as changes in rainfall, which may be important for the palm diversity in the protected area. However, we know little about the population dynamics of these palms, how persistent their populations are or if they behave similarly over long time frames. In 1998, populations of five co-occurring species of Pinanga at several locations in the Ulu Temburong National Park were documented. This project aimed to undertake a comprehensive resurvey of the original five Pinanga palm species populations in order to assess if they showed similar population changes across sites and species after two decades. Overall, most species maintained their population size in the surveyed region but not consistently among sites, and one species significantly declined in abundance. There was considerable variation in population growth rate (R) within and among species and sites that was significantly correlated with density and the percentage of multi-stemmed plants. There was evidence of pulsed recruitment in some species and or sites rather than steady or exponential patterns of population growth.

Diversity of the MHC class II DRB gene in the wolverine (Carnivora: Mustelidae: Gulo gulo) in Finland

The wolverine (Gulo gulo) in Finland has undergone significant population declines in the past. Since major histocompatibility complex (MHC) genes encode proteins involved in pathogen recognition, the diversity of these genes provides insights into the immunological fitness of regional populations. We sequenced 862 amplicons (242 bp) of MHC class II DRB exon 2 from 32 Finnish wolverines and identified 11 functional alleles and three pseudogenes. A molecular phylogenetic analysis indicated trans-species polymorphism, and PAML and MEME analyses indicated positive selection, suggesting that the Finnish wolverine DRB genes have evolved under balancing and positive selection. In contrast to DRB gene analyses in other species, allele frequencies in the Finnish wolverines clearly indicated the existence of two regional subpopulations, congruent with previous studies based on neutral genetic markers. In the Finnish wolverine, rapid population declines in the past have promoted genetic drift, resulting in a lower genetic diversity of DRB loci, including fewer alleles and positively selected sites, than other mustelid species analyzed previously. Our data suggest that the MHC region in the Finnish wolverine population was likely affected by a recent bottleneck.

Sustained positive consequences of genetic rescue of fitness and behavioural traits in inbred populations of Drosophila melanogaster

One solution to alleviate the detrimental genetic effects associated with reductions in population size and fragmentation is to introduce immigrants from other populations. While the effects of this genetic rescue on fitness traits are fairly well known, it is less clear to what extent inbreeding depression and subsequent genetic rescue affect behavioural traits. In this study, replicated crosses between inbred lines of Drosophila melanogaster were performed in order to investigate the effects of inbreeding and genetic rescue on egg-to-adult viability and negative geotaxis behaviour-a locomotor response used to measure, e.g. the effects of physiological ageing. Transgenerational effects of outcrossing were investigated by examining the fitness consequences in both the F₁ and F₄  generation. The majority of inbred lines showed evidence for inbreeding depression for both egg-to-adult viability and behavioural performance (95% and 66% of lines, respectively), with inbreeding depression being more pronounced for viability compared with the locomotor response. Subsequent outcrossing with immigrants led to an alleviation of the negative effects for both viability and geotaxis response resulting in inbred lines being similar to the outbred controls, with beneficial effects persisting from F₁ to F₄ . Overall, the results clearly show that genetic rescue can provide transgenerational rescue of small, inbred populations by rapidly improving population fitness components. Thus, we show that even the negative effects of inbreeding on behaviour, similar to that of neurodegeneration associated with physiological ageing, can be reversed by genetic rescue.

Variation in intraspecific demography drives localised concordance but species-wide discordance in response to past climatic change

BACKGROUND

Understanding how species biology may facilitate resilience to climate change remains a critical factor in detecting and protecting species at risk of extinction. Many studies have focused on the role of particular ecological traits in driving species responses, but less so on demographic history and levels of standing genetic variation. Additionally, spatial variation in the interaction of demographic and adaptive factors may further complicate prediction of species responses to environmental change. We used environmental and genomic datasets to reconstruct the phylogeographic histories of two ecologically similar and largely co-distributed freshwater fishes, the southern (Nannoperca australis) and Yarra (N. obscura) pygmy perches, to assess the degree of concordance in their responses to Plio-Pleistocene climatic changes. We described contemporary genetic diversity, phylogenetic histories, demographic histories, and historical species distributions across both species, and statistically evaluated the degree of concordance in co-occurring populations.

RESULTS

Marked differences in contemporary genetic diversity, historical distribution changes and historical migration were observed across the species, with a distinct lack of genetic diversity and historical range expansion suggested for N. obscura. Although several co-occurring populations within a shared climatic refugium demonstrated concordant demographic histories, idiosyncratic population size changes were found at the range edges of the more spatially restricted species. Discordant responses between species were associated with low standing genetic variation in peripheral populations. This might have hindered adaptive potential, as documented in recent demographic declines and population extinctions for the two species.

CONCLUSION

Our results highlight both the role of spatial scale in the degree of concordance in species responses to climate change, and the importance of standing genetic variation in facilitating range shifts. Even when ecological traits are similar between species, long-term genetic diversity and historical population demography may lead to discordant responses to ongoing and future climate change.

Limiting factors for queen conch (Lobatus gigas) reproduction: A simulation-based evaluation

Queen conch are among the most economically, socially, and culturally important fishery resources in the Caribbean. Despite a multitude of fisheries management measures enacted across the region, populations are depleted and failing to recover. It is believed that queen conch are highly susceptible to depensatory processes, impacting reproductive success and contributing to the lack of recovery. We developed a model of reproductive dynamics to evaluate how variations in biological factors such as population density, movement speeds, rest periods between mating events, scent tracking, visual perception of conspecifics, sexual facilitation, and barriers to movement affect reproductive success and overall reproductive output. We compared simulation results to empirical observations of mating and spawning frequencies from conch populations in the central Bahamas and Florida Keys. Our results confirm that low probability of mate finding associated with decreased population density is the primary driver behind observed breeding behavior in the field, but is insufficient to explain observed trends. Specifically, sexual facilitation coupled with differences in movement speeds and ability to perceive conspecifics may explain the observed lack of mating at low densities and differences between mating frequencies in the central Bahamas and Florida Keys, respectively. Our simulations suggest that effective management strategies for queen conch should aim to protect high-density reproductive aggregations and critical breeding habitats.

A continent-wide high genetic load in African buffalo revealed by clines in the frequency of deleterious alleles, genetic hitchhiking and linkage disequilibrium

A high genetic load can negatively affect population viability and increase susceptibility to diseases and other environmental stressors. Prior microsatellite studies of two African buffalo (Syncerus caffer) populations in South Africa indicated substantial genome-wide genetic load due to high-frequency occurrence of deleterious alleles. The occurrence of these alleles, which negatively affect male body condition and bovine tuberculosis resistance, throughout most of the buffalo's range were evaluated in this study. Using available microsatellite data (2-17 microsatellite loci) for 1676 animals from 34 localities (from 25°S to 5°N), we uncovered continent-wide frequency clines of microsatellite alleles associated with the aforementioned male traits. Frequencies decreased over a south-to-north latitude range (average per-locus Pearson r = -0.22). The frequency clines coincided with a multilocus-heterozygosity cline (adjusted R2 = 0.84), showing up to a 16% decrease in southern Africa compared to East Africa. Furthermore, continent-wide linkage disequilibrium (LD) at five linked locus pairs was detected, characterized by a high fraction of positive interlocus associations (0.66, 95% CI: 0.53, 0.77) between male-deleterious-trait-associated alleles. Our findings suggest continent-wide and genome-wide selection of male-deleterious alleles driven by an earlier observed sex-chromosomal meiotic drive system, resulting in frequency clines, reduced heterozygosity due to hitchhiking effects and extensive LD due to male-deleterious alleles co-occurring in haplotypes. The selection pressures involved must be high to prevent destruction of allele-frequency clines and haplotypes by LD decay. Since most buffalo populations are stable, these results indicate that natural mammal populations, depending on their genetic background, can withstand a high genetic load.

Genetic purging in captive endangered ungulates with extremely low effective population sizes

Inbreeding threatens the survival of small populations by producing inbreeding depression, but also exposes recessive deleterious effects in homozygosis allowing for genetic purging. Using inbreeding-purging theory, we analyze early survival in four pedigreed captive breeding programs of endangered ungulates where population growth was prioritized so that most adult females were allowed to contribute offspring according to their fitness. We find evidence that purging can substantially reduce inbreeding depression in Gazella cuvieri (with effective population size N_e = 14) and Nanger dama (N_e = 11). No purging is detected in Ammotragus lervia (N_e = 4), in agreement with the notion that drift overcomes purging under fast inbreeding, nor in G. dorcas (N_e = 39) where, due to the larger population size, purging is slower and detection is expected to require more generations. Thus, although smaller populations are always expected to show smaller fitness (as well as less adaptive potential) than larger ones due to higher homozygosis and deleterious fixation, our results show that a substantial fraction of their inbreeding load and inbreeding depression can be purged when breeding contributions are governed by natural selection. Since management strategies intended to maximize the ratio from the effective to the actual population size tend to reduce purging, the search for a compromise between these strategies and purging could be beneficial in the long term. This could be achieved either by allowing some level of random mating and some role of natural selection in determining breeding contributions, or by undertaking reintroductions into the wild at the earliest opportunity.

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.

Whole-genome sequencing of Schistosoma mansoni reveals extensive diversity with limited selection despite mass drug administration

Control and elimination of the parasitic disease schistosomiasis relies on mass administration of praziquantel. Whilst these programmes reduce infection prevalence and intensity, their impact on parasite transmission and evolution is poorly understood. Here we examine the genomic impact of repeated mass drug administration on Schistosoma mansoni populations with documented reduced praziquantel efficacy. We sequenced whole-genomes of 198 S. mansoni larvae from 34 Ugandan children from regions with contrasting praziquantel exposure. Parasites infecting children from Lake Victoria, a transmission hotspot, form a diverse panmictic population. A single round of treatment did not reduce this diversity with no apparent population contraction caused by long-term praziquantel use. We find evidence of positive selection acting on members of gene families previously implicated in praziquantel action, but detect no high frequency functionally impactful variants. As efforts to eliminate schistosomiasis intensify, our study provides a foundation for genomic surveillance of this major human parasite.

Wildlife conservation in a fragmented landscape: the Eurasian red squirrel on the Isle of Wight

Island populations may have a higher extinction risk due to reduced genetic diversity and need to be managed effectively in order to reduce the risk of biodiversity loss. The Eurasian red squirrels (Sciurus vulgaris) in the south of England only survive on three islands (the Isle of Wight, Brownsea and Furzey islands), with the Isle of Wight harbouring the largest population in the region. Fourteen microsatellites were used to determine the genetic structure of red squirrel populations on the Isle of Wight, as well as their relatedness to other populations of the species. Our results demonstrated that squirrels on these islands were less genetically diverse than those in Continental mainland populations, as would be expected. It also confirmed previous results from mitochondrial DNA which indicated that the squirrels on the Isle of Wight were relatively closely related to Brownsea island squirrels in the south of England. Importantly, our findings showed that genetic mixing between squirrels in the east and west of the Isle of Wight was very limited. Given the potential deleterious effects of small population size on genetic health, landscape management to encourage dispersal of squirrels between these populations should be a priority.

Demographic and ecogeographic factors limit wild grapevine spread at the southern edge of its distribution range

The spatial distribution of plants is constrained by demographic and ecogeographic factors that determine the range and abundance of the species. Wild grapevine (Vitis vinifera ssp. sylvestris) is distributed from Switzerland in the north to Israel in the south. However, little is known about the ecogeographic constraints of this species and its genetic and phenotypic characteristics, especially at the southern edge of its distribution range in the Levant region. In this study, we explore the population structure of southern Levantine wild grapevines and the correlation between demographic and ecogeographic characteristics. Based on our genetic analysis, the wild grapevine populations in this region can be divided into two major subgroups in accordance with a multivariate spatial and ecogeographical clustering model. The identified subpopulations also differ in morphological traits, mainly leaf hairiness which may imply adaptation to environmental stress. The findings suggest that the Upper Jordan River population was spread to the Sea of Galilee area and that a third smaller subpopulation at the south of the Golan Heights may represent a distinguished gene pool or a recent establishment of a new population. A spatial distribution model indicated that distance to water sources, Normalized difference vegetation index, and precipitation are the main environmental factors constraining V. v. sylvestris distribution at its southern distribution range. These factors in addition to limited gene flow between populations prevent further spread of wild grapevines southwards to semi-arid regions.

Host Genetic Diversity and Infectious Diseases. Focus on Wild Boar, Red Deer and Tuberculosis

Host genetic diversity tends to limit disease spread in nature and buffers populations against epidemics. Genetic diversity in wildlife is expected to receive increasing attention in contexts related to disease transmission and human health. Ungulates such as wild boar (Sus scrofa) and red deer (Cervus elaphus) are important zoonotic hosts that can be precursors to disease emergence and spread in humans. Tuberculosis is a zoonotic disease with relevant consequences and can present high prevalence in wild boar and red deer populations. Here, we review studies on the genetic diversity of ungulates and determine to what extent these studies consider its importance on the spread of disease. This assessment also focused on wild boar, red deer, and tuberculosis. We found a disconnection between studies treating genetic diversity and those dealing with infectious diseases. Contrarily, genetic diversity studies in ungulates are mainly concerned with conservation. Despite the existing disconnection between studies on genetic diversity and studies on disease emergence and spread, the knowledge gathered in each discipline can be applied to the other. The bidirectional applications are illustrated in wild boar and red deer populations from Spain, where TB is an important threat for wildlife, livestock, and humans.

Salmonella Infection in Turtles: A Risk for Staff Involved in Wildlife Management?

Simple Summary

The aim of this study was to investigate the occurrence of non-typhoidal Salmonella in the turtles housed in a regional wildlife rescue centre of Apulia, in southern Italy, to assess the presence of Salmonella serovars that may represent a risk for operators involved in wildlife management. Sixty-nine tortoises, of which 36 were males and 33 were females, belonging to different species (Testudo hermanni hermanni, T. h. boettgeri, T. graeca, and T. marginata) were tested. All the turtles were adults (34 between 6 and 10 years of age and 35 more than 10 years of age). Salmonella was statistically detected more frequently in T. hermanni hermanni. No differences of the infection prevalence related to animal gender or age were found. Two different species, S. enterica and S. bongori, three S. enterica subspecies (enterica, diarizonae, salamae), and five different serovars (Hermannswerder, Abony, Ferruch, Richmond, Vancouver) within the group S. enterica subspecies enterica were identified. Two Salmonella types with different combinations were simultaneously found in specimens of T. h. hermanni. Most of the detected Salmonella types may represent a potential risk for operators in wildlife rescue centres.

Abstract

Monitoring of infections that may be transmitted to humans by animals in wildlife rescue centres is very important in order to protect the staff engaged in rehabilitation practices. Salmonella may be a natural inhabitant of the intestinal tract of turtles, rarely causing disease. This may represent a potential risk for humans, increasing the sanitary risk for operators in wildlife rescue centres. In this paper, the occurrence of non-typhoidal Salmonella among terrestrial turtles housed in a wildlife rescue centre in Southern Italy was investigated, in order to assess the serovars more frequently carried by turtles and identify those that may represent a risk for operators involved in wildlife management. Sixty-nine adult turtles (Testudo hermanni hermanni, T. h. boettgeri, T. graeca, and T. marginata) were tested. Detection and serotyping of Salmonella strains were performed according to ISO 6579-1 and ISO/TR 6579-3:2013, respectively. The distribution of Salmonella spp. was significantly higher in T. hermanni hermanni than in other species, independent of the age and gender of the animals. Two different Salmonella species, S. enterica and S. bongori, three S. enterica subspecies (enterica, diarizonae, salamae), and five different serovars (Hermannswerder, Abony, Ferruch, Richmond, Vancouver) within the group S. enterica subspecies enterica were identified. Different combinations of Salmonella types were simultaneously found in specimens of T. h. hermanni. Most of detected Salmonella types may represent a potential risk for public health. Adopting correct animal husbandry procedures and informing on potential sanitary risks may be useful for minimising the risk of transmission of Salmonella to workers involved in wildlife management.

Genetic signatures of population bottlenecks, relatedness, and inbreeding highlight recent and novel conservation concerns in the Egyptian vulture

The assessment of temporal variation in genetic features can be particularly informative on the factors behind demography and viability of wildlife populations and species. We used molecular methods to evaluate neutral genetic variation, relatedness, bottlenecks, and inbreeding in a declining population of Egyptian vulture (Neophron percnopterus) in central Spain. The results show that the genetic diversity remained relatively stable over a period of twelve years despite the decline in census and effective population sizes in the last decades. A relatively high proportion of nestlings from different and distant territories showed high relatedness in each study year. We also found support for an increasing impact of severe recent (contemporary) rather than distant (historical) past demographic bottlenecks, and the first evidence of inbred mating between full siblings coinciding with lethal malformations in offspring. The inbred nestling with feather malformations was positive to beak and feather disease virus recorded for the first time in this species. These results alert on recent and novel threats potentially affecting health and reducing the adaptive potential of individuals in this threatened species.

An isolated white-tailed deer (Odocoileus virginianus) population on St. John, US Virgin Islands shows low inbreeding and comparable heterozygosity to other larger populations

This is the first study to document the genetic diversity of the white-tailed deer population on St. John, US Virgin Islands. The island population was founded by a small number of animals, has very limited hunting or predation, and recently experienced a reduction in size following an extended drought in 2015. DNA samples were collected from hair from 23 anesthetized adult deer (13 males, 10 females) ranging in age from 1 to 8 years (3.36 ± 1.9 years) and also from fecal DNA samples, for a total of 42 individuals analyzed for genetic diversity. The St. John deer data set averaged 4.19 alleles per marker and demonstrates the second lowest number of alleles (A) when compared to other populations of Odocoileus virginianus (4.19). Heterozygosity was similar to the other studies (0.54) with little evidence of inbreeding. To explain the level of heterozygosity and level of inbreeding within the St. John population, three hypotheses are proposed, including the effect of intrinsic biological traits within the population, a recent infusion of highly heterogeneous loci from North American populations, and a consistent level of immigration from a nearby island. Additional work is needed to further understand the genetic history of the St. John and regional deer populations.

Polyandry provides reproductive and genetic benefits in colonising populations

Polyandry, when females mate with more than one male, is theorised to play an important role in successful colonisation of new habitats. In addition to possible benefits from sexual selection, even mild polyandry could facilitate colonisation by protecting against inbreeding and reducing the costs of mating with incompatible or infertile males. Here, we measure the importance of mild polyandry for population viability and reproductive fitness following experimental founder events into a higher‐temperature regime. Using colonisation experiments with the model beetle Tribolium castaneum, in which females can produce offspring for up to 140 days following a single mating, we founded more than 100 replicate populations using single females that had been given the opportunity to mate with either one or two males and then tracked their subsequent population dynamics. Following population viability and fitness across 10 generations, we found that extinction rates were significantly lower in populations founded by females given polyandrous opportunities to mate with two males (9%) compared to populations founded by monogamous females (34%). In addition, populations founded by females that had been provided with opportunities to store sperm from two different males showed double the median productivity following colonisation compared to monogamous‐founded populations. Notably, we identified short‐term and longer‐term benefits to post‐colonisation populations from double‐mating, with results suggesting that polyandry acts to both protect against mating with incompatible males through the founder event, and reduce inbreeding depression as the colonisation proceeds for 10 generations. Our results therefore show that even mild polyandry provides both reproductive and genetic benefits for colonising populations.

Polygenic basis for adaptive morphological variation in a threatened Aotearoa | New Zealand bird, the hihi (Notiomystis cincta)

To predict if a threatened species can adapt to changing selective pressures, it is crucial to understand the genetic basis of adaptive traits, especially in species historically affected by severe bottlenecks. We estimated the heritability of three hihi (Notiomystis cincta) morphological traits known to be under selection (nestling tarsus length, body mass and head-bill length) using 523 individuals and 39 699 single nucleotide polymorphisms (SNPs) from a 50 K Affymetrix SNP chip. We then examined the genetic architecture of the traits via chromosome partitioning analyses and genome-wide association scans (GWAS). Heritabilities estimated using pedigree relatedness or genomic relatedness were low. For tarsus length, the proportion of genetic variance explained by each chromosome was positively correlated with its size, and more than one chromosome explained significant variation for body mass and head-bill length. Finally, GWAS analyses suggested many loci of small effect contributing to trait variation for all three traits, although one locus (an SNP within an intron of the transcription factor HEY2) was tentatively associated with tarsus length. Our findings suggest a polygenic nature for the morphological traits, with many small effect size loci contributing to the majority of the variation, similar to results from many other wild populations. However, the small effective population size, polygenic architecture and already low heritabilities suggest that both the total response and rate of response to selection are likely to be limited in hihi.

Survival of the cheapest: how proteome cost minimization drives evolution

Darwin's theory of evolution emphasized that positive selection of functional proficiency provides the fitness that ultimately determines the structure of life, a view that has dominated biochemical thinking of enzymes as perfectly optimized for their specific functions. The 20th-century modern synthesis, structural biology, and the central dogma explained the machinery of evolution, and nearly neutral theory explained how selection competes with random fixation dynamics that produce molecular clocks essential e.g. for dating evolutionary histories. However, quantitative proteomics revealed that selection pressures not relating to optimal function play much larger roles than previously thought, acting perhaps most importantly via protein expression levels. This paper first summarizes recent progress in the 21st century toward recovering this universal selection pressure. Then, the paper argues that proteome cost minimization is the dominant, underlying 'non-function' selection pressure controlling most of the evolution of already functionally adapted living systems. A theory of proteome cost minimization is described and argued to have consequences for understanding evolutionary trade-offs, aging, cancer, and neurodegenerative protein-misfolding diseases.

Museomics identifies genetic erosion in two butterfly species across the 20th century in Finland

Erosion of biodiversity generated by anthropogenic activities has been studied for decades and in many areas at the species level, using taxa monitoring. In contrast, genetic erosion within species has rarely been tracked, and is often studied by inferring past population dynamics from contemporaneous estimators. An alternative to such inferences is the direct examination of past genes, by analysing museum collection specimens. While providing direct access to genetic variation over time, historical DNA is usually not optimally preserved, and it is necessary to apply genotyping methods based on hybridization-capture to unravel past genetic variation. In this study, we apply such a method (i.e., HyRAD), to large time series of two butterfly species in Finland, and present a new bioinformatic pipeline, namely PopHyRAD, that standardizes and optimizes the analysis of HyRAD data at the within-species level. In the localities for which the data retrieved have sufficient power to accurately examine genetic dynamics through time, we show that genetic erosion has increased across the last 100 years, as revealed by signatures of allele extinctions and heterozygosity decreases, despite local variations. In one of the two butterflies (Erebia embla), isolation by distance also increased through time, revealing the effect of greater habitat fragmentation over time.

Adaptive genetic diversity and evidence of population genetic structure in the endangered Sierra Madre Sparrow (Xenospiza baileyi)

The magnitude and distribution of genetic diversity through space and time can provide useful information relating to evolutionary potential and conservation status in threatened species. In assessing genetic diversity in species that are of conservation concern, several studies have focused on the use of Toll-like receptors (TLRs). TLRs are innate immune genes related to pathogen resistance, and polymorphisms may reflect not only levels of functional diversity, but may also be used to assess genetic diversity within and among populations. Here, we combined four potentially adaptive markers (TLRs) with one mitochondrial (COI) marker to evaluate genetic variation in the endangered Sierra Madre Sparrow (Xenospiza baileyi). This species offers an ideal model to investigate population and evolutionary genetic processes that may be occurring in a habitat restricted endangered species with disjunct populations (Mexico City and Durango), the census sizes of which differ by an order of magnitude. TLRs diversity in the Sierra Madre Sparrow was relatively high, which was not expected given its two small, geographically isolated populations. Genetic diversity was different (but not significantly so) between the two populations, with less diversity seen in the smaller Durango population. Population genetic structure between populations was due to isolation and different selective forces acting on different TLRs; population structure was also evident in COI. Reduction of genetic diversity in COI was observed over 20 years in the Durango population, a result likely caused by habitat loss, a factor which may be the main cause of diversity decline generally. Our results provide information related to the ways in which adaptive variation can be altered by demographic changes due to human-mediated habitat alterations. Furthermore, our findings may help to guide conservation schemes for both populations and their restricted habitat.

Severe population bottleneck and cranial morphology change in the Mednyi Island subspecies of Arctic fox Vulpes lagopus (Carnivora: Canidae)

Arctic foxes, Vulpes lagopus living on Mednyi Island suffered a drastic decline in population size in the late 1970s due to an outbreak of mange epizootic. This dramatic fall in numbers rendered the subspecies endangered, and the concomitant loss of variability resulted in a population bottleneck. Here, we investigate whether differences in cranial morphology between Mednyi Island Arctic foxes and Bering Island Arctic foxes could be attributed to the severe population bottleneck suffered by the Mednyi population in the 1970s. We used morphometric traits as proxies for genetic data to provide estimates of FST. Results show higher FST estimates for the Mednyi population than for the Bering population, which we interpret as a bottleneck signature. FST results also indicate a pattern of divergence between the two populations consistent with random genetic drift. Bottleneck detection is critical for the interpretation of the demographic history of the endangered Mednyi Island Arctic fox, with consequences for conservation management.

Population modelling and genetics of a critically endangered Madagascan palm Tahina spectabilis

Madagascar is home to 208 indigenous palm species, almost all of them endemic and >80% of which are endangered. We undertook complete population census and sampling for genetic analysis of a relatively recently discovered giant fan palm, the Critically Endangered Tahina spectablis in 2008 and 2016. Our 2016 study included newly discovered populations and added to our genetic study. We incorporated these new populations into species distribution niche model (SDM) and projected these onto maps of the region. We developed population matrix models based on observed demographic data to model population change and predict the species vulnerability to extinction by undertaking population viability analysis (PVA). We investigated the potential conservation value of reintroduced planted populations within the species potential suitable habitat. We found that the population studied in 2008 had grown in size due to seedling regeneration but had declined in the number of reproductively mature plants, and we were able to estimate that the species reproduces and dies after approximately 70 years. Our models suggest that if the habitat where it resides continues to be protected the species is unlikely to go extinct due to inherent population decline and that it will likely experience significant population growth after approximately 80 years due to the reproductive and life cycle attributes of the species. The newly discovered populations contain more genetic diversity than the first discovered southern population which is genetically depauperate. The species appears to demonstrate a pattern of dispersal leading to isolated founder plants which may eventually lead to population development depending on local establishment opportunities. The conservation efforts currently put in place including the reintroduction of plants within the species potential suitable habitat if maintained are thought likely to enable the species to sustain itself but it remains vulnerable to anthropogenic impacts.

The role of demographic history and selection in shaping genetic diversity of the Galápagos penguin (Spheniscus mendiculus)

Although many studies have documented the effects of demographic bottlenecks on the genetic diversity of natural populations, there is conflicting evidence of the roles that genetic drift and selection may play in driving changes in genetic variation at adaptive loci. We analyzed genetic variation at microsatellite and mitochondrial loci in conjunction with an adaptive MHC class II locus in the Galápagos penguin (Spheniscus mendiculus), a species that has undergone serial demographic bottlenecks associated with El Niño events through its evolutionary history. We compared levels of variation in the Galápagos penguin to those of its congener, the Magellanic penguin (Spheniscus magellanicus), which has consistently maintained a large population size and thus was used as a non-bottlenecked control. The comparison of neutral and adaptive markers in these two demographically distinct species allowed assessment of the potential role of balancing selection in maintaining levels of MHC variation during bottleneck events. Our analysis suggests that the lack of genetic diversity at both neutral and adaptive loci in the Galápagos penguin likely resulted from its restricted range, relatively low abundance, and history of demographic bottlenecks. The Galápagos penguin revealed two MHC alleles, one mitochondrial haplotype, and six alleles across five microsatellite loci, which represents only a small fraction of the diversity detected in Magellanic penguins. Despite the decreased genetic diversity in the Galápagos penguin, results revealed signals of balancing selection at the MHC, which suggest that selection can mitigate some of the effects of genetic drift during bottleneck events. Although Galápagos penguin populations have persisted for a long time, increased frequency of El Niño events due to global climate change, as well as the low diversity exhibited at immunological loci, may put this species at further risk of extinction.

Comparative phylogeography of mainland and insular species of Neotropical molossid bats (Molossus)

Historical events, habitat preferences, and geographic barriers might result in distinct genetic patterns in insular versus mainland populations. Comparison between these two biogeographic systems provides an opportunity to investigate the relative role of isolation in phylogeographic patterns and to elucidate the importance of evolution and demographic history in population structure. Herein, we use a genotype-by-sequencing approach (GBS) to explore population structure within three species of mastiff bats (Molossus molossus, M. coibensis, and M. milleri), which represent different ecological histories and geographical distributions in the genus. We tested the hypotheses that oceanic straits serve as barriers to dispersal in Caribbean bats and that isolated island populations are more likely to experience genetic drift and bottlenecks in comparison with highly connected ones, thus leading to different phylogeographic patterns. We show that population structures vary according to general habitat preferences, levels of population isolation, and historical fluctuations in climate. In our dataset, mainland geographic barriers played only a small role in isolation of lineages. However, oceanic straits posed a partial barrier to the dispersal for some populations within some species (M. milleri), but do not seem to disrupt gene flow in others (M. molossus). Lineages on distant islands undergo genetic bottlenecks more frequently than island lineages closer to the mainland, which have a greater exchange of haplotypes.

Genomic variation predicts adaptive evolutionary responses better than population bottleneck history

The relationship between population size, inbreeding, loss of genetic variation and evolutionary potential of fitness traits is still unresolved, and large-scale empirical studies testing theoretical expectations are surprisingly scarce. Here we present a highly replicated experimental evolution setup with 120 lines of Drosophila melanogaster having experienced inbreeding caused by low population size for a variable number of generations. Genetic variation in inbred lines and in outbred control lines was assessed by genotyping-by-sequencing (GBS) of pooled samples consisting of 15 males per line. All lines were reared on a novel stressful medium for 10 generations during which body mass, productivity, and extinctions were scored in each generation. In addition, we investigated egg-to-adult viability in the benign and the stressful environments before and after rearing at the stressful conditions for 10 generations. We found strong positive correlations between levels of genetic variation and evolutionary response in all investigated traits, and showed that genomic variation was more informative in predicting evolutionary responses than population history reflected by expected inbreeding levels. We also found that lines with lower genetic diversity were at greater risk of extinction. For viability, the results suggested a trade-off in the costs of adapting to the stressful environments when tested in a benign environment. This work presents convincing support for long-standing evolutionary theory, and it provides novel insights into the association between genetic variation and evolutionary capacity in a gradient of diversity rather than dichotomous inbred/outbred groups.

Spatial differences in genetic diversity and northward migration suggest genetic erosion along the boreal caribou southern range limit and continued range retraction

With increasing human activities and associated landscape changes, distributions of terrestrial mammals become fragmented. These changes in distribution are often associated with reduced population sizes and loss of genetic connectivity and diversity (i.e., genetic erosion) which may further diminish a species' ability to respond to changing environmental conditions and lead to local population extinctions. We studied threatened boreal caribou (Rangifer tarandus caribou) populations across their distribution in Ontario/Manitoba (Canada) to assess changes in genetic diversity and connectivity in areas of high and low anthropogenic activity. Using data from >1,000 caribou and nine microsatellite loci, we assessed population genetic structure, genetic diversity, and recent migration rates using a combination of network and population genetic analyses. We used Bayesian clustering analyses to identify population genetic structure and explored spatial and temporal variation in those patterns by assembling networks based on R ST and F ST as historical and contemporary genetic edge distances, respectively. The Bayesian clustering analyses identified broad-scale patterns of genetic structure and closely aligned with the R ST network. The F ST network revealed substantial contemporary genetic differentiation, particularly in areas presenting contemporary anthropogenic disturbances and habitat fragmentation. In general, relatively lower genetic diversity and greater genetic differentiation were detected along the southern range limit, differing from areas in the northern parts of the distribution. Moreover, estimation of migration rates suggested a northward movement of animals away from the southern range limit. The patterns of genetic erosion revealed in our study suggest ongoing range retraction of boreal caribou in central Canada.

The Genomic Footprints of the Fall and Recovery of the Crested Ibis

Human-induced environmental change and habitat fragmentation pose major threats to biodiversity and require active conservation efforts to mitigate their consequences. Genetic rescue through translocation and the introduction of variation into imperiled populations has been argued as a powerful means to preserve, or even increase, the genetic diversity and evolutionary potential of endangered species [1-4]. However, factors such as outbreeding depression [5, 6] and a reduction in available genetic diversity render the success of such approaches uncertain. An improved evaluation of the consequence of genetic restoration requires knowledge of temporal changes to genetic diversity before and after the advent of management programs. To provide such information, a growing number of studies have included small numbers of genomic loci extracted from historic and even ancient specimens [7, 8]. We extend this approach to its natural conclusion, by characterizing the complete genomic sequences of modern and historic population samples of the crested ibis (Nipponia nippon), an endangered bird that is perhaps the most successful example of how conservation effort has brought a species back from the brink of extinction. Though its once tiny population has today recovered to >2,000 individuals [9], this process was accompanied by almost half of ancestral loss of genetic variation and high deleterious mutation load. We furthermore show how genetic drift coupled to inbreeding following the population bottleneck has largely purged the ancient polymorphisms from the current population. In conclusion, we demonstrate the unique promise of exploiting genomic information held within museum samples for conservation and ecological research.

Genetic Diversity among Clouded Salamanders (Hynobius nebulosus) in Shiga Prefecture

The Japanese clouded salamander (Hynobius nebulosus) is a lentic-breeding species distributed throughout western Japan. Threats, such as habitat loss, have led to it being categorized as a vulnerable species. To explore the phylogeographic features and population differentiation among clouded salamanders in Shiga prefecture, we analyzed sequences of the mitochondrial cytochrome b gene. DNA samples were collected from 29 distinct breeding sites, and 53 cytochrome b haplotypes were identified. On the basis of comparison of the composition and frequency of haplotypes in each breeding site, salamanders in each habitat appeared to have distinct characteristics. Significant genetic differentiation was observed in 93.3% of possible pairs of habitats in Shiga prefecture, and 67.7% of habitat pairs within the same locality were found to be significantly different. These results suggest that the salamanders' poor locomotion combined with topographic effects may have contributed to the diversity of locally distributed salamanders in Shiga. Phylogenetic analysis showed that haplotypes of H. nebulosus in Shiga can be divided into five groups (the Nagahama-Maibara group, Otsu group, Konan group 1, Konan group 2, and Takashima group), each with a distinct geographical distribution. Haplotypes of the Otsu group, however, were exceptionally widely distributed. The results of the present study will contribute to the future of H. nebulosus conservation management in Shiga.

Successful re-establishment of a rabbit population from embryos vitrified 15 years ago: The importance of biobanks in livestock conservation

Genetic resource banks (GRB) are a valuable tool for maintaining genetic variability and preserving breeds from pathogens or catastrophe, enabling us to assess and correct breeding schemes, minimizing the impact of genetic drift and facilitating diffusion. This study tests their efficiency in re-establishing two extinct populations of a synthetic rabbit line selected for daily weight gain, using vitrified embryos from two generations (18th and 36th) separated by 15 years of genetic selection. The effect of long-term storage of vitrified embryos in liquid nitrogen was also evaluated. A total of 516 vitrified embryos using the same protocol were transferred into 54 recipients. The embryos had been maintained in liquid nitrogen during 2 different periods, (i) 1 year (301 embryos and 26 transfers, 36th generation) and (ii) 15 years (259 embryos and 28 transfers, 18th generation). A total of 80.0% (8/10 to 18th) and 60.0% (9/15 to 36th) of the foundational sire families were eventually re-established. Over approximately one year, animals within each population were crossed to produce the next generation and re-establish the original population size. Our study demonstrated that our GRB of embryos vitrified 15 years ago is a successful strategy to re-establish rabbit populations to continue the breeding programme.