Increased inbreeding and strong kinship structure in Taxus baccata estimated from both AFLP and SSR data

Habitat fragmentation can have severe genetic consequences for trees, such as increased inbreeding and decreased effective population size. In effect, local populations suffer from reduction of genetic variation, and thus loss of adaptive capacity, which consequently increases their risk of extinction. In Europe, Taxus baccata is among a number of tree species experiencing strong habitat fragmentation. However, there is little empirical data on the population genetic consequences of fragmentation for this species. This study aimed to characterize local genetic structure in two natural remnants of English yew in Poland based on both amplified fragment length polymorphism (AFLP) and microsatellite (SSR) markers. We introduced a Bayesian approach that estimates the average inbreeding coefficient using AFLP (dominant) markers. Results showed that, in spite of high dispersal potential (bird-mediated seed dispersal and wind-mediated pollen dispersal), English yew populations show strong kinship structure, with a spatial extent of 50–100 m, depending on the population. The estimated inbreeding levels ranged from 0.016 to 0.063, depending on the population and marker used. Several patterns were evident: (1) AFLP markers showed stronger kinship structure than SSRs; (2) AFLP markers provided higher inbreeding estimates than SSRs; and (3) kinship structure and inbreeding were more pronounced in denser populations regardless of the marker used. Our results suggest that, because both kinship structure and (bi-parental) inbreeding exist in populations of English yew, gene dispersal can be fairly limited in this species. Furthermore, at a local scale, gene dispersal intensity can be more limited in a dense population.     

Evaluating sample size to estimate genetic management metrics in the genomics era

Inbreeding and relationship metrics among and within populations are useful measures for genetic management of wild populations, but accuracy and precision of estimates can be influenced by the number of individual genotypes analysed. Biologists are confronted with varied advice regarding the sample size necessary for reliable estimates when using genomic tools. We developed a simulation framework to identify the optimal sample size for three widely used metrics to enable quantification of expected variance and relative bias of estimates and a comparison of results among populations. We applied this approach to analyse empirical genomic data for 30 individuals from each of four different free‐ranging Rocky Mountain bighorn sheep (Ovis canadensis canadensis) populations in Montana and Wyoming, USA, through cross‐species application of an Ovine array and analysis of approximately 14,000 single nucleotide polymorphisms (SNPs) after filtering. We examined intra‐ and interpopulation relationships using kinship and identity by state metrics, as well as F_ST between populations. By evaluating our simulation results, we concluded that a sample size of 25 was adequate for assessing these metrics using the Ovine array to genotype Rocky Mountain bighorn sheep herds. However, we conclude that a universal sample size rule may not be able to sufficiently address the complexities that impact genomic kinship and inbreeding estimates. Thus, we recommend that a pilot study and sample size simulation using R code we developed that includes empirical genotypes from a subset of populations of interest would be an effective approach to ensure rigour in estimating genomic kinship and population differentiation.     

Quantifying gene flow from spatial genetic structure data in a metapopulation of Chamaecrista fasciculata (Leguminosae)

Abstract An extensive allozyme survey was conducted within a natural "meta" population of the native North American annual legume, Chamaecrista fasciculata (Leguminosae) to quantify genetic structure at different spatial scales. Gene flow was then estimated by a recently developed indirect method based on a continuous population model, using pairwise kinship coefficients between individuals. The indirect estimates of gene flow, quantified in terms of neighborhood size, with an average value on the order of 150 individuals, were concordant among different spatial scales (subpopulation, population, metapopulation). This gene-flow value lies within the range of direct estimates previously documented from observations of pollen and seed dispersal for the same metapopulation. Monte Carlo simulations using the direct measures of gene flow as parameters further demonstrated that the observed spatial pattern of allozyme variation was congruent with a model of isolation by distance. Combining previously published estimates of pollen dispersal distances with kinship coefficients from this study, we quantified biparental inbreeding relative to either a single subpopulation or the whole metapopulation. At the level of a neighborhood, little biparental inbreeding was observed and most departure from Hardy-Weinberg genotypic proportions was explained by self-fertilization, whereas both selfing and biparental inbreeding contributed to nonrandom mating at the metapopulation level. Gene flow was also estimated from indirect methods based on a discontinuous population structure model. We discuss these results with respect to the effect of a patchy population structure on estimation of gene flow.     

On the expected relationship between inbreeding,fitness, and extinction

We assessed the expected relationship between the level and the cost of inbreeding, measured either in terms of fitness, inbreeding depression or probability of extinction. First, we show that the assumption of frequent, slightly deleterious mutations do agree with observations and experiments, on the contrary to the assumption of few, moderately deleterious mutations. For the same inbreeding coefficient, populations can greatly differ in fitness according to the following: (i) population size; larger populations show higher fitness (ii) the history of population size; in a population that recovers after a bottleneck, higher inbreeding can lead to higher fitness and (iii) population demography; population growth rate and carrying capacity determine the relationship between inbreeding and extinction. With regards to the relationship between inbreeding depression and inbreeding coefficient, the population size that minimizes inbreeding depression depends on the level of inbreeding: inbreeding depression can even decrease when population size increases. It is therefore clear that to infer the costs of inbreeding, one must know both the history of inbreeding (e.g. past bottlenecks) and population demography.     

A metapopulation perspective on genetic diversity and differentiation in partially self-fertilizing plants

Abstract.— Partial self-fertilization is common in higher plants. Mating system variation is known to have important consequences for how genetic variation is distributed within and among populations. Selfing is known to reduce effective population size, and inbreeding species are therefore expected to have lower levels of genetic variation than comparable out crossing taxa. However, several recent empirical studies have shown that reductions in genetic diversity within populations of inbreeding species are far greater than the expected reductions based on the reduced effective population size. Two different processes have been argued to cause these patterns, selective sweeps (or hitchhiking) and background selection. Both are expected to be most effective in reducing genetic variation in regions of low recombination rates. Selfing is known to reduce the effective recombination rate, and inbreeding taxa are thus thought to be particularly vulnerable to the effects of hitchhiking or background selection. Here I propose a third explanation for the lower-than-expected levels of genetic diversity within populations of selfing species; recurrent extinctions and recolonizations of local populations, also known as metapopulation dynamics. I show that selfing in a metapopulation setting can result in large reductions in genetic diversity within populations, far greater than expected based the lower effective population size inbreeding species is expected to have. The reason for this depends on an interaction between selfing and pollen migration.     

A threefold genetic allee effect: population size affects cross-compatibility, inbreeding depression and drift load in the self-incompatible Ranunculus reptans

A decline in population size can lead to the loss of allelic variation, increased inbreeding, and the accumulation of genetic load through drift. We estimated the fitness consequences of these processes in offspring of controlled within-population crosses from 13 populations of the self-incompatible, clonal plant Ranunculus reptans. We used allozyme allelic richness as a proxy for long-term population size, which was positively correlated with current population size. Crosses between plants of smaller populations were less likely to be compatible. Inbreeding load, assessed as the slope of the relationship between offspring performance and parental kinship coefficients, was not related to population size, suggesting that deleterious mutations had not been purged from small populations. Offspring from smaller populations were on average more inbred, so inbreeding depression in clonal fitness was higher in small populations. We estimated variation in drift load from the mean fitness of outbred offspring and found enhanced drift load affecting female fertility within small populations. We conclude that self-incompatibility systems do not necessarily prevent small populations from suffering from inbreeding depression and drift load and may exacerbate the challenge of finding suitable mates.     

Heterosis and inbreeding depression in bottlenecked populations: a test in the hermaphroditic freshwater snail Lymnaea stagnalis

Small population size is expected to induce heterosis, due to the random fixation and accumulation of mildly deleterious mutations, whereas within‐population inbreeding depression should decrease due to increased homozygosity. Population bottlenecks, although less effective, may have similar consequences. We tested this hypothesis in the self‐fertile freshwater snail Lymnaea stagnalis, by subjecting experimental populations to a single bottleneck of varied magnitude. Although patterns were not strong, heterosis was significant in the most severely bottlenecked populations, under stressful conditions. This was mainly due to hatching rate, suggesting that early acting and highly deleterious alleles were involved. Although L. stagnalis is a preferential outcrosser, inbreeding depression was very low and showed no clear relationship with bottleneck size. In the less reduced populations, inbreeding depression for hatching success increased under high inbreeding. This may be consistent with the occurence of synergistic epistasis between fitness loci, which may contribute to favour outcrossing in L. stagnalis.     

Genetic implications of mating structure in a Caribbean isolate.

Mating patterns in the population of St. Barthélémy, French West Indies, are examined to evaluate the expected genetic consequences of consanguinity avoidance and of failure to mate due to celibacy and emigration. Nearly 60% of all individuals reaching mating age on "St. Bart" never contribute to the gene pool of succeeding generations. This greatly reduces the effective population size and, therefore, increases genetic drift or random inbreeding. The consequent loss of heterozygosity is partially countered by the fact that the individuals who fail to reproduce within the population tend to be more closely related to the population than are those who do reproduce. This nonrandom failure to mate results in nonrandom inbreeding, which reduces total pedigree inbreeding. Total inbreeding is also reduced by consanguinity avoidance. However, the effect of avoiding close consanguineous matings is, to some extent, reversed by an apparent preference for more distant relationships, especially with second cousins. Generally similar results are found for the related Northside French population of St. Thomas. U.S. Virgin Islands. A recently developed mathematical model predicts the expected effect on heterozygosity of partial or complete avoidance of, or preference for, any combination of consanguineous matings, and also takes into account nonrandom failure to mate. Application of this model to these two populations suggests that under some circumstances nonrandomness in celibacy and emigration may have a greater influence on genotype distributions than does consanguinity avoidance.     

Immigration prevents inbreeding in a growing colony of a long‐lived and philopatric seabird

Mating between close relatives can have deleterious effects on reproductive success or offspring fitness, which should favour the evolution of active or passive inbreeding avoidance mechanisms. In birds, evidence for active inbreeding avoidance by kin‐discriminative mate choice is scarce; many studies describe random mating in relation to kinship and thus support passive inbreeding avoidance by natal dispersal. However, most studies were conducted in island populations of short‐lived passerines with fast alternation of generations. In this study, we present inbreeding estimates based on pedigree data from a 16‐year study in a coastal colony of Common Terns Sterna hirundo, a long‐lived seabird with delayed sexual maturation and low rates of extra‐pair paternity. Incestuous mating was rare (four of 2387 pairs), even if partially accounting for incomplete pedigrees. Although the average relatedness of observed pairs was lower than would be expected from random pairing, the inbreeding coefficient did not differ from random mating. Hence, we found no clear evidence for active inbreeding avoidance by kin‐discriminative mate choice, and the low level of inbreeding seems to be related to the high immigration rate in the colony and thus to be maintained passively by dispersal.     

Pedigree analysis and inbreeding effects over morphological traits in Campolina horse population

Genetic improvement, without control of inbreeding, can go to loss of genetic variability, reducing the potential for genetic gains in the domestic populations. The aim of this study was to analyze the population structure and the inbreeding depression in Campolina horses. Phenotype information from 43 465 individuals was analyzed, data provided by the Campolina Breeders Association. A pedigree file containing 107 951 horses was used to connected the phenotyped individuals. The inbreeding coefficient was performed by use of the diagonal of the relationship matrix and the genealogical parameters were computed using proper softwares. The effective population size was estimated based on the rate of inbreeding and census information, and the stratification of the population was verified by the average relationship coefficient between animals born in different regions of Brazil. The effects of inbreeding on morphological traits were made by the use of inbreeding coefficient as a covariate in the model of random regression. The inbreeding coefficient increased from 1990 on, impacting effective population size and, consequently, shrinking genetic variability. The paternal inbreeding was greater than maternal, which may be attributed to the preference for inbred animals in reproduction. The average genetic relationship coefficient of animals born in different states was lower than individuals born within the same state. The increase in the inbreeding coefficient was negatively associated with all studied traits, showing the importance to avoid genetic losses in the long term. Although results do not indicate a severe narrowing of the population until the present date, the average relationship coefficient shows signs of increase, which could cause a drastic reduction in genetic variability if inbred mating is not successfully controlled in the Campolina horse population.     

Trends in inbreeding, isonymy, and repeated pairs of surnames in the Valserine Valley, French Jura, 1763-1972

The estimates of inbreeding derived from pedigrees and frequency of isonymous marriages (i.e., between persons of the same surname) are compared using genealogical and isonymic information from 4,899 marriages recorded between 1763 and 1972 in 4 rural villages of the French Jura region (a mountainous area near the Swiss border). Before the second half of the 20th century, the two kinds of estimates show a different temporal evolution. The mean inbreeding coefficient based on pedigrees increases between 1763 and 1852 and reaches a maximum between 1853 and 1882 (alpha = 0.0028), with a very low percentage (< 1%) the result of remote kinship. The mean inbreeding coefficient based on isonymy is always higher, with a maximum observed between 1793 and 1822 (F = 0.0200), and it remains roughly the same between 1763 and 1882 (F = 0.0150), with a high percentage resulting from a random component (Fr), a consequence of the small population size and genetic drift. After 1883, the 2 mean coefficients decrease. This discordance is largely explained by the poor quality, for the first periods, of the genealogical data base, which ignores the more remote links of kinship, justifying the use of the model of Crow and Mange (1965) to explore consanguinity during the more ancient periods. The temporal evolution of the repeated pairs of surnames index (RP) confirms the recent evolution of the marital structure of the valley. Moreover, it appears that isonymous marriages and repeated and unique pairs of surnames constitute 3 distinct matrimonial groups characterized by both a different mean coefficient of inbreeding (alpha) and a different rate of endogamy.     

The genetic effective and adult census size of an Australian population of tiger prawns (Penaeus esculentus)

This study compares estimates of the census size of the spawning population with genetic estimates of effective current and long-term population size for an abundant and commercially important marine invertebrate, the brown tiger prawn (Penaeus esculentus). Our aim was to focus on the relationship between genetic effective and census size that may provide a source of information for viability analyses of naturally occurring populations. Samples were taken in 2001, 2002 and 2003 from a population on the east coast of Australia and temporal allelic variation was measured at eight polymorphic microsatellite loci. Moments-based and maximum-likelihood estimates of current genetic effective population size ranged from 797 to 1304. The mean long-term genetic effective population size was 9968. Although small for a large population, the effective population size estimates were above the threshold where genetic diversity is lost at neutral alleles through drift or inbreeding. Simulation studies correctly predicted that under these experimental conditions the genetic estimates would have non-infinite upper confidence limits and revealed they might be overestimates of the true size. We also show that estimates of mortality and variance in family size may be derived from data on average fecundity, current genetic effective and census spawning population size, assuming effective population size is equivalent to the number of breeders. This work confirms that it is feasible to obtain accurate estimates of current genetic effective population size for abundant Type III species using existing genetic marker technology.     

A model for genetic relationship in areally continous plant populations.

Representations are based on plant populations, continuously distributed over their habitats according to specified density patterns. Migration of genetic material takes place via pollen and seed dispersal. Monoecious plants with arbitrary rates of self-fertilization and dioecious plants are considered. The model was constructed with the intention of determining coefficients of inbreeding and kinship for all locations within the seed population after its dispersal over the habitat, assuming the respective genetic relationships of the parental generation to be known. To display the consequences of single components hidden in the general result, the following specifications have been treated: finite population size combined with random dispersal of seed, equilibrium states for hypothetically infinite population size with “limited” dispersal of pollen and seed, random dispersal of pollen, and random dispersal of seed.     

The imprecision of heterozygosity-fitness correlations hinders the detection of inbreeding and inbreeding depression in a threatened species

In nonpedigreed wild populations, inbreeding depression is often quantified through the use of heterozygosity-fitness correlations (HFCs), based on molecular estimates of relatedness. Although such correlations are typically interpreted as evidence of inbreeding depression, by assuming that the marker heterozygosity is a proxy for genome-wide heterozygosity, theory predicts that these relationships should be difficult to detect. Until now, the vast majority of empirical research in this area has been performed on generally outbred, nonbottlenecked populations, but differences in population genetic processes may limit extrapolation of results to threatened populations. Here, we present an analysis of HFCs, and their implications for the interpretation of inbreeding, in a free-ranging pedigreed population of a bottlenecked species: the endangered takahe (Porphyrio hochstetteri). Pedigree-based inbreeding depression has already been detected in this species. Using 23 microsatellite loci, we observed only weak evidence of the expected relationship between multilocus heterozygosity and fitness at individual life-history stages (such as survival to hatching and fledging), and parameter estimates were imprecise (had high error). Furthermore, our molecular data set could not accurately predict the inbreeding status of individuals (as 'inbred' or 'outbred', determined from pedigrees), nor could we show that the observed HFCs were the result of genome-wide identity disequilibrium. These results may be attributed to high variance in heterozygosity within inbreeding classes. This study is an empirical example from a free-ranging endangered species, suggesting that even relatively large numbers (>20) of microsatellites may give poor precision for estimating individual genome-wide heterozygosity. We argue that pedigree methods remain the most effective method of quantifying inbreeding in wild populations, particularly those that have gone through severe bottlenecks.     

Genetic kinship among an isolated Adi tribe of Arunachal Pradesh: isonymy in the Adi Panggi

Isolated tribes in remote areas are important for genetic studies, and one such little known subtribe of the Adi tribe, namely, the Adi Panggi (Pangi) of the Upper Siang District of Arunachal Pradesh, India, was studied for surname distribution to deduce the deviation from random mating and genetic kinship between villages. The estimates of homonymy (homozygosity) vary between villages; husbands show wider variation (0.009 to 0.23) than wives (0.005 to 0.054). The remote villages of Sumsing and Sibum and Geku Town show lower entropy among husbands' surnames than among Panggi wives. The highest equivalent surname number was found among Sibum husbands (9.9), Panggi wives (12.6), and Panggi and non-Panggi wives (13.5). The estimates of unbiased random isonymy among husbands and wives together show the smallest values in Sibum (0.05) and the highest values in Sumsing and Ramku (0.16). The random and nonrandom components of the inbreeding coefficient show avoidance of inbreeding among the Panggi villages (-0.012 to -0.27) except in Sibum (0.012). Genetic kinship between villages based on the Mij distance shows different clusters of villages among husbands and wives. Both the Panggi wives and the Panggi and non-Panggi wives show a similar pattern of clustering between villages. The wide homonymy variation between villages among the patrilocal Adi Panggi indicates differential genetic kinetics among husbands and wives, avoidance of inbreeding, and female-oriented differential gene flow with little effect on the overall intervillage genetic kinship.     

No evidence for inbreeding avoidance in a natural population of song sparrows (Melospiza melodia)

We studied mate choice and inbreeding avoidance a natural population of song sparrows (Melospiza melodia) on Mandarte Island, Canada. Inbreeding occurred regularly: 59% all matings were between known relatives. We tested for inbreeding avoidance by comparing the observed levels of inbreeding to those expected if mate choice had been random with respect to relatedness. Independent of our assumptions about the availability of mates in the random mating model, we found that the expected and observed distributions of inbreeding coefficients were similar, as was the expected and observed frequency of close (f >/= 0.125) inbreeding. Furthermore, there was no difference in relatedness observed pairs and those that would have resulted had birds mated instead with their nearest neighbors. The only evidence to suggest any inbreeding avoidance was a reduced rate of parent-offspring matings as compared to one random mating model but not the other. Hence, despite substantial inbreeding depression in this population, we found little evidence for inbreeding avoidance through mate choice. We present a simple model to suggest that variation in inbreeding avoidance behaviors in birds may arise from differences in survival rates: in species with low survival rates, the costs of forfeiting matings to avoid inbreeding may exceed the costs of inbreeding.     

Genetic and demographic dynamics of small populations of Silene latifolia

Small local populations of Silene alba, a short-lived herbaceous plant, were sampled in 1994 and again in 1999. Sampling included estimates of population size and genetic diversity, as measured at six polymorphic allozyme loci. When averaged across populations, there was very little change between samples (about three generations) in population size, measures of within-population genetic diversity such as number of alleles or expected heterozygosity, or in the apportionment of genetic diversity within and among populations as measured by F(st). However, individual populations changed considerably, both in terms of numbers of individuals and genetic composition. Some populations doubled in size between samples, while others shrank by more than 75%. Similarly, expected heterozygosity and allele number increased by more than two-fold in individual populations and decreased by more than three-fold in others. When population-specific change in number and change in measures of genetic diversity were considered together, significant positive correlations were found between the demographic and genetic variables. It is speculated that some populations were released from the demographic consequences of inbreeding depression by gene flow.     

Interpopulation variation in inbreeding is primarily driven by tolerance of mating with relatives in a spermcasting invertebrate

The degree to which individuals inbreed is a fundamental aspect of population biology shaped by both passive and active processes. Yet, the relative influences of random and non-random mating on the overall magnitude of inbreeding are not well characterized for many taxa. We quantified variation in inbreeding among qualitatively accessible and isolated populations of a sessile marine invertebrate (the colonial ascidian Lissoclinum verrilli) in which hermaphroditic colonies cast sperm into the water column for subsequent uptake and internal fertilization. We compared estimates of inbreeding to simulations predicting random mating within sites to evaluate if levels of inbreeding were (1) less than expected because of active attempts to limit inbreeding, (2) as predicted by genetic subdivision and passive inbreeding tolerance, or (3) greater than simulations due to active attempts to promote inbreeding via self-fertilization or a preference for related mates. We found evidence of restricted gene flow and significant differences in the genetic diversity of L. verrilli colonies among sites, indicating that on average colonies were weakly related in accessible locations, but their levels of relatedness matched that of first cousins or half-siblings on isolated substrates. Irrespective of population size, progeny arrays revealed variation in the magnitude of inbreeding across sites that tracked with the mean relatedness of conspecifics. Biparental reproduction was confirmed in most offspring (86%) and estimates of total inbreeding largely overlapped with simulations of random mating, suggesting that interpopulation variation in mother–offspring resemblance was primarily due to genetic subdivision and passive tolerance of related mates. Our results highlight the influence of demographic isolation on the genetic composition of populations, and support theory predicting that tolerance of biparental inbreeding, even when mates are closely related, may be favoured under a broad set of ecological and evolutionary conditions.     

Local extinction and recolonization, species effective population size, and modern human origins

A primary objection from a population genetics perspective to a multiregional model of modern human origins is that the model posits a large census size, whereas genetic data suggest a small effective population size. The relationship between census size and effective size is complex, but arguments based on an island model of migration show that if the effective population size reflects the number of breeding individuals and the effects of population subdivision, then an effective population size of 10,000 is inconsistent with the census size of 500,000 to 1,000,000 that has been suggested by archeological evidence. However, these models have ignored the effects of population extinction and recolonization, which increase the expected variance among demes and reduce the inbreeding effective population size. Using models developed for population extinction and recolonization, we show that a large census size consistent with the multiregional model can be reconciled with an effective population size of 10,000, but genetic variation among demes must be high, reflecting low interdeme migration rates and a colonization process that involves a small number of colonists or kin-structured colonization. Ethnographic and archeological evidence is insufficient to determine whether such demographic conditions existed among Pleistocene human populations, and further work needs to be done. More realistic models that incorporate isolation by distance and heterogeneity in extinction rates and effective deme sizes also need to be developed. However, if true, a process of population extinction and recolonization has interesting implications for human demographic history.