Gene genealogies in geographically structured populations.

Population genetics theory has dealt only with the spatial or geographic pattern of degrees of relatedness or genetic similarity separately for each point in time. However, a frequent goal of experimental studies is to infer migration patterns that occurred in the past or over extended periods of time. To fully understand how a present geographic pattern of genetic variation reflects one in the past, it is necessary to build genealogy models that directly relate the two. For the first time, space-time probabilities of identity by descent and coalescence probabilities are formulated and characterized in this article. Formulations for general migration processes are developed and applied to specific types of systems. The results can be used to determine the level of certainty that genes found in present populations are descended from ancient genes in the same population or nearby populations vs. geographically distant populations. Some parameter combinations result in past populations that are quite distant geographically being essentially as likely to contain ancestors of genes at a given population as the past population located at the same place. This has implications for the geographic point of origin of ancestral, "Eve," genes. The results also form the first model for emerging "space-time" molecular genetic data.     

Evolutionary and statistical properties of three genetic distances

Many genetic distances have been developed to summarize allele frequency differences between populations. I review the evolutionary and statistical properties of three popular genetic distances: DS, DA, and theta;, using computer simulation of two simple evolutionary histories: an isolation model of population divergence and an equilibrium migration model. The effect of effective population size, mutation rate, and mutation mechanism upon the parametric value between pairs of populations in these models explored, and the unique properties of each distance are described. The effect of these evolutionary parameters on study design is also investigated and similar results are found for each genetic distance in each model of evolution: large sample sizes are warranted when populations are relatively genetically similar; and loci with more alleles produce better estimates of genetic distance.     

Dissecting genetic networks underlying complex phenotypes: the theoretical framework

Great progress has been made in genetic dissection of quantitative trait variation during the past two decades, but many studies still reveal only a small fraction of quantitative trait loci (QTLs), and epistasis remains elusive. We integrate contemporary knowledge of signal transduction pathways with principles of quantitative and population genetics to characterize genetic networks underlying complex traits, using a model founded upon one-way functional dependency of downstream genes on upstream regulators (the principle of hierarchy) and mutual functional dependency among related genes (functional genetic units, FGU). Both simulated and real data suggest that complementary epistasis contributes greatly to quantitative trait variation, and obscures the phenotypic effects of many 'downstream' loci in pathways. The mathematical relationships between the main effects and epistatic effects of genes acting at different levels of signaling pathways were established using the quantitative and population genetic parameters. Both loss of function and "co-adapted" gene complexes formed by multiple alleles with differentiated functions (effects) are predicted to be frequent types of allelic diversity at loci that contribute to the genetic variation of complex traits in populations. Downstream FGUs appear to be more vulnerable to loss of function than their upstream regulators, but this vulnerability is apparently compensated by different FGUs of similar functions. Other predictions from the model may account for puzzling results regarding responses to selection, genotype by environment interaction, and the genetic basis of heterosis.     

Migration and isolation effects on the Rapanui population (Easter island) through the analysis of STRs on the Y chromosome

Cultural and biological data suggests the Polynesian origin of the Rapanui population, although the presence of foreign genes in the native population, as a result of admixture with Europeans in the last two centuries has also to be considered. In order to estimate the genetic affinities of the present-day inhabitants of Easter Island and the nearby populations, we used seven polymorphisms of the Y chromosome. However we want to estimate the grade of admixture on the genetic structure that was brought from foreigners within the last two centuries upon the more geographically isolated populations in the world. The preliminary results showed the presence of 18 haplotypes analyzed on 30 male samples. The analysis of the allelic frequency showed a distribution typical of the Polynesian populations. Available data in literature, even with some differences probably due to either the founder effect or historical and ecological events that created sudden demographic variations on the island population. The phylogentic analysis of the haplotypes obtained through Network Median Joining showed two different cluster of haplotypes, of which one represents about 64% of the present haplotypes on Easter Island, which are characterized from the presence of the allele DYS19^*16, very frequent in the Pacific populations. The other cluster is characterized from the presence of the allele DYS19^*14, absent within the populations in the Pacific and with reasonable high frequency within the European populations and South American. Most probably the two clusters are the product of several colonizations that Easter Island had endured from the time of the Chilian and European Colonies. It was demonstrated in fact, that the arrival in 1914 of 50 German and English prisoners would have left a considerable genetic impact on the population of Rapanui, which during this period was of small size.     

Reinforcement and divergence under assortative mating

Traits that cause assortative mating such as the flowering time in plants and body size in animals can produce reproductive isolation between hybridizing populations. Can selection against unfit hybrids cause two populations to diverge in their mean values for these kinds of traits? Here I present a haploid analytical model of one population that receives gene flow from another. The partial pre-zygotic isolation between the two populations is caused by assortative mating for a trait that is influenced by any number of genes with additive effects. The post-zygotic isolation is caused by selection against genetic incompatibilities that can involve any form of selection on individual genes and gene combinations (epistasis). The analysis assumes that the introgression rate and selection coefficients are small. The results show that the assortment trait mean will not diverge from the immigrants unless there is direct selection on the trait favouring it to do so or there are genes of very large effect. The amount of divergence at equilibrium is determined by a balance between direct selection on the assortment trait and introgression from the other population. Additional selection against hybrid genetic incompatibilities reduces the effective migration rate and allows greater divergence. The role of assortment in speciation is discussed in the light of these results.     

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.     

Inter‐island divergence within Drosophila mauritiana,a species of the D. simulans complex: Past history and/or speciation in progress?

Speciation with gene flow may be more common than generally thought, which makes detailed understanding of the extent and pattern of genetic divergence between geographically isolated populations useful. Species of the Drosophila simulans complex provide a good model for speciation and evolutionary studies, and hence understanding their population genetic structure will increase our understanding of the context in which speciation has occurred. Here, we describe genetic diversity and genetic differentiation of two distant populations of D. mauritiana (Mauritius and Rodrigues Islands) at mitochondrial and nuclear loci. We surveyed the two populations for their mitochondrial haplotypes, eight nuclear genes and 18 microsatellite loci. A new mitochondrial type is fixed in the Rodrigues population of D. mauritiana. The two populations are highly differentiated, their divergence appears relatively ancient (100,000 years) compared to the origin of the species, around 0.25MYA, and they exhibit very limited gene flow. However, they have similar levels of divergence from their sibling, D. simulans. Both nuclear genes and microsatellites revealed contrasting demographic histories between the two populations, expansion for the Mauritius population and stable population size for the Rodrigues Island population. The discovery of pronounced geographic structure within D. mauritiana combined to genetic structuring and low gene flow between the two island populations illuminates the evolutionary history of the species and clearly merits further attention in the broad context of speciation.     

Extinction of populations due to inbreeding depression with demographic disturbances

The process of population extinction due to inbreeding depression with constant demographic disturbances every generation is analysed using a population genetic and demographic model. The demographic disturbances introduced into the model represent loss of population size that is induced by any kind of human activities, e.g. through hunting and destruction of habitats. The genetic heterozygosity among recessive deleterious genes and the population size are assumed to be in equilibrium before the demographic disturbances start. The effects of deleterious mutations are represented by decreases in the growth rate and carrying capacity of a population. Numerical simulations indicate rapid extinction due to synergistic interaction between inbreeding depression and declining population size for realistic ranges of per-locus mutation rate, equilibrium population size, intrinsic rate of population growth, and strength of demographic disturbances. Large populations at equilibrium are more liable to extinction when disturbed due to inbreeding depression than small populations. This is a consequence of the fact that large populations maintain more recessive deleterious mutations than small populations. The rapid extinction predicted in the present study indicates the importance of the demographic history of a population in relation to extinction due to inbreeding depression.     

Sex ratio rather than population size affects genetic diversity in Antennaria dioica

• Habitat fragmentation and small population size can lead to genetic erosion in threatened plant populations. Classical theory implies that dioecy can counteract genetic erosion as it decreases the magnitude of inbreeding and genetic drift due to obligate outcrossing. However, in small populations, sex ratios may be strongly male‐ or female‐biased, leading to substantial reductions in effective population size. This may theoretically result in a unimodal relationship between sex ratios and genetic diversity; yet, empirical studies on this relationship are scarce.
• Using AFLP markers, we studied genetic diversity, structure and differentiation in 14 highly fragmented Antennaria dioica populations from the Central European lowlands. Our analyses focused on the relationship between sex ratio, population size and genetic diversity.
• Although most populations were small (mean: 35.5 patches), genetic diversity was moderately high. We found evidence for isolation‐by‐distance, but overall differentiation of the populations was rather weak. Females dominated 11 populations, which overall resulted in a slightly female‐biased sex ratio (61.5%). There was no significant relationship between population size and genetic diversity. The proportion of females was not unimodally but positively linearly related to genetic diversity.
• The high genetic diversity and low genetic differentiation suggest that A. dioica has been widely distributed in the Central European lowlands in the past, while fragmentation occurred only in the last decades. Sex ratio has more immediate consequences on genetic diversity than population size. An increasing proportion of females can increase genetic diversity in dioecious plants, probably due to a higher amount of sexual reproduction.

     

Variance of genetic distance and correlation of heterozygosity between populations under the pressure of stepwise mutation

Using the stepwise mutation model of Ohta and Kimura (1973), formulas are developed for the correlation of heterozygosity and the variance of genetic distance between two finite populations. Studied in detail is the case where the sizes of the two descendant populations are equal to that of the ancestral population and the mutation rate is the same for all loci. Numerical computations are carried out by using the present formulas and those of Li and Nei (1975Genet. Res.25) for the infinite-allele model. The results are as follows: The correlation of heterozygosity decreases with time faster for the stepwise mutation model than for the infinite-allele model. However, the relationships between the correlation of heterozygosity and the normalized genetic identity for the two models are very similar, if the average heterozygosities of the two populations are around 0.20 or less. On the other hand, the variance of genetic distance for the stepwise mutation model may become considerably smaller than that for the infinite-allele model, if the average heterozygosities of the two populations are larger than 0.05. The ratio of the standard deviation to the mean is, however, very large for the stepwise mutation model as well as the infinite-allele model.     

On the size distribution of private microsatellite alleles

Private microsatellite alleles tend to be found in the tails rather than in the interior of the allele size distribution. To explain this phenomenon, we have investigated the size distribution of private alleles in a coalescent model of two populations, assuming the symmetric stepwise mutation model as the mode of microsatellite mutation. For the case in which four alleles are sampled, two from each population, we condition on the configuration in which three distinct allele sizes are present, one of which is common to both populations, one of which is private to one population, and the third of which is private to the other population. Conditional on this configuration, we calculate the probability that the two private alleles occupy the two tails of the size distribution. This probability, which increases as a function of mutation rate and divergence time between the two populations, is seen to be greater than the value that would be predicted if there was no relationship between privacy and location in the allele size distribution. In accordance with the prediction of the model, we find that in pairs of human populations, the frequency with which private microsatellite alleles occur in the tails of the allele size distribution increases as a function of genetic differentiation between populations.     

WAITING FOR SPECIATION: THE EFFECT OF POPULATION SUBDIVISION ON THE TIME TO SPECIATION

We study the time required for speciation in a species that is divided into small versus large populations. Following Dobzhansky and Muller, we assume that hybrid sterility or inviability is caused by “complementary genes,” that is, by the accumulation of genes that cause sterility or inviability when brought together in hybrids but that have no deleterious effect on their normal species genetic background. When divergence between populations is caused by genetic drift, we show that the time to speciation is independent of population subdivision: speciation occurs just as quickly in a species split into a few large populations as into many small populations. When divergence is driven by natural selection, however, the time to speciation is very sensitive to population subdivision and speciation occurs most rapidly when a species is split into two large populations. These results contradict several popular intuitions about the effect of population size on speciation.     

Functional genomics of life history variation in a butterfly metapopulation

In fragmented landscapes, small populations frequently go extinct and new ones are established with poorly understood consequences for genetic diversity and evolution of life history traits. Here, we apply functional genomic tools to an ecological model system, the well-studied metapopulation of the Glanville fritillary butterfly. We investigate how dispersal and colonization select upon existing genetic variation affecting life history traits by comparing common-garden reared 2-day adult females from new populations with those from established older populations. New-population females had higher expression of abdomen genes involved in egg provisioning and thorax genes involved in the maintenance of flight muscle proteins. Physiological studies confirmed that new-population butterflies have accelerated egg maturation, apparently regulated by higher juvenile hormone titer and angiotensin converting enzyme mRNA, as well as enhanced flight metabolism. Gene expression varied between allelic forms of two metabolic genes (Pgi and Sdhd), which themselves were associated with differences in flight metabolic rate, population age and population growth rate. These results identify likely molecular mechanisms underpinning life history variation that is maintained by extinction-colonization dynamics in metapopulations.     

Relationships between Population Size, Genetic Diversity and Fitness Components in the Rare Plant Dictamnus albus in Central Germany

An increasing number of studies support the hypothesis that smaller populations face a higher risk of extinction, and declining population sizes are therefore one of the focal points in plant conservation. In small populations, loss of genetic diversity is often related to reduced reproductive fitness. For the rare Dictamnus albus in Central Germany, an earlier study had already confirmed a significant correlation between population size and genetic diversity. In order to assess whether these variables correlate with fitness components, plant height; flower, fruit and seed production; and germination were studied in a total of 11 populations of different size. In the seven populations that were sampled over two consecutive years, differences among populations and among years were tested using a Two-Way ANOVA. Co-linearity among variables was assessed using principal component analysis (PCA), followed by calculating correlations between ordination axes and both population size and genetic diversity. Plant height and flower number were uncorrelated to the other variables and, together with germination, did not show any correlation to either population size or genetic diversity. However, both size and genetic diversity of populations correlated significantly with other PCA axes that reflected reproductive components such as fruit number, seed number, seed fruit ratio, and seed mass. Our results support the idea that reproduction is hampered in small populations and raise concerns over the loss of genetic diversity in D. albus.     

Mapping quantitative trait loci in multiple populations of Arabidopsis thaliana identifies natural allelic variation for trichome density

The majority of biological traits are genetically complex. Mapping the quantitative trait loci (QTL) that determine these phenotypes is a powerful means for estimating many parameters of the genetic architecture for a trait and potentially identifying the genes responsible for natural variation. Typically, such experiments are conducted in a single mapping population and, therefore, have only the potential to reveal genomic regions that are polymorphic between the progenitors of the population. What remains unclear is how well the QTL identified in any one mapping experiment characterize the genetics that underlie natural variation in traits. Here we provide QTL mapping data for trichome density from four recombinant inbred mapping populations of Arabidopsis thaliana. By aligning the linkage maps for these four populations onto a common physical map, the results from each experiment were directly compared. Seven of the nine QTL identified are population specific while two were mapped in all four populations. Our results show that many lineage-specific alleles that either increase or decrease trichome density persist in natural populations and that most of this genetic variation is additive. More generally, these findings suggest that the use of multiple populations holds great promise for better understanding the genetic architecture of natural variation.     

Combining genetic structure and demographic analyses to estimate persistence in endangered Key deer (<Emphasis Type="Italic">Odocoileus virginianus clavium</Emphasis>)

Recent improvements in genetic analyses have paved the way in using molecular data to answer questions regarding evolutionary history, genetic structure, and demography. Key deer are a federally endangered subspecies assumed to be genetically unique, homogeneous, and have a female-biased population of approximately 900 deer. We used 985 bp of the mitochondrial cytochrome b gene and 12 microsatellite loci to test two hypotheses: (1) that Key deer are isolated and have reduced diversity compared to mainland deer and (2) that isolation of the Florida Keys has led to a small population size and a high risk of extinction. Our results indicate that Key deer are indeed genetically isolated from mainland white-tailed deer and that there is a lack of genetic substructure between islands. While Key deer exhibit reduced levels of genetic diversity compared to their mainland counterparts, they contain enough diversity to uniquely identify individual deer. Based on genetic identification, we estimated a census size of around 1000 individuals with a heavily skewed female-biased adult sex ratio. Furthermore, our genetic and contemporary demographic data were used to generate a species persistence model of the Key deer. Sensitivity tests within the population viability analysis brought to light the importance of fetal sex ratio and female survival as the primary factors at risk of driving the subspecies to extinction. This study serves as a prime example of how persistence models can be used to evaluate population viability in natural populations of endangered organisms.     

Genetic estimates of contemporary effective population size: what can they tell us about the importance of genetic stochasticity for wild population persistence?

Genetic stochasticity due to small population size contributes to population extinction, especially when population fragmentation disrupts gene flow. Estimates of effective population size ( N _e) can therefore be informative about population persistence, but there is a need for an assessment of their consistency and informative relevance. Here we review the body of empirical estimates of N _e for wild populations obtained with the temporal genetic method and published since Frankham's (1995 ) review. Theoretical considerations have identified important sources of bias for this analytical approach, and we use empirical data to investigate the extent of these biases. We find that particularly model selection and sampling require more attention in future studies.
We report a median unbiased N _e estimate of 260 (among 83 studies) and find that this median estimate tends to be smaller for populations of conservation concern, which may therefore be more sensitive to genetic stochasticity. Furthermore, we report a median N _e/ N ratio of 0.14, and find that this ratio may actually be higher for small populations, suggesting changes in biological interactions at low population abundances. We confirm the role of gene flow in countering genetic stochasticity by finding that N _e correlates strongest with neutral genetic metrics when populations can be considered isolated. This underlines the importance of gene flow for the estimation of N _e, and of population connectivity for conservation in general. Reductions in contemporary gene flow due to ongoing habitat fragmentation will likely increase the prevalence of genetic stochasticity, which should therefore remain a focal point in the conservation of biodiversity.     

Complex genetic architecture of population differences in adult lifespan of a beetle: nonadditive inheritance, gender differences, body size and a large maternal effect

Evolutionary responses to selection can be complicated when there is substantial nonadditivity, which limits our ability to extrapolate from simple models of selection to population differentiation and speciation. Studies of Drosophila melanogaster indicate that lifespan and the rate of senescence are influenced by many genes that have environment- and sex-specific effects. These studies also demonstrate that interactions among alleles (dominance) and loci (epistasis) are common, with the degree of interaction differing between the sexes and among environments. However, little is known about the genetic architecture of lifespan or mortality rates for organisms other than D. melanogaster. We studied genetic architecture of differences in lifespan and shapes of mortality curves between two populations of the seed beetle, Callosobruchus maculatus (South India and Burkina Faso populations). These two populations differ in various traits (such as body size and adult lifespan) that have likely evolved via host-specific selection. We found that the genetic architecture of lifespan differences between populations differs substantially between males and females; there was a large maternal effect on male lifespan (but not on female lifespan), and substantial dominance of long-life alleles in females (but not males). The large maternal effect in males was genetically based (there was no significant cytoplasmic effect) likely due to population differences in maternal effects genes that influence lifespan of progeny. Rearing host did not affect the genetic architecture of lifespan, and there was no evidence that genes on the Y-chromosome influence the population differences in lifespan. Epistatic interactions among loci were detectable for the mortality rate of both males and females, but were detectable for lifespan only after controlling for body size variation among lines. The detection of epistasis, dominance, and sex-specific genetic effects on C. maculatus lifespan is consistent with results from line cross and quantitative trait locus studies of D. melanogaster.     

Huge populations and old species of Costa Rican and Panamanian dirt frogs inferred from mitochondrial and nuclear gene sequences

Molecular genetic data were used to investigate population sizes and ages of Eleutherodactylus (Anura: Leptodactylidae), a species-rich group of small leaf-litter frogs endemic to Central America. Population genetic structure and divergence was investigated for four closely related species surveyed across nine localities in Costa Rica and Panama. DNA sequence data were collected from a mitochondrial gene (ND2) and a nuclear gene (c-myc). Phylogenetic analyses yielded concordant results between loci, with reciprocal monophyly of mitochondrial DNA haplotypes for all species and of c-myc haplotypes for three of the four species. Estimates of genetic differentiation among populations (FST) based upon mitochondrial data were always higher than nuclear-based FST estimates, even after correcting for the expected fourfold lower effective population size (Ne) of the mitochondrial genome. Comparing within-population variation and the relative mutation rates of the two genes revealed that the Ne of the mitochondrial genome was 15-fold lower than the estimate of the nuclear genome based on c-myc. Nuclear FST estimates were approximately 0 for the most proximal pairs of populations, but ranged from 0.5 to 1.0 for all other pairs, even within the same nominal species. The nuclear locus yielded estimates of Ne within localities on the order of 105. This value is two to three orders of magnitude larger than any previous Ne estimate from frogs, but is nonetheless consistent with published demographic data. Applying a molecular clock model suggested that morphologically indistinguishable populations within one species may be 107 years old. These results demonstrate that even a geologically young and dynamic region of the tropics can support very old lineages that harbour great levels of genetic diversity within populations. The association of high nucleotide diversity within populations, large divergence between populations, and high species diversity is also discussed in light of neutral community models.     

Fluctuating Selection,Egg Banks and Population Genetic Structure in Cyclically Parthenogenetic Species

Associations between neutral genetic markers and genes under selection have been suggested to explain the population genetic structure of neutral genes in cyclically parthenogenetic freshwater invertebrates. A simulation model was constructed in order to analyse the extrapolated consequences of observed fluctuations in genotype frequencies in Daphnia, in the presence of egg banks. When of sufficient depth and magnitude, egg banks in combination with fluctuating selection were shown to maintain genetic variation within populations indefinitely. The level of equilibrium diversity increased with the depth and magnitude of the banks, and with intensity of selection. The same threshold was responsible for genetic differentiation between populations, which was independent of migration rate, and which was attained very rapidly following initial Hardy–Weinberg equilibrium. In the absence of selection, egg banks increased the effective size of local populations, thereby decreasing genetic differentiation at migration-drift equilibrium. These results suggest that egg banks are crucial to the genetic structure in the presence of fluctuating selective pressures, but more data are needed if this knowledge is to be used in an improved general understanding of the genetic structure of cyclically parthenogenetic species.