Exact moment calculations for genetic models with migration, mutation, and drift

Using properties of moment stationarity we develop exact expressions for the mean and covariance of allele frequencies at a single locus for a set of populations subject to drift, mutation, and migration. Some general results can be obtained even for arbitrary mutation and migration matrices, for example: (1) Under quite general conditions, the mean vector depends only on mutation rates, not on migration rates or the number of populations. (2) Allele frequencies covary among all pairs of populations connected by migration. As a result, the drift, mutation, migration process is not ergodic when any finite number of populations is exchanging genes. in addition, we provide closed form expressions for the mean and covariance of allele frequencies in Wright's finite-island model of migration under several simple models of mutation, and we show that the correlation in allele frequencies among populations can be very large for realistic rates of mutation unless an enormous number of populations are exchanging genes. As a result, the traditional diffusion approximation provides a poor approximation of the stationary distribution of allele frequencies among populations. Finally, we discuss some implications of our results for measures of population structure based on Wright's F-statistics.

Stay or drift? Queen acceptance in the ant Formica paralugubris

The acceptance of new queens in ant colonies has profound effects on colony kin structure and inclusive fitness of workers. Therefore, it is important to study the recognition and discrimination behaviour of workers towards reproductive individuals entering established colonies. We examined the acceptance rate of queens in populations of the highly polygynous ant F. paralugubris, where the genetic differentiation among nests and discrimination ability among workers suggest that workers might reject foreign queens. We experimentally introduced young queens in their natal nest and in foreign nests. Surprisingly, the survival rate of mated queens did not differ significantly when introduced in a foreign male-producing nest, a foreign female-producing nest, or the natal nest. Moreover, the survival of virgin queens in their natal nest was twice the one of mated queens, suggesting that mating status plays an important role for acceptance. The results indicate that other factors than queen discrimination by workers are implicated in the limited long-distance gene flow between nests in these populations.

The effect of host plant and isolation on the genetic structure of phytophagous insects: A preliminary study on a bruchid beetle

Genetic differentiation is a consequence of the combination of drift and restriction in gene flow between populations due to barriers to dispersal, or selection against individuals resulting from inter-population matings In phytophagous insects, local adaptation to different kinds of host plants can sometimes lead to reproductive isolation and thus to genetic structuring, or even to speciation Acanthoscelides. obtectus Say is a bean bruchid specialized on beans of the Phaseolus vulgaris group, attacking both wild and domesticated forms of P vulgaris., and P coccineus This study reveals that the genetic structure of populations of this bruchid is explained mainly by their geographical location and is not related to a particular kind (wild or domesticated) of bean In contrast, the species of bean might have led, to some extent, to genetic structuring in these bruchids, although our sampling is too limited to address such process unambiguously. If confirmed, it would corroborate preliminary results found for the parasitoid species that attack Acanthoscelides species, which might show a genetic structure depending on the species of host plant

quantiNemo: an individual-based program to simulate quantitative traits with explicit genetic architecture in a dynamic metapopulation.

quantiNemo is an individual-based, genetically explicit stochastic simulation program. It was developed to investigate the effects of selection, mutation, recombination and drift on quantitative traits with varying architectures in structured populations connected by migration and located in a heterogeneous habitat. quantiNemo is highly flexible at various levels: population, selection, trait(s) architecture, genetic map for QTL and/or markers, environment, demography, mating system, etc. quantiNemo is coded in C++ using an object-oriented approach and runs on any computer platform. Availability: Executables for several platforms, user's manual, and source code are freely available under the GNU General Public License at http://www2.unil.ch/popgen/softwares/quantinemo.

Linking genetic and species diversities in a community of freshwater gastropods

Summary Biodiversity is usually studied through species or genetic diversities. To date, these two levels of diversity have remained the independent .fields of investigations of community ecologists and population geneticists. However, recent joint analyses of species and genetic diversities have suggested that common processes may underlie the two levels. Positive correlations between species diversity and genetic diversity may arise when the effects of drift and migration overwhelm selective effects. The first goal of this thesis was to make a joint investigation of the patterns of species and genetic diversity in a community of freshwater gastropods living in a floodplain habitat. The second goal was to determine, as far as possible, the relative influences of the processes underlying the patterns observed at each level. In chapter 2 we investigate the relative influences of different evolutionary forces in shaping the genetic structure of Radix balthica populations. Results revealed that the structure inferred using quantitative traits was lower or equal to the one inferred using neutral molecular markers. Consequently, the pattern of structure observed could be only due to random drift, possibly to uniform selection, but definitely not to selection for local optima. In chapter 3, we analyze the temporal variation of species and genetic diversities in five localities. An extended period of drought occurred at the end of the study period leading to decay of both species and genetic diversities. This parallel loss of diversity following a natural perturbation highlighted the role sometimes predominant of random drift over selection on patterns of biodiversity in a floodplain habitat. In chapter 4, we compare the spatial genetic structures of two sympatric species: Radix balthica and Planorbis carinatus. We found that R. balthica populations are weakly structured and have moderate to high values of gene diversity. In contrast, P. carinatus populations are highly structured and poorly diverse. Then we measured correlations between various indices of species and genetic diversity using genetic data .from the two species. We found only one significant correlation: between species richness and gene diversity of P. carinatus. This result highlights the .need to use genetic date from more than one species to infer correlations between species and genetic diversities. Overall, this thesis provided new insights into the common processes underlying patterns of species and genetic diversity. Résumé La biodiversité est généralement étudiée au niveau de la diversité génétique ou spécifique. Ces deux niveaux sont restés jusqu'à maintenant les domaines d'investigation séparés des généticiens des populations et des écologistes des communautés. Cependant, des analyses conjointes des diversités génétique et spécifique ont récemment suggéré que des processus similaires pouvaient influencer ces deux niveaux. Des corrélations positives entre les diversités génétique et spécifique pourraient être dues aux effets de migration et de dérive qui dominent les effets sélectifs. Le premier but de cette thèse était de faire une étude conjointe des diversités génétique et spécifique dans une communauté de gastéropodes d'eau douce. Le second objectif était de déterminer les influences relatives des différents processus liés à chaque niveau de diversité. Dans le chapitre 2 nous cherchons à déterminer quelles forces évolutives influencent la structure génétique de quatre populations de Radix balthica. La structure mesurée sur des traits quantitatifs s'est révélée être plus faible ou égale à celle mesurée avec des marqueurs moléculaires neutres. La structure observée pourrait ainsi être due uniquement à la dérive génétique, potentiellement à la sélection uniforme, mais en aucun cas à la sélection locale pour différents optima. Dans le chapitre 3 nous analysons la variation temporelle des diversités génétique et spécifique dans cinq localités. Une récente période de sécheresse a causé une diminution parallèle des deux niveaux de diversité. Cette perturbation à mis en évidence le rôle parfois prépondérant de la dérive par rapport à celui de la sélection dans le déterminisme de la biodiversité dans un écosytème alluvial. Dans le chapitre 4, nous comparons la structure génétique spatiale de deux espèces vivant en sympatrie : Radix balthica et Planorbis carinatus. Les populations de R. balthica sont peu structurées et présentent un niveau de diversité relativement élevé alors que celles de P. carinatus sont fortement structurées et peu diversifiées. Nous avons ensuite mesuré différentes corrélations entre les diversités génétique et spécifique, mais la seule relation significative a été trouvée entre la richesse spécifique et la diversité génétique de P. carinatus. Ainsi, cette thèse a permis de découvrir de nouveaux aspects des processus qui influencent en parallèle la diversité aux niveaux génétique et spécifique.

Genetic analysis of complex demographic scenarios: Spatially expanding populations of the cane toad, Bufo marinus

Inferring the spatial expansion dynamics of invading species from molecular data is notoriously difficult due to the complexity of the processes involved. For these demographic scenarios, genetic data obtained from highly variable markers may be profitably combined with specific sampling schemes and information from other sources using a Bayesian approach. The geographic range of the introduced toad Bufo marinus is still expanding in eastern and northern Australia, in each case from isolates established around 1960. A large amount of demographic and historical information is available on both expansion areas. In each area, samples were collected along a transect representing populations of different ages and genotyped at 10 microsatellite loci. Five demographic models of expansion, differing in the dispersal pattern for migrants and founders and in the number of founders, were considered. Because the demographic history is complex, we used an approximate Bayesian method, based on a rejection-regression algorithm. to formally test the relative likelihoods of the five models of expansion and to infer demographic parameters. A stepwise migration-foundation model with founder events was statistically better supported than other four models in both expansion areas. Posterior distributions supported different dynamics of expansion in the studied areas. Populations in the eastern expansion area have a lower stable effective population size and have been founded by a smaller number of individuals than those in the northern expansion area. Once demographically stabilized, populations exchange a substantial number of effective migrants per generation in both expansion areas, and such exchanges are larger in northern than in eastern Australia. The effective number of migrants appears to be considerably lower than that of founders in both expansion areas. We found our inferences to be relatively robust to various assumptions on marker. demographic, and historical features. The method presented here is the only robust, model-based method available so far, which allows inferring complex population dynamics over a short time scale. It also provides the basis for investigating the interplay between population dynamics, drift, and selection in invasive species.

Genetic structure in a tropical lek-breeding bird, the blue manakin (Chiroxiphia caudata) in the Brazilian Atlantic Forest

Determining the genetic structure of tropical bird populations is important for assessing potential genetic effects of future habitat fragmentation and for testing hypotheses about evolutionary mechanisms promoting diversification. Here we used 10 microsatellite DNA loci to describe levels of genetic differentiation for five populations of the lek-mating blue manakin (Chiroxiphia caudata), sampled along a 414-km transect within the largest remaining continuous tract of the highly endangered Atlantic Forest habitat in southeast Brazil. We found small but significant levels of differentiation between most populations. F-ST values varied from 0.0 to 0.023 (overall F-ST = 0.012) that conformed to a strong isolation by distance relationship, suggesting that observed levels of differentiation are a result of migration-drift equilibrium. N(e)m values estimated using a coalescent-based method were small (<= 2 migrants per generation) and close to the minimum level required to maintain genetic similarity between populations. An implication of these results is that if future habitat fragmentation reduces dispersal between populations to even a small extent, then individual populations may undergo a loss of genetic diversity due to an increase in the relative importance of drift, since inbreeding effective population sizes are relatively small (N-e similar to 1000). Our findings also demonstrate that population structuring can occur in a tropical bird in continuous habitat in the absence of geographical barriers possibly due to behavioural features of the species.

Contemporary and historic factors influence differently genetic differentiation and diversity in a tropical palm

Population genetics theory predicts loss in genetic variability because of drift and inbreeding in isolated plant populations; however, it has been argued that long-distance pollination and seed dispersal may be able to maintain gene flow, even in highly fragmented landscapes. We tested how historical effective population size, historical migration and contemporary landscape structure, such as forest cover, patch isolation and matrix resistance, affect genetic variability and differentiation of seedlings in a tropical palm (Euterpe edulis) in a human-modified rainforest. We sampled 16 sites within five landscapes in the Brazilian Atlantic forest and assessed genetic variability and differentiation using eight microsatellite loci. Using a model selection approach, none of the covariates explained the variation observed in inbreeding coefficients among populations. The variation in genetic diversity among sites was best explained by historical effective population size. Allelic richness was best explained by historical effective population size and matrix resistance, whereas genetic differentiation was explained by matrix resistance. Coalescence analysis revealed high historical migration between sites within landscapes and constant historical population sizes, showing that the genetic differentiation is most likely due to recent changes caused by habitat loss and fragmentation. Overall, recent landscape changes have a greater influence on among-population genetic variation than historical gene flow process. As immediate restoration actions in landscapes with low forest amount, the development of more permeable matrices to allow the movement of pollinators and seed dispersers may be an effective strategy to maintain microevolutionary processes.

Population amalgamation and genetic variation: observations on artificially agglomerated tribal populations of Central and South America.

The interpretation of data on genetic variation with regard to the relative roles of different evolutionary factors that produce and maintain genetic variation depends critically on our assumptions concerning effective population size and the level of migration between neighboring populations. In humans, recent population growth and movements of specific ethnic groups across wide geographic areas mean that any theory based on assumptions of constant population size and absence of substructure is generally untenable. We examine the effects of population subdivision on the pattern of protein genetic variation in a total sample drawn from an artificial agglomerate of 12 tribal populations of Central and South America, analyzing the pooled sample as though it were a single population. Several striking findings emerge. (1) Mean heterozygosity is not sensitive to agglomeration, but the number of different alleles (allele count) is inflated, relative to neutral mutation/drift/equilibrium expectation. (2) The inflation is most serious for rare alleles, especially those which originally occurred as tribally restricted "private" polymorphisms. (3) The degree of inflation is an increasing function of both the number of populations encompassed by the sample and of the genetic divergence among them. (4) Treating an agglomerated population as though it were a panmictic unit of long standing can lead to serious biases in estimates of mutation rates, selection pressures, and effective population sizes. Current DNA studies indicate the presence of numerous genetic variants in human populations. The findings and conclusions of this paper are all fully applicable to the study of genetic variation at the DNA level as well.

Contemporary genetic, historical genetic and morphological data sets for the Yellow-tufted Honeyeater Lichenostomus melanops

Understanding the evolutionary history of threatened populations can improve their conservation management. Re-establishment of past but recent gene flow could re-invigorate threatened populations and replenish genetic diversity, necessary for population persistence. One of the four nominal subspecies of the common yellow-tufted honeyeater, Lichenostomus melanops cassidix, is critically endangered despite substantial conservation efforts over 55 years. Using a combination of morphometric, genetic and modelling approaches we tested for its evolutionary distinctiveness and conservation merit. We confirmed that cassidix has at least one morphometric distinction. It also differs genetically from the other subspecies in allele frequencies but not phylogenetically, implying that its evolution was recent. Modelling historical distribution supported the lack of vicariance and suggested a possibility of gene flow among subspecies at least since the late Pleistocene. Multi-locus coalescent analyses indicated that cassidix diverged from its common ancestor with neighbouring subspecies gippslandicus sometime from the mid-Pleistocene to the Holocene, and that it has the smallest historical effective population size of all subspecies. It appears that cassidix diverged from its ancestor with gippslandicus through a combination of drift and local selection. From patterns of genetic subdivision on two spatial scales and morphological variation we concluded that cassidix, gippslandicus and (melanops + meltoni) are diagnosable as subspecies. Low genetic diversity and effective population size of cassidix may translate to low genetic fitness and evolutionary potential, thus managed gene flow from gippslandicus is recommended for its recovery.

Microsatellite variability and genetic distances.

We analyze the within- and between-population dynamics of the distribution of the number of repeats at multiple microsatellite DNA loci subject to stepwise mutation. Analytical expressions for moments up to the fourth order within a locus and the variance of between-locus variance at mutation-drift equilibrium have been obtained. These statistics may be used to test the appropriateness of the one-step mutation model and to detect between-locus variation in the mutation rate. Published data are compatible with the one-step mutation model, although they do not reject the two-step model. Using both multinomial sampling and diffusion approximations for the analysis of the genetic distance introduced by Goldstein et al. [Goldstein, D. B., Linares, A. R., Cavalli-Sforza, L. L. & Feldman, M. W. (1995) Proc. Natl. Acad. Sci. USA 92, 6723-6727], we show that this distance follows a chi 2 distribution with degrees of freedom equal to the number of loci when there is no variation in mutation rates among the loci. In the presence of such variation, the variance of the distance is obtained. We conclude that the number of microsatellite loci required for the construction of phylogenetic trees with reliable branch lengths may be several hundred. Also, mutations that change repeat scores by several units, even though extremely rare, may dramatically influence estimates of population parameters.

Genetic analysis of complex demographic scenarios: Spatially expanding populations of the cane toad, Bufo marinus

Inferring the spatial expansion dynamics of invading species from molecular data is notoriously difficult due to the complexity of the processes involved. For these demographic scenarios, genetic data obtained from highly variable markers may be profitably combined with specific sampling schemes and information from other sources using a Bayesian approach. The geographic range of the introduced toad Bufo marinus is still expanding in eastern and northern Australia, in each case from isolates established around 1960. A large amount of demographic and historical information is available on both expansion areas. In each area, samples were collected along a transect representing populations of different ages and genotyped at 10 microsatellite loci. Five demographic models of expansion, differing in the dispersal pattern for migrants and founders and in the number of founders, were considered. Because the demographic history is complex, we used an approximate Bayesian method, based on a rejection-regression algorithm. to formally test the relative likelihoods of the five models of expansion and to infer demographic parameters. A stepwise migration-foundation model with founder events was statistically better supported than other four models in both expansion areas. Posterior distributions supported different dynamics of expansion in the studied areas. Populations in the eastern expansion area have a lower stable effective population size and have been founded by a smaller number of individuals than those in the northern expansion area. Once demographically stabilized, populations exchange a substantial number of effective migrants per generation in both expansion areas, and such exchanges are larger in northern than in eastern Australia. The effective number of migrants appears to be considerably lower than that of founders in both expansion areas. We found our inferences to be relatively robust to various assumptions on marker. demographic, and historical features. The method presented here is the only robust, model-based method available so far, which allows inferring complex population dynamics over a short time scale. It also provides the basis for investigating the interplay between population dynamics, drift, and selection in invasive species.

Fisheries stocks from an ecological perspective: Disentangling ecological connectivity from genetic interchange

The concept of a stock of fish as a management unit has been around for well over a hundred years, and this has formed the basis for fisheries science. Methods for delimiting stocks have advanced considerably over recent years, including genetic, telemetric, tagging, geochemical and phenotypic information. In parallel with these developments, concepts in population ecology such as meta-population dynamics and connectivity have advanced. The pragmatic view of stocks has always accepted some mixing during spawning, feeding and/or larval drift. Here we consider the mismatch between ecological connectivity of a matrix of populations typically focussed on demographic measurements, and genetic connectivity of populations that focus on genetic exchange detected using modern molecular approaches. We suggest that from an ecological-connectivity perspective populations can be delimited as management units if there is limited exchange during recruitment or via migration in most years. From a genetic-connectivity perspective such limited exchange can maintain panmixia. We use case-studies of species endangered by overexploitation and/or habitat degradation to show how current methods of stock delimitation can help in managing populations and in conservation.

Fisheries stocks from an ecological perspective: Disentangling ecological connectivity from genetic interchange

The concept of a stock of fish as a management unit has been around for well over a hundred years, and this has formed the basis for fisheries science. Methods for delimiting stocks have advanced considerably over recent years, including genetic, telemetric, tagging, geochemical and phenotypic information. In parallel with these developments, concepts in population ecology such as meta-population dynamics and connectivity have advanced. The pragmatic view of stocks has always accepted some mixing during spawning, feeding and/or larval drift. Here we consider the mismatch between ecological connectivity of a matrix of populations typically focussed on demographic measurements, and genetic connectivity of populations that focus on genetic exchange detected using modern molecular approaches. We suggest that from an ecological-connectivity perspective populations can be delimited as management units if there is limited exchange during recruitment or via migration in most years. From a genetic-connectivity perspective such limited exchange can maintain panmixia. We use case-studies of species endangered by overexploitation and/or habitat degradation to show how current methods of stock delimitation can help in managing populations and in conservation.