Alps, genes, and chromosomes: their role in the formation of species in the Sorex araneus group (Mammalia, Insectivora), as inferred from two hybrid zones.

During the Pleistocene glaciations, the Alps were an efficient barrier to gene flow between isolated populations, often leading to allopatric speciation. Afterwards, the Alps strongly influenced the post-glacial recolonization of Europe and represent a major suture zone between differentiated populations. Two hybrid zones in the Swiss and French Alps between genetically and chromosomally well-differentiated species-the Valais shrew, Sorex antinorii, and the common shrew, S. araneus-were studied karyotypically and by analyzing the distribution of seven microsatellite loci. In the center of the Haslital hybrid zone the two species coexist over a distance of 900 m. Hybrid karyotypes, among them the most complex known in Sorex, are rare. F-statistics based on microsatellite data revealed a strong heterozygote deficit only in the center of the zone, due to the sympatric distribution of the two species with little hybridization between them. Structuring within the species (both F(IS) and F(ST)) was low. An hierarchical analysis showed a high level of interspecific differentiation. Results were compared with those previously reported in another hybrid zone located at Les Houches in the French Alps. Genetic structuring within and between species was comparable in both hybrid zones, although chromosomal incompatibilities are more important in Haslital, where a linkage block of the race-specific chromosomes should additionally impede gene flow. Evidence for a more restricted gene flow in Haslital comes from the genetically intermediate hybrid karyotypes, whereas in Les Houches, hybrid karyotypes are genetically identical to individuals of the pure karyotypic races. Genic and chromosomal introgression was observed in Les Houches, but not in Haslital. The possible influence of a river, separating the two species at Les Houches, on gene flow is discussed.

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.

Dispersal and Inbreeding Avoidance.

Using a game-theoretical approach, we investigate the dispersal patterns expected if inbreeding avoidance were the only reason for dispersal. The evolutionary outcome is always complete philopatry by one sex. The rate of dispersal by the other sex depends on patch size and mating system, as well as inbreeding and dispersal costs. If such costs are sex independent, then two stable equilibria coexist (male or female philopatry), with symmetric domains of attraction. Which sex disperses is determined entirely by history, genetic drift, and gene flow. An asymmetry in costs makes one domain of attraction extend at the expense of the other. In such a case, the dispersing sex might also be, paradoxically, the one that incurs the higher dispersal costs. As asymmetry increases, one equilibrium eventually disappears, which may result in a sudden evolutionary shift in the identity of the dispersing sex. Our results underline the necessity to control for phylogenetic relationships (e.g., through the use of independent-comparisons methods) when investigating empirical trends in dispersal. Our model also makes quantitative predictions on the rate of dispersal by the dispersing sex and suggests that inbreeding avoidance may only rarely be the sole reason for dispersal.

Local adaptation maintains clinal variation in melanin-based coloration of European barn owls (Tyto alba).

Ecological parameters vary in space, and the resulting heterogeneity of selective forces can drive adaptive population divergence. Clinal variation represents a classical model to study the interplay of gene flow and selection in the dynamics of this local adaptation process. Although geographic variation in phenotypic traits in discrete populations could be remainders of past adaptation, maintenance of adaptive clinal variation requires recurrent selection. Clinal variation in genetically determined traits is generally attributed to adaptation of different genotypes to local conditions along an environmental gradient, although it can as well arise from neutral processes. Here, we investigated whether selection accounts for the strong clinal variation observed in a highly heritable pheomelanin-based color trait in the European barn owl by comparing spatial differentiation of color and of neutral genes among populations. Barn owl's coloration varies continuously from white in southwestern Europe to reddish-brown in northeastern Europe. A very low differentiation at neutral genetic markers suggests that substantial gene flow occurs among populations. The persistence of pronounced color differentiation despite this strong gene flow is consistent with the hypothesis that selection is the primary force maintaining color variation among European populations. Therefore, the color cline is most likely the result of local adaptation.

Variable queen number in ant colonies: no impact on queen turnover, inbreeding, and population genetic differentiation in the ant Formica selysi.

Variation in queen number alters the genetic structure of social insect colonies, which in turn affects patterns of kin-selected conflict and cooperation. Theory suggests that shifts from single- to multiple-queen colonies are often associated with other changes in the breeding system, such as higher queen turnover, more local mating, and restricted dispersal. These changes may restrict gene flow between the two types of colonies and it has been suggested that this might ultimately lead to sympatric speciation. We performed a detailed microsatellite analysis of a large population of the ant Formica selysi, which revealed extensive variation in social structure, with 71 colonies headed by a single queen and 41 by multiple queens. This polymorphism in social structure appeared stable over time, since little change in the number of queens per colony was detected over a five-year period. Apart from queen number, single- and multiple-queen colonies had very similar breeding systems. Queen turnover was absent or very low in both types of colonies. Single- and multiple-queen colonies exhibited very small but significant levels of inbreeding, which indicates a slight deviation from random mating at a local scale and suggests that a small proportion of queens mate with related males. For both types of colonies, there was very little genetic structuring above the level of the nest, with no sign of isolation by distance. These similarities in the breeding systems were associated with a complete lack of genetic differentiation between single- and multiple-queen colonies, which provides no support for the hypothesis that change in queen number leads to restricted gene flow between social forms. Overall, this study suggests that the higher rates of queen turnover, local mating, and population structuring that are often associated with multiple-queen colonies do not appear when single- and multiple-queen colonies still coexist within the same population, but build up over time in populations consisting mostly of multiple-queen colonies.

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.

Mate availability and male dispersal in the Argentine ant Linepithema humile (Mayr) (=Iridomyrmex humilis)

In the Argentine ant Linepithema humile (=Iridomyrmex humilis) only males disperse whereas female sexuals (unmated winged queens) stay in their mother nest where they mate. This study investigated (1) whether dispersing males are accepted into foreign colonies, (2) whether they can mate with resident female sexuals, and (3) whether the propensity of males to disperse is affected by the expectation of mating in their mother nest. Field experiments demonstrated that males were accepted into foreign colonies only when these colonies contained female sexuals or queen pupae. Before and after the time of (sic) female sexuals, workers attacked and killed most of the foreign males. Laboratory experiments snowed that males that successfully enter foreign colonies can mate with resident female sexuals. The propensity of males to disperse was significantly influenced by the presence of female sexuals in their nest. Males were more likely to fly out from colonies containing no female sexuals than from those with them. These results are consistent with males preferentially dispersing when there is little or no opportunity to mate in their mother nest. Thus there are two mating strategies available for males: staying in their mother nest when an opportunity to mate arises or dispersing and attempting to mate in a foreign nest when there are no female sexuals in their mother nest. This latter behaviour could mediate gene flow between colonies and account for the lack of significant inbreeding previously documented in this species.

Genetic structure of the Lesser Gymnure (Genus Hylomys) in SE-Asia: evidence for two species

The Lesser Gymnure is a small galericine Insectivore living in the mainland forests of Southeast Asia, including the major islands of Sumatra, Java and Borneo. It is the only representative of the supposed monospecific genus Hylomys which is morphologically rather constant over its geographic range. Only marginal subspecific differentiation is currently recognized based in slight pelage colour variations. We report here the results of 35 gene loci revealed by protein electrophoresis on 23 specimens sampled over most of Southeast Asia. They were compared with outgroup, Erinaceaus europeaus, member of a distinct subfamily. The surveyed populations of Lesser Gymnures clearly group themselves into distinct taxa, one of which seems restricted to Sumatra, while the other occupies the whole geographic range. The genetic distance between these groups is two times greater than the divergence observed within groups: it is of the same order of magnitude as what is usually reported for congeneric mammal species, which supports their specific distinction. The lack of gene flow is also demonstrated by several diagnostic loci defining unambiguously each species. Both are only distantly related to the outgroup, a result which is consistant with their actual classification into two distinct subfamilies (Erinaceae and Galericinae). Concordant genetic and geographic subdivision of the widespread species further suggest that eustatic sea level changes during the Pleistocene produced predictable patterns in species differentiation.

Multilocus phylogeography of the common lizard Zootoca vivipara at the Ibero-Pyrenean suture zone reveals lowland barriers and high-elevation introgression.

BACKGROUND: The geographic distribution of evolutionary lineages and the patterns of gene flow upon secondary contact provide insight into the process of divergence and speciation. We explore the evolutionary history of the common lizard Zootoca vivipara (= Lacerta vivipara) in the Iberian Peninsula and test the role of the Pyrenees and the Cantabrian Mountains in restricting gene flow and driving lineage isolation and divergence. We also assess patterns of introgression among lineages upon secondary contact, and test for the role of high-elevation trans-mountain colonisations in explaining spatial patterns of genetic diversity. We use mtDNA sequence data and genome-wide AFLP loci to reconstruct phylogenetic relationships among lineages, and measure genetic structure RESULTS: The main genetic split in mtDNA corresponds generally to the French and Spanish sides of the Pyrenees as previously reported, in contrast to genome-wide AFLP data, which show a major division between NW Spain and the rest. Both types of markers support the existence of four distinct and geographically congruent genetic groups, which are consistent with major topographic barriers. Both datasets reveal the presence of three independent contact zones between lineages in the Pyrenean region, one in the Basque lowlands, one in the low-elevation mountains of the western Pyrenees, and one in the French side of the central Pyrenees. The latter shows genetic evidence of a recent, high-altitude trans-Pyrenean incursion from Spain into France. CONCLUSIONS: The distribution and age of major lineages is consistent with a Pleistocene origin and a role for both the Pyrenees and the Cantabrian Mountains in driving isolation and differentiation of Z. vivipara lineages at large geographic scales. However, mountain ranges are not always effective barriers to dispersal, and have not prevented a recent high-elevation trans-Pyrenean incursion that has led to asymmetrical introgression among divergent lineages. Cytonuclear discordance in patterns of genetic structure and introgression at contact zones suggests selection may be involved at various scales. Suture zones are important areas for the study of lineage formation and speciation, and our results show that biogeographic barriers can yield markedly different phylogeographic patterns in different vertebrate and invertebrate taxa.

Comparative population genetic structure in a plant-pollinator/seed predator system.

Comparative analyses of spatial genetic structure of populations of plants and the insects they interact with provide an indication of how gene flow, natural selection and genetic drift may jointly influence the distribution of genetic variation and potential for local co-adaptation for interacting species. Here, we analysed the spatial scale of genetic structure within and among nine populations of an interacting species pair, the white campion Silene latifolia and the moth Hadena bicruris, along a latitudinal gradient across Northern/Central Europe. This dioecious, short-lived perennial plant inhabits patchy, often disturbed environments. The moth H. bicruris acts both as its pollinator and specialist seed predator that reproduces by laying eggs in S. latifolia flowers. We used nine microsatellite markers for S. latifolia and eight newly developed markers for H. bicruris. We found high levels of inbreeding in most populations of both plant and pollinator/seed predator. Among populations, significant genetic structure was observed for S. latifolia but not for its pollinator/seed predator, suggesting that despite migration among populations of H. bicruris, pollen is not, or only rarely, carried over between populations, thus maintaining genetic structure among plant populations. There was a weak positive correlation between genetic distances of S. latifolia and H. bicruris. These results indicate that while significant structure of S. latifolia populations creates the potential for differentiation at traits relevant for the interaction with the pollinator/seed predator, substantial gene flow in H. bicruris may counteract this process in at least some populations.

Wheat alleles introgress into selfing wild relatives: empirical estimates from approximate Bayesian computation in Aegilops triuncialis.

Extensive gene flow between wheat (Triticum sp.) and several wild relatives of the genus Aegilops has recently been detected despite notoriously high levels of selfing in these species. Here, we assess and model the spread of wheat alleles into natural populations of the barbed goatgrass (Aegilops triuncialis), a wild wheat relative prevailing in the Mediterranean flora. Our sampling, based on an extensive survey of 31 Ae. triuncialis populations collected along a 60 km × 20 km area in southern Spain (Grazalema Mountain chain, Andalousia, totalling 458 specimens), is completed with 33 wheat cultivars representative of the European domesticated pool. All specimens were genotyped with amplified fragment length polymorphism with the aim of estimating wheat admixture levels in Ae. triuncialis populations. This survey first confirmed extensive hybridization and backcrossing of wheat into the wild species. We then used explicit modelling of populations and approximate Bayesian computation to estimate the selfing rate of Ae. triuncialis along with the magnitude, the tempo and the geographical distance over which wheat alleles introgress into Ae. triuncialis populations. These simulations confirmed that extensive introgression of wheat alleles (2.7 × 10(-4) wheat immigrants for each Ae. triuncialis resident, at each generation) into Ae. triuncialis occurs despite a high selfing rate (Fis ≈ 1 and selfing rate = 97%). These results are discussed in the light of risks associated with the release of genetically modified wheat cultivars in Mediterranean agrosystems.

Avian Evolutionary Genomics: Studies of Ficedula Flycatchers

In this thesis, different genetic tools are used to investigate both natural variation and speciation in the Ficedula flycatcher system: pied (Ficedula hypoleuca) and collared (F. albicollis) flycatchers. The molecular evolution of a gene involved in postnatal body growth, GH, has shown high degree of conservation at the mature protein between birds and mammals, whereas the variation observed in its signal peptide seems to be adaptive in pied flycatcher (I & II). Speciation is the process by which reproductive barriers to gene flow evolve between populations, and understanding the mechanisms involved in pre- and post-zygotic isolation have been investigated in Ficedula flycatchers. The Z chromosome have been suggested to be the hotspot for genes involved in speciation, thus sequencing of 13 Z-linked coding genes from the two species in allopatry and sympatry have been conducted (III). Surprisingly, the majority of Z-linked genes seemed to be highly conserved, suggesting instead a potential involvement of regulatory regions. Previous studies have shown that genes involved in hybrid fitness, female preferences and male plumage colouration are sex-linked. Hence, three pigmentation genes have been investigated: MC1R, AGRP, and TYRP1. Of these three genes, TYRP1 was identified as a strong candidate to be associated with black-brown plumage variation in sympatric populations, and hence is a strong candidate for a gene contributing to pre-zygotic isolation (IV). In sympatric areas, where pied and collared flycatchers have overlapping breeding areas, hybridization sometimes occurs leading to the production of unfit hybrids. By using a proteomic approach a novel expression pattern in hybrids was revealed compared to the parental species (V) and differentially expressed proteins subsequently identified by sequence similarity (VI). In conclusion, the Z chromosome appears to play an important role in flycatcher speciation, but probably not at the coding level. In addition the novel expression patterns might give new insights into the maladaptive hybrids.

Effects of colonization asymmetries on metapopulation persistence.

Ocean currents, prevailing winds, and the hierarchical structures of river networks are known to create asymmetries in re-colonization between habitat patches. The impacts of such asymmetries on metapopulation persistence are seldom considered, especially rarely in theoretical studies. Considering three classical models (the island, the stepping stone and the distance-dependent model), we explore how metapopulation persistence is affected by (i) asymmetry in dispersal strength, in which the colonization rate between two patches differs in direction, and (ii) asymmetry in connectivity, in which the overall colonization pattern displays asymmetry (circulating or dendritic networks). Viability can be drastically reduced when directional bias in dispersal strength is higher than 25%. Re-colonization patterns that allow for strong local connectivity provide the highest persistence compared to systems that allow circulation. Finally, asymmetry has relatively weak effects when metapopulations maintain strong general connectivity.

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.

Phylogeography and dispersal in the velvet gecko (Oedura lesueurii), and potential implications for conservation of an endangered snake (Hoplocephalus bungaroides).

ABSTRACT: BACKGROUND: To conserve critically endangered predators, we also need to conserve the prey species upon which they depend. Velvet geckos (Oedura lesueurii) are a primary prey for the endangered broad-headed snake (Hoplocephalus bungaroides), which is restricted to sandstone habitats in southeastern Australia. We sequenced the ND2 gene from 179 velvet geckos, to clarify the lizards' phylogeographic history and landscape genetics. We also analysed 260 records from a longterm (3-year) capture-mark-recapture program at three sites, to evaluate dispersal rates of geckos as a function of locality, sex and body size. RESULTS: The genetic analyses revealed three ancient lineages in the north, south and centre of the species' current range. Estimates of gene flow suggest low dispersal rates, constrained by the availability of contiguous rocky habitat. Mark-recapture records confirm that these lizards are highly sedentary, with most animals moving < 30 m from their original capture site even over multi-year periods. CONCLUSION: The low vagility of these lizards suggests that they will be slow to colonise vacant habitat patches; and hence, efforts to restore degraded habitats for broad-headed snakes may need to include translocation of lizards.