Whole-genome sequencing of allotetraploid bermudagrass reveals the origin of Cynodon and candidate genes for salt tolerance

Bermudagrass (Cynodon dactylon) is a globally distributed, extensively used warm-season turf and forage grass with high tolerance to salinity and drought stress in alkaline environments. However, the origin of the species and genetic mechanisms for salinity tolerance in the species are basically unknown. Accordingly, we set out to study evolution divergence events in the Cynodon genome and to identify genes for salinity tolerance. We developed a 604.0 Mb chromosome-level polyploid genome sequence for bermudagrass 'A12359' (n = 18). The C. dactylon genome comprises 2 complete sets of homoeologous chromosomes, each with approximately 30 000 genes, and most genes are conserved as syntenic pairs. Phylogenetic study showed that the initial Cynodon species diverged from Oropetium thomaeum approximately 19.7-25.4 million years ago (Mya), the A and B subgenomes of C. dactylon diverged approximately 6.3-9.1 Mya, and the bermudagrass polyploidization event occurred 1.5 Mya on the African continent. Moreover, we identified 82 candidate genes associated with seven agronomic traits using a genome-wide association study, and three single-nucleotide polymorphisms were strongly associated with three salt resistance genes: RAP2-2, CNG channels, and F14D7.1. These genes may be associated with enhanced bermudagrass salt tolerance. These bermudagrass genomic resources, when integrated, may provide fundamental insights into evolution of diploid and tetraploid genomes and enhance the efficacy of comparative genomics in studying salt tolerance in Cynodon.

Does one subgenome become dominant in the formation and evolution of a polyploid?

BACKGROUND

Polyploids are common in flowering plants and they tend to have more expanded ranges of distributions than their diploid progenitors. Possible mechanisms underlying polyploid success have been intensively investigated. Previous studies showed that polyploidy generates novel changes and that subgenomes in allopolyploid species often differ in gene number, gene expression levels and levels of epigenetic alteration. It is widely believed that such differences are the results of conflicts among the subgenomes. These differences have been treated by some as subgenome dominance, and it is claimed that the magnitude of subgenome dominance increases in polyploid evolution.

SCOPE

In addition to changes which occurred during evolution, differences between subgenomes of a polyploid species may also be affected by differences between the diploid donors and changes which occurred during polyploidization. The variable genome components in many plant species are extensive, which would result in exaggerated differences between a subgenome and its progenitor when a single genotype or a small number of genotypes are used to represent a polyploid or its donors. When artificially resynthesized polyploids are used as surrogates for newly formed genotypes which have not been exposed to evolutionary selection, differences between diploid genotypes available today and those involved in the formation of the natural polyploid genotypes must also be considered.

CONCLUSIONS

Contrary to the now widely held views that subgenome biases in polyploids are the results of conflicts among the subgenomes and that one of the parental subgenomes generally retains more genes which are more highly expressed, available results show that subgenome biases mainly reflect legacy from the progenitors and that they can be detected before the completion of polyploidization events. Further, there is no convincing evidence that the magnitudes of subgenome biases have significantly changed during evolution for any of the allopolyploid species assessed.

Chromosome-scale genome assembly of the hexaploid Taiwanese goosefoot 'djulis' (Chenopodium formosanum)

Djulis (Chenopodium formosanum Koidz.) is a crop grown since antiquity in Taiwan. It is a BCD-genome hexaploid (2n = 6x = 54) domesticated form of lambsquarters (C. album L.) and a relative of the allotetraploid (AABB) C. quinoa. As with quinoa, djulis seed contains a complete protein profile and many nutritionally important vitamins and minerals. While still sold locally in Taiwanese markets, its traditional culinary uses are being lost as diets of younger generations change. Moreover, indigenous Taiwanese peoples who have long safeguarded djulis are losing their traditional farmlands. We used PacBio sequencing and Hi-C-based scaffolding to produce a chromosome-scale, reference-quality assembly of djulis. The final genome assembly spans 1.63 Gb in 798 scaffolds, with 97.8% of the sequence contained in 27 scaffolds representing the nine haploid chromosomes of each sub-genome of the species. BUSCO results indicated that 98.5% of the conserved orthologous genes for Viridiplantae are complete within the assembled genome, with 92.9% duplicated, as expected for a polyploid. A total of 67.8% of the assembly is repetitive, with the most common repeat being Gypsy long terminal repeat retrotransposons, which had significantly expanded in the B sub-genome. Gene annotation using Iso-Seq data from multiple tissues identified 75,056 putative gene models. Comparisons to quinoa showed strong patterns of synteny which allowed for the identification of homoeologous chromosomes, and sub-genome-specific sequences were used to assign homoeologs to each sub-genome. These results represent the first hexaploid genome assembly and the first assemblies of the C and D genomes of the Chenopodioideae subfamily.

Genome-Wide Duplication of Allotetraploid Brassica napus Produces Novel Characteristics and Extensive Ploidy Variation in Self-Pollinated Progeny

Whole genome duplications (WGDs) have played a major role in angiosperm species evolution. Polyploid plants have undergone multiple cycles of ancient WGD events during their evolutionary history. However, little attention has been paid to the additional WGD of the existing allopolyploids. In this study, we explored the influences of additional WGD on the allopolyploid Brassica napus. Compared to tetraploid B. napus, octoploid B. napus (AAAACCCC, 2n = 8x =76) showed significant differences in phenotype, reproductive ability and the ploidy of self-pollinated progeny. Genome duplication also altered a key reproductive organ feature in B. napus, that is, increased the number of pollen apertures. Unlike autopolyploids produced from the diploid Brassica species, the octoploid B. napus produced from allotetraploid B. napus had a relatively stable meiotic process, high pollen viability and moderate fertility under self-pollination conditions, indicating that sub-genomic interactions may be important for the successful establishment of higher-order polyploids. Doubling the genome of B. napus provided us with an opportunity to gain insight into the flexibility of the Brassica genomes. The genome size of self-pollinated progeny of octoploid B. napus varied greatly, and was accompanied by extensive genomic instability, such as aneuploidy, mixed-ploidy and mitotic abnormality. The octoploid B. napus could go through any of genome reduction, equilibrium or expansion in the short-term, thus providing a novel karyotype library for the Brassica genus. Our results reveal the short-term evolutionary consequences of recurrent polyploidization events, and help to deepen our understanding of polyploid plant evolution.

Expression patterns and functional divergence of homologous genes accompanied by polyploidization in cotton (Gossypium hirsutum L.)

Naturally allotetraploid cotton has been widely used as an ideal model to investigate gene expression remodeling as a consequence of polyploidization. However, the global gene pattern variation during early fiber development was unknown. In this study, through RNA-seq technology, we comprehensively investigated the expression patterns of homologous genes between allotetraploid cotton (G. hirsutum) and its diploid progenitors (G. arboreum and G. raimondii) at the fiber early development stage. In tetraploid cotton, genes showed expression level dominance (ELD) bias toward the A genome. This phenomenon was explained by the up-/downregulation of the homologs from the nondominant progenitor (D genome). Gene ontology (GO) enrichment results indicated that the ELD-A genes might be a prominent cause responsible for fiber property change through regulating the fatty acid biosynthesis/metabolism and microtubule procession, and the ELD-D genes might be involved in transcription regulation and stress inducement. In addition, the number and proportion of completely A- and D-subfunctionalized gene were similar at different fiber development stages. However, for neofunctionalization, the number and proportion of reactivated D-derived genes were greater than those of A at 3 and 5 DPA. Eventually, we found that some homologous genes belonging to several specific pathways might create novel asymmetric transcripts between two subgenomes during polyploidization and domestication process, further making the fiber property meet the human demands. Our study identified determinate pathways and their involved genes between allotetraploid cotton and their progenitors at early fiber development stages, providing new insights into the mechanism of cotton fiber evolution.

Evolutionary history and genetic diversity of apomictic allopolyploids in Hieracium s.str.: morphological versus genomic features

PREMISE

The origin of allopolyploids is believed to shape their evolutionary potential, ecology, and geographical ranges. Morphologically distinct apomictic types sharing the same parental species belong to the most challenging groups of polyploids. We evaluated the origins and variation of two triploid taxa (Hieracium pallidiflorum, H. picroides) presumably derived from the same diploid parental pair (H. intybaceum, H. prenanthoides).

METHODS

We used a suite of approaches ranging from morphological, phylogenetic (three unlinked molecular markers), and cytogenetic analyses (in situ hybridization) to genome size screening and genome skimming.

RESULTS

Genotyping proved the expected parentage of all analyzed accessions of H. pallidiflorum and H. picroides and revealed that nearly all of them originated independently. Genome sizes and genome dosage largely corresponded to morphology, whereas the maternal origin of the allopolyploids had no discernable effect. Polyploid accessions of both parental species usually contained genetic material from other species. Given the phylogenetic distance of the parents, their chromosomes appeared only weakly differentiated in genomic in situ hybridization (GISH), as well as in overall comparisons of the repetitive fraction of their genomes. Furthermore, the repeatome of a phylogenetically more closely related species (H. umbellatum) differed significantly more.

CONCLUSIONS

We proved (1) multiple origins of hybridogeneous apomicts from the same diploid parental taxa, and (2) allopolyploid origins of polyploid accessions of the parental species. We also showed that the evolutionary dynamics of very fast evolving markers such as satellite DNA or transposable elements does not necessarily follow patterns of speciation.

Genome skimming provides new insight into the relationships in Ludwigia section Macrocarpon, a polyploid complex

PREMISE OF THE STUDY

Interpreting relationships within groups containing polyploids, which are frequent in angiosperms, can be greatly assisted by genomic techniques. In this study, we used a genome-skimming approach to investigate the evolutionary relationships and origins of polyploids in the monophyletic group, Ludwigia section Macrocarpon (Onagraceae), which includes diploid, tetraploid, and hexaploid taxa.

METHODS

We sampled all known taxa and ploidy levels in the section and conducted shotgun sequencing. We assembled plastomes, mitochondrial sequences, and completed nuclear ribosomal regions, reconstructed phylogenies, and conducted comparative genomic analyses for plastomes to gain insights into the relationships among studied taxa.

KEY RESULTS

Within the section, results showed that the South American diploid taxa L. bonariensis and L. lagunae were closely related. We reported the first chromosome count (2n = 4× = 32) for L. neograndiflora, which is closely related to the two South American diploid taxa, although its exact origin remains unclear. The samples of the widespread, polyploid taxon L. octovalvis do not form a monophyletic group. Both tetraploid and hexaploid L. octovalvis lineages have originated more than once. At least one tetraploid in the L. octovalvis lineage may have been involved in the origins of hexaploids. One or more extinct/unsampled intermediate tetraploids in the L. octovalvis lineages had also likely been involved in the origins of hexaploids.

CONCLUSIONS

Genome skimming provided important insights into the complex evolutionary relationships within sect. Macrocarpon, but additional sampling and data from single-copy nuclear regions are necessary to further elucidate the origins of the polyploids in this section.

Isolation and characterization of microsatellite loci in Sorbus porrigentiformis and cross-amplification in S. aria and S. rupicola (Rosaceae)

PREMISE OF THE STUDY

Southwestern Britain is an emblematic hotspot of polyploid diversity of whitebeams (Sorbus aria agg.; Rosaceae) with ca. 30 polyploid endemic species. The tetraploid S. porrigentiformis is postulated as one of the parents of most of these endemics, along with the sexual diploid S. aria s. str. and the tetraploid S. rupicola.

METHODS AND RESULTS

We isolated 16 nuclear microsatellite loci from S. porrigentiformis and characterized them on 45 trees representing the three putative parental species. Eleven loci were polymorphic, and eight of them exhibited species-specific alleles. Allele numbers ranged from one to 11, and observed heterozygosity ranged from 0.40 to 1.00. The intraspecific levels of variation were very low, in agreement with the facultative apomictic reproduction hypothesized for this species.

CONCLUSIONS

The species-specific alleles will be useful for tracing the origin of the narrowly distributed Sorbus taxa. In addition, the assessment of diversity levels will help design a conservation strategy for the polyploid complex.

Evolutionary Conservation and Divergence of Gene Coexpression Networks in Gossypium (Cotton) Seeds

The cotton genus (Gossypium) provides a superior system for the study of diversification, genome evolution, polyploidization, and human-mediated selection. To gain insight into phenotypic diversification in cotton seeds, we conducted coexpression network analysis of developing seeds from diploid and allopolyploid cotton species and explored network properties. Key network modules and functional associations were identified related to seed oil content and seed weight. We compared species-specific networks to reveal topological changes, including rewired edges and differentially coexpressed genes, associated with speciation, polyploidy, and cotton domestication. Network comparisons among species indicate that topologies are altered in addition to gene expression profiles, indicating that changes in transcriptomic coexpression relationships play a role in the developmental architecture of cotton seed development. The global network topology of allopolyploids, especially for domesticated G. hirsutum, resembles the network of the A-genome diploid more than that of the D-genome parent, despite its D-like phenotype in oil content. Expression modifications associated with allopolyploidy include coexpression level dominance and transgressive expression, suggesting that the transcriptomic architecture in polyploids is to some extent a modular combination of that of its progenitor genomes. Among allopolyploids, intermodular relationships are more preserved between two different wild allopolyploid species than they are between wild and domesticated forms of a cultivated cotton, and regulatory connections of oil synthesis-related pathways are denser and more closely clustered in domesticated vs. wild G. hirsutum. These results demonstrate substantial modification of genic coexpression under domestication. Our work demonstrates how network inference informs our understanding of the transcriptomic architecture of phenotypic variation associated with temporal scales ranging from thousands (domestication) to millions (speciation) of years, and by polyploidy.

Multiple hybridization events in Cardamine (Brassicaceae) during the last 150 years: revisiting a textbook example of neoallopolyploidy

BACKGROUND AND AIMS

Recently formed allopolyploid species represent excellent subjects for exploring early stages of polyploid evolution. The hexaploid Cardamine schulzii was regarded as one of the few nascent allopolyploid species formed within the past ∼150 years that presumably arose by autopolyploidization of a triploid hybrid, C. × insueta; however, the most recent investigations have shown that it is a trigenomic hybrid. The aims of this study were to explore the efficiency of progenitor-specific microsatellite markers in detecting the hybrid origins and genome composition of these two allopolyploids, to estimate the frequency of polyploid formation events, and to outline their evolutionary potential for long-term persistence and speciation.

METHODS

Flow-cytometric ploidy-level screening and genotyping by progenitor-specific microsatellite markers (20 microsatellite loci) were carried out on samples focused on hybridizing populations at Urnerboden, Switzerland, but also including comparative material of the parental species from other sites in the Alps and more distant areas.

KEY RESULTS

It was confirmed that hybridization between the diploids C. amara and C. rivularis auct. gave rise to triploid C. × insueta, and it is inferred that this has occurred repeatedly. Evidence is provided that C. schulzii comprises three parental genomes and supports its origin from hybridization events between C. × insueta and the locally co-occurring hypotetraploid C. pratensis, leading to two cytotypes of C. schulzii: hypopentaploid and hypohexaploid. Each cytotype of C. schulzii is genetically uniform, suggesting their single origins.

CONCLUSIONS

Persistence of C. schulzii has presumably been achieved only by perennial growth and clonal reproduction. This contrasts with C. × insueta, in which multiple origins and occasional sexual reproduction have generated sufficient genetic variation for long-term survival and evolutionary success. This study illustrates a complex case of recurrent hybridization and polyploidization events, and highlights the role of triploids that promoted the origin of trigenomic hybrids.

Diversity and endemism in deglaciated areas: ploidy, relative genome size and niche differentiation in the Galium pusillum complex (Rubiaceae) in Northern and Central Europe

BACKGROUND AND AIMS

Plants endemic to areas covered by ice sheets during the last glaciation represent paradigmatic examples of rapid speciation in changing environments, yet very few systems outside the harsh arctic zone have been comprehensively investigated so far. The Galium pusillum aggregate (Rubiaceae) is a challenging species complex that exhibits a marked differentiation in boreal parts of Northern Europe. As a first step towards understanding its evolutionary history in deglaciated regions, this study assesses cytological variation and ecological preferences of the northern endemics and compares the results with corresponding data for species occurring in neighbouring unglaciated parts of Central and Western Europe.

METHODS

DNA flow cytometry was used together with confirmatory chromosome counts to determine ploidy levels and relative genome sizes in 1158 individuals from 181 populations. A formalized analysis of habitat preferences was applied to explore niche differentiation among species and ploidy levels.

KEY RESULTS

The G. pusillum complex evolved at diploid and tetraploid levels in Northern Europe, in contrast to the high-polyploid evolution of most other northern endemics. A high level of eco-geographic segregation was observed between different species (particularly along gradients of soil pH and competition) which is unusual for plants in deglaciated areas and most probably contributes to maintaining species integrity. Relative monoploid DNA contents of the species from previously glaciated regions were significantly lower than those of their counterparts from mostly unglaciated Central Europe, suggesting independent evolutionary histories.

CONCLUSIONS

The aggregate of G. pusillum in Northern Europe represents an exceptional case with a geographically vicariant and ecologically distinct diploid/tetraploid species endemic to formerly glaciated areas. The high level of interspecific differentiation substantially widens our perception of the evolutionary dynamics and speciation rates in the dramatically changing environments of Northern Europe.

The evolutionary history of the white-rayed species of Melampodium (Asteraceae) involved multiple cycles of hybridization and polyploidization

PREMISE OF THE STUDY

Polyploidy plays an important role in race differentiation and eventually speciation. Underlying mechanisms include chromosomal and genomic changes facilitating reproductive isolation and/or stabilization of hybrids. A prerequisite for studying these processes is a sound knowledge on the origin of polyploids. A well-suited group for studying polyploid evolution consists of the three species of Melampodium ser. Leucantha (Asteraceae): M. argophyllum, M. cinereum, and M. leucanthum.

METHODS

The origin of polyploids was inferred using network and tree-based phylogenetic analyses of several plastid and nuclear DNA sequences and of fingerprint data (AFLP). Genome evolution was assessed via genome size measurements, karyotype analysis, and in situ hybridization of ribosomal DNA.

KEY RESULTS

Tetraploid cytotypes of the phylogenetically distinct M. cinereum and M. leucanthum had, compared to the diploid cytotypes, doubled genome sizes and no evidence of gross chromosomal rearrangements. Hexaploid M. argophyllum constituted a separate lineage with limited intermixing with the other species, except in analyses from nuclear ITS. Its genome size was lower than expected if M. cinereum and/or M. leucanthum were involved in its origin, and no chromosomal rearrangements were evident.

CONCLUSIONS

Polyploids in M. cinereum and M. leucanthum are of recent autopolyploid origin in line with the lack of significant genomic changes. Hexaploid M. argophyllum also appears to be of autopolyploid origin against the previous hypothesis of an allopolyploid origin involving the other two species, but some gene flow with the other species in early phases of differentiation cannot be excluded.