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.

Multivariate QST-FST comparisons: a neutrality test for the evolution of the g matrix in structured populations.

Neutrality tests in quantitative genetics provide a statistical framework for the detection of selection on polygenic traits in wild populations. However, the existing method based on comparisons of divergence at neutral markers and quantitative traits (Q(st)-F(st)) suffers from several limitations that hinder a clear interpretation of the results with typical empirical designs. In this article, we propose a multivariate extension of this neutrality test based on empirical estimates of the among-populations (D) and within-populations (G) covariance matrices by MANOVA. A simple pattern is expected under neutrality: D = 2F(st)/(1 - F(st))G, so that neutrality implies both proportionality of the two matrices and a specific value of the proportionality coefficient. This pattern is tested using Flury's framework for matrix comparison [common principal-component (CPC) analysis], a well-known tool in G matrix evolution studies. We show the importance of using a Bartlett adjustment of the test for the small sample sizes typically found in empirical studies. We propose a dual test: (i) that the proportionality coefficient is not different from its neutral expectation [2F(st)/(1 - F(st))] and (ii) that the MANOVA estimates of mean square matrices between and among populations are proportional. These two tests combined provide a more stringent test for neutrality than the classic Q(st)-F(st) comparison and avoid several statistical problems. Extensive simulations of realistic empirical designs suggest that these tests correctly detect the expected pattern under neutrality and have enough power to efficiently detect mild to strong selection (homogeneous, heterogeneous, or mixed) when it is occurring on a set of traits. This method also provides a rigorous and quantitative framework for disentangling the effects of different selection regimes and of drift on the evolution of the G matrix. We discuss practical requirements for the proper application of our test in empirical studies and potential extensions.

Genetic differentiation in two European tree frog (Hyla arborea) metapopulations in contrasted landscapes of Western Switzerland

The survival of threatened species as the European tree frog (Hyla arborea) is strongly dependent on the genetic variability within populations, as well as gene flow between them. In Switzerland, only two sectors in its western part still harbour metapopulations. The first is characterised by a very heterogeneous and urbanized landscape, while the second is characterised by a uninterrupted array of suitable habitats. In this study, six microsatellite loci were used to establish levels of genetic differentiation among the populations from the two different locations. The results show that the metapopulations have: (i) weak levels of genetic differentiation (FST within metapopulation ≈ 0.04), (ii) no difference in levels of genetic structuring between them, (iii) significant (p = 0.019) differences in terms of genetic diversity (Hs) and observed heterozygozity (Ho), the metapopulation located in a disturbed landscape showing lower values. Our results suggest that even if the dispersal of H. arborea among contiguous ponds seems to be efficient in areas of heterogeneous landscape, a loss of genetic diversity can occur.

Parental influences on pathogen resistance in brown trout embryos and effects of outcrossing within a river network.

Phenotypic plasticity can increase tolerance to heterogeneous environments but the elevations and slopes of reaction norms are often population specific. Disruption of locally adapted reaction norms through outcrossing can lower individual viability. Here, we sampled five genetically distinct populations of brown trout (Salmo trutta) from within a river network, crossed them in a full-factorial design, and challenged the embryos with the opportunistic pathogen Pseudomonas fluorescens. By virtue of our design, we were able to disentangle effects of genetic crossing distance from sire and dam effects on early life-history traits. While pathogen infection did not increase mortality, it was associated with delayed hatching of smaller larvae with reduced yolk sac reserves. We found no evidence of a relationship between genetic distance (W, FST) and the expression of early-life history traits. Moreover, hybrids did not differ in phenotypic means or reaction norms in comparison to offspring from within-population crosses. Heritable variation in early life-history traits was found to remain stable across the control and pathogen environments. Our findings show that outcrossing within a rather narrow geographical scale can have neutral effects on F1 hybrid viability at the embryonic stage, i.e. at a stage when environmental and genetic effects on phenotypes are usually large.

Genetic consequences of the ice ages on nurseries of the bat Myotis myotis: a mitochondrial and nuclear survey.

Analyses of mitochondrial DNA (mtDNA) control region polymorphism and of variation at 10 nuclear microsatellite loci were used to investigate the mechanisms and genetic consequences of postglacial expansion of Myotis myotis in Europe. Initial sampling consisted of 480 bats genotyped in 24 nursery colonies arranged along a transect of approximately 3000 km. The phylogeographical survey based on mtDNA sequences revealed the existence of major genetic subdivisions across this area, with several suture zones between haplogroups. Such zones of secondary contact were found in the Alps and Rhodopes, whereas other potential barriers to gene flow, like the Pyrenees, did not coincide with genetic discontinuities. Areas of population admixture increased locally the genetic diversity of colonies, which confounded the northward decrease in nucleotide diversity predicted using classical models of postglacial range expansion. However, when analyses were restricted to a subset of 15 nurseries originating from a single presumed glacial refugium, mtDNA polymorphism did indeed support a northwards decrease in diversity. Populations were also highly structured (PhiST = 0.384). Conversely, the same subset of colonies showed no significant latitudinal decrease in microsatellite diversity and much less population structure (FST = 0.010), but pairwise genetic differentiation at these nuclear markers was strongly correlated with increasing geographical distance. Together, this evidence suggests that alleles carried via male bats have maintained enough nuclear gene flow to counteract the effects of recurrent bottlenecks generally associated with recolonization processes. As females are highly philopatric, we argue that the maternally transmitted mtDNA marker better reflects the situation of past, historical gene flow, whereas current levels of gene flow are better reflected by microsatellite markers.

Strong reproductive barriers in a narrow hybrid zone of West-Mediterranean green toads (Bufo viridis subgroup) with Plio-Pleistocene divergence

Background One key question in evolutionary biology deals with the mode and rate at which reproductive isolation accumulates during allopatric speciation. Little is known about secondary contacts of recently diverged anuran species. Here we conduct a multi-locus field study to investigate a contact zone between two lineages of green toads with an estimated divergence time of 2.7 My, and report results from preliminary experimental crosses. Results The Sicilian endemic Bufo siculus and the Italian mainland-origin B. balearicus form a narrow hybrid zone east of Mt. Etna. Despite bidirectional mtDNA introgression over a ca. 40 km North-South cline, no F1 hybrids could be found, and nuclear genomes display almost no admixture. Populations from each side of the contact zone showed depressed genetic diversity and very strong differentiation (FST = 0.52). Preliminary experimental crosses point to a slightly reduced fitness in F1 hybrids, a strong hybrid breakdown in backcrossed offspring (F1 x parental, with very few reaching metamorphosis) and a complete and early mortality in F2 (F1 x F1). Conclusion Genetic patterns at the contact zone are molded by drift and selection. Local effective sizes are reduced by the geography and history of the contact zone, B. balearicus populations being at the front wave of a recent expansion (late Pleistocene). Selection against hybrids likely results from intrinsic genomic causes (disruption of coadapted sets of genes in backcrosses and F2-hybrids), possibly reinforced by local adaptation (the ranges of the two taxa roughly coincide with the borders of semiarid and arid climates). The absence of F1 in the field might be due to premating isolation mechanisms. Our results, show that these lineages have evolved almost complete reproductive isolation after some 2.7 My of divergence, contrasting sharply with evidence from laboratory experiments that some anuran species may still produce viable F1 offspring after > 20 My of divergence.

Geographic variation in sex-chromosome differentiation in the common frog (Rana temporaria).

In sharp contrast with birds and mammals, sex-determination systems in ectothermic vertebrates are often highly dynamic and sometimes multifactorial. Both environmental and genetic effects have been documented in common frogs (Rana temporaria). One genetic linkage group, mapping to the largest pair of chromosomes and harbouring the candidate sex-determining gene Dmrt1, associates with sex in several populations throughout Europe, but association varies both within and among populations. Here, we show that sex association at this linkage group differs among populations along a 1500-km transect across Sweden. Genetic differentiation between sexes is strongest (FST  = 0.152) in a northern-boreal population, where male-specific alleles and heterozygote excesses (FIS  = -0.418 in males, +0.025 in females) testify to a male-heterogametic system and lack of X-Y recombination. In the southernmost population (nemoral climate), in contrast, sexes share the same alleles at the same frequencies (FST  = 0.007 between sexes), suggesting unrestricted recombination. Other populations show intermediate levels of sex differentiation, with males falling in two categories: some cluster with females, while others display male-specific Y haplotypes. This polymorphism may result from differences between populations in the patterns of X-Y recombination, co-option of an alternative sex-chromosome pair, or a mixed sex-determination system where maleness is controlled either by genes or by environment depending on populations or families. We propose approaches to test among these alternative models, to disentangle the effects of climate and phylogeography on the latitudinal trend, and to sort out how this polymorphism relates to the 'sexual races' described in common frogs in the 1930s.

Microsatellite variation reveals low genetic subdivision in a chromosome race of Sorex araneus (Mammalia, Insectivora)

Two hundred and forty-five individuals of the common shrew (Sorex araneus, Insectivora, Mammalia) from 24 sampling localities situated in four different valleys of the western European Alps were genotyped for six microsatellite loci. Allelic variability ranged from 3 to 32 different alleles at a single locus and the average gene diversity over all loci was 0.69. An analysis for F and R statistics revealed weak genetic population subdivision (Fst = 0.032; Rst = 0.016). This suggests considerable gene flow and little phylogeographic structure within and between valleys. We tested whether a stepwise mutation model (SMM) better explained variation at the microsatellite loci than an infinite allele model (IAM). No trend in favor of either model was detected.

The effect of innovation and sex-specific migration on neutral cultural differentiation

Studies of behaviour are increasingly focusing on acquisition of traits through cultural inheritance. Comparison of patterns of spatial population structure (FST) between neutral genetic loci and behavioural or cultural traits can been used to test hypotheses about demography, life history, and the mechanisms of inheritance/transmission of these traits in humans, chimpanzees and other animals. Here, we develop analytical expectations to show how FST in cultural traits can differ strongly from that measured at neutral genetic markers if migration is largely restricted to one sex but social learning is predominantly modelled on the other (e.g. males migrate, females serve as models for cultural traits), if one individual is the learning model for many, or if rates of innovation (individual learning) are high or rates of social learning are low. We discuss how comparisons of FST between genetic loci and behavioural traits can be applied to evaluate the importance of innovation in shaping patterns of cultural differentiation, as even low rates of innovation can considerably reduce FST, relative to observed structure at neutral genetic loci. Our results also suggest that differentiation in neutral cultural traits should occur over much smaller scales in species with male migration and female enculturation (or the reverse).

High quantitative and no molecular differentiation of a freshwater snail (Galba truncatula) between temporary and permanent water habitats.

We investigate the variation in quantitative and molecular traits in the freshwater snail Galba truncatula, from permanent and temporary water habitats. Using a common garden experiment, we measured 20 quantitative traits and molecular variation using seven microsatellites in 17 populations belonging to these two habitats. We estimated trait means in each habitat. We also estimated the distributions of overall genetic quantitative variation (QST), and of molecular variation (FST), within and between habitats. Overall, we observed a lack of association between molecular and quantitative variance. Among habitats, we found QST>FST, an indication of selection for different optima. Individuals from temporary water habitat matured older, at a larger size and were less fecund than individuals from permanent water habitat. We discuss these findings in the light of several theories for life-history traits evolution.

Testing the effects of genetic crossing distance on embryo survival within a metapopulation of brown trout (Salmo trutta)

Predicting progeny performance from parental genetic divergence can potentially enhance the efficiency of supportive breeding programmes and facilitate risk assessment. Yet, experimental testing of the effects of breeding distance on offspring performance remains rare, especially in wild populations of vertebrates. Recent studies have demonstrated that embryos of salmonid fish are sensitive indicators of additive genetic variance for viability traits. We therefore used gametes of wild brown trout (Salmo trutta) from five genetically distinct populations of a river catchment in Switzerland, and used a full factorial design to produce over 2,000 embryos in 100 different crosses with varying genetic distances (FST range 0.005-0.035). Customized egg capsules allowed recording the survival of individual embryos until hatching under natural field conditions. Our breeding design enabled us to evaluate the role of the environment, of genetic and nongenetic parental contributions, and of interactions between these factors, on embryo viability. We found that embryo survival was strongly affected by maternal environmental (i.e. non-genetic) effects and by the microenvironment, i.e. by the location within the gravel. However, embryo survival was not predicted by population divergence, parental allelic dissimilarity, or heterozygosity, neither in the field nor under laboratory conditions. Our findings suggest that the genetic effects of inter-population hybridization within a genetically differentiated meta-population can be minor in comparison to environmental effects.

Heritability and quantitative genetic divergence of serotiny, a fire-persistence plant trait.

BACKGROUND AND AIMS: Although it is well known that fire acts as a selective pressure shaping plant phenotypes, there are no quantitative estimates of the heritability of any trait related to plant persistence under recurrent fires, such as serotiny. In this study, the heritability of serotiny in Pinus halepensis is calculated, and an evaluation is made as to whether fire has left a selection signature on the level of serotiny among populations by comparing the genetic divergence of serotiny with the expected divergence of neutral molecular markers (QST-FST comparison). METHODS: A common garden of P. halepensis was used, located in inland Spain and composed of 145 open-pollinated families from 29 provenances covering the entire natural range of P. halepensis in the Iberian Peninsula and Balearic Islands. Narrow-sense heritability (h(2)) and quantitative genetic differentiation among populations for serotiny (QST) were estimated by means of an 'animal model' fitted by Bayesian inference. In order to determine whether genetic differentiation for serotiny is the result of differential natural selection, QST estimates for serotiny were compared with FST estimates obtained from allozyme data. Finally, a test was made of whether levels of serotiny in the different provenances were related to different fire regimes, using summer rainfall as a proxy for fire regime in each provenance. KEY RESULTS: Serotiny showed a significant narrow-sense heritability (h(2)) of 0·20 (credible interval 0·09-0·40). Quantitative genetic differentiation among provenances for serotiny (QST = 0·44) was significantly higher than expected under a neutral process (FST = 0·12), suggesting adaptive differentiation. A significant negative relationship was found between the serotiny level of trees in the common garden and summer rainfall of their provenance sites. CONCLUSIONS: Serotiny is a heritable trait in P. halepensis, and selection acts on it, giving rise to contrasting serotiny levels among populations depending on the fire regime, and supporting the role of fire in generating genetic divergence for adaptive traits.

Testing differentiation in diploid populations.

We examine the power of different exact tests of differentiation for diploid populations. Since there is not necessarily random mating within populations, the appropriate hypothesis to construct exact tests is that of independent sampling of genotypes. There are two categories of tests, FST-estimator tests and goodness of fit tests. In this latter category, we distinguish "allelic statistics", which account for the nature of alleles within genotypes, from "genotypic statistics" that do not. We show that the power of FST-estimator tests and of allelic goodness of fit tests are similar when sampling is balanced, and higher than the power of genotypic goodness of fit tests. When sampling is unbalanced, the most powerful tests are shown to belong to the allelic goodness of fit group.

Local adaptation and life history traits evolution in heterogeneous environment : an empirical example in the freshwater snail, Galba truncatula

Abstract: To understand the processes of evolution, biologists are interested in the ability of a population to respond to natural or artificial selection. The amount of genetic variation is often viewed as the main factor allowing a species to answer to selection. Many theories have thus focused on the maintenance of genetic variability. Ecologists and population geneticists have long-suspected that the structure of the environment is connected to the maintenance of diversity. Theorists have shown that diversity can be permanently and stably maintained in temporal and spatial varying environment in certain conditions. Moreover, varying environments have been also theoretically demonstrated to cause the evolution of divergent life history strategies in the different niches constituting the environment. Although there is a huge number of theoretical studies selection and on life history evolution in heterogeneous environments, there is a clear lack of empirical studies. The purpose of this thesis was to. empirically study the evolutionary consequences of a heterogeneous environment in a freshwater snail Galba truncatula. Indeed, G. truncatula lives in two habitat types according the water availability. First, it can be found in streams or ponds which never completely dry out: a permanent habitat. Second, G. truncatula can be found in pools that freeze during winter and dry during summer: a temporary habitat. Using a common garden approach, we empirically demonstrated local adaptation of G. truncatula to temporary and permanent habitats. We used at first a comparison of molecular (FST) vs. quantitative (QST) genetic differentiation between temporary and permanent habitats. To confirm the pattern QST> FST between habitats suggesting local adaptation, we then tested the desiccation resistance of individuals from temporary and permanent habitats. This study confirmed that drought resistance seemed to be the main factor selected between habitats, and life history traits linked to the desiccation resistance were thus found divergent between habitats. However, despite this evidence of selection acting on mean values of traits between habitats, drift was suggested to be the main factor responsible of variation in variances-covariances between populations. At last, we found life history traits variation of individuals in a heterogeneous environment varying in parasite prevalence. This thesis empirically demonstrated the importance of heterogeneous environments in local adaptation and life history evolution and suggested that more experimental studies are needed to investigate this topic. Résumé: Les biologistes se sont depuis toujours intéressés en l'aptitude d'une population à répondre à la sélection naturelle. Cette réponse dépend de la quantité de variabilité génétique présente dans cette population. Plus particulièrement, les théoriciens se sont penchés sur la question du maintient de la variabilité génétique au sein d'environnements hétérogènes. Ils ont alors démontré que, sous certaines conditions, la diversité génétique peut se maintenir de manière stable et permanente dans des environnements variant au niveau spatial et temporel. De plus, ces environments variables ont été démontrés comme responsable de divergence de traits d'histoire de vie au sein des différentes niches constituant l'environnement. Cependant, malgré ce nombre important d'études théoriques portant sur la sélection et l'évolution des traits d'histoire de vie en environnement hétérogène, les études empiriques sont plus rares. Le but de cette thèse était donc d'étudier les conséquences évolutives d'un environnement hétérogène chez un esgarcot d'eau douce Galba truncatula. En effet, G. truncatula est trouvé dans deux types d'habitats qui diffèrent par leur niveau d'eau. Le premier, l'habitat temporaire, est constitué de flaques d'eau qui peuvent s'assécher pendant l'été et geler pendant l'hiver. Le second, l'habitat permanent, correspond à des marres ou à des ruisseaux qui ont un niveau d'eau constant durant toute l'année. Utilisant une approche expérimentale de type "jardin commun", nous avons démontré l'adaptation locale des individus à leur type d'habitat, permanent ou temporaire. Nous avons utilisé l'approche Fsr/QsT qui compare la différentiation génétique moléculaire avec la différentiation génétique quantitative entre les 2 habitats. Le phénomène d'adapation locale démontré par QsT > FsT, a été testé experimentalement en mesurant la résistance à la dessiccation d'individus d'habitat temporaire et permanent. Cette étude confirma que la résistance à la sécheresse a été sélectionné entre habitats et que les traits responsables de cette resistance sont différents entre habitats. Cependant si la sélection agit sur la valeur moyenne des traits entre habitats, la dérive génétique semble être le responsable majeur de la différence de variances-covariances entre populations. Pour finir, une variation de traits d'histoire de vie a été trouvée au sein d'un environnement hétérogène constitué de populations variants au niveau de leur taux de parasitisme. Pour conclure, cette thèse a donc démontré l'importance d'un environnement hétérogène sur l'adaptation locale et l'évolution des traits d'histoire de vie et suggère que plus d'études empiriques sur le sujet sont nécessaires.

Evolutionary aspects of population structure for molecular and quantitative traits in the freshwater snail Radix balthica.

Detecting the action of selection in natural populations can be achieved using the QST-FST comparison that relies on the estimation of FST with neutral markers, and QST using quantitative traits potentially under selection. QST higher than FST suggests the action of directional selection and thus potential local adaptation. In this article, we apply the QST-FST comparison to four populations of the hermaphroditic freshwater snail Radix balthica located in a floodplain habitat. In contrast to most studies published so far, we did not detect evidence of directional selection for local optima for any of the traits we measured: QST calculated using three different methods was never higher than FST. A strong inbreeding depression was also detected, indicating that outcrossing is probably predominant over selfing in the studied populations. Our results suggest that in this floodplain habitat, local adaptation of R. balthica populations may be hindered by genetic drift, and possibly altered by uneven gene flow linked to flood frequency.