Local genetic structure in a clonal dioecious angiosperm

Conservation genomics of Agave tequilana Weber var. azul: low genetic differentiation and heterozygote excess in the tequila agave from Jalisco, Mexico

Background

Genetic diversity is fundamental for the survival of species. In particular, in a climate change scenario, it is crucial that populations maintain genetic diversity so they can adapt to novel environmental conditions. Genetic diversity in wild agaves is usually high, with low genetic differentiation among populations, in part maintained by the agave pollinators such as the nectarivorous bats. In cultivated agaves, patterns of genetic diversity vary according to the intensity of use, management, and domestication stage. In Agave tequilana Weber var. azul (A. tequilana thereafter), the plant used for tequila production, clonal propagation has been strongly encouraged. These practices may lead to a reduction in genetic diversity.

Methods

We studied the diversity patterns with genome-wide SNPs, using restriction site associated DNA sequencing in cultivated samples of A. tequilana from three sites of Jalisco, Mexico. For one locality, seeds were collected and germinated in a greenhouse. We compared the genomic diversity, levels of inbreeding, genetic differentiation, and connectivity among studied sites and between adults and juvenile plants.

Results

Agave tequilana presented a genomic diversity of H_T = 0.12. The observed heterozygosity was higher than the expected heterozygosity. Adults were more heterozygous than juveniles. This could be a consequence of heterosis or hybrid vigor. We found a shallow genetic structure (average paired F_ST = 0.0044). In the analysis of recent gene flow, we estimated an average migration rate among the different populations of m = 0.25. In particular, we found a population that was the primary source of gene flow and had greater genomic diversity (H_E and H_O ), so we propose that this population should continue to be monitored as a potential genetic reservoir.

Discussion

Our results may be the consequence of more traditional management in the studied specific region of Jalisco. Also, the exchange of seeds or propagules by producers and the existence of gene flow due to occasional sexual reproduction may play an important role in maintaining diversity in A. tequilana. For populations to resist pests, to continue evolving and reduce their risk of extinction under a climate change scenario, it is necessary to maintain genetic diversity. Under this premise we encourage to continue acting in conservation programs for this species and its pollinators.

2b-RAD Genotyping of the Seagrass Cymodocea nodosa Along a Latitudinal Cline Identifies Candidate Genes for Environmental Adaptation

Plant populations distributed along broad latitudinal gradients often show patterns of clinal variation in genotype and phenotype. Differences in photoperiod and temperature cues across latitudes influence major phenological events, such as timing of flowering or seed dormancy. Here, we used an array of 4,941 SNPs derived from 2b-RAD genotyping to characterize population differentiation and levels of genetic and genotypic diversity of three populations of the seagrass Cymodocea nodosa along a latitudinal gradient extending across the Atlantic-Mediterranean boundary (i.e., Gran Canaria—Canary Islands, Faro—Portugal, and Ebro Delta—Spain). Our main goal was to search for potential outlier loci that could underlie adaptive differentiation of populations across the latitudinal distribution of the species. We hypothesized that such polymorphisms could be related to variation in photoperiod-temperature regime occurring across latitudes. The three populations were clearly differentiated and exhibited diverse levels of clonality and genetic diversity. Cymodocea nodosa from the Mediterranean displayed the highest genotypic richness, while the Portuguese population had the highest clonality values. Gran Canaria exhibited the lowest genetic diversity (as observed heterozygosity). Nine SNPs were reliably identified as outliers across the three sites by two different methods (i.e., BayeScan and pcadapt), and three SNPs could be associated to specific protein-coding genes by screening available C. nodosa transcriptomes. Two SNPs-carrying contigs encoded for transcription factors, while the other one encoded for an enzyme specifically involved in the regulation of flowering time, namely Lysine-specific histone demethylase 1 homolog 2. When analyzing biological processes enriched within the whole dataset of outlier SNPs identified by at least one method, “regulation of transcription” and “signalling” were among the most represented. Our results highlight the fundamental importance signal integration and gene-regulatory networks, as well as epigenetic regulation via DNA (de)methylation, could have for enabling adaptation of seagrass populations along environmental gradients.

Partitioning resilience of a marine foundation species into resistance and recovery trajectories

The resilience of an ecological unit encompasses resistance during adverse conditions and the capacity to recover. We adopted a 'resistance-recovery' framework to experimentally partition the resilience of a foundation species (the seagrass Cymodocea nodosa). The shoot abundances of nine seagrass meadows were followed before, during and after simulated light reduction conditions. We determined the significance of ecological, environmental and genetic drivers on seagrass resistance (% of shoots retained during the light deprivation treatments) and recovery (duration from the end of the perturbed state back to initial conditions). To identify whether seagrass recovery was linearly related to prior resistance, we then established the connection between trajectories of resistance and recovery. Finally, we assessed whether recovery patterns were affected by biological drivers (production of sexual products-seeds-and asexual propagation) at the meadow-scale. Resistance to shading significantly increased with the genetic diversity of the meadow and seagrass recovery was conditioned by initial resistance during shading. A threshold in resistance (here, at a ca. 70% of shoot abundances retained during the light deprivation treatments) denoted a critical point that considerably delays seagrass recovery if overpassed. Seed densities, but not rhizome elongation rates, were higher in meadows that exhibited large resistance and quick recovery, which correlated positively with meadow genetic diversity. Our results highlight the critical role of resistance to a disturbance for persistence of a marine foundation species. Estimation of critical trade-offs between seagrass resistance and recovery is a promising field of research to better manage impacts on seagrass meadows.

Fine-scale genetic structure and flowering output of the seagrass Enhalus acoroides undergoing disturbance

Seagrass are under great stress in the tropical coast of Asia, where Enhalus acoroides is frequently the dominant species with a large food web. Here, we investigate the question of the fine-scale genetic structure of this ecologically important foundation species, subject to severe anthropogenic disturbance in China. The genetic structure will illuminate potential mechanisms for population dynamics and sustainability, which are critical for preservation of biodiversity and for decision-making in management and restoration. We evaluated the fine-scale spatial genetic structure (SGS) and flowering output of E. acoroides, and indirectly estimated the relative importance of sexual versus asexual reproduction for population persistence using spatial autocorrelation analysis. Results reveal high clonal diversity for this species, as predicted from its high sexual reproduction output. The stronger Sp statistic at the ramet-level compared with genet-level indicates that clonality increases the SGS pattern for E. acoroides. Significant SGS at the genet-level may be explained by the aggregated dispersal of seed/pollen cohorts. The estimated gene dispersal variance suggests that dispersal mediated by sexual reproduction is more important than clonal growth in this study area. The ongoing anthropogenic disturbance will negatively affect the mating pattern and the SGS patterns in the future due to massive death of shoots, and less frequency of sexual reproduction.

Carbon economy of Mediterranean seagrasses in response to thermal stress

Increased plant mortality in temperate seagrass populations has been recently observed after summer heatwaves, although the underlying causes of plant death are yet unknown. The potential energetic constrains resulting from anomalous thermal events could be the reason that triggered seagrass mortality, as demonstrated for benthic invertebrates. To test this hypothesis, the carbon balance of Posidonia oceanica and Cymodocea nodosa plants from contrasting thermal environments was investigated during a simulated heatwave, by analyzing their photosynthetic performance, carbon balance (ratio photosynthesis:respiration), carbohydrates content, growth and mortality. Both species were able to overcome and recover from the thermal stress produced by the six-week exposure to temperatures 4 °C above mean summer levels, albeit plants from cold waters were more sensitive to warming than plants from warm waters as reflected by their inability to maintain their P:R ratio unaltered. The strategies through which plants tend to preserve their energetic status varied depending on the biology of the species and the thermal origin of plants. These included respiratory homeostasis (P. oceanica warm-plants), carbon diversion from growth to respiration (C. nodosa cold-plants) or storage (P. oceanica warm-plants) and changes in biomass allocation (C. nodosa warm-plants). Findings suggest an important geographic heterogeneity in the overall response of Mediterranean seagrasses to warming with potential negative impacts on the functions and services offered by seagrass meadows including among others their capacity for carbon sequestration and carbon export to adjacent ecosystems.

Coping with potential bi-parental inbreeding: limited pollen and seed dispersal and large genets in the dioecious marine angiosperm Thalassia testudinum

The high prevalence of dioecy in marine angiosperms or seagrasses (>50% of all species) is thought to enforce cross-fertilization. However, seagrasses are clonal plants, and they may still be subject to sibling-mating or bi-parental inbreeding if the genetic neighborhood is smaller than the size of the genets. We tested this by determining the genetic neighborhoods of the dioecious seagrass Thalassia testudinum at two sites (Back-Reef and Mid-Lagoon) in Puerto Morelos Reef Lagoon, Mexico, by measuring dispersal of pollen and seeds in situ, and by fine-scale spatial autocorrelation analysis with eight polymorphic microsatellite DNA markers. Prevalence of inbreeding was verified by estimating pairwise kinship coefficients; and by analysing the genotypes of seedlings grown from seeds in mesocosms. Average dispersal of pollen was 0.3-1.6 m (max. 4.8 m) and of seeds was 0.3-0.4 m (max. 1.8 m), resulting in a neighborhood area of 7.4 m2 (range 3.4-11.4 m2) at Back-Reef and 1.9 (range 1.87-1.92 m2) at Mid-Lagoon. Neighborhood area (Na) derived from spatial autocorrelation was 0.1-20.5 m2 at Back-Reef and 0.1-16.9 m2 at Mid-Lagoon. Maximal extensions of the genets, in 19 × 30 m plots, were 19.2 m (median 7.5 m) and 10.8 m (median 4.8 m) at Back-Reef and Mid-Lagoon. There was no indication of deficit or excess of heterozygotes nor were coefficients of inbreeding (FIS) significant. The seedlings did not show statistically significant deficit of heterozygotes (except for 1 locus at Back-Reef). Contrary to our expectations, we did not find evidence of bi-parental inbreeding in this dioecious seagrass with large genets but small genetic neighborhoods. Proposed mechanisms to avoid bi-parental inbreeding are possible selection against homozygotes during fecundation or ovule development. Additionally, the genets grew highly dispersed (aggregation index Ac was 0.09 and 0.10 for Back-Reef and Mid-Lagoon, respectively); such highly dispersed guerrilla-like clonal growth form likely increases the probability of crossing between different potentially unrelated genets.

Complex spatial clonal structure in the macroalgae Fucus radicans with both sexual and asexual recruitment

In dioecious species with both sexual and asexual reproduction, the spatial distribution of individual clones affects the potential for sexual reproduction and local adaptation. The seaweed Fucus radicans, endemic to the Baltic Sea, has separate sexes, but new attached thalli may also form asexually. We mapped the spatial distribution of clones (multilocus genotypes, MLGs) over macrogeographic (>500 km) and microgeographic (<100 m) scales in the Baltic Sea to assess the relationship between clonal spatial structure, sexual recruitment, and the potential for natural selection. Sexual recruitment was predominant in some areas, while in others asexual recruitment dominated. Where clones of both sexes were locally intermingled, sexual recruitment was nevertheless low. In some highly clonal populations, the sex ratio was strongly skewed due to dominance of one or a few clones of the same sex. The two largest clones (one female and one male) were distributed over 100–550 km of coast and accompanied by small and local MLGs formed by somatic mutations and differing by 1–2 mutations from the large clones. Rare sexual events, occasional long‐distance migration, and somatic mutations contribute new genotypic variation potentially available to natural selection. However, dominance of a few very large (and presumably old) clones over extensive spatial and temporal scales suggested that either these have superior traits or natural selection has only been marginally involved in the structuring of genotypes.

Causes and consequences of contrasting genetic structure in sympatrically growing and closely related species

Gene flow, natural selection and genetic drift are processes that play a major role in shaping the genetic structure of natural populations. In addition, genetic structures of individual populations are strongly correlated with their geographical position within the species distribution area. The highest levels of genetic variation are usually found in the centre of a species' distribution and tend to decrease beyond that point. Additionally, narrowly endemic taxa are expected to be characterized by lower levels of genetic variation than their widespread congeners. To understand the historical circumstances that shape populations of sympatric and closely related taxa, microsatellite markers were used, while populations of the three closely related and sympatric Mediterranean Salvia species (S. officinalis L., S. fruticosa Mill. and S. brachyodon Vandas) served as a study model. In the populations of widespread S. officinalis, located in the central parts of this species' distribution area, no population genetic disturbances were detected. The narrow endemic S. brachyodon showed heterozygote excess, clonal reproduction and a genetic bottleneck. Because the genetic bottleneck was likely caused by the disappearance of suitable open-type habitats, the recent wildfire that cleared the terrain probably saved the S. brachyodon population from gradual deterioration and extinction. At the same time, clonal reproduction could serve as a valuable mechanism in the preservation of genetic variability. The results of the disjunct S. fruticosa population indicated heterozygote deficiency, inbreeding, hybridization with S. officinalis and population expansion. The hybridization with S. officinalis along with the abandonment of the agro-pastoral system are likely the main drivers of the strong expansion of S. fruticosa in the studied location. As many relevant findings and conclusions regarding historical and contemporary demography of individual populations or species can be reached only through their comparison with closely related taxa, this study demonstrates the importance and advantages of such a multi-species approach.

Is the degree of clonality of forest herbs dependent on gap age? Using fingerprinting approaches to assess optimum successional stages for montane forest herbs

Using molecular fingerprinting (amplified fragment length polymorphism [AFLP] method), we explored the potential of small-scale population analysis for understanding colonization patterns of herb layer species in forests after canopy disturbance. We investigated three common forest understorey species with different life forms (Trientalis europaea, Calamagrostis villosa, and Vaccinium myrtillus) in the Harz Mountains in Germany in three different gap age classes and undisturbed forest. For two of them (T. europaea and C. villosa), we analyzed clone sizes and clonal structure. We hypothesized that clone sizes depend on age since gap formation and are affected by light availability. Mean patch sizes of V. myrtillus, T. europaea, and C. villosa formed were 3.7 m², 27.9 m², and 40.6 m², respectively. Trientalis europaea and C. villosa patches consisted mostly of more than one genet. Largest clone sizes of T. europaea were encountered in gaps of intermediate successional age (15–60 years, averaged minimum estimation of clone sizes: 6.56 m²) whereas clone size of C. villosa was found to be independent from gap age and had a mean minimum clone size of 0.49 m². In both species, clone size was positively related to light availability. Additionally, there was a positive relationship between clone size and ramet density for T. europaea and C. villosa. Genetic variation was higher within populations of T. europaea and C. villosa than among populations. Trientalis europaea was the only species with a clear genetic isolation by distance, pointing at an equilibrium between gene flow and genetic drift. In conclusion, we showed that forest canopy gap dynamics clearly affect the small-scale structure of populations of understorey plants. Species with high lateral growth rates, such as T. europaea offer the possibility to serve as “ecological clock” for dating ecological processes.

Genetic variability of an endangered Bromeliaceae species (Pitcairnia albiflos) from the Brazilian Atlantic rainforest

Pitcairnia albiflos is a Bromeliaceae species endemic to Brazil that has been included as data-deficient in the extinction risk list of Brazilian flora. We analyzed genetic variability in P. albiflos populations using RAPD markers to investigate population structure and reproductive mechanisms and also to evaluate the actual extinction risk level of this species. Leaves of 56 individuals of P. albiflos from three populations were collected: Urca Hill (UH, 20 individuals), Chacrinha State Park (CSP, 24 individuals) and Tijuca National Park (TNP, 12 individuals). The RAPD technique was effective in characterizing the genetic diversity in the P. albiflos populations since it was possible to differentiate the populations and to identify exclusive bands for at least two of them. Even if there is low genetic diversity among them (CSP-UH = 0.463; CSP-TNP = 0.440; UH-TNP = 0.524), the populations seem to be isolated according to the low genetic diversity observed within them (H(pop) CSP = 0.060; H(pop) UH = 0.042; H(pop) TNP = 0.130). This fact might be the result of clonal and self-reproduction predominance and also from environmental degradation around the collection areas. Consequently, it would be important to protect all populations both in situ and ex situ to prevent the decrease of genetic variability. The low genetic variability among individuals of the same population confirms the inclusion of this species as critically endangered in the risk list for Brazilian flora.

Clonal and spatial genetic structure within populations of a coastal plant, Carex kobomugi (Cyperaceae)

Clarification of clonal growth pattern is critical for understanding the population dynamics and reproductive system evolution of clonal plant species. The contribution of clonality to the spatial genetic structure (SGS) within populations is also an important issue. I examined the spatial distribution of genetic variability within two populations of the coastal plant Carex kobomugi using seven microsatellite loci. Genotyping of 226 and 140 ramets within 14 × 40 m and 14 × 34 m plots on two populations revealed 36 and 33 multilocus genotypes, respectively. To quantify the extent of intermingling among clones, for each genet, I calculated the dominance of ramets belonging to a particular genet within a spatial range of the genet. Furthermore, I analyzed spatial distribution of genotypes within 2 × 2 m and 1 × 2 m quadrats using second-order spatial statistics. These analyses indicated that clones are highly intermingled, suggesting a low level of spatial interaction among clones. Spatial autocorrelation analysis of kinship coefficient including all pairs of ramets showed significantly stronger SGS than analysis considering only pairs between different genets. I conclude that clonal propagation largely contributes to SGS at a fine scale.

Dynamics of distribution and performance of ramets constructing genets: a demographic-genetic study in a clonal plant, Convallaria keiskei

BACKGROUND AND AIMS

In clonal plants producing vegetative offspring, performance at the genet level as well as at the ramet level should be investigated in order to understand the entire picture of the population dynamics and the life history characteristics. In this study, demography, including reproduction and survival, the growth patterns and the spatial distributions of ramets within genets of the clonal herb Convallaria keiskei were explored.

METHODS

Vegetative growth, flowering and survival of shoots whose genets were identified using microsatellite markers were monitored in four study plots for 3 years (2003-2005). The size structures of ramets in genets and their temporal shifts were then analysed. Their spatial distributions were also examined.

KEY RESULTS

During the census, 274 and 149 ramets were mapped in two 1 x 2 m plots, and 83 and 94 ramets in two 2 x 2 m quadrats. Thirty-eight genotypes were identified from 580 samples. Each plot included 5-18 genets, and most ramets belonged to the predominant genet(s) in each plot. Shoots foliated yearly for several years, but flowering ramets did not have an inflorescence the next year. A considerable number of new clonal offspring persistently appeared, forming a bell-shaped curve of the size structure of ramets in each genet. Comparing the structures modelled by the normal distributions suggested variation among ramets belonging to a single genet and variation among genets. Furthermore, spatial analyses revealed clumped and distant distributions of ramet pairs in a genet, in which the distant patterns corresponded to the linearly elongating clonal growth pattern of this species.

CONCLUSION

Characteristics of ramet performances such as flowering and recruitment of clonal offspring, in addition to growth, played a large part in the regulation of genet dynamics and distribution, which were different among the studied genets. These might be characteristics particularly relevant to clonal life histories.

Population genetic structure of native versus naturalized sympatric shrub willows (Salix; Salicaceae)

Vegetative propagation of an introduced species can contribute significantly to its ability to spread and become naturalized, potentially in competition with native species. This study focused on the naturalization of a willow shrub, Salix purpurea, which was introduced to the United States from Europe and is commonly sympatric with the native shrub willow, S. eriocephala. Both species are capable of vegetative and sexual reproduction, but little is known about their relative frequency, nor the impact of clonal propagation on population-level genetic diversity. We analyzed genotypes at several microsatellite loci in 993 individuals belonging to 30 subpopulations of S. eriocephala and 28 subpopulations of S. purpurea in areas of sympatry across three watersheds to compare their genetic diversity and genetic structure. Our results revealed six subpopulations of S. purpurea containing plants with identical multilocus genotypes, while clonal individuals were rare among S. eriocephala populations. These species are dioecious with relatively high levels of heterozygosity, but S. eriocephala had much higher allelic diversity and genotypic diversity than did S. purpurea. These results strongly suggest that vegetative propagation has contributed to the naturalization of S. purpurea and has resulted in higher levels of genetic differentiation among S. purpurea populations than among native S. eriocephala populations.

Clonal architecture in an intertidal bed of the dwarf eelgrass Zostera noltii in the Northern Wadden Sea: persistence through extreme physical perturbation and the importance of a seed bank

Genotypic structure and temporal dynamics of the dwarf seagrass, Zostera noltii, were studied in an intertidal meadow that has persisted since prior to 1936 near the Wadden Sea island of Sylt. Samples were collected from two 10 × 10 m plots separated by 250 m from May 2002 to June 2005 and from four 1 × 1 m plots from June 2003 to September 2004. All the samples were genotyped with nine microsatellite loci. No genotypes were shared between the plots separated by 250 m. Genetic diversity was higher in the Wadden Sea than in the other regions of its geographic range. The average clone size (genets) (SD) in the two plots was 1.38 (0.26) and 1.46 (0.4) m², respectively, with a range up to 9 m² and <20% persisted for >4 years. A high genetic and genotypic diversity was maintained by annual recruitment of seedlings despite a dramatic decrease in ramet density that coincided with the severe heat stress event of 2003. Fine-scale (1 m²) analysis suggested that extensive loss of seagrass cover precluded space competition among the genets, while a persistent seed bank prevented local extinction. Long-term persistence of Z. noltii meadows in the intertidal Wadden Sea was achieved by high genet turnover and frequent seedling recruitment from a seed bank, in contrast to the low diversity observed in large and long-living clones of Z. noltii and other seagrasses in subtidal habitats.

Clonal diversity of Clintonia udensis Trautv. et Mey. populations and its correlation with ecological factors

The clonal diversity of Clintonia udensis Trautv. et Mey. was detected by ISSR markers among 16 populations, and its correlation with ecological factors was analyzed as well in this work. Results showed that individuals (clonal ramets) per genotype were 1.12 and 1.149 at population and species levels, respectively, and that the 16 populations were all multiclonal. The detected genotypes were localized, without exception, within populations but demonstrated relatively high clonal differentiation among populations. The clonal diversity of the studied populations was high, with the average Simpson's index of 0.975, while the genets showed a clonal architecture of "guerilla". The population genetic diversities revealed by genet were consistent with those by ramet, further confirming their genetic differentiation among populations. And its genotype diversity within populations probably resulted largely from the frequent seedling regeneration and self-compatibility. In addition, the correlation analysis further revealed that, among the ecological factors, Simpson's index of C. udensis had a significant positive correlation (P<0.05) with pH values in the soil but not others.

Clonal extent, apical dominance and networking features in the phalanx angiosperm Zostera noltii Hornem

Disaggregating seagrass meadows and studying its components separately (clones, ramets, shoots) can provide us insights on meadow dynamics and growth patterns. The clonal growth, dependent upon clonal rules may regulate and impose constraints to plant architecture and, therefore, determine how individual clones evolve into the environment. In order to investigate the relationship between clonal growth rules and clone architecture, the belowground network architecture of single-clones of the seagrass Zostera noltii was studied. Networks were traced in situ after washing out the overlying sediment, and network characteristics were measured using digital analysis: area covered by clone, total rhizome length, type of rhizomatic axes (main, secondary, tertiary, quaternary), number and length of the internodes, branching angles and branching frequencies. This approach revealed that Z. noltii is able to develop into large clones integrating up to 300 internodes, 676 cm of rhizome, 208 shoots and 4,300 cm2 of plant area. Internodal length depended on both, the distance to the apical shoot (time effect) and the axes type (apical dominance effect). However, average branching angle was independent of axis type (average 58.3 ± 0.75), but varied significantly depending on the distance from the apical shoot. This average branching angle allows Z. noltii maximize the rate of centrifugal expansion, maintaining a high density in colonized areas to produce close stands but also minimizing the investment in belowground biomass and ramets overlapping. The clonal architecture of Z. noltii seems to be regulated by the interaction of both, apical dominance strength and clonal integration distance. Moreover, clonal growth rules and growth pattern seem to constrain clonality through (clonal) plant architecture regulations (i.e. branching is restricted in secondary axes, similar average branching angles regardless the axes, the higher the distance to the apex the higher the number of internodes in secondary axes, shorter internodes in secondary and tertiary axes). Future research efforts should focus on how these complex relationships between apical dominance and clonal integration interact to elucidate the temporal (seasonal) and spatial scales of both processes and the outcome at the plant architectural level.

Heterozygote excess in a self-incompatible and partially clonal forest tree species -- Prunus avium L

Wild cherry (Prunus avium L.), a partially asexual self-incompatible forest tree, shows heterozygote excess, which is a poorly studied phenomenon. In three natural populations, we found significant heterozygote excess at almost all investigated loci (eight microsatellites and markers for the self-incompatibility locus). We examined four hypotheses to account for this observed heterozygote excess. First, negative F(IS) can result from a lack of selfed progeny in small populations of outcrossing species. A second explanation for negative F(IS) is selection during the life cycle of the most heterozygous individuals. A third explanation is negative assortative mating when reproduction occurs between individuals bearing phenotypes more dissimilar than by chance. The last explanation for negative F(IS) relies on asexual reproduction. Expectations for each hypothesis were tested using empirical data. Patterns of F(IS) differed among loci. Nevertheless, our experimental results did not confirm the small sample size hypothesis. Although one locus is probably under a hitch-hiking effect from the SI locus, we rejected the effect of the self-incompatibility locus for the genome as a whole. Similarly, although one locus showed a clear pattern consistent with the selection of heterozygous individuals, the heterosis effect over the whole genome was rejected. Finally, our results revealed that clonality probably explains significant negative F(IS) in wild cherry populations when considering all individuals. More theoretical effort is needed to develop expectations and hypotheses, and test them in the case of species combining self-incompatibility and partially asexual reproduction.