Comparative Transcriptomics of Strawberries (Fragaria spp.) Provides Insights into Evolutionary Patterns

The relationship between niche breadth and phylogenetic characteristics of eight species of rhubarb on the Qinghai-Tibet Plateau, Asia

The relationship between spatial distribution and phylogeny has been widely debated in recent decades. To understand biogeographic and evolutionary history relationships and to explore the interspecific similarities and phylogenetic correlations of niche characteristics, we collected and recorded all distribution points for eight species of rhubarb on the Qinghai-Tibet Plateau, used different methods to describe the ecological niche, and explored the relationship between phylogeny, ecological niche, and distribution range. The results reveal that: (1) the ranges of optimal environmental variables for species with close kinship are not exactly the same, ecologically similar species are not necessarily sister species, and the overlap of rhubarb has no significant correlation with phylogeny. Therefore, the impact of ecological dimensions on species formation is greater than that of geographical latitude for the eight species of rhubarb. (2) Among the eight species of rhubarb, the breadth of ecological niche is positively correlated with the current suitable habitat area and negatively correlated with fluctuations in future suitable habitat area. In the future, except for Rheum tanguticum and Rheum palmatum, the suitable planting areas for the other six species of rhubarb will decrease as greenhouse gas emissions concentrations and time increase. Therefore, species with smaller ecological niches are at a greater risk of habitat loss compared to species with larger ecological niches. (3) In both existing and future distribution prediction models of rhubarb, we observed that both the widely distributed Rheum spiciforme and the niche narrow Rheum nobile, all eight species of rhubarb are present in the Hengduan Mountains, based on our analysis, we propose that the Hengduan Mountains should be regarded as a priority conservation area for rhubarb, to preserve the species' biodiversity. Our study lays the groundwork for identifying evolutionary trends in ecological specialization.

Mosaic environment-driven evolution of the deep-sea mussel Gigantidas platifrons bacterial endosymbiont

BACKGROUND

The within-species diversity of symbiotic bacteria represents an important genetic resource for their environmental adaptation, especially for horizontally transmitted endosymbionts. Although strain-level intraspecies variation has recently been detected in many deep-sea endosymbionts, their ecological role in environmental adaptation, their genome evolution pattern under heterogeneous geochemical environments, and the underlying molecular forces remain unclear.

RESULTS

Here, we conducted a fine-scale metagenomic analysis of the deep-sea mussel Gigantidas platifrons bacterial endosymbiont collected from distinct habitats: hydrothermal vent and methane seep. Endosymbiont genomes were assembled using a pipeline that distinguishes within-species variation and revealed highly heterogeneous compositions in mussels from different habitats. Phylogenetic analysis separated the assemblies into three distinct environment-linked clades. Their functional differentiation follows a mosaic evolutionary pattern. Core genes, essential for central metabolic function and symbiosis, were conserved across all clades. Clade-specific genes associated with heavy metal resistance, pH homeostasis, and nitrate utilization exhibited signals of accelerated evolution. Notably, transposable elements and plasmids contributed to the genetic reshuffling of the symbiont genomes and likely accelerated adaptive evolution through pseudogenization and the introduction of new genes.

CONCLUSIONS

The current study uncovers the environment-driven evolution of deep-sea symbionts mediated by mobile genetic elements. Its findings highlight a potentially common and critical role of within-species diversity in animal-microbiome symbioses. Video Abstract.

Transposon signatures of allopolyploid genome evolution

Hybridization brings together chromosome sets from two or more distinct progenitor species. Genome duplication associated with hybridization, or allopolyploidy, allows these chromosome sets to persist as distinct subgenomes during subsequent meioses. Here, we present a general method for identifying the subgenomes of a polyploid based on shared ancestry as revealed by the genomic distribution of repetitive elements that were active in the progenitors. This subgenome-enriched transposable element signal is intrinsic to the polyploid, allowing broader applicability than other approaches that depend on the availability of sequenced diploid relatives. We develop the statistical basis of the method, demonstrate its applicability in the well-studied cases of tobacco, cotton, and Brassica napus, and apply it to several cases: allotetraploid cyprinids, allohexaploid false flax, and allooctoploid strawberry. These analyses provide insight into the origins of these polyploids, revise the subgenome identities of strawberry, and provide perspective on subgenome dominance in higher polyploids.

Stem transcriptome screen for selection in wild and cultivated pitahaya (Selenicereus undatus): an epiphytic cactus with edible fruit

Dragon fruit, pitahaya or pitaya are common names for the species in the Hylocereus group of Selenicereus that produce edible fruit. These Neotropical epiphytic cacti are considered promising underutilized crops and are currently cultivated around the world. The most important species, S. undatus, has been managed in the Maya domain for centuries and is the focus of this article. Transcriptome profiles from stems of wild and cultivated plants of this species were compared. We hypothesized that differences in transcriptomic signatures could be associated with genes related to drought stress. De novo transcriptome assembly and the analysis of differentially expressed genes (DEGs) allowed us to identify a total of 9,203 DEGs in the Hunucmá cultivar relative of wild Mozomboa plants. Of these, 4,883 represent up-regulated genes and 4,320, down-regulated genes. Additionally, 6,568 DEGs were identified from a comparison between the Umán cultivar and wild plants, revealing 3,286 up-regulated and 3,282 down-regulated genes. Approximately half of the DEGs are shared by the two cultivated plants. Differences between the two cultivars that were collected in the same region could be the result of differences in management. Metabolism was the most representative functional category in both cultivars. The up-regulated genes of both cultivars formed a network related to the hormone-mediated signaling pathway that includes cellular responses to auxin stimulus and to hormone stimulus. These cellular reactions have been documented in several cultivated plants in which drought-tolerant cultivars modify auxin transport and ethylene signaling, resulting in a better redistribution of assimilates.

Divergence in interspecific and intersubspecific gene expression between two closely related horseshoe bats (Rhinolophus)

Closely related species have been used as representative systems to investigate the genetic mechanisms involved in the early stages of species differentiation. Previous studies have indicated that variation in gene expression might be a sensitive indicator of initial species divergence, although the role of expression divergence, and especially that associated with phenotypic variation remained relatively undefined. For three organs (cochlea, brain, and liver) from two closely related bat species (Rhinolophus siamensis and R. episcopus), the interspecific and intersubspecific gene expression profiles were compared using transcriptomics in this study. Striking organ specificity of expression was observed, and expression profiles exhibited similarities between cochlea and brain tissues. Numerous differentially expressed genes (DEGs) were identified for each organ in the interspecific comparison (cochlea/brain/liver: 1,069/647/692) and intersubspecific comparison (608/528/368). Functional enrichment analysis indicated vital variation in expression related to the immune system, ion activities, neuronal function, and multisensory system regulation in both comparisons. DEGs relevant to the variation in echolocation calls (RF) were found, and some of them were involved in the pivotal patterns of expression variation. The regulation of immune, ion channel, neural activity, and sophisticated sensory functions at the expression level might be key mechanisms in the early species divergence of bats, and the expression variation related to acoustical signal could have played a crucial part. This study expands our knowledge of gene expression and patterns of variation for three key organs to echolocation at both the interspecific and intersubspecific levels. Further, the framework described here provides insight into the genetic basis of phenotypic variation during the incipient stage of species differentiation.

Evolutionary history and pan-genome dynamics of strawberry (Fragaria spp.)

Strawberry (Fragaria spp.) has emerged as a model system for various fundamental and applied research in recent years. In total, the genomes of five different species have been sequenced over the past 10 y. Here, we report chromosome-scale reference genomes for five strawberry species, including three newly sequenced species' genomes, and genome resequencing data for 128 additional accessions to estimate the genetic diversity, structure, and demographic history of key Fragaria species. Our analyses obtained fully resolved and strongly supported phylogenies and divergence times for most diploid strawberry species. These analyses also uncovered a new diploid species (Fragaria emeiensis Jia J. Lei). Finally, we constructed a pan-genome for Fragaria and examined the evolutionary dynamics of gene families. Notably, we identified multiple independent single base mutations of the MYB10 gene associated with white pigmented fruit shared by different strawberry species. These reference genomes and datasets, combined with our phylogenetic estimates, should serve as a powerful comparative genomic platform and resource for future studies in strawberry.

Genetic Variability and Structure of Fragaria nilgerrensis Schlecht. Germplasm in Sichuan Province

Fragaria nilgerrensis Schlecht. (wild strawberry) is widely distributed in Southwest China, characterized by stress tolerance and the fruits of a notable peach aroma. So far there is only limited knowledge of variability and genetic structure in this species. Using AFLP markers, we investigated the genetic variability of 37 plants of F. nilgerrensis sampled in six main mountain areas of Sichuan Province and analyzed their genetic structure. Genetic similarity according to Nei and Li was used for cluster analysis based on UPGMA method and Agglomerative Hierarchical Clustering. Stratification of plants into more distinctive genetic groups was determined using Bayesian structure analysis. Six primer combinations produced a total of 1302 fragments of which 818 (62.8%) were polymorphic. Bayesian analysis showed the 37 plants of F. nilgerrensis grouped into five distinctive genetic groups. Most of the plants from the same mountain area clustered into the same genetic group, indicating each area as an area with the unique genetic profile. The genetic parameters analyzed here indicate a huge genetic variability of F. nilgerrensis in Sichuan Province. Our results provide reference data for surveying and protecting the germplasm resources of F. nilgerrensis that could be used in strawberry breeding programs.

Evaluation of genetic diversity and population structure of Fragaria nilgerrensis using EST-SSR markers

Fragaria nilgerrensis is a diploid wild strawberry widely distributed in Southwest China. Its white color and "peach-like" fragrance of fruits are valuable characters for the genetic improvement of cultivated strawberry plants. Its strong biotic and abiotic resistance and tolerance also enable it to survive in different habitats in the field. In this study, we evaluated the level of genetic variation within and between 16 populations with 169 individuals of F. nilgerrensis using 16 newly developed EST-SSR (expressed sequence tag-simple sequence repeats) markers. The results show that the genetic diversity of this species was high, based on Nei's genetic diversity (0.26) and polymorphic loci (0.41), although it is self-compatible and has clonal propagation. Significant genetic differentiation among populations was also detected by AMOVA analysis (Fst = 0.34), which could be indicative of little gene flow (Nm = 0.43) in F. nilgerrensis. The phylogenetic tree indicates that most of individuals from the same population have clustered together. These populations were not grouped based on the geographical distance, consistent with the Mantel test result (R² = 0.0063, P > 0.05). All the populations were assigned into two ancestral groups, with some individuals admixed, suggesting ancestral gene flow had occurred between these two groups. Our developed EST-SSR markers as well as the genetic diversity and population structure analysis of F. nilgerrensis are important for genetic improvement in the breeding process. Moreover, the populations that contain high genetic diversity would be a priority for collection and conservation.

Tracing the Diploid Ancestry of the Cultivated Octoploid Strawberry

The commercial strawberry, Fragaria × ananassa, is a recent allo-octoploid that is cultivated worldwide. However, other than Fragaria vesca, which is universally accepted one of its diploid ancestors, its other early diploid progenitors remain unclear. Here, we performed comparative analyses of the genomes of five diploid strawberries, F. iinumae, F. vesca, F. nilgerrensis, F. nubicola, and F. viridis, of which the latter three are newly sequenced. We found that the genomes of these species share highly conserved gene content and gene order. Using an alignment-based approach, we show that F. iinumae and F. vesca are the diploid progenitors to the octoploid F. × ananassa, whereas the other three diploids that we analyzed in this study are not parental species. We generated a fully resolved, dated phylogeny of Fragaria, and determined that the genus arose ∼6.37 Ma. Our results effectively resolve conflicting hypotheses regarding the putative diploid progenitors of the cultivated strawberry, establish a reliable backbone phylogeny for the genus, and provide genetic resources for molecular breeding.

Contrasting dynamics in abscisic acid metabolism in different Fragaria spp. during fruit ripening and identification of the enzymes involved

Abscisic acid (ABA) is a key hormone in non-climacteric Fragaria spp, regulating multiple physiological processes throughout fruit ripening. Its concentration increases during ripening, and it promotes fruit (receptacle) development. However, its metabolism in the fruit is largely unknown. We analyzed the concentrations of ABA and its catabolites at different developmental stages of strawberry ripening in diploid and octoploid genotypes and identified two functional ABA-glucosyltransferases (FvUGT71A49 and FvUGT73AC3) and two regiospecific ABA-8'-hydroxylases (FaCYP707A4a and FaCYP707A1/3). ABA-glucose ester content increased during ripening in diploid F. vesca varieties but decreased in octoploid F.×ananassa. Dihydrophaseic acid content increased throughout ripening in all analyzed receptacles, while 7'-hydroxy-ABA and neo-phaseic acid did not show significant changes during ripening. In the studied F. vesca varieties, the receptacle seems to be the main tissue for ABA metabolism, as the concentration of ABA and its metabolites in the receptacle was generally 100 times higher than in achenes. The accumulation patterns of ABA catabolites and transcriptomic data from the literature show that all strawberry fruits produce and metabolize considerable amounts of the plant hormone ABA during ripening, which is therefore a conserved process, but also illustrate the diversity of this metabolic pathway which is species, variety, and tissue dependent.

Molecular dating for phylogenies containing a mix of populations and species by using Bayesian and RelTime approaches

Simultaneous molecular dating of population and species divergences is essential in many biological investigations, including phylogeography, phylodynamics and species delimitation studies. In these investigations, multiple sequence alignments consist of both intra- and interspecies samples (mixed samples). As a result, the phylogenetic trees contain interspecies, interpopulation and within-population divergences. Bayesian relaxed clock methods are often employed in these analyses, but they assume the same tree prior for both inter- and intraspecies branching processes and require specification of a clock model for branch rates (independent vs. autocorrelated rates models). We evaluated the impact of a single tree prior on Bayesian divergence time estimates by analysing computer-simulated data sets. We also examined the effect of the assumption of independence of evolutionary rate variation among branches when the branch rates are autocorrelated. Bayesian approach with coalescent tree priors generally produced excellent molecular dates and highest posterior densities with high coverage probabilities. We also evaluated the performance of a non-Bayesian method, RelTime, which does not require the specification of a tree prior or a clock model. RelTime's performance was similar to that of the Bayesian approach, suggesting that it is also suitable to analyse data sets containing both populations and species variation when its computational efficiency is needed.

The high-quality genome of diploid strawberry (Fragaria nilgerrensis) provides new insights into anthocyanin accumulation

Fragaria nilgerrensis is a wild diploid strawberry species endemic to east and southeast region in Asia and provides a rich source of genetic variations for strawberry improvement. Here, we present a chromosome-scale assembly of F. nilgerrensis using single-molecule real-time (SMRT) Pacific Biosciences sequencing and chromosome conformation capture (Hi-C) genome scaffolding. The genome assembly size was 270.3 Mb, with a contig N50 of ∼8.5 Mb. A total of 28 780 genes and 117.2 Mb of transposable elements were annotated for this genome. Next, detailed comparative genomics with the high-quality F. vesca reference genome was conducted to obtain the difference among transposable elements, SNPs, Indels, and so on. The genome size of F. nilgerrensis was enhanced by around 50 Mb relatively to F. vesca, which is mainly due to expansion of transposable elements. In comparison with the F. vesca genome, we identified 4 561 825 SNPs, 846 301 Indels, 4243 inversions, 35 498 translocations and 10 099 relocations. We also found a marked expansion of genes involved in phenylpropanoid biosynthesis, starch and sucrose metabolism, cyanoamino acid metabolism, plant-pathogen interaction, brassinosteroid biosynthesis and plant hormone signal transduction in F. nilgerrensis, which may account for its specific phenotypes and considerable environmental adaptability. Interestingly, we found sequence variations in the upstream regulatory region of FnMYB10, a core transcriptional activator of anthocyanin biosynthesis, resulted in the low expression level of the FnMYB10 gene, which is likely responsible for white fruit phenotype of F. nilgerrensis. The high-quality F. nilgerrensis genome will be a valuable resource for biological research and comparative genomics research.

FaMYB9 is involved in the regulation of C6 volatile biosynthesis in strawberry

The large-scale untargeted proteomic and metabolomic studies were conducted in strawberry (Fragaria × ananassa) cv. Akihime fruit at five developmental stages. We found that some C6 volatiles highly contributed to the enrichment of volatiles at the red stage of strawberry fruit. We found that 12 genes involved in LOX pathway for volatile biosynthesis showed multiple patterns in protein abundance during fruit development and ripening, and 9 out of the 12 genes exhibited a significant increase in their relative expression levels at the red stage of fruit. We also found that the MYB9 gene (FaMYB9) expression level was positively correlated with the content of C6 volatiles (R = 0.989) and with the relative expression level and protein abundance of FaLOX5 at different strawberry fruit developmental stages (R = 0.954). The interaction between FaMYB9 and FaLOX5 was detected by yeast two-hybrid, co-immunoprecipitation (Co-IP), bimolecular fluorescence complementation (BiFC), and immunofluorescence (IF) analyses. Transient silencing of FaMYB9 delayed the fruit development and ripening, resulting in a significant decrease in the contents of C6 volatiles, while overexpression of FaMYB9 increased the fruit development and ripening and the contents of C6 volatiles in Akihime fruit. Therefore, FaMYB9 is positively involved in C6 volatile biosynthesis.

Genome Synteny Has Been Conserved Among the Octoploid Progenitors of Cultivated Strawberry Over Millions of Years of Evolution

Allo-octoploid cultivated strawberry (Fragaria × ananassa) originated through a combination of polyploid and homoploid hybridization, domestication of an interspecific hybrid lineage, and continued admixture of wild species over the last 300 years. While genes appear to flow freely between the octoploid progenitors, the genome structures and diversity of the octoploid species remain poorly understood. The complexity and absence of an octoploid genome frustrated early efforts to study chromosome evolution, resolve subgenomic structure, and develop a single coherent linkage group nomenclature. Here, we show that octoploid Fragaria species harbor millions of subgenome-specific DNA variants. Their diversity was sufficient to distinguish duplicated (homoeologous and paralogous) DNA sequences and develop 50K and 850K SNP genotyping arrays populated with co-dominant, disomic SNP markers distributed throughout the octoploid genome. Whole-genome shotgun genotyping of an interspecific segregating population yielded 1.9M genetically mapped subgenome variants in 5,521 haploblocks spanning 3,394 cM in F. chiloensis subsp. lucida, and 1.6M genetically mapped subgenome variants in 3,179 haploblocks spanning 2,017 cM in F. × ananassa. These studies provide a dense genomic framework of subgenome-specific DNA markers for seamlessly cross-referencing genetic and physical mapping information and unifying existing chromosome nomenclatures. Using comparative genomics, we show that geographically diverse wild octoploids are effectively diploidized, nearly completely collinear, and retain strong macro-synteny with diploid progenitor species. The preservation of genome structure among allo-octoploid taxa is a critical factor in the unique history of garden strawberry, where unimpeded gene flow supported its origin and domestication through repeated cycles of interspecific hybridization.

Mining and characterization of novel EST-SSR markers of Parrotia subaequalis (Hamamelidaceae) from the first Illumina-based transcriptome datasets

Parrotia subaequalis is an endangered Tertiary relict tree from eastern China. Despite its important ecological and horticultural value, no transcriptomic data and limited molecular markers are currently available in this species. In this study, we first performed high-throughput transcriptome sequencing of two individuals representing the northernmost (TX) and southernmost (SJD) population of P. subaequalis on the Illumina HiSeq 2500 platform. We gathered a total of 69,135 unigenes for P. subaequalis (TX) and 84,009 unigenes for P. subaequalis (SJD). From two unigenes datasets, 497 candidate polymorphic novel expressed sequence tag-simple sequence repeats (EST-SSRs) were identified using CandiSSR. Among these repeats, di-nucleotide repeats were the most abundant repeat type (62.78%) followed by tri-, tetra- and hexa-nucleotide repeats. We then randomly selected 54 primer pairs for polymorphism validation, of which 27 (50%) were successfully amplified and showed polymorphisms in 96 individuals from six natural populations of P. subaequalis. The average number of alleles per locus and the polymorphism information content values were 3.70 and 0.343; the average observed and expected heterozygosity were 0.378 and 0.394. A relatively high level of genetic diversity (H_T = 0.393) and genetic differentiation level (F_ST = 0.171) were surveyed, indicating P. subaequalis maintained high levels of species diversity in the long-term evolutionary history. Additionally, a high level of cross-transferability (92.59%) was displayed in five congeneric Hamamelidaceae species. Therefore, these new transcriptomic data and novel polymorphic EST-SSR markers will pinpoint genetic resources and facilitate future studies on population genetics and molecular breeding of P. subaequalis and other Hamamelidaceae species.

Functional trait divergence and trait plasticity confer polyploid advantage in heterogeneous environments

Polyploidy, or whole-genome duplication often with hybridization, is common in eukaryotes and is thought to drive ecological and evolutionary success, especially in plants. The mechanisms of polyploid success in ecologically relevant contexts, however, remain largely unknown. We conducted an extensive test of functional trait divergence and plasticity in conferring polyploid fitness advantage in heterogeneous environments, by growing clonal replicates of a worldwide genotype collection of six allopolyploid and five diploid wild strawberry (Fragaria) taxa in three climatically different common gardens. Among leaf functional traits, we detected divergence in trait means but not plasticities between polyploids and diploids, suggesting that increased genomic redundancy in polyploids does not necessarily translate into greater trait plasticity in response to environmental change. Across the heterogeneous garden environments, however, polyploids exhibited fitness advantage, which was conferred by both trait means and adaptive trait plasticities, supporting a 'jack-and-master' hypothesis for polyploids. Our findings elucidate essential ecological mechanisms underlying polyploid adaptation to heterogeneous environments, and provide an important insight into the prevalence and persistence of polyploid plants.

The Strawberry Tales: Size Matters

The history of the strawberry is one of ambitious kings and cunning spies, overseas adventures, and religious symbolism. Here, we will tell these stories, illustrated by relevant paintings and drawings, to give insight into the domestication of the strawberry. Furthermore, we will briefly discuss the complex strawberry genomes and aroma.

Plastid genomes reveal recurrent formation of allopolyploid Fragaria

PREMISE OF THE STUDY

Recurrent formation of polyploid taxa is a common observation in many plant groups. Haploid, cytoplasmic genomes like the plastid genome can be used to overcome the problem of homeologous genes and recombination in polyploid taxa. Fragaria (Rosaceae) contains several octo- and decaploid species. We use plastome sequences to infer the plastid ancestry of these taxa with special focus on the decaploid Fragaria cascadensis.

METHODS

We used genome skimming of 96 polyploid Fragaria samples on a single Illumina HiSeq 3000 lane to obtain whole plastome sequences. These sequences were used for phylogenetic reconstructions and dating analyses. Ploidy of all samples was inferred with flow cytometry, and plastid inheritance was examined in a controlled cross of F. cascadensis.

KEY RESULTS

The plastid genome phylogeny shows that only the octoploid F. chiloensis is monophyletic, all other polyploid taxa were supported to be para- or polyphyletic. The decaploid Fragaria cascadensis has biparental plastid inheritance and four different plastid donors. Diversification of the F. cascadensis clades occurred in the last 230,000 years. The southern part of its distribution range harbors considerably higher genetic diversity, suggestive of a potential refugium.

CONCLUSIONS

Fragaria cascadensis had at least four independent origins from parents with different plastomes. In contrast, para- and polyphyletic taxa of the octoploid Fragaria species are best explained by incomplete lineage sorting and/or hybridization. Biogeographic patterns in F. cascadensis are probably a result of range shift during the last glacial maximum.

Microevolution of the VQ gene family in six species of Fragaria

VQ motif-containing proteins play crucial roles in plant growth, development, and stress responses. However, no information of VQ motif-containing proteins has been studied at the microevolutionary level in species of Fragaria. In this study, a total of 19, 21, 23, 23, 23, and 25 genes containing the VQ motif were identified from the genomes of F. nipponica, F. iinumae, F. orientalis, F. vesca, F. nubicola, and F. x ananassa, respectively. We classified the VQ genes into 15 clades with grapevine VQ genes, which indicated that at least 15 ancient VQ genes existed before the divergence of the six studied species of Fragaria. Phylogenetic analysis indicated that 28 gene duplication events have occurred in the evolutionary process of the six species of Fragaria. Structural analysis showed that most of the VQ genes have no introns and that VQ proteins in each clade have a similar motif composition. The majority of gene pairs had K_a/K_s ratios less than 1, which illustrated that most of the VQ genes underwent purifying selection in the six species of Fragaria. Four types of cis-elements in promoters of VQ genes were detected, which is an important basis for further studies about plant stress responses. Furthermore, the expression analysis of FvVQ genes indicated that these genes are expressed differentially in the examined organs and tissues. The identification of VQ genes and the analysis of VQ gene duplication and polyploidization events in the six species of Fragaria provide important information on the evolutionary fate of VQ genes during the divergence of the six species of Fragaria.

Present-day sympatry belies the evolutionary origin of a high-order polyploid

Disentangling the evolutionary histories of polyploids, especially those with high ploidies, can reveal fundamental processes in speciation. Despite occurring frequently during evolution, the origins of many extant polyploid plant species remain largely unknown. By integrating linkage mapping, polyploid phylogeny and sex-determining region (SDR) in a unified framework, we statistically evaluated evolutionary hypotheses concerning the origin of a recently recognized decaploid strawberry (Fragaria cascadensis). The maximum-likelihood phylogenies and topology tests across homeologous groups consistently rejected the seemingly parsimonious hypothesis of 'contemporary sympatric speciation' via hybridization between octoploid and diploid congeners. Instead, most chromosomes supported 'ancient hybrid speciation' between a maternal octoploid progenitor ancestral to extant octoploid strawberries and a paternal, extinct Fragaria iinumae-like diploid progenitor, probably in Beringia during the Pleistocene. The absence of a shared SDR between the decaploid and other Fragaria is also consistent with an older origin rather than a recent hybrid origin in situ. Our study reveals a long evolutionary history of the decaploid despite its recent discovery, and highlights the pitfalls of inferring polyploid origins from niche/range alone or combined with morphology. It can serve as an exemplary starting step towards building much-needed model systems of established polyploids that have been, and remain to be, recognized.

Genotyping-by-sequencing enables linkage mapping in three octoploid cultivated strawberry families

Genotyping-by-sequencing (GBS) was used to survey genome-wide single-nucleotide polymorphisms (SNPs) in three biparental strawberry (Fragaria × ananassa) populations with the goal of evaluating this technique in a species with a complex octoploid genome. GBS sequence data were aligned to the F. vesca ‘Fvb’ reference genome in order to call SNPs. Numbers of polymorphic SNPs per population ranged from 1,163 to 3,190. Linkage maps consisting of 30–65 linkage groups were produced from the SNP sets derived from each parent. The linkage groups covered 99% of the Fvb reference genome, with three to seven linkage groups from a given parent aligned to any particular chromosome. A phylogenetic analysis performed using the POLiMAPS pipeline revealed linkage groups that were most similar to ancestral species F. vesca for each chromosome. Linkage groups that were most similar to a second ancestral species, F. iinumae, were only resolved for Fvb 4. The quantity of missing data and heterogeneity in genome coverage inherent in GBS complicated the analysis, but POLiMAPS resolved F. × ananassa chromosomal regions derived from diploid ancestor F. vesca.