Integrating Genetic Data and Demographic Modeling to Facilitate Conservation of Small,Isolated Mountain Goat Populations

Acquisition of field data and analytical methods needed for conservation and management of wildlife populations represent significant challenges, particularly for species that inhabit landscapes that are difficult to access or species that persist in small, isolated populations. In such instances, integrating diverse and complementary data streams, such as genetic and non-genetic data, can advance our understanding of population dynamics and associated management implications. We examined how genetic and morphologic data can be used to articulate population structure of a low-density, peninsular population of mountain goats (Oreamnos americanus) on the Cleveland Peninsula, Alaska, USA, and surrounding areas, 2005–2018. We then use a population demographic modeling approach to examine how the use of population structure information influences sustainable harvest quotas, as compared to a panmictic, null model. Specifically, we conducted extensive field sampling of genetic (n = 446) and morphologic (i.e., horn length, n = 371) data to characterize population structure. We conducted demographic analyses and examined harvest modeling scenarios using a sex- and age-specific matrix population modeling approach. Genetic and morphologic data analyses suggested peninsular subpopulations were demographically isolated, relative to surrounding mainland populations. Specifically, genetic structuring was evident and followed an isolation-by-distance, stepping-stone pattern indicating limited interchange, low effective population sizes, and reduced genetic diversity along a peninsular extremity to mainland gradient. Harvest modeling indicated that overharvest would likely occur if the panmictic, null model was used to guide harvest because the smallest genetically defined population at the peninsular extremity was too small to permit any level of sustainable harvest. Our analyses illustrate the importance of using genetic and morphologic data, in combination with demographic modeling, to quantitatively delineate population boundaries and dynamics for ensuring viability of small, isolated populations. © 2020 The Wildlife Society.     

TEMPORAL FLUCTUATIONS IN DEMOGRAPHIC PARAMETERS AND THE GENETIC VARIANCE AMONG POPULATIONS

Contrary to assumptions commonly made in the study of population genetics, the demographic properties of many populations are not always constant. Important characteristics of populations such as migration rate and population size may vary in time and space. Moreover, local populations often come and go; the rate of extinction and the properties of colonization may also vary. In this paper, the approach to equilibrium following a disturbance in the genetic variance among populations is described. The rate of migration is shown to be critical in determining the extent to which extinction and recolonization affects genetic differentiation. Perturbations and variations through time and space in demographic parameters such as population size and migration rate are shown to be important in determining the partitioning of genetic variance. Equations are given to predict the average through time of genetic differentiation among populations in the event of a single disturbance or in constant fluctuations in the pertinent demographic parameters. In general, these fluctuations increase the F_ST of a species. Spatial demographic variation affects F_STmuch more than temporal variation. These demographic properties make some species unsuitable for the empirical analysis of migration with indirect genetic measures. Demographic instability may play a large role in the evolution of genetic variation.     

Microsatellite DNA analysis of northern pike (Esox lucius L.) populations: insights into the genetic structure and demographic history of a genetically depauperate species

The northern pike Esox lucius L. is a freshwater fish exhibiting pronounced population subdivision and low genetic variability. However, there is limited knowledge on phylogeographical patterns within the species, and it is not known whether the low genetic variability reflects primarily current low effective population sizes or historical bottlenecks. We analysed six microsatellite loci in ten populations from Europe and North America. Genetic variation was low, with the average number of alleles within populations ranging from 2.3 to 4.0 per locus. Genetic differentiation among populations was high (overall θ_ST = 0.51; overall ρ_ST = 0.50). Multidimensional scaling analysis of genetic distances between populations and spatial analysis of molecular variance suggested a single phylogeographical race within the sampled populations from northern Europe, whereas North American and southern European populations were highly distinct. A population from Ireland was monomorphic at all loci, presumably reflecting founder events associated with introduction of the species to the island in the sixteenth century. Bayesian analysis of demographic parameters showed differences in θ (a product of effective population size and mutation rate) among populations from large and small water bodies, but the relative differences in θ were smaller than expected, which could reflect population subdivision within the larger water bodies. Finally, the analyses showed drastic population declines on a time scale of several thousand years within European populations, which we ascribe to either glacial bottlenecks or postglacial founder events.  © 2005 The Linnean Society of London, Biological Journal of the Linnean Society , 2005, 84 , 91–101.     

Possible Secondary Population-Level Effects of Selective Harvest of Adult Male Muskoxen

Selective harvest regimes are often focused on males resulting in skewed sex-ratios, and for many ungulate species this strategy is sustainable. However, muskoxen (Ovibos moschatus) are very social and mature bulls (≥4 years old), particularly prime-age bulls (6–10 years old), play important roles in predator defense and recruitment. A year-round social structure incorporating large males into mixed-sex groups could make this species more susceptible to the effects of selective harvest if population composition and sex-ratios influence overall survival and reproductive success. Using detailed data collected on the muskox population occupying the Seward Peninsula, Alaska during 2002–2012, we formulated the hypothesis that the selective harvest of mature bulls may be related to documented changes in population composition and growth rates in this species. In addition, we reviewed existing published information from two other populations in Alaska, the Cape Thompson and Northeastern populations, to compare population growth rates among the three areas under differential harvest rates relative to our hypothesis. We found that on the Seward Peninsula, mature bull:adult cow ratios declined 4–12%/year and short-yearling:adult cow ratios (i.e., recruitment) declined 8–9%/year in the most heavily harvested areas. Growth rates in all 3 populations decreased disproportionately after increases in the number of bulls harvested, and calf:cow ratios declined in the Northeastern population as harvest increased. While lack of appropriate data prevented us from excluding other potential causes such as density dependent effects and changes in predator densities, our results did align with our hypothesis, suggesting that in the interest of conservation, harvest of mature males should be restricted until causal factors can be more definitively identified. If confirmed by additional research, our findings would have important implications for harvest management and conservation of muskoxen and other ungulate species with similar life-histories.     

Demography of the Pacific walrus (Odobenus rosmarus divergens): 1974–2006

Global climate change may fundamentally alter population dynamics of many species for which baseline population parameter estimates are imprecise or lacking. Historically, the Pacific walrus is thought to have been limited by harvest, but it may become limited by global warming‐induced reductions in sea ice. Loss of sea ice, on which walruses rest between foraging bouts, may reduce access to food, thus lowering vital rates. Rigorous walrus survival rate estimates do not exist, and other population parameter estimates are out of date or have well‐documented bias and imprecision. To provide useful population parameter estimates we developed a Bayesian, hidden process demographic model of walrus population dynamics from 1974 through 2006 that combined annual age‐specific harvest estimates with five population size estimates, six standing age structure estimates, and two reproductive rate estimates. Median density independent natural survival was high for juveniles (0.97) and adults (0.99), and annual density dependent vital rates rose from 0.06 to 0.11 for reproduction, 0.31 to 0.59 for survival of neonatal calves, and 0.39 to 0.85 for survival of older calves, concomitant with a population decline. This integrated population model provides a baseline for estimating changing population dynamics resulting from changing harvests or sea ice.     

Molecular ecological approaches to studying the evolutionary impact of selective harvesting in wildlife

Harvesting of wildlife populations by humans is usually targeted by sex, age or phenotypic criteria, and is therefore selective. Selective harvesting has the potential to elicit a genetic response from the target populations in several ways. First, selective harvesting may affect population demographic structure (age structure, sex ratio), which in turn may have consequences for effective population size and hence genetic diversity. Second, wildlife-harvesting regimes that use selective criteria based on phenotypic characteristics (e.g. minimum body size, horn length or antler size) have the potential to impose artificial selection on harvested populations. If there is heritable genetic variation for the target characteristic and harvesting occurs before the age of maturity, then an evolutionary response over time may ensue. Molecular ecological techniques offer ways to predict and detect genetic change in harvested populations, and therefore have great utility for effective wildlife management. Molecular markers can be used to assess the genetic structure of wildlife populations, and thereby assist in the prediction of genetic impacts by delineating evolutionarily meaningful management units. Genetic markers can be used for monitoring genetic diversity and changes in effective population size and breeding systems. Tracking evolutionary change at the phenotypic level in the wild through quantitative genetic analysis can be made possible by genetically determined pedigrees. Finally, advances in genome sequencing and bioinformatics offer the opportunity to study the molecular basis of phenotypic variation through trait mapping and candidate gene approaches. With this understanding, it could be possible to monitor the selective impacts of harvesting at a molecular level in the future. Effective wildlife management practice needs to consider more than the direct impact of harvesting on population dynamics. Programs that utilize molecular genetic tools will be better positioned to assess the long-term evolutionary impact of artificial selection on the evolutionary trajectory and viability of harvested populations.     

Microsatellite analysis of a population crash and bottleneck in the Mauna Kea silversword, Argyroxiphium sandwicense ssp. sandwicense (Asteraceae), and its implications for reintroduction

The Mauna Kea silversword, Argyroxiphium sandwicense ssp. sandwicense, has experienced both a severe population crash associated with an increase in alien ungulate populations on Mauna Kea, and a population bottleneck associated with reintroduction. In this paper, we address the genetic consequences of both demographic events using eight microsatellite loci. The population crash was not accompanied by a significant reduction in number of alleles or heterozygosity. However, the population bottleneck was accompanied by significant reductions in observed number of alleles, effective number of alleles, and expected heterozygosity, though not in observed heterozygosity. The effective size of the population bottleneck was calculated using both observed heterozygosities and allele frequency variances. Both methods corroborated the historical census size of the population bottleneck of at most three individuals. The results suggest that: (i) small populations, even those that result from severe reductions in historical population size and extent, are not necessarily genetically depauperate; and (ii) species reintroduction plans need to be conceived and implemented carefully, with due consideration to the genetic impact of sampling for reintroduction.     

Demographic Model of the Swiss Cattle Population for the Years 2009-2011 Stratified by Gender,Age and Production Type

Demographic composition and dynamics of animal and human populations are important determinants for the transmission dynamics of infectious disease and for the effect of infectious disease or environmental disasters on productivity. In many circumstances, demographic data are not available or of poor quality. Since 1999 Switzerland has been recording cattle movements, births, deaths and slaughter in an animal movement database (AMD). The data present in the AMD offers the opportunity for analysing and understanding the dynamic of the Swiss cattle population. A dynamic population model can serve as a building block for future disease transmission models and help policy makers in developing strategies regarding animal health, animal welfare, livestock management and productivity. The Swiss cattle population was therefore modelled using a system of ordinary differential equations. The model was stratified by production type (dairy or beef), age and gender (male and female calves: 0–1 year, heifers and young bulls: 1–2 years, cows and bulls: older than 2 years). The simulation of the Swiss cattle population reflects the observed pattern accurately. Parameters were optimized on the basis of the goodness-of-fit (using the Powell algorithm). The fitted rates were compared with calculated rates from the AMD and differed only marginally. This gives confidence in the fitted rates of parameters that are not directly deductible from the AMD (e.g. the proportion of calves that are moved from the dairy system to fattening plants).     

Patterns and causes of demographic variation in a harvested moose population: evidence for the effects of climate and density-dependent drivers

1.?Better understanding of the mechanisms affecting demographic variation in ungulate populations is needed to support sustainable management of harvested populations. While studies of moose Alces alces L. populations have previously explored temporal variation in demographic processes, managers responsible for populations that span large heterogeneous landscapes would benefit from an understanding of how demography varies across biogeographical gradients in climate and other population drivers. Evidence of thresholds in population response to manageable and un-manageable drivers could aid resource managers in identifying limits to the magnitude of sustainable change. 2.?Generalized additive models (GAMs) were used to evaluate the relative importance of population density, habitat abundance, summer and winter climatic conditions, primary production, and harvest intensity in explaining spatial variation in moose vital rates in Ontario, Canada. Tree regression was used to test for thresholds in the magnitudes of environmental predictor variables that significantly affected population vital rates. 3.?Moose population growth rate was negatively related to moose density and positively related to the abundance of mixed deciduous habitat abundant in forage. Calf recruitment was negatively related to a later start of the growing season and calf harvest. The ratio of bulls to cows was related to male harvest and hunter access, and thresholds were evident in predictor variables for all vital rate models. 4.?Findings indicate that the contributions of density-dependent and independent factors can vary depending on the scale of population process. The importance of density dependence and habitat supply to low-density ungulate populations was evident, and management strategies for ungulates may be improved by explicitly linking forest management and harvest. Findings emphasize the importance of considering summer climatic influences to ungulate populations, as recruitment in moose was more sensitive to the timing of vegetation green-up than winter severity. The efficacy of management decisions for harvested ungulates may require regional shifts in targets where populations span bioclimatic gradients. The use of GAMs in combination with recursive partitioning was demonstrated to be an informative analytical framework that captured nonlinear relationships common in natural processes and thresholds that are relevant to population management in diverse systems.     

Population size,habitat fragmentation,and the nature of adaptive variation in a stream fish

Whether and how habitat fragmentation and population size jointly affect adaptive genetic variation and adaptive population differentiation are largely unexplored. Owing to pronounced genetic drift, small, fragmented populations are thought to exhibit reduced adaptive genetic variation relative to large populations. Yet fragmentation is known to increase variability within and among habitats as population size decreases. Such variability might instead favour the maintenance of adaptive polymorphisms and/or generate more variability in adaptive differentiation at smaller population size. We investigated these alternative hypotheses by analysing coding-gene, single-nucleotide polymorphisms associated with different biological functions in fragmented brook trout populations of variable sizes. Putative adaptive differentiation was greater between small and large populations or among small populations than among large populations. These trends were stronger for genetic population size measures than demographic ones and were present despite pronounced drift in small populations. Our results suggest that fragmentation affects natural selection and that the changes elicited in the adaptive genetic composition and differentiation of fragmented populations vary with population size. By generating more variable evolutionary responses, the alteration of selective pressures during habitat fragmentation may affect future population persistence independently of, and perhaps long before, the effects of demographic and genetic stochasticity are manifest.     

The effects of cyclic dynamics and mating system on the effective size of an island mouflon population

The Haute Island mouflon (Ovis aries) population is isolated on one small (6.5 km2) island of the remote Kerguelen archipelago. Given a promiscuous mating system, a cyclic demography and a strong female-biased sex ratio after population crashes, we expected a low effective population size (Ne). We estimated Ne using demographic and temporal genetic approaches based on genetic information at 25 microsatellite loci from 62 and 58 mouflons sampled in 1988 and 2003, respectively. Genetic Ne estimates were higher than expected, varying between 104 and 250 depending on the methods used. Both demographic and genetic approaches show the Haute Island Ne is buffered against population crashes. The unexpectedly high Ne likely results from the cyclic winter crashes that allow young males to reproduce, limiting the variance of male reproductive success. Based on individual-based simulations, we suggest that despite a strongly female-biased sex ratio, the effects of the mating system on the effective population size more closely resemble random mating or weak polygyny.     

Disentangling the impact of demographic factors on population differentiation of an endangered freshwater crayfish (Austropotamobius pallipes) using population density and microsatellite data

1. Habitat fragmentation of stream ecosystems often results in decreased connectivity between populations and lower population sizes. Hence, understanding how habitat fragmentation affects genetic erosion is important for the preservation of freshwater biodiversity, in particular, as small populations suffer from loss of genetic diversity through genetic drift and loss of fitness because of inbreeding, increasing the risk of extinction. 2. Here, we assess the impact of demographic factors on population differentiation in the endangered freshwater crayfish Austropotamobius pallipes by analysing population genetic structure, estimating effective population sizes and comparing levels of polymorphism at five microsatellite loci with density estimates of 10 populations within a small French catchment that has become progressively confined to headwaters over the last six decades. 3. Levels of expected heterozygosity and allelic richness per population were relatively low (0.214–0.396 and 1.6–2.6, respectively). We found strong genetic differentiation between these geographically close populations (F_ST = 0.283), with weak statistical evidence for a pattern of isolation by distance. Estimates of effective population size were low (<150) in most populations, but potentially reached several thousands in three populations. 4. Population density and allelic richness were strongly positively correlated. A robust relationship between population density and heterozygosity values was also noted, but only after discarding two populations for which significant genetic signatures of a recent bottleneck were found; these two populations displayed high expected heterozygosity compared with a very low density. Populations with the highest densities of individuals had the highest effective population size estimates and vice versa. 5. Our results clearly show the importance of demographic factors and genetic drift on A. pallipes populations. Furthermore, analysis of genetic and population density data is a pragmatic and efficient approach to corroborate inferences from genetic data and can be particularly useful in the identification of populations experiencing a bottleneck and therefore in conservation genetics studies aiming at identifying priority populations for conservation.     

Estimating the growth of a newly established moose population using reproductive value

Estimating the population growth rate and environmental stochasticity of long-lived species is difficult because annual variation in population size is influenced by temporal autocorrelations caused by fluctuations in the age-structure. Here we use the dynamics of the reproductive value to estimate the long-term growth rate s and the environmental variance of a moose population that recently colonized the island of Vega in northern Norway. We show that the population growth rate was high (ŝ=0.26). The major stochastic influences on the population dynamics were due to demographic stochasticity, whereas the environmental variance was not significantly different from 0. This supports the suggestion that population growth rates of polytocous ungulates are high, and that demographic stochasticity must be assessed when estimating the growth of small ungulate populations.     

Genetic effective size of a wild primate population: influence of current and historical demography

A comprehensive assessment of the determinants of effective population size (N(e)) requires estimates of variance in lifetime reproductive success and past changes in census numbers. For natural populations, such information can be best obtained by combining longitudinal data on individual life histories and genetic marker-based inferences of demographic history. Independent estimates of the variance effective size (N(ev), obtained from life-history data) and the inbreeding effective size (N((eI), obtained from genetic data) provide a means of disentangling the effects of current and historical demography. The purpose of this study was to assess the demographic determinants of N(e) in one of the most intensively studied natural populations of a vertebrate species: the population of savannah baboons (Papio cynocephalus) in the Amboseli Basin, southern Kenya. We tested the hypotheses that N(eV) < N < N(eI) (where N = population census number) due to a recent demographic bottleneck. N(eV) was estimated using a stochastic demographic model based on detailed life-history data spanning a 28-year period. Using empirical estimates of age-specific rates of survival and fertility for both sexes, individual-based simulations were used to estimate the variance in lifetime reproductive success. The resultant values translated into an N(eV)/N estimate of 0.329 (SD = 0.116, 95% CI = 0.172-0.537). Historical N(eI), was estimated from 14-locus microsatellite genotypes using a coalescent-based simulation model. Estimates of N(eI) were 2.2 to 7.2 times higher than the contemporary census number of the Amboseli baboon population. In addition to the effects of immigration, the disparity between historical N(eI) and contemporary N is likely attributable to the time lag between the recent drop in census numbers and the rate of increase in the average probability of allelic identity-by-descent. Thus, observed levels of genetic diversity may primarily reflect the population's prebottleneck history rather than its current demography.     

The diversity of population responses to environmental change

The current extinction and climate change crises pressure us to predict population dynamics with ever‐greater accuracy. Although predictions rest on the well‐advanced theory of age‐structured populations, two key issues remain poorly explored. Specifically, how the age‐dependency in demographic rates and the year‐to‐year interactions between survival and fecundity affect stochastic population growth rates. We use inference, simulations and mathematical derivations to explore how environmental perturbations determine population growth rates for populations with different age‐specific demographic rates and when ages are reduced to stages. We find that stage‐ vs. age‐based models can produce markedly divergent stochastic population growth rates. The differences are most pronounced when there are survival‐fecundity‐trade‐offs, which reduce the variance in the population growth rate. Finally, the expected value and variance of the stochastic growth rates of populations with different age‐specific demographic rates can diverge to the extent that, while some populations may thrive, others will inevitably go extinct.     

Effective size of a fluctuating age-structured population

Previous theories on the effective size of age-structured populations assumed a constant environment and, usually, a constant population size and age structure. We derive formulas for the variance effective size of populations subject to fluctuations in age structure and total population size produced by a combination of demographic and environmental stochasticity. Haploid and monoecious or dioecious diploid populations are analyzed. Recent results from stochastic demography are employed to derive a two-dimensional diffusion approximation for the joint dynamics of the total population size, N, and the frequency of a selectively neutral allele, p. The infinitesimal variance for p, multiplied by the generation time, yields an expression for the effective population size per generation. This depends on the current value of N, the generation time, demographic stochasticity, and genetic stochasticity due to Mendelian segregation, but is independent of environmental stochasticity. A formula for the effective population size over longer time intervals incorporates deterministic growth and environmental stochasticity to account for changes in N.     

Two measures of effective population size for graphs

Effective population size is a key parameter in population ecology because it allows prediction of the dynamics of genetic variation and the rate of genetic drift and inbreeding. It is important for the definition of "nearly neutral" mutations and, hence, has consequences for the fixation or extinction probabilities of advantageous and deleterious mutations. As graph-based population models become increasingly popular for studying evolution in spatially or socially structured populations, a neutral theory for evolution on graphs is called for. Here, we derive formulae for two alternative measures of effective population size, the variance effective and inbreeding effective size of general unweighted and undirected graphs. We show how these two quantities relate to each other and we derive effective sizes for the complete graph the cycle and bipartite graphs. For one-dimensional lattices and small-world graphs, we estimate the inbreeding effective size using simulations. The presented method is suitable for any structured population of haploid individuals with overlapping generations.     

Scaling of the mean and variance of population dynamics under fluctuating regimes

Theoretical ecologists have long sought to understand how the persistence of populations depends on the interactions between exogenous (biotic and abiotic) and endogenous (e.g., demographic and genetic) drivers of population dynamics. Recent work focuses on the autocorrelation structure of environmental perturbations and its effects on the persistence of populations. Accurate estimation of extinction times and especially determination of the mechanisms affecting extinction times is important for biodiversity conservation. Here we examine the interaction between environmental fluctuations and the scaling effect of the mean population size with its variance. We investigate how interactions between environmental and demographic stochasticity can affect the mean time to extinction, change optimal patch size dynamics, and how it can alter the often-assumed linear relationship between the census size and the effective population size. The importance of the correlation between environmental and demographic variation depends on the relative importance of the two types of variation. We found the correlation to be important when the two types of variation were approximately equal; however, the importance of the correlation diminishes as one source of variation dominates. The implications of these findings are discussed from a conservation and eco-evolutionary point of view.     

Reproductive success and effective population size in woodrats (Neotoma macrotis)

Discrepancies between the census size and the genetically effective size of populations (N(e)) can be caused by a number of behavioural and demographic factors operating within populations. Specifically, strong skew in male reproductive success, as would be expected in a polygynous mating system, could cause a substantial decrease in N(e) relative to census size. Because the mating system of Neotoma macrotis had previously been described as one nearing harem polygyny, I examined the distribution of reproductive success and genetic variation within a population of this species. Combining genetic data and three years of field observations, I show that variance in reproductive success does not deviate from poisson expectations within either sex and variance in success is similar between the sexes. Furthermore, both males and females had multiple partners across litters in addition to some evidence of multiple paternity within litters. Despite a lack of strong skew in reproductive success, an estimate of N(e) based on a number of demographic parameters suggests that the ratio of N(e)/N in this population is 0.48. Although the ratio of N(e)/N suggests that the population is experiencing higher rates of genetic drift than would be expected based on census size alone, the population maintains high levels of genetic diversity. Estimates of neighbourhood size and patterns of recruitment to the study site suggest that immigration plays an important role in this population and may contribute to the maintenance of high levels of genetic diversity.     

Do hypotheses from short‐term studies hold in the long‐term? An empirical test

Abstract. 1. A sequence of population estimates for two now-extinct populations of Euphydryas editha bayensis is presented. After removing biased sampling days, estimates of demographic parameters from the long-term data were used to test five hypotheses built from studies of shorter duration. Such tests of short-term conclusions are rare.
2. The long-term demographic parameters include sex ratio, mortality, dispersal, mean flight date, and duration of flight season. The two populations differed with respect to sex ratio and mean flight date, and sexes differed with respect to mortality and dispersal.
3. Three of the five hypotheses were supported directly or indirectly by patterns in the parameters. These hypotheses predict that dynamics are asynchronous over the long term, that larval mortality, not adult abundance and mortality, is the primary determinant of changes in population size, and that topography mediates larval mortality.
4. Two hypotheses were not supported or supported only in part. Flight phenology differed between the study populations as predicted, but flight order was opposite that expected from the topographic composition of each habitat. Variability in sex ratio and the occurrence of female-biased ratios in the habitat of one of the populations also suggest that previous observations of sex ratio are not generalisable.
5. Populations were extremely volatile over the study period. Removal of biased sampling days did not change basic trends or fluctuations in the data. This volatility suggests that E. editha populations residing in similar habitats may risk immediate extinction.