Inbreeding in reintroduced populations: the effects of early reintroduction history and contemporary processes

Maintaining genetic variation and minimizing inbreeding are central goals of conservation genetics. It is therefore crucial to understand the important population parameters that affect inbreeding, particularly in reintroduction programs. Using data from 41 reintroduced Alpine ibex (Capra ibex ibex) populations we estimated inbreeding since the beginning of reintroductions using population-specific Fst, and inbreeding over the last few generations with contemporary effective population sizes. Total levels of inbreeding since reintroduction of ibex were, on average, close to that from one generation of half-sib mating. Contemporary effective population sizes did not reflect total inbreeding since reintroduction, but 16% of variation in contemporary effective population sizes among populations was due to variation in current population sizes. Substantial variation in inbreeding levels among populations was explained by founder group sizes and the harmonic mean population sizes since founding. This study emphasizes that, in addition to founder group sizes, early population growth rates are important parameters determining inbreeding levels in reintroduced populations.     

Group-structured genetic models in analyses of the population and behavioral ecology of poikilothermic vertebrates

Estimates of gene correlations among individuals within and among populations are frequently derived from statistical analyses of genetic data (e.g., F statistics). These measures can be important tools in molecular ecology and conservation, and offer important insights into population breeding structure. Using recently derived theory developed for group-structured populations, we show that fixation indices, when combined with basic population ecological and demographic data can be used to investigate population mating systems and to predict dispersal rates, trajectories and asymptotic levels of fixation indices, and effective population size. Four case studies of poikilothermic vertebrates are used to demonstrate the broad utility of evolutionary and ecological inferences afforded by group-structured models.     

Effective Size and F-Statistics of Subdivided Populations. I. Monoecious Species with Partial Selfing

Assuming discrete generations and autosomal inheritance involving genes that do not affect viability or reproductive ability, we have derived recurrence equations for the inbreeding coefficient and coancestry between individuals within and among subpopulations for a subdivided monoecious population with arbitrary distributions of male and female gametes per family, variable pollen and seed migration rates, and partial selfing. From the equations, formulas for effective size and expressions for F-statistics are obtained. For the special case of a single unsubdivided population, our equations reduce to the simple expressions derived by previous authors. It is shown that population structure (subdivision and migration) is important in determining the inbreeding coefficient and effective size. Failure to recognize internal structures of populations may lead to considerable bias in predicting effective size. Inbreeding coefficient, coancestry between individuals within and among subpopulations accrue at different and variable rates over initial generations before they converge to the same asymptotic rate of increase. For a given population, the smaller the pollen and seed migration rates, the more generations are required to attain the asymptotic rate and the larger the asymptotic effective size. The equations presented herein can be used for the study of evolutionary biology and conservation genetics.     

Effective Size and F-Statistics of Subdivided Populations. II. Dioecious Species

Assuming discrete generations and autosomal inheritance involving genes that do not affect viability or reproductive ability, we have derived recurrence equations for the inbreeding coefficient and coancestry between individuals within and among subpopulations for a subdivided monoecious population with arbitrary distributions of male and female gametes per family, variable pollen and seed migration rates, and partial selfing. From the equations, formulas for effective size and expressions for F-statistics are obtained. For the special case of a single unsubdivided population, our equations reduce to the simple expressions derived by previous authors. It is shown that population structure (subdivision and migration) is important in determining the inbreeding coefficient and effective size. Failure to recognize internal structures of populations may lead to considerable bias in predicting effective size. Inbreeding coefficient, coancestry between individuals within and among subpopulations accrue at different and variable rates over initial generations before they converge to the same asymptotic rate of increase. For a given population, the smaller the pollen and seed migration rates, the more generations are required to attain the asymptotic rate and the larger the asymptotic effective size. The equations presented herein can be used for the study of evolutionary biology and conservation genetics.     

Effective Sizes for Subdivided Populations

Many derivations of effective population sizes have been suggested in the literature; however, few account for the breeding structure and none can readily be expanded to subdivided populations. Breeding structures influence gene correlations through their effects on the number of breeding individuals of each sex, the mean number of progeny per female, and the variance in the number of progeny produced by males and females. Additionally, hierarchical structuring in a population is determined by the number of breeding groups and the migration rates of males and females among such groups. This study derives analytical solutions for effective sizes that can be applied to subdivided populations. Parameters that encapsulate breeding structure and subdivision are utilized to derive the traditional inbreeding and variance effective sizes. Also, it is shown that effective sizes can be determined for any hierarchical level of population structure for which gene correlations can accrue. Derivations of effective sizes for the accumulation of gene correlations within breeding groups (coancestral effective size) and among breeding groups (intergroup effective size) are given. The results converge to traditional, single population measures when similar assumptions are applied. In particular, inbreeding and intergroup effective sizes are shown to be special cases of the coancestral effective size, and intergroup and variance effective sizes will be equal if the population census remains constant. Instantaneous solutions for effective sizes, at any time after gene correlation begins to accrue, are given in terms of traditional F statistics or transition equations. All effective sizes are shown to converge upon a common asymptotic value when breeding tactics and migration rates are constant. The asymptotic effective size can be expressed in terms of the fixation indices and the number of breeding groups; however, the rate of approach to the asymptote is dependent upon dispersal rates. For accurate assessment of effective sizes, initial, instantaneous or asymptotic, the expressions must be applied at the lowest levels at which migration among breeding groups is nonrandom. Thus, the expressions may be applicable to lineages within socially structured populations, fragmented populations (if random exchange of genes prevails within each population), or combinations of intra- and interpopulation discontinuities of gene flow. Failure to recognize internal structures of populations may lead to considerable overestimates of inbreeding effective size, while usually underestimating variance effective size.     

Prevalence of multiple mating by female common dormice, <Emphasis Type="Italic">Muscardinus avellanarius</Emphasis>

Mating behaviour is an important component of species’ life histories. Knowledge of natural patterns of mating can lead also to more effective management strategies for populations of conservation concern. Despite a high conservation profile many aspects of the biology of the common dormouse (Muscardinus avellanarius) remain unknown, potentially limiting present conservation efforts. We determine the mating behaviour of M. avellanarius at two woodland sites in the UK: (1) Bontuchel (a natural population in Wales) and (2) Wych (a population in England that was established by reintroducing captive-bred animals) by genotyping mothers and litters at a panel of 10 microsatellite loci. Adult female body weight positively correlates with litter size and no apparent reproductive skew was evident. We found that multiple mating by female dormice is prevalent at both sites, with litters containing three or more offspring sired by multiple fathers; moreover, multiple mating is adopted by released animals even after a period of captive breeding where females are mated singly or as a breeding pair. We also present evidence for low proportion of fathers identified in our samples that probably related to unsampled individuals and/or larger than anticipated population sizes. This first report of mating behaviour in M. avellanarius highlights the role of genetic studies to uncover species’ reproductive behaviours and include these data for conservation management.     

Influence of Gene Flow and Breeding Tactics on Gene Diversity within Populations

Expressions describing the accumulation of gene correlations within and among lineages and individuals of a population are derived. The model permits different migration rates by males and females and accounts for various breeding tactics within lineages. The resultant equations enable calculation of the probabilistic quantities for the fixation indices, rates of loss of genetic variation, accumulation of inbreeding, and coefficients of relationship for the population at any generation. All fixation indices were found to attain asymptotic values rapidly despite the consistent loss of genetic variation and accumulation of inbreeding within the population. The time required to attain asymptotic values, however, was prolonged when gene flow among lineages was relatively low (less than 20%). The degree of genetic differentiation among breeding groups, inbreeding coefficients, and gene correlations within lineages were found to be primarily functions of breeding tactics within groups rather than gene flow among groups. Thus, the asymptotic value of S. Wright's island model is not appropriate for describing genetic differences among groups within populations. An alternative solution is provided that under limited conditions will reduce to the original island model. The evolution of polygynous breeding tactics appears to be more favorable for promoting intragroup gene correlations than modification of migration rates. Inbreeding and variance effective sizes are derived for populations that are structured by different migration and breeding tactics. Processes that reduce the inbreeding effective population size result in a concomitant increase in variance effective population size.     

Reconstructing recent divergence: evaluating nonequilibrium population structure in New Zealand chinook salmon

Newly established or perturbed populations are often the focus of conservation concerns but they pose special challenges for population genetics because drift?migration equilibrium is unlikely. To advance our understanding of the evolution of such populations, we investigated structure and gene flow among populations of chinook salmon that formed via natural straying following introduction to New Zealand in the early 1900s. We examined 11 microsatellite loci from samples collected in several sites and years to address two questions: (i) what population differentiation has arisen in the ≈ 30 generations since salmon were introduced to New Zealand, relative to temporal variation within populations; and (ii) what are the approximate effective population sizes and amounts of gene flow in these populations? These questions are routinely addressed in studies of indigenous populations, but less often in the case of new populations and rarely with consideration of equilibrium assumptions. We show that despite the recent introduction, continued gene flow and high temporal variability among samples, detectable population structure has arisen among the New Zealand populations, consistent with their colonization pattern and isolation by geographical distance. Furthermore, we use simple individual‐based simulations and estimates of effective population sizes to estimate the effective gene flow among drainages under likely nonequilibrium conditions. Similar methodology may be broadly applicable to other studies of population structure and phenotypic evolution under similar nonequilibrium, high gene flow conditions.     

Effective population size is strongly correlated with breeding pond size in the endangered California tiger salamander, <Emphasis Type="Italic">Ambystoma californiense</Emphasis>

Maintaining genetic diversity and population viability in endangered and threatened species is a primary concern of conservation biology. Genetic diversity depends on population connectivity and effective population size (N_e), both of which are often compromised in endangered taxa. While the importance of population connectivity and gene flow has been well studied, investigating effective population sizes in natural systems has received far less attention. However, N_e plays a prominent role in the maintenance of genetic diversity, the prevention of inbreeding depression, and in determining the probability of population persistence. In this study, we examined the relationship between breeding pond characteristics and N_e in the endangered California tiger salamander, Ambystoma californiense. We sampled 203 individuals from 10 breeding ponds on a local landscape, and used 11 polymorphic microsatellite loci to quantify genetic structure, gene flow, and effective population sizes. We also measured the areas of each pond using satellite imagery and classified ponds as either hydrologically-modified perennial ponds or naturally occurring vernal pools, the latter of which constitute the natural breeding habitat for A. californiense. We found no correlation between pond area and heterozygosity or allelic diversity, but we identified a strong positive relationship between breeding pond area and N_e, particularly for vernal pools. Our results provide some of the first empirical evidence that variation in breeding habitat can be associated with differences in N_e and suggest that a more complete understanding of the environmental features that influence N_e is an important component of conservation genetics and management.     

gesp: A computer program for modelling genetic effective population size,inbreeding and divergence in substructured populations

The genetically effective population size (N_e) is of key importance for quantifying rates of inbreeding and genetic drift and is often used in conservation management to set targets for genetic viability. The concept was developed for single, isolated populations and the mathematical means for analysing the expected N_e in complex, subdivided populations have previously not been available. We recently developed such analytical theory and central parts of that work have now been incorporated into a freely available software tool presented here. gesp (Genetic Effective population size, inbreeding and divergence in Substructured Populations) is R‐based and designed to model short‐ and long‐term patterns of genetic differentiation and effective population size of subdivided populations. The algorithms performed by gesp allow exact computation of global and local inbreeding and eigenvalue effective population size, predictions of genetic divergence among populations (G_ST) as well as departures from random mating (F_IS, F_IT) while varying (i) subpopulation census and effective size, separately or including trend of the global population size, (ii) rate and direction of migration between all pairs of subpopulations, (iii) degree of relatedness and divergence among subpopulations, (iv) ploidy (haploid or diploid) and (v) degree of selfing. Here, we describe gesp and exemplify its use in conservation genetics modelling.     

Relationships between migration rates and landscape resistance assessed using individual-based simulations

Linking landscape effects on gene flow to processes such as dispersal and mating is essential to provide a conceptual foundation for landscape genetics. It is particularly important to determine how classical population genetic models relate to recent individual-based landscape genetic models when assessing individual movement and its influence on population genetic structure. We used classical Wright-Fisher models and spatially explicit, individual-based, landscape genetic models to simulate gene flow via dispersal and mating in a series of landscapes representing two patches of habitat separated by a barrier. We developed a mathematical formula that predicts the relationship between barrier strength (i.e., permeability) and the migration rate (m) across the barrier, thereby linking spatially explicit landscape genetics to classical population genetics theory. We then assessed the reliability of the function by obtaining population genetics parameters (m, F(ST) ) using simulations for both spatially explicit and Wright-Fisher simulation models for a range of gene flow rates. Next, we show that relaxing some of the assumptions of the Wright-Fisher model can substantially change population substructure (i.e., F(ST) ). For example, isolation by distance among individuals on each side of a barrier maintains an F(ST) of ～0.20 regardless of migration rate across the barrier, whereas panmixia on each side of the barrier results in an F(ST) that changes with m as predicted by classical population genetics theory. We suggest that individual-based, spatially explicit modelling provides a general framework to investigate how interactions between movement and landscape resistance drive population genetic patterns and connectivity across complex landscapes.     

Genebank biodiversity assessments regarding optimal sample size and seed harvesting techniques for the regeneration of carrot accessions

Small, finite populations are particularly vulnerable to diversity loss during regeneration. The regeneration of a highly outbreeding open-pollinated variety relies on estimated effective population size, via the measurement of temporal change in allele frequencies. Using appropriate estimators for dominant gene markers, effective sizes were calculated for five sizes of a mating population and two seed harvesting procedures. We have shown that, in the case of carrot regeneration, 70 equally harvested plants should provide an effective size (N _e) of at least 50 plants. This value seems sufficient to limit genetic drift and to preserve an efficient level of genetic diversity within the collection. The efficiency of balanced samples (made of an equal number of seeds per plant) is compared to that of bulk samples (seeds randomly chosen among the total seed lot coming from all the plants).     

Microsatellite variation within and among populations of Oryza officinalis (Poaceae), an endangered wild rice from China

Oryza officinalis Wall. ex Watt. is an agriculturally important but seriously endangered species of wild rice. To obtain more accurate estimates of population structure for improved conservation planning of the species, genetic variability at 14 microsatellite DNA loci was examined in population samples covering most of the species' range in China. Considerable genetic variability (overall Na = 1.886, P = 62%, HO = 0.056, HE = 0.216, and HS = 0.277) was detected at the 14 loci in 442 individuals of the 18 natural populations. The evaluation of partitioning of genetic variability (FST = 0.442) suggested high genetic differentiation among the Chinese O. officinalis populations. An overall value of Nm = 0.316 suggested limited gene flow occurred among the sampled populations. Most of the populations showed heterozygote deficits in tests of Hardy-Weinberg equilibrium and significantly positive FIS values. This could be due to some inbreeding occurring in this predominantly outcrossing species. For effective in situ conservation and restoration genetics, maintenance of significant historical processes is particularly important, including high outbreeding, considerable gene flow, and large population effective sizes. The high FST values detected among populations in this study are instructive for adopting a conservation plan that includes representative populations with the greatest genetic variation for either in situ conservation management or germplasm collection expeditions.     

The mating system and gene dynamics of plateau pikas

Evolutionary theory suggests that mating systems should have substantial effects on gene dynamics of local populations. In polygynous species, local 'breeding groups' may produce significant genetic structure, due to genetic differences among groups, and rate of loss of genetic variation from such populations may be considerably slowed. We examined possible influences of the variable mating system and family group structure on genetic properties of a population of plateau pikas (Ochotona curzoniae). Pika gene dynamics were examined via F-statistics and effective population sizes (N(e)), calculated from genetic correlations within and among individuals and families. Genetic correlations were estimated from mating patterns, population demography, and dispersal patterns. Substantial genetic structure within the population was indicated by a strongly positive F(LS). Genetic influence of natal dispersal out of pika families was indicated by a strongly negative inbreeding statistic (F(IL)=-0.34). Effective size of the population was not greatly different from the census population, whereas a traditional estimate of effective size of the population was much lower, indicating that the family structure of the pikas results in a slowed loss of genetic variation over time. Thus, even though mating patterns of plateau pikas were variable, family structure had a strong influence on pika gene dynamics.     

Bayes estimation of species divergence times and ancestral population sizes using DNA sequences from multiple loci

The effective population sizes of ancestral as well as modern species are important parameters in models of population genetics and human evolution. The commonly used method for estimating ancestral population sizes, based on counting mismatches between the species tree and the inferred gene trees, is highly biased as it ignores uncertainties in gene tree reconstruction. In this article, we develop a Bayes method for simultaneous estimation of the species divergence times and current and ancestral population sizes. The method uses DNA sequence data from multiple loci and extracts information about conflicts among gene tree topologies and coalescent times to estimate ancestral population sizes. The topology of the species tree is assumed known. A Markov chain Monte Carlo algorithm is implemented to integrate over uncertain gene trees and branch lengths (or coalescence times) at each locus as well as species divergence times. The method can handle any species tree and allows different numbers of sequences at different loci. We apply the method to published noncoding DNA sequences from the human and the great apes. There are strong correlations between posterior estimates of speciation times and ancestral population sizes. With the use of an informative prior for the human-chimpanzee divergence date, the population size of the common ancestor of the two species is estimated to be approximately 20,000, with a 95% credibility interval (8000, 40,000). Our estimates, however, are affected by model assumptions as well as data quality. We suggest that reliable estimates have yet to await more data and more realistic models.     

Efficient genetic markers for population biology

Population genetics has come of age. Three important components have come together: efficient techniques to examine informative segments of DNA, statistics to analyse DNA data and the availability of easy-to-use computer packages. Single-locus genetic markers and those that produce gene genealogies yield information that is truly comparable among studies. These markers answer biological questions most efficiently and also contribute to much broader investigations of evolutionary, population and conservation biology. For these reasons, single-locus and genealogical markers should be the focus of the intensive genetic data collection that has begun owing to the power of genetics in population biology.     

普通野生稻小种群的交配系统与遗传多样性

小种群的遗传动态是保育遗传学关注的核心问题之一,而种群遗传动态又与交配系统密切相关.普通野生稻(Oryza rufipogon Griff.)是具有重要经济价值的濒危物种,目前其种群规模都较小,研究其小种群交配系统与遗传变异性对普通野生稻的保护具有重要意义.运用7对SSR引物,对采自江西东乡普通野生稻小种群的36份种茎和其中20个家系共计601份子代进行了分析.结果显示:该种群的表观异交率为0.318,多位点法估计(MLTR)的多位点异交率为0.481;50%以上的子代共享亲本,非随机交配明显;东乡普通野生稻种群交配系统属于混合交配类型.比较亲本和子代种群的遗传变异性显示:子代种群比亲本种群遗传变异性更丰富;子代种群的杂合子不足与种群变小自交比例上升有关;而亲本种群杂合子过剩可能与杂合基因型的选择优势有关.这些结果说明创造条件扩大种群规模对普通野生稻的原生境保护显得尤为重要.     

Population genetics meets behavioral ecology

Populations are often composed of more than just randomly mating subpopulations - many organisms from social groups with distinct patterns of mating and dispersal. Such patterns have recieved much attention in behavioral ecology, yet theories of population genetics rarely take social structures into account. Consequently, population geneticists often report high levels of apparent in breeding and concomitantly low efective sizes, even for species that avoid mating between close kin. Recently, a view of gene dynamics has been introduced that takes dispersal and social structure into account. Accounting for social structure in population genetics leads to a different perspective on how genetic variation is partitoned and the rate at which genic diversity is lost in natural populations - a view that is more consistent with observed behaviors for the minimization of inbreeding.     

Impact of mating design on selection response in Brassica rapa L.

The impact of four mating designs on selection response for leaf area was assessed at four different population sizes, using fast-cycling Brassica rapa L. Mating designs were either balanced (partial diallel or pair mating) or unbalanced (factorial mating designs with either one or two testers). When balanced, the mating designs required different numbers of crossings for the same number of parents: the partial diallel design, in the configuration retained here, required three times as many crossings as pair mating. Population sizes were 4, 8, 16, and 32. The percentage of selected individuals was kept constant at 25%. Despite an average estimated heritability around 0.4, the overall response to selection after five generations was fairly weak in all three replicates. For a given population size, selection response was larger under balanced mating designs than under unbalanced ones. There was no difference among balanced mating designs. Both results indicate that effective population size is more important than population size or the number of crossings in maintaining genetic gain.     

群体遗传结构中的基因流

群体遗传结构上的差异是遗传多样性的一种重要体现,对群体遗传结构的研究已有较久的历史,而其中的基因流研究近些年来越来越受到重视。它对群体遗传学、进化生物学、保护生物学、生态学有着极其重要的作用。虽然传统的群体遗传学能估测基因流大小,但它的精确性还有很大局限性。随着生物技术的进步,对基因流的研究逐渐向分子水平过渡,应用蛋白质电泳技术、分子标记技术(RAPD、RFLP、VNTR、ISSR、DNA测序等)方法对群体间基因流的流动水平进行了深入细致的研究。通过综述群体遗传结构的几种模式:陆岛模式、海岛模式、阶石模式、距离隔离模式、层次模式,以及在群体遗传结构的几种模式基础上的基因流的研究方法、作用、地位和近些年来研究者的研究成果,并指出了这些方法的局限性。