The use of neutral and non-neutral SSRs to analyse the genetic structure of a Tunisian collection of Medicago truncatula lines and to reveal associations with eco-environmental variables.

In this study, we investigated the genetic diversity of a collection of 136 Medicago truncatula lines from 10 Tunisian natural populations collected in well-defined locations and in various ecological conditions of soil, salinity and water availability. The genetic diversity was evaluated using a set of 18 microsatellites (SSRs), representing the 8 chromosomes of M. truncatula. A neutrality test showed that 7 SSRs were non-neutral with evidence of balancing selection. The 11 neutral SSRs revealed a geographical pooling with the Tunisian Dorsale axis restricting migration of alleles. The 7 non-neutral alleles demonstrate a correlation with rainfall, altitude and salinity environmental variables suggesting that these SSRs are linked to genes involved in water use efficiency, resistance to salinity or adaptation to altitude, and that there is local adaptation of M. truncatula to these variables. This demonstrates that the choice of so-called neutral markers should be carefully evaluated in population genetic studies. This study illustrates the genetic diversity occurring in natural Tunisian populations of M. truncatula and describes the first collection of this species dedicated to natural variation involved in adaptation to the environment.

How multilocus genotypic pattern helps to understand the history of selfing populations: a case study in Medicago truncatula.

The occurrence of populations exhibiting high genetic diversity in predominantly selfing species remains a puzzling question, since under regular selfing genetic diversity is expected to be depleted at a faster rate than under outcrossing. Fine-scale population genetics approaches may help to answer this question. Here we study a natural population of the legume Medicago truncatula in which both the fine-scale spatial structure and the selfing rate are characterized using three different methods. Selfing rate estimates were very high ( approximately 99%) irrespective of the method used. A clear pattern of isolation by distance reflecting small seed dispersal distances was detected. Combining genotypic data over loci, we could define 34 multilocus genotypes. Among those, six highly inbred genotypes (lines) represented more than 75% of the individuals studied and harboured all the allelic variation present in the population. We also detected a large set of multilocus genotypes resembling recombinant inbred lines between the most frequent lines occurring in the population. This finding illustrates the importance of rare recombination in redistributing available allelic diversity into new genotypic combinations. This study shows how multilocus and fine-scale spatial analyses may help to understand the population history of self-fertilizing species, especially to make inferences about the relative role of foundation/migration and recombination events in such populations.

Analysis of population structure in autotetraploid species.

Population structure parameters commonly used for diploid species are reexamined for the particular case of tetrasomic inheritance (autotetraploid species). Recurrence equations that describe the evolution of identity probabilities for neutral genes in an "island model" of population structure are derived assuming tetrasomic inheritance. The expected equilibrium value of FST is computed. In contrast to diploids, the correlation of genes between individuals within populations with respect to genes between populations (FST) may vary among loci due to the particular segregation patterns expected under tetrasomic inheritance and is consequently inappropriate for estimating demographic parameters in such populations. We thus define a new parameter (rho) and derive its relationship with Nm. This relationship is shown to be independent from both the selfing rate and the proportion of double reduction. Finally, the statistical procedure required to evaluate these parameters using data on gene frequencies distribution among autotetraploid populations is developed.

Inbreeding depression in small populations of self-incompatible plants.

Self-incompatibility (SI) is a widespread mechanism that prevents inbreeding in flowering plants. In many species, SI is controlled by a single locus (the S locus) where numerous alleles are maintained by negative frequency-dependent selection. Inbreeding depression, the decline in fitness of selfed individuals compared to outcrossed ones, is an essential factor in the evolution of SI systems. Conversely, breeding systems influence levels of inbreeding depression. Little is known about the joint effect of SI and drift on inbreeding depression. Here we studied, using a two-locus model, the effect of SI (frequency-dependent selection) on a locus subject to recurrent deleterious mutations causing inbreeding depression. Simulations were performed to assess the effect of population size and linkage between the two loci on the level of inbreeding depression and genetic load. We show that the sheltering of deleterious alleles linked to the S locus strengthens inbreeding depression in small populations. We discuss the implications of our results for the evolution of SI systems.

Mutation load depending on variance in reproductive success and mating system.

Equalization of reproductive success of individuals, although it results in an increase of effective population size, leads also to an increase of the mutation load. The magnitude of this increase depends highly on the mode of fitness interactions between deleterious mutations, and is higher in the case of inbreeding. Recommended practices in conservation genetics must be evaluated in regards to these differing consequences of an increase of effective population size. To keep a balance between retaining genetic variability and minimizing the increase of the mutation load, equalization of reproductive success of a set of individuals rather than of every individual might be more advantageous.

AMaCAID: a useful tool for Accurate Marker Choice for Accession Identification and Discrimination.

AMaCAID is an R program designed to analyse multilocus genotypic patterns in large samples. It allows (i) the computation of the number and frequency of the different multilocus patterns available in a molecular data set and (ii) the analysis of discriminatory power of each combination of k markers among n available. It thus enables the identification of the minimum number of markers required to distinguish all the observed genotypes and the subset of markers that maximize the number of distinct genotypes. AMaCAID can be used with any kind of molecular markers, on data sets mixing different kinds of markers, but also on qualitative characters like morphological or taxonomic traits. AMaCAID has been built primarily to select subsets of markers for identifying accessions and monitoring their genetic stability during regeneration cycles in an ex situ genebank. It can, however, also be used to screen any kind of data set that characterizes a set of individuals or species (e.g. taxonomic or phylogenetic studies) for discrimination purposes. The size of the assayed sample has no limitation, but the program only performs computations on all combinations of markers when there are less than 25 markers. For larger number of markers/characters, it is possible to ask AMaCAID to screen a large but limited number of combinations of markers. We apply AMaCAID to three data sets involving either molecular or taxonomic data and give some results on the computing time of the program with respect to the size of the data set.

Effective population size associated with self-fertilization: lessons from temporal changes in allele frequencies in the selfing annual Medicago truncatula.

Despite its significance in evolutionary and conservation biology, few estimates of effective population size (N(e)) are available in plant species. Self-fertilization is expected to affect N(e), through both its effect on homozygosity and population dynamics. Here, we estimated N(e) using temporal variation in allele frequencies for two contrasted populations of the selfing annual Medicago truncatula: a large and continuous population and a subdivided population. Estimated N(e) values were around 5-10% of the population census size suggesting that other factors than selfing must contribute to variation in allele frequencies. Further comparisons between monolocus allelic variation and changes in the multilocus genotypic composition of the populations show that the local dynamics of inbred lines can play an important role in the fluctuations of allele frequencies. Finally, comparing N(e) estimates and levels of genetic variation suggest that H(e) is a poor estimator of the contemporaneous variance effective population size.

Marker-based investigation of inbreeding depression in the endangered species Brassica insularis.

Various methods have been proposed to estimate inbreeding depression and to assess its consequences for natural populations. As an alternative to controlled crosses, the use of molecular markers has allowed direct investigation of inbreeding depression in natural populations, but usually suffers from low statistical power. Here, we investigated the effect of inbreeding depression on survival in two populations of the rare species Brassica insularis, using both controlled crosses and a marker-based approach. We compare the respective merits of the two approaches for studying inbreeding depression. We also use information from the molecular markers to dissect in detail patterns of inbreeding depression in this species. A posteriori, we find that combining the approaches was not necessary to obtain simple point estimates of inbreeding depression. However, using molecular markers may give insight into the genetic basis of inbreeding depression, such as the occurrence of epistatic interactions among deleterious alleles or purging.

Domestication history in the Medicago sativa species complex: inferences from nuclear sequence polymorphism.

DNA sequence polymorphism carries genealogical information and allows for testing hypotheses on selection and population history, especially through coalescent-based analysis. Understanding the evolutionary forces at work in plant domestication and subsequent selection is of critical importance for the management of genetic resources. In this study, we surveyed DNA sequence diversity at two assumed neutral nuclear loci in the wild-domesticated species complex of alfalfa (Medicago sativa L.). A high level of polymorphism was detected. The domesticated pool contains on average 31% less diversity than the wild pool, but with a high heterogeneity among loci. Coalescent simulations of the domestication process showed that this result cannot be explained by assuming a constant population size but is rather consistent with a demographic bottleneck during domestication. A very low level of divergence was detected between the wild and the domesticated forms as well as between the related subspecies of the M. sativa species complex. However, the originality of the Spanish wild populations, already observed based on mitochondrial DNA polymorphism, was confirmed. These results, together with patterns of intrapopulation polymorphism, suggest that nuclear sequence polymorphism could be a promising tool, complementary to mitochondrial DNA and phenotypic evaluations, to investigate historical demographic and evolutionary processes.

Spatial structure of genetic variation and primary succession in the pioneer tree species Antirhea borbonica on La Réunion.

In habitats where colonization and extinction are recurrent, the distribution of gene frequencies among patches of suitable habitat may reflect the age structure of different populations. In this study, we quantify population genetic structure for a pioneer tree species, Antirhea borbonica, in a chrono-sequence of primary succession on the lava flows of the Piton de La Fournaise volcano (La Réunion). Using microsatellite loci and amplified fragment length polymorphism (AFLP) markers, we quantified genetic variation within and among populations for early- and late-succession populations in a landscape where extinction and recolonization are recurrent (the 'Grand Brûlé') and for late-succession populations in a more stable landscape. This study produced three main results. First, we detected no evidence that founder events increase genetic differentiation among colonizing populations; F(ST) values among early- and among late-succession populations were similar. Second, we found no evidence for isolation by distance; genetic distance was not correlated with spatial distance within and among populations. Third, F(IS) values are consistently high in all populations, despite the fact that A. borbonica populations are functionally close to dioecy and thus expected to have an outcrossing mating system. Multiple colonization events from different sources may limit differentiation among young populations and spatial isolation may enhance differentiation among late-succession populations. Ecological processes acting during colonization may create the conditions for spatial aggregation within pioneer populations, and thus contribute to the high F(IS) values.

Spatial effects and rare outcrossing events in Medicago truncatula (Fabaceae).

In order to elucidate the mechanisms underlying the large amount of RAPD polymorphism found in 1990 in a population of the selfing annual Medicago truncatula GAERTN. (Fabaceae), we have analysed most of the individuals (n = 363) from the same population 6 years later using microsatellite loci. We confirm the result of the earlier study, namely that this population is very polymorphic and highly subdivided, with approximately 37% of the variance distributed among subpopulations, only 50 m apart one from another. We use standard F-statistics analyses, linkage disequilibria, minimum spanning network, multilocus assignment tests and spatial autocorrelation analyses to test the hypotheses that spatial structure and outcrossing events are involved in maintaining the large amount of genetic diversity at the level of each subpopulation. Interestingly, fine-scale spatial structure could be observed in only one subpopulation suggesting that other mechanisms are acting elsewhere. To the best of our knowledge, this is the first study of fine spatial genetic structure in a predominantly selfing species.

Differentiation between natural and cultivated populations of Medicago sativa (Leguminosae) from Spain: analysis with random amplified polymorphic DNA (RAPD) markers and comparison to allozymes.

The conservation of a crop's wild relatives as genetic resources requires an understanding of the way genetic diversity is maintained in their populations, notably the effect of crop-to-wild gene flow. In this study, the amount of differentiation between natural and cultivated populations of Medicago sativa was analysed using random amplified polymorphic DNA (RAPD) markers and an extension of the AMOVA procedure adapted to autotetraploid organisms. Simulations of structured populations were performed to test whether AMOVA provides estimates of population structure in autotetraploids that can be directly compared to those obtained for allozyme data. Simulations showed that straight phi-statistics allow a good estimation of population differentiation when unbiased allelic frequencies are used to correct the conditional expectations of squared genetic distances. But such unbiased estimates can not be practically guaranteed, and population structure is notably overestimated when some populations are fixed for the presence of amplified fragments. However, removing fixed loci from the data set improves the statistical power of the test for population structure. The genetic variation of 15 natural and six cultivated populations of M. sativa was analysed at 25 RAPD loci and compared to estimates computed with allozymes on the same set of populations. Although RAPD markers revealed less within-population genetic diversity than allozymes, the quantitative and qualitative patterns of population structure were in full agreement with allozymes. This confirmed the conclusions drawn from the allozymic survey: crop-to-wild gene flow occurred in many locations, but some other mechanisms opposed cultivated traits to be maintained into natural populations.

Mitochondrial DNA diversity and male sterility in natural populations of Daucus carota ssp carota.

Mitochondrial variability was investigated in natural populations of wild carrot (Daucus carota ssp carota) in different regions: South of France, Greece, and various sites in the Mediterranean Basin and Asia. Total DNA was digested with two restriction endonucleases (EcoRV and HindIII) and probed with three mitochondrial DMA-specific genes (coxI, atp6, and coxII). Twenty-five different mitochondrial types were found in 80 analyzed individuals. Thirteen mitotypes were found among the 7 French populations studied. On average, 4.4 different mitotypes were observed per population, and these mitotypes were well-distributed among the populations. All of the mitochondrial types were specific to a single region. However, the proportion of shared restriction fragments between 2 mitotypes from different regions was not particularly lower than that which occurred among mitotypes from a single region. On the basis of the sexual phenotype [male-sterile (MS) or hermaphrodite] of the plants studied in situ and that of their progeny, 2 mitotypes were found to be highly associated with male sterility. Eighty percent of the plants bearing these mitotypes were MS in situ, and all of these plants produced more than 30% MS plants in their progeny. This association with male sterility was consistent in several populations, suggesting an association with a cytoplasmic male-sterility system. Moreover, these two mitotypes had very similar mitochondrial DNA restriction patterns and were well-differentiated from the other mitotypes observed in wild plants and also from those observed in the two CMS types already known in the cultivated carrot. This suggests that they correspond to a third cytoplasmic sterility.

Evidence for gene flow between wild and cultivated Medicago sativa (Leguminosae)based on allozyme markers andquantitative traits.

Genetic differentiation between co-occurring crops and their wild relatives will be greatly modified by crop-to-weed gene flow and variation between human and natural selective pressures. The maintenance of original morphological features in most natural populations of Medicago sativa in Spain questions the relative extent of these antagonistic forces. In this paper, we measured and compared the pattern of population differentiation within and among the wild and cultivated gene pool with respect to both allozymes and quantitative traits. Patterns of diversity defined three kinds of natural populations. First, some populations were intermediate with respect to both allozymes and quantitative traits. This suggests that crop-to-weed gene flow may have created hybrid populations in some locations. Second, some populations were different from all the cultivated landraces with respect to both allozymes and quantitative traits. This probably results from variable gene flow in space and in time, due to demographic stochasticity in either natural or cultivated populations. Third, differentiation from cultivated landraces was only achieved for the quantitative traits but not for allozymes in two populations. This suggests that natural selection in some locations may oppose gene flow to establish cultivated traits into the natural introgressed populations.

How mitochondrial DNA diversity can help to understand the dynamics of wild-cultivated complexes. The case of Medicago sativa in Spain.

In order to clarify the relationships (genetic exchange and shared ancestry) between natural and cultivated populations of alfalfa (Medicago sativa L.) in Spain, we investigated the patterns of mitochondrial DNA variation (characterized through restriction fragment length polymorphism) for 248 individuals in seven natural and six cultivated populations of this species. Mitochondrial variation was evidenced in both natural and cultivated populations of M. sativa. Among the seven mitotypes identified in the species, two were specific of the natural populations, a result attesting the fact that the Spanish wild form of M. sativa is an original genetic pool compared to the cultivated one. Other mitotypes were observed in both natural and cultivated populations, suggesting the occurrence of gene flow through seeds from cultivated towards natural populations. Comparisons with previously gathered nuclear and phenotypic data give insights into the different evolutionary forces acting on the different kinds of Spanish natural populations examined so far.

Inferences from mitochondrial DNA patterns on the domestication history of alfalfa (Medicago sativa).

The evolutionary history of a wild-cultivated complex is the outcome of different factors, among which are the domestication of the cultivated form, its geographical and demographic expansion, as well as gene flow between natural and cultivated populations. To clarify this history for alfalfa (Medicago sativa), we adopted a phylogeographical approach based on a sample of natural and cultivated populations of the M. sativa species complex. This sample was characterized for mtDNA variation through restriction fragment length polymorphisms. Twenty-two mitotypes were identified in the whole data set (155 individuals). In the wild pool from the presumed area of origin (Near East to Central Asia), the diversity was high but the absence of geographical differentiation hinders a more precise location of the centre of domestication. Within the cultivated alfalfa, the geographical structure suggests strongly the existence of at least two independent routes of dissemination of alfalfa from its centre of origin. In particular, original mitotypes detected in some regions suggest that wild populations not sampled in this study, and located outside the presumed area of domestication, may have contributed to the diversity of the cultivated pool. Finally, analysis of the mtDNA variation shows that natural populations collected in Spain constitute an endemic wild pool now introgressed partially by cultivated alfalfa.