Recurrent allopolyploidy and its implications for conservation in vascular plants: a commentary on 'Population genomics of the Isoetes appalachiana (Isoetaceae) complex supports a "diploids-first" approach to conservation'

This article comments on:

David Wickell, Jacob Landis, Elizabeth Zimmer and Fay-Wei Li, Population genomics of the Isoetes appalachiana (Isoetaceae) complex supports a ‘diploids-first’ approach to conservation, Annals of Botany, Volume 133, Issue 2, 01 February 2024, Pages 261–272, https://doi.org/10.1093/aob/mcad180

Divergent mechanisms of reduced growth performance in Betula ermanii saplings from high-altitude and low-latitude range edges

The reduced growth performance of individuals from range edges is a common phenomenon in various taxa, and considered to be an evolutionary factor that limits the species' range. However, most studies did not distinguish between two mechanisms that can lead to this reduction: genetic load and adaptive selection to harsh conditions. To address this lack of understanding, we investigated the climatic and genetic factors underlying the growth performance of Betula ermanii saplings transplanted from 11 populations including high-altitude edge and low-latitude edge population. We estimated the climatic position of the populations within the overall B. ermanii's distribution, and the genetic composition and diversity using restriction-site associated DNA sequencing, and measured survival, growth rates and individual size of the saplings. The high-altitude edge population (APW) was located below the 95% significance interval for the mean annual temperature range, but did not show any distinctive genetic characteristics. In contrast, the low-latitude edge population (SHK) exhibited a high level of linkage disequilibrium, low genetic diversity, a distinct genetic composition from the other populations, and a high relatedness coefficient. Both APW and SHK saplings displayed lower survival rates, heights and diameters, while SHK saplings also exhibited lower growth rates than the other populations' saplings. The low heights and diameters of APW saplings was likely the result of adaptive selection to harsh conditions, while the low survival and growth rates of SHK saplings was likely the result of genetic load. Our findings shed light on the mechanisms underlying the reduced growth performance of range-edge populations.

The Himalayan uplift and evolution of aquatic biodiversity across Asia: Snowtrout (Cyprininae: Schizothorax) as a test case

Global biodiversity hotspots are often remote, tectonically active areas undergoing climatic fluctuations, such as the Himalaya Mountains and neighboring Qinghai-Tibetan Plateau (QTP). They provide biogeographic templates upon which endemic biodiversity can be mapped to infer diversification scenarios. Yet, this process can be somewhat opaque for the Himalaya, given substantial data gaps separating eastern and western regions. To help clarify, we evaluated phylogeographic and phylogenetic hypotheses for a widespread fish (Snowtrout: Cyprininae; Schizothorax) by sequencing 1,140 base pair of mtDNA cytochrome-b (cytb) from Central Himalaya samples (Nepal: N = 53; Bhutan: N = 19), augmented with 68 GenBank sequences (N = 60 Schizothorax/N = 8 outgroups). Genealogical relationships (N = 132) were analyzed via maximum likelihood (ML), Bayesian (BA), and haplotype network clustering, with clade divergence estimated via TimeTree. Snowtrout seemingly originated in Central Asia, dispersed across the QTP, then into Bhutan via southward-flowing tributaries of the east-flowing Yarlung-Tsangpo River (YLTR). Headwaters of five large Asian rivers provided dispersal corridors from Central into eastern/southeastern Asia. South of the Himalaya, the YLTR transitions into the Brahmaputra River, facilitating successive westward colonization of Himalayan drainages first in Bhutan, then Nepal, followed by far-western drainages subsequently captured by the (now) westward-flowing Indus River. Two distinct Bhutanese phylogenetic groups were recovered: Bhutan-1 (with three subclades) seemingly represents southward dispersal from the QTP; Bhutan-2 apparently illustrates northward colonization from the Lower Brahmaputra. The close phylogenetic/phylogeographic relationships between the Indus River (Pakistan) and western tributaries of the Upper Ganges (India/Nepal) potentially implicate an historic, now disjunct connection. Greater species-divergences occurred across rather than within-basins, suggesting vicariance as a driver. The Himalaya is a component of the Earth's largest glacial reservoir (i.e., the "third-pole") separate from the Arctic/Antarctic. Its unique aquatic biodiversity must be defined and conserved through broad, trans-national collaborations. Our study provides an initial baseline for this process.

Untangling the hedge: Genetic diversity in clonally and sexually transmitted genomes of European wild roses, Rosa L

While European wild roses are abundant and widely distributed, their morphological taxonomy is complicated and ambiguous. In particular, the polyploid Rosa section Caninae (dogroses) is characterised by its unusual meiosis, causing simultaneous clonal and sexual transmission of sub-genomes. This hemisexual reproduction, which often co-occurs with vegetative reproduction, defies the standard definition of species boundaries. We analysed seven highly polymorphic microsatellite loci, scored for over 2 600 Rosa samples of differing ploidy, collected across Europe within three independent research projects. Based on their morphology, these samples had been identified as belonging to 21 dogrose and five other native rose species. We quantified the degree of clonality within species and at individual sampling sites. We then compared the genetic structure within our data to current rose morpho-systematics and searched for hemisexually co-inherited sets of alleles at individual loci. We found considerably fewer copies of identical multi-locus genotypes in dogroses than in roses with regular meiosis, with some variation recorded among species. While clonality showed no detectable geographic pattern, some genotypes appeared to be more widespread. Microsatellite data confirmed the current classification of subsections, but they did not support most of the generally accepted dogrose microspecies. Under canina meiosis, we found co-inherited sets of alleles as expected, but could not distinguish between sexually and clonally inherited sub-genomes, with only some of the detected allele combinations being lineage-specific.

Genetic Covariation Between the Vertically Transmitted Endophyte Epichloë canadensis and Its Host Canada Wildrye

Symbiotic mutualisms are thought to be stabilized by correlations between the interacting genotypes which may be strengthened via vertical transmission and/or reduced genetic variability within each species. Vertical transmission, however, may weaken interactions over time as the endosymbionts would acquire mutations that could not be purged. Additionally, temporal variation in a conditional mutualism could create genetic variation and increased variation in the interaction outcome. In this study, we assessed genetic variation in both members of a symbiosis, the endosymbiotic fungal endophyte Epichloë canadensis and its grass host Canada wildrye (Elymus canadensis). Both species exhibited comparable levels of diversity, mostly within populations rather than between. There were significant differences between populations, although not in the same pattern for the two species, and the differences were not correlated with geographic distance for either species. Interindividual genetic distance matrices for the two species were significantly correlated, although all combinations of discriminant analysis of principle components (DAPC) defined multilocus genotype groups were found suggesting that strict genotype matching is not necessary. Variation in interaction outcome is common in grass/endophyte interactions, and our results suggest that the accumulation of mutations overtime combined with temporal variation in selection pressures increasing genetic variation in the symbiosis may be the cause.

Double reduction estimation and equilibrium tests in natural autopolyploid populations

Many bioinformatics pipelines include tests for equilibrium. Tests for diploids are well studied and widely available, but extending these approaches to autopolyploids is hampered by the presence of double reduction, the comigration of sister chromatid segments into the same gamete during meiosis. Though a hindrance for equilibrium tests, double reduction rates are quantities of interest in their own right, as they provide insights about the meiotic behavior of autopolyploid organisms. Here, we develop procedures to (i) test for equilibrium while accounting for double reduction, and (ii) estimate the double reduction rate given equilibrium. To do so, we take two approaches: a likelihood approach, and a novel U-statistic minimization approach that we show generalizes the classical equilibrium χ2 test in diploids. For small sample sizes and uncertain genotypes, we further develop a bootstrap procedure based on our U-statistic to test for equilibrium. We validate our methods on both simulated and real data.

Meeting in Liguria: hybridisation between Apennine endemic Euphorbia barrelieri and western Mediterranean E. nicaeensis led to the allopolyploid origin of E. ligustica

The Mediterranean Basin is renowned for its extremely rich biota and is considered as one of the 25 Global Biodiversity Hotspots, but its diversity is not homogeneously distributed. Outstanding in the number of (endemic) species are the Ligurian Alps (Italy). At the foot of the Ligurian Alps, little above the Mediterranean Sea, a disjunct occurrence of Italian endemic Euphorbia barrelieri was reported. Using an array of integrative methods ranging from cytogenetic (chromosome number and relative genome size estimation), over phylogenetic approaches (plastid, ITS and RAD sequencing) to multivariate morphometrics we disentangled the origin of these populations that were shown to be tetraploid. We performed phylogenetic analyses of the nuclear ITS and plastid regions of a broad taxonomic sampling of Euphorbia sect. Pithyusa to identify possible species involved in the origin of the tetraploid populations and then applied various analyses of RADseq data to identify the putative parental species. Our results have shown that the Ligurian populations of E. barrelieri are of allotetraploid origin that involved E. barrelieri and western Mediterranean E. nicaeensis as parental species. We thus describe a new species, E. ligustica, and hypothesise that its adaptation to similar environments in which E. barrelieri occurs, triggered development of similar morphology, whereas its genetic composition appears to be closer to that of E. nicaeensis. Our study emphasises the importance of polyploidisation for plant diversification, highlights the value of the Ligurian Alps as a hotspot of biodiversity and endemism and underlines the importance of integrative taxonomic approaches in uncovering cryptic diversity.

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

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

Population Genomic Analysis of Diploid-Autopolyploid Species

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

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

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

De Novo Transcriptome Assembly and EST-SSR Marker Development and Application in Chrysosplenium macrophyllum

Chrysosplenium macrophyllum Oliv., belonging to the family Saxifragaceae, is a traditional and unique Chinese herbal medicine. However, the lack of adequate molecular markers has hampered the progress regarding population genetics and evolution within this species. In this research, we used the DNBSEQ-T7 Sequencer (MGI) sequencing assay to analyze the transcriptome profiles of C. macrophyllum. SSR markers were developed on the basis of transcriptomic sequences and further validated on C. macrophyllum and other Chrysosplenium species. The genetic diversity and structure of the 12 populations were analyzed by using polymorphic expressed sequence tag simple sequence repeat (EST-SSR) markers. A potential pool of 3127 non-redundant EST-SSR markers were identified for C. macrophyllum in this study. The developed EST-SSR markers had high amplification rates and cross-species transferability in Chrysosplenium. Our results also showed that the natural populations of C. macrophyllum had a high level of genetic diversity. Genetic distance, principal component analysis, and popular structure analysis revealed that all 60 samples clustered into two major groups that were consistent with their geographical origins. This study provided a batch of highly polymorphic EST-SSR molecular markers that were developed via transcriptome sequencing. These markers will be of great significance for the study of the genetic diversity and evolutionary history of C. macrophyllum and other Chrysosplenium species.

A first look at sea-lavenders genomics - can genome wide SNP information tip the scales of controversy in the Limonium vulgare species complex?

BACKGROUND

Sea-lavenders (Limonium Mill., Plumbaginaceae) are a cosmopolitan group of diploid and polyploid plants often adapted to extreme saline environments, with a mostly Tethyan distribution, occurring in the Mediterranean, Irano-Turanian, Euro-Siberian and in the New World. The halophylic Limonium vulgare polyploid complex in particular, presents a large distribution throughout extreme salt-marsh habitats and shows little morphological but high taximetric variation, frequently blurring species delimitation. In this work we pursue three main goals: assert whether SNP data from polyploid individuals has the resolution to distinguish the seven sampled species, to better understand how genetically structured Limonium vulgare is, and attempt to identify specific molecular mechanisms for the differentiation between L. maritimum and L. vulgare. For this purpose, 95 individuals were genotyped using Genotyping by Sequencing (GBS), which were assembled as two independent datasets using IPYRAD. All analyses performed downstream of assembly were fully automated. Phylogenetic inference, PCA, and admixture plots were used to infer answers to the study's main goals.

RESULTS

Close to 10,000 SNPs were obtained for each dataset. Phylogenetic analyses reveal that polyploid data can be used to infer species relationships. Population structure analyses suggest a genetically structured L. vulgare. A set of 34 SNPs were found to be fully segregated between L. vulgare and L. maritimum, two of which are potentially linked to proteins that might be involved in the speciation process.

CONCLUSION

Despite polyploid data analyses shortcomings, GBS generated SNPs have the resolution to discern all seven included species. Limonium vulgare revealed pronounced genetic structure along a geographical north-south cline. L. maritimum always appears as a distinct genetic entity. Segregated SNPs between L. vulgare and L. maritimum indicate salinity response and morphological trait control genes as potentially interesting to follow up for studying these species' divergence process.

Population Structure of White Sturgeon (Acipenser transmontanus) in the Columbia River Inferred from Single-Nucleotide Polymorphisms

White sturgeon (Acipenser transmontanus) are the largest freshwater fish in North America, with reproducing populations in the Sacramento-San Joaquin, Fraser, and Columbia River Basins. Of these, the Columbia River is the largest, but it is also highly fragmented by hydroelectric dams, and many segments are characterized by declining abundance and persistent recruitment failure. Efforts to conserve and supplement these fish requires an understanding of their spatial genetic structure. Here, we assembled a large set of samples from throughout the Columbia River Basin, along with representative collections from adjacent basins, and genotyped them using a panel of 325 single-nucleotide markers. Results from individual- and group-based analyses of these data indicate that white sturgeon in the uppermost Columbia River Basin, in the Kootenai and upper Snake Rivers, are the most distinct, while the remaining populations downstream in the basin can be described as a genetic gradient consistent with an isolation-by-distance effect. Notably, the population in the lowest reaches of the Columbia River is more distinct from the middle or upper reaches than from outside basins, and suggests historically a higher or more recent gene exchange through coastal routes than with populations in the interior Columbia Basin. Nonetheless, proximal reaches were generally only marginally or non-significantly divergent, suggesting that transplanting larvae or juveniles from nearby sources poses relatively little risk of outbreeding depression. Indeed, we inferred examples of dispersal between reaches via close-kin mark-recapture and genetic mark-recapture that indicate movement between nearby reaches is not unusual. Samples from the Kootenai and upper Snake Rivers exhibited notably lower genetic diversity than the remaining samples as a result of population bottlenecks, genetic drift, and/or historical divergence. Conservation actions, such as supplementation, are underway to maintain population viability and will require balanced efforts to increase demographic abundance while maintaining genetic diversity.

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.

Weak population differentiation and high diversity in Salsola tragus in the inland Pacific Northwest, USA

BACKGROUND

Salsola tragus is a widespread and problematic weed of semi-arid wheat production globally, and in the inland Pacific Northwest region of the USA. The species exhibits high levels of phenotypic diversity across its range and, at least in California USA, previous work has described cryptic diversity comprising a multi-species complex. Such cryptic diversity could suggest the potential for a differential response to management inputs between groups, and have important implications for the spread of herbicide resistance or other adaptive traits within populations. We used a genotyping-by-sequencing approach to characterize the population structure of S. tragus in the inland Pacific Northwest.

RESULTS

Our results indicated that the population in this region is comprised of a single, tetraploid species (S. tragus sensu latu) with weak population structure on a regional scale. Isolation-by-distance appears to be the primary pattern of structure, but an independent set of weakly differentiated clusters of unknown origin were also apparent, along with a mixed mating system and high levels of largely unstructured genetic diversity.

CONCLUSIONS

Despite considerable phenotypic variability within S. tragus in the region, agronomic weed managers can likely consider it as a single entity across the region, rather than a collection of cryptic subgroups with possible differential responses to management inputs or agroecosystem conditions. A lack of strong barriers to migration and gene flow mean that adaptive traits, such as herbicide resistance, can be expected to spread rapidly through populations across the region. © 2022 Society of Chemical Industry.

Pollinator movement activity influences genetic diversity and differentiation of spatially isolated populations of clonal forest herbs

In agricultural landscapes, forest herbs live in small, spatially isolated forest patches. For their long-term survival, their populations depend on animals as genetic linkers that provide pollen- or seed-mediated gene flow among different forest patches. However, whether insect pollinators serve as genetic linkers among spatially isolated forest herb populations in agricultural landscapes remains to be shown. Here, we used population genetic methods to analyze: (A) the genetic diversity and genetic differentiation of populations of two common, slow-colonizing temperate forest herb species [Polygonatum multiflorum (L.) All. and Anemone nemorosa L.] in spatially isolated populations within three agricultural landscapes in Germany and Sweden and (B) the movement activity of their most relevant associated pollinator species, i.e., the bumblebee Bombus pascuorum (Scopoli, 1,763) and the hoverfly Melanostoma scalare (Fabricus, 1,794), respectively, which differ in their mobility. We tested whether the indicated pollinator movement activity affected the genetic diversity and genetic differentiation of the forest herb populations. Bumblebee movement indicators that solely indicated movement activity between the forest patches affected both genetic diversity and genetic differentiation of the associated forest herb P. multiflorum in a way that can be explained by pollen-mediated gene flow among the forest herb populations. In contrast, movement indicators reflecting the total movement activity at a forest patch (including within-forest patch movement activity) showed unexpected effects for both plant-pollinator pairs that might be explained by accelerated genetic drift due to enhanced sexual reproduction. Our integrated approach revealed that bumblebees serve as genetic linkers of associated forest herb populations, even if they are more than 2 km apart from each other. No such evidence was found for the forest associated hoverfly species which showed significant genetic differentiation among forest patches itself. Our approach also indicated that a higher within-forest patch movement activity of both pollinator species might enhance sexual recruitment and thus diminishes the temporal buffer that clonal growth provides against habitat fragmentation effects.

Phenotypic variability and genetic diversity analysis of cultivated potatoes in China

Phenotypic evaluation and molecular biotechnology are both important in the identification and utilization of crop germplasm resources. In this study, the phenotypic variation and genetic diversity of 149 main potato cultivars in China were investigated with 12 phenotypic traits and 24 SSR markers. The coefficient of variation of 12 phenotypic traits ranged from 12.11% to 156.93%. The results of SSR markers exhibited a relatively high level of genetic variation (Na =5.458 ± 1.499, Ne =3.300 ± 1.087, I =1.397 ± 0.298, Ho =0.797 ± 0.178, He = 0.660 ± 0.117, and PIC=0.702 ± 0.087). Population structure and phylogenetic tree analysis divided the varieties into three subgroups. The results indicated that ninety percent of the molecular variance was attributed to within-group differences, and the remaining 10% was attributed to variation among groups. Consistent with previous report, alleles of the STI032 marker were significantly associated with tuber starch content and growth period traits in the population. The results of this study could facilitate the utilization of potato germplasm resources, molecular genetic breeding and improvement.

Genotyping-by-Sequencing and Morphology Revealed the Role of Polyploidization and Hybridization in the Diversification of the Centaurea aspera L. Complex of Section Seridia (Juss.) DC. (Asteraceae)

Hybridization and polyploidy are major drivers of plant evolution. In Centaurea (Asteraceae), both mechanisms are frequent and lead to reticulate evolutions. However, in the Western Mediterranean section, Seridia studies are scarce. In this section, Centaurea aspera forms a complex including four European diploid and one Moroccan autotetraploid subspecies, an allopolyploid, and hybrids among them. Here, we aimed to delimit the different taxa, identify any introgressions, and discuss their evolutionary history. Samples of all taxa were analysed using 1688 SNPs obtained through GBS and were morphologically characterized. Three genetically well-differentiated clusters were observed, corresponding to the allopolyploid C. seridis, the diploid C. aspera and the cryptic autotetraploid C. aspera ssp. gentilii, which is proposed to be considered as a species. Centaurea seridis showed a high isolation by distance, a greater morphological variability, and a lack of interspecific gene flow. Diploid and autotetraploid C. aspera individuals were morphologically similar, and some introgressions were detected in Southern Spain, where new forms may promote diversification. This gene flow might have taken place during the Messinian and before autopolyploidization occurred in Morocco. In the C. aspera complex, current interspecific barriers are strong, while polyploidization may provide a better adaptation to drier environments.

Asymptotic tests for Hardy-Weinberg equilibrium in hexaploids

Hexaploids, a group of organisms containing three complete sets of chromosomes in a single nucleus, are of utmost importance to evolutionary studies and breeding programs. Many studies have focused on hexaploid linkage analysis and QTL mapping in controlled crosses, but little methodology has been developed to reveal how hexaploids diversify and evolve in natural populations. We formulate a general framework for studying the pattern of genetic variation in autohexaploid populations through testing deviation from Hardy-Weinberg equilibrium (HWE) at individual molecular markers. We confirm that hexaploids cannot reach exact HWE but can approach asymptotic HWE at 8-9 generations of random mating. We derive a statistical algorithm for testing HWE and the occurrence of double reduction for autopolyploids, a phenomenon that affects population variation during long evolutionary processes. We perform computer simulation to validate the statistical behavior of our test procedure and demonstrate its usefulness by analyzing a real data set for autohexaploid chrysanthemum. When extended to allohexaploids, our test procedure will provide a generic tool for illustrating the genome structure of hexaploids in the quest to infer their evolutionary status and design association studies of complex traits.

Genetic Diversity and Population Structure of Psidium Species from Restinga: A Coastal and Disturbed Ecosystem of the Brazilian Atlantic Forest

The Atlantic Forest is one of the most threatened biomes in the world. Here, we use a common set of microsatellite markers to assess the genetic diversity and population structure of three species from the genus Psidium (P. guajava, P. macahense, and P. guineense), located in a disturbed environment of the Atlantic Forest, the restinga, in Espírito Santo, Brazil. Psidium guajava populations presented the highest number of alleles (95) followed by P. guineense (81) and P. macahense (68). The genetic variability was high (P. guajava = 0.71; P. guineense = 0.74; P. macahense = 0.63), with greater variation within populations (72 to 84%) than among populations (15 to 27%), reflecting elevated values of genetic differentiation (P. guajava, F_ST: 0.15; P. macahense, Ø_ST: 0.27; P. guineense, Ø_ST: 0.21). The populations were clustered into two main groups and considered moderately structured. This is the first report of genetic studies and evidence of polyploidy to P. macahense. Our results may provide information that can be used in management and conservation strategies, to preserve the diversity of Psidium populations.

The transitivity of the Hardy-Weinberg law

The Hardy-Weinberg law is shown to be transitive in the sense that a multi-allelic polymorphism that is in equilibrium will retain its equilibrium status if any allele together with its corresponding genotypes is deleted from the population. Similarly, the transitivity principle also applies if alleles are joined, which leads to the summation of allele frequencies and their corresponding genotype frequencies. These basic polymorphism properties are intuitive, but they have apparently not been formalized or investigated. This article provides a straightforward proof of the transitivity principle, and its usefulness in genetic data analysis is explored, using high-quality autosomal microsatellite databases from the US National Institute of Standards and Technology. We address the reduction of multi-allelic polymorphisms to variants with fewer alleles, two in the limit. Equilibrium test results obtained with the original and reduced polymorphisms are generally observed to be coherent, in particular when results obtained with length-based and sequence-based microsatellites are compared. We exploit the transitivity principle in order to identify disequilibrium-related alleles, and show its usefulness for detecting population substructure and genotyping problems that relate to null alleles and allele imbalance.

Characterizing Tetraploid Populations of Actinidia chinensis for Kiwifruit Genetic Improvement

Understanding genetic diversity and structure in natural populations and their suitable habitat response to environmental changes is critical for the protection and utilization of germplasm resources. We evaluated the genetic diversity and structure of 24 A. chinensis populations using simple sequence repeat (SSR) molecular markers. The potential suitable distribution of tetraploid A. chinensis estimated under the current climate and predicted for the future climate was generated with ecological niche modeling (ENM). The results indicated that the polyploid populations of A.chinensis have high levels of genetic diversity and that there are distinct eastern and western genetic clusters. The population structure of A. chinensis can be explained by an isolation-by-distance model. The results also revealed that potentially suitable areas of tetraploids will likely be gradually lost and the habitat will likely be increasingly fragmented in the future. This study provides an extensive overview of tetraploid A. chinensis across its distribution range, contributing to a better understanding of its germplasm resources. These results can also provide the scientific basis for the protection and sustainable utilization of kiwifruit wild resources.

A generalist-specialist trade-off between switchgrass cytotypes impacts climate adaptation and geographic range

Polyploidy, which occurs in roughly half of all flowering plants and an even higher percentage of grasses, is thought to be a major driver of adaptation. Higher numbers of copies of each gene in polyploid genomes can increase genetic diversity, which could drive shifts in habitat preference, adaptability, and fitness. To test the effects of increased ploidy, we compared genomic diversity, environmental niche, and fitness responses across climatic gradients between tetraploid and octoploid switchgrass. We found that the octoploids contained novel combinations of the ancestral tetraploid genetic diversity, which was linked to the expansion of switchgrass into unsuitable habitats for tetraploid populations. Our experiments revealed evidence of niche divergence, differential fitness, and a generalist–specialist trade-off between cytotypes.

Polyploidy results from whole-genome duplication and is a unique form of heritable variation with pronounced evolutionary implications. Different ploidy levels, or cytotypes, can exist within a single species, and such systems provide an opportunity to assess how ploidy variation alters phenotypic novelty, adaptability, and fitness, which can, in turn, drive the development of unique ecological niches that promote the coexistence of multiple cytotypes. Switchgrass, Panicum virgatum, is a widespread, perennial C4 grass in North America with multiple naturally occurring cytotypes, primarily tetraploids (4×) and octoploids (8×). Using a combination of genomic, quantitative genetic, landscape, and niche modeling approaches, we detect divergent levels of genetic admixture, evidence of niche differentiation, and differential environmental sensitivity between switchgrass cytotypes. Taken together, these findings support a generalist (8×)–specialist (4×) trade-off. Our results indicate that the 8× represent a unique combination of genetic variation that has allowed the expansion of switchgrass’ ecological niche and thus putatively represents a valuable breeding resource.

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

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

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.

Shorebirds' Longer Migratory Distances Are Associated With Larger ADCYAP1 Microsatellites and Greater Morphological Complexity of Hippocampal Astrocytes

For the epic journey of autumn migration, long-distance migratory birds use innate and learned information and follow strict schedules imposed by genetic and epigenetic mechanisms, the details of which remain largely unknown. In addition, bird migration requires integrated action of different multisensory systems for learning and memory, and the hippocampus appears to be the integration center for this task. In previous studies we found that contrasting long-distance migratory flights differentially affected the morphological complexity of two types of hippocampus astrocytes. Recently, a significant association was found between the latitude of the reproductive site and the size of the ADCYAP1 allele in long distance migratory birds. We tested for correlations between astrocyte morphological complexity, migratory distances, and size of the ADCYAP1 allele in three long-distance migrant species of shorebird and one non-migrant. Significant differences among species were found in the number and morphological complexity of the astrocytes, as well as in the size of the microsatellites of the ADCYAP1 gene. We found significant associations between the size of the ADCYAP1 microsatellites, the migratory distances, and the degree of morphological complexity of the astrocytes. We suggest that associations between astrocyte number and morphological complexity, ADCYAP1 microsatellite size, and migratory behavior may be part of the adaptive response to the migratory process of shorebirds.

A Tracing Model for the Evolutionary Equilibrium of Octoploids

Testing Hardy-Weinberg equilibrium (HWE) is a fundamental approach for inferring population diversity and evolution, but its application to octoploids containing eight chromosome sets has not well been justified. We derive a mathematical model to trace how genotype frequencies transmit from parental to offspring generations in the natural populations of autooctoploids. We find that octoploids, including autooctolpoids undergoing double reduction, attach asymptotic HWE (aHWE) after 15 generations of random mating, in a contrast to diploids where one generation can assure exact equilibrium and, also, different from tetraploids that use 5 generations to reach aHWE. We develop a statistical procedure for testing aHWE in octoploids and apply it to analyze a real data set from octoploid switchgrass distributed in two ecologically different regions, demonstrating the usefulness of the test procedure. Our model provides a tool for studying the population genetic diversity of octoploids, inferring their evolutionary history, and identifying the ecological relationship of octoploid-genome structure with environmental adaptation.

The activation of gene expression and alternative splicing in the formation and evolution of allopolyploid Brassica napus

Allopolyploids contain two or more sets of subgenomes. To establish a compatible relationship between subgenomes, a series of gene expression changes occurred in allopolyploids. What evolutionary changes of transcripts have taken place in Brassica napus during the early establishment and subsequent evolution was a fascinating scientific question. Here, we study this issue using a set of materials (natural, resynthesized B. napus and their progenitors/parents) by long-read RNA sequencing technology. The results showed that more genes were up-regulated in resynthesized B. napus compared with its two parents, and more up-regulated expressed genes were observed in natural B. napus compared with resynthesized B. napus. The presence of up-regulation genes in organism may help it adapt to the influence of "genomic shock" and cope with natural environment. Isoforms are produced from precursor mRNAs by alternative splicing (AS) events, and more than 60% of novel isoforms were identified in all materials, which could improve the reference genome information of B. napus. We found that the isoform numbers, the number of genes potentially involved in AS and alternative polyadenylation increased in B. napus after evolution, which may involve in the adaptation of plants to natural environment. In addition, all identified isoforms were functional annotated by searching 7 databases. In general, this study could improve our overall understanding of the full-length transcriptome of B. napus, and help us recognize the significant gene expression changes and isoform abundance changes occurred in allopolyploid B. napus during evolution.

When do autopolyploids need poly-sequencing data?

The sequencing depth required to genotype autopolyploid populations is a very controversial topic. Different studies have adopted variable depth values without a clear guide on the optimal sequencing depth value. Many studies suggest high depth thresholds for different ploidies that may not be practical and substantially increase the overall genotyping cost for different projects. However, such conservative thresholds may not be required to achieve the most common research goals. In fact, some recent reports in the field of quantitative genetics found that much lower sequencing depth thresholds could achieve the same accuracy as high depth thresholds. In this manuscript, I discuss when researchers need to use stringent sequencing depth thresholds and when they can use more relaxed ones. I support my argument by calculating the probabilities of sampling different homologues at a given sequencing depth. I also discuss the uses and the uncertainty in calculating a continuous allelic dosage as the proportion of sequencing reads that hold the alternative allele, which is becoming a common method now in quantitative genetics to replace discrete dosage estimation.

Parallel introgression, not recurrent emergence, explains apparent elevational ecotypes of polyploid Himalayan snowtrout

The recurrence of similar evolutionary patterns within different habitats often reflects parallel selective pressures acting upon either standing or independently occurring genetic variation to produce a convergence of phenotypes. This interpretation (i.e. parallel divergences within adjacent streams) has been hypothesized for drainage-specific morphological 'ecotypes' observed in polyploid snowtrout (Cyprinidae: Schizothorax). However, parallel patterns of differential introgression during secondary contact are a viable alternative hypothesis. Here, we used ddRADseq (N = 35 319 de novo and N = 10 884 transcriptome-aligned SNPs), as derived from Nepali/Bhutanese samples (N = 48 each), to test these competing hypotheses. We first employed genome-wide allelic depths to derive appropriate ploidy models, then a Bayesian approach to yield genotypes statistically consistent under the inferred expectations. Elevational 'ecotypes' were consistent in geometric morphometric space, but with phylogenetic relationships at the drainage level, sustaining a hypothesis of independent emergence. However, partitioned analyses of phylogeny and admixture identified subsets of loci under selection that retained genealogical concordance with morphology, suggesting instead that apparent patterns of morphological/phylogenetic discordance are driven by widespread genomic homogenization. Here, admixture occurring in secondary contact effectively 'masks' previous isolation. Our results underscore two salient factors: (i) morphological adaptations are retained despite hybridization and (ii) the degree of admixture varies across tributaries, presumably concomitant with underlying environmental or anthropogenic factors.

Comparative Transcriptome Analysis Revealed Genes Involved in Sexual and Polyploid Growth Dimorphisms in Loach (Misgurnus anguillicaudatus)

Simple Summary

Misgurnus anguillicaudatus not only exhibits sexual size dimorphism, but also shows polyploid size dimorphism. Here, we performed comparative transcriptome integration analysis of multiple tissues of diploid and tetraploid M. anguillicaudatus of both sexes. We found that differences in energy metabolism and steroid hormone synthesis levels may be the main causes of sexual and polyploidy growth dimorphisms of M. anguillicaudatus. Fast-growing M. anguillicaudatus (tetraploids, females) have higher levels of energy metabolism and lower steroid hormone synthesis and fatty acid degradation abilities than slow-growing M. anguillicaudatus (diploids, males).

Abstract

Sexual and polyploidy size dimorphisms are widespread phenomena in fish, but the molecular mechanisms remain unclear. Loach (Misgurnus anguillicaudatus) displays both sexual and polyploid growth dimorphism phenomena, and are therefore ideal models to study these two phenomena. In this study, RNA-seq was used for the first time to explore the differentially expressed genes (DEGs) between both sexes of diploid and tetraploid loaches in four tissues (brain, gonad, liver, and muscle). Results showed that 21,003, 17, and 1 DEGs were identified in gonad, liver, and muscle tissues, respectively, between females and males in both diploids and tetraploids. Regarding the ploidy levels, 4956, 1496, 2187, and 1726 DEGs were identified in the brain, gonad, liver, and muscle tissues, respectively, between tetraploids and diploids of the same sex. When both sexual and polyploid size dimorphisms were considered simultaneously in the four tissues, only 424 DEGs were found in the gonads, indicating that these gonadal DEGs may play an important regulatory role in regulating sexual and polyploid size dimorphisms. Regardless of the sex or ploidy comparison, the significant DEGs involved in glycolysis/gluconeogenesis and oxidative phosphorylation pathways were upregulated in faster-growing individuals, while steroid hormone biosynthesis-related genes and fatty acid degradation and elongation-related genes were downregulated. This suggests that fast-growing loaches (tetraploids, females) have higher energy metabolism levels and lower steroid hormone synthesis and fatty acid degradation abilities than slow-growing loaches (diploids, males). Our findings provide an archive for future systematic research on fish sexual and polyploid dimorphisms.

Alignment-free methods for polyploid genomes: Quick and reliable genetic distance estimation

Polyploid genomes pose several inherent challenges to population genetic analyses. While alignment-based methods are fundamentally limited in their applicability to polyploids, alignment-free methods bypass most of these limits. We investigated the use of Mash, a k-mer analysis tool that uses the MinHash method to reduce complexity in large genomic data sets, for basic population genetic analyses of polyploid sequences. We measured the degree to which Mash correctly estimated pairwise genetic distance in simulated haploid and polyploid short-read sequences with various levels of missing data. Mash-based estimates of genetic distance were comparable to alignment-based estimates, and were less impacted by missing data. We also used Mash to analyse publicly available short-read data for three polyploid and one diploid species, then compared Mash results to published results. For both simulated and real data, Mash accurately estimated pairwise genetic differences for polyploids as well as diploids as much as 476 times faster than alignment-based methods, though we found that Mash genetic distance estimates could be biased by per-sample read depth. Mash may be a particularly useful addition to the toolkit of polyploid geneticists for rapid confirmation of alignment-based results and for basic population genetics in reference-free systems or those with only poor-quality sequence data available.

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

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

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

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

MultiGWAS: An integrative tool for Genome Wide Association Studies in tetraploid organisms

The genome-wide association studies (GWASs) are essential to determine the genetic bases of either ecological or economic phenotypic variation across individuals within populations of the model and nonmodel organisms. For this research question, the GWAS replication testing different parameters and models to validate the results' reproducibility is common. However, straightforward methodologies that manage both replication and tetraploid data are still missing. To solve this problem, we designed the MultiGWAS, a tool that does GWAS for diploid and tetraploid organisms by executing in parallel four software packages, two designed for polyploid data (GWASpoly and SHEsis) and two designed for diploid data (GAPIT and TASSEL). MultiGWAS has several advantages. It runs either in the command line or in a graphical interface; it manages different genotype formats, including VCF. Moreover, it allows control for population structure, relatedness, and several quality control checks on genotype data. Besides, MultiGWAS can test for additive and dominant gene action models, and, through a proprietary scoring function, select the best model to report its associations. Finally, it generates several reports that facilitate identifying false associations from both the significant and the best-ranked association Single Nucleotide Polymorphisms (SNPs) among the four software packages. We tested MultiGWAS with public tetraploid potato data for tuber shape and several simulated data under both additive and dominant models. These tests demonstrated that MultiGWAS is better at detecting reliable associations than using each of the four software packages individually. Moreover, the parallel analysis of polyploid and diploid software that only offers MultiGWAS demonstrates its utility in understanding the best genetic model behind the SNP association in tetraploid organisms. Therefore, MultiGWAS probed to be an excellent alternative for wrapping GWAS replication in diploid and tetraploid organisms in a single analysis environment.

Patterns of population genomic diversity in the invasive Japanese knotweed species complex

PREMISE

Invasive species are expected to undergo a reduction in genetic diversity due to founder effects, which should limit their ability to adapt to new habitats. Still, many invasive species achieve widespread distributions and dense populations. This paradox of invasions could potentially be overcome through multiple introductions or hybridization, both of which increase genetic diversity. We conducted a population genomics study of Japanese knotweed (Reynoutria japonica), which is a polyploid, clonally reproducing invasive species that has been notoriously successful worldwide despite supposedly low genetic diversity.

METHODS

We used genotyping by sequencing to collect 12,912 SNP markers from 88 samples collected at 38 locations across North America for the species complex. We used alignment-free k-mer hashing analysis in addition to traditional population genetic analyses to account for the challenges of genotyping polyploids.

RESULTS

Genotypes conformed to three genetic clusters, likely representing Japanese knotweed, giant knotweed, and hybrid bohemian knotweed. We found that, contrary to previous findings, the Japanese knotweed cluster had substantial genetic diversity, though it had no apparent genetic structure across the landscape. In contrast, giant knotweed and hybrids showed distinct population groups. We did not find evidence of isolation by distance in the species complex, likely reflecting the stochastic introduction history of this species complex.

CONCLUSIONS

The results indicate that clonal invasive species can show substantial genetic diversity and can be successful at colonizing a variety of habitats without showing evidence of local adaptation or genetic structure.

Low dispersal and ploidy differences in a grass maintain photosynthetic diversity despite gene flow and habitat overlap

Geographical isolation facilitates the emergence of distinct phenotypes within a single species, but reproductive barriers or selection are needed to maintain the polymorphism after secondary contact. Here, we explore the processes that maintain intraspecific variation of C₄ photosynthesis, a complex trait that results from the combined action of multiple genes. The grass Alloteropsis semialata includes C₄ and non-C₄ populations, which have coexisted as a polyploid series for more than 1 million years in the miombo woodlands of Africa. Using population genomics, we show that there is genome-wide divergence for the photosynthetic types, but the current geographical distribution does not reflect a simple habitat displacement scenario as the genetic clusters overlap, being occasionally mixed within a given habitat. Despite evidence of recurrent introgression between non-C₄ and C₄ groups, in both diploids and polyploids, the distinct genetic lineages retain their identity, potentially because of selection against hybrids. Coupled with strong isolation by distance within each genetic group, this selection created a geographical mosaic of photosynthetic types. Diploid C₄ and non-C₄ types never grew together, and the C₄ type from mixed populations constantly belonged to the hexaploid lineage. By limiting reproductive interactions between photosynthetic types, the ploidy difference probably allows their co-occurrence, reinforcing the functional diversity within this species. Together, these factors enabled the persistence of divergent physiological traits of ecological importance within a single species despite gene flow and habitat overlap.

Computational characterization of double reduction in autotetraploid natural populations

Population genetic theory has been well developed for diploid species, but its extension to study genetic diversity, variation and evolution in autopolyploids, a class of polyploids derived from the genome doubling of a single ancestral species, requires the incorporation of multisomic inheritance. Double reduction, which is characteristic of autopolyploidy, has long been believed to shape the evolutionary consequence of organisms in changing environments. Here, we develop a computational model for testing and estimating double reduction and its genomic distribution in autotetraploids. The model is implemented with the expectation-maximization (EM) algorithm to dissect unobservable allelic recombinations among multiple chromosomes, enabling the simultaneous estimation of allele frequencies and double reduction in natural populations. The framework fills an important gap in the population genetic theory of autopolyploids.

Model-based genotype and ancestry estimation for potential hybrids with mixed-ploidy

Non-random mating among individuals can lead to spatial clustering of genetically similar individuals and population stratification. This deviation from panmixia is commonly observed in natural populations. Consequently, individuals can have parentage in single populations or involving hybridization between differentiated populations. Accounting for this mixture and structure is important when mapping the genetics of traits and learning about the formative evolutionary processes that shape genetic variation among individuals and populations. Stratified genetic relatedness among individuals is commonly quantified using estimates of ancestry that are derived from a statistical model. Development of these models for polyploid and mixed-ploidy individuals and populations has lagged behind those for diploids. Here, we extend and test a hierarchical Bayesian model, called entropy, which can use low-depth sequence data to estimate genotype and ancestry parameters in autopolyploid and mixed-ploidy individuals (including sex chromosomes and autosomes within individuals). Our analysis of simulated data illustrated the trade-off between sequencing depth and genome coverage and found lower error associated with low-depth sequencing across a larger fraction of the genome than with high-depth sequencing across a smaller fraction of the genome. The model has high accuracy and sensitivity as verified with simulated data and through analysis of admixture among populations of diploid and tetraploid Arabidopsis arenosa.

Genetic Diversity and Connectivity in Plant Species Differing in Clonality and Dispersal Mechanisms in Wetland Island Habitats

In plants, long-distance dispersal is both attenuated and directed by specific movement vectors, including animals, wind, and/or water. Hence, movement vectors partly shape metapopulation genetic patterns that are, however, also influenced by other life-history traits such as clonal growth. We studied the relationship between area, isolation, plant-species richness, reproduction, and dispersal mechanisms with genetic diversity and divergence in 4 widespread wetland plant-species in a total of 20 island-like kettle-hole habitats surrounded by an intensive agricultural landscape. Our results showed that genetic parameters reflect the reproduction strategies with the highest genetic diversity being observed in the non-clonal, outcrossing Oenanthe aquatica compared to the clonal Lycopus europaeus, Typha latifolia, and Phragmites australis. Lycopus showed a positive relationship between genetic diversity and kettle-hole area, but a negative relationship with the number of neighboring kettle holes (less isolation). Genetic diversity increased with plant-species richness in the clonal species Phragmites and Lycopus; while it decreased in the non-clonal Oenanthe. Finally, genetic divergence and, therefore, connectivity differed between alternative dispersal strategies, where wind-dispersed Typha and Phragmites had a higher gene flow between the analyzed kettle holes compared with the insect-pollinated, hydrochorous Lycopus and Oenanthe. Our study provides information on genetic patterns related to reproduction and dispersal mechanisms of 4 common wetland species contributing to the understanding of the functioning of plant metacommunities occurring in kettle holes embedded in agricultural landscapes.

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

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

Estimation of Molecular Pairwise Relatedness in Autopolyploid Crops

A suitable pairwise relatedness estimation is key to genetic studies. Several methods are proposed to compute relatedness in autopolyploids based on molecular data. However, unlike diploids, autopolyploids still need further studies considering scenarios with many linked molecular markers with known dosage. In this study, we provide guidelines for plant geneticists and breeders to access trustworthy pairwise relatedness estimates. To this end, we simulated populations considering different ploidy levels, meiotic pairings patterns, number of loci and alleles, and inbreeding levels. Analysis were performed to access the accuracy of distinct methods and to demonstrate the usefulness of molecular marker in practical situations. Overall, our results suggest that at least 100 effective biallelic molecular markers are required to have good pairwise relatedness estimation if methods based on correlation is used. For this number of loci, current methods based on multiallelic markers show lower performance than biallelic ones. To estimate relatedness in cases of inbreeding or close relationships (as parent-offspring, full-sibs, or half-sibs) is more challenging. Methods to estimate pairwise relatedness based on molecular markers, for different ploidy levels or pedigrees were implemented in the AGHmatrix R package.

Comparative analysis of genomic- and EST-SSRs in European plum (Prunus domestica L.): implications for the diversity analysis of polyploids

Simple sequence repeats (SSRs) are among the most useful DNA markers in plant science. The aim of this study was to compare the features and usefulness of genomic SSR (gSSR) and EST-SSRs in European plum (Prunus domestica L.), an economically important, hexaploid stone fruit crop globally cultivated to produce fleshy fruits and derived foodstuff. The analysis of an ample set of morphologically diverse varieties indicated that gSSRs and EST-SSRs provide different estimates of some of the locus-based indicators of diversity. Moreover, the two classes of SSRs gave different, weakly correlated, estimations of distance-based parameters with gSSRs being more powerful for discriminating purposes. The two SSR classes provide complementary information in European plum, making the contribution of EST-SSRs useful not only as non-neutral markers. The differences between SSR classes are discussed considering the neutral and non-neutral evolution, and the polyploidy and asexual propagation of the cultivated tree varieties.

Patterns of Genetic Diversification in the Invasive Hybrid Plant Pathogen Phytophthora × alni and Its Parental Species P. uniformis

In pathogenic fungi and oomycetes, interspecific hybridization may lead to the formation of new species having a greater impact on natural ecosystems than the parental species. From the early 1990s, a severe alder (Alnus spp.) decline due to an unknown Phytophthora species was observed in several European countries. Genetic analyses revealed that the disease was caused by the triploid hybrid P. × alni, which originated in Europe from the hybridization of P. uniformis and P. × multiformis. Here, we investigated the population structure of P. × alni (158 isolates) and P. uniformis (85 isolates) in several European countries using microsatellite markers. Our analyses confirmed the genetic structure previously observed in other European populations, with P. uniformis populations consisting of at most two multilocus genotypes (MLGs) and P. × alni populations dominated by MLG Pxa-1. The genetic structure of P. × alni populations in the Czech Republic, Hungary and Sweden seemed to reflect the physical isolation of river systems. Most rare P. × alni MLGs showed a loss of heterozygosity (LOH) at one or a few microsatellite loci compared with other MLGs. This LOH may allow a stabilization within the P. × alni genome or a rapid adaptation to stress situations. Alternatively, alleles may be lost because of random genetic drift in small, isolated populations, with no effect on fitness of P. × alni. Additional studies would be necessary to confirm these patterns of population diversification and to better understand the factors driving it.