Analysis of the genetic diversity and structure across a wide range of germplasm reveals prominent gene flow in apple at the European level

Abstract Background: The amount and structure of genetic diversity in dessert apple germplasm conserved at a European level is mostly unknown, since all diversity studies conducted in Europe until now have been performed on regional or national collections. Here, we applied a common set of 16 SSR markers to genotype more than 2,400 accessions across 14 collections representing three broad European geographic regions (North+East, West and South) with the aim to analyze the extent, distribution and structure of variation in the apple genetic resources in Europe. Results: A Bayesian model-based clustering approach showed that diversity was organized in three groups, although these were only moderately differentiated (FST=0.031). A nested Bayesian clustering approach allowed identification of subgroups which revealed internal patterns of substructure within the groups, allowing a finer delineation of the variation into eight subgroups (FST=0.044). The first level of stratification revealed an asymmetric division of the germplasm among the three groups, and a clear association was found with the geographical regions of origin of the cultivars. The substructure revealed clear partitioning of genetic groups among countries, but also interesting associations between subgroups and breeding purposes of recent cultivars or particular usage such as cider production. Additional parentage analyses allowed us to identify both putative parents of more than 40 old and/or local cultivars giving interesting insights in the pedigree of some emblematic cultivars. Conclusions: The variation found at group and sub-group levels may reflect a combination of historical processes of migration/selection and adaptive factors to diverse agricultural environments that, together with genetic drift, have resulted in extensive genetic variation but limited population structure. The European dessert apple germplasm represents an important source of genetic diversity with a strong historical and patrimonial value. The present work thus constitutes a decisive step in the field of conservation genetics. Moreover, the obtained data can be used for defining a European apple core collection useful for further identification of genomic regions associated with commercially important horticultural traits in apple through genome-wide association studies.

Genetic structure and systematic relationships within the Ophrys fuciflora aggregate (Orchidaceae: Orchidinae): high diversity in Kent and a wind-induced discontinuity bisecting the Adriatic

A recent phylogenetic study based on multiple datasets is used as the framework for a more detailed examination of one of the ten molecularly circumscribed groups identified, the Ophrys fuciflora aggregate. The group is highly morphologically variable, prone to phenotypic convergence, shows low levels of sequence divergence and contains an unusually large proportion of threatened taxa, including the rarest Ophrys species in the UK. The aims of this study were to (a) circumscribe minimum resolvable genetically distinct entities within the O. fuciflora aggregate, and (b) assess the likelihood of gene flow between genetically and geographically distinct entities at the species and population levels. Fifty-five accessions sampled in Europe and Asia Minor from the O. fuciflora aggregate were studied using the AFLP genetic fingerprinting technique to evaluate levels of infraspecific and interspecific genetic variation and to assess genetic relationships between UK populations of O. fuciflora s.s. in Kent and in their continental European and Mediterranean counterparts. The two genetically and geographically distinct groups recovered, one located in England and central Europe and one in south-eastern Europe, are incongruent with current species delimitation within the aggregate as a whole and also within O. fuciflora s.s. Genetic diversity is higher in Kent than in the rest of western and central Europe. Gene flow is more likely to occur between populations in closer geographical proximity than those that are morphologically more similar. Little if any gene flow occurs between populations located in the south-eastern Mediterranean and those dispersed throughout the remainder of the distribution, revealing a genetic discontinuity that runs north-south through the Adriatic. This discontinuity is also evident in other clades of Ophrys and is tentatively attributed to the long-term influence of prevailing winds on the long-distance distribution of pollinia and especially seeds. A cline of gene flow connects populations from Kent and central and southern Europe; these individuals should therefore be considered part of an extensive meta-population. Gene flow is also evident among populations from Kent, which appear to constitute a single metapopulation. They show some evidence of hybridization, and possibly also introgression, with O. apifera.

The role of satellite image-processing for national-scale estimates of gene flow from genetically modified crops: rapeseed in the UK as a model.

1. There is concern over the possibility of unwanted environmental change following transgene movement from genetically modified (GM) rapeseed Brassica napus to its wild and weedy relatives. 2. The aim of this research was to develop a remote sensing-assisted methodology to help quantify gene flow from crops to their wild relatives over wide areas. Emphasis was placed on locating sites of sympatry, where the frequency of gene flow is likely to be highest, and on measuring the size of rapeseed fields to allow spatially explicit modelling of wind-mediated pollen-dispersal patterns. 3. Remote sensing was used as a tool to locate rapeseed fields, and a variety of image-processing techniques was adopted to facilitate the compilation of a spatially explicit profile of sympatry between the crop and Brassica rapa. 4. Classified satellite images containing rapeseed fields were first used to infer the spatial relationship between donor rapeseed fields and recipient riverside B. rapa populations. Such images also have utility for improving the efficiency of ground surveys by identifying probable sites of sympatry. The same data were then also used for the calculation of mean field size. 5. This paper forms a companion paper to Wilkinson et al. (2003), in which these elements were combined to produce a spatially explicit profile of hybrid formation over the UK. The current paper demonstrates the value of remote sensing and image processing for large-scale studies of gene flow, and describes a generic method that could be applied to a variety of crops in many countries. 6. Synthesis and applications. The decision to approve or prevent the release of a GM cultivar is made at a national rather than regional level. It is highly desirable that data relating to the decision-making process are collected at the same scale, rather than relying on extrapolation from smaller experiments designed at the plot, field or even regional scale. It would be extremely difficult and labour intensive to attempt to carry out such large-scale investigations without the use of remote-sensing technology. This study used rapeseed in the UK as a model to demonstrate the value of remote sensing in assembling empirical information at a national level.

Landscape gene flow, coexistence and threshold effect: The case of genetically modified herbicide tolerant oilseed rape (Brassica napus)

Globally there have been a number of concerns about the development of genetically modified crops many of which relate to the implications of gene flow at various levels. In Europe these concerns have led the European Union (EU) to promote the concept of 'coexistence' to allow the freedom to plant conventional and genetically modified (GM) varieties but to minimise the presence of transgenic material within conventional crops. Should a premium for non-GM varieties emerge on the market, the presence of transgenes would generate a 'negative externality' to conventional growers. The establishment of maximum tolerance level for the adventitious presence of GM material in conventional crops produces a threshold effect in the external costs. The existing literature suggests that apart from the biological characteristics of the plant under consideration (e.g. self-pollination rates, entomophilous species, anemophilous species, etc.), gene flow at the landscape level is affected by the relative size of the source and sink populations and the spatial arrangement of the fields in the landscape. In this paper, we take genetically modified herbicide tolerant oilseed rape (GM HT OSR) as a model crop. Starting from an individual pollen dispersal function, we develop a spatially explicit numerical model in order to assess the effect of the size of the source/sink populations and the degree of spatial aggregation on the extent of gene flow into conventional OSR varieties under two alternative settings. We find that when the transgene presence in conventional produce is detected at the field level, the external cost will increase with the size of the source area and with the level of spatial disaggregation. on the other hand when the transgene presence is averaged among all conventional fields in the landscape (e.g. because of grain mixing before detection), the external cost will only depend on the relative size of the source area. The model could readily be incorporated into an economic evaluation of policies to regulate adoption of GM HT OSR. (c) 2007 Elsevier B.V. All rights reserved.

Gene flow, risk assessment and the environmental release of transgenic plants

The release of genetically modified plants is governed by regulations that aim to provide an assessment of potential impact on the environment. One of the most important components of this risk assessment is an evaluation of the probability of gene flow. In this review, we provide an overview of the current literature on gene flow from transgenic plants, providing a framework of issues for those considering the release of a transgenic plant into the environment. For some plants gene flow from transgenic crops is well documented, and this information is discussed in detail in this review. Mechanisms of gene flow vary from plant species to plant species and range from the possibility of asexual propagation, short- or long-distance pollen dispersal mediated by insects or wind and seed dispersal. Volunteer populations of transgenic plants may occur where seed is inadvertently spread during harvest or commercial distribution. If there are wild populations related to the transgenic crop then hybridization and eventually introgression in the wild may occur, as it has for herbicide resistant transgenic oilseed rape (Brassica napus). Tools to measure the amount of gene flow, experimental data measuring the distance of pollen dispersal, and experiments measuring hybridization and seed survivability are discussed in this review. The various methods that have been proposed to prevent gene flow from genetically modified plants are also described. The current "transgenic traits'! in the major crops confer resistance to herbicides and certain insects. Such traits could confer a selective advantage (an increase in fitness) in wild plant populations in some circumstances, were gene flow to occur. However, there is ample evidence that gene flow from crops to related wild species occurred before the development of transgenic crops and this should be taken into account in the risk assessment process.

The population genetic structure of clonal organisms generated by exponentially bounded and fat-tailed dispersal

Long distance dispersal (LDD) plays an important role in many population processes like colonization, range expansion, and epidemics. LDD of small particles like fungal spores is often a result of turbulent wind dispersal and is best described by functions with power-law behavior in the tails ("fat tailed"). The influence of fat-tailed LDD on population genetic structure is reported in this article. In computer simulations, the population structure generated by power-law dispersal with exponents in the range of -2 to -1, in distinct contrast to that generated by exponential dispersal, has a fractal structure. As the power-law exponent becomes smaller, the distribution of individual genotypes becomes more self-similar at different scales. Common statistics like G(ST) are not well suited to summarizing differences between the population genetic structures. Instead, fractal and self-similarity statistics demonstrated differences in structure arising from fat-tailed and exponential dispersal. When dispersal is fat tailed, a log-log plot of the Simpson index against distance between subpopulations has an approximately constant gradient over a large range of spatial scales. The fractal dimension D-2 is linearly inversely related to the power-law exponent, with a slope of similar to -2. In a large simulation arena, fat-tailed LDD allows colonization of the entire space by all genotypes whereas exponentially bounded dispersal eventually confines all descendants of a single clonal lineage to a relatively small area.

Spontaneous gene flow from rapeseed (Brassica napus) to wild Brassica oleracea

Research on the environmental risks of gene flow from genetically modified ( GM) crops to wild relatives has traditionally emphasized recipients yielding most hybrids. For GM rapeseed ( Brassica napus), interest has centred on the 'frequently hybridizing' Brassica rapa over relatives such as Brassica oleracea, where spontaneous hybrids are unreported in the wild. In two sites, where rapeseed and wild B. oleracea grow together, we used flow cytometry and crop-specific microsatellite markers to identify one triploid F-1 hybrid, together with nine diploid and two near triploid introgressants. Given the newly discovered capacity for spontaneous introgression into B. oleracea, we then surveyed associated flora and fauna to evaluate the capacity of both recipients to harm cohabitant species with acknowledged conservational importance. Only B. oleracea occupies rich communities containing species afforded legislative protection; these include one rare micromoth species that feeds on B. oleracea and warrants further assessment. We conclude that increased attention should now focus on B. oleracea and similar species that yield few crop-hybrids, but possess scope to affect rare or endangered associates.

The role of satellite image-processing for national-scale estimates of gene flow from genetically modified crops: rapeseed in the UK as a model

1.There is concern over the possibility of unwanted environmental change following transgene movement from genetically modified (GM) rapeseed Brassica napus to its wild and weedy relatives. 2. The aim of this research was to develop a remote sensing-assisted methodology to help quantify gene flow from crops to their wild relatives over wide areas. Emphasis was placed on locating sites of sympatry, where the frequency of gene flow is likely to be highest, and on measuring the size of rapeseed fields to allow spatially explicit modelling of wind-mediated pollen-dispersal patterns. 3. Remote sensing was used as a tool to locate rapeseed fields, and a variety of image-processing techniques was adopted to facilitate the compilation of a spatially explicit profile of sympatry between the crop and Brassica rapa. 4. Classified satellite images containing rapeseed fields were first used to infer the spatial relationship between donor rapeseed fields and recipient riverside B. rapa populations. Such images also have utility for improving the efficiency of ground surveys by identifying probable sites of sympatry. The same data were then also used for the calculation of mean field size. 5. This paper forms a companion paper to Wilkinson et al. (2003), in which these elements were combined to produce a spatially explicit profile of hybrid formation over the UK. The current paper demonstrates the value of remote sensing and image processing for large-scale studies of gene flow, and describes a generic method that could be applied to a variety of crops in many countries. 6.Synthesis and applications. The decision to approve or prevent the release of a GM cultivar is made at a national rather than regional level. It is highly desirable that data relating to the decision-making process are collected at the same scale, rather than relying on extrapolation from smaller experiments designed at the plot, field or even regional scale. It would be extremely difficult and labour intensive to attempt to carry out such large-scale investigations without the use of remote-sensing technology. This study used rapeseed in the UK as a model to demonstrate the value of remote sensing in assembling empirical information at a national level.

Genetic structure and linkage disequilibrium in landrace populations of barley in Sardinia

Multilocus digenic linkage disequilibria (LD) and their population structure were investigated in eleven landrace populations of barley (Hordeum vulgare ssp. vulgare L.) in Sardinia, using 134 dominant simple-sequence amplified polymorphism markers. The analysis of molecular variance for these markers indicated that the populations were partially differentiated (F ST = 0.18), and clustered into three geographic areas. Consistent with this population pattern, STRUCTURE analysis allocated individuals from a bulk of all populations into four genetic groups, and these groups also showed geographic patterns. In agreement with other molecular studies in barley, the general level of LD was low (13 % of locus pairs, with P < 0.01) in the bulk of 337 lines, and decayed steeply with map distance between markers. The partitioning of multilocus associations into various components indicated that genetic drift and founder effects played a major role in determining the overall genetic makeup of the diversity in these landrace populations, but that epistatic homogenising or diversifying selection was also present. Notably, the variance of the disequilibrium component was relatively high, which implies caution in the pooling of barley lines for association studies. Finally, we compared the analyses of multilocus structure in barley landrace populations with parallel analyses in both composite crosses of barley on the one hand and in natural populations of wild barley on the other. Neither of these serves as suitable mimics of landraces in barley, which require their own study. Overall, the results suggest that these populations can be exploited for LD mapping if population structure is controlled.

Genetic evidence of frequent long-distance recruitment in a vertebrate-dispersed tree

The importance of dispersal for the maintenance of biodiversity, while long-recognized, has remained unresolved. We used molecular markers to measure effective dispersal in a natural population of the vertebrate-dispersed Neotropical tree, Simarouba amara (Simaroubaceae) by comparing the distances between maternal parents and their offspring and comparing gene movement via seed and pollen in the 50 ha plot of the Barro Colorado Island forest, Central Panama. In all cases (parent-pair, mother-offspring, father-offspring, sib-sib) distances between related pairs were significantly greater than distances to nearest possible neighbours within each category. Long-distance seedling establishment was frequent: 74% of assigned seedlings established > 100 m from the maternal parent [mean = 392 +/- 234.6 m (SD), range = 9.3-1000.5 m] and pollen-mediated gene flow was comparable to that of seed [mean = 345.0 +/- 157.7 m (SD), range 57.6-739.7 m]. For S. amara we found approximately a 10-fold difference between distances estimated by inverse modelling and mean seedling recruitment distances (39 m vs. 392 m). Our findings have important implications for future studies in forest demography and regeneration, with most seedlings establishing at distances far exceeding those demonstrated by negative density-dependent effects.

Conserving marginal populations of the foodplant (Impatiens noli-tangere) of an endangered moth (Eustroma reticulatum) in a changing climate

Impatiens noli-tangere is scarce in the UK and probably only native to the Lake District and Wales. It is the sole food plant for the endangered moth Eustroma reticulattum. Significant annual fluctuations in the size of I. noli-tangere populations endanger the continued presence of E. reticulatum in the UK. In this study, variation in population size was monitored across native populations of L noli-tangere in the English Lake District and Wales. In 1998, there was a crash in the population size of all metapopulations in the Lake District but not of those found in Wales. A molecular survey of the genetic affinities of samples in 1999 from both regions and a reference population from Switzerland was performed using AFLP and ISSR analyses. The consensus UPGMA dendrogram and a PCO scatter plot revealed clear differentiation between the populations of L noli-tangere in Wales and those in the Lake District. Most of the genetic variation in the UK (H-T= 0.064) was partitioned between (G(ST) = 0.455) rather than within (H-S = 0.034) regions, inferring little gene flow occurs between regions. There was similar bias towards differentiation between metapopulations in Wales, again consistent with low levels of interpopulation gene flow. This contrasts with far lower levels of differentiation in the Lake District which suggests modest rates of gene flow may occur between populations. It is concluded that in the event of local extinction of sites or populations, reintroductions should be restricted to samples collected from the same region. We then surveyed climatic variables to identify those most likely to cause local extinctions. Climatic correlates of population size were sought from two Lake District metapopulations situated close to a meteorological station. A combination of three climatic variables common to both sites explained 81-84% of the variation in plant number between 1990 and 2001. Projected trends for these climatic variables were used in a Monte Carlo simulation which suggested an increased risk of I. noli-tangere population crashes by 2050 at Coniston Water. but not at Derwentwater. Implications of these findings for practical conservation strategies are explored. (C) 2003 Elsevier Ltd. All rights reserved.

An approximate Bayesian computation approach to overcome biases that arise when using amplified fragment length polymorphism markers to study population structure

There is great interest in using amplified fragment length polymorphism (AFLP) markers because they are inexpensive and easy to produce. It is, therefore, possible to generate a large number of markers that have a wide coverage of species genotnes. Several statistical methods have been proposed to study the genetic structure using AFLP's but they assume Hardy-Weinberg equilibrium and do not estimate the inbreeding coefficient, F-IS. A Bayesian method has been proposed by Holsinger and colleagues that relaxes these simplifying assumptions but we have identified two sources of bias that can influence estimates based on these markers: (i) the use of a uniform prior on ancestral allele frequencies and (ii) the ascertainment bias of AFLP markers. We present a new Bayesian method that avoids these biases by using an implementation based on the approximate Bayesian computation (ABC) algorithm. This new method estimates population-specific F-IS and F-ST values and offers users the possibility of taking into account the criteria for selecting the markers that are used in the analyses. The software is available at our web site (http://www-leca.uif-grenoble.fi-/logiciels.htm). Finally, we provide advice on how to avoid the effects of ascertainment bias.