Estimation of allele frequencies in polyploids under certain patterns of inheritance

Relative Genetic Homogeneity within a Phenotypically Diverse group: The Case of Lake Tana Labeobarbus (Cyprinidae) Species Flock, Ethiopia

The Lake Tana Labeobarbus species flock represents one of the world's most famous examples of lacustrine species radiations. Previous studies of this group have resulted in the description of at least 15 species based on their differences in functional morphology and definition of two clades (lacustrine and riverine spawning clades) based on life history traits. A total of 166 fish representing 14 Labeobarbus species were genotyped using 10 lineage-specific hexaploid microsatellite loci. Six of these loci were developed for this study based on DNA sequence contigs derived from a microsatellite-enriched genomic library of Labeobarbus intermedius from Lake Tana; the remaining four loci were obtained from a previous study. The genotypes of the 10 loci were analyzed to examine genetic diversity and population structure within Lake Tana Labeobarbus. Overall mean allelic richness (NA) was 17.6 alleles per locus and observed (Ho) and expected (He) heterozygosities were 0.84 ± 0.14 and 0.73 ± 0.09, respectively, across all Lake Tana Labeobarbus samples examined. Our analyses reveal that there is little genetic differentiation among species (FST = 0.020-0.099; only 10 of 91 species comparisons were significant), but moderate differentiation (FST = 0.11, p < 0.05) between lacustrine and riverine spawning populations. Relative to previous phylogenetic hypotheses, our phenetic analysis employing the R-based Analysis of Phylogenetics and Evolution (APE) program seems to perform marginally better in revealing lineages within Lake Tana Labeobarbus. Herein, our results are compared to a previous microsatellite-based study of the same populations.

Genetic and epigenetic variation separately contribute to range expansion and local metalliferous habitat adaptation during invasions of Chenopodium ambrosioides into China

BACKGROUND AND AIMS

Invasive plants often colonize wide-ranging geographical areas with various local microenvironments. The specific roles of epigenetic and genetic variation during such expansion are still unclear. Chenopodium ambrosioides is a well-known invasive alien species in China that can thrive in metalliferous habitats. This study aims to comprehensively understand the effects of genetic and epigenetic variation on the successful invasion of C. ambrosioides.

METHODS

We sampled 367 individuals from 21 heavy metal-contaminated and uncontaminated sites with a wide geographical distribution in regions of China. We obtained environmental factors of these sampling sites, including 13 meteorological factors and the contents of four heavy metals in soils. Microsatellite markers were used to investigate the demographic history of C. ambrosioides populations in China. We also analysed the effect of epigenetic variation on metalliferous microhabitat adaptation using methylation-sensitive amplified polymorphism (MSAP) markers. A common garden experiment was conducted to compare heritable phenotypic variations among populations.

KEY RESULTS

Two distinct genetic clusters that diverged thousands of years ago were identified, suggesting that the eastern and south-western C. ambrosioides populations in China may have originated from independent introduction events without recombination. Genetic variation was shown to be a dominant determinant of phenotypic differentiation relative to epigenetic variation, and further affected the geographical distribution pattern of invasive C. ambrosioides. The global DNA unmethylation level was reduced in metalliferous habitats. Dozens of methylated loci were significantly associated with the heavy metal accumulation trait of C. ambrosioides and may contribute to coping with metalliferous microenvironments.

CONCLUSIONS

Our study of C. ambrosioides highlighted the dominant roles of genetic variation in large geographical range expansion and epigenetic variation in local metalliferous habitat adaptation.

Assessment of the Genetic Distinctiveness and Uniformity of Pre-Basic Seed Stocks of Italian Ryegrass Varieties

Lolium multiflorum Lam., commonly known as Italian ryegrass, is a forage grass mostly valued for its high palatability and digestibility, along with its high productivity. However, Italian ryegrass has an outbreeding nature and therefore has high genetic heterogeneity within each variety. Consequently, the exclusive use of morphological descriptors in the existing varietal identification and registration process based on the Distinctness, Uniformity, and Stability (DUS) test results in an inadequately precise assessment. The primary objective of this work was to effectively test whether the uniformity observed at the phenological level within each population of Italian ryegrass was confirmed at the genetic level through an SSR marker analysis. In this research, using 12 polymorphic SSR loci, we analyzed 672 samples belonging to 14 different Italian ryegrass commercial varieties to determine the pairwise genetic similarity (GS), verified the distribution of genetic diversity within and among varieties, and investigated the population structure. Although the fourteen commercial varieties did not show elevated genetic differentiation, with only 13% of the total variation attributable to among-cultivar genetic variation, when analyzed as a core, each variety constitutes a genetic cluster on its own, resulting in distinct characteristics from the others, except for two varieties. In this way, by combining a genetic tool with the traditional morphological approach, we were able to limit biases linked to the environmental effect of field trials, assessing the real source of diversity among varieties and concretely answering the key requisites of the Plant Variety Protection (PVP) system.

Biogeographic Population Structure of Chimeric Blades of Porphyra in the Northeast Atlantic Reveals Southern Rich Gene Pools, Introgression and Cryptic Plasticity

The genus Porphyra sensu lato (Bangiaceae, Rhodophyta), an important seaweed grown in aquaculture, is the most genetically diverse group of the Class Bangiophyceae, but has poorly understood genetic variability linked to complex evolutionary processes. Genetic studies in the last decades have largely focused on resolving gene phylogenies; however, there is little information on historical population biogeography, structure and gene flow in the Bangiaceae, probably due to their cryptic nature, chimerism and polyploidy, which render analyses challenging. This study aims to understand biogeographic population structure in the two abundant Porphyra species in the Northeast Atlantic: Porphyra dioica (a dioecious annual) and Porphyra linearis (protandrous hermaphroditic winter annual), occupying distinct niches (seasonality and position on the shore). Here, we present a large-scale biogeographic genetic analysis across their distribution in the Northeast Atlantic, using 10 microsatellites and cpDNA as genetic markers and integrating chimerism and polyploidy, including simulations considering alleles derived from different ploidy levels and/or from different genotypes within the chimeric blade. For P. linearis, both markers revealed strong genetic differentiation of north-central eastern Atlantic populations (from Iceland to the Basque region of Northeast Iberia) vs. southern populations (Galicia in Northwest Iberia, and Portugal), with higher genetic diversity in the south vs. a northern homogenous low diversity. For. P. dioica, microsatellite analyses also revealed two genetic regions, but with weaker differentiation, and cpDNA revealed little structure with all the haplotypes mixed across its distribution. The southern cluster in P. linearis also included introgressed individuals with cpDNA from P. dioica and a winter form of P. dioica occurred spatially intermixed with P. linearis. This third entity had a similar morphology and seasonality as P. linearis but genomes (either nuclear or chloroplast) from P. dioica. We hypothesize a northward colonization from southern Europe (where the ancestral populations reside and host most of the gene pool of these species). In P. linearis recently established populations colonized the north resulting in homogeneous low diversity, whereas for P. dioica the signature of this colonization is not as obvious due to hypothetical higher gene flow among populations, possibly linked to its reproductive biology and annual life history.

Training Set Construction for Genomic Prediction in Auto-Tetraploids: An Example in Potato

Training set construction is an important prerequisite to Genomic Prediction (GP), and while this has been studied in diploids, polyploids have not received the same attention. Polyploidy is a common feature in many crop plants, like for example banana and blueberry, but also potato which is the third most important crop in the world in terms of food consumption, after rice and wheat. The aim of this study was to investigate the impact of different training set construction methods using a publicly available diversity panel of tetraploid potatoes. Four methods of training set construction were compared: simple random sampling, stratified random sampling, genetic distance sampling and sampling based on the coefficient of determination (CDmean). For stratified random sampling, population structure analyses were carried out in order to define sub-populations, but since sub-populations accounted for only 16.6% of genetic variation, there were negligible differences between stratified and simple random sampling. For genetic distance sampling, four genetic distance measures were compared and though they performed similarly, Euclidean distance was the most consistent. In the majority of cases the CDmean method was the best sampling method, and compared to simple random sampling gave improvements of 4-14% in cross-validation scenarios, and 2-8% in scenarios with an independent test set, while genetic distance sampling gave improvements of 5.5-10.5% and 0.4-4.5%. No interaction was found between sampling method and the statistical model for the traits analyzed.

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.

Genetic structure of coast redwood (Sequoia sempervirens [D. Don] Endl.) populations in and outside of the natural distribution range based on nuclear and chloroplast microsatellite markers

Coast redwood (Sequoia sempervirens) naturally growing in southern Oregon and northern California is one of the few conifer tree species that are polyploid. Despite its unique ecological and economic importance, its population genetic structure is still insufficiently studied. To obtain additional data on its population genetic structure we genotyped 317 samples collected from populations in California (data set C) and 144 trees growing in a provenance trial in France (data set F) using 12 nuclear (five random nuclear genomic nSSRs and seven expressed sequence tag EST-SSRs) and six chloroplast (cpSSRs) microsatellite or simple sequence repeat (SSR) markers, respectively. These data sets were also used as reference to infer the origin of 147 coast redwood trees growing in Germany (data set G). Coast redwood was introduced to Europe, including Germany as an ornamental species, decades ago. Due to its fast growth and high timber quality, it could be considered as a potential commercial timber species, especially in perspective to climate warming that makes more regions in Germany suitable for its growing. The well performing trees in colder Germany could be potential frost resistant genotypes, but their genetic properties and origin are mostly unknown. Within the natural range in southern Oregon and northern California, only two relatively weak clusters were identified, one northern and one southern, separated by the San Francisco Bay. High genetic diversity, but low differentiation was found based on the 12 nuclear SSR markers for all three data sets F, C and G. We found that investigated 147 German trees represented only 37 different genotypes. They showed genetic diversity at the level less than diversity observed within the natural range in the northern or southern cluster, but more similar to the diversity observed in the southern cluster. It was difficult to assign German trees to the original single native populations using the six cpSSR markers, but rather to either the northern or southern cluster. The high number of haplotypes found in the data sets based on six cpSSR markers and low genetic differentiation based on 12 nuclear SSRs found in this study helps us study and better understand population genetic structure of this complex polyploid tree and supports the selection of potential genotypes for German forestry.

Sequencing depth and genotype quality: accuracy and breeding operation considerations for genomic selection applications in autopolyploid crops

KEY MESSAGE

Polypoid crop breeders can balance resources between density and sequencing depth, dosage information and fewer highly informative SNPs recommended, non-additive models and QTL advantages on prediction dependent on trait architecture. The autopolyploid nature of potato and sweetpotato ensures a wide range of meiotic configurations and linkage phases leading to complex gene-action and pose problems in genotype data quality and genomic selection analyses. We used a 315-progeny biparental F1 population of hexaploid sweetpotato and a diversity panel of 380 tetraploid potato, genotyped using different platforms to answer the following questions: (i) do polyploid crop breeders need to invest more for additional sequencing depth? (ii) how many markers are required to make selection decisions? (iii) does considering non-additive genetic effects improve predictive ability (PA)? (iv) does considering dosage or quantitative trait loci (QTL) offer significant improvement to PA? Our results show that only a small number of highly informative single nucleotide polymorphisms (SNPs; ≤ 1000) are adequate for prediction in the type of populations we analyzed. We also show that considering dosage information and models considering only additive effects had the best PA for most traits, while the comparative advantage of considering non-additive genetic effects and including known QTL in the predictive model depended on trait architecture. We conclude that genomic selection can help accelerate the rate of genetic gains in potato and sweetpotato. However, application of genomic selection should be considered as part of optimizing the entire breeding program. Additionally, since the predictions in the current study are based on single populations, further studies on the effects of haplotype structure and inheritance on PA should be studied in actual multi-generation breeding populations.

Manipulation of cytosine methylation does not remove latitudinal clines in two invasive goldenrod species in Central Europe

Invasive species frequently differentiate phenotypically in novel environments within a few generations, often even with limited genetic variation. For the invasive plants Solidago canadensis and S. gigantea, we tested whether such differentiation might have occurred through heritable epigenetic changes in cytosine methylation. In a 2-year common-garden experiment, we grew plants from seeds collected along a latitudinal gradient in their non-native Central European range to test for trait differentiation and whether differentiation disappeared when seeds were treated with the demethylation agent zebularine. Microsatellite markers revealed no population structure along the latitudinal gradient in S. canadensis, but three genetic clusters in S. gigantea. Solidago canadensis showed latitudinal clines in flowering phenology and growth. In S. gigantea, the number of clonal offspring decreased with latitude. Although zebularine had a significant effect on early growth, probably through effects on cytosine methylation, latitudinal clines remained (or even got stronger) in plants raised from seeds treated with zebularine. Thus, our experiment provides no evidence that epigenetic mechanisms by selective cytosine methylation contribute to the observed phenotypic differentiation in invasive goldenrods in Central Europe.

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.

genodive version 3.0: Easy-to-use software for the analysis of genetic data of diploids and polyploids

genodive version 3.0 is a user-friendly program for the analysis of population genetic data. This version presents a major update from the previous version and now offers a wide spectrum of different types of analyses. genodive has an intuitive graphical user interface that allows direct manipulation of the data through transformation, imputation of missing data, and exclusion and inclusion of individuals, population and/or loci. Furthermore, genodive seamlessly supports 15 different file formats for importing or exporting data from or to other programs. One major feature of genodive is that it supports both diploid and polyploid data, up to octaploidy (2n = 8x) for some analyses, but up to hexadecaploidy (2n = 16x) for other analyses. The different types of analyses offered by genodive include multiple statistics for estimating population differentiation (φ_ST , F_ST , F'_ST , G_ST , G'_ST , G''_ST , D_est , R_ST , ρ), analysis of molecular variance-based K-means clustering, Hardy-Weinberg equilibrium, hybrid index, population assignment, clone assignment, Mantel test, Spatial Autocorrelation, 23 ways of calculating genetic distances, and both principal components and principal coordinates analyses. A unique feature of genodive is that it can also open data sets with nongenetic variables, for example environmental data or geographical coordinates that can be included in the analysis. In addition, genodive makes it possible to run several external programs (lfmm, structure, instruct and vegan) directly from its own user interface, avoiding the need for data reformatting and use of the command line. genodive is available for computers running Mac OS X 10.7 or higher and can be downloaded freely from: http://www.patrickmeirmans.com/software.

Genotypic Frequencies at Equilibrium for Polysomic Inheritance Under Double-Reduction

Polyploids are organisms whose genomes consist of more than two complete sets of chromosomes. Both autopolyploids and allopolyploids may display polysomic inheritance. A peculiarity of polysomic inheritance is multivalent formation during meiosis resulting in double-reduction, which occurs when sister chromatid fragments segregate into the same gamete. Double-reduction can result in gametes carrying identical-by-descent alleles and slightly increasing homozygosity. This will cause the genotypic frequencies to deviate from expected values and will thus bias the results of standard population genetic analytical methods used in molecular ecology and selective breeding. In this study, we extend existing double-reduction models to account for any even level of ploidy, and derive the symbolic expressions for genotypic frequencies via two methods. Inbreeding coefficients and heterozygosity under double-reduction and inbreeding are also calculated. Numerical solutions obtained by computer simulations are compared with analytical solutions predicted by the model to validate the model.

polyRAD: Genotype Calling with Uncertainty from Sequencing Data in Polyploids and Diploids

Low or uneven read depth is a common limitation of genotyping-by-sequencing (GBS) and restriction site-associated DNA sequencing (RAD-seq), resulting in high missing data rates, heterozygotes miscalled as homozygotes, and uncertainty of allele copy number in heterozygous polyploids. Bayesian genotype calling can mitigate these issues, but previously has only been implemented in software that requires a reference genome or uses priors that may be inappropriate for the population. Here we present several novel Bayesian algorithms that estimate genotype posterior probabilities, all of which are implemented in a new R package, polyRAD. Appropriate priors can be specified for mapping populations, populations in Hardy-Weinberg equilibrium, or structured populations, and in each case can be informed by genotypes at linked markers. The polyRAD software imports read depth from several existing pipelines, and outputs continuous or discrete numerical genotypes suitable for analyses such as genome-wide association and genomic prediction.

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.

Comparison of population genetic structures of the plant Silene stellata and its obligate pollinating seed predator moth Hadena ectypa

BACKGROUND AND AIMS

Population genetic structures and patterns of gene flow of interacting species provide important insights into the spatial scale of their interactions and the potential for local co-adaptation. We analysed the genetic structures of the plant Silene stellata and the nocturnal moth Hadena ectypa. Hadena ectypa acts as one of the important pollinators of S. stellata as well as being an obligate seed parasite on the plant. Although H. ectypa provides a substantial pollination service to S. stellata, this system is largely considered parasitic due to the severe seed predation by the Hadena larvae. Previous research on this system has found variable interaction outcomes across space, indicating the potential for a geographical selection mosaic.

METHODS

Using 11 microsatellite markers for S. stellata and nine markers for H. ectypa, we analysed the population genetic structure and the patterns and intensity of gene flow within and among three local populations in the Appalachians.

KEY RESULTS

We found no spatial genetic structure in the moth populations, while significant differentiation was detected among the local plant populations. Additionally, we observed that gene flow rates among H. ectypa populations were more uniform and that the mean gene flow rate in H. ectypa was twice as large as that in S. stellata.

CONCLUSIONS

Our results suggest that although the moths move frequently among populations, long-distance pollen carryover only happens occasionally. The difference in gene flow rates between S. stellata and H. ectypa could prevent strict local co-adaptation. Furthermore, higher gene flow rates in H. ectypa could also increase resistance of the local S. stellata populations to the parasitic effect of H. ectypa and therefore help to stabilize the Silene-Hadena interaction dynamics.

Different patterns of colonization of Oxalis alpina in the Sky Islands of the Sonoran desert via pollen and seed flow

Historical factors such as climatic oscillations during the Pleistocene epoch have dramatically impacted species distributions. Studies of the patterns of genetic structure in angiosperm species using molecular markers with different modes of inheritance contribute to a better understanding of potential differences in colonization and patterns of gene flow via pollen and seeds. These markers may also provide insights into the evolution of reproductive systems in plants. Oxalis alpina is a tetraploid, herbaceous species inhabiting the Sky Island region of the southwestern United States and northern Mexico. Our main objective in this study was to analyze the influence of climatic oscillations on the genetic structure of O. alpina and the impact of these oscillations on the evolutionary transition from tristylous to distylous reproductive systems. We used microsatellite markers and compared our results to a previous study using chloroplast genetic markers. The phylogeographic structure inferred by both markers was different, suggesting that intrinsic characteristics including the pollination system and seed dispersal have influenced patterns of gene flow. Microsatellites exhibited low genetic structure, showed no significant association between geographic and genetic distances, and all individual genotypes were assigned to two main groups. In contrast, chloroplast markers exhibited a strong association between geographic and genetic distance, had higher levels of genetic differentiation, and were assigned to five groups. Both types of DNA markers showed evidence of a northward expansion as a consequence of climate warming occurring in the last 10,000 years. The data from both types of markers support the hypothesis for several independent transitions from tristyly to distyly.

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.

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.

Ecological distributions, phenological isolation, and genetic structure in sympatric and parapatric populations of the Larrea tridentata polyploid complex

PREMISE OF THE STUDY

Polyploidy is widely recognized as a mechanism of diversification. Contributions of polyploidy to specific pre- and postzygotic barriers-and classifications of polyploid speciation as "ecological" vs. "non-ecological"-are more contentious. Evaluation of these issues requires comprehensive studies that test ecological characteristics of cytotypes as well as the coincidence of genetic structure with cytotype distributions.

METHODS

We investigated a classical example of autopolyploid speciation, Larrea tridentata, at multiple areas of cytotype co-occurrence. Habitat and phenological differences were compared between diploid, tetraploid, and hexaploid populations on the basis of edaphic, community composition, and flowering time surveys. Frequency of hybridization between diploids and tetraploids was investigated using a diploid-specific chloroplast DNA (cpDNA) marker; genetic structure for all cytotypes was assessed using amplified fragment length polymorphisms (AFLPs).

KEY RESULTS

Across contact zones, we found cytotypes in habitats distinguished by soil and vegetation. We observed modest differences in timing and production of flowers, indicating a degree of assortative mating that was asymmetric between cytotypes. Nonetheless, cpDNA analyses in diploid-tetraploid contact zones suggested that ∼5% of tetraploid plants had hybrid origins involving unilateral sexual polyploidization. Genetic structure of AFLPs largely coincided with cytotype distributions in diploid-tetraploid contact zones. In contrast, there was little structure in areas of contact between tetraploids and hexaploids, suggesting intercytotype gene flow or recurrent hexaploid formation.

CONCLUSIONS

Diploid, tetraploid, and hexaploid cytotypes of L. tridentata are segregated by environmental distributions and flowering phenology in contact zones, with diploid and tetraploid populations having corresponding differences in genetic structure.

Transcriptome-derived evidence supports recent polyploidization and a major phylogeographic division in Trithuria submersa (Hydatellaceae, Nymphaeales)

Relatively little is known about species-level genetic diversity in flowering plants outside the eudicots and monocots, and it is often unclear how to interpret genetic patterns in lineages with whole-genome duplications. We addressed these issues in a polyploid representative of Hydatellaceae, part of the water-lily order Nymphaeales. We examined a transcriptome of Trithuria submersa for evidence of recent whole-genome duplication, and applied transcriptome-derived microsatellite (expressed-sequence tag simple-sequence repeat (EST-SSR)) primers to survey genetic variation in populations across its range in mainland Australia. A transcriptome-based Ks plot revealed at least one recent polyploidization event, consistent with fixed heterozygous genotypes representing underlying sets of homeologous loci. A strong genetic division coincides with a trans-Nullarbor biogeographic boundary. Patterns of 'allelic' variation (no more than two variants per EST-SSR genotype) and recently published chromosomal evidence are consistent with the predicted polyploidization event and substantial homozygosity underlying fixed heterozygote SSR genotypes, which in turn reflect a selfing mating system. The Nullarbor Plain is a barrier to gene flow between two deep lineages of T. submersa that may represent cryptic species. The markers developed here should also be useful for further disentangling species relationships, and provide a first step towards future genomic studies in Trithuria.

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.

Accounting for genotype uncertainty in the estimation of allele frequencies in autopolyploids

Despite the increasing opportunity to collect large-scale data sets for population genomic analyses, the use of high-throughput sequencing to study populations of polyploids has seen little application. This is due in large part to problems associated with determining allele copy number in the genotypes of polyploid individuals (allelic dosage uncertainty-ADU), which complicates the calculation of important quantities such as allele frequencies. Here, we describe a statistical model to estimate biallelic SNP frequencies in a population of autopolyploids using high-throughput sequencing data in the form of read counts. We bridge the gap from data collection (using restriction enzyme based techniques [e.g. GBS, RADseq]) to allele frequency estimation in a unified inferential framework using a hierarchical Bayesian model to sum over genotype uncertainty. Simulated data sets were generated under various conditions for tetraploid, hexaploid and octoploid populations to evaluate the model's performance and to help guide the collection of empirical data. We also provide an implementation of our model in the R package polyfreqs and demonstrate its use with two example analyses that investigate (i) levels of expected and observed heterozygosity and (ii) model adequacy. Our simulations show that the number of individuals sampled from a population has a greater impact on estimation error than sequencing coverage. The example analyses also show that our model and software can be used to make inferences beyond the estimation of allele frequencies for autopolyploids by providing assessments of model adequacy and estimates of heterozygosity.

Genetic Diversity and Spatial Genetic Structure of the Grassland Perennial Saxifraga granulata along Two River Systems

Due to changes in land use, the natural habitats of an increasing number of plant species have become more and more fragmented. In landscapes that consist of patches of suitable habitat, the frequency and extent of long-distance seed dispersal can be expected to be an important factor determining local genetic diversity and regional population structure of the remaining populations. In plant species that are restricted to riparian habitats, rivers can be expected to have a strong impact on the dynamics and spatial genetic structure of populations as they may enable long-distance seed dispersal and thus maintain gene flow between fragmented populations. In this study, we used polymorphic microsatellite markers to investigate the genetic diversity and the spatial genetic structure of 28 populations of Saxifraga granulata along two rivers in central Belgium. We hypothesized that rivers might be essential for gene flow among increasingly isolated populations of this species. Genetic diversity was high (HS = 0.68), which to a certain extent can be explained by the octoploid nature of S. granulata in the study area. Populations along the Dijle and Demer rivers were also highly differentiated (G"ST = 0.269 and 0.164 and DEST = 0.190 and 0.124, respectively) and showed significant isolation-by-distance, indicating moderate levels of gene flow primarily between populations that are geographically close to each other. Along the river Demer population genetic diversity was higher upstream than downstream, suggesting that seed dispersal via the water was not the primary mode of dispersal. Overall, these results indicate that despite increasing fragmentation populations along both rivers were highly genetically diverse. The high ploidy level and longevity of S. granulata have most likely buffered negative effects of fragmentation on genetic diversity and the spatial genetic structure of populations in riparian grasslands.

Population genetics of autopolyploids under a mixed mating model and the estimation of selfing rate

Nowadays, the population genetics analysis of autopolyploid species faces many difficulties due to (i) limited development of population genetics tools under polysomic inheritance, (ii) difficulties to assess allelic dosage when genotyping individuals and (iii) a form of inbreeding resulting from the mechanism of 'double reduction'. Consequently, few data analysis computer programs are applicable to autopolyploids. To contribute bridging this gap, this article first derives theoretical expectations for the inbreeding and identity disequilibrium coefficients under polysomic inheritance in a mixed mating model. Moment estimators of these coefficients are proposed when exact genotypes or just markers phenotypes (i.e. allelic dosage unknown) are available. This led to the development of estimators of the selfing rate based on adult genotypes or phenotypes and applicable to any even-ploidy level. Their statistical performances and robustness were assessed by numerical simulations. Contrary to inbreeding-based estimators, the identity disequilibrium-based estimator using phenotypes is robust (absolute bias generally < 0.05), even in the presence of double reduction, null alleles or biparental inbreeding due to isolation by distance. A fairly good precision of the selfing rate estimates (root mean squared error < 0.1) is already achievable using a sample of 30-50 individuals phenotyped at 10 loci bearing 5-10 alleles each, conditions reachable using microsatellite markers. Diallelic markers (e.g. SNP) can also perform satisfactorily in diploids and tetraploids but more polymorphic markers are necessary for higher ploidy levels. The method is implemented in the software SPAGeDi and should contribute to reduce the lack of population genetics tools applicable to autopolyploids.

orchard: Paternity program for autotetraploid species

Advances in molecular marker technology have provided new opportunities to study the population genetics of polyploid taxa. Paternity analysis using microsatellite markers can be used in detection of gene flow between individuals and populations, in mating system analysis, to identify factors that influence fecundity and fertility, to identify behaviour of parent-offspring relationships and in the analysis of the reproductive success of different ecological groups. As there is no specific program for carrying out paternity analysis in tetraploid species, specialized software was designed for the assignment of paternity for autotetraploid species. orchard is a novel implementation of exclusion and likelihood statistics for carrying out paternity analysis of autotetraploids. First, the program performs an exclusion method, and then, a likelihood statistic is used with nonexcluded candidate fathers. Optional features include estimation of allele dosage of known mother trees and the estimation of pollen flow distances. orchard was tested using a data set of microsatellite data of Dipteryx odorata, a tetraploid Amazonian tree species.

Genetic analysis of resistance to Pseudomonas syringae pv. actinidiae (Psa) in a kiwifruit progeny test: an application of generalised linear mixed models (GLMMs)

Linear Mixed models (LMMs) that incorporate genetic and spatial covariance structures have been used for many years to estimate genetic parameters and to predict breeding values in animal and plant breeding. Although the theoretical aspects for extending LMM to generalised linear mixed models (GLMMs) have been around for some time, suitable software has been developed only within the last decade or so. The GLIMMIX procedure in SAS® is becoming popular for fitting GLMMs in various disciplines. Applications of GLMMs to genetic analysis have been limited, probably because of the complexity of the models used. This is particularly so for Proc GLIMMIX because, unlike ASReml software, it is not specifically tailored for analysis of breeding data and some pre-procedure coding is necessary. Binary data that fits the GLMM framework is commonly encountered in breeding experiments, such as when evaluating individuals for resistance by observing the presence or absence of disease. Bacterial canker (Psa) caused by Pseudomonas syringae pv. actinidiae is a serious disease of kiwifruit in New Zealand and other kiwifruit-producing countries. Data from a progeny test trial was available to identify parents with high breeding values for resistance. We successfully applied the GLIMMIX procedure for this purpose. Heritability for resistance was moderate, and we identified two parents and their family as having high potential for Psa resistance breeding. There are several potential pitfalls when using GLMMs with binary data and these are briefly discussed.

PollenCALC: software for estimation of pollen compatibility of self-incompatible allo- and autotetraploid species

Background

Self-incompatibility (SI) is a biological mechanism to avoid inbreeding in allogamous plants. In grasses, this mechanism is controlled by a two-locus system (S-Z). Calculation of male and female gamete frequencies is complex for tetraploid species. We are not aware of any software available for predicting pollen haplotype frequencies and pollen compatibility in tetraploid species.

Results

PollenCALC is a software tool written in C++ programming language that can predict pollen compatibility percentages for polyploid species with a two-locus (S, Z) self-incompatibility system. The program predicts pollen genotypes and frequencies based on defined meiotic parameters for allo- or autotetraploid species with a gametophytic S-Z SI system. These predictions can be used to obtain expected values for for diploid and for (allo- or autotetraploidy SI grasses.

Conclusion

The information provided by this calculator can be used to predict compatibility of pair-crosses in plant breeding applications, to analyze segregation distortion for S and Z genes, as well as linked markers in mapping populations, hypothesis testing of the number of S and Z alleles in a pair cross, and the underlying genetic model.

The validation of a 15 STR multiplex PCR for Cannabis species

Trade and acquisition of Cannabis drugs are illegal in many countries worldwide; nevertheless, crimes related with these drugs are a major problem for the investigative authorities. With this manuscript, we want to introduce a 15 short tandem repeat (STR) Cannabis marker set that can be amplified in one PCR reaction. This multiplex PCR is specific to Cannabis species and combines highly informative STR markers. The 15 STR multiplex is easy to use and was validated according to common laboratory quality standards. Due to the fact that a lot of Cannabis plants are cultivated by clonal propagation and may show aneuploidy, polyploidy or multiple gene loci, it is not possible to apply biostatistics that follow the Hardy-Weinberg law. However, this multiplex will help the police to trace back trade routes of drug syndicates or dealers and it can help to link Cannabis plants to a crime scene.

Interspecific and interploidal gene flow in Central European Arabidopsis (Brassicaceae)

BACKGROUND

Effects of polyploidisation on gene flow between natural populations are little known. Central European diploid and tetraploid populations of Arabidopsis arenosa and A. lyrata are here used to study interspecific and interploidal gene flow, using a combination of nuclear and plastid markers.

RESULTS

Ploidal levels were confirmed by flow cytometry. Network analyses clearly separated diploids according to species. Tetraploids and diploids were highly intermingled within species, and some tetraploids intermingled with the other species, as well. Isolation with migration analyses suggested interspecific introgression from tetraploid A. arenosa to tetraploid A. lyrata and vice versa, and some interploidal gene flow, which was unidirectional from diploid to tetraploid in A. arenosa and bidirectional in A. lyrata.

CONCLUSIONS

Interspecific genetic isolation at diploid level combined with introgression at tetraploid level indicates that polyploidy may buffer against negative consequences of interspecific hybridisation. The role of introgression in polyploid systems may, however, differ between plant species, and even within the small genus Arabidopsis, we find very different evolutionary fates when it comes to introgression.

Evaluating the utility of microsatellites for investigations of autopolyploid taxa

Autopolyploid taxa present numerous challenges for population genetic analyses due to difficulties determining allele dosage. Dosage ambiguity hinders accurate assessment of allele frequencies, multilocus genotypes (MLGTs), as well as levels and patterns of clonality. The pervasiveness of polyploidy in the evolutionary history of plant taxa makes this a recurring problem. Whereas diploidization of loci may occur over time, duplication of at least some loci is still frequently evident. Fortunately, with high-quality allozyme gels, it is possible to accurately infer allele dosage and, thus, determine exact MLGTs. However, accurately assessing dosage of microsatellite peaks is nearly impossible when studying wild populations with a large number of alleles per locus. Even if precise knowledge of genotypes is not required, for comparable numbers of alleles per locus and loci, the number of "phenotypes" is always lower with microsatellites than allozymes due to the inability to assess allele dosage. Microsatellite loci typically have more alleles per locus relative to allozymes although fewer loci are generally employed. Here, we present a mathematical model for comparing the relative utility of simple sequence repeat (SSR) versus allozyme markers to discriminate MLGTs. For example, the average plant allozyme study (2.6 alleles per locus, 10 polymorphic loci) has better discriminating power than SSR markers with 10 alleles at each of 3 loci, 9 alleles at 4 loci, 6 alleles at 5 loci, 5 alleles at 6 loci, and 4 alleles at 8 loci, demonstrating the value of assessing the relative discriminating power of these markers.

Characterization of microsatellite loci and reliable genotyping in a polyploid plant, Mercurialis perennis (Euphorbiaceae)

For many applications in population genetics, codominant simple sequence repeats (SSRs) may have substantial advantages over dominant anonymous markers such as amplified fragment length polymorphisms (AFLPs). In high polyploids, however, allele dosage of SSRs cannot easily be determined and alleles are not easily attributable to potentially diploidized loci. Here, we argue that SSRs may nonetheless be better than AFLPs for polyploid taxa if they are analyzed as effectively dominant markers because they are more reliable and more precise. We describe the transfer of SSRs developed for diploid Mercurialis huetii to the clonal dioecious M. perennis. Primers were tested on a set of 54 male and female plants from natural decaploid populations. Eight of 65 tested loci produced polymorphic fragments. Binary profiles from 4 different scoring routines were used to define multilocus lineages (MLLs). Allowing for fragment differences within 1 MLL, all analyses revealed the same 14 MLLs without conflicting with merigenet, sex, or plot assignment. For semiautomatic scoring, a combination of as few as 2 of the 4 most polymorphic loci resulted in unambiguous discrimination of clones. Our study demonstrates that microsatellite fingerprinting of polyploid plants is a cost efficient and reliable alternative to AFLPs, not least because fewer loci are required than for diploids.

High homologous gene conservation despite extreme autopolyploid redundancy in sugarcane

Modern sugarcane (Saccharum spp.) is the leading sugar crop and a primary energy crop. It has the highest level of 'vertical' redundancy (2n=12x=120) of all polyploid plants studied to date. It was produced about a century ago through hybridization between two autopolyploid species, namely S. officinarum and S. spontaneum. In order to investigate the genome dynamics in this highly polyploid context, we sequenced and compared seven hom(oe)ologous haplotypes (bacterial artificial chromosome clones). Our analysis revealed a high level of gene retention and colinearity, as well as high gene structure and sequence conservation, with an average sequence divergence of 4% for exons. Remarkably, all of the hom(oe)ologous genes were predicted as being functional (except for one gene fragment) and showed signs of evolving under purifying selection, with the exception of genes within segmental duplications. By contrast, transposable elements displayed a general absence of colinearity among hom(oe)ologous haplotypes and appeared to have undergone dynamic expansion in Saccharum, compared with sorghum, its close relative in the Andropogonea tribe. These results reinforce the general trend emerging from recent studies indicating the diverse and nuanced effect of polyploidy on genome dynamics.

Radiative evolution of polyploid races of the Iberian carnation Dianthus broteri (Caryophyllaceae)

*The micro-evolutionary mechanisms that drive large-scale radiations are not completely understood, partly because of a shortage of population-level studies aimed at identifying putative causes of rapid evolutionary change. The Dianthus broteri complex, representing the largest polyploid series known to date for any species in the genus (2x, 4x, 6x and 12x cytotypes), belongs to a lineage that was recently found to have diversified at unusually rapid rates. *We used a combination of genome sequencing (internal transcribed spacer (ITS), plus chloroplast DNA (cpDNA) regions trnH-psbA, psbA-trnK and trnK-matK) and amplified fragment length polymorphism (AFLP) fingerprinting in 25 populations to infer the evolutionary history of extant polyploid races. *The haplotype, ribotype and AFLP reconstructions showed a star-shaped arrangement suggesting a pattern of radiative evolution. The major, widespread haplotype occurred at all ploidy levels, whereas 20 minor haplotypes were restricted to single populations and cytotypes. In addition, AFLP analyses retrieved well-supported cytogeographic groups: six clades were clearly differentiated in terms of ploidy level and geography. Molecular data indicate that gene flow among different cytotypes is rare or nonexistent. *Our study supports a scenario of rapid diversification in carnations in which autopolyploidy and allopolyploidy, in interaction with geography and/or isolation, have played prominent roles.

Isolation of polymorphic microsatellite markers for the tetraploid Dipteryx odorata, an intensely exploited Amazonian tree species

Dipteryx odorata is an intensely exploited Amazonian tree legume. Microsatellite markers were developed to study the genetic structure, gene flow and reproductive biology of D. odorata. Eight highly polymorphic microsatellite markers were isolated from enriched repeat libraries screened for microsatellite repeats. An average of 16 alleles and 0.964 phenotype diversity per locus were found in 76 individuals from the Tapajos National Forest, in the state of Pará in the Brazilian Amazon.

Development and characterization of microsatellite markers for analysis of population differentiation in the tree legume Acacia koa (Fabaceae: Mimosoideae) in the Hawaiian Islands

The aim of this research was to develop and use microsatellite markers to characterize the high-value timber tree Acacia koa (koa), which is endemic to the Hawaiian Islands. Genomic DNA fragments of 300-1000 bp were cloned and sequenced following enrichment for microsatellite motifs by PCR using 7 oligonucleotide repeat primers in separate reactions. Among 96 sequences analyzed, 63 contained unique microsatellite motifs flanked by variable sequences. A dual PCR method involving a primer walking step was used to develop 15 primer pairs. Another 16 primer pairs were developed directly from the variable sequences on both sides of the microsatellite motifs. These 31 primer pairs were tested on 172 koa plants representing 11 populations collected from 4 of the major Hawaiian Islands. Nine of the primers that identified polymorphic microsatellite loci and 3 that detected unique alleles exclusively in some populations were used for genetic diversity studies of koa. Cluster analysis and multidimensional scaling of the allelic phenotype data revealed that koa from Kauai formed a distinct group separate from koa of the neighboring islands of Oahu, Maui, and Hawaii. The oldest of the four islands, Kauai, also had the most diverse populations of koa.

Genetic diversity and differentiation processes in the ploidy series of Olea europaea L.: a multiscale approach from subspecies to insular populations

Geographical isolation and polyploidization are central concepts in plant evolution. The hierarchical organization of archipelagos in this study provides a framework for testing the evolutionary consequences for polyploid taxa and populations occurring in isolation. Using amplified fragment length polymorphism and simple sequence repeat markers, we determined the genetic diversity and differentiation patterns at three levels of geographical isolation in Olea europaea: mainland-archipelagos, islands within an archipelago, and populations within an island. At the subspecies scale, the hexaploid ssp. maroccana (southwest Morocco) exhibited higher genetic diversity than the insular counterparts. In contrast, the tetraploid ssp. cerasiformis (Madeira) displayed values similar to those obtained for the diploid ssp. guanchica (Canary Islands). Geographical isolation was associated with a high genetic differentiation at this scale. In the Canarian archipelago, the stepping-stone model of differentiation suggested in a previous study was partially supported. Within the western lineage, an east-to-west differentiation pattern was confirmed. Conversely, the easternmost populations were more related to the mainland ssp. europaea than to the western guanchica lineage. Genetic diversity across the Canarian archipelago was significantly correlated with the date of the last volcanic activity in the area/island where each population occurs. At the island scale, this pattern was not confirmed in older islands (Tenerife and Madeira), where populations were genetically homogeneous. In contrast, founder effects resulted in low genetic diversity and marked genetic differentiation among populations of the youngest island, La Palma.