The effects of dominance, regular inbreeding and sampling design on Q(ST), an estimator of population differentiation for quantitative traits.

To test whether quantitative traits are under directional or homogenizing selection, it is common practice to compare population differentiation estimates at molecular markers (F(ST)) and quantitative traits (Q(ST)). If the trait is neutral and its determinism is additive, then theory predicts that Q(ST) = F(ST), while Q(ST) > F(ST) is predicted under directional selection for different local optima, and Q(ST) < F(ST) is predicted under homogenizing selection. However, nonadditive effects can alter these predictions. Here, we investigate the influence of dominance on the relation between Q(ST) and F(ST) for neutral traits. Using analytical results and computer simulations, we show that dominance generally deflates Q(ST) relative to F(ST). Under inbreeding, the effect of dominance vanishes, and we show that for selfing species, a better estimate of Q(ST) is obtained from selfed families than from half-sib families. We also compare several sampling designs and find that it is always best to sample many populations (>20) with few families (five) rather than few populations with many families. Provided that estimates of Q(ST) are derived from individuals originating from many populations, we conclude that the pattern Q(ST) > F(ST), and hence the inference of directional selection for different local optima, is robust to the effect of nonadditive gene actions.

Population differentiation and its consequences in the white campion "Silene latifolia" (Caryophyllaceae): a functional analysis

Summary In his theory On the Origin of Species by Means of Natural Selection (1859), Darwin describes evolution as a gradual change in population over time and that natural selection is a process that caused evolution. Because quantitative variation in species is partly influenced by several genes and thus heritable, association between levels of genetic variation at neutral markers and at quantitative traits and their partitioning within and among populations are important to study mechanisms that drive evolution in populations. Most studies addressing quantitative variation in plants focused on morphological and life history traits but not in traits affecting reproductive success. The aim of this thesis is to better understand how patterns of variation for neutral molecular markers and phenotypic traits drive the evolution of reproduction and defensive mechanisms in six European populations of Silene latifolia, a dioecious plant species. We found evidence for extremely high within and between population variation at six microsatellite loci and at most quantitative traits studied in plants grown under standardized conditions (morphology, life history and reproductive traits). Interestingly, there was clinal variation between age at first flowering and latitude. This pattern is likely due to natural selection since differentiation of this trait was high, heritable and probably higher than differentiation at neutral markers. Our study focused on sex specific selective pressures: mechanisms of intersexual coadaptation and defence mechanism against the seed predator Hadena bicruris. To address divergence at reproductive traits, we studied male and female population of origin effects and in particular pollen competitive ability on male post-pollination success in the study populations with within and between populations crosses. We crossed the same female plant with pollen from a male within the same population of origin and pollen from two males from two distinct populations, using a fixed tester male as a competitor. Additionally, we conducted control crosses with pollen from each male as a single donor. We analysed paternity success of each competitor with two microsatellite loci, seed set and offspring fitness. Male population of origin showed significant among-population variation for siring success at pollen competition. In vitro pollen germination rate showed heritable variation among populations and was positively correlated to siring success. Local or foreign pollen did not have a consistent advantage. Furthermore, female population of origin affected the outcome of pollen competition in some populations. There was no difference of seed set or offspring fitness in within/ between population crosses. This suggests that reproductive divergence may occur via pollen competition in Silene latifolia. The specialist seed predator Hadena bicruris may also induce divergence between populations. We tested potential constitutive and induced defence mechanisms against the specialist predator Hadena bicruris. Because fruit wall thickness is smaller in the invasive range (Northern America) were the moth is absent, this suggests that a thicker fruit wall is a potentially defensive trait against larval attack, and that relaxed selection in the absence of the seed predator has resulted in an evolutionary loss of this defence in the invasive range. Fruit wall thickness was different among three populations. Experimental exposure to moth eggs increased fruit abortion. Fruits built after attack on exposed plants did not have thicker fruit walls compared to fruits on non-exposed plants. Furthermore, fruits with thicker fruit walls were not less profitable, nor did they require longer handling time when exposed to larvae, suggesting no defensive role of fruit wall thickness. Our results show that there is high molecular and phenotypic variation in Silene latifolia and that traits potentially involved in reproductive success both for intra-specific (between sexes) and inter-specific interactions are heritable. Different selective forces may thus interact and cause differential evolution of geographically separated Silene latifolia populations in Europe, leading to the observed differentiation. Résumé Dans sa théorie de l'évolution, L'origine des espèces, ch. 4 (1859), Darwin décrit l'évolution comme un processus continu au cours du temps à l'intérieur de populations et que la sélection naturelle en est le moteur. La variation quantitative est en partie déterminée par plusieurs gènes, donc transmissible à la descendance. Associer le niveau de variation génétique à des marqueurs neutres au niveau de la variation à des traits quantitatifs, ainsi que la répartition à l'intérieur et entre les populations d'une espèce donnée de cette variation, sont importants dans la compréhension des forces évolutives. La plupart des études scientifiques sur la variation quantitative chez les plantes se sont intéressées à la morphologie et à la phénologie mais pas aux caractères impliqués dans le succès reproducteur. L'objectif de cette thèse est de mieux comprendre comment la répartition de la variation à des marqueurs neutres et des caractères quantitatifs influence l'évolution de la reproduction et des mécanismes de défense dans six populations Européennes de l'espèce dioïque Silene latifolia. Nous avons mis en évidence une grande diversité intra et inter-population à six loci microsatellites ainsi qu'à la plupart des caractères quantitatifs mesurés (morphologie, phénologie et traits reproducteurs) sur des plantes cultivées dans des conditions standardisées. Un résultat intéressant est la présence d'un cline latitudinal pour l'âge à la floraison. Ceci est probablement une conséquence de la sélection naturelle, puisque ce caractère est différencié entre les populations étudiées, héritable et que la différenciation de ce trait est supérieure à la différenciation des marqueurs neutres étudiés. Notre étude a ensuite porté plus précisément sur les pressions de sélection spécifiques aux sexes : la coadaptation entre les sexes et les mécanismes de défense contre l'insecte granivore Hadena bicruris. Afin d'évaluer la divergence sur les traits reproducteurs, nous avons étudié les effets des populations d'origine des mâles et des femelles et en particulier le succès reproducteur des mâles après pollinisation à l'aide de croisements inter et intra-population. Nous avons pollinisé la même femelle avec du pollen provenant d'un mâle de la même population ainsi qu'avec le pollen de deux mâles provenant de deux autres populations en situation de compétition avec un pollen provenant d'une population test. Des croisements contrôle ont été réalisés avec les mêmes mâles en pollinisation pure. Nous avons évalué le succès reproducteur de chaque mâle à l'aide d'analyses de paternité ainsi que la production de graines et la fitness de la descendance. L'origine du mâle avait un effet sur la paternité. Le taux de croissance in vitro du pollen est un caractère héritable et a eu un effet positif sur le succès reproducteur. De plus, l'origine de la femelle avait un effet sur le succès des mâles en compétition dans certaines populations. Nos résultats suggèrent qu'une divergence reproductive chez Silene latifolia pourrait apparaître suite à la compétition pollinique. Nous avons ensuite testé des mécanismes potentiels de défense constitutive et induite contre l'herbivore spécialiste Hadena bicruris, un papillon nocturne qui pourrait aussi jouer un rôle dans la différenciation des populations. L'épaisseur des fruits étant plus faible dans les régions où la plante est invasive (Amérique du Nord) et où l'insecte est absent, ce trait pourrait jouer un rôle défensif. Une pression de sélection plus faible causée par l'absence de l'herbivore aurait abouti à une perte de cette défense dans ces régions. Nous avons montré que l'épaisseur du fruit est variable selon les populations. L'infestation artificielle de fruit par l'insecte induit l'abscission sélective des fruits. Les fruits produits après une infestation n'étaient pas plus épais que les fruits issus de plantes non infestées. De plus, les fruits épais n'étaient pas moins nutritifs et ne causaient pas de perte de temps pour la prédation pour les larves, ce qui suggère que l'épaisseur des fruits ne joue pas un rôle défensif. Nos résultats montrent que plusieurs pressions de sélection interviennent et interagissent dans l'évolution de populations distantes, provoquant la divergence des populations Européennes de l'espèce Silene latifolia.

Extreme population differentiation in vulnerable slavemaking ant with fragmented distribution

Understanding levels of population differentiation and inbreeding are important issues in conservation biology, especially for social Hymenoptera with fragmented and small population sizes. Isolated populations are more vulnerable to genetic loss and extinction than those with extended continuous distributions. However, small populations are not always a consequence of a recent reduction of their habitat. Thus, determining the history of population isolation and current patterns of genetic variation of a species is crucial for its conservation. Rossomyrmex minuchae is a slave-making ant with patchy distribution in South Eastern Spain and is classified as vulnerable by the IUCN. In contrast, the other three known species of the genus are presumed to show more uniform distributions. Here we investigate the genetic diversity and population structure of R. minuchae and compare it with that found in two other species of the genus: R. anatolicus and R. quandratinodum. We conclude that although genetic diversity of R. minuchae is low, there is no evidence of a recent bottleneck, suggesting a gradual and natural fragmentation process. We also show extreme population differentiation at nuclear and mitochondrial markers, and isolation by distance at a local scale. Despite some evidence for inbreeding and low genetic variation within populations, we found almost no diploid males, a finding which contrasts with that expected in inbred Hymenoptera with single locus complementary sex determination. This could mean that sex is determined by another mechanism. We argue that continued low population size means that detrimental effects of inbreeding and low genetic variation are likely in the future. We suggest that a policy of artificial gene flow aimed at increasing within population variation is considered as a management option.

Statistical properties of population differentiation estimators under stepwise mutation in a finite island model.

Microsatellite loci mutate at an extremely high rate and are generally thought to evolve through a stepwise mutation model. Several differentiation statistics taking into account the particular mutation scheme of the microsatellite have been proposed. The most commonly used is R(ST) which is independent of the mutation rate under a generalized stepwise mutation model. F(ST) and R(ST) are commonly reported in the literature, but often differ widely. Here we compare their statistical performances using individual-based simulations of a finite island model. The simulations were run under different levels of gene flow, mutation rates, population number and sizes. In addition to the per locus statistical properties, we compare two ways of combining R(ST) over loci. Our simulations show that even under a strict stepwise mutation model, no statistic is best overall. All estimators suffer to different extents from large bias and variance. While R(ST) better reflects population differentiation in populations characterized by very low gene-exchange, F(ST) gives better estimates in cases of high levels of gene flow. The number of loci sampled (12, 24, or 96) has only a minor effect on the relative performance of the estimators under study. For all estimators there is a striking effect of the number of samples, with the differentiation estimates showing very odd distributions for two samples.

On the transition of genetic differentiation from isolation to panmixia: What we can learn from GST and D.

Population genetic differentiation characterizes the repartition of alleles among populations. It is commonly thought that genetic differentiation measures, such as GST and D, should be near zero when allele frequencies are close to their expected value in panmictic populations, and close to one when they are close to their expected value in isolated populations. To analyse those properties, we first derive analytically a reference function f of known parameters that describes how important features of genetic differentiation (e.g. gene diversity, proportion of private alleles, frequency of the most common allele) are close to their expected panmictic and isolation value. We find that the behaviour of function f differs according to three distinct mutation regimes defined by the scaled mutation rate and the number of populations. Then, we compare GST and D to f, and demonstrate that their signal of differentiation strongly depends on the mutation regime. In particular, we show that D captures well the variations of genetic diversity when mutation is weak, otherwise it overestimates it when panmixia is not met. GST detects population differentiation when mutation is intermediate but has a low sensitivity to the variations of genetic diversity when mutation is weak. When mutation is strong the domain of sensitivity of both measures are altered. Finally, we also point out the importance of the number of populations on genetic differentiation measures, and provide recommendations for the use of GST and D.

The population genetics of clonal and partially clonal diploids.

The consequences of variable rates of clonal reproduction on the population genetics of neutral markers are explored in diploid organisms within a subdivided population (island model). We use both analytical and stochastic simulation approaches. High rates of clonal reproduction will positively affect heterozygosity. As a consequence, nearly twice as many alleles per locus can be maintained and population differentiation estimated as F(ST) value is strongly decreased in purely clonal populations as compared to purely sexual ones. With increasing clonal reproduction, effective population size first slowly increases and then points toward extreme values when the reproductive system tends toward strict clonality. This reflects the fact that polymorphism is protected within individuals due to fixed heterozygosity. Contrarily, genotypic diversity smoothly decreases with increasing rates of clonal reproduction. Asexual populations thus maintain higher genetic diversity at each single locus but a lower number of different genotypes. Mixed clonal/sexual reproduction is nearly indistinguishable from strict sexual reproduction as long as the proportion of clonal reproduction is not strongly predominant for all quantities investigated, except for genotypic diversities (both at individual loci and over multiple loci).

Genetic differentiation in Silene dioica metapopulations: estimation of spatiotemporal effects in a successional plant species

Silene dioica is a diploid, dioecious, perennial, insect-pollinated herb and part of the deciduous phase of primary succession in Skeppsvik Archipelago, Gulf of Bothnia, Sweden. These islands are composed of material deposited and left underwater by melting ice at the end of the last ice age. A rapid and relatively constant rate of land uplift of 0.9 cm per year continually creates new islands available for colonization by plants. Because the higher deposits appear first, islands differ in age. Because it is possible to estimate the ages of islands and populations of plant species belonging to early stages of succession, the genetic dynamics occurring within an age-structured metapopulation can be investigated in this archipelago. Fifty-two island populations of S. dioica of known ages, sizes, and distances from each other were studied through electrophoretic data. A number of factors increase the degree of genetic differentiation among these island populations relative to an island model at equilibrium. Newly founded populations were more differentiated than those of intermediate age, which suggests that colonization dynamics increase genetic variance among populations. The very old populations, which decrease in size as they approach extinction, were more differentiated than intermediate-aged populations. Isolation by distance occurs in this system. Colonizers are likely to come from more than one source, and the migrant pool model best explains colonization events in the archipelago. Degree of environmental exposure also affects population differentiation.

Species-specific population structure in rock-specialized sympatric cichlid species in Lake Tanganyika, East Africa.

Species richness and geographical phenotypic variation in East African lacustrine cichlids are often correlated with ecological specializations and limited dispersal. This study compares mitochondrial and microsatellite genetic diversity and structure among three sympatric rock-dwelling cichlids of Lake Tanganyika, Eretmodus cyanostictus, Tropheus moorii, and Ophthalmotilapia ventralis. The species represent three endemic, phylogenetically distinct tribes (Eretmodini, Tropheini, and Ectodini), and display divergent ecomorphological and behavioral specialization. Sample locations span both continuous, rocky shoreline and a potential dispersal barrier in the form of a muddy bay. High genetic diversity and population differentiation were detected in T. moorii and E. cyanostictus, whereas much lower variation and structure were found in O. ventralis. In particular, while a 7-km-wide muddy bay curtails dispersal in all three species to a similar extent, gene flow along mostly continuous habitat appeared to be controlled by distance in E. cyanostictus, further restricted by site philopatry and/or minor habitat discontinuities in T. moorii, and unrestrained in O. ventralis. In contrast to the general pattern of high gene flow along continuous shorelines in rock-dwelling cichlids of Lake Malawi, our study identifies differences in population structure among stenotopic Lake Tanganyika species. The amount of genetic differentiation among populations was not related to the degree of geographical variation of body color, especially since more phenotypic variation is observed in O. ventralis than in the genetically highly structured E. cyanostictus.

Effects of selection and drift on G matrix evolution in a heterogeneous environment: a multivariate Qst-Fst Test with the freshwater snail Galba truncatula.

Unraveling the effect of selection vs. drift on the evolution of quantitative traits is commonly achieved by one of two methods. Either one contrasts population differentiation estimates for genetic markers and quantitative traits (the Q(st)-F(st) contrast) or multivariate methods are used to study the covariance between sets of traits. In particular, many studies have focused on the genetic variance-covariance matrix (the G matrix). However, both drift and selection can cause changes in G. To understand their joint effects, we recently combined the two methods into a single test (accompanying article by Martin et al.), which we apply here to a network of 16 natural populations of the freshwater snail Galba truncatula. Using this new neutrality test, extended to hierarchical population structures, we studied the multivariate equivalent of the Q(st)-F(st) contrast for several life-history traits of G. truncatula. We found strong evidence of selection acting on multivariate phenotypes. Selection was homogeneous among populations within each habitat and heterogeneous between habitats. We found that the G matrices were relatively stable within each habitat, with proportionality between the among-populations (D) and the within-populations (G) covariance matrices. The effect of habitat heterogeneity is to break this proportionality because of selection for habitat-dependent optima. Individual-based simulations mimicking our empirical system confirmed that these patterns are expected under the selective regime inferred. We show that homogenizing selection can mimic some effect of drift on the G matrix (G and D almost proportional), but that incorporating information from molecular markers (multivariate Q(st)-F(st)) allows disentangling the two effects.

Sex-biased dispersal in a migratory bat: a characterization using sex-specific demographic parameters.

We studied the noctule bat (Nyctalus noctula), in which the mitochondrial F(ST) is about 10 times that revealed by nuclear markers, to address two questions. We first verified whether random dispersal of one sex is compatible with highly contrasted mitochondrial and nuclear population structures. Using computer simulations, we then assessed the power of multilocus population differentiation tests when the expected population structure departs only slightly from panmixia. Using an island model with sex-specific demographic parameters, we found that random male dispersal is consistent with the population structure observed in the noctule. However, other parameter combinations are also compatible with the data. We computed the minimum sex bias in dispersal (at least 69% of the dispersing individuals are males), a result that would not be available if we had used more classical population genetic models. The power of multilocus population differentiation tests was unexpectedly high, the tests being significant in almost 100% of the replicates, although the observed population structure infered from nuclear markers was extremely low (F(ST) = 0.6%).

How accurate is the current picture of human genetic variation?

Our understanding of the distribution of worldwide human genomic diversity has greatly increased over recent years thanks to the availability of large data sets derived from short tandem repeats (STRs), insertion deletion polymorphisms (indels) and single nucleotide polymorphisms (SNPs). A concern, however, is that the current picture of worldwide human genomic diversity may be inaccurate because of biases in the selection process of genetic markers (so-called 'ascertainment bias'). To evaluate this problem, we first compared the distribution of genomic diversity between these three types of genetic markers in the populations from the HGDP-CEPH panel for evidence of bias or incongruities. In a second step, using a very relaxed set of criteria to prevent the intrusion of bias, we developed a new set of unbiased STR markers and compared the results against those from available panels. Contrarily to recent claims, our results show that the STR markers suffer from no discernible bias, and can thus be used as a baseline reference for human genetic diversity and population differentiation. The bias on SNPs is moderate compared to that on the set of indels analysed, which we recommend should be avoided for work describing the distribution of human genetic diversity or making inference on human settlement history.

Evolution in subdivided plant populations: concepts, recent advances and future directions.

SUMMARY: Research into the evolution of subdivided plant populations has long involved the study of phenotypic variation across plant geographic ranges and the genetic details underlying that variation. Genetic polymorphism at different marker loci has also allowed us to infer the long- and short-term histories of gene flow within and among populations, including range expansions and colonization-extinction dynamics. However, the advent of affordable genome-wide sequences for large numbers of individuals is opening up new possibilities for the study of subdivided populations. In this review, we consider what the new tools and technologies may allow us to do. In particular, we encourage researchers to look beyond the description of variation and to use genomic tools to address new hypotheses, or old ones afresh. Because subdivided plant populations are complex structures, we caution researchers away from adopting simplistic interpretations of their data, and to consider the patterns they observe in terms of the population genetic processes that have given rise to them; here, the genealogical framework of the coalescent will continue to be conceptually and analytically useful.

HIERFSTAT, a package for R to compute and test hierarchical F-statistics

The package HIERFSTAT for the statistical software R, created by the R Development Core Team, allows the estimate of hierarchical F-statistics from a hierarchy with any numbers of levels. In addition, it allows testing the statistical significance of population differentiation for these different levels, using a generalized likelihood-ratio test. The package HIERFSTAT is available at http://www.unil.ch/popgen/softwares/hierfstat.htm.

Divergence with gene flow and fine-scale phylogeographical structure in the wedge-billed woodcreeper, Glyphorynchus spirurus, a Neotropical rainforest bird.

Determining the relative roles of vicariance and selection in restricting gene flow between populations is of central importance to the evolutionary process of population divergence and speciation. Here we use molecular and morphological data to contrast the effect of isolation (by mountains and geographical distance) with that of ecological factors (altitudinal gradients) in promoting differentiation in the wedge-billed woodcreeper, Glyphorynchus spirurus, a tropical forest bird, in Ecuador. Tarsus length and beak size increased relative to body size with altitude on both sides of the Andes, and were correlated with the amount of moss on tree trunks, suggesting the role of selection in driving adaptive divergence. In contrast, molecular data revealed a considerable degree of admixture along these altitudinal gradients, suggesting that adaptive divergence in morphological traits has occurred in the presence of gene flow. As suggested by mitochondrial DNA sequence data, the Andes act as a barrier to gene flow between ancient subspecific lineages. Genome-wide amplified fragment length polymorphism markers reflected more recent patterns of gene flow and revealed fine-scale patterns of population differentiation that were not detectable with mitochondrial DNA, including the differentiation of isolated coastal populations west of the Andes. Our results support the predominant role of geographical isolation in driving genetic differentiation in G. spirurus, yet suggest the role of selection in driving parallel morphological divergence along ecological gradients.