The effects of inheritance in tetraploids on genetic diversity and population divergence

GenAPoPop 1.0: A user-friendly software to analyse genetic diversity and structure from partially clonal and selfed autopolyploid organisms

Autopolyploidy is quite common in most clades of eukaryotes. The emergence of sequence-based genotyping methods with individual and marker tags now enables confident allele dosage, overcoming the main obstacle to the democratization of the population genetic approaches when studying ecology and evolution of autopolyploid populations and species. Reproductive modes, including clonality, selfing and allogamy, have deep consequences on the ecology and evolution of population and species. Analysing genetic diversity and its dynamics over generations is one efficient way to infer the relative importance of clonality, selfing and allogamy in populations. GenAPoPop is a user-friendly solution to compute the specific corpus of population genetic indices, including indices about genotypic diversity, needed to analyse partially clonal, selfed and allogamous polysomic populations genotyped with confident allele dosage. It also easily provides the posterior probabilities of quantitative reproductive modes in autopolyploid populations genotyped at two-time steps and a graphical representation of the minimum spanning trees of the genetic distances between polyploid individuals, facilitating the interpretation of the genetic coancestry between individuals in hierarchically structured populations. GenAPoPop complements the previously existing solutions, including SPAGEDI and POLYGENE, to use genotypings to study the ecology and evolution of autopolyploid populations. It was specially developed with a simple graphical interface and workflow, and comes with a simulator to facilitate practical courses and teaching of population genetics for autopolyploid populations.

Demographic history inference and the polyploid continuum

Polyploidy is an important generator of evolutionary novelty across diverse groups in the Tree of Life, including many crops. However, the impact of whole-genome duplication depends on the mode of formation: doubling within a single lineage (autopolyploidy) versus doubling after hybridization between two different lineages (allopolyploidy). Researchers have historically treated these two scenarios as completely separate cases based on patterns of chromosome pairing, but these cases represent ideals on a continuum of chromosomal interactions among duplicated genomes. Understanding the history of polyploid species thus demands quantitative inferences of demographic history and rates of exchange between subgenomes. To meet this need, we developed diffusion models for genetic variation in polyploids with subgenomes that cannot be bioinformatically separated and with potentially variable inheritance patterns, implementing them in the dadi software. We validated our models using forward SLiM simulations and found that our inference approach is able to accurately infer evolutionary parameters (timing, bottleneck size) involved with the formation of auto- and allotetraploids, as well as exchange rates in segmental allotetraploids. We then applied our models to empirical data for allotetraploid shepherd's purse (Capsella bursa-pastoris), finding evidence for allelic exchange between the subgenomes. Taken together, our model provides a foundation for demographic modeling in polyploids using diffusion equations, which will help increase our understanding of the impact of demography and selection in polyploid lineages.

Alternative Evolutionary Pathways in Paspalum Involving Allotetraploidy, Sexuality, and Varied Mating Systems

The genetic systems of Paspalum species have not been extensively studied. We analyzed the ploidy, reproductive mode, mating system, and fertility of four Paspalum species-Paspalum durifolium, Paspalum ionanthum, Paspalum regnellii, and Paspalum urvillei. An analysis of 378 individuals from 20 populations of northeastern Argentina was conducted. All populations of the four Paspalum species were pure tetraploid and had a sexual and stable reproductive mode. However, some populations of P. durifolium and P. ionanthum showed low levels of apospory. Populations of P. durifolium and P. ionanthum had low seed sets under self-pollination but were fertile under open pollination, showing that self-incompatibility likely caused self-sterility. In contrast, populations of P. regnellii or P. urvillei showed no evidence of apospory, and seed sets in both self- and open pollination conditions were high, suggesting that they are self-compatible due to the absence of pollen-pistil molecular incompatibility mechanisms. The evolutionary origin of the four Paspalum species could explain these differences. This study supplies valuable insights into the genetic systems of Paspalum species, which could have implications for their conservation and management.

Population structure and hybridization under contemporary and future climates in a heteroploid foundational shrub species (Artemisia tridentata)

Current and past climatic changes can shift plant climatic niches, which may cause spatial overlap or separation between related taxa. The former often leads to hybridization and introgression, which may generate novel variation and influence the adaptive capacity of plants. An additional mechanism facilitating adaptations to novel environments and an important evolutionary driver in plants is polyploidy as the result of whole genome duplication. Artemisia tridentata (big sagebrush) is a landscape-dominating foundational shrub in the western United States which occupies distinct ecological niches, exhibiting diploid and tetraploid cytotypes. Tetraploids have a large impact on the species' landscape dominance as they occupy a preponderance of the arid spectrum of A. tridentata range. Three distinct subspecies are recognized, which co-occur in ecotones - the transition zone between two or more distinct ecological niches - allowing for hybridization and introgression. Here we assess the genomic distinctiveness and extent of hybridization among subspecies at different ploidies under both contemporary and predicted future climates. We sampled five transects throughout the western United States where a subspecies overlap was predicted using subspecies-specific climate niche models. Along each transect, we sampled multiple plots representing the parental and the potential hybrid habitats. We performed reduced representation sequencing and processed the data using a ploidy-informed genotyping approach. Population genomic analyses revealed distinct diploid subspecies and at least two distinct tetraploid gene pools, indicating independent origins of the tetraploid populations. We detected low levels of hybridization (2.5%) between the diploid subspecies, while we found evidence for increased admixture between ploidy levels (18%), indicating hybridization has an important role in the formation of tetraploids. Our analyses highlight the importance of subspecies co-occurrence within these ecotones to maintain gene exchange and potential formation of tetraploid populations. Genomic confirmations of subspecies in the ecotones support the subspecies overlap predicted by the contemporary climate niche models. However, future mid-century projections of subspecies niches predict a substantial loss in range and subspecies overlap. Thus, reductions in hybridization potential could affect new recruitment of genetically variable tetraploids that are vital to this species' ecological role. Our results underscore the importance of ecotone conservation and restoration.

Learning to tango with four (or more): the molecular basis of adaptation to polyploid meiosis

Polyploidy, which arises from genome duplication, has occurred throughout the history of eukaryotes, though it is especially common in plants. The resulting increased size, heterozygosity, and complexity of the genome can be an evolutionary opportunity, facilitating diversification, adaptation and the evolution of functional novelty. On the other hand, when they first arise, polyploids face a number of challenges, one of the biggest being the meiotic pairing, recombination and segregation of the suddenly more than two copies of each chromosome, which can limit their fertility. Both for developing polyploidy as a crop improvement tool (which holds great promise due to the high and lasting multi-stress resilience of polyploids), as well as for our basic understanding of meiosis and plant evolution, we need to know both the specific nature of the challenges polyploids face, as well as how they can be overcome in evolution. In recent years there has been a dramatic uptick in our understanding of the molecular basis of polyploid adaptations to meiotic challenges, and that is the focus of this review.

Genetic differentiation and genetic structure of mixed-ploidy Camellia hainanica populations

Camellia hainanica, which is common in China's Hainan Province, is an important woody olive tree species. Due to many years of geographic isolation, C. hainanica has not received the attention it deserves, which limits the exploitation of germplasm resources. Therefore, it is necessary to study population genetic characteristics for further utilization and conservation of C. hainanica. In this study, 96 individuals in six wild Camellia hainanica populations were used for ploidy analysis of the chromosome number, and the genetic diversity and population structure were investigated using 12 pairs of SSR primers. The results show complex ploidy differentiation in C. hainanica species. The ploidy of wild C. hainanica includes tetraploid, pentaploid, hexaploid, heptaploid, octoploid and decaploid species. Genetic analysis shows that genetic diversity and genetic differentiation among populations are low. Populations can be divided into two clusters based on their genetic structure, which matches their geographic location. Finally, to further maintain the genetic diversity of C. hainanica, ex-situ cultivation and in-situ management measures should be considered to protect it in the future.

Inferring the Demographic History and Inheritance Mode of Tetraploid Species Using ABC

Genomic patterns of diversity and divergence are impacted by certain life history traits, reproductive systems, and demographic history. The latter is characterized by fluctuations in population sizes over time, as well as by temporal patterns of introgression. For a given organism, identifying a demographic history that deviates from the standard neutral model allows a better understanding of its evolution but also helps to reduce the risk of false positives when screening for molecular targets of natural selection. Tetraploid organisms and beyond have demographic histories that are complicated by the mode of polyploidization, the mode of inheritance, and different scenarios of gene flow between sub-genomes and diploid parental species. Here we provide guidelines for experimenters wishing to address these issues through a flexible statistical framework: approximate Bayesian computation (ABC). The emphasis is on the general philosophy of the approach to encourage future users to exploit the enormous flexibility of ABC beyond the limitations imposed by generalist data analysis pipelines.

Inference of Polyploid Origin and Inheritance Mode from Population Genomic Data

Whole-genome duplications yield varied chromosomal pairing patterns, ranging from strictly bivalent to multivalent, resulting in disomic and polysomic inheritance modes. In the bivalent case, homeologous chromosomes form pairs, where in a multivalent pattern all copies are homologous and are therefore free to pair and recombine. As sufficient sequencing data is more readily available than high-quality cytological assessments of meiotic behavior or population genetic assessment of allelic segregation, especially for non-model organisms, bioinformatics approaches to infer origins and inheritance modes of polyploids using short-read sequencing data are attractive. Here we describe two such approaches, where the first is based on distributions of allelic read depth at heterozygous sites within an individual, as the expectations of such distributions are different for disomic and polysomic inheritance modes. The second approach is more laborious and based on a phylogenetic assessment of partially phased haplotypes of a polyploid in comparison to the closest diploid relatives. We discuss the sources of deviations from expected inheritance patterns, advantages and pitfalls of both methods, effects of mating types on the performance of the methods, and possible future developments.

Population Genomic Analysis of Diploid-Autopolyploid Species

This chapter outlines an empirical analysis of genome-wide single-nucleotide polymorphism (SNP) variation and its underlying drivers among multiple natural populations within a diploid-autopolyploid species. The aim is to reconstruct the genetic structure among natural populations of varying ploidy and infer footprints of selection in these populations, framed around specific questions that are typically encountered when analyzing a mixed-ploidy data set,e.g., addressing the relevance of natural whole-genome duplication for speciation and adaptation. We briefly review the options for the analysis of polyploid population genomic data involving variant calling, population structure, demographic history inference, and selection scanning approaches. Further, we provide suggestions for methods and associated software, possible caveats, and examples of their application to mixed-ploidy and autopolyploid data sets.

Multiple introductions, polyploidy and mixed reproductive strategies are linked to genetic diversity and structure in the most widespread invasive plant across Southern Ocean archipelagos

Biological invasions in remote areas that experience low human activity provide unique opportunities to elucidate processes responsible for invasion success. Here we study the most widespread invasive plant species across the isolated islands of the Southern Ocean, the annual bluegrass, Poa annua. To analyse geographical variation in genome size, genetic diversity and reproductive strategies, we sampled all major sub-Antarctic archipelagos in this region and generated microsatellite data for 470 individual plants representing 31 populations. We also estimated genome sizes for a subset of individuals using flow cytometry. Occasional events of island colonization are expected to result in high genetic structure among islands, overall low genetic diversity and increased self-fertilization, but we show that this is not the case for P. annua. Microsatellite data indicated low population genetic structure and lack of isolation by distance among the sub-Antarctic archipelagos we sampled, but high population structure within each archipelago. We identified high levels of genetic diversity, low clonality and low selfing rates in sub-Antarctic P. annua populations (contrary to rates typical of continental populations). In turn, estimates of selfing declined in populations as genetic diversity increased. Additionally, we found that most P. annua individuals are probably tetraploid and that only slight variation exists in genome size across the Southern Ocean. Our findings suggest multiple independent introductions of P. annua into the sub-Antarctic, which promoted the establishment of genetically diverse populations. Despite multiple introductions, the adoption of convergent reproductive strategies (outcrossing) happened independently in each major archipelago. The combination of polyploidy and a mixed reproductive strategy probably benefited P. annua in the Southern Ocean by increasing genetic diversity and its ability to cope with the novel environmental conditions.

A population-level statistic for assessing Mendelian behavior of genotyping-by-sequencing data from highly duplicated genomes

Background

Given the economic and environmental importance of allopolyploids and other species with highly duplicated genomes, there is a need for methods to distinguish paralogs, i.e. duplicate sequences within a genome, from Mendelian loci, i.e. single copy sequences that pair at meiosis. The ratio of observed to expected heterozygosity is an effective tool for filtering loci but requires genotyping to be performed first at a high computational cost, whereas counting the number of sequence tags detected per genotype is computationally quick but very ineffective in inbred or polyploid populations. Therefore, new methods are needed for filtering paralogs.

Results

We introduce a novel statistic, H_ind/H_E, that uses the probability that two reads sampled from a genotype will belong to different alleles, instead of observed heterozygosity. The expected value of H_ind/H_E is the same across all loci in a dataset, regardless of read depth or allele frequency. In contrast to methods based on observed heterozygosity, it can be estimated and used for filtering loci prior to genotype calling. In addition to filtering paralogs, it can be used to filter loci with null alleles or high overdispersion, and identify individuals with unexpected ploidy and hybrid status. We demonstrate that the statistic is useful at read depths as low as five to 10, well below the depth needed for accurate genotype calling in polyploid and outcrossing species.

Conclusions

Our methodology for estimating H_ind/H_E across loci and individuals, as well as determining reasonable thresholds for filtering loci, is implemented in polyRAD v1.6, available at https://github.com/lvclark/polyRAD. In large sequencing datasets, we anticipate that the ability to filter markers and identify problematic individuals prior to genotype calling will save researchers considerable computational time.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-022-04635-9.

Little hope for the polyploid endemic Pyrenean Larkspur ( Delphinium montanum ): Evidences from population genomics and Ecological Niche Modeling

Species endemic to restricted geographical ranges represent a particular conservation issue, be it for their heritage interest. In a context of global change, this is particularly the case for plants which belong to high‐mountain ecosystems and, because of their ecological requirements, are doomed to survive or disappear on their “sky islands”. The Pyrenean Larkspur (Delphinium montanum, Ranunculaceae) is endemic to the Eastern part of the Pyrenees (France and Spain). It is now only observable at a dozen of localities and some populations show signs of decline, such as a recurrent lack of flowering. Implementing population genomics approach (e.g., RAD‐seq like) is particularly useful to understand genomic patterns of diversity and differentiation in order to provide recommendations in term of conservation. However, it remains challenging for species such as D. montanum that are autotetraploid with a large genome size (1C‐value >10 pg) as most methods currently available were developed for diploid species. A Bayesian framework able to call genotypes with uncertainty allowed us to assess genetic diversity and population structure in this system. Our results show evidence for inbreeding (mean G_IS = 0.361) within all the populations and substantial population structure (mean G_ST = 0.403) at the metapopulation level. In addition to a lack of connectivity between populations, spatial projections of Ecological Niche Modeling (ENM) analyses under different climatic scenarios predict a dramatic decrease of suitable habitat for D. montanum in the future. Based on these results, we discuss the relevance and feasibility of different conservation measures.

High-throughput sequencing-based microsatellite genotyping for polyploids to resolve allele dosage uncertainty and improve analyses of genetic diversity, structure and differentiation: A case study of the hexaploid Camellia oleifera

Conventional microsatellite (simple sequence repeat, SSR) genotyping methods cannot accurately identify polyploid genotypes leading to allele dosage uncertainty, introducing biases in population genetic analysis. Here, a new SSR genotyping method was developed to directly infer accurate polyploid genotypes. The frequency distribution of SSR sequences was obtained based on deep-coverage high-throughput sequencing data. Corrections were performed accounting for the "stutter peak" and amplification efficiency of SSR sequences. Perl scripts and an online SSR genotyping tool "SSRSeq" were provided to process the sequencing data and output genotypes with corrected allele dosages. Hexaploid Camellia oleifera is the dominant woody oilseed crop in China. Understanding the geographical pattern of genetic variation in wild C. oleifera is essential for the conservation and utilization of genetic resources. Six wild C. oleifera populations were sampled across geographical ranges in subtropical evergreen broadleaf forests of China. Using 35 SSR markers, the high-throughput sequencing-based SSRSeq method was applied to obtain accurate hexaploid genotypes of wild C. oleifera. The results demonstrated that the new method could resolve allele dosage uncertainty and considerably improve genetic diversity, structure and differentiation analyses for polyploids. The genetic variation patterns of wild C. oleifera across geographical ranges agree with the "central-marginal hypothesis", stating that genetic diversity is high in the central population and declines from the central to the peripheral populations, and genetic differentiation increases from the centre to the periphery. This method and findings can facilitate the utilization of wild C. oleifera genetic resources for the breeding of cultivated C. oleifera.

Spatial and Ecological Drivers of Genetic Structure in Greek Populations of Alkanna tinctoria (Boraginaceae), a Polyploid Medicinal Herb

Background and Aims: Quantifying genetic variation is fundamental to understand a species' demographic trajectory and its ability to adapt to future changes. In comparison with diploids, however, genetic variation and factors fostering genetic divergence remain poorly studied in polyploids due to analytical challenges. Here, by employing a ploidy-aware framework, we investigated the genetic structure and its determinants in polyploid Alkanna tinctoria (Boraginaceae), an ancient medicinal herb that is the source of bioactive compounds known as alkannin and shikonin (A/S). From a practical perspective, such investigation can inform biodiversity management strategies. Methods: We collected 14 populations of A. tinctoria within its main distribution range in Greece and genotyped them using restriction site-associated DNA sequencing. In addition, we included two populations of A. sieberi. By using a ploidy-aware genotype calling based on likelihoods, we generated a dataset of 16,107 high-quality SNPs. Classical and model-based analysis was done to characterize the genetic structure within and between the sampled populations, complemented by genome size measurements and chromosomal counts. Finally, to reveal the drivers of genetic structure, we searched for associations between allele frequencies and spatial and climatic variables. Key Results: We found support for a marked regional structure in A. tinctoria along a latitudinal gradient in line with phytogeographic divisions. Several analyses identified interspecific admixture affecting both mainland and island populations. Modeling of spatial and climatic variables further demonstrated a larger contribution of neutral processes and a lesser albeit significant role of selection in shaping the observed genetic structure in A. tinctoria. Conclusion: Current findings provide evidence of strong genetic structure in A. tinctoria mainly driven by neutral processes. The revealed natural genomic variation in Greek Alkanna can be used to further predict variation in A/S production, whereas our bioinformatics approach should prove useful for the study of other non-model polyploid species.

Type and intensity of surrounding human land use, not local environment, shape genetic structure of a native grassland plant

Landscape heterogeneity can shape genetic structure and functional connectivity of populations. When this heterogeneity imposes variable costs of moving across the landscape, populations can be structured according to a pattern of "isolation by resistance" (IBR). At the same time, divergent local environmental filters can limit gene flow, creating an alternative pattern of "isolation by environment" (IBE). Here, we evaluate IBR and IBE in the insect-pollinated, biennial plant Sabatia angularis (L.) Pursh (Gentianaceae) across serpentine grasslands in the fragmented landscape of SE Pennsylvania, USA using ~4500 neutral SNP loci. Specifically, we test the extent to which radical alteration of the landscape matrix by humans has fundamentally altered the cost of movement, imprinting a pattern of IBR dictated by land use type and intensity, and the potential for IBE in relation to a gradient of heavy metal concentrations found in serpentine soil. We reveal a strong signal of IBR and a weak signal of IBE across sites, indicating the greater importance of the landscape matrix in shaping genetic structure of S. angularis populations in the study region. Based on Circuitscape and least cost path approaches, we find that both low- and high-intensity urbanization resist gene flow by orders of magnitude greater than "natural" habitats, although resistance to low-intensity urbanization weakens at larger spatial scales. While cropland presents a substantially lower barrier than urban development, cumulative human land use surrounding populations predicts within-population genetic diversity and inbreeding in S. angularis. Our results emphasize the role of forest buffers and corridors in facilitating gene flow between serpentine grassland patches and averting local extinction of plant populations.

Skipping sex: A nonrecombinant genomic assemblage of complementary reproductive modules

The unusual occurrence and developmental diversity of asexual eukaryotes remain a puzzle. De novo formation of a functioning asexual genome requires a unique assembly of sets of genes or gene states to disrupt cellular mechanisms of meiosis and gametogenesis, and to affect discrete components of sexuality and produce clonal or hemiclonal offspring. We highlight two usually overlooked but essential conditions to understand the molecular nature of clonal organisms, that is, a nonrecombinant genomic assemblage retaining modifiers of the sexual program, and a complementation between altered reproductive components. These subtle conditions are the basis for physiologically viable and genetically balanced transitions between generations. Genomic and developmental evidence from asexual animals and plants indicates the lack of complementation of molecular changes in the sexual reproductive program is likely the main cause of asexuals' rarity, and can provide an explanatory frame for the developmental diversity and lability of developmental patterns in some asexuals as well as for the discordant time to extinction estimations.

Spatial, climate and ploidy factors drive genomic diversity and resilience in the widespread grass Themeda triandra

Global climate change poses a significant threat to natural communities around the world, with many plant species showing signs of climate stress. Grassland ecosystems are not an exception, with climate change compounding contemporary pressures such as habitat loss and fragmentation. In this study, we assess the climate resilience of Themeda triandra, a foundational species and the most widespread plant in Australia, by assessing the relative contributions of spatial, environmental and ploidy factors to contemporary genomic variation. Reduced-representation genome sequencing on 472 samples from 52 locations was used to test how the distribution of genomic variation, including ploidy polymorphism, supports adaptation to hotter and drier climates. We explicitly quantified isolation by distance (IBD) and isolation by environment (IBE) and predicted genomic vulnerability of populations to future climates based on expected deviation from current genomic composition. We found that a majority (54%) of genomic variation could be attributed to IBD, while an additional 22% (27% when including ploidy information) could be explained by two temperature and two precipitation climate variables demonstrating IBE. Ploidy polymorphisms were common within populations (31/52 populations), indicating that ploidy mixing is characteristic of T. triandra populations. Genomic vulnerabilities were found to be heterogeneously distributed throughout the landscape, and our analysis suggested that ploidy polymorphism, along with other factors linked to polyploidy, reduced vulnerability to future climates by 60% (0.25-0.10). Our data suggests that polyploidy may facilitate adaptation to hotter climates and highlight the importance of incorporating ploidy in adaptive management strategies to promote the resilience of this and other foundation species.

Conservation of the Threatened Species, Pulsatilla vulgaris Mill. (Pasqueflower), is Aided by Reproductive System and Polyploidy

Population loss due to habitat disturbance is a major concern in biodiversity conservation. Here we investigate the genetic causes of the demographic decline observed in English populations of Pulsatilla vulgaris and the consequences for conservation. Using 10 nuclear microsatellite markers, we compare genetic variation in wild populations with restored and seed-regenerated populations (674 samples). Emergence of genetic structure and loss of allelic variation in natural populations are not as evident as expected from demographic trends. Restored populations show genetic variation comparable to their source populations and, in general, to the wild ones. Genetic homogeneity is observed in regeneration trials, although some alleles not captured in source populations are detected. We infer that polyploidy, longevity, and clonal reproduction have provided P. vulgaris with the standing genetic variation necessary to make the species resilient to the effects of demographic decline, suggesting that the use of multiple sources for reintroduction may be beneficial to mimic natural gene flow and the availability of multiple allele copies typical of polyploid species.

Hidden Hybridization and Habitat Differentiation in a Mediterranean Macrophyte, the Euryhaline Genus Ruppia

In many aquatic plant taxa, classification based on morphology has always been difficult. Molecular markers revealed that the complexity in several of these aquatic taxa could be addressed to recurrent hybridization events and cryptic species diversity. The submerged macrophyte genus Ruppia is one of these aquatic genera with a complex taxonomy due to the absence of clear distinguishable traits and several hybridization events. Two species co-exist throughout Europe, R. maritima and R. spiralis (previously known as R. cirrhosa), but recent molecular studies also found several indications of hybridization, introgression and chloroplast capture between these species. However, the full extent and frequency of hybridization and introgression in this genus has not been studied so far, nor is it clear how these hybrid lineages can co-exist locally with their parental species. In this paper, we wanted to detect whether a single coastal wetland where both species co-exist can act as a Ruppia hybrid zone. As a case study, we chose the Camargue, a Mediterranean coastal wetland that harbors a wide diversity in aquatic habitats, especially in terms of salinity and hydro-regime. We sampled several Ruppia populations within this wetland. To identify each sample and reconstruct the local genetic structure of the two parental species and their hybrids, we used both chloroplast and nuclear microsatellite markers. Afterward, we tested whether different species had different habitat preferences. Our results confirmed that R. maritima and R. spiralis are two strongly divergent species with different reproductive ecologies and different habitat preferences. This prevents frequent hybridization and consequently we could not detect any trace of a recent hybridization event. However, we found several populations of later-generation hybrids, including a population of R. maritima x hybrid backcrosses. The hybrid populations occupy a different habitat and are genetically distinct from their parental species, although they tend to be morphological similar to parental R. maritima. Although local hybridization and introgression in Ruppia is less frequent than we expected, the taxonomy of Ruppia is complicated due to ancient hybridizations and several back-crossings.

Population genetic structures of two ecologically distinct species Betula platyphylla and B. ermanii inferred based on nuclear and chloroplast DNA markers

Climatic oscillations during the last glacial maximum (LGM) significantly affected the distribution patterns and genetic structure of extant plants. Northeast China (NEC) is a major biodiversity center in East Asia, and the influence of historical climate change on NEC populations is critical for understanding species responses to future climate change. However, only a few phylogeographic studies of cool temperate deciduous tree species have been conducted in the area, and results are inconsistent for species with different niches or distribution areas. We employed multiple chloroplast and nuclear markers to investigate the genetic structure of two ecologically contrasting species, Betula platyphylla and B. ermanii, in NEC. Rare haplotypes were identified in the chloroplast genome of these species, and both exhibited high levels of nucleotide diversity based on a fragment of the nuclear gene G3PDH and microsatellites. Moreover, significant phylogeographic structure was detected for B. platyphylla, suggesting that these populations had recolonized from independent glacial refuges, whereas no genetic structure was found for B. ermanii.

OPEN RESEARCH BADGES

The nSSR datasets used in the current study and the table of pairwise FST (below diagonal) and its standardized F'ST (above diagonal) among 25 populations based on seven SSRs are available from the Dryad (DOI: https://doi.org/10.5061/dryad.230d176). Sequences generated from this study were deposited in GenBank under Accession nos. KY199568-KY200162 and MK819541-MK819970.

Role of ploidy in colonization of alpine habitats in natural populations of Arabidopsis arenosa

BACKGROUND AND AIMS

Polyploidy is an important driver of plant diversification and adaptation to novel environments. As a consequence of genome doubling, polyploids often exhibit greater colonizing ability or occupy a wider ecological niche than diploids. Although elevation has been traditionally considered as a key driver structuring ploidy variation, we do not know if environmental and phenotypic differentiation among ploidy cytotypes varies along an elevational gradient. Here, we tested for the consequences of genome duplication on genetic diversity, phenotypic variation and habitat preferences on closely related diploid and tetraploid populations that coexist along approx. 2300 m of varying elevation.

METHODS

We sampled and phenotyped 45 natural diploid and tetraploid populations of Arabidopsis arenosa in one mountain range in Central Europe (Western Carpathians) and recorded abiotic and biotic variables at each collection site. We inferred genetic variation, population structure and demographic history in a sub-set of 29 populations genotyped for approx. 36 000 single nucleotide polymorphisms.

KEY RESULTS

We found minor effects of polyploidy on colonization of alpine stands and low genetic differentiation between the two cytotypes, mirroring recent divergence of the polyploids from the local diploid lineage and repeated reticulation events among the cytotypes. This pattern was corroborated by the absence of ecological niche differentiation between the two cytotypes and overall phenotypic similarity at a given elevation.

CONCLUSIONS

The case of A. arenosa contrasts with previous studies that frequently showed clear niche differentiation between cytotypes. Our work stresses the importance of considering genetic structure and past demographic processes when interpreting the patterns of ploidy distributions, especially in species that underwent recent polyploidization events.

Genetic diversity and phenotypic variation within hatchery-produced oyster cohorts predict size and success in the field

The rapid growth of the aquaculture industry to meet global seafood demand offers both risks and opportunities for resource management and conservation. In particular, hatcheries hold promise for stock enhancement and restoration, yet cultivation practices may lead to enhanced variation between populations at the expense of variation within populations, with uncertain implications for performance and resilience. To date, few studies have assessed how production techniques impact genetic diversity and population structure, as well as resultant trait variation in and performance of cultivated offspring. We collaborated with a commercial hatchery to produce multiple cohorts of the eastern oyster (Crassostrea virginica) from field-collected broodstock using standard practices. We recorded key characteristics of the broodstock (male : female ratio, effective population size), quantified the genetic diversity of the resulting cohorts, and tested their trait variation and performance across multiple field sites and experimental conditions. Oyster cohorts produced under the same conditions in a single hatchery varied almost twofold in genetic diversity. In addition, cohort genetic diversity was a significant positive predictor of oyster performance traits, including initial size and survival in the field. Oyster cohorts produced in the hatchery had lower within-cohort genetic variation and higher among-cohort genetic structure than adults surveyed from the same source sites. These findings are consistent with "sweepstakes reproduction" in oysters, even when manually spawned. A readily measured characteristic of broodstock, the ratio of males to females, was positively correlated with within-cohort genetic diversity of the resulting offspring. Thus, this metric may offer a tractable way both to meet short-term production goals for seafood demand and to ensure the capacity of hatchery-produced stock to achieve conservation objectives, such as the recovery of self-sustaining wild populations.

Fernández de Castro A, Seguí J, Traveset A, Sosa PA. Alpine species in dynamic insular ecosystems through time: conservation genetics and niche shift estimates of the endemic and vulnerable Viola cheiranthifolia

BACKGROUND AND AIMS

Alpine oceanic ecosystems are considered amongst the most ephemeral and restricted habitats, with a biota highly vulnerable to climate changes and disturbances. As an example of an alpine insular endemic, the past and future population genetic structure and diversity, and the future distribution of Viola cheiranthifolia (Violaceae), endemic to Tenerife (Canary Islands), were estimated. The main goals were to predict distribution changes of this alpine oceanic plant under climate change, and to assist in actions for its conservation.

METHODS

To perform population genetic analysis, 14 specific microsatellite markers and algorithms which considered the polyploid condition of V. cheiranthifolia were employed. The niche modelling approach incorporated temperature gradients, topography and snow cover maps. Models were projected into climate change scenarios to assess the extent of the altitudinal shifts of environmental suitability. Finally, simulations were performed to predict whether the environmental suitability loss will affect the genetic diversity of populations.

KEY RESULTS

Viola cheiranthifolia presents short dispersal capacity, moderate levels of genetic diversity and a clear population genetic structure divided into two main groups (Teide and Las Cañadas Wall), showing signs of recolonization dynamics after volcanic eruptions. Future estimates of the distribution of the study populations also showed that, despite being extremely vulnerable to climate change, the species will not lose all its potential area in the next decades. The simulations to estimate genetic diversity loss show that it is correlated to suitability loss, especially in Las Cañadas Wall.

CONCLUSIONS

The low dispersal capacity of V. cheiranthifolia, coupled with herbivory pressure, mainly from rabbits, will make its adaptation to future climate conditions in this fragile alpine ecosystem difficult. Conservation actions should be focused on herbivore control, population reinforcement and surveillance of niche shifts, especially in Guajara, which represents the oldest isolated population and a genetic reservoir for the species.

The Site Frequency/Dosage Spectrum of Autopolyploid Populations

The Site Frequency Spectrum (SFS) and the heterozygosity of allelic variants are among the most important summary statistics for population genetic analysis of diploid organisms. We discuss the generalization of these statistics to populations of autopolyploid organisms in terms of the joint Site Frequency/Dosage Spectrum and its expected value for autopolyploid populations that follow the standard neutral model. Based on these results, we present estimators of nucleotide variability from High-Throughput Sequencing (HTS) data of autopolyploids and discuss potential issues related to sequencing errors and variant calling. We use these estimators to generalize Tajima's D and other SFS-based neutrality tests to HTS data from autopolyploid organisms. Finally, we discuss how these approaches fail when the number of individuals is small. In fact, in autopolyploids there are many possible deviations from the Hardy–Weinberg equilibrium, each reflected in a different shape of the individual dosage distribution. The SFS from small samples is often dominated by the shape of these deviations of the dosage distribution from its Hardy–Weinberg expectations.

Insights into population genetics and evolution of polyploids and their ancestors

We have developed the first comprehensive simulator for polyploid genomes (PolySim) and demonstrated its value by performing large-scale simulations to examine the effect of different population parameters on the evolution of polyploids. PolySim is unlimited in terms of ploidy, population size or number of simulated loci. Our process considered the evolution of polyploids from diploid ancestors, polysomic inheritance, inbreeding, recombination rate change in polyploids and gene flow from lower to higher ploidies. We compared the number of segregating single nucleotide polymorphisms, minor allele frequency, heterozygosity, R2 and average kinship relatedness between different simulated scenarios, and to real data from polyploid species. As expected, allotetraploid populations showed no difference from their ancestral diploids when population size remained constant and there was no gene flow or multivalent (MV) pairing between subgenomes. Autotetraploid populations showed significant differences from their ancestors for most parameters and diverged from their ancestral populations faster than allotetraploids. Autotetraploids can have significantly higher heterozygosity, relatedness and extended linkage disequilibrium compared with allotetraploids. Interestingly, autotetraploids were more sensitive to increasing selfing rate and decreasing population size. MV formation can homogenize allotetraploid subgenomes, but this homogenization requires a higher MV rate than previously proposed. Our results can be considered as the first building block to understand polyploid population evolutionary dynamics. PolySim can be used to simulate a wide variety of polyploid organisms that mimic empirical populations, which, in combination with quantitative genetics tools, can be used to investigate the power of genomewide association, genomic selection or breeding programme designs in these species.

Competition for Chiasma Formation Between Identical and Homologous (But Not Identical) Chromosomes in Synthetic Autotetraploids of Arabidopsis thaliana

Polyploid organisms provide additional opportunities to study meiosis in a more complex context since more than two potential homologous chromosomes are available. When the chromosome complement of a diploid individual is duplicated, each chromosome is accompanied by one identical and two homologous chromosomes within the same nucleus. In this situation, a competition in pairing/synapsis/chiasma formation between identical and homologous (but not necessarily identical) chromosomes can occur. Several studies have been conducted in different species to address whether there are preferences in crossover formation between identical rather than homologous chromosomes. In this study, multivalent and chiasma frequencies were cytologically analyzed in synthetic autotetraploids of Arabidopsis thaliana including the accessions Col, Ler, and the Col/Ler hybrid. Fluorescence in situ hybridization was conducted to identify each chromosome at metaphase I. The new Col and Ler tetraploids showed high multivalent frequencies, exceeding the theoretical 66.66% expected on a simple random end-pairing model, thus indicating that there are more than two autonomous synaptic sites per chromosome despite their small size. However, a significant excess of bivalent pairs was found in the Col/Ler hybrid, mainly due to the contribution of chromosomes 2 and 3. The mean chiasma frequencies of the three artificial autotetraploids were about twofold the corresponding mean cell chiasma frequencies of their diploid counterparts. The relative contribution of each chromosome to the total chiasma frequency was similar in the three genotypes, with the exception of a lower contribution of chromosome 3 in the hybrid. Preferences for chiasma formation between identical and homologous chromosomes were analyzed in Col/Ler 4x, taking advantage of the cytological differences between the accessions: variations in the size of the 45S rDNA region on the short arm of chromosome 2 and changes in the size and localization of the 5S rDNA region in chromosome 3. We observed a different behavior of chromosomes 2 and 3, i.e., random chiasma formation between identical and homologous chromosomes 2, and preferences for chiasma formation between homologous chromosomes 3. Hence, our results reveal the existence of chromosome-specific mechanisms responsible for these preferences.

A two-locus model of selection in autotetraploids: Chromosomal gametic disequilibrium and selection for an adaptive epistatic gene combination

In this paper, we present a two-locus model of selection for an autotetraploid population. We also investigate a measure of disequilibrium that occurs between homologous chromosomes in the diploid gametes of autotetraploids, namely chromosomal gametic disequilibrium. We apply the model and measure of disequilibrium to compare how an adaptive epistatic gene combination is inherited and selected for in an autotetraploid versus diploid population. Autotetraploids are expected to have higher genomic mutation and recombination rates relative to diploids, due to a greater ploidy level. These two processes can work in opposition in terms of selection for adaptive epistatic gene combinations. While a higher genomic mutation rate can generate the alleles that confer an epistatic combination more quickly, a higher recombination rate is expected to break the combination down more quickly. We show that chromosomal gametic disequilibrium in autotetraploids can potentially compensate for less linkage disequilibrium in autotetraploids. We also explore how double reduction affects the inheritance of and selection for an epistatic gene combination. Over all, our analysis provides theoretical evidence that adaptive epistatic combinations can be selected for more efficiently in autotetraploids versus diploids. This may provide insight into empirical work that finds epistasis has a role in causing population differentiation between autotetraploid plant populations.

Population assignment in autopolyploids

Understanding the patterns of contemporary gene dispersal within and among populations is of critical importance to population genetics and in managing populations for conservation. In contrast to diploids, there are few studies of gene dispersal in autopolyploids, in part due to complex polysomic inheritance and genotype ambiguity. Here we develop a novel approach for population assignment for codominant markers for autotetraploids and autohexaploids. This method accounts for polysomic inheritance, unreduced gametes and unknown allele dosage. It can also utilise information regarding the origin and genotype of one parent for population assignment of maternal or paternal parents. Using simulations, we demonstrate that our approach achieves high levels of accuracy for assignment even when population divergence is low (F_ST~0.06) and with only 12 microsatellite loci. We also show that substantially higher accuracy is achieved when known maternal information is utilised, regardless of whether allele dosage is known. Although this novel method exhibited near identical levels of accuracy to Structure when population divergence was high, it performed substantially better for most parameters at moderate (F_ST=0.06) to low levels of divergence (F_ST=0.03). These methods fill an important gap in the toolset for autopolyploids and pave the way for investigating contemporary gene dispersal in a widespread group of organisms.

How petals change their spots: cis-regulatory re-wiring in Clarkia (Onagraceae)

A long-standing question in evolutionary developmental biology is how new traits evolve. Although most floral pigmentation studies have focused on how pigment intensity and composition diversify, few, if any, have explored how a pattern element can shift position. In the present study, we examine the genetic changes underlying shifts in the position of petal spots in Clarkia. Comparative transcriptome analyses were used to identify potential candidate genes responsible for spot formation. Co-segregation analyses in F₂ individuals segregating for different spot positions, quantitative PCR, and pyrosequencing, were used to confirm the role of the candidate gene in determining spot position. Transient expression assays were used to identify the expression domain of different alleles. An R2R3Myb transcription factor (CgMyb1) activated spot formation, and different alleles of CgMyb1 were expressed in different domains, leading to spot formation in different petal locations. Reporter assays revealed that promoters from different alleles determine different locations of expression. The evolutionary shift in spot position is due to one or more cis-regulatory changes in the promoter of CgMyb1, indicating that shifts in pattern element position can be caused by changes in a single gene, and that cis-regulatory rewiring can be used to alter the relative position of an existing character.

Allo-allo-triploid Sphagnum × falcatulum: single individuals contain most of the Holantarctic diversity for ancestrally indicative markers

BACKGROUND AND AIMS

Allopolyploids exhibit both different levels and different patterns of genetic variation than are typical of diploids. However, scant attention has been given to the partitioning of allelic information and diversity in allopolyploids, particularly that among homeologous monoploid components of the hologenome. Sphagnum × falcatulum is a double allopolyploid peat moss that spans a considerable portion of the Holantarctic. With monoploid genomes from three ancestral species, this organism exhibits a complex evolutionary history involving serial inter-subgeneric allopolyploidizations.

METHODS

Studying populations from three disjunct regions [South Island (New Zealand); Tierra de Fuego archipelago (Chile, Argentina); Tasmania (Australia)], allelic information for five highly stable microsatellite markers that differed among the three (ancestral) monoploid genomes was examined. Using Shannon information and diversity measures, the holoploid information, as well as the information within and among the three component monoploid genomes, was partitioned into separate components for individuals within and among populations and regions, and those information components were then converted into corresponding diversity measures.

KEY RESULTS

The majority (76 %) of alleles detected across these five markers are most likely to have been captured by hybridization, but the information within each of the three monoploid genomes varied, suggesting a history of recurrent allopolyploidization between ancestral species containing different levels of genetic diversity. Information within individuals, equivalent to the information among monoploid genomes (for this dataset), was relatively stable, and represented 83 % of the grand total information across the Holantarctic, with both inter-regional and inter-population diversification each accounting for about 5 % of the total information.

CONCLUSIONS

Sphagnum × falcatulum probably inherited the great majority of its genetic diversity at these markers by reticulation, rather than by subsequent evolutionary radiation. However, some post-hybridization genetic diversification has become fixed in at least one regional population. Methodology allowing statistical analysis of any ploidy level is presented.

Congruent population structure across paralogous and nonparalogous loci in Salish Sea chum salmon (Oncorhynchus keta)

Whole-genome duplications are major evolutionary events with a lasting impact on genome structure. Duplication events complicate genetic analyses as paralogous sequences are difficult to distinguish; consequently, paralogs are often excluded from studies. The effects of an ancient whole-genome duplication (approximately 88 MYA) are still evident in salmonids through the persistence of numerous paralogous gene sequences and partial tetrasomic inheritance. We use restriction site-associated DNA sequencing on 10 collections of chum salmon from the Salish Sea in the USA and Canada to investigate genetic diversity and population structure in both tetrasomic and rediploidized regions of the genome. We use a pedigree and high-density linkage map to identify paralogous loci and to investigate genetic variation across the genome. By applying multivariate statistical methods, we show that it is possible to characterize paralogous loci and that they display similar patterns of population structure as the diploidized portion of the genome. We find genetic associations with the adaptively important trait of run-timing in both sets of loci. By including paralogous loci in genome scans, we can observe evolutionary signals in genomic regions that have routinely been excluded from population genetic studies in other polyploid-derived species.

Genomic diversity guides conservation strategies among rare terrestrial orchid species when taxonomy remains uncertain

BACKGROUND AND AIMS

Species are often used as the unit for conservation, but may not be suitable for species complexes where taxa are difficult to distinguish. Under such circumstances, it may be more appropriate to consider species groups or populations as evolutionarily significant units (ESUs). A population genomic approach was employed to investigate the diversity within and among closely related species to create a more robust, lineage-specific conservation strategy for a nationally endangered terrestrial orchid and its relatives from south-eastern Australia.

METHODS

Four putative species were sampled from a total of 16 populations in the Victorian Volcanic Plain (VVP) bioregion and one population of a sub-alpine outgroup in south-eastern Australia. Morphological measurements were taken in situ along with leaf material for genotyping by sequencing (GBS) and microsatellite analyses.

KEY RESULTS

Species could not be differentiated using morphological measurements. Microsatellite and GBS markers confirmed the outgroup as distinct, but only GBS markers provided resolution of population genetic structure. The nationally endangered Diuris basaltica was indistinguishable from two related species ( D. chryseopsis and D. behrii ), while the state-protected D. gregaria showed genomic differentiation.

CONCLUSIONS

Genomic diversity identified among the four Diuris species suggests that conservation of this taxonomically complex group will be best served by considering them as one ESU rather than separately aligned with species as currently recognized. This approach will maximize evolutionary potential among all species during increased isolation and environmental change. The methods used here can be applied generally to conserve evolutionary processes for groups where taxonomic uncertainty hinders the use of species as conservation units.

Autopolyploidy genome duplication preserves other ancient genome duplications in Atlantic salmon (Salmo salar)

Salmonids (e.g. Atlantic salmon, Pacific salmon, and trouts) have a long legacy of genome duplication. In addition to three ancient genome duplications that all teleosts are thought to share, salmonids have had one additional genome duplication. We explored a methodology for untangling these duplications from each other to better understand them in Atlantic salmon. In this methodology, homeologous regions (paralogous/duplicated genomic regions originating from a whole genome duplication) from the most recent genome duplication were assumed to have duplicated genes at greater density and have greater sequence similarity. This assumption was used to differentiate duplicated gene pairs in Atlantic salmon that are either from the most recent genome duplication or from earlier duplications. From a comparison with multiple vertebrate species, it is clear that Atlantic salmon have retained more duplicated genes from ancient genome duplications than other vertebrates--often at higher density in the genome and containing fewer synonymous mutations. It may be that polysomic inheritance is the mechanism responsible for maintaining ancient gene duplicates in salmonids. Polysomic inheritance (when multiple chromosomes pair during meiosis) is thought to be relatively common in salmonids compared to other vertebrate species. These findings illuminate how genome duplications may not only increase the number of duplicated genes, but may also be involved in the maintenance of them from previous genome duplications as well.

Very high MHC Class IIB diversity without spatial differentiation in the mediterranean population of greater Flamingos

Background

Selective pressure from pathogens is thought to shape the allelic diversity of major histocompatibility complex (MHC) genes in vertebrates. In particular, both local adaptation to pathogens and gene flow are thought to explain a large part of the intraspecific variation observed in MHC allelic diversity. To date, however, evidence that adaptation to locally prevalent pathogens maintains MHC variation is limited to species with limited dispersal and, hence, reduced gene flow. On the one hand high gene flow can disrupt local adaptation in species with high dispersal rates, on the other hand such species are much more likely to experience spatial variation in pathogen pressure, suggesting that there may be intense pathogen mediated selection pressure operating across breeding sites in panmictic species. Such pathogen mediated selection pressure operating across breeding sites should therefore be sufficient to maintain high MHC diversity in high dispersing species in the absence of local adaptation mechanisms. We used the Greater Flamingo, Phoenicopterus roseus, a long-lived colonial bird showing a homogeneous genetic structure of neutral markers at the scale of the Mediterranean region, to test the prediction that higher MHC allelic diversity with no population structure should occur in large panmictic populations of long-distance dispersing birds than in other resident species.

Results

We assessed the level of allelic diversity at the MHC Class IIB exon 2 from 116 individuals born in four different breeding colonies of Greater Flamingo in the Mediterranean region. We found one of the highest allelic diversity (109 alleles, 2 loci) of any non-passerine avian species investigated so far relative to the number of individuals and loci genotyped. There was no evidence of population structure between the four major Mediterranean breeding colonies.

Conclusion

Our results suggest that local adaptation at MHC Class IIB in Greater Flamingos is constrained by high gene flow and high MHC diversity appears to be maintained by population wide pathogen-mediated selection rather than local pathogen-mediated selection. Further understanding of how pathogens vary across space and time will be crucial to further elucidate the mechanisms maintaining MHC diversity in species with large panmictic populations and high dispersal rates.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-0905-3) contains supplementary material, which is available to authorized users.

Don't throw the baby out with the bathwater: identifying and mapping paralogs in salmonids

Many eukaryotic genomes contain a large fraction of gene duplicates (or paralogs) as a result of ancient or recent whole-genome duplications (Ohno 1970; Jaillon et al. 2004; Kellis et al. 2004). Identifying paralogs with NGS data is a pervasive problem in both ancient polyploids and neopolyploids. Likewise, paralogs are often treated as a nuisance that has to be detected and removed (Everett et al. 2012). In this issue of Molecular Ecology Resources, Waples et al. (2015) show that exclusion might not be necessary and how we may miss out on important genomic information in doing so. They present a novel statistical approach to detect paralogs based on the segregation of RAD loci in haploid offspring and test their method by constructing linkage maps with and without these duplicated loci in chum salmon, Oncorhynchus keta (Fig.1). Their linkage map including the resolved paralogs shows that these are mostly located in the distal regions of several linkage groups. Particularly intriguing is their finding that these homoeologous regions appear impoverished in transposable elements (TE). Given the role that TE play in genome remodelling, it is noteworthy that these elements are of low abundance in regions showing residual tetrasomic inheritance. This raises the question whether re-diploidization is constrained in these regions and whether they might have a role to play in salmonid speciation. This study provides an original approach to identifying duplicated loci in species with a pedigree, as well as providing a dense linkage map for chum salmon, and interesting insights into the retention of gene duplicates in an ancient polyploid.

Inheritance in doubled-diploid clementine and comparative study with SDR unreduced gametes of diploid clementine

Tetraploid clementine displays mainly tetrasomic inheritance. Genetic structures of 2n SDR and 2 × gametes from DD clementine are complementary and will guides triploids citrus breeding strategies. Triploid breeding is developed worldwide to create new seedless cultivars. Citrus triploid hybrids can be recovered from 2x × 2x sexual hybridizations as a consequence of the formation of unreduced gametes (2n), or from 4x × 2x interploid hybridizations in which tetraploid parents used are most often doubled-diploid (DD). Here we have analyzed the inheritance in doubled-diploid clementine and compared the genetic structures of gametes of DD clementine with SDR unreduced gametes of diploid clementine. Parental heterozygosity restitution (PHR) with DD parents depends on the rate of preferential chromosome pairing and thus the proportion of disomic versus tetrasomic segregations. Doubled-diploid clementine largely exhibited tetrasomic segregation. However, three linkage groups had intermediate segregation and one had a tendency for disomy. Significant doubled reduction rates (DR) rates were observed in six of the nine LGs. Differences of PHR between 2n SDR and 2x DD gametes were highest in the centromeric region and progressively decreased toward the distal regions where they were not significant. Over all markers, PHR was lower (two-thirds) in SDR 2n gametes than in DD-derived diploid gametes. The two strategies appear complementary in terms of genotypic variability. Interploid 4x × 2x hybridization is potentially more efficient for developing new cultivars that are phenotypically closer to the diploid parent of the DD than sexual hybridization through SDR 2n gametes. Conversely, 2x × 2x triploidisation has the potential to produce novel products with characteristics for market segmentation strategies.

Genetic differentiation and admixture between sibling allopolyploids in the Dactylorhiza majalis complex

Allopolyploidization often happens recurrently, but the evolutionary significance of its iterative nature is not yet fully understood. Of particular interest are the gene flow dynamics and the mechanisms that allow young sibling polyploids to remain distinct while sharing the same ploidy, heritage and overlapping distribution areas. By using eight highly variable nuclear microsatellites, newly reported here, we investigate the patterns of divergence and gene flow between 386 polyploid and 42 diploid individuals, representing the sibling allopolyploids Dactylorhiza majalis s.s. and D. traunsteineri s.l. and their parents at localities across Europe. We make use in our inference of the distinct distribution ranges of the polyploids, including areas in which they are sympatric (that is, the Alps) or allopatric (for example, Pyrenees with D. majalis only and Britain with D. traunsteineri only). Our results show a phylogeographic signal, but no clear genetic differentiation between the allopolyploids, despite the visible phenotypic divergence between them. The results indicate that gene flow between sibling Dactylorhiza allopolyploids is frequent in sympatry, with potential implications for the genetic patterns across their entire distribution range. Limited interploidal introgression is also evidenced, in particular between D. incarnata and D. traunsteineri. Altogether the allopolyploid genomes appear to be porous for introgression from related diploids and polyploids. We conclude that the observed phenotypic divergence between D. majalis and D. traunsteineri is maintained by strong divergent selection on specific genomic areas with strong penetrance, but which are short enough to remain undetected by genotyping dispersed neutral markers.