Double trouble: taxonomy and definitions of polyploidy

Polyploid plants take cytonuclear perturbations in stride

Hybridization in plants is often accompanied by nuclear genome doubling (allopolyploidy), which has been hypothesized to perturb interactions between nuclear and organellar (mitochondrial and plastid) genomes by creating imbalances in the relative copy number of these genomes and producing genetic incompatibilities between maternally derived organellar genomes and the half of the allopolyploid nuclear genome from the paternal progenitor. Several evolutionary responses have been predicted to ameliorate these effects, including selection for changes in protein sequences that restore cytonuclear interactions; biased gene retention/expression/conversion favoring maternal nuclear gene copies; and fine-tuning of relative cytonuclear genome copy numbers and expression levels. Numerous recent studies, however, have found that evolutionary responses are inconsistent and rarely scale to genome-wide generalities. The apparent robustness of plant cytonuclear interactions to allopolyploidy may reflect features that are general to allopolyploids such as the lack of F2 hybrid breakdown under disomic inheritance, and others that are more plant-specific, including slow sequence divergence in organellar genomes and preexisting regulatory responses to changes in cell size and endopolyploidy during development. Thus, cytonuclear interactions may only rarely act as the main barrier to establishment of allopolyploid lineages, perhaps helping to explain why allopolyploidy is so pervasive in plant evolution.

Defining autopolyploidy: Cytology, genetics, and taxonomy

Autopolyploidy is taxonomically defined as the presence of more than two copies of each genome within an organism or species, where the genomes present must all originate within the same species. Alternatively, "genetic" or "cytological" autopolyploidy is defined by polysomic inheritance: random pairing and segregation of the four (or more) homologous chromosomes present, with no preferential pairing partners. In this review, we provide an overview of methods used to categorize species as taxonomic and cytological autopolyploids, including both modern and obsolete cytological methods, marker-segregation-based and genomics methods. Subsequently, we also investigated how frequently polysomic inheritance has been reliably documented in autopolyploids. Pure or predominantly polysomic inheritance was documented in 39 of 43 putative autopolyploid species where inheritance data was available (91%) and in seven of eight synthetic autopolyploids, with several cases of more mixed inheritance within species. We found no clear cases of autopolyploids with disomic inheritance, which was likely a function of our search methodology. Interestingly, we found seven species with purely polysomic inheritance and another five species with partial or predominant polysomic inheritance that appear to be taxonomic allopolyploids. Our results suggest that observations of polysomic inheritance can lead to relabeling of taxonomically allopolyploid species as autopolyploid and highlight the need for further cytogenetic and genomic investigation into polyploid origins and inheritance types.

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.

Autopolyploidy, Allopolyploidy, and Phylogenetic Networks with Horizontal Arcs

Polyploidization is an evolutionary process by which a species acquires multiple copies of its complete set of chromosomes. The reticulate nature of the signal left behind by it means that phylogenetic networks offer themselves as a framework to reconstruct the evolutionary past of species affected by it. The main strategy for doing this is to first construct a so-called multiple-labelled tree and to then somehow derive such a network from it. The following question therefore arises: How much can be said about that past if such a tree is not readily available? By viewing a polyploid dataset as a certain vector which we call a ploidy (level) profile, we show that among other results, there always exists a phylogenetic network in the form of a beaded phylogenetic tree with additional arcs that realizes a given ploidy profile. Intriguingly, the two end vertices of almost all of these additional arcs can be interpreted as having co-existed in time thereby adding biological realism to our network, a feature that is, in general, not enjoyed by phylogenetic networks. In addition, we show that our network may be viewed as a generator of ploidy profile space, a novel concept similar to phylogenetic tree space that we introduce to be able to compare phylogenetic networks that realize one and the same ploidy profile. We illustrate our findings in terms of a publicly available Viola dataset.

Neutral processes underlying the macro eco-evolutionary dynamics of mixed-ploidy systems

Polyploidy, i.e. the occurrence of multiple sets of chromosomes, is regarded as an important phenomenon in plant ecology and evolution, with all flowering plants likely having a polyploid ancestry. Owing to genome shock, minority cytotype exclusion and reduced fertility, polyploids emerging in diploid populations are expected to face significant challenges to successful establishment. Their establishment and persistence are often explained by possible fitness or niche differences that would relieve the competitive pressure with diploid progenitors. Experimental evidence for such advantages is, however, not unambiguous, and considerable niche overlap exists among most polyploid species and their diploid counterparts. Here, we develop a neutral spatially explicit eco-evolutionary model to understand whether neutral processes can explain the eco-evolutionary patterns of polyploids. We present a general mechanism for polyploid establishment by showing that sexually reproducing organisms assemble in space in an iterative manner, reducing frequency-dependent mating disadvantages and overcoming potential reduced fertility issues. Moreover, we construct a mechanistic theoretical framework that allows us to understand the long-term evolution of mixed-ploidy populations and show that our model is remarkably consistent with recent phylogenomic estimates of species extinctions in the Brassicaceae family.

Morphological and environmental differentiation as prezygotic reproductive barriers between parapatric and allopatric Campanula rotundifolia agg. cytotypes

BACKGROUND AND AIMS

Reproductive isolation and local establishment are necessary for plant speciation. Polyploidy, the possession of more than two complete chromosome sets, creates a strong postzygotic reproductive barrier between diploid and tetraploid cytotypes. However, this barrier weakens between polyploids (e.g. tetraploids and hexaploids). Reproductive isolation may be enhanced by cytotype morphological and environmental differentiation. Moreover, morphological adaptations to local conditions contribute to plant establishment. However, the relative contributions of ploidy level and the environment to morphology have generally been neglected. Thus, the extent of morphological variation driven by ploidy level and the environment was modelled for diploid, tetraploid and hexaploid cytotypes of Campanula rotundifolia agg. Cytotype distribution was updated, and morphological and environmental differentiation was tested in the presence and absence of natural contact zones.

METHODS

Cytotype distribution was assessed from 231 localities in Central Europe, including 48 localities with known chromosome counts, using flow cytometry. Differentiation in environmental niche and morphology was tested for cytotype pairs using discriminant analyses. A structural equation model was used to explore the synergies between cytotype, environment and morphology.

KEY RESULTS

Tremendous discrepancies were revealed between the reported and detected cytotype distribution. Neither mixed-ploidy populations nor interploidy hybrids were detected in the contact zones. Diploids had the broadest environmental niche, while hexaploids had the smallest and specialized niche. Hexaploids and spatially isolated cytotype pairs differed morphologically, including allopatric tetraploids. While leaf and shoot morphology were influenced by environmental conditions and polyploidy, flower morphology depended exclusively on the cytotype.

CONCLUSIONS

Reproductive isolation mechanisms vary between cytotypes. While diploids and polyploids are isolated postzygotically, the environmental niche shift is essential between higher polyploids. The impact of polyploidy and the environment on plant morphology implies the adaptive potential of polyploids, while the exclusive relationship between flower morphology and cytotype highlights the role of polyploidy in reproductive isolation.

Studying Whole-Genome Duplication Using Experimental Evolution of Chlamydomonas

In this chapter, we present the use of Chlamydomonas reinhardtii in experiments designed to study the evolutionary impacts of whole genome duplication. We shortly introduce the algal species and depict why it is an excellent model for experimental evolution. Subsequently, we discuss the most relevant steps and methods in the design of a ploidy-related Chlamydomonas experiment. These steps include strain selection, ploidy determination, different methods of making diplo- and polyploid Chlamydomonas cells, replication, culturing conditions, preservation, and the ways to quantify phenotypic and genotypic change.

Constraining Whole-Genome Duplication Events in Geological Time

The timing of whole-genome duplication (WGD) events is crucial to understanding their role in evolution and underpins many hypotheses linking WGD to increased diversity and complexity. As such, means of estimating the timing of the WGD events relative to their macroevolutionary outcomes are of considerable importance. Molecular clock methods facilitate direct estimation of the absolute timing of WGD events, integrating information on the rate of sequence evolution between species while accommodating the uncertainty inherent to the fossil record. We present an explanation of the best practice for constructing fossil calibrations and estimating the age of WGD events via molecular clock methods in the program MCMCtree, with an example dataset based on a well-characterized WGD event within the flowering dogwoods (Cornus). The approach presented herein allows for the estimation of the age of WGD events and subsequent speciation events, allowing the relationship between WGD and the macroevolutionary outcomes to be explored. In our example, we show that in the case of flowering dogwoods, the WGD event long predates the end-Cretaceous mass extinction and that the two events may be independent.

Brachypodium: 20 years as a grass biology model system; the way forward?

It has been 20 years since Brachypodium distachyon was suggested as a model grass species, but ongoing research now encompasses the entire genus. Extensive Brachypodium genome sequencing programmes have provided resources to explore the determinants and drivers of population diversity. This has been accompanied by cytomolecular studies to make Brachypodium a platform to investigate speciation, polyploidisation, perenniality, and various aspects of chromosome and interphase nucleus organisation. The value of Brachypodium as a functional genomic platform has been underscored by the identification of key genes for development, biotic and abiotic stress, and cell wall structure and function. While Brachypodium is relevant to the biofuel industry, its impact goes far beyond that as an intriguing model to study climate change and combinatorial stress.

Parental legacy versus regulatory innovation in salt stress responsiveness of allopolyploid cotton (Gossypium) species

Polyploidy provides an opportunity for evolutionary innovation and species diversification, especially under stressful conditions. In allopolyploids, the conditional dynamics of homoeologous gene expression can be either inherited from ancestral states pre-existing in the parental diploids or novel upon polyploidization, the latter potentially permitting a wider range of phenotypic responses to stresses. To gain insight into regulatory mechanisms underlying the diversity of salt resistance in Gossypium species, we compared global transcriptomic responses to modest salinity stress in two allotetraploid (AD-genome) cotton species, Gossypium hirsutum and G. mustelinum, relative to their model diploid progenitors (A-genome and D-genome). Multivariate and pairwise analyses of salt-responsive changes revealed a profound alteration of gene expression for about one third of the transcriptome. Transcriptional responses and associated functional implications of salt acclimation varied across species, as did species-specific coexpression modules among species and ploidy levels. Salt responsiveness in both allopolyploids was strongly biased toward the D-genome progenitor. A much lower level of transgressive downregulation was observed in the more salt-tolerant G. mustelinum than in the less tolerant G. hirsutum. By disentangling inherited effects from evolved responses, we show that expression biases that are not conditional upon salt stress approximately equally reflect parental legacy and regulatory novelty upon allopolyploidization, whereas stress-responsive biases are predominantly novel, or evolved, in allopolyploids. Overall, our work suggests that allopolyploid cottons acquired a wide range of stress response flexibility relative to their diploid ancestors, most likely mediated by complex suites of duplicated genes and regulatory factors.

Genomic constitution, allopolyploidy, and evolutionary proposal for Cynodon Rich. based on GISH

Polyploidy is the main mechanism for chromosome number variation in Cynodon. Taxonomic boundaries are difficult to define and, although phylogenetic studies indicate that some species are closely related, the degree of genomic similarity remains unknown. Furthermore, the Cynodon species classification as auto or allopolyploids is still controversial. Thus, this study aimed to investigate the genomic constitution in diploid and polyploid species using different approaches of genomic in situ hybridization (GISH). To better understand the hybridization events, we also investigated the occurrence of unreduced gametes in C. dactylon diploid pollen grains. We suggest a genomic nomenclature of diploid species as DD, D¹D¹, and D²D² for C. dactylon, C. incompletus, and C. nlemfuensis, and DDD²D² and DD²D¹D¹ for the segmental allotetraploids of Cynodon dactylon and C. transvaalensis, respectively. Furthermore, an evolutionary proposal was built based on our results and previous data from other studies, showing possible crosses that may have occurred between Cynodon species.

Global Patterns of subgenome evolution in organelle-targeted genes of six allotetraploid angiosperms

Whole-genome duplications (WGDs) are a prominent process of diversification in eukaryotes. The genetic and evolutionary forces that WGD imposes on cytoplasmic genomes are not well understood, despite the central role that cytonuclear interactions play in eukaryotic function and fitness. Cellular respiration and photosynthesis depend on successful interaction between the 3,000+ nuclear-encoded proteins destined for the mitochondria or plastids and the gene products of cytoplasmic genomes in multi-subunit complexes such as OXPHOS, organellar ribosomes, Photosystems I and II, and Rubisco. Allopolyploids are thus faced with the critical task of coordinating interactions between the nuclear and cytoplasmic genes that were inherited from different species. Because the cytoplasmic genomes share a more recent history of common descent with the maternal nuclear subgenome than the paternal subgenome, evolutionary “mismatches” between the paternal subgenome and the cytoplasmic genomes in allopolyploids might lead to the accelerated rates of evolution in the paternal homoeologs of allopolyploids, either through relaxed purifying selection or strong directional selection to rectify these mismatches. We report evidence from six independently formed allotetraploids that the subgenomes exhibit unequal rates of protein-sequence evolution, but we found no evidence that cytonuclear incompatibilities result in altered evolutionary trajectories of the paternal homoeologs of organelle-targeted genes. The analyses of gene content revealed mixed evidence for whether the organelle-targeted genes are lost more rapidly than the non-organelle-targeted genes. Together, these global analyses provide insights into the complex evolutionary dynamics of allopolyploids, showing that the allopolyploid subgenomes have separate evolutionary trajectories despite sharing the same nucleus, generation time, and ecological context.

Chromosome painting reveals genomic structure of three polyploid species of Ipomoea

Cultivated sweetpotato (Ipomoea batatas (L.) Lam.) from the family Convolvulaceae is a hexaploid species with 2n = 6x = 90, and has been controversial regarding its nature as an autopolyploid arising within a species or allopolyploid forming between species. Here, we developed oligonucleotide-based painting probes for two chromosomes of I.nil, a model diploid Ipomoea species. Using these probes, we revealed pairing behavior of homoeologous chromosomes in I. batatas and its two possible polyploid ancestral species, tetraploid I. tabascana (2n = 4x = 60) and hexaploid I. trifida (2n = 6x = 90). Chromosome painting analysis revealed a high percentage of quadrivalent formation in zygotene-pachytene cells of I. tabascana, which supported that I. tabascana was an autotetraploid likely derived by doubling of structurally similar and homologous genomes rather than a hybrid between I. batatas and I. trifida (2x). The high frequency of hexavalent/bivalent and tetravalent pairing was observed in I. trifida (6x) and I. batatas. However, the percentage of hexavalent pairing in I. trifida (6x) was far higher than that in I. batatas. Thus, present results tended to support that I. trifida (6x) was an autohexaploid, while I. batatas was more likely to be a segmental allohexaploid.

Evidence for hexasomic inheritance in Chrysanthemum morifolium Ramat. based on analysis of EST-SSR markers

Chrysanthemums (Chrysanthemum morifolium Ramat.) are ornamental flowers, which are famous worldwide. The mode of inheritance has great implications for the genetic analysis of polyploid species. However, genetic analysis of chrysanthemum has been hampered because of its controversial inheritance mode (disomic or hexasomic). To classify the inheritance mode of chrysanthemums, an analysis of three approaches was carried out in an F₁ progeny of 192 offspring using 223 expressed sequence tag-simple sequence repeat (EST-SSR) markers. The analysis included segregation analysis, the ratio of simplex marker alleles linked in coupling to repulsion, as well as the transmission and segregation patterns of EST-SSR marker alleles. After segregation analysis, 204 marker alleles fit hexasomic inheritance and 150 marker alleles fit disomic inheritance, showing that marker alleles were inherited predominantly in a hexasomic manner. Furthermore, the results of the analysis of allele configuration and segregation behavior of five EST-SSR markers also suggested random pairing of chromosomes. Additionally, the ratio of simplex marker alleles linked in coupling to repulsion was 1:0, further supporting hexasomic inheritance. Therefore, it could be inferred that chrysanthemum is a complete or near-complete hexasome.

High-Throughput Genomic Data Reveal Complex Phylogenetic Relationships in Stylosanthes Sw (Leguminosae)

Allopolyploidy is widely present across plant lineages. Though estimating the correct phylogenetic relationships and origin of allopolyploids may sometimes become a hard task. In the genus Stylosanthes Sw. (Leguminosae), an important legume crop, allopolyploidy is a key speciation force. This makes difficult adequate species recognition and breeding efforts on the genus. Based on comparative analysis of nine high-throughput sequencing (HTS) samples, including three allopolyploids (S. capitata Vogel cv. “Campo Grande,” S. capitata “RS024” and S. scabra Vogel) and six diploids (S. hamata Taub, S. viscosa (L.) Sw., S. macrocephala M. B. Ferreira and Sousa Costa, S. guianensis (Aubl.) Sw., S. pilosa M. B. Ferreira and Sousa Costa and S. seabrana B. L. Maass & 't Mannetje) we provide a working pipeline to identify organelle and nuclear genome signatures that allowed us to trace the origin and parental genome recognition of allopolyploids. First, organelle genomes were de novo assembled and used to identify maternal genome donors by alignment-based phylogenies and synteny analysis. Second, nuclear-derived reads were subjected to repetitive DNA identification with RepeatExplorer2. Identified repeats were compared based on abundance and presence on diploids in relation to allopolyploids by comparative repeat analysis. Third, reads were extracted and grouped based on the following groups: chloroplast, mitochondrial, satellite DNA, ribosomal DNA, repeat clustered- and total genomic reads. These sets of reads were then subjected to alignment and assembly free phylogenetic analyses and were compared to classical alignment-based phylogenetic methods. Comparative analysis of shared and unique satellite repeats also allowed the tracing of allopolyploid origin in Stylosanthes, especially those with high abundance such as the StyloSat1 in the Scabra complex. This satellite was in situ mapped in the proximal region of the chromosomes and made it possible to identify its previously proposed parents. Hence, with simple genome skimming data we were able to provide evidence for the recognition of parental genomes and understand genome evolution of two Stylosanthes allopolyploids.

Reciprocal allopolyploid grasses (Festuca × Lolium) display stable patterns of genome dominance

Allopolyploidization entailing the merger of two distinct genomes in a single hybrid organism, is an important process in plant evolution and a valuable tool in breeding programs. Newly established hybrids often experience massive genomic perturbations, including karyotype reshuffling and gene expression modifications. These phenomena may be asymmetric with respect to the two progenitors, with one of the parental genomes being "dominant." Such "genome dominance" can manifest in several ways, including biased homoeolog gene expression and expression level dominance. Here we employed a k-mer-based approach to study gene expression in reciprocal Festuca pratensis Huds. × Lolium multiflorum Lam. allopolyploid grasses. Our study revealed significantly more genes where expression mimicked that of the Lolium parent compared with the Festuca parent. This genome dominance was heritable to successive generation and its direction was only slightly modified by environmental conditions and plant age. Our results suggest that Lolium genome dominance was at least partially caused by its more efficient trans-acting gene expression regulatory factors. Unraveling the mechanisms responsible for propagation of parent-specific traits in hybrid crops contributes to our understanding of allopolyploid genome evolution and opens a way to targeted breeding strategies.

Patterns and Processes of Diploidization in Land Plants

Most land plants are now known to be ancient polyploids that have rediploidized. Diploidization involves many changes in genome organization that ultimately restore bivalent chromosome pairing and disomic inheritance, and resolve dosage and other issues caused by genome duplication. In this review, we discuss the nature of polyploidy and its impact on chromosome pairing behavior. We also provide an overview of two major and largely independent processes of diploidization: cytological diploidization and genic diploidization/fractionation. Finally, we compare variation in gene fractionation across land plants and highlight the differences in diploidization between plants and animals. Altogether, we demonstrate recent advancements in our understanding of variation in the patterns and processes of diploidization in land plants and provide a road map for future research to unlock the mysteries of diploidization and eukaryotic genome evolution.

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.

Gradual polyploid genome evolution revealed by pan-genomic analysis of Brachypodium hybridum and its diploid progenitors

Our understanding of polyploid genome evolution is constrained because we cannot know the exact founders of a particular polyploid. To differentiate between founder effects and post polyploidization evolution, we use a pan-genomic approach to study the allotetraploid Brachypodium hybridum and its diploid progenitors. Comparative analysis suggests that most B. hybridum whole gene presence/absence variation is part of the standing variation in its diploid progenitors. Analysis of nuclear single nucleotide variants, plastomes and k-mers associated with retrotransposons reveals two independent origins for B. hybridum, ~1.4 and ~0.14 million years ago. Examination of gene expression in the younger B. hybridum lineage reveals no bias in overall subgenome expression. Our results are consistent with a gradual accumulation of genomic changes after polyploidization and a lack of subgenome expression dominance. Significantly, if we did not use a pan-genomic approach, we would grossly overestimate the number of genomic changes attributable to post polyploidization evolution.

In Support of Winge's Theory of "Hybridization Followed by Chromosome Doubling"

Polyploidy-or chromosome doubling-plays a significant role in plant speciation and evolution. Much of the existing evidence indicates that fusion of unreduced (or 2n) gametes is the major pathway responsible for polyploid formation. In the early 1900s, a theory was put forward that the mechanism of "hybridization followed by chromosome doubling" would enable the survival and development of the hybrid zygote by providing each chromosome with a homolog with which to pair. However, to date there is only scant empirical evidence supporting this theory. In our previous study, interspecific-interploid crosses between the tetraploid Hylocereus megalanthus, as the female parent, and the diploid H. undatus, as the male parent, yielded only allopentaploids, allohexaploids, and 5x-and 6x-aneuploids instead of the expected allotriploids. No viable hybrids were obtained from the reciprocal cross. Since H. undatus underwent normal meiosis with regular pairing in the pollen mother cells and only reduced pollen grains were observed, the allohexaploids obtained supported the concept of "chromosome doubling." In this work, we report ploidy level, fruit morphology, and pollen viability and diameter in a group of putative hybrids obtained from an embryo rescue procedure following controlled H. megalanthus × H. undatus crosses, with the aim to elucidate, for the first time, the timing and developmental stage of the chromosome doubling. As in our previous report, no triploids were obtained, but tetraploids, pentaploids, hexaploids, and 5x- and 6x-aneuploids were found in the regenerated plants. The tetraploids exhibited the morphological features of the maternal parent and could not be considered true hybrids. Based on our previous studies, we can assume that the pentaploids were a result of a fertilization event between one unreduced (2n) female gamete from the tetraploid H. megalanthus and a normal (n) haploid male gamete from H. undatus. All the allohexaploids obtained from the embryo rescue technique where those that regenerated from fertilized ovules 10 days after pollination (at the pro-embryo stage), showing that the chromosome doubling event occurred at a very early development stage, i.e., at the zygote stage or shortly after zygote formation. These allohexaploids thus constitute empirical evidence of "hybridization followed by chromosome doubling."

Comparatively Barcoded Chromosomes of Brachypodium Perennials Tell the Story of Their Karyotype Structure and Evolution

The Brachypodium genus is an informative model system for studying grass karyotype organization. Previous studies of a limited number of species and reference chromosomes have not provided a comprehensive picture of the enigmatic phylogenetic relationships in the genus. Comparative chromosome barcoding, which enables the reconstruction of the evolutionary history of individual chromosomes and their segments, allowed us to infer the relationships between putative ancestral karyotypes of extinct species and extant karyotypes of current species. We used over 80 chromosome-specific BAC (bacterial artificial chromosome) clones derived from five reference chromosomes of B. distachyon as probes against the karyotypes of twelve accessions representing five diploid and polyploid Brachypodium perennials. The results showed that descending dysploidy is common in Brachypodium and occurs primarily via nested chromosome fusions. Brachypodiumdistachyon was rejected as a putative ancestor for allotetraploid perennials and B. stacei for B. mexicanum. We propose two alternative models of perennial polyploid evolution involving either the incorporation of a putative x = 5 ancestral karyotype with different descending dysploidy patterns compared to B. distachyon chromosomes or hybridization of two x = 9 ancestors followed by genome doubling and descending dysploidy. Details of the karyotype structure and evolution in several Brachypodium perennials are revealed for the first time.

Chromosome painting in meiosis reveals pairing of specific chromosomes in polyploid Solanum species

Analysis of chromosome pairing has been an important tool to assess the genetic similarity of homologous and homoeologous chromosomes in polyploids. However, it is technically challenging to monitor the pairing of specific chromosomes in polyploid species, especially for plant species with a large number of small chromosomes. We developed oligonucleotide-based painting probes for four different potato chromosomes. We demonstrate that these probes are robust enough to monitor a single chromosome throughout the prophase I of meiosis in polyploid Solanum species. Cultivated potato (Solanum tuberosum, 2n = 4x = 48) is an autotetraploid. We demonstrate that the four copies of each potato chromosome pair as a quadrivalent in 66-78% of the meiotic cells at the pachytene stage. Solanum demissum (2n = 6x = 72) is a hexaploid and has been controversial regarding its nature as an autopolyploid or allopolyploid. Interestingly, no hexavalent pairing was observed in meiosis. Instead, we observed three independent bivalents in 83-98% of the meiotic cells at late diakinesis and early metaphase I for the four chromosomes. These results suggest that S. demissum has evolved into a cytologically stable state with predominantly bivalent pairing in meiosis.

Tools for Genetic Studies in Experimental Populations of Polyploids

Polyploid organisms carry more than two copies of each chromosome, a condition rarely tolerated in animals but which occurs relatively frequently in the plant kingdom. One of the principal challenges faced by polyploid organisms is to evolve stable meiotic mechanisms to faithfully transmit genetic information to the next generation upon which the study of inheritance is based. In this review we look at the tools available to the research community to better understand polyploid inheritance, many of which have only recently been developed. Most of these tools are intended for experimental populations (rather than natural populations), facilitating genomics-assisted crop improvement and plant breeding. This is hardly surprising given that a large proportion of domesticated plant species are polyploid. We focus on three main areas: (1) polyploid genotyping; (2) genetic and physical mapping; and (3) quantitative trait analysis and genomic selection. We also briefly review some miscellaneous topics such as the mode of inheritance and the availability of polyploid simulation software. The current polyploid analytic toolbox includes software for assigning marker genotypes (and in particular, estimating the dosage of marker alleles in the heterozygous condition), establishing chromosome-scale linkage phase among marker alleles, constructing (short-range) haplotypes, generating linkage maps, performing genome-wide association studies (GWAS) and quantitative trait locus (QTL) analyses, and simulating polyploid populations. These tools can also help elucidate the mode of inheritance (disomic, polysomic or a mixture of both as in segmental allopolyploids) or reveal whether double reduction and multivalent chromosomal pairing occur. An increasing number of polyploids (or associated diploids) are being sequenced, leading to publicly available reference genome assemblies. Much work remains in order to keep pace with developments in genomic technologies. However, such technologies also offer the promise of understanding polyploid genomes at a level which hitherto has remained elusive.

Ten steps to get started in Genome Assembly and Annotation

As a part of the ELIXIR-EXCELERATE efforts in capacity building, we present here 10 steps to facilitate researchers getting started in genome assembly and genome annotation. The guidelines given are broadly applicable, intended to be stable over time, and cover all aspects from start to finish of a general assembly and annotation project.

Intrinsic properties of genomes are discussed, as is the importance of using high quality DNA. Different sequencing technologies and generally applicable workflows for genome assembly are also detailed. We cover structural and functional annotation and encourage readers to also annotate transposable elements, something that is often omitted from annotation workflows. The importance of data management is stressed, and we give advice on where to submit data and how to make your results Findable, Accessible, Interoperable, and Reusable (FAIR).

Phylogenomics of the olive tree (Olea europaea) reveals the relative contribution of ancient allo- and autopolyploidization events

BACKGROUND

Polyploidization is one of the major evolutionary processes that shape eukaryotic genomes, being particularly common in plants. Polyploids can arise through direct genome doubling within a species (autopolyploidization) or through the merging of genomes from distinct species after hybridization (allopolyploidization). The relative contribution of both mechanisms in plant evolution is debated. Here we used phylogenomics to dissect the tempo and mode of duplications in the genome of the olive tree (Olea europaea), one of the first domesticated Mediterranean fruit trees.

RESULTS

Our results depict a complex scenario involving at least three past polyploidization events, of which two-at the bases of the family Oleaceae and the tribe Oleeae, respectively-are likely to be the result of ancient allopolyploidization. A more recent polyploidization involves specifically the olive tree and relatives.

CONCLUSION

Our results show the power of phylogenomics to distinguish between allo- and auto polyploidization events and clarify the contributions of duplications in the evolutionary history of the olive tree.

Karyotypic variation and pollen stainability in resynthesized allopolyploids Tragopogon miscellus and T. mirus

PREMISE OF THE STUDY

Polyploidy has extensively shaped the evolution of plants, but the early stages of polyploidy are still poorly understood. The neoallopolyploid species Tragopogon mirus and T. miscellus are both characterized by widespread karyotypic variation, including frequent aneuploidy and intergenomic translocations. Our study illuminates the origins and early impacts of this variation by addressing two questions: How quickly does karyotypic variation accumulate in Tragopogon allopolyploids following whole-genome duplication (WGD), and how does the fertility of resynthesized Tragopogon allopolyploids evolve shortly after WGD?

METHODS

We used genomic in situ hybridization and lactophenol-cotton blue staining to estimate the karyotypic variation and pollen stainability, respectively, of resynthesized T. mirus and T. miscellus during the first five generations after WGD.

KEY RESULTS

Widespread karyotypic variation developed quickly in synthetics and resembled that of naturally occurring T. mirus and T. miscellus by generation S₄ . Pollen stainability in resynthesized allopolyploids was consistently lower than that of natural T. mirus and T. miscellus, as well as their respective diploid progenitor species. Logistic regression showed that mean pollen stainability increased slightly over four generations in resynthesized T. mirus but remained at equivalent levels in T. miscellus.

CONCLUSIONS

Our results clarify some of the changes that occur in T. mirus and T. miscellus immediately following their origin, most notably the rapid onset of karyotypic variation within these species immediately following WGD.

Conclusive evidence for hexasomic inheritance in chrysanthemum based on analysis of a 183 k SNP array

BACKGROUND

Cultivated chrysanthemum is an outcrossing hexaploid (2n = 6× = 54) with a disputed mode of inheritance. In this paper, we present a single nucleotide polymorphism (SNP) selection pipeline that was used to design an Affymetrix Axiom array with 183 k SNPs from RNA sequencing data (1). With this array, we genotyped four bi-parental populations (with sizes of 405, 53, 76 and 37 offspring plants respectively), and a cultivar panel of 63 genotypes. Further, we present a method for dosage scoring in hexaploids from signal intensities of the array based on mixture models (2) and validation of selection steps in the SNP selection pipeline (3). The resulting genotypic data is used to draw conclusions on the mode of inheritance in chrysanthemum (4), and to make an inference on allelic expression bias (5).

RESULTS

With use of the mixture model approach, we successfully called the dosage of 73,936 out of 183,130 SNPs (40.4%) that segregated in any of the bi-parental populations. To investigate the mode of inheritance, we analysed markers that segregated in the large bi-parental population (n = 405). Analysis of segregation of duplex x nulliplex SNPs resulted in evidence for genome-wide hexasomic inheritance. This evidence was substantiated by the absence of strong linkage between markers in repulsion, which indicated absence of full disomic inheritance. We present the success rate of SNP discovery out of RNA sequencing data as affected by different selection steps, among which SNP coverage over genotypes and use of different types of sequence read mapping software. Genomic dosage highly correlated with relative allele coverage from the RNA sequencing data, indicating that most alleles are expressed according to their genomic dosage.

CONCLUSIONS

The large population, genotyped with a very large number of markers, is a unique framework for extensive genetic analyses in hexaploid chrysanthemum. As starting point, we show conclusive evidence for genome-wide hexasomic inheritance.

Partial preferential chromosome pairing is genotype dependent in tetraploid rose

It has long been recognised that polyploid species do not always neatly fall into the categories of auto- or allopolyploid, leading to the term 'segmental allopolyploid' to describe everything in between. The meiotic behaviour of such intermediate species is not fully understood, nor is there consensus as to how to model their inheritance patterns. In this study we used a tetraploid cut rose (Rosa hybrida) population, genotyped using the 68K WagRhSNP array, to construct an ultra-high-density linkage map of all homologous chromosomes using methods previously developed for autotetraploids. Using the predicted bivalent configurations in this population we quantified differences in pairing behaviour among and along homologous chromosomes, leading us to correct our estimates of recombination frequency to account for this behaviour. This resulted in the re-mapping of 25 695 SNP markers across all homologues of the seven rose chromosomes, tailored to the pairing behaviour of each chromosome in each parent. We confirmed the inferred differences in pairing behaviour among chromosomes by examining repulsion-phase linkage estimates, which also carry information about preferential pairing and recombination. Currently, the closest sequenced relative to rose is Fragaria vesca. Aligning the integrated ultra-dense rose map with the strawberry genome sequence provided a detailed picture of the synteny, confirming overall co-linearity but also revealing new genomic rearrangements. Our results suggest that pairing affinities may vary along chromosome arms, which broadens our current understanding of segmental allopolyploidy.