Genetic diversity and population structure in the endangered tree Hopea hainanensis (Dipterocarpaceae) on Hainan Island, China

Hopea hainanensis Merrill & Chun (Dipterocarpaceae) is an endangered tree species restricted to Hainan Island, China and a small part of Northern Vietnam. On Hainan Island, it is an important indicator species for tropical forests. However, because of its highly valued timber, H. hainanensis has suffered from overexploitation, leading to a sharp population decline. To facilitate the conservation of this species, genetic diversity and population structure were assessed using 12 SSR markers for 10 populations sampled across Hainan Island. Compared to non-threatened Hopea species, H. hainanensis exhibited reduced overall genetic diversity and increased population differentiation (AMOVA: F_ST = 0.23). Bayesian model-based clustering and principal coordinate analysis consistently assigned H. hainanensis individuals into three genetic groups, which were found to be widespread and overlapping geographically. A Mantel test found no correlation between genetic and geographical distances (r = 0.040, p = 0.418). The observed genetic structure suggests that long-distance gene flow occurred among H. hainanensis populations prior to habitat fragmentation. A recent population bottleneck was revealed, which may cause rapid loss of genetic diversity and increased differentiation across populations. Based on these findings, appropriate strategies for the long-term conservation of the endangered species H. hainanensis are proposed.

Contrasted histories of organelle and nuclear genomes underlying physiological diversification in a grass species

C₄ photosynthesis evolved multiple times independently in angiosperms, but most origins are relatively old so that the early events linked to photosynthetic diversification are blurred. The grass Alloteropsis semialata is an exception, as this species encompasses C₄ and non-C₄ populations. Using phylogenomics and population genomics, we infer the history of dispersal and secondary gene flow before, during and after photosynthetic divergence in A. semialata. We further analyse the genome composition of individuals with varied ploidy levels to establish the origins of polyploids in this species. Detailed organelle phylogenies indicate limited seed dispersal within the mountainous region of origin and the emergence of a C₄ lineage after dispersal to warmer areas of lower elevation. Nuclear genome analyses highlight repeated secondary gene flow. In particular, the nuclear genome associated with the C₄ phenotype was swept into a distantly related maternal lineage probably via unidirectional pollen flow. Multiple intraspecific allopolyploidy events mediated additional secondary genetic exchanges between photosynthetic types. Overall, our results show that limited dispersal and isolation allowed lineage divergence, with photosynthetic innovation happening after migration to new environments, and pollen-mediated gene flow led to the rapid spread of the derived C₄ physiology away from its region of origin.

Genome-wide genotyping elucidates the geographical diversification and dispersal of the polyploid and clonally propagated yam (Dioscorea alata)

BACKGROUND AND AIMS

Inferring the diffusion history of many human-dispersed species is still not straightforward due to unresolved past human migrations. The centre of diversification and routes of migration of the autopolyploid and clonally propagated greater yam, Dioscorea alata, one of the oldest edible tubers, remain unclear. Here, we address yam demographic and dispersal history using a worldwide sample.

METHODS

We characterized genome-wide patterns of genetic variation using genotyping by sequencing 643 greater yam accessions spanning four continents. First, we disentangled the polyploid and clonal components of yam diversity using allele frequency distribution and identity by descent approaches. We then addressed yam geographical origin and diffusion history with a model-based coalescent inferential approach.

KEY RESULTS

Diploid genotypes were more frequent than triploids and tetraploids worldwide. Genetic diversity was generally low and clonality appeared to be a main factor of diversification. The most likely evolutionary scenario supported an early divergence of mainland Southeast Asian and Pacific gene pools with continuous migration between them. The genetic make-up of triploids and tetraploids suggests that they have originated from these two regions before westward yam migration. The Indian Peninsula gene pool gave origin to the African gene pool, which was later introduced to the Caribbean region.

CONCLUSIONS

Our results are congruent with the hypothesis of independent domestication origins of the two main Asian and Pacific gene pools. The low genetic diversity and high clonality observed suggest a strong domestication bottleneck followed by thousands of years of widespread vegetative propagation and polyploidization. Both processes reduced the extent of diversity available for breeding, and this is likely to threaten future adaptation.

Conservation genomics and pollination biology of an endangered, edaphic-endemic, octoploid herb: El Dorado bedstraw (Galium californicum subsp. sierrae; Rubiaceae)

El Dorado bedstraw (Galium californicum subsp. sierrae) is a federally endangered dioecious, octoploid, perennial herb found only in the Pine Hill region of El Dorado County, CA, USA. Like many species of Galium, El Dorado bedstraw is capable of both sexual and asexual reproduction, spreading via stem-layering as well as seeds. El Dorado bedstraw is also dioecious, and thus dependent on pollinators to transfer pollen from male to female stems. The capacity for asexual reproduction has conservation implications for this plant, due to the potential for populations to become dominated by a small number of clones in the absence of recruitment from seeds. No previous work has examined either the population genetics or pollination biology of this plant. Here, double-digest restriction site-associated DNA sequencing was used to develop a genetic dataset for a sample of El Dorado bedstraw (12 individuals from each of seven locations). Genomic data was used to calculate population genetic statistics and quantify the degree to which clonality affects the sampled populations. Visual observation of insect visitors at every sampling location was used to assess the potential for pollen transfer within and among locations. A total of 23 clonal colonies were detected across 82 successfully sequenced stems, consisting of an average of 2.4 individuals (range: 2–6). Significant isolation by distance among locations was detected using a Mantel test. Insect pollinators were from eleven families, consisting mainly of small species with weak flight. It is recommended that clonality and small-scale population differentiation be taken into account in conservation measures.

Population Genetics of Polyploid Complex Psidium cattleyanum Sabine (Myrtaceae): Preliminary Analyses Based on New Species-Specific Microsatellite Loci and Extension to Other Species of the Genus

Polyploidy is a phenomenon that alters the genetic diversity of populations and has been reported as one of the most important evolutionary forces for plant diversification. The Psidium cattleyanum complex comprises a group of wild populations with several ploidy levels reported in the literature. The multiple cytotypes, associated with its wide distribution area, make this species a potential key model for understanding evolutionary processes related to polyploidization. In this study, we isolated and characterized nuclear microsatellite markers of P. cattleyanum and tested their transferability to other nine species of the genus. We performed a preliminary analysis of genetic diversity and population structure in three populations of P. cattleyanum. The three populations analyzed had different chromosome numbers, being polyploid cytotypes (2n = 6x = 66, 2n = 7x = 77 and 2n = 8x = 88). We designed 46 primer pairs and successfully amplified 37 markers, from which the 10 best were selected for analysis. Considering both the PIC and DP values, most of markers were highly informative. The new SSR markers were used to assess the levels of genetic diversity of the populations and detected one population with predominance of sexual reproduction. DAPC analysis pointed the formation of three groups, which corresponded to the populations analyzed. The markers were successfully amplified in related species, with some species presenting 80% transferability. By producing this panel of polymorphic microsatellites, we contribute to the understanding evolution in groups of natural polyploids for future studies.

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.

Genomic insight into the developmental history of southern highbush blueberry populations

Interspecific hybridization is a common breeding approach for introducing novel traits and genetic diversity to breeding populations. Southern highbush blueberry (SHB) is a blueberry cultivar group that has been intensively bred over the last 60 years. Specifically, it was developed by multiple interspecific crosses between northern highbush blueberry [NHB, Vaccinium corymbosum L. (2n = 4x = 48)] and low-chill Vaccinium species to expand the geographic limits of highbush blueberry production. In this study, we genotyped polyploid blueberries, including 105 SHB, 17 NHB, and 10 rabbiteye blueberry (RE) (Vaccinium virgatum Aiton), from the accessions planted at Poplarville, Mississippi, and accessions distributed in Japan, based on the double-digest restriction site-associated DNA sequencing. The genome-wide SNP data clearly indicated that RE cultivars were genetically distinct from SHB and NHB cultivars, whereas NHB and SHB were genetically indistinguishable. The population structure results appeared to reflect the differences in the allele selection strategies that breeders used for developing germplasm adapted to local climates. The genotype data implied that there are no or very few genomic segments that were commonly introgressed from low-chill Vaccinium species to the SHB genome. Principal component analysis-based outlier detection analysis found a few loci associated with a variable that could partially differentiate NHB and SHB. These SNP loci were detected in Mb-scale haplotype blocks and may be close to the functional genes related to SHB development. Collectively, the data generated in this study suggest a polygenic adaptation of SHB to the southern climate, and may be relevant for future population-scale genome-wide analyses of blueberry.

Secondary origin, hybridization and sexual reproduction in a diploid-tetraploid contact zone of the facultatively apomictic orchid Zygopetalum mackayi

The production of triploids and apomictic reproduction are important processes for polyploid establishment and cytotype coexistence, but we know little about the interaction between triploids and facultatively apomictic plants. To bridge this gap, we studied the pollen-dependent, facultatively apomictic orchid Zygopetalum mackayi from high-elevation outcrops of southeast Brazil. We described the nature of the contact between Z. mackayi cytotypes and patterns of genetic diversity and structure based on eight microsatellite markers and 155 individuals of pure tetraploid, pure diploid and mixed cytotype populations. Our results revealed high values of genetic and genotypic diversity within all populations of Z. mackayi. Each cytotype emerged as a genetic distinct cluster, combining individuals from different populations. Triploids clustered in an intermediate position between diploids and tetraploids. Most genetic variance is associated with individuals within populations and genetic differentiation is high among populations. Mixed cytotype populations of Z. mackayi originate from secondary contact. Triploids are hybrids between diploids and tetraploids and likely act as a bridge. Our results point to the predominance of sexual reproduction in all populations but do not corroborate previous basic chromosome number for this species. Polyploidy rather than facultative apomixis may explain the larger geographic distribution of tetraploids of Z. mackayi.

Phylogenomic Relationships and Evolution of Polyploid Salix Species Revealed by RAD Sequencing Data

Polyploidy is common in the genus Salix. However, little is known about the origin, parentage and genomic composition of polyploid species because of a lack of suitable molecular markers and analysis tools. We established a phylogenomic framework including species of all described sections of Eurasian shrub willows. We analyzed the genomic composition of seven polyploid willow species in comparison to putative diploid parental species to draw conclusions on their origin and the effects of backcrossing and post-origin evolution. We applied recently developed programs like SNAPP, HyDe, and SNiPloid to establish a bioinformatic pipeline for unravelling the complexity of polyploid genomes. RAD sequencing revealed 23,393 loci and 320,010 high quality SNPs for the analysis of relationships of 35 species of Eurasian shrub willows (Salix subg. Chamaetia/Vetrix). Polyploid willow species appear to be predominantly of allopolyploid origin. More ancient allopolyploidization events were observed for two hexaploid and one octoploid species, while our data suggested a more recent allopolyploid origin for the included tetraploids and identified putative parental taxa. SNiPloid analyses disentangled the different genomic signatures resulting from hybrid origin, backcrossing, and secondary post-origin evolution in the polyploid species. Our RAD sequencing data demonstrate that willow genomes are shaped by ancient and recent reticulate evolution, polyploidization, and post-origin divergence of 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.

Diving in uncharted waters: An updated genetics toolkit highlights the challenges of polyploidy in landscape genomics analyses

As molecular ecologists, we have by necessity become adept at working across computational platforms. A diverse community of scientists has developed a broad array of analytical resources spanning command line to graphical user interface across Linux, Mac, and Windows environments and a dizzying array of program-specific input formats. In light of this, we often explore our data like free divers - filling our lungs with air and descending for a short period of time into one part of our data set before resurfacing, reformatting, and preparing for our next analysis. In this issue of Molecular Ecology Resources, Meirmans (2020) presents an updated version of GenoDive, a program with a toolkit that provides users with the opportunity to stay a while and delve deeper into the diverse portfolio of information provided by a genomic data set. The comprehensive nature of GenoDive coupled with its unique capability to handle both diploid and polyploid data also provides an opportunity to reflect on the unevenness of resources available for the analysis of polyploid versus diploid data. Since new updates include the addition of plug-ins for genotype-environment association analyses, we limit the observations presented here to the common tools used for landscape genomics analyses.

The transferability of microsatellite loci from a homoploid to a polyploid hybrid complex: an example from fine-leaved Festuca species (Poaceae)

Background

Microsatellite loci, or single sequence repeats (SSR), are widely used as powerful markers in population genetics. They represent an attractive tool for studying plants such as grasses, whose evolution is driven by hybridisation and polyploidization. However, the development of microsatellite markers has been challenging and time-consuming, especially for non-model organisms lacking available genome-wide sequence data. One straightforward and low-cost approach is to transfer the SSR loci developed for one species, or complex, to another closely-related one. This work evaluates the transferability of microsatellite loci from homoploid to allopolyploid complexes of fine-leaved Festuca species and to assess their use in two new species. The studied complex (F. amethystina-F. tatrae) is a useful model for research on the local adaptability of grasses with different ploidy levels. Since both species can be considered as rare or threatened (F. tatrae-as a mountain and narrow endemic species and F. amethystina-a mountain species with relict lowland populations), any tool enabling studies on genetic diversity and population genetics, such as SSR markers, could also be very useful in a conservation context.

Methods

The ploidy level within populations was estimated using flow cytometry. One diploid and one tetraploid population of F. amethystina and a diploid population of F. tatrae were chosen to test the transferability of SSR loci. Because our work describes the transfer of SSR nuclear markers designed originally for F. gautieri, a phylogenetic tree was prepared based on the ITS marker to assess the genetic distance between the studied complexes. The PCR products were separated on a high-resolution agarose gel, intended for SSR marker analysis. Appropriate solutions for the allotetraploid population and whole mixed-ploidy complex were implemented.

Results

Flow cytometry confirmed earlier data regarding DNA content in the investigated species and cytotypes. The phylogenetic ITS tree indicated a small genetic distance between F. gautieri complexes and the studied species. Ten microsatellite markers were successfully transferred. All markers were polymorphic. In total, 163 different alleles were scored from the 10 SSR loci. PCoA of accessions revealed well-separated groups corresponding to studied populations. Over 60% of the total variance is explained by differentiation within populations and one third among them.

Conclusions

The transferred markers are valid tools for the study of population genetics and inheritance relationships within cytotypes and species and between them. The presented markers can be used to study inbreeding depression in the Festuca species, and variations in the degrees of genetic diversity between different cytotypes in mountain and lowland areas. Our findings can also be applied to study conservation strategies for ensuring biodiversity at the genetic level in polyploid complexes.

Genotyping-by-sequencing and ecological niche modeling illuminate phylogeography, admixture, and Pleistocene range dynamics in quaking aspen (Populus tremuloides)

Populus tremuloides is the widest-ranging tree species in North America and an ecologically important component of mesic forest ecosystems displaced by the Pleistocene glaciations. Using phylogeographic analyses of genome-wide SNPs (34,796 SNPs, 183 individuals) and ecological niche modeling, we inferred population structure, ploidy levels, admixture, and Pleistocene range dynamics of P. tremuloides, and tested several historical biogeographical hypotheses. We found three genetic lineages located mainly in coastal-Cascades (cluster 1), east-slope Cascades-Sierra Nevadas-Northern Rockies (cluster 2), and U.S. Rocky Mountains through southern Canadian (cluster 3) regions of the P. tremuloides range, with tree graph relationships of the form ((cluster 1, cluster 2), cluster 3). Populations consisted mainly of diploids (86%) but also small numbers of triploids (12%) and tetraploids (1%), and ploidy did not adversely affect our genetic inferences. The main vector of admixture was from cluster 3 into cluster 2, with the admixture zone trending northwest through the Rocky Mountains along a recognized phenotypic cline (Utah to Idaho). Clusters 1 and 2 provided strong support for the "stable-edge hypothesis" that unglaciated southwestern populations persisted in situ since the last glaciation. By contrast, despite a lack of clinal genetic variation, cluster 3 exhibited "trailing-edge" dynamics from niche suitability predictions signifying complete northward postglacial expansion. Results were also consistent with the "inland dispersal hypothesis" predicting postglacial assembly of Pacific Northwestern forest ecosystems, but rejected the hypothesis that Pacific-coastal populations were colonized during outburst flooding from glacial Lake Missoula. Overall, congruent patterns between our phylogeographic and ecological niche modeling results and fossil pollen data demonstrate complex mixtures of stable-edge, refugial locations, and postglacial expansion within P. tremuloides. These findings confirm and refine previous genetic studies, while strongly supporting a distinct Pacific-coastal genetic lineage of quaking aspen.

Molecular diversity and genetic structure of Saccharum complex accessions

Sugarcane is an important crop for food and energy security, providing sucrose and bioethanol from sugar content and bioelectricity from lignocellulosic bagasse. In order to evaluate the diversity and genetic structure of the Brazilian Panel of Sugarcane Genotypes (BPSG), a core collection composed by 254 accessions of the Saccharum complex, eight TRAP markers anchored in sucrose and lignin metabolism genes were evaluated. A total of 584 polymorphic fragments were identified and used to investigate the genetic structure of BPSG through analysis of molecular variance (AMOVA), principal components analysis (PCA), a Bayesian method using STRUCTURE software, genetic dissimilarity and phylogenetic tree. AMOVA showed a moderate genetic differentiation between ancestors and improved accessions, 0.14, and the molecular variance was higher within populations than among populations, with values of 86%, 95% and 97% when constrasting improved with ancestors, foreign with ancestors and improved with foreign, respectively. The PCA approach suggests clustering in according with evolutionary and Brazilian breeding sugarcane history, since improved accessions from older generations were positioned closer to ancestors than improved accessions from recent generations. This result was also confirmed by STRUCTURE analysis and phylogenetic tree. The Bayesian method was able to separate ancestors of the improved accessions while the phylogenetic tree showed clusters considering the family relatedness within three major clades; the first being composed mainly by ancestors and the other two mainly by improved accessions. This work can contribute to better management of the crosses considering functional regions of the sugarcane genome.

Patterns of genetic diversity in North Africa: Moroccan-Algerian genetic split in Juniperus thurifera subsp. africana

Juniperus thurifera is a key element of the forest communities in arid and semi-arid areas of the western Mediterranean. Previous genetic and morphological investigations suggested that Algerian populations are genetically more similar to European than to Moroccan populations and advocated their recognition at the variety rank. We aimed to investigate the spatial genetic structure in J. thurifera to verify the distinct character of the Algerian population in terms of the genetic breaks reported among several North African taxa. We also modelled species distributions since the Eemian to recognise the impact of past climatic changes on the current pattern of diversity and predict possible changes in species distribution in the future. Species-specific microsatellites were used in the analysis of 11 populations from Algeria, Morocco and Europe. We revealed the significant genetic distinctiveness of the Algerian populations from the Moroccan and European stands that may have important taxonomic and conservation implications. The diversity pattern revealed for J. thurifera reflects the east-west genetic splits reported among some North African plant and animal taxa and suggests an impact of shared historical processes. Additionally, modelling of the distribution allowed us to identify possible glacial refugia and their impact on the modern pattern of differentiation in J. thurifera. Reduction of species occurrence, especially in the European domain, is likely according to the future projections of the species distribution.

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.

Development and characterization of a duplex PCR assay in Chinese sturgeon (Acipenser sinensis) for genetic analysis

Chinese sturgeon (Acipenser sinensis) has been listed as a critically endangered species on the IUCN Red List and is an endemic fish of China. Five sets of duplex polymerase chain reactions (PCR) assays were developed with 10 tetranucleotide microsatellites for Chinese sturgeon. The size of CS57, ZHX43, ZHX69, AS105, ZHX51, AS074, ZHX2, AS078, AS026 and AS073 products in 184 Chinese sturgeon individuals ranged from 257-305, 191-241, 251-285, 172-244, 236-260, 169-209, 194-234, 92-176, 165-257 and 120-164, respectively. The observed allele number of the 10 microsatellites ranged from 7 to 16, and the total number of alleles was 106. The number of alleles per individual in CS57, ZHX43, AS105, AS074, AS078 and AS026 was 1-4. The number of alleles per individual in ZHX69, ZHX51, ZHX2 and AS073 was 2-4. The mean number of alleles per locus per individual ranged from 2.01-3.76. The expected heterozygosity (H_E), observed heterozygosity (H_O), polymorphic information content (PIC) and Shannon-Weiner diversity index (H') ranged from 0.582 to 0.899, from 0.676 to 1, from 0.518 to 0.886 and from 1.034 to 2.34, respectively. Despite many advantages, the use of microsatellites as genetic analysis tools can be limited by the cost of the associated experiment. To solve this problem, this set of five duplex PCRs will provide tools that are more helpful, less expensive and less time consuming than others used for genetic analyses in Chinese sturgeon.

So Closely Related and Yet So Different: Strong Contrasts Between the Evolutionary Histories of Species of the Cardamine pratensis Polyploid Complex in Central Europe

Recurrent polyploid formation and weak reproductive barriers between independent polyploid lineages generate intricate species complexes with high diversity and reticulate evolutionary history. Uncovering the evolutionary processes that formed their present-day cytotypic and genetic structure is a challenging task. We studied the species complex of Cardamine pratensis, composed of diploid endemics in the European Mediterranean and diploid-polyploid lineages more widely distributed across Europe, focusing on the poorly understood variation in Central Europe. To elucidate the evolution of Central European populations we analyzed ploidy level and genome size variation, genetic patterns inferred from microsatellite markers and target enrichment of low-copy nuclear genes (Hyb-Seq), and environmental niche differentiation. We observed almost continuous variation in chromosome numbers and genome size in C. pratensis s.str., which is caused by the co-occurrence of euploid and dysploid cytotypes, along with aneuploids, and is likely accompanied by inter-cytotype mating. We inferred that the polyploid cytotypes of C. pratensis s.str. are both of single and multiple, spatially and temporally recurrent origins. The tetraploid Cardamine majovskyi evolved at least twice in different regions by autopolyploidy from diploid Cardamine matthioli. The extensive genome size and genetic variation of Cardamine rivularis reflects differentiation induced by the geographic isolation of disjunct populations, establishment of triploids of different origins, and hybridization with sympatric C. matthioli. Geographically structured genetic lineages identified in the species under study, which are also ecologically divergent, are interpreted as descendants from different source populations in multiple glacial refugia. The postglacial range expansion was accompanied by substantial genetic admixture between the lineages of C. pratensis s.str., which is reflected by diffuse borders in their contact zones. In conclusion, we identified an interplay of diverse processes that have driven the evolution of the species studied, including allopatric and ecological divergence, hybridization, multiple polyploid origins, and genetic reshuffling caused by Pleistocene climate-induced range dynamics.

Phylogenomic Relationships of Diploids and the Origins of Allotetraploids in Dactylorhiza (Orchidaceae)

Disentangling phylogenetic relationships proves challenging for groups that have evolved recently, especially if there is ongoing reticulation. Although they are in most cases immediately isolated from diploid relatives, sets of sibling allopolyploids often hybridize with each other, thereby increasing the complexity of an already challenging situation. Dactylorhiza (Orchidaceae: Orchidinae) is a genus much affected by allopolyploid speciation and reticulate phylogenetic relationships. Here, we use genetic variation at tens of thousands of genomic positions to unravel the convoluted evolutionary history of Dactylorhiza. We first investigate circumscription and relationships of diploid species in the genus using coalescent and maximum likelihood methods, and then group 16 allotetraploids by maximum affiliation to their putative parental diploids, implementing a method based on genotype likelihoods. The direction of hybrid crosses is inferred for each allotetraploid using information from maternally inherited plastid RADseq loci. Starting from age estimates of parental taxa, the relative ages of these allotetraploid entities are inferred by quantifying their genetic similarity to the diploids and numbers of private alleles compared with sibling allotetraploids. Whereas northwestern Europe is dominated by young allotetraploids of postglacial origins, comparatively older allotetraploids are distributed further south, where climatic conditions remained relatively stable during the Pleistocene glaciations. Our bioinformatics approach should prove effective for the study of other naturally occurring, nonmodel, polyploid plant complexes.

Two-Locus Local Adaptation by Additive or Epistatic Gene Combinations in Autotetraploids Versus Diploids

In this article, we present a theoretical comparison of local adaptation between diploid and autotetraploid populations when fitness is determined by either additive or epistatic interactions between alleles at 2 loci. A continent-island model of local adaptation is derived, with 1-way migration from the continent to the island and distinct genotypes adaptive on the continent versus the island. The meiotic component of the model accounts for multivalent formation and the processes of chromosomal gametic disequilibrium and double reduction, which are unique to autotetraploids. Both the adaptability and efficiency of adaptation are investigated, where adaptability asks whether a population adapts and efficiency is the rate of adaptation. With an additive genetic basis to fitness, diploids experience better adaptability and efficiency than autotetraploids. With epistasis, our results indicate a limited parameter space in which autotetraploids have greater adaptability than diploids, but results indicate an interesting difference between adaptability and efficiency of adaptation. Oftentimes, diploids exhibit greater adaptability whereas autotetraploids exhibit greater efficiency of adaptation. These findings provide evidence for the advantage of epistasis within autotetraploids when efficiency of adaptation is of interest. Although autotetraploids are more efficient, under the same conditions and at equilibrium, diploid populations often have higher mean local fitness. Overall, the most ideal situation for autotetraploid local adaptation compared to diploids is when epistasis is strong, mutation is weak, recombination is high, selection is strong, deleterious selection is additive, chromosomal gametic disequilibrium is present, and double reduction is absent.

Allelic Variants for Candidate Nitrogen Fixation Genes Revealed by Sequencing in Red Clover (Trifolium pratense L.)

Plant–rhizobia symbiosis can activate key genes involved in regulating nodulation associated with biological nitrogen fixation (BNF). Although the general molecular basis of the BNF process is frequently studied, little is known about its intraspecific variability and the characteristics of its allelic variants. This study’s main goals were to describe phenotypic and genotypic variation in the context of nitrogen fixation in red clover (Trifolium pretense L.) and identify variants in BNF candidate genes associated with BNF efficiency. Acetylene reduction assay validation was the criterion for selecting individual plants with particular BNF rates. Sequences in 86 key candidate genes were obtained by hybridization-based sequence capture target enrichment of plants with alternative phenotypes for nitrogen fixation. Two genes associated with BNF were identified: ethylene response factor required for nodule differentiation (EFD) and molybdate transporter 1 (MOT1). In addition, whole-genome population genotyping by double-digest restriction-site-associated sequencing (ddRADseq) was performed, and BNF was evaluated by the natural ¹⁵N abundance method. Polymorphisms associated with BNF and reflecting phenotype variability were identified. The genetic structure of plant accessions was not linked to BNF rate of measured plants. Knowledge of the genetic variation within BNF candidate genes and the characteristics of genetic variants will be beneficial in molecular diagnostics and breeding of red clover.

Genetic diversity and structure of baobab (Adansonia digitata L.) in southeastern Kenya

Baobab (Adansonia digitata L.) is an iconic tree of African savannahs. Its multipurpose character and nutritional composition of fruits and leaves offer high economic and social potential for local communities. There is an urgent need to characterize the genetic diversity of the Kenyan baobab populations in order to facilitate further conservation and domestication programmes. This study aims at documenting the genetic diversity and structure of baobab populations in southeastern Kenya. Leaf or bark samples were collected from 189 baobab trees in seven populations distributed in two geographical groups, i.e. four inland and three coastal populations. Nine microsatellite loci were used to assess genetic diversity. Overall, genetic diversity of the species was high and similarly distributed over the populations. Bayesian clustering and principal coordinate analysis congruently divided the populations into two distinct clusters, suggesting significant differences between inland and coastal populations. The genetic differentiation between coastal and inland populations suggests a limited possibility of gene flow between these populations. Further conservation and domestications studies should take into consideration thegeographical origin of trees and more attention should be paid to morphological characterization of fruits and leaves of the coastal and inland populations to understand the causes and the impact of the differentiation.

Review: Molecular techniques to assess genetic variation within and between Panax ginseng and Panax quinquefolius

A variety of methods have been used to examine genetic differences in P. ginseng and P. quinquefolius. They have shown genetic differences within populations of P. ginseng (within and between elite cultivars, landraces and wild accessions), within populations of P. quinquefolius (within and between wild and cultivated accessions) and between P. ginseng and P. quinquefolius as well as other Panax species. Some examples of their applications have been to show that some elite cultivars are not uniform, there are possible founder effects in certain populations, there has been the spread of cultivated types into wild populations, relative diversity differs between different populations and identification of the source and purity of commercial samples. More work in the use of molecular markers for ginseng are needed, however, particularly the use of Next Generation Sequencing. Potential applications are the use of sequence analysis for genetic selection, breeding to develop new cultivars and providing traceability from field to consumer. Research on molecular markers in ginseng has lagged compared to other crops probably because of less of an emphasis on breeding for cultivar development and relatively small areas of production. The many potential benefits for ginseng production have yet to be realized.

Selfing ability and drift load evolve with range expansion

Colonization at expanding range edges often involves few founders, reducing effective population size. This process can promote the evolution of self-fertilization, but implicating historical processes as drivers of trait evolution is often difficult and requires an explicit model of biogeographic history. In plants, contemporary limits to outcrossing are often invoked as evolutionary drivers of self-fertilization, but historical expansions may shape mating system diversity, with leading-edge populations evolving elevated selfing ability. In a widespread plant, Campanula americana, we identified a glacial refugium in the southern Appalachian Mountains from spatial patterns of genetic drift among 24 populations. Populations farther from this refugium have smaller effective sizes and fewer rare alleles. They also displayed elevated heterosis in among-population crosses, reflecting the accumulation of deleterious mutations during range expansion. Although populations with elevated heterosis had reduced segregating mutation load, the magnitude of inbreeding depression lacked geographic pattern. The ability to self-fertilize was strongly positively correlated with the distance from the refugium and mutation accumulation-a pattern that contrasts sharply with contemporary mate and pollinator limitation. In this and other species, diversity in sexual systems may reflect the legacy of evolution in small, colonizing populations, with little or no relation to the ecology of modern populations.

Multiple auto- and allopolyploidisations marked the Pleistocene history of the widespread Eurasian steppe plant Astragalus onobrychis (Fabaceae)

The Eurasian steppes occupy a significant portion of the worldwide land surface and their biota have been affected by specific past range dynamics driven by ice ages-related climatic fluctuations. The dynamic alterations in conditions during the Pleistocene often triggered reticulate evolution and whole genome duplication events. Employing genomic, genetic and cytogenetic tools as well as morphometry we investigate the intricate evolution of Astragalus onobrychis, a widespread Eurasian steppe plant with diploid, tetraploid and octoploid cytotypes. To analyse the heteroploid RADseq dataset we employ both genotype-based and genotype-free methods that result in highly consistent results, and complement our inference with information from the plastid ycf1 region. We uncover a complex and reticulate evolutionary history, including at least one auto-tetraploidization event and two allo-octoploidization events; one of them involved also genetic contributions from other species, most likely A. goktschaicus. The present genetic structure points to the existence of four main clades within A. onobrychis, which only partly correspond to different ploidies. Time-calibrated diffusion models suggest that diversification within A. onobrychis was associated with ice age-related climatic fluctuations during the last million years. We finally argue for the usefulness of uniparentally inherited plastid markers, even in the genomic era, especially when investigating heteroploid systems.

STRUCTURE is more robust than other clustering methods in simulated mixed-ploidy populations

Analysis of population genetic structure has become a standard approach in population genetics. In polyploid complexes, clustering analyses can elucidate the origin of polyploid populations and patterns of admixture between different cytotypes. However, combining diploid and polyploid data can theoretically lead to biased inference with (artefactual) clustering by ploidy. We used simulated mixed-ploidy (diploid-autotetraploid) data to systematically compare the performance of k-means clustering and the model-based clustering methods implemented in STRUCTURE, ADMIXTURE, FASTSTRUCTURE and INSTRUCT under different scenarios of differentiation and with different marker types. Under scenarios of strong population differentiation, the tested applications performed equally well. However, when population differentiation was weak, STRUCTURE was the only method that allowed unbiased inference with markers with limited genotypic information (co-dominant markers with unknown dosage or dominant markers). Still, since STRUCTURE was comparatively slow, the much faster but less powerful FASTSTRUCTURE provides a reasonable alternative for large datasets. Finally, although bias makes k-means clustering unsuitable for markers with incomplete genotype information, for large numbers of loci (>1000) with known dosage k-means clustering was superior to FASTSTRUCTURE in terms of power and speed. We conclude that STRUCTURE is the most robust method for the analysis of genetic structure in mixed-ploidy populations, although alternative methods should be considered under some specific conditions.

Dissimilarity of individual microsatellite profiles under different mutation models: Empirical approach

Microsatellites (simple sequence repeats, SSRs) still remain popular molecular markers for studying neutral genetic variation. Two alternative models outline how new microsatellite alleles evolve. Infinite alleles model (IAM) assumes that all possible alleles are equally likely to result from a mutation, while stepwise mutation model (SMM) describes microsatellite evolution as stepwise adding or subtracting single repeat units. Genetic relationships between individuals can be analyzed in higher precision when assuming the SMM scenario with allele size differences as a proxy of genetic distance. If population structure is not predetermined in advance, an empirical data analysis usually includes (a) estimating proximity between individual SSR profiles with a selected dissimilarity measure and (b) determining putative genetic structure of a given set of individuals using methods of clustering and/or ordination for the obtained dissimilarity matrix. We developed new dissimilarity indices between SSR profiles of haploid, diploid, or polyploid organisms assuming different mutation models and compared the performance of these indices for determining genetic structure with population data and with simulations. More specifically, we compared SMM with a constant or variable mutation rate at different SSR loci to IAM using data from natural populations of a freshwater bryozoan Cristatella mucedo (diploid), wheat leaf rust Puccinia triticina (dikaryon), and wheat powdery mildew Blumeria graminis (monokaryon). We show that inferences about population genetic structure are sensitive to the assumed mutation model. With simulations, we found that Bruvo's distance performs generally poorly, while the new metrics are capturing the differences in the genetic structure of the populations.

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

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

pSBVB: A Versatile Simulation Tool To Evaluate Genomic Selection in Polyploid Species

Genomic Selection (GS) is the procedure whereby molecular information is used to predict complex phenotypes and it is standard in many animal and plant breeding schemes. However, only a small number of studies have been reported in horticultural crops, and in polyploid species in particular. In this paper, we have developed a versatile forward simulation tool, called polyploid Sequence Based Virtual Breeding (pSBVB), to evaluate GS strategies in polyploids; pSBVB is an efficient gene dropping software that can simulate any number of complex phenotypes, allowing a very flexible modeling of phenotypes suited to polyploids. As input, it takes genotype data from the founder population, which can vary from single nucleotide polymorphisms (SNP) chips up to sequence, a list of causal variants for every trait and their heritabilities, and the pedigree. Recombination rates between homeologous chromosomes can be specified, so that both allo- and autopolyploid species can be considered. The program outputs phenotype and genotype data for all individuals in the pedigree. Optionally, it can produce several genomic relationship matrices that consider exact or approximate genotype values. pSBVB can therefore be used to evaluate GS strategies in polyploid species (say varying SNP density, genetic architecture or population size, among other factors), or to optimize experimental designs for association studies. We illustrate pSBVB with SNP data from tetraploid potato and partial sequence data from octoploid strawberry, and we show that GS is a promising breeding strategy for polyploid species but that the actual advantage critically depends on the underlying genetic architecture. Source code, examples and a complete manual are freely available in GitHub https://github.com/lauzingaretti/pSBVB.

Sexual intraspecific recombination but not de novo origin governs the genesis of new apomictic genotypes in Potentilla puberula (Rosaceae)

Apomixis - asexual reproduction via seeds - might arise de novo following polyploidisation events, or via reproductive transfer of apomixis. Both processes can be obtained within species or via hybridisation. We aimed to determine the origin of apomictic genotypes in Potentilla puberula, a rosaceous species showing reproductive differentiation with ploidy: sexual tetraploids and apomictic penta- to octoploids, which regularly co-occur in sympatry. The study is based on 726 individuals, comprising all cytotypes, collected from 138 populations in the Eastern European Alps. We established relationships of cytotypes based on AFLP fingerprinting and cpDNA sequencing to test (1) whether the apomicts are of recurrent allopolyploid origin or originated from within the species via autopolyploidy, and (2) whether there are indications for reproductive transfer versus de novo origin of apomixis. Three principal pathways were identified which explain the origin of new apomictic genotypes, all involving at least one apomictic parent and thus compatible with the idea of reproductive transfer of the apomictic trait to the progeny: (1) self-fertilisation of unreduced egg cells in apomicts; (2) cross-fertilisation among apomicts; and (3) occasionally, heteroploid crosses among sexuals and apomicts. Autopolyploids derived from tetraploid sexuals were repeatedly observed, but did not express apomixis. Finally, our results suggest no role of other species in the origin of extant apomictic genotypes of P. puberula, although local hybrids with P. crantzii were identified. In conclusion, our results show that the formation of new apomictic genotypes required a genetic contribution from at least one apomictic parent. This finding is in accordance with the idea that apomixis is inheritable in P. puberula. On the contrary, lack of apomixis in penta- and hexaploids derived from sexual backgrounds did not support the hypothesis of a de novo origin of apomixis. Relatively high frequency of remnant sexuality in the apomicts involving different cytological pathways of seed formation can explain their high cytological and genotypic diversity. Finally, lack of global introgression from a third taxon is in support of P. puberula as a concise, although highly diverse, species.

Integrating Networks, Phylogenomics, and Population Genomics for the Study of Polyploidy

Duplication events are regarded as sources of evolutionary novelty, but our understanding of general trends for the long-term trajectory of additional genomic material is still lacking. Organisms with a history of whole genome duplication (WGD) offer a unique opportunity to study potential trends in the context of gene retention and/or loss, gene and network dosage, and changes in gene expression. In this review, we discuss the prevalence of polyploidy across the tree of life, followed by an overview of studies investigating genome evolution and gene expression. We then provide an overview of methods in network biology, phylogenomics, and population genomics that are critical for advancing our understanding of evolution post-WGD, highlighting the need for models that can accommodate polyploids. Finally, we close with a brief note on the importance of random processes in the evolution of polyploids with respect to neutral versus selective forces, ancestral polymorphisms, and the formation of autopolyploids versus allopolyploids.

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.