Search

Across the tree of life, many organisms are able to reproduce clonally, via vegetative spread, budding or parthenogenesis. In population genetic analyses of clonally reproducing organisms, it is common practice to retain only a single representative per multilocus genotype. Though this practice of clone correction is widespread, the theoretical justification behind it has been very little studied. Here, I use individual-based simulations to study the effect of clone correction on the estimation of the genetic summary statistics H O , H S , F IS , F ST , F'' ST and D est . The simulations follow the standard finite island model, consisting of a set of populations connected by gene flow, but with a variable rate of sexual versus asexual reproduction. The results of the simulations show that by itself, the inclusion of replicated genotypes does not lead to a deviation in the values of the summary statistics, except when the rate of sexual reproduction is less than about one in thousand. However, clone correction can introduce a strong deviation in the values of most of the statistics, when compared to a scenario of full sexual reproduction. For H S and F IS , this deviation can be informative about the process of asexual reproduction, but for F ST , F'' ST and D est , clone correction can lead to incorrect conclusions. I therefore argue that clone correction is not strictly necessary, but can in some cases be insightful. However, when clone correction is applied, it is imperative that results for both the corrected and uncorrected data are presented., (© 2024 The Author(s). Molecular Ecology Resources published by John Wiley & Sons Ltd.)

Despite a long presence in the contiguous United States (US), the distribution of invasive wild pigs (Sus scrofa × domesticus) has expanded rapidly since the 1980s, suggesting a more recent evolutionary shift towards greater invasiveness. Contemporary populations of wild pigs represent exoferal hybrid descendants of domestic pigs and European wild boar, with such hybridization expected to enrich genetic diversity and increase the adaptive potential of populations. Our objective was to characterize how genetic enrichment through hybridization increases the invasiveness of populations by identifying signals of selection and the ancestral origins of selected loci. Our study focused on invasive wild pigs within Great Smoky Mountains National Park, which represents a hybrid population descendent from the admixture of established populations of feral pigs and an introduction of European wild boar to North America. Accordingly, we genotyped 881 wild pigs with multiple high-density single-nucleotide polymorphism (SNP) arrays. We found 233 markers under putative selection spread over 79 regions across 16 out of 18 autosomes, which contained genes involved in traits affecting feralization. Among these, genes were found to be related to skull formation and neurogenesis, with two genes, TYRP1 and TYR, also encoding for crucial melanogenesis enzymes. The most common haplotypes associated with regions under selection for the Great Smoky Mountains population were also common among other populations throughout the region, indicating a key role of putatively selective variants in the fitness of invasive populations. Interestingly, many of these haplotypes were absent among European wild boar reference genotypes, indicating feralization through genetic adaptation., (© 2024 The Authors. Molecular Ecology published by John Wiley & Sons Ltd.)

Analyzing autopolyploid genetic data still presents numerous challenges due to, e.g., missing dosage information of genotypes and the presence of multiple ploidy levels within species or populations, but also because the choice of software is limited when compared to what is available for diploid data. However, over the last years, the number of software programs that can deal with polyploid data is slowly increasing. The software GENODIVE is one of the most widely used programs for the analysis of polyploid genetic data, presenting a wide array of different methods. In this chapter, I outline several frequently used types of population genetic analyses and explain how these apply to polyploid data, including possible pitfalls and biases. I then explain how GENODIVE approaches these analyses and whether and how it can overcome possible biases. Specifically, I focus on analyses of genetic diversity, Hardy-Weinberg equilibrium, quantifying population differentiation, clustering, and calculation of genetic distances. GENODIVE can be downloaded freely from http://www.patrickmeirmans.com/software ., (© 2023. The Author(s), under exclusive license to Springer Science+Business Media, LLC, part of Springer Nature.)

Understanding the impact of historical and demographic processes on genetic variation is essential for devising conservation strategies and predicting responses to climate change. Recolonization after Pleistocene glaciations is expected to leave distinct genetic signatures, characterised by lower genetic diversity in previously glaciated regions. Populations' positions within species ranges also shape genetic variation, following the central-marginal paradigm dictating that peripheral populations are depauperate, sparse and isolated. However, the general applicability of these patterns and relative importance of historical and demographic factors remains unknown. Here, we analysed the distribution of genetic variation in 91 native species of North American plants by coupling microsatellite data and species distribution modelling. We tested the contributions of historical climatic shifts and the central-marginal hypothesis on genetic diversity and structure on the whole data set and across subsets based on taxonomic groups and growth forms. Decreased diversity was found with increased distance from potential glacial refugia, coinciding with the expected make-up of postglacially colonised localities. At the range periphery, lower genetic diversity, higher inbreeding levels and genetic differentiation were reported, following the assumptions of the central-marginal hypothesis. History and demography were found to have approximately equal importance in shaping genetic variation., (© 2021 The Authors. Molecular Ecology published by John Wiley & Sons Ltd.)

The genus Porphyra sensu lato (Bangiaceae, Rhodophyta), an important seaweed grown in aquaculture, is the most genetically diverse group of the Class Bangiophyceae, but has poorly understood genetic variability linked to complex evolutionary processes. Genetic studies in the last decades have largely focused on resolving gene phylogenies; however, there is little information on historical population biogeography, structure and gene flow in the Bangiaceae, probably due to their cryptic nature, chimerism and polyploidy, which render analyses challenging. This study aims to understand biogeographic population structure in the two abundant Porphyra species in the Northeast Atlantic: Porphyra dioica (a dioecious annual) and Porphyra linearis (protandrous hermaphroditic winter annual), occupying distinct niches (seasonality and position on the shore). Here, we present a large-scale biogeographic genetic analysis across their distribution in the Northeast Atlantic, using 10 microsatellites and cpDNA as genetic markers and integrating chimerism and polyploidy, including simulations considering alleles derived from different ploidy levels and/or from different genotypes within the chimeric blade. For P. linearis , both markers revealed strong genetic differentiation of north-central eastern Atlantic populations (from Iceland to the Basque region of Northeast Iberia) vs. southern populations (Galicia in Northwest Iberia, and Portugal), with higher genetic diversity in the south vs. a northern homogenous low diversity. For. P. dioica , microsatellite analyses also revealed two genetic regions, but with weaker differentiation, and cpDNA revealed little structure with all the haplotypes mixed across its distribution. The southern cluster in P. linearis also included introgressed individuals with cpDNA from P. dioica and a winter form of P. dioica occurred spatially intermixed with P. linearis . This third entity had a similar morphology and seasonality as P. linearis but genomes (either nuclear or chloroplast) from P. dioica . We hypothesize a northward colonization from southern Europe (where the ancestral populations reside and host most of the gene pool of these species). In P. linearis recently established populations colonized the north resulting in homogeneous low diversity, whereas for P. dioica the signature of this colonization is not as obvious due to hypothetical higher gene flow among populations, possibly linked to its reproductive biology and annual life history., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Varela-Álvarez, Meirmans, Guiry and Serrão.)

Many sexual-asexual complexes show a distinct pattern where the asexuals have larger and more northerly ranges than closely related sexuals. A prime candidate to explain this so-called "geographical parthenogenesis" is ecological niche divergence between the sexuals and asexuals. Modern niche modelling techniques allow testing niche divergence by directly comparing the niches of sexuals and asexuals. In this study, I use such techniques to perform range-wide tests of whether nine bioclimatic variables, including annual mean temperature and annual precipitation, contribute to geographical parthenogenesis in two dandelion taxa: Taraxacum section Ruderalia and Taraxacum section Erythrosperma, which are both comprised of sexual diploids and asexual triploids. For both sections, I found evidence of niche divergence, though the exact nature of this divergence was different for the two sections. In section Ruderalia, the sexuals preferred warmer and wetter conditions, whereas in section Erythrosperma, the sexuals preferred dryer conditions. Using Species Distribution Modelling, consistent differences between the sexuals and asexuals were found when looking at the niche determinants: the variables that are most important for modelling the distribution. Furthermore, and in contrast with theoretical expectations that predict that the sexuals should have a wider niche, in section Erythrosperma the asexuals were found to have a wider niche than the sexuals. In conclusion, differences in niche optima, niche determinants, and niche width all contribute to the pattern of geographical parthenogenesis of these two dandelion taxa. However, the results also indicate that the exact causation of geographical parthenogenesis is not uniform across taxa., (© 2021 The Authors. Journal of Evolutionary Biology published by John Wiley & Sons Ltd on behalf of European Society for Evolutionary Biology.)

genodive version 3.0 is a user-friendly program for the analysis of population genetic data. This version presents a major update from the previous version and now offers a wide spectrum of different types of analyses. genodive has an intuitive graphical user interface that allows direct manipulation of the data through transformation, imputation of missing data, and exclusion and inclusion of individuals, population and/or loci. Furthermore, genodive seamlessly supports 15 different file formats for importing or exporting data from or to other programs. One major feature of genodive is that it supports both diploid and polyploid data, up to octaploidy (2n = 8x) for some analyses, but up to hexadecaploidy (2n = 16x) for other analyses. The different types of analyses offered by genodive include multiple statistics for estimating population differentiation (φ ST , F ST , F' ST , G ST , G' ST , G'' ST , D est , R ST , ρ), analysis of molecular variance-based K-means clustering, Hardy-Weinberg equilibrium, hybrid index, population assignment, clone assignment, Mantel test, Spatial Autocorrelation, 23 ways of calculating genetic distances, and both principal components and principal coordinates analyses. A unique feature of genodive is that it can also open data sets with nongenetic variables, for example environmental data or geographical coordinates that can be included in the analysis. In addition, genodive makes it possible to run several external programs (lfmm, structure, instruct and vegan) directly from its own user interface, avoiding the need for data reformatting and use of the command line. genodive is available for computers running Mac OS X 10.7 or higher and can be downloaded freely from: http://www.patrickmeirmans.com/software., (© 2020 The Authors. Molecular Ecology Resources published by John Wiley & Sons Ltd.)

Analysis of population genetic structure has become a standard approach in population genetics. In polyploid complexes, clustering analyses can elucidate the origin of polyploid populations and patterns of admixture between different cytotypes. However, combining diploid and polyploid data can theoretically lead to biased inference with (artefactual) clustering by ploidy. We used simulated mixed-ploidy (diploid-autotetraploid) data to systematically compare the performance of k-means clustering and the model-based clustering methods implemented in STRUCTURE, ADMIXTURE, FASTSTRUCTURE and INSTRUCT under different scenarios of differentiation and with different marker types. Under scenarios of strong population differentiation, the tested applications performed equally well. However, when population differentiation was weak, STRUCTURE was the only method that allowed unbiased inference with markers with limited genotypic information (co-dominant markers with unknown dosage or dominant markers). Still, since STRUCTURE was comparatively slow, the much faster but less powerful FASTSTRUCTURE provides a reasonable alternative for large datasets. Finally, although bias makes k-means clustering unsuitable for markers with incomplete genotype information, for large numbers of loci (>1000) with known dosage k-means clustering was superior to FASTSTRUCTURE in terms of power and speed. We conclude that STRUCTURE is the most robust method for the analysis of genetic structure in mixed-ploidy populations, although alternative methods should be considered under some specific conditions.

Studying the genetic population structure of species can reveal important insights into several key evolutionary, historical, demographic, and anthropogenic processes. One of the most important statistical tools for inferring genetic clusters is the program STRUCTURE. Recently, several papers have pointed out that STRUCTURE may show a bias when the sampling design is unbalanced, resulting in spurious joining of underrepresented populations and spurious separation of overrepresented populations. Suggestions to overcome this bias include subsampling and changing the ancestry model, but the performance of these two methods has not yet been tested on actual data. Here, I use a data set of 12 high-alpine plant species to test whether unbalanced sampling affects the STRUCTURE inference of population differentiation between the European Alps and the Carpathians. For four of the 12 species, subsampling of the Alpine populations-to match the sample size between the Alps and the Carpathians-resulted in a drastically different clustering than the full data set. On the other hand, STRUCTURE results with the alternative ancestry model were indistinguishable from the results with the default model. Based on these results, the subsampling strategy seems a more viable approach to overcome the bias than the alternative ancestry model. However, subsampling is only possible when there is an a priori expectation of what constitute the main clusters. Though these results do not mean that the use of STRUCTURE should be discarded, it does indicate that users of the software should be cautious about the interpretation of the results when sampling is unbalanced.

Why and how sexual reproduction is maintained in natural populations, the so-called "queen of problems," is a key unanswered question in evolutionary biology. Recent efforts to solve the problem of sex have often emphasized results generated from laboratory settings. Here, we use a survey of representative "sex in the wild" literature to review and synthesize the outcomes of empirical studies focused on natural populations. Especially notable results included relatively strong support for mechanisms involving niche differentiation and a near absence of attention to adaptive evolution. Support for a major role of parasites is largely confined to a single study system, and only three systems contribute most of the support for mutation accumulation hypotheses. This evidence for taxon specificity suggests that outcomes of particular studies should not be more broadly extrapolated without extreme caution. We conclude by suggesting steps forward, highlighting tests of niche differentiation mechanisms in both laboratory and nature, and empirical evaluation of adaptive evolution-focused hypotheses in the wild. We also emphasize the value of leveraging the growing body of genomic resources for nonmodel taxa to address whether the clearance of harmful mutations and spread of beneficial variants in natural populations proceeds as expected under various hypotheses for sex., (© 2018 The Author(s). Evolution © 2018 The Society for the Study of Evolution.)

The Lesser Antillean Iguana (Iguana delicatissima) is an endangered species threatened by habitat loss and hybridization with non-native Green Iguanas (Iguana iguana). Iguana delicatissima has been extirpated on several islands, and the Green Iguana has invaded most islands with extant populations. Information is essential to protect this species from extinction. We collected data on 293 iguanas including 17 juveniles from St. Eustasius, one of the few remaining I. delicatissima strongholds. Genetic data were leveraged to test for hybridization presence with the Green Iguana using both mitochondrial and nuclear genes, including 16 microsatellite loci. The microsatellites were also analyzed to estimate genetic diversity, population structure, and effective population size. Using molecular and morphological data, we identified 286 I. delicatissima individuals captured during our first fieldwork effort, and 7 non-native iguanas captured during a second effort, showing hybridization occurs within this population. Comparing homologous microsatellites used in studies on Dominica and Chancel, the I. delicatissima population on St. Eustatius has extremely low genetic diversity (HO = 0.051; HE = 0.057), suggesting this population is genetically depauperate. Furthermore, there is significant evidence for inbreeding (FIS = 0.12) and weak spatial genetic structure (FST = 0.021, P = 0.002) within this population. Besides immediate threats including hybridization, this population's low genetic diversity, presence of physiological abnormalities and low recruitment could indicate presence of inbreeding depression that threatens its long-term survival. We conclude there is a continued region-wide threat to I. delicatissima and highlight the need for immediate conservation action to stop the continuing spread of Green Iguanas and to eliminate hybridization from St. Eustatius.

Though polyploidy is an important aspect of the evolutionary genetics of both plants and animals, the development of population genetic theory of polyploids has seriously lagged behind that of diploids. This is unfortunate since the analysis of polyploid genetic data-and the interpretation of the results-requires even more scrutiny than with diploid data. This is because of several polyploidy-specific complications in segregation and genotyping such as tetrasomy, double reduction, and missing dosage information. Here, we review the theoretical and statistical aspects of the population genetics of polyploids. We discuss several widely used types of inferences, including genetic diversity, Hardy-Weinberg equilibrium, population differentiation, genetic distance, and detecting population structure. For each, we point out how the statistical approach, expected result, and interpretation differ between different ploidy levels. We also discuss for each type of inference what biases may arise from the polyploid-specific complications and how these biases can be overcome. From our overview, it is clear that the statistical toolbox that is available for the analysis of genetic data is flexible and still expanding. Modern sequencing techniques will soon be able to overcome some of the current limitations to the analysis of polyploid data, though the techniques are lagging behind those available for diploids. Furthermore, the availability of more data may aggravate the biases that can arise, and increase the risk of false inferences. Therefore, simulations such as we used throughout this review are an important tool to verify the results of analyses of polyploid genetic data.

Many species suffer from anthropogenic habitat fragmentation. The resulting small and isolated populations are more prone to extinction due to, amongst others, genetic erosion, inbreeding depression and Allee-effects. Genetic rescue can help mitigate such problems, but might result in outbreeding depression. We evaluated offspring fitness after selfing and outcrossing within and among three very small and isolated remnant populations of the heterostylous plant Primula vulgaris . We used greenhouse-grown offspring from these populations to test several fitness components. One population was fixed for the pin-morph, and was outcrossed with another population in the field to obtain seeds. Genetic diversity of parent and offspring populations was studied using microsatellites. Morph and population-specific heterosis, inbreeding and outbreeding depression were observed for fruit and seed set, seed weight and cumulative fitness. Highest fitness was observed in the field-outcrossed F1-population, which also showed outbreeding depression following subsequent between-population (back)crossing. Despite outbreeding depression, fitness was still relatively high. Inbreeding coefficients indicated that the offspring were more inbred than their parent populations. Offspring heterozygosity and inbreeding coefficients correlated with observed fitness. One population is evolving homostyly, showing a thrum morph with an elongated style and high autonomous fruit and seed set. This has important implications for conservation strategies such as genetic rescue, as the mating system will be altered by the introduction of homostyles.

Hybridization between species is known to greatly affect their genetic diversity and, therefore, their evolution. Also, within species, there may be genetic clusters between which gene flow is limited, which may impact natural selection. However, few studies have looked simultaneously at the influence of among-species and within-species gene flow. Here, we study the influence of hybridization between Populus balsamifera and Populus trichocarpa on population structure and adaptation in P. balsamifera. We did this by sampling a total of 1517 individuals from across the ranges of these two species, and by genotyping them using a combination of 93 nuclear and 17 cpDNA SNPs. We found that hybridization is mostly limited to the contact zone where the species' distributions overlap. Within P. balsamifera, we found multiple levels of population structure. Interestingly, the border between the Eastern and Central clusters is very sharp, whereas the border between the Central and Western clusters is diffuse. Outlier analysis revealed that three loci associated with the sharp border were also associated with climate. We hypothesize that the observed clusters derive from three refugia during the Pleistocene ice ages. Between the Central and Western clusters, post-glacial long-distance gene flow has led to the diffusion of their border. In the Eastern cluster, we hypothesize that endogenous genomic barriers have developed, leading to the sharp border and a spurious climate association. We conclude that the large-scale genetic structure of P. balsamifera is mostly shaped by historical factors and the influence of interspecific hybridization is limited., (© 2017 European Society For Evolutionary Biology. Journal of Evolutionary Biology © 2017 European Society For Evolutionary Biology.)

Accurately detecting signatures of local adaptation using genetic-environment associations (GEAs) requires controlling for neutral patterns of population structure to reduce the risk of false positives. However, a high degree of collinearity between climatic gradients and neutral population structure can greatly reduce power, and the performance of GEA methods in such case is rarely evaluated in empirical studies. In this study, we attempted to disentangle the effects of local adaptation and isolation by environment (IBE) from those of isolation by distance (IBD) and isolation by colonization from glacial refugia (IBC) using range-wide samples in two white pine species. For this, SNPs from 168 genes, including 52 candidate genes for growth and phenology, were genotyped in 133 and 61 populations of Pinus strobus and P. monticola , respectively. For P. strobus and using all 153 SNPs, climate (IBE) did not significantly explained among-population variation when controlling for IBD and IBC in redundancy analyses (RDAs). However, 26 SNPs were significantly associated with climate in single-locus GEA analyses (Bayenv2 and LFMM), suggesting that local adaptation took place in the presence of high gene flow. For P. monticola , we found no evidence of IBE using RDAs and weaker signatures of local adaptation using GEA and F ST outlier tests, consistent with adaptation via phenotypic plasticity. In both species, the majority of the explained among-population variation (69 to 96%) could not be partitioned between the effects of IBE, IBD, and IBC. GEA methods can account differently for this confounded variation, and this could explain the small overlap of SNPs detected between Bayenv2 and LFMM. Our study illustrates the inherent difficulty of taking into account neutral structure in natural populations and the importance of sampling designs that maximize climatic variation, while minimizing collinearity between climatic gradients and neutral structure.

We used flow cytometry, chromosome counting and AFLP markers to investigate gene flow from the crop plant oilseed rape, Brassica napus (AACC) to wild B. rapa (AA) in the Netherlands. From 89 B. napus source populations investigated, all near cropping fields or at transhipment sites, only 19 contained a B. rapa population within a 2.5-km radius. During our survey we found only three populations with F1 hybrids (AAC), as recognized by their nine extra chromosomes and by flow cytometry. These hybrids were all collected in mixed populations where the two species grew in close proximity. Populations with F1 hybrids were not close to crops, but instead were located on road verges with highly disturbed soils, in which both species were probably recruited from the soil seed bank. Many plants in the F2, BC1 or higher backcrosses are expected to carry one to eight C chromosomes. However, these plants were not observed among the hybrids. We further investigated introgression with molecular markers (AFLP) and compared sympatric B. rapa populations (near populations of B. napus) with control populations of B. rapa (no B. napus within at least 7 km). We found no difference between sympatric and control populations in the number of C markers in B. rapa, nor did we find that these sympatric populations closely resembled B. napus. Our data show that hybrids occur but also suggest no recent introgression of alleles from the crop plant B. napus into wild B. rapa in the Dutch populations studied., (© 2014 German Botanical Society and The Royal Botanical Society of the Netherlands.)

The estimation of migration rates using molecular markers is an important aspect of many population genetic studies. Several different methods are available for estimating migration, but most of these make multiple limiting assumptions. One method that is relatively free from assumptions is BAYESASS, which uses assignment methods in a Bayesian framework. However, when tested using simulated data, this method was found to have problems with the convergence of the Markov chain Monte Carlo. Here, I perform a literature study to test whether these convergence problems are also present when BAYESASS is used to estimate migration rates from empirical data. A review of 100 studies that have used BAYESASS shows that this is indeed the case. The estimated proportions of nonmigrants were mostly either close to 2/3 or 1, indicating that the MCMC tends to get trapped near the bounds of the prior distribution. In addition, I found that the quality of the inference was negatively affected by the number of sampled populations, but increased with increasing numbers of sampled individuals and with the strength of the population structure as measured by FST . Based on these results, I give several recommendations that should help to reduce problems when using BAYESASS with empirical data. Most importantly, I argue that researchers should be more realistic about inferences of migration rates and that BAYESASS will give optimal results only when the experiment has been especially designed around its use., (© 2013 John Wiley & Sons Ltd.)

Background and Aims: Transgene introgression from crops into wild relatives may increase the resistance of wild plants to herbicides, insects, etc. The chance of transgene introgression depends not only on the rate of hybridization and the establishment of hybrids in local wild populations, but also on the metapopulation dynamics of the wild relative. The aim of the study was to estimate gene flow in a metapopulation for assessing and managing the risks of transgene introgression., Methods: Wild carrots (Daucus carota) were sampled from 12 patches in a metapopulation. Eleven microsatellites were used to genotype wild carrots. Genetic structure was estimated based on the FST statistic. Contemporary (over the last several generations) and historical (over many generations) gene flow was estimated with assignment and coalescent methods, respectively., Key Results: The genetic structure in the wild carrot metapopulation was moderate (FST = 0·082) and most of the genetic variation resided within patches. A pattern of isolation by distance was detected, suggesting that most of the gene flow occurred between neighbouring patches (≤1 km). The mean contemporary gene flow was 5 times higher than the historical estimate, and the correlation between them was very low. Moreover, the contemporary gene flow in roadsides was twice that in a nature reserve, and the correlation between contemporary and historical estimates was much higher in the nature reserve. Mowing of roadsides may contribute to the increase in contemporary gene flow. Simulations demonstrated that the higher contemporary gene flow could accelerate the process of transgene introgression in the metapopulation., Conclusions: Human disturbance such as mowing may alter gene flow patterns in wild populations, affecting the metapopulation dynamics of wild plants and the processes of transgene introgression in the metapopulation. The risk assessment and management of transgene introgression and the control of weeds need to take metapopulation dynamics into consideration.

Polyploids are traditionally classified into allopolyploids and autopolyploids, based on their evolutionary origin and their disomic or multisomic mode of inheritance. Over the past decade it has become increasingly clear that there is a continuum between disomic and multisomic inheritance, with the rate of tetrasomy differing among species and among chromosomes within species. Here, we use a simple population genetic model to study the impact of the mode of inheritance on the genetic diversity and population divergence of tetraploids. We found that under almost strict disomic inheritance the tetraploid genome is divided into two separate subgenomes, such as found in classical allopolyploids. In those cases, assuming full tetrasomy in the analysis of polyploid genetic data will lead to an important bias in estimates of genetic diversity and population divergence. However, we found that even a low rate of allele exchange between the two subgenomes, at about one event per generation, is sufficient to homogenise the allele frequencies over the subgenomes, and the estimates become essentially unbiased. The inbreeding coefficient F(IS) can then be used to detect whether the estimates of diversity and divergence will be biased when full multisomy is assumed. Finally, we found that different summary statistics for measuring the strength of population differentiation are differentially affected by a deviation from full tetrasomy. Our model results provide several useful guidelines for the analysis of polyploid data, helping researchers to determine when their inferences are biased and which summary statistics to use.

Introgression is the permanent incorporation of genes from the genome of one population into another. This can have severe consequences, such as extinction of endemic species, or the spread of transgenes. Quantification of the risk of introgression is an important component of genetically modified crop regulation. Most theoretical introgression studies aimed at such quantification disregard one or more of the most important factors concerning introgression: realistic genetical mechanisms, repeated invasions and stochasticity. In addition, the use of linkage as a risk mitigation strategy has not been studied properly yet with genetic introgression models. Current genetic introgression studies fail to take repeated invasions and demographic stochasticity into account properly, and use incorrect measures of introgression risk that can be manipulated by arbitrary choices. In this study, we present proper methods for risk quantification that overcome these difficulties. We generalize a probabilistic risk measure, the so-called hazard rate of introgression, for application to introgression models with complex genetics and small natural population sizes. We illustrate the method by studying the effects of linkage and recombination on transgene introgression risk at different population sizes.

Determining the genetic structure of populations is becoming an increasingly important aspect of genetic studies. One of the most frequently used methods is the calculation of F-statistics using an Analysis of Molecular Variance (AMOVA). However, this has the drawback that the population hierarchy has to be known a priori. Therefore, the population structure is often based on the results of a clustering analysis. Here I show how these two steps, clustering and calculation of F-statistics, can be combined in a single analysis. I do this by showing how the AMOVA framework is theoretically related to the widely used method of K-means clustering and can be used for the clustering of populations into groups. Simulations were used to show that the method performed very well both under random mating and under nonrandom mating. However, when the migration rates were high, the results were better under random mating than under predominant selfing or clonal reproduction. Two summary statistics were tested for estimating the number of clusters. Overall, pseudo-F showed the better performance, but BIC is better for detecting whether any significant structure is present. The results show that the AMOVA-based K-means clustering is useful for clustering population genetic data. Programs to perform the clustering can be downloaded from www.patrickmeirmans.com/software.

The genetic population structure of many species is characterised by a pattern of isolation by distance (IBD): due to limited dispersal, individuals that are geographically close tend to be genetically more similar than individuals that are far apart. Despite the ubiquity of IBD in nature, many commonly used statistical tests are based on a null model that is completely non-spatial, the Island model. Here, I argue that patterns of spatial autocorrelation deriving from IBD present a problem for such tests as it can severely bias their outcome. I use simulated data to illustrate this problem for two widely used types of tests: tests of hierarchical population structure and the detection of loci under selection. My results show that for both types of tests the presence of IBD can indeed lead to a large number of false positives. I therefore argue that all analyses in a study should take the spatial dependence in the data into account, unless it can be shown that there is no spatial autocorrelation in the allele frequency distribution that is under investigation. Thus, it is urgent to develop additional statistical approaches that are based on a spatially explicit null model instead of the non-spatial Island model., (© 2012 Blackwell Publishing Ltd.)

Understanding the maintenance of sexual reproduction constitutes a difficult problem for evolutionary biologists because of the immediate costs that sex seems to incur. Typically, general benefits to sex and recombination are investigated that might outweigh these costs. However, several factors can strongly influence the complex balance between costs and benefits of sex; these include constraints on the evolution of asexuality, ecological differentiation, and certain lif-history traits. We review these factors and their empirical support for the first time in a unified framework and find that they can reduce the costs of sex, circumvent them, or make them inapplicable. These factors can even tip the scales to a net benefit for sex. The reviewed factors affect species and species groups differently, and we conclude consequently that understanding the maintenance of sex could turn out to be more species-specific than commonly assumed. Interestingly, our study suggests that, in some species, no general benefits to sex and recombination might be needed to understand the maintenance of sex, as in our case study of dandelions.

A plant species' genetic population structure is the result of a complex combination of its life history, ecological preferences, position in the ecosystem and historical factors. As a result, many different statistical methods exist that measure different aspects of species' genetic structure. However, little is known about how these methods are interrelated and how they are related to a species' ecology and life history. In this study, we used the IntraBioDiv amplified fragment length polymorphisms data set from 27 high-alpine species to calculate eight genetic summary statistics that we jointly correlate to a set of six ecological and life-history traits. We found that there is a large amount of redundancy among the calculated summary statistics and that there is a significant association with the matrix of species traits. In a multivariate analysis, two main aspects of population structure were visible among the 27 species. The first aspect is related to the species' dispersal capacities and the second is most likely related to the species' postglacial recolonization of the Alps. Furthermore, we found that some summary statistics, most importantly Mantel's r and Jost's D, show different behaviour than expected based on theory. We therefore advise caution in drawing too strong conclusions from these statistics., (© 2011 Blackwell Publishing Ltd.)

Background: The genomic architecture of bud phenology and height growth remains poorly known in most forest trees. In non model species, QTL studies have shown limited application because most often QTL data could not be validated from one experiment to another. The aim of our study was to overcome this limitation by basing QTL detection on the construction of genetic maps highly-enriched in gene markers, and by assessing QTLs across pedigrees, years, and environments., Results: Four saturated individual linkage maps representing two unrelated mapping populations of 260 and 500 clonally replicated progeny were assembled from 471 to 570 markers, including from 283 to 451 gene SNPs obtained using a multiplexed genotyping assay. Thence, a composite linkage map was assembled with 836 gene markers.For individual linkage maps, a total of 33 distinct quantitative trait loci (QTLs) were observed for bud flush, 52 for bud set, and 52 for height growth. For the composite map, the corresponding numbers of QTL clusters were 11, 13, and 10. About 20% of QTLs were replicated between the two mapping populations and nearly 50% revealed spatial and/or temporal stability. Three to four occurrences of overlapping QTLs between characters were noted, indicating regions with potential pleiotropic effects. Moreover, some of the genes involved in the QTLs were also underlined by recent genome scans or expression profile studies.Overall, the proportion of phenotypic variance explained by each QTL ranged from 3.0 to 16.4% for bud flush, from 2.7 to 22.2% for bud set, and from 2.5 to 10.5% for height growth. Up to 70% of the total character variance could be accounted for by QTLs for bud flush or bud set, and up to 59% for height growth., Conclusions: This study provides a basic understanding of the genomic architecture related to bud flush, bud set, and height growth in a conifer species, and a useful indicator to compare with Angiosperms. It will serve as a basic reference to functional and association genetic studies of adaptation and growth in Picea taxa. The putative QTNs identified will be tested for associations in natural populations, with potential applications in molecular breeding and gene conservation programs. QTLs mapping consistently across years and environments could also be the most important targets for breeding, because they represent genomic regions that may be least affected by G × E interactions.

MARLIN is a software to create, run, analyse, and visualize spatially explicit population genetic simulations. It provides an intuitive user interface with which the geographical layout of a metapopulation can be drawn by hand or loaded from a map. Furthermore, the interface allows easy selection of the many different simulation settings. MARLIN then uses the program QuantiNemo to run the simulation in the background. When simulations are finished, MARLIN directly analyses and plots the results, thereby greatly simplifying the simulation workflow. This combination of simulation and analysis makes MARLIN ideal for teaching and for scientists who are interested in doing simulations without having to learn command-line operations. MARLIN is available for computers running Mac OS X and can be downloaded from: http://www.patrickmeirmans.com/software., (© 2010 Blackwell Publishing Ltd.)

Although F(ST) is widely used as a measure of population structure, it has been criticized recently because of its dependency on within-population diversity. This dependency can lead to difficulties in interpretation and in the comparison of estimates among species or among loci and has led to the development of two replacement statistics, F'(ST) and D. F'(ST) is the normal F(ST) standardized by the maximum value it can obtain, given the observed within-population diversity. D uses a multiplicative partitioning of diversity, based on the effective number of alleles rather than on the expected heterozygosity. In this study, we review the relationships between the three classes of statistics (F(ST), F'(ST) and D), their estimation and their properties. We illustrate the relationships between the statistics using a data set of estimates from 84 species taken from the last 4 years of Molecular Ecology. As with F(ST), unbiased estimators are available for the two new statistics D and F'(ST). Here, we develop a new unbiased F'(ST) estimator based on G(ST), which we call G''(ST). However, F'(ST) can be calculated using any F(ST) estimator for which the maximum value can be obtained. As all three statistics have their advantages and their drawbacks, we recommend continued use of F(ST) in combination with either F'(ST) or D. In most cases, F'(ST) would be the best choice among the latter two as it is most suited for inferences of the influence of demographic processes such as genetic drift and migration on genetic population structure., (© 2010 Blackwell Publishing Ltd.)

Premise of the Study: Poplars and their hybrids are seen as important candidates for bioenergy initiatives. However, many concerns have been raised about large-scale plantations of new poplar cultivars. The deployment of such plants with novel traits brings the risk of potential spread of novel genome regions (including exotic genes, transgenes, or other heritable modifications) into natural populations of related species. The possibility of introgression is especially high in poplars because reproductive barriers between species are weak. Knowledge of the frequency of hybridization between cultivated trees and natural populations is one important step in the risk-assessment process. •, Methods: We studied the rate of spontaneous hybridization from two sexually mature poplar plantations into adjacent natural populations of Populus deltoides and P. balsamifera. The two plantations, both in eastern Canada, contain many different complex hybrid clones with components from exotic species, mostly P. nigra, P. trichocarpa, and P. maximowiczii. We analyzed 12 species-specific single nucleotide polymorphisms from six different genes in 5373 offspring sampled from the natural populations. •, Results: Contributions from all three exotics were found in the offspring, confirming low reproductive barriers among poplar species in these sections. The frequency of hybrid offspring varied among pollen donors, recipient populations, and years. •, Conclusions: The remarkably high rate of hybridization that was found in the smallest natural population sampled suggests that small peripheral populations carry a higher risk of introgression. These results could be used as a starting point for developing regulatory guidelines for the introduction of plants with novel traits.

As the evolutionary significance of hybridization is largely dictated by its extent beyond the first generation, we broadly surveyed patterns of introgression across a sympatric zone of two native poplars (Populus balsamifera, Populus deltoides) in Quebec, Canada within which European exotic Populus nigra and its hybrids have been extensively planted since the 1800s. Single nucleotide polymorphisms (SNPs) that appeared fixed within each species were characterized by DNA-sequencing pools of pure individuals. Thirty-five of these diagnostic SNPs were employed in a high-throughput assay that genotyped 635 trees of different age classes, sampled from 15 sites with various degrees of anthropogenic disturbance. The degree of admixture within sampled trees was then assessed through Bayesian clustering of genotypes. Hybrids were present in seven of the populations, with 2.4% of all sampled trees showing spontaneous admixture. Sites with hybrids were significantly more disturbed than pure stands, while hybrids comprised both immature juveniles and trees of reproductive age. All three possible F1s were detected. Advanced-generation hybrids were consistently biased towards P. balsamifera regardless of whether hybridization had occurred with P. deltoides or P. nigra. Gene exchange between P. deltoides and P. nigra was not detected beyond the F1 generation; however, detection of a trihybrid demonstrates that even this apparent reproductive isolation does not necessarily result in an evolutionary dead end. Collectively, results demonstrate the natural fertility of hybrid poplars and suggest that introduced genes could potentially affect the genetic integrity of native trees, similar to that arising from introgression between natives.

Sexual reproduction is both extremely costly and extremely common relative to asexuality, indicating that it must confer profound benefits. This in turn points to major disadvantages of asexual reproduction, which is usually given as an explanation for why almost all asexual lineages are apparently quite short-lived. However, a growing body of evidence suggests that some asexual lineages are actually quite old. Insight into why sex is so common may come from understanding why asexual lineages persist in some places or taxa but not others. Here, we review the distribution of asexual lineage ages estimated from a diverse array of taxa, and we discuss our results in light of the main mutational and environmental hypotheses for sex. Along with strengthening the case for wide variation in asexual lineage age and the existence of many old asexual taxa, we also found that the distribution of asexual lineage age estimates follows a surprisingly regular distribution, to the extent that asexual taxa viewed as "scandalously" ancient merely fall on the high end of this distribution. We interpret this result to mean that similar mechanisms may determine asexual lineage age across eukaryotic taxa. We also derive some qualitative predictions for asexual lineage age under different theories for sex and discuss empirical evidence for these predictions. Ultimately, we were limited in the extent to which we could use these data to make inferences about the maintenance of sex by the absence of both clear theoretical expectations and estimates of key parameters.

Most models on introgression from genetically modified (GM) plants have focused on small spatial scales, modelling gene flow from a field containing GM plants into a single adjacent population of a wild relative. Here, we present a model to study the effect of introgression from multiple plantations into the whole metapopulation of the wild relative. The most important result of the model is that even very low levels of introgression and selection can lead to a high probability that the transgene goes to fixation in the metapopulation. Furthermore, the overall frequency of the transgene in the metapopulation, after a certain number of generations of introgression, depends on the population dynamics. If there is a high rate of migration or a high rate of population turnover, the overall transgene frequency is much higher than with lower rates. However, under an island model of population structure, this increased frequency has only a very small effect on the probability of fixation of the transgene. Considering these results, studies on the potential ecological risks of introgression from GM plants should look not only at the rate of introgression and selection acting on the transgene, but also at the metapopulation dynamics of the wild relative.

Comparison of population structure between studies can be difficult, because the value of the often-used FST-statistic depends on the amount of genetic variation within populations. Recently, a standardized measure of genetic differentiation was developed based on GST, which addressed this problem, though no method was provided to estimate this standardized measure without bias. Here I present a method to estimate a standardized measure of population differentiation based on the analysis of molecular variance framework. One advantage of the method is that it can be readily expanded to include different hierarchical levels in the tested population structure.

Male reproductive output, pollen in plants and sperm in animals has been shown to constitute a substantial cost for many organisms. In parthenogenetic hermaphrodites, selection is therefore expected to reduce the allocation of resources to male reproductive output. However, sustained production of pollen or sperm has been observed in numerous asexual hermaphrodites. We studied the widespread production of pollen by triploid asexual dandelions, Taraxacum sect. Ruderalia, comparing rare male sterile individuals with pollen producing asexuals. We found that individuals can show plasticity in the production of pollen, but that it is nevertheless possible to distinguish between (facultatively) male sterile asexuals and male fertile asexuals. Based on evidence from genetic markers and crosses, we conclude that the male sterility in asexual dandelions is caused by nuclear genes, in contrast to the cytoplasmically inherited male sterility previously found in sexual dandelions. Male sterile lineages did not produce more seeds per flower head, heavier seeds or seeds that were more viable. However, male sterile plants did produce more seed heads and hence more seeds than pollen producing ones, indicating that they were able to reallocate resources toward seed production. Considering the difference in seed production, it remains puzzling that not more asexual dandelions are male sterile.

Ecological differentiation is widely seen as an important factor enabling the stable coexistence of closely related plants of different ploidy levels. We studied ecological and genetic differentiation between co-occurring sexual diploid and apomictic triploid Taraxacum section Ruderalia by analysing spatial patterns both in the distribution of cytotypes and in the distribution of genetic variation within and between the cytotypes. A significant relationship between ploidy level and elevation was found. This mode of ecological differentiation however, was not sufficient to explain the significant spatial structure in the distribution of diploids and triploids within the population. Strong congruence was found between the spatial genetic patterns within the diploids and within the triploids. We argue that this congruence is an indication of gene flow between neighbouring plants of different ploidy levels.