Genetic differentiation and habitat preference of flowering-time variants within Gymnadenia conopsea

Historical biogeography and local adaptation explain population genetic structure in a widespread terrestrial orchid

BACKGROUND AND AIMS

Historical changes in environmental conditions and colonization-extinction dynamics have a direct impact on the genetic structure of plant populations. However, understanding how past environmental conditions influenced the evolution of species with high gene flow is challenging when signals for genetic isolation and adaptation are swamped by gene flow. We investigated the spatial distribution and genetic structure of the widespread terrestrial orchid Epipactis helleborine to identify glacial refugia, characterize postglacial population dynamics and assess its adaptive potential.

METHODS

Ecological niche modelling was used to locate possible glacial refugia and postglacial recolonization opportunities of E. helleborine. A large single-nucleotide polymorphism (SNP) dataset obtained through genotyping by sequencing was used to define population genetic diversity and structure and to identify sources of postglacial gene flow. Outlier analyses were used to elucidate how adaptation to the local environment contributed to population divergence.

KEY RESULTS

The distribution of climatically suitable areas was restricted during the Last Glacial Maximum to the Mediterranean, south-western Europe and small areas in the Alps and Carpathians. Within-population genetic diversity was high in E. helleborine (mean expected heterozygosity, 0.373 ± 0.006; observed heterozygosity, 0.571 ± 0.012; allelic richness, 1.387 ± 0.007). Italy and central Europe are likely to have acted as important genetic sources during postglacial recolonization. Adaptive SNPs were associated with temperature, elevation and precipitation.

CONCLUSIONS

Forests in the Mediterranean and Carpathians are likely to have acted as glacial refugia for Epipactis helleborine. Postglacial migration northwards and to higher elevations resulted in the dispersal and diversification of E. helleborine in central Europe and Italy, and to geographical isolation and divergent adaptation in Greek and Italian populations. Distinguishing adaptive from neutral genetic diversity allowed us to conclude that E. helleborine has a high adaptive potential to climate change and demonstrates that signals of adaptation and historical isolation can be identified even in species with high gene flow.

Extensive pollinator sharing does not promote character displacement in two orchid congeners

Pollinator sharing between close relatives can be costly and can promote pollination niche partitioning and floral divergence. This should be reflected by a higher species divergence in sympatry than in allopatry. We tested this hypothesis in two orchid congeners with overlapping distributions and flowering times. We characterized floral traits and pollination niches and quantified pollen limitation in 15 pure and mixed populations, and we measured phenotypic selection on floral traits and performed controlled crosses in one mixed site. Most floral traits differed between species, yet pollinator sharing was extensive. Only the timing of scent emission diverged more in mixed sites than in pure sites, and this was not mirrored by the timing of pollinator visitation. We did not detect divergent selection on floral traits. Seed production was pollen limited in most populations but not more severely in mixed sites than in pure sites. Interspecific crosses produced the same or a higher proportion of viable seeds than intraspecific crosses. The two orchid species attract the same pollinator species despite showing divergent floral traits. However, this does not promote character displacement, implying a low cost of pollinator sharing. Our results highlight the importance of characterizing both traits and ecological niches in character displacement studies.

Genome size variability of the population of Gymnadenia conopsea (Orchidaceae) in the Novosibirsk Region (Russia)

Gymnadenia conopsea (Fragrant orchid) is a widespread model species for the study of polymorphism in the European part of Russia and abroad. Flow cytometry refers to rapid methods effective for assessing the structure of a population by DNA size estimation. However, this method is only available in habitats with close proximity to the place of the experiment, because this indicator is determined only in fresh leaves. A small population of G. conopsea in the industrial area adjacent to the Novosibirsk urban agglomeration in the Novosibirsk Region (NR) was studied. Only diploids were detected in the G. conopsea cenopopulation occurred in the industrial zone of the NR. In this study we pointed on significance of genome size variation, which correlates with the adaptive evolution of the G. conopsea in West Siberia. The revealed absence of minor cytotypes (individuals with a higher level of ploidy) in G. conopsea population reduces intraspecific and intrapopulation diversity of the species.

Divergent selection on flowering phenology but not on floral morphology between two closely related orchids

Closely related species often differ in traits that influence reproductive success, suggesting that divergent selection on such traits contribute to the maintenance of species boundaries. Gymnadenia conopsea ss. and Gymnadenia densiflora are two closely related, perennial orchid species that differ in (a) floral traits important for pollination, including flowering phenology, floral display, and spur length, and (b) dominant pollinators. If plant-pollinator interactions contribute to the maintenance of trait differences between these two taxa, we expect current divergent selection on flowering phenology and floral morphology between the two species. We quantified phenotypic selection via female fitness in one year on flowering start, three floral display traits (plant height, number of flowers, and corolla size) and spur length, in six populations of G. conopsea s.s. and in four populations of G. densiflora. There was indication of divergent selection on flowering start in the expected direction, with selection for earlier flowering in two populations of the early-flowering G. conopsea s.s. and for later flowering in one population of the late-flowering G. densiflora. No divergent selection on floral morphology was detected, and there was no significant stabilizing selection on any trait in the two species. The results suggest ongoing adaptive differentiation of flowering phenology, strengthening this premating reproductive barrier between the two species. Synthesis: This study is among the first to test whether divergent selection on floral traits contribute to the maintenance of species differences between closely related plants. Phenological isolation confers a substantial potential for reproductive isolation, and divergent selection on flowering time can thus greatly influence reproductive isolation and adaptive differentiation.

Restriction-site associated DNA sequencing supports a sister group relationship of Nigritella and Gymnadenia (Orchidaceae)

The orchid genus Nigritella is closely related to Gymnadenia and has from time to time been merged with the latter. Although Nigritella is morphologically distinct, it has been suggested that the separating characters are easily modifiable and subject to rapid evolutionary change. So far, molecular phylogenetic studies have either given support for the inclusion of Nigritella in Gymnadenia, or for their separation as different genera. To resolve this issue, we analysed data obtained from Restriction-site associated DNA sequencing, RADseq, which provides a large number of SNPs distributed across the entire genome. To analyse samples of different ploidies, we take an analytical approach of building a reduced genomic reference based on de novo RADseq loci reconstructed from diploid accessions only, which we further use to map and call variants across both diploid and polyploid accessions. We found that Nigritella is distinct from Gymnadenia forming a well-supported separate clade, and that genetic diversity within Gymnadenia is high. Within Gymnadenia, taxa characterized by an ITS-E ribotype (G. conopsea s.str. (early flowering) and G. odoratissima), are divergent from taxa characterized by ITS-L ribotype (G. frivaldii, G. densiflora and late flowering G. conopsea). Gymnigritella runei is confirmed to have an allopolyploid origin from diploid Gymnadenia conopsea and tetraploid N. nigra ssp. nigra on the basis of RADseq data. Within Nigritella the aggregation of polyploid members into three clear-cut groups as suggested by allozyme and nuclear microsatellite data was further supported.

Habitat preference and flowering-time variation contribute to reproductive isolation between diploid and autotetraploid Anacamptis pyramidalis

Tetraploid lineages are typically reproductively isolated from their diploid ancestors by post-zygotic isolation via triploid sterility. Nevertheless, polyploids often also exhibit ecological divergence that could contribute to reproductive isolation from diploid ancestors. In this study, we disentangled the contribution of different forms of reproductive isolation between sympatric diploid and autotetraploid individuals of the food-deceptive orchid Anacamptis pyramidalis by quantifying the strength of seven reproductive barriers: three prepollination, one post-pollination prezygotic and three post-zygotic. The overall reproductive isolation between the two cytotypes was found very high, with a preponderant contribution of two prepollination barriers, that is phenological and microhabitat differences. Although the contribution of post-zygotic isolation (triploid sterility) is confirmed in our study, these results highlight that prepollination isolation, not necessarily involving pollinator preference, can represent a strong component of reproductive isolation between different cytotypes. Thus, in the context of polyploidy as quantum speciation, that generates reproductive isolation via triploid sterility, ecological divergence can strengthen the reproductive isolation between cytotypes, reducing the waste of gametes in low fitness interploidy crosses and thus favouring the initial establishment of the polyploid lineage. Under this light, speciation by polyploidy involves ecological processes and should not be strictly considered as a nonecological form of speciation.

Effects of population structure on pollen flow, clonality rates and reproductive success in fragmented Serapias lingua populations

BACKGROUND

Fragmentation of habitats by roads, railroads, fields, buildings and other human activities can affect population size, pollination success, sexual and asexual reproduction specially in plants showing pollinator limitation, such as Mediterranean orchids. In this study, we assessed pollen flow, selfing rates, vegetative reproduction and female reproductive success and their correlations with habitat characters in nine fragmented subpopulations of Serapias lingua. To improve understanding of population structure effects on plant biology, we examined genetic differentiation among populations, pollen flow, selfing rates and clonal reproduction using nuclear microsatellite markers.

RESULTS

Smaller populations showed a significant heterozygote deficit occurred at all five nuclear microsatellite loci, the coefficient of genetic differentiation among populations was 0.053 and pairwise FST was significantly correlated with the geographical distance between populations. Paternity analysis of seeds showed that most pollen flow occurred within a population and there was a positive correlation between percentage of received pollen and distance between populations. The fruit production rate varied between 5.10 % and 20.30 % and increased with increasing population size, while the percentage of viable seeds (78-85 %) did not differ significantly among populations. The extent of clonality together with the clonal and sexual reproductive strategies varied greatly among the nine populations and correlated with the habitats where they occur. The small, isolated populations tended to have high clonal diversity and low fruit production, whereas the large populations with little disturbance were prone to have reductions in clonal growth and increased sexual reproduction.

CONCLUSIONS

We found that clonality offers an advantage in small and isolated populations of S. lingua, where clones may have a greater ability to persist than sexually reproducing individuals.

Floral isolation is the major reproductive barrier between a pair of rewarding orchid sister species

The crucial role of reproductive isolation in speciation has long been recognized; however, a limited number of studies quantify different isolation barriers and embed reproductive isolation in a phylogenetic context. In this study, we investigate reproductive isolation between the often sympatrically occurring orchid species, Gymnadenia conopsea and G. odoratissima. We examine the phylogenetic relationship between the two species and analyse floral isolation, fruit set and seed viability from interspecies crosses, as well as the ploidy level. Additionally, we quantify interspecies differences in floral signals and morphology. The results suggest that the two species have a sister-species relationship. In terms of reproductive isolation, we found complete floral isolation between the two species, but little to no post-pollination isolation; the species also mostly had the same ploidy level in the studied populations. We also show clear distinctions in floral signals, as well as in floral size and spur length. We propose that respective adaptation to short- vs. long-tongued pollinators was the driver of speciation in the here studied Gymnadenia species. Our study supports the key role of floral isolation in orchid speciation and shows that floral isolation is not restricted to highly specialized pollination systems, but can also occur between species with less specialized pollination.

Minority cytotypes in European populations of the Gymnadenia conopsea complex (Orchidaceae) greatly increase intraspecific and intrapopulation diversity

BACKGROUND AND AIMS

Patterns of ploidy variation among and within populations can provide valuable insights into the evolutionary mechanisms shaping the dynamics of plant systems showing ploidy diversity. Whereas data on majority ploidies are, by definition, often sufficiently extensive, much less is known about the incidence and evolutionary role of minority cytotypes.

METHODS

Ploidy and proportions of endoreplicated genome were determined using DAPI (4',6-diamidino-2-phenylindole) flow cytometry in 6150 Gymnadenia plants (fragrant orchids) collected from 141 populations in 17 European countries. All widely recognized European species, and several taxa of less certain taxonomic status were sampled within Gymnadenia conopsea sensu lato.

KEY RESULTS

Most Gymnadenia populations were taxonomically and/or ploidy heterogeneous. Two majority (2x and 4x) and three minority (3x, 5x and 6x) cytotypes were identified. Evolution largely proceeded at the diploid level, whereas tetraploids were much more geographically and taxonomically restricted. Although minority ploidies constituted <2 % of the individuals sampled, they were found in 35 % of populations across the entire area investigated. The amount of nuclear DNA, together with the level of progressively partial endoreplication, separated all Gymnadenia species currently widely recognized in Europe.

CONCLUSIONS

Despite their low frequency, minority cytotypes substantially increase intraspecific and intrapopulation ploidy diversity estimates for fragrant orchids. The cytogenetic structure of Gymnadenia populations is remarkably dynamic and shaped by multiple evolutionary mechanisms, including both the ongoing production of unreduced gametes and heteroploid hybridization. Overall, it is likely that the level of ploidy heterogeneity experienced by most plant species/populations is currently underestimated; intensive sampling is necessary to obtain a holistic picture.

Remarkable coexistence of multiple cytotypes of the Gymnadenia conopsea aggregate (the fragrant orchid): evidence from flow cytometry

BACKGROUND AND AIMS

One of the prerequisites for polyploid research in natural systems is knowledge of the geographical distribution of cytotypes. Here inter- and intrapopulational ploidy diversity was examined in the Gymnadenia conopsea aggregate in central Europe and potential explanations and evolutionary consequences of the observed spatial patterns investigated.

METHODS

DAPI flow cytometry supplemented by confirmatory chromosome counts was used to determine ploidy in 3581 samples of the G. conopsea aggregate from 43 populations. The fine-scale spatial pattern of cytotype distribution (intra- and interploidy associations) was analysed with univariate and bivariate K-functions.

KEY RESULTS

Gymnadenia tissues undergo a progressively partial endoreplication, which accounts for about 60 % and 75 % of the total genome in G. conopsea and G. densiflora, respectively. Flow cytometric profiles are therefore species-specific and can be used as a marker for rapid and reliable species recognition. Two majority (4x, 8x) and three minority (6x, 10x, 12x) cytotypes were found, often in mixed-ploidy populations (harbouring up to all five different ploidy levels). The scarcity of the minority cytotypes (about 2·7 %) suggests the existence of strong pre- or postzygotic mating barriers. Spatial structure was observed in plots of populations with the highest cytotype variation, including clumping of individuals of the same ploidy and negative association between tetra- and octoploids.

CONCLUSIONS

The remarkable ploidy coexistence in the G. conopsea aggregate has reshaped our perception of intrapopulational ploidy diversity under natural conditions. This system offers unique opportunities for studying processes governing the formation and establishment of polyploids and assessing the evolutionary significance of the various pre- and postzygotic mating barriers that maintain this ploidy mixture.

Fungi from the roots of the common terrestrial orchid Gymnadenia conopsea

The fungal community associated with the terrestrial photosynthetic orchid Gymnadenia conopsea was characterized through PCR-amplification directly from root extracted DNA and cloning of the PCR products. Six populations in two geographically distinct regions in Germany were investigated. New ITS-primers amplifying a wide taxonomic range including Basidiomycetes and Ascomycetes revealed a high taxonomic and ecological diversity of fungal associates, including typical orchid mycorrhizas of the Tulasnellaceae and Ceratobasidiaceae as well as several ectomycorrhizal taxa of the Pezizales. The wide spectrum of potential mycorrhizal partners may contribute to this orchid's ability to colonize different habitat types with their characteristic microbial communities. The fungal community of G. conopsea showed a clear spatial structure. With 43 % shared taxa the species composition of the two regions showed only little overlap. Regardless of regions, populations were highly variable concerning taxon richness, varying between 5 and 14 taxa per population. The spatial structure and the continuous presence of mycorrhizal taxa on the one hand and the low specificity towards certain fungal taxa on the other hand suggest that the fungal community associated with G. conopsea is determined by multiple factors. In this context, germination as well as pronounced morphological and genetic differentiation within G. conopsea deserve attention as potential factors affecting the composition of the fungal community.

Frequency-dependent selection and the evolution of assortative mating

A long-standing goal in evolutionary biology is to identify the conditions that promote the evolution of reproductive isolation and speciation. The factors promoting sympatric speciation have been of particular interest, both because it is notoriously difficult to prove empirically and because theoretical models have generated conflicting results, depending on the assumptions made. Here, we analyze the conditions under which selection favors the evolution of assortative mating, thereby reducing gene flow between sympatric groups, using a general model of selection, which allows fitness to be frequency dependent. Our analytical results are based on a two-locus diploid model, with one locus altering the trait under selection and the other locus controlling the strength of assortment (a "one-allele" model). Examining both equilibrium and nonequilibrium scenarios, we demonstrate that whenever heterozygotes are less fit, on average, than homozygotes at the trait locus, indirect selection for assortative mating is generated. While costs of assortative mating hinder the evolution of reproductive isolation, they do not prevent it unless they are sufficiently great. Assortative mating that arises because individuals mate within groups (formed in time or space) is most conducive to the evolution of complete assortative mating from random mating. Assortative mating based on female preferences is more restrictive, because the resulting sexual selection can lead to loss of the trait polymorphism and cause the relative fitness of heterozygotes to rise above homozygotes, eliminating the force favoring assortment. When assortative mating is already prevalent, however, sexual selection can itself cause low heterozygous fitness, promoting the evolution of complete reproductive isolation (akin to "reinforcement") regardless of the form of natural selection.

Hybrids and fruit set in a mixed flowering-time population of Gymnadenia conopsea (Orchidaceae)

We have recently found that the morphologically determined subspecies Gymnadenia conopsea ssp conopsea in Sweden includes early and late flowering individuals. We were interested in the interactions between the flowering time groups; if there were gene flow between them and if so this was detrimental or advantageous. A spatially mixed population of early and late flowering individuals was studied using three microsatellite loci. We measured patterns in genetic differentiation and inferred occurrence of hybridisation and introgression. Variation in flowering time, fertility and relative and absolute fruit set was measured. The pattern of introgression between flowering-time groups differed between loci. In two of the three investigated loci, allele separation was distinct between early and late flowering plants and one genetically obvious hybrid was infertile. In the third locus, several alleles were shared between the two flowering time variants. The degree of introgression was associated to fruit set failure, which was higher in the late flowering plants and lower in early flowering plants. A small group of early flowering individuals with somewhat delayed flowering compared to the main group was genetically distinct and had lower relative and absolute fruit set. This group was not genetically intermediate, but rather constituting an independent group, with lower fruit set possibly caused by absence of pollinators. There seem to be a strong barrier against introgression into the late flowering group which is kept genetically distinct and less diverse. The early flowering group is diverse, includes two subgroups and seems to benefit from gene flow.

Local-scale genetic structure in the peatmoss Sphagnum fuscum

Sphagnum (peatmoss) dominates huge areas of the Northern Hemisphere and acts as a significant carbon sink on a global scale, yet little is known about the genetic structure of Sphagnum populations. We investigated genetic structure within a population of the common peatmoss Sphagnum fuscum, to assess local patterns of genetic diversity and the spatial extent of clones. One hundred seventeen shoots were sampled from five transects in Fuglmyra, central Norway, and sequenced for three anonymous DNA regions. Five neighbourhood patches were marked along each transect, and from each patch, five stems were sampled for molecular analyses. Seventeen haplotypes could be distinguished and two major groups of haplotypes differed by 12 mutational steps. The two major haplotype groups differed significantly in microhabitat association along the distance to groundwater table and the pH gradients, indicating microhabitat differentiation. The haplotypes within these groups were all genetically similar, differing by one or two mutations. The most common haplotype occurred in four transects separated by 250-m distance. Most of the molecular variation in the population was found among transects, and within patches. Large dominating clones within each transect resulted in low variation explained by the among-patch-within-transect component of spatial structure. Mutation appears to account for a larger proportion of the population variation than recombination. Within the population, vegetative growth and asexual reproduction from gametophyte fragments dominate as the main reproductive mode.

Bayesian inference of evolutionary history from chloroplast microsatellites in the cosmopolitan weed Capsella bursa-pastoris (Brassicaceae)

Besides showing an extraordinary degree of phenotypic variability, Capsella bursa-pastoris (Brassicaceae) is also one of the world's most common plant species and a serious weed in many countries. We have employed a coalescent-based Bayesian analysis of chloroplast microsatellite data to infer demographic and evolutionary parameters of this species. Two different demographic models applied to data from seven chloroplast microsatellite loci among 59 accessions show that the effective population size of C. bursa-pastoris is very small indicating a rapid expansion of the species, a result that is in accordance with fossil and historical data. Against this background, analysis of flowering time variation among accessions suggests that ecotypic differentiation in flowering time has occurred recently in the species' history. Finally, our results also indicate that mononucleotide repeat loci in the chloroplast genome can deteriorate in relatively short periods of evolutionary time.

SSR analysis of the Medicago truncatula SARDI core collection reveals substantial diversity and unusual genotype dispersal throughout the Mediterranean basin

The world's oldest and largest Medicago truncatula collection is housed at the South Australian Research and Development Institute (SARDI). We used six simple sequence repeat (SSR) loci to analyse the genetic diversity and relationships between randomly selected individuals from 192 accessions in the core collection. M. truncatula is composed of three subspecies (ssp.): ssp. truncatula, ssp. longeaculeata, and ssp. tricycla. Analysis at the level of six SSR loci supports the concept of ssp. tricycla, all the samples of which showed unique alleles at two loci. Contingency Chi-squared tests were significant between ssp. tricycla and ssp. truncatula at four loci, suggesting a barrier to gene flow between these subspecies. In accessions defined as ssp. longeaculeata, no unique allelic distribution or diagnostic sizes were observed, suggesting this apparent ssp. is a morphological variant of ssp. truncatula. The data also suggest M. truncatula that exhibits unusually wide genotype dispersal throughout its native Mediterranean region, possibly due to animal and trade-related movements. Our results showed the collection to be highly diverse, exhibiting an average of 25 SSR alleles per locus, with over 90% of individuals showing discrete genotypes. The rich diversity of the SARDI collection provides an invaluable resource for studying natural allelic variation of M. truncatula. To efficiently exploit the variation in the SARDI collection, we have defined a subset of accessions (n = 61) that maximises the diversity.

Population structure attributable to reproductive time: isolation by time and adaptation by time

Many populations are composed of a mixture of individuals that reproduce at different times, and these times are often heritable. Under these conditions, gene flow should be limited between early and late reproducers, even within populations having a unimodal temporal distribution of reproductive activity. This temporal restriction on gene flow might be called "isolation by time" (IBT) to acknowledge its analogy with isolation by distance (IBD). IBD and IBT are not exactly equivalent, however, owing to differences between dispersal in space and dispersal in time. We review empirical studies of natural populations that provide evidence for IBT based on heritabilities of reproductive time and on molecular genetic differences associated with reproductive time. When IBT is present, variation in selection through the reproductive season may lead to adaptive temporal variation in phenotypic traits [adaptation by time (ABT)]. We introduce a novel theoretical model that shows how ABT increases as (i) selection on the trait increases; (ii) environmental influences on reproductive time decrease; (iii) the heritability of reproductive time increases; and (iv) the temporal distribution of reproductive activity becomes increasingly uniform. We then review empirical studies of natural populations that provide evidence for ABT by documenting adaptive temporal clines in phenotypic traits. The best evidence for IBT and ABT currently comes from salmonid fishes and flowering plants, but we expect that future work will show these processes are more widespread.