Diverse genome organization following 13 independent mesopolyploid events in Brassicaceae contrasts with convergent patterns of gene retention

Biased Retention of Environment-Responsive Genes Following Genome Fractionation

The molecular underpinnings and consequences of cycles of whole-genome duplication (WGD) and subsequent gene loss through subgenome fractionation remain largely elusive. Endogenous drivers, such as transposable elements (TEs), have been postulated to shape genome-wide dominance and biased fractionation, leading to a conserved least-fractionated (LF) subgenome and a degenerated most-fractionated (MF) subgenome. In contrast, the role of exogenous factors, such as those induced by environmental stresses, has been overlooked. In this study, a chromosome-scale assembly of the alpine buckler mustard (Biscutella laevigata; Brassicaceae) that underwent a WGD event about 11 million years ago is coupled with transcriptional responses to heat, cold, drought, and herbivory to assess how gene expression is associated with differential gene retention across the MF and LF subgenomes. Counteracting the impact of TEs in reducing the expression and retention of nearby genes across the MF subgenome, dosage balance is highlighted as a main endogenous promoter of the retention of duplicated gene products under purifying selection. Consistent with the "turn a hobby into a job" model, about one-third of environment-responsive duplicates exhibit novel expression patterns, with one copy typically remaining conditionally expressed, whereas the other copy has evolved constitutive expression, highlighting exogenous factors as a major driver of gene retention. Showing uneven patterns of fractionation, with regions remaining unbiased, but with others showing high bias and significant enrichment in environment-responsive genes, this mesopolyploid genome presents evolutionary signatures consistent with an interplay of endogenous and exogenous factors having driven gene content following WGD-fractionation cycles.

Polyploids of Brassicaceae: Genomic Insights and Assembly Strategies

The Brassicaceae family is distinguished by its inclusion of high-value crops such as cabbage, broccoli, mustard, and wasabi, all noted for their glucosinolates. In this family, many polyploidy species are distributed and shaped by numerous whole-genome duplications, independent genome doublings, and hybridization events. The evolutionary trajectory of the family is marked by enhanced diversification and lineage splitting after paleo- and meso-polyploidization, with discernible remnants of whole-genome duplications within their genomes. The recent neopolyploidization events notably increased the proportion of polyploid species within the family. Although sequencing efforts for the Brassicaceae genome have been robust, accurately distinguishing sub-genomes remains a significant challenge, frequently complicating the assembly process. Assembly strategies include comparative analyses with ancestral species and examining k-mers, long terminal repeat retrotransposons, and pollen sequencing. This review comprehensively explores the unique genomic characteristics of the Brassicaceae family, with a particular emphasis on polyploidization events and the latest strategies for sequencing and assembly. This review will significantly improve our understanding of polyploidy in the Brassicaceae family and assist in future genome assembly methods.

Complementing model species with model clades

Model species continue to underpin groundbreaking plant science research. At the same time, the phylogenetic resolution of the land plant tree of life continues to improve. The intersection of these 2 research paths creates a unique opportunity to further extend the usefulness of model species across larger taxonomic groups. Here we promote the utility of the Arabidopsis thaliana model species, especially the ability to connect its genetic and functional resources, to species across the entire Brassicales order. We focus on the utility of using genomics and phylogenomics to bridge the evolution and diversification of several traits across the Brassicales to the resources in Arabidopsis, thereby extending scope from a model species by establishing a "model clade." These Brassicales-wide traits are discussed in the context of both the model species Arabidopsis and the family Brassicaceae. We promote the utility of such a "model clade" and make suggestions for building global networks to support future studies in the model order Brassicales.

The evolution of ephemeral flora in Xinjiang, China: insights from plastid phylogenomic analyses of Brassicaceae

BACKGROUND

The ephemeral flora of northern Xinjiang, China, plays an important role in the desert ecosystems. However, the evolutionary history of this flora remains unclear. To gain new insights into its origin and evolutionary dynamics, we comprehensively sampled ephemeral plants of Brassicaceae, one of the essential plant groups of the ephemeral flora.

RESULTS

We reconstructed a phylogenetic tree using plastid genomes and estimated their divergence times. Our results indicate that ephemeral species began to colonize the arid areas in north Xinjiang during the Early Miocene and there was a greater dispersal of ephemeral species from the surrounding areas into the ephemeral community of north Xinjiang during the Middle and Late Miocene, in contrast to the Early Miocene or Pliocene periods.

CONCLUSIONS

Our findings, together with previous studies, suggest that the ephemeral flora originated in the Early Miocene, and species assembly became rapid from the Middle Miocene onwards, possibly attributable to global climate changes and regional geological events.

Interploidy gene flow via a 'pentaploid bridge' and ploidy reduction in Cystopteris fragilis fern complex (Cystopteridaceae: Polypodiales)

KEY MESSAGE

Our results indicate the existence of interploidy gene flow in Cystopteris fragilis, resulting in sexual triploid and diploid gametophytes from pentaploid parents. Similar evolutionary dynamics might operate in other fern complexes and need further investigation. Polyploidization and hybridization are a key evolutionary processes in ferns. Here, we outline an interploidy gene flow pathway operating in the polyploid Cystopteris fragilis complex. The conditions necessary for the existence of this pathway were tested. A total of 365 C. fragilis individuals were collected covering representatives of all three predominant ploidy levels (tetraploid, pentaploid, and hexaploid), cultivated, had their ploidy level estimated by flow cytometry, and their spores collected. The spores, as well as gametophytes and sporophytes established from them, were analysed by flow cytometry. Spore abortion rate was also estimated. In tetraploids, we observed the formation of unreduced (tetraploid) spores (ca 2%). Collected pentaploid individuals indicate ongoing hybridization between ploidy levels. Pentaploids formed up to 52% viable spores, ca 79% of them reduced, i.e. diploid and triploid. Reduced spores formed viable gametophytes, and, in the case of triploids, filial hexaploid sporophytes, showing evidence of sexual reproduction. Some tetraploid sporophytes reproduce apomictically (based on uniform ploidy of their metagenesis up to filial sporophytes). Triploid and diploid gametophytes from pentaploid parents are able to mate among themselves, or with "normal" reduced gametophytes from the sexual tetraploid sporophytes (the dominant ploidy level in the sporophytes in this populations), to produce tetraploid, pentaploid, and hexaploid sporophytes, allowing for geneflow from the pentaploids to both the tetraploid and hexaploid populations. Similar evolutionary dynamics might operate in other fern complexes and need further investigation.

Plant biology: Phylogenomics of mustards and their relatives

A new nuclear phylogeny for the large family Brassicaceae opens the way for advanced comparative studies of adaptation, development, coevolution, hybridization, and diversification in this crucial group, which is the source for many of the genomic resources now used across the flowering plants.

Global Brassicaceae phylogeny based on filtering of 1,000-gene dataset

The mustard family (Brassicaceae) is a scientifically and economically important family, containing the model plant Arabidopsis thaliana and numerous crop species that feed billions worldwide. Despite its relevance, most phylogenetic trees of the family are incompletely sampled and often contain poorly supported branches. Here, we present the most complete Brassicaceae genus-level family phylogenies to date (Brassicaceae Tree of Life or BrassiToL) based on nuclear (1,081 genes, 319 of the 349 genera; 57 of the 58 tribes) and plastome (60 genes, 265 genera; all tribes) data. We found cytonuclear discordance between the two, which is likely a result of rampant hybridization among closely and more distantly related lineages. To evaluate the impact of such hybridization on the nuclear phylogeny reconstruction, we performed five different gene sampling routines, which increasingly removed putatively paralog genes. Our cleaned subset of 297 genes revealed high support for the tribes, whereas support for the main lineages (supertribes) was moderate. Calibration based on the 20 most clock-like nuclear genes suggests a late Eocene to late Oligocene origin of the family. Finally, our results strongly support a recently published new family classification, dividing the family into two subfamilies (one with five supertribes), together representing 58 tribes. This includes five recently described or re-established tribes, including Arabidopsideae, a monogeneric tribe accommodating Arabidopsis without any close relatives. With a worldwide community of thousands of researchers working on Brassicaceae and its diverse members, our new genus-level family phylogeny will be an indispensable tool for studies on biodiversity and plant biology.

The meso-octoploid Heliophila variabilis genome sheds a new light on the impact of polyploidization and diploidization on the diversity of the Cape flora

Although the South African Cape flora is one of the most remarkable biodiversity hotspots, its high diversity has not been associated with polyploidy. Here, we report the chromosome-scale genome assembly of an ephemeral cruciferous species Heliophila variabilis (~334 Mb, n = 11) adapted to South African semiarid biomes. Two pairs of differently fractionated subgenomes suggest an allo-octoploid origin of the genome at least 12 million years ago. The ancestral octoploid Heliophila genome (2n = 8x = ~60) has probably originated through hybridization between two allotetraploids (2n = 4x = ~30) formed by distant, intertribal, hybridization. Rediploidization of the ancestral genome was marked by extensive reorganization of parental subgenomes, genome downsizing, and speciation events in the genus Heliophila. We found evidence for loss-of-function changes in genes associated with leaf development and early flowering, and over-retention and sub/neofunctionalization of genes involved in pathogen response and chemical defense. The genomic resources of H. variabilis will help elucidate the role of polyploidization and genome diploidization in plant adaptation to hot arid environments and origin of the Cape flora. The sequenced H. variabilis represents the first chromosome-scale genome assembly of a meso-octoploid representative of the mustard family.

Patterns, mechanisms, and consequences of homoeologous exchange in allopolyploid angiosperms: a genomic and epigenomic perspective

Allopolyploids result from hybridization between different evolutionary lineages coupled with genome doubling. Homoeologous chromosomes (chromosomes with common shared ancestry) may undergo recombination immediately after allopolyploid formation and continue over successive generations. The outcome of this meiotic pairing behavior is dynamic and complex. Homoeologous exchanges (HEs) may lead to the formation of unbalanced gametes, reduced fertility, and selective disadvantage. By contrast, HEs could act as sources of novel evolutionary substrates, shifting the relative dosage of parental gene copies, generating novel phenotypic diversity, and helping the establishment of neo-allopolyploids. However, HE patterns vary among lineages, across generations, and even within individual genomes and chromosomes. The causes and consequences of this variation are not fully understood, though interest in this evolutionary phenomenon has increased in the last decade. Recent technological advances show promise in uncovering the mechanistic basis of HEs. Here, we describe recent observations of the common patterns among allopolyploid angiosperm lineages, underlying genomic and epigenomic features, and consequences of HEs. We identify critical research gaps and discuss future directions with far-reaching implications in understanding allopolyploid evolution and applying them to the development of important phenotypic traits of polyploid crops.

Genome assembly of the Brassicaceae diploid Orychophragmus violaceus reveals complex whole-genome duplication and evolution of dihydroxy fatty acid metabolism

Orychophragmus violaceus is a Brassicaceae species widely cultivated in China, particularly as a winter cover crop in northern China because of its low-temperature tolerance and low water demand. Recently, O. violaceus has also been cultivated as a potential industrial oilseed crop because of its abundant 24-carbon dihydroxy fatty acids (diOH-FAs), which contribute to superior high-temperature lubricant properties. In this study, we performed de novo assembly of the O. violaceus genome. Whole-genome synteny analysis of the genomes of its relatives demonstrated that O. violaceus is a diploid that has undergone an extra whole-genome duplication (WGD) after the Brassicaceae-specific α-WGD event, with a basic chromosome number of x = 12. Formation of diOH-FAs is hypothesized to have occurred after the WGD event. Based on the genome and the transcriptome data from multiple stages of seed development, we predicted that OvDGAT1-1 and OvDGAT1-2 are candidate genes for the regulation of diOH-FA storage in O. violaceus seeds. These results may greatly facilitate the development of heat-tolerant and eco-friendly plant-based lubricants using O. violaceus seed oil and improve our understanding of the genomic evolution of Brassicaceae.

The genome of Orychophragmus violaceus provides genomic insights into the evolution of Brassicaceae polyploidization and its distinct traits

Orychophragmus violaceus, referred to as "eryuelan" (February orchid) in China, is an early-flowering ornamental plant. The high oil content and abundance of unsaturated fatty acids in O. violaceus seeds make it a potential high-quality oilseed crop. Here, we generated a whole-genome assembly for O. violaceus using Nanopore and Hi-C sequencing technologies. The assembled genome of O. violaceus was ∼1.3 Gb in size, with 12 pairs of chromosomes. Through investigation of ancestral genome evolution, we determined that the genome of O. violaceus experienced a tetraploidization event from a diploid progenitor with the translocated proto-Calepineae karyotype. Comparisons between the reconstructed subgenomes of O. violaceus identified indicators of subgenome dominance, indicating that subgenomes likely originated via allotetraploidy. O. violaceus was phylogenetically close to the Brassica genus, and tetraploidy in O. violaceus occurred approximately 8.57 million years ago, close in time to the whole-genome triplication of Brassica that likely arose via an intermediate tetraploid lineage. However, the tetraploidization in Orychophragmus was independent of the hexaploidization in Brassica, as evidenced by the results from detailed phylogenetic analyses and comparisons of the break and fusion points of ancestral genomic blocks. Moreover, identification of multi-copy genes regulating the production of high-quality oil highlighted the contributions of both tetraploidization and tandem duplication to functional innovation in O. violaceus. These findings provide novel insights into the polyploidization evolution of plant species and will promote both functional genomic studies and domestication/breeding efforts in O. violaceus.

The genome sequence of Hirschfeldia incana, a new Brassicaceae model to improve photosynthetic light-use efficiency

Photosynthesis is a key process in sustaining plant and human life. Improving the photosynthetic capacity of agricultural crops is an attractive means to increase their yields. While the core mechanisms of photosynthesis are highly conserved in C₃ plants, these mechanisms are very flexible, allowing considerable diversity in photosynthetic properties. Among this diversity is the maintenance of high photosynthetic light-use efficiency at high irradiance as identified in a small number of exceptional C₃ species. Hirschfeldia incana, a member of the Brassicaceae family, is such an exceptional species, and because it is easy to grow, it is an excellent model for studying the genetic and physiological basis of this trait. Here, we present a reference genome of H. incana and confirm its high photosynthetic light-use efficiency. While H. incana has the highest photosynthetic rates found so far in the Brassicaceae, the light-saturated assimilation rates of closely related Brassica rapa and Brassica nigra are also high. The H. incana genome has extensively diversified from that of B. rapa and B. nigra through large chromosomal rearrangements, species-specific transposon activity, and differential retention of duplicated genes. Duplicated genes in H. incana, B. rapa, and B. nigra that are involved in photosynthesis and/or photoprotection show a positive correlation between copy number and gene expression, providing leads into the mechanisms underlying the high photosynthetic efficiency of these species. Our work demonstrates that the H. incana genome serves as a valuable resource for studying the evolution of high photosynthetic light-use efficiency and enhancing photosynthetic rates in crop species.

Positive selection and heat-response transcriptomes reveal adaptive features of the Brassicaceae desert model, Anastatica hierochuntica

Plant adaptation to a desert environment and its endemic heat stress is poorly understood at the molecular level. The naturally heat-tolerant Brassicaceae species Anastatica hierochuntica is an ideal extremophyte model to identify genetic adaptations that have evolved to allow plants to tolerate heat stress and thrive in deserts. We generated an A. hierochuntica reference transcriptome and identified extremophyte adaptations by comparing Arabidopsis thaliana and A. hierochuntica transcriptome responses to heat, and detecting positively selected genes in A. hierochuntica. The two species exhibit similar transcriptome adjustment in response to heat and the A. hierochuntica transcriptome does not exist in a constitutive heat 'stress-ready' state. Furthermore, the A. hierochuntica global transcriptome as well as heat-responsive orthologs, display a lower basal and higher heat-induced expression than in A. thaliana. Genes positively selected in multiple extremophytes are associated with stomatal opening, nutrient acquisition, and UV-B induced DNA repair while those unique to A. hierochuntica are consistent with its photoperiod-insensitive, early-flowering phenotype. We suggest that evolution of a flexible transcriptome confers the ability to quickly react to extreme diurnal temperature fluctuations characteristic of a desert environment while positive selection of genes involved in stress tolerance and early flowering could facilitate an opportunistic desert lifestyle.

Genome diploidization associates with cladogenesis, trait disparity, and plastid gene evolution

Angiosperm genome evolution was marked by many clade-specific whole-genome duplication events. The Microlepidieae is one of the monophyletic clades in the mustard family (Brassicaceae) formed after an ancient allotetraploidization. Postpolyploid cladogenesis has resulted in the extant c. 17 genera and 60 species endemic to Australia and New Zealand (10 species). As postpolyploid genome diploidization is a trial-and-error process under natural selection, it may proceed with different intensity and be associated with speciation events. In Microlepidieae, different extents of homoeologous recombination between the two parental subgenomes generated clades marked by slow ("cold") versus fast ("hot") genome diploidization. To gain a deeper understanding of postpolyploid genome evolution in Microlepidieae, we analyzed phylogenetic relationships in this tribe using complete chloroplast sequences, entire 35S rDNA units, and abundant repetitive sequences. The four recovered intra-tribal clades mirror the varied diploidization of Microlepidieae genomes, suggesting that the intrinsic genomic features underlying the extent of diploidization are shared among genera and species within one clade. Nevertheless, even congeneric species may exert considerable morphological disparity (e.g. in fruit shape), whereas some species within different clades experience extensive morphological convergence despite the different pace of their genome diploidization. We showed that faster genome diploidization is positively associated with mean morphological disparity and evolution of chloroplast genes (plastid-nuclear genome coevolution). Higher speciation rates in perennials than in annual species were observed. Altogether, our results confirm the potential of Microlepidieae as a promising subject for the analysis of postpolyploid genome diploidization in Brassicaceae.

Analysis of flavonol regulator evolution in the Brassicaceae reveals MYB12, MYB111 and MYB21 duplications and MYB11 and MYB24 gene loss

Background

Flavonols are the largest subgroup of flavonoids, possessing multiple functions in plants including protection against ultraviolet radiation, antimicrobial activities, and flower pigmentation together with anthocyanins. They are of agronomical and economical importance because the major off-taste component in rapeseed protein isolates is a flavonol derivative, which limits rapeseed protein use for human consumption. Flavonol production in Arabidopsis thaliana is mainly regulated by the subgroup 7 (SG7) R2R3-MYB transcription factors MYB11, MYB12, and MYB111. Recently, the SG19 MYBs MYB21, MYB24, and MYB57 were shown to regulate flavonol accumulation in pollen and stamens. The members of each subgroup are closely related, showing gene redundancy and tissue-specific expression in A. thaliana. However, the evolution of these flavonol regulators inside the Brassicaceae, especially inside the Brassiceae, which include the rapeseed crop species, is not fully understood.

Results

We studied the SG7 and SG19 MYBs in 44 species, including 31 species of the Brassicaceae, by phylogenetic analyses followed by synteny and gene expression analyses. Thereby we identified a deep MYB12 and MYB111 duplication inside the Brassicaceae, which likely occurred before the divergence of Brassiceae and Thelypodieae. These duplications of SG7 members were followed by the loss of MYB11 after the divergence of Eruca vesicaria from the remaining Brassiceae species. Similarly, MYB21 experienced duplication before the emergence of the Brassiceae tribe, where the gene loss of MYB24 is also proposed to have happened. The members of each subgroup revealed frequent overlapping spatio-temporal expression patterns in the Brassiceae member B. napus, which are assumed to compensate for the loss of MYB11 and MYB24 in the analysed tissues.

Conclusions

We identified a duplication of MYB12, MYB111, and MYB21 inside the Brassicaceae and MYB11 and MYB24 gene loss inside the tribe Brassiceae. We propose that polyploidization events have shaped the evolution of the flavonol regulators in the Brassicaceae, especially in the Brassiceae.

Idahoa and Subularia: Hidden polyploid origins of two enigmatic genera of crucifers

PREMISE

The monotypic Idahoa (I. scapigera) and the bispecific Subularia (S. aquatica and S. monticola) belong to Brassicaceae with unclear phylogenetic relationships and no tribal assignment. To fill this knowledge gap, we investigated these species and their closest relatives by combining cytogenomic and phylogenomic methods.

METHODS

We used whole plastome sequences in maximum likelihood and Bayesian inference analyses. We tested the phylogenetic informativeness of shared genomic repeats. We combined nuclear gene tree reconciliation and comparative chromosome painting (CCP) to examine the occurrence of past whole-genome duplications (WGDs).

RESULTS

The plastid data set corroborated the sister relationship between Idahoa and Subularia within the crucifer Lineage V but failed to resolve consistent topologies using both inference methods. The shared repetitive sequences provided conflicting pwhylogenetic signals. CCP analysis unexpectedly revealed that Idahoa (2n = 16) has a diploidized mesotetraploid genome, whereas two Subularia species (2n = 28 and 30) have diploidized mesoctoploid genomes. Several ancient allopolyploidy events have also been detected in closely related taxa (Chamira circaeoides, Cremolobeae, Eudemeae, and Notothlaspideae).

CONCLUSIONS

Our results suggest that the contentious phylogenetic placement of Idahoa and Subularia is best explained by two WGDs involving one or more shared parental genomes. The newly identified mesopolyploid genomes highlight the challenges of studying plant clades with complex polyploidy histories and provide a better framework for understanding genome evolution in the crucifer family.

Understanding the Role of PIN Auxin Carrier Genes under Biotic and Abiotic Stresses in Olea europaea L

The PIN-FORMED (PIN) proteins represent the most important polar auxin transporters in plants. Here, we characterized the PIN gene family in two olive genotypes, the Olea europaea subsp. europaea var. sylvestris and the var. europaea (cv. ‘Farga’). Twelve and 17 PIN genes were identified for vars. sylvestris and europaea, respectively, being distributed across 6 subfamilies. Genes encoding canonical OePINs consist of six exons, while genes encoding non-canonical OePINs are composed of five exons, with implications at protein specificities and functionality. A copia-LTR retrotransposon located in intron 4 of OePIN2b of var. europaea and the exaptation of partial sequences of that element as exons of the OePIN2b of var. sylvestris reveals such kind of event as a driving force in the olive PIN evolution. RNA-seq data showed that members from the subfamilies 1, 2, and 3 responded to abiotic and biotic stress factors. Co-expression of OePINs with genes involved in stress signaling and oxidative stress homeostasis were identified. This study highlights the importance of PIN genes on stress responses, contributing for a holistic understanding of the role of auxins in plants.

Evolutionary footprints of a cold relic in a rapidly warming world

With accelerating global warming, understanding the evolutionary dynamics of plant adaptation to environmental change is increasingly urgent. Here, we reveal the enigmatic history of the genus Cochlearia (Brassicaceae), a Pleistocene relic that originated from a drought-adapted Mediterranean sister genus during the Miocene. Cochlearia rapidly diversified and adapted to circum-Arctic regions and other cold-characterized habitat types during the Pleistocene. This sudden change in ecological preferences was accompanied by a highly complex, reticulate polyploid evolution, which was apparently triggered by the impact of repeated Pleistocene glaciation cycles. Our results illustrate that two early diversified Arctic-alpine diploid gene pools contributed differently to the evolution of this young polyploid genus now captured in a cold-adapted niche. Metabolomics revealed central carbon metabolism responses to cold in diverse species and ecotypes, likely due to continuous connections to cold habitats that may have facilitated widespread adaptation to alpine and subalpine habitats, and which we speculate were coopted from existing drought adaptations. Given the growing scientific interest in the adaptive evolution of temperature-related traits, our results provide much-needed taxonomic and phylogenomic resolution of a model system as well as first insights into the origins of its adaptation to cold.

Genome structure and evolution in the cruciferous tribe Thlaspideae (Brassicaceae)

Whole-genome duplications (WGDs) and chromosome rearrangements (CRs) play the key role in driving the diversification and evolution of plant lineages. Although the direct link between WGDs and plant diversification is well documented, relatively few studies focus on the evolutionary significance of CRs. The cruciferous tribe Thlaspideae represents an ideal model system to address the role of large-scale chromosome alterations in genome evolution, as most Thlaspideae species share the same diploid chromosome number (2n = 2x = 14). Here we constructed the genome structure in 12 Thlaspideae species, including field pennycress (Thlaspi arvense) and garlic mustard (Alliaria petiolata). We detected and precisely characterized genus- and species-specific CRs, mostly pericentric inversions, as the main genome-diversifying drivers in the tribe. We reconstructed the structure of seven chromosomes of an ancestral Thlaspideae genome, identified evolutionary stable chromosomes versus chromosomes prone to CRs, estimated the rate of CRs, and uncovered an allohexaploid origin of garlic mustard from diploid taxa closely related to A. petiolata and Parlatoria cakiloidea. Furthermore, we performed detailed bioinformatic analysis of the Thlaspideae repeatomes, and identified repetitive elements applicable as unique species- and genus-specific barcodes and chromosome landmarks. This study deepens our general understanding of the evolutionary role of CRs, particularly pericentric inversions, in plant genome diversification, and provides a robust base for follow-up whole-genome sequencing efforts.

Chromosome identification in oil palm (Elaeis guineensis) using in situ hybridization with massive pools of single copy oligonucleotides and transferability across Arecaceae species

Chromosome identification is essential for linking sequence and chromosomal maps, verifying sequence assemblies, showing structural variations and tracking inheritance or recombination of chromosomes and chromosomal segments during evolution and breeding programs. Unfortunately, identification of individual chromosomes and chromosome arms has been a major challenge for some economically important crop species with a near-continuous chromosome size range and similar morphology. Here, we developed oligonucleotide-based chromosome-specific probes that enabled us to establish a reference chromosome identification system for oil palm (Elaeis guineensis Jacq., 2n = 32). Massive oligonucleotide sequence pools were anchored to individual chromosome arms using dual and triple fluorescent in situ hybridization (EgOligoFISH). Three fluorescently tagged probe libraries were developed to contain, in total 52,506 gene-rich single-copy 47-mer oligonucleotides spanning each 0.2-0.5 Mb across strategically placed chromosome regions. They generated 19 distinct FISH signals and together with rDNA probes enabled identification of all 32 E. guineensis chromosome arms. The probes were able to identify individual homoeologous chromosome regions in the related Arecaceae palm species: American oil palm (Elaeis oleifera), date palm (Phoenix dactylifera) and coconut (Cocos nucifera) showing the comparative organization and concerted evolution of genomes in the Arecaceae. The oligonucleotide probes developed here provide a valuable approach to chromosome arm identification and allow tracking chromosome transfer in hybridization and breeding programs in oil palm, as well as comparative studies within Arecaceae.

Can functional genomic diversity provide novel insights into mechanisms of community assembly? A pilot study from an invaded alpine streambed

An important focus of community ecology, including invasion biology, is to investigate functional trait diversity patterns to disentangle the effects of environmental and biotic interactions. However, a notable limitation is that studies usually rely on a small and easy-to-measure set of functional traits, which might not immediately reflect ongoing ecological responses to changing abiotic or biotic conditions, including those that occur at a molecular or physiological level. We explored the potential of using the diversity of expressed genes-functional genomic diversity (FGD)-to understand ecological dynamics of a recent and ongoing alpine invasion. We quantified FGD based on transcriptomic data measured for 26 plant species occurring along adjacent invaded and pristine streambeds. We used an RNA-seq approach to summarize the overall number of expressed transcripts and their annotations to functional categories, and contrasted this with functional trait diversity (FTD) measured from a suite of characters that have been traditionally considered in plant ecology. We found greater FGD and FTD in the invaded community, independent of differences in species richness. However, the magnitude of functional dispersion was greater from the perspective of FGD than from FTD. Comparing FGD between congeneric alien-native species pairs, we did not find many significant differences in the proportion of genes whose annotations matched functional categories. Still, native species with a greater relative abundance in the invaded community compared with the pristine tended to express a greater fraction of genes at significant levels in the invaded community, suggesting that changes in FGD may relate to shifts in community composition. Comparisons of diversity patterns from the community to the species level offer complementary insights into processes and mechanisms driving invasion dynamics. FGD has the potential to illuminate cryptic changes in ecological diversity, and we foresee promising avenues for future extensions across taxonomic levels and macro-ecosystems.

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.

Linked by Ancestral Bonds: Multiple Whole-Genome Duplications and Reticulate Evolution in a Brassicaceae Tribe

Pervasive hybridization and whole-genome duplications (WGDs) influenced genome evolution in several eukaryotic lineages. Although frequent and recurrent hybridizations may result in reticulate phylogenies, the evolutionary events underlying these reticulations, including detailed structure of the ancestral diploid and polyploid genomes, were only rarely reconstructed. Here, we elucidate the complex genomic history of a monophyletic clade from the mustard family (Brassicaceae), showing contentious relationships to the early-diverging clades of this model plant family. Genome evolution in the crucifer tribe Biscutelleae (∼60 species, 5 genera) was dominated by pervasive hybridizations and subsequent genome duplications. Diversification of an ancestral diploid genome into several divergent but crossable genomes was followed by hybridizations between these genomes. Whereas a single genus (Megadenia) remained diploid, the four remaining genera originated by allopolyploidy (Biscutella, Lunaria, Ricotia) or autopolyploidy (Heldreichia). The contentious relationships among the Biscutelleae genera, and between the tribe and other early diverged crucifer lineages, are best explained by close genomic relatedness among the recurrently hybridizing ancestral genomes. By using complementary cytogenomics and phylogenomics approaches, we demonstrate that the origin of a monophyletic plant clade can be more complex than a parsimonious assumption of a single WGD spurring postpolyploid cladogenesis. Instead, recurrent hybridization among the same and/or closely related parental genomes may phylogenetically interlink diploid and polyploid genomes despite the incidence of multiple independent WGDs. Our results provide new insights into evolution of early-diverging Brassicaceae lineages and elucidate challenges in resolving the contentious relationships within and between land plant lineages with pervasive hybridization and WGDs.

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

PREMISE

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

A story from the Miocene: Clock-dated phylogeny of Sisymbrium L. (Sisymbrieae, Brassicaceae)

Morphological variability and imprecise generic boundaries have hindered systematic, taxonomical, and nomenclatural studies of Sisymbrium L. (Brassicaceae, Sisymbrieae DC.). The members of this almost exclusively Old-World genus grow mostly on highly porous substrates across open steppe, semidesert, or ruderal habitats in the temperate zone of the Northern Hemisphere and African subtropics. The present study placed the biological history of Sisymbrium L. into time and space and rendered the tribus Sisymbrieae as monotypic. Five nuclear-encoded and three chloroplast-encoded loci of approximately 85% of all currently accepted species were investigated. Several accessions per species covering their whole distribution range allowed for a more representative assessment of intraspecific genetic diversity. In the light of fossil absence, the impact of different secondary calibration methods and taxon sets on time spans was tested, and we showed that such a combinatorial nested dating approach is beneficial. Multigene phylogeny accompanied with a time divergence estimation analysis placed the onset and development of this tribus into the western Irano-Turanian floristic region during the Miocene. Continuous increase in continentality and decrease in temperatures promoted the diversity of the Sisymbrieae, which invaded the open grasslands habitats in Eurasia, Mediterranean, and South Africa throughout the Pliocene and Pleistocene. Our results support the assumption of the Irano-Turanian region as a biodiversity reservoir for adjacent regions.

Pilot RNA-seq data from 24 species of vascular plants at Harvard Forest

PREMISE

Large-scale projects such as the National Ecological Observatory Network (NEON) collect ecological data on entire biomes to track climate change. NEON provides an opportunity to launch community transcriptomic projects that ask integrative questions in ecology and evolution. We conducted a pilot study to investigate the challenges of collecting RNA-seq data from diverse plant communities.

METHODS

We generated >650 Gbp of RNA-seq for 24 vascular plant species representing 12 genera and nine families at the Harvard Forest NEON site. Each species was sampled twice in 2016 (July and August). We assessed transcriptome quality and content with TransRate, BUSCO, and Gene Ontology annotations.

RESULTS

Only modest differences in assembly quality were observed across multiple k-mers. On average, transcriptomes contained hits to >70% of loci in the BUSCO database. We found no significant difference in the number of assembled and annotated transcripts between diploid and polyploid transcriptomes.

DISCUSSION

We provide new RNA-seq data sets for 24 species of vascular plants in Harvard Forest. Challenges associated with this type of study included recovery of high-quality RNA from diverse species and access to NEON sites for genomic sampling. Overcoming these challenges offers opportunities for large-scale studies at the intersection of ecology and genomics.

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.

The Genome Sequence of Alpine Megacarpaea delavayi Identifies Species-Specific Whole-Genome Duplication

Megacarpaea delavayi (Brassicaceae), a plant found the high mountains of southwest China at high altitudes (3000–4800 m), is used as a vegetable or medicine. Here, we report a draft genome for this species. The assembly genome of M. delavayi is 883 Mb, and 61.59% of the genome is composed of repeat sequences. Annotation of the genome identified a total of 41,114 protein-coding genes. We found that M. delavayi experienced an independent whole-genome duplication (WGD), paralleling those independent WGDs in Iberis, Biscutella, and Anastatica in the early Miocene. Phylogenetic analyses based on the single-copy genes confirmed the position of the genus Megacarpaea within the expanded lineage II of the family and resolved its basal divergence to a subclade consisting of Anastatica, Iberis, and Biscutella. Species-specific and fast-evolving genes in M. delavayi are mainly involved in “DNA repair” and “response to UV-B radiation.” These genetic changes may together help this species survive in high-altitude environments. The reference genome reported here provides a valuable resource for studying adaptation of this and other alpine plants to the high-altitude habitats.

Phylogeny and multiple independent whole-genome duplication events in the Brassicales

PREMISE

Whole-genome duplications (WGDs) are prevalent throughout the evolutionary history of plants. For example, dozens of WGDs have been phylogenetically localized across the order Brassicales, specifically, within the family Brassicaceae. A WGD event has also been identified in the Cleomaceae, the sister family to Brassicaceae, yet its placement, as well as that of WGDs in other families in the order, remains unclear.

METHODS

Phylo-transcriptomic data were generated and used to infer a nuclear phylogeny for 74 Brassicales taxa. Genome survey sequencing was also performed on 66 of those taxa to infer a chloroplast phylogeny. These phylogenies were used to assess and confirm relationships among the major families of the Brassicales and within Brassicaceae. Multiple WGD inference methods were then used to assess the placement of WGDs on the nuclear phylogeny.

RESULTS

Well-supported chloroplast and nuclear phylogenies for the Brassicales and the putative placement of the Cleomaceae-specific WGD event Th-ɑ are presented. This work also provides evidence for previously hypothesized WGDs, including a well-supported event shared by at least two members of the Resedaceae family, and a possible event within the Capparaceae.

CONCLUSIONS

Phylogenetics and the placement of WGDs within highly polyploid lineages continues to be a major challenge. This study adds to the conversation on WGD inference difficulties by demonstrating that sampling is especially important for WGD identification and phylogenetic placement. Given its economic importance and genomic resources, the Brassicales continues to be an ideal group for assessing WGD inference methods.

Nested whole-genome duplications coincide with diversification and high morphological disparity in Brassicaceae

Angiosperms have become the dominant terrestrial plant group by diversifying for ~145 million years into a broad range of environments. During the course of evolution, numerous morphological innovations arose, often preceded by whole genome duplications (WGD). The mustard family (Brassicaceae), a successful angiosperm clade with ~4000 species, has been diversifying into many evolutionary lineages for more than 30 million years. Here we develop a species inventory, analyze morphological variation, and present a maternal, plastome-based genus-level phylogeny. We show that increased morphological disparity, despite an apparent absence of clade-specific morphological innovations, is found in tribes with WGDs or diversification rate shifts. Both are important processes in Brassicaceae, resulting in an overall high net diversification rate. Character states show frequent and independent gain and loss, and form varying combinations. Therefore, Brassicaceae pave the way to concepts of phylogenetic genome-wide association studies to analyze the evolution of morphological form and function.

The Genome of Microthlaspi erraticum (Brassicaceae) Provides Insights Into the Adaptation to Highly Calcareous Soils

Microthlaspi erraticum is widely distributed in temperate Eurasia, but restricted to Ca²⁺-rich habitats, predominantly on white Jurassic limestone, which is made up by calcium carbonate, with little other minerals. Thus, naturally occurring Microthlaspi erraticum individuals are confronted with a high concentration of Ca²⁺ ions while Mg²⁺ ion concentration is relatively low. As there is a competitive uptake between these two ions, adaptation to the soil condition can be expected. In this study, it was the aim to explore the genomic consequences of this adaptation by sequencing and analysing the genome of Microthlaspi erraticum. Its genome size is comparable with other diploid Brassicaceae, while more genes were predicted. Two Mg²⁺ transporters known to be expressed in roots were duplicated and one showed a significant degree of positive selection. It is speculated that this evolved due to the pressure to take up Mg²⁺ ions efficiently in the presence of an overwhelming amount of Ca²⁺ ions. Future studies on plants specialized on similar soils and affinity tests of the transporters are needed to provide unequivocal evidence for this hypothesis. If verified, the transporters found in this study might be useful for breeding Brassicaceae crops for higher yield on Ca²⁺-rich and Mg²⁺ -poor soils.

Genomic Blocks in Aethionema arabicum Support Arabideae as Next Diverging Clade in Brassicaceae

The tribe Aethionemeae is sister to all other crucifers, making it a crucial group for unraveling genome evolution and phylogenetic relationships within the crown group Brassicaceae. In this study, we extend the analysis of Brassicaceae genomic blocks (GBs) to Aethionema whereby we identified unique block boundaries shared only with the tribe Arabideae. This was achieved using bioinformatic methods to analyze synteny between the recently updated genome sequence of Aethionema arabicum and other high-quality Brassicaceae genome sequences. We show that compared to the largely conserved genomic structure of most non-polyploid Brassicaceae lineages, GBs are highly rearranged in Aethionema. Furthermore, we detected similarities between the genomes of Aethionema and Arabis alpina, in which also a high number of genomic rearrangements compared to those of other Brassicaceae was found. These similarities suggest that tribe Arabideae, a clade showing conflicting phylogenetic position between studies, may have diverged before diversification of the other major lineages, and highlight the potential of synteny information for phylogenetic inference.

Has the Polyploid Wave Ebbed?

There was a wave of whole genome duplications (WGD) during and subsequent to the K-Pg interface, which was followed by an increase in the proportion of species that were polyploid. I consider why this wave of polyploid speciation has continued to rise through the divergent evolution of polyploid lineages, and through rounds of homoploid and heteroploid chromosomal change. I also consider why the polyploid speciation wave is likely to rise in the next millennium. I propose that the speed of polyploid genesis through ploidal increase and through diversification among polyploids likely will be greater than the speed of diploid speciation. The increase in polyploid diversity is expected to lag well behind episodes of WGD, owing to the very long period required for species diversification either by lineage splitting or additional rounds of polyploidy, in addition to the long period of genomic adjustment to higher ploidal levels in neopolyploids.

Temporal patterns of diversification in Brassicaceae demonstrate decoupling of rate shifts and mesopolyploidization events

BACKGROUND AND AIMS

Whole-genome duplication (WGD) events are considered important driving forces of diversification. At least 11 out of 52 Brassicaceae tribes had independent mesopolyploid WGDs followed by diploidization processes. However, the association between mesopolyploidy and subsequent diversification is equivocal. Herein we show the results from a family-wide diversification analysis on Brassicaceae, and elaborate on the hypothesis that polyploidization per se is a fundamental driver in Brassicaceae evolution.

METHODS

We established a time-calibrated chronogram based on whole plastid genomes comprising representative Brassicaceae taxa and published data spanning the entire Rosidae clade. This allowed us to set multiple calibration points and anchored various Brassicaceae taxa for subsequent downstream analyses. All major splits among Brassicaceae lineages were used in BEAST analyses of 48 individually analysed tribes comprising 2101 taxa in total using the internal transcribed spacers of nuclear ribosomal DNA. Diversification patterns were investigated on these tribe-wide chronograms using BAMM and were compared with family-wide data on genome size variation and species richness.

KEY RESULTS

Brassicaceae diverged 29.9 million years ago (Mya) during the Oligocene, and the majority of tribes started diversification in the Miocene with an average crown group age of about 12.5 Mya. This matches the cooling phase right after the Mid Miocene climatic optimum. Significant rate shifts were detected in 12 out of 52 tribes during the Mio- and Pliocene, decoupled from preceding mesopolyploid WGDs. Among the various factors analysed, the combined effect of tribal crown group age and net diversification rate (speciation minus extinction) is likely to explain sufficiently species richness across Brassicaceae tribes.

CONCLUSIONS

The onset of the evolutionary splits among tribes took place under cooler and drier conditions. Pleistocene glacial cycles may have contributed to the maintenance of high diversification rates. Rate shifts are not consistently associated with mesopolyploid WGD. We propose, therefore, that WGDs in general serve as a constant 'pump' for continuous and high species diversification.

An extraordinarily stable karyotype of the woody Populus species revealed by chromosome painting

The karyotype represents the basic genetic make-up of a eukaryotic species. Comparative cytogenetic analysis of related species based on individually identified chromosomes has been conducted in only a few plant groups and not yet in woody plants. We have developed a complete set of 19 chromosome painting probes based on the reference genome of the model woody plant Populus trichocarpa. Using sequential fluorescence in situ hybridization we were able to identify all poplar chromosomes in the same metaphase cells, which led to the development of poplar karyotypes based on individually identified chromosomes. We demonstrate that five Populus species, belonging to five different sections within Populus, have maintained a remarkably conserved karyotype. No inter-chromosomal structural rearrangements were observed on any of the 19 chromosomes among the five species. Thus, the chromosomal synteny in Populus has been remarkably maintained after nearly 14 million years of divergence. We propose that the karyotypes of woody species are more stable than those of herbaceous plants since it may take a longer period of time for woody plants to fix chromosome number or structural variants in natural populations.

Evolution of Tandem Repeats Is Mirroring Post-polyploid Cladogenesis in Heliophila (Brassicaceae)

The unigeneric tribe Heliophileae encompassing more than 100 Heliophila species is morphologically the most diverse Brassicaceae lineage. The tribe is endemic to southern Africa, confined chiefly to the southwestern South Africa, home of two biodiversity hotspots (Cape Floristic Region and Succulent Karoo). The monospecific Chamira (C. circaeoides), the only crucifer species with persistent cotyledons, is traditionally retrieved as the closest relative of Heliophileae. Our transcriptome analysis revealed a whole-genome duplication (WGD) ∼26.15-29.20 million years ago, presumably preceding the Chamira/Heliophila split. The WGD was then followed by genome-wide diploidization, species radiations, and cladogenesis in Heliophila. The expanded phylogeny based on nuclear ribosomal DNA internal transcribed spacer (ITS) uncovered four major infrageneric clades (A-D) in Heliophila and corroborated the sister relationship between Chamira and Heliophila. Herein, we analyzed how the diploidization process impacted the evolution of repetitive sequences through low-coverage whole-genome sequencing of 15 Heliophila species, representing the four clades, and Chamira. Despite the firmly established infrageneric cladogenesis and different ecological life histories (four perennials vs. 11 annual species), repeatome analysis showed overall comparable evolution of genome sizes (288-484 Mb) and repeat content (25.04-38.90%) across Heliophila species and clades. Among Heliophila species, long terminal repeat (LTR) retrotransposons were the predominant components of the analyzed genomes (11.51-22.42%), whereas tandem repeats had lower abundances (1.03-12.10%). In Chamira, the tandem repeat content (17.92%, 16 diverse tandem repeats) equals the abundance of LTR retrotransposons (16.69%). Among the 108 tandem repeats identified in Heliophila, only 16 repeats were found to be shared among two or more species; no tandem repeats were shared by Chamira and Heliophila genomes. Six "relic" tandem repeats were shared between any two different Heliophila clades by a common descent. Four and six clade-specific repeats shared among clade A and C species, respectively, support the monophyly of these two clades. Three repeats shared by all clade A species corroborate the recent diversification of this clade revealed by plastome-based molecular dating. Phylogenetic analysis based on repeat sequence similarities separated the Heliophila species to three clades [A, C, and (B+D)], mirroring the post-polyploid cladogenesis in Heliophila inferred from rDNA ITS and plastome sequences.

The story of promiscuous crucifers: origin and genome evolution of an invasive species, Cardamine occulta (Brassicaceae), and its relatives

BACKGROUND AND AIMS

Cardamine occulta (Brassicaceae) is an octoploid weedy species (2n = 8x = 64) originated in Eastern Asia. It has been introduced to other continents including Europe and considered to be an invasive species. Despite its wide distribution, the polyploid origin of C. occulta remained unexplored. The feasibility of comparative chromosome painting (CCP) in crucifers allowed us to elucidate the origin and genome evolution in Cardamine species. We aimed to investigate the genome structure of C. occulta in comparison with its tetraploid (2n = 4x = 32, C. kokaiensis and C. scutata) and octoploid (2n = 8x = 64, C. dentipetala) relatives.

METHODS

Genomic in situ hybridization (GISH) and large-scale CCP were applied to uncover the parental genomes and chromosome composition of the investigated Cardamine species.

KEY RESULTS

All investigated species descended from a common ancestral Cardamine genome (n = 8), structurally resembling the Ancestral Crucifer Karyotype (n = 8), but differentiated by a translocation between chromosomes AK6 and AK8. Allotetraploid C. scutata originated by hybridization between two diploid species, C. parviflora and C. amara (2n = 2x = 16). By contrast, C. kokaiensis has an autotetraploid origin from a parental genome related to C. parviflora. Interestingly, octoploid C. occulta probably originated through hybridization between the tetraploids C. scutata and C. kokaiensis. The octoploid genome of C. dentipetala probably originated from C. scutata via autopolyploidization. Except for five species-specific centromere repositionings and one pericentric inversion post-dating the polyploidization events, the parental subgenomes remained stable in the tetra- and octoploids.

CONCLUSIONS

Comparative genome structure, origin and evolutionary history was reconstructed in C. occulta and related species. For the first time, whole-genome cytogenomic maps were established for octoploid plants. Post-polyploid evolution in Asian Cardamine polyploids has not been associated with descending dysploidy and intergenomic rearrangements. The combination of different parental (sub)genomes adapted to distinct habitats provides an evolutionary advantage to newly formed polyploids by occupying new ecological niches.

Interspecies association mapping links reduced CG to TG substitution rates to the loss of gene-body methylation

Comparative genomics can unravel the genetic basis of species differences; however, successful reports on quantitative traits are still scarce. Here we present genome assemblies of 31 so-far unassembled Brassicaceae plant species and combine them with 16 previously published assemblies to establish the Brassicaceae Diversity Panel. Using a new interspecies association strategy for quantitative traits, we found a so-far unknown association between the unexpectedly high variation in CG to TG substitution rates in genes and the absence of CHROMOMETHYLASE3 (CMT3) orthologues. Low substitution rates were associated with the loss of CMT3, while species with conserved CMT3 orthologues showed high substitution rates. Species without CMT3 also lacked gene-body methylation (gbM), suggesting an evolutionary trade-off between the unknown function of gbM and low substitution rates in Brassicaceae, possibly due to low mutability of non-methylated cytosines.

A Target Capture-Based Method to Estimate Ploidy From Herbarium Specimens

Whole genome duplication (WGD) events are common in many plant lineages, but the ploidy status and possible occurrence of intraspecific ploidy variation are unknown for most species. Standard methods for ploidy determination are chromosome counting and flow cytometry approaches. While flow cytometry approaches typically use fresh tissue, an increasing number of studies have shown that recently dried specimens can be used to yield ploidy data. Recent studies have started to explore whether high-throughput sequencing (HTS) data can be used to assess ploidy levels by analyzing allelic frequencies from single copy nuclear genes. Here, we compare different approaches using a range of yam (Dioscorea) tissues of varying ages, drying methods and quality, including herbarium tissue. Our aims were to: (1) explore the limits of flow cytometry in estimating ploidy level from dried samples, including herbarium vouchers collected between 1831 and 2011, and (2) optimize a HTS-based method to estimate ploidy by considering allelic frequencies from nuclear genes obtained using a target-capture method. We show that, although flow cytometry can be used to estimate ploidy levels from herbarium specimens collected up to fifteen years ago, success rate is low (5.9%). We validated our HTS-based estimates of ploidy using 260 genes by benchmarking with dried samples of species of known ploidy (Dioscorea alata, D. communis, and D. sylvatica). Subsequently, we successfully applied the method to the 85 herbarium samples analyzed with flow cytometry, and successfully provided results for 91.7% of them, comprising species across the phylogenetic tree of Dioscorea. We also explored the limits of using this HTS-based approach for identifying high ploidy levels in herbarium material and the effects of heterozygosity and sequence coverage. Overall, we demonstrated that ploidy diversity within and between species may be ascertained from historical collections, allowing the determination of polyploidization events from samples collected up to two centuries ago. This approach has the potential to provide insights into the drivers and dynamics of ploidy level changes during plant evolution and crop domestication.

The repetitive DNA landscape in Avena (Poaceae): chromosome and genome evolution defined by major repeat classes in whole-genome sequence reads

BACKGROUND

Repetitive DNA motifs - not coding genetic information and repeated millions to hundreds of times - make up the majority of many genomes. Here, we identify the nature, abundance and organization of all the repetitive DNA families in oats (Avena sativa, 2n = 6x = 42, AACCDD), a recognized health-food, and its wild relatives.

RESULTS

Whole-genome sequencing followed by k-mer and RepeatExplorer graph-based clustering analyses enabled assessment of repetitive DNA composition in common oat and its wild relatives' genomes. Fluorescence in situ hybridization (FISH)-based karyotypes are developed to understand chromosome and repetitive sequence evolution of common oat. We show that some 200 repeated DNA motifs make up 70% of the Avena genome, with less than 20 families making up 20% of the total. Retroelements represent the major component, with Ty3/Gypsy elements representing more than 40% of all the DNA, nearly three times more abundant than Ty1/Copia elements. DNA transposons are about 5% of the total, while tandemly repeated, satellite DNA sequences fit into 55 families and represent about 2% of the genome. The Avena species are monophyletic, but both bioinformatic comparisons of repeats in the different genomes, and in situ hybridization to metaphase chromosomes from the hexaploid species, shows that some repeat families are specific to individual genomes, or the A and D genomes together. Notably, there are terminal regions of many chromosomes showing different repeat families from the rest of the chromosome, suggesting presence of translocations between the genomes.

CONCLUSIONS

The relatively small number of repeat families shows there are evolutionary constraints on their nature and amplification, with mechanisms leading to homogenization, while repeat characterization is useful in providing genome markers and to assist with future assemblies of this large genome (c. 4100 Mb in the diploid). The frequency of inter-genomic translocations suggests optimum strategies to exploit genetic variation from diploid oats for improvement of the hexaploid may differ from those used widely in bread wheat.

Brassicaceae flowers: diversity amid uniformity

The mustard family Brassicaceae, which includes the model plant Arabidopsis thaliana, exhibits morphological stasis and significant uniformity of floral plan. Nonetheless, there is untapped diversity in almost every aspect of floral morphology in the family that lends itself to comparative study, including organ number, shape, form, and color. Studies on the genetic basis of morphological diversity, enabled by extensive genetic tools and genomic resources and the close phylogenetic distance among mustards, have revealed a mosaic of conservation and divergence in numerous floral traits. Here I review the morphological diversity of the flowers of Brassicaceae and discuss studies addressing the underlying genetic and developmental mechanisms shaping floral diversity. To put flowers in the context of the floral display, I describe diversity in inflorescence morphology and the variation that exists in the structures preceding the floral organs. Reconstructing the floral morphospace in Brassicaceae coupled with next-generation sequencing data and unbiased approaches to interrogate gene function in species throughout the mustard phylogeny offers promising ways to understand how developmental mechanisms originate and diversify.

Resolving the backbone of the Brassicaceae phylogeny for investigating trait diversity

The Brassicaceae family comprises c. 4000 species including economically important crops and the model plant Arabidopsis thaliana. Despite their importance, the relationships among major lineages in the family remain unresolved, hampering comparative research. Here, we inferred a Brassicaceae phylogeny using newly generated targeted enrichment sequence data of 1827 exons (> 940 000 bases) representing 63 species, as well as sequenced genome data of 16 species, together representing 50 of the 52 currently recognized Brassicaceae tribes. A third of the samples were derived from herbarium material, facilitating broad taxonomic coverage of the family. Six major clades formed successive sister groups to the rest of Brassicaceae. We also recovered strong support for novel relationships among tribes, and resolved the position of 16 taxa previously not assigned to a tribe. The broad utility of these phylogenetic results is illustrated through a comparative investigation of genome-wide expression signatures that distinguish simple from complex leaves in Brassicaceae. Our study provides an easily extendable dataset for further advances in Brassicaceae systematics and a timely higher-level phylogenetic framework for a wide range of comparative studies of multiple traits in an intensively investigated group of plants.

Hybridization-facilitated genome merger and repeated chromosome fusion after 8 million years

The small genus Ricotia (nine species, Brassicaceae) is confined to the eastern Mediterranean. By comparative chromosome painting and a dated multi-gene chloroplast phylogeny, we reconstructed the origin and subsequent evolution of Ricotia. The ancestral Ricotia genome originated through hybridization between two older genomes with n = 7 and n = 8 chromosomes, respectively, on the Turkish mainland during the Early Miocene (c. 17.8 million years ago, Ma). Since then, the allotetraploid (n = 15) genome has been altered by two independent descending dysploidies (DD) to n = 14 in Ricotia aucheri and the Tenuifolia clade (2 spp.). By the Late Miocene (c. 10 Ma), the latter clade started to evolve in the most diverse Ricotia core clade (6 spp.), the process preceded by a DD event to n = 13. It is noteworthy that this dysploidy was mediated by a unique chromosomal rearrangement, merging together the same two chromosomes as were merged during the origin of a fusion chromosome within the paternal n = 7 genome c. 20 Ma. This shows that within a time period of c. 8 Myr genome evolution can repeat itself and that structurally very similar chromosomes may originate repeatedly from the same ancestral chromosomes by different pathways (end-to-end translocation versus nested chromosome insertion).

Multiple large-scale gene and genome duplications during the evolution of hexapods

Polyploidy or whole genome duplication (WGD) is a major contributor to genome evolution and diversity. Although polyploidy is recognized as an important component of plant evolution, it is generally considered to play a relatively minor role in animal evolution. Ancient polyploidy is found in the ancestry of some animals, especially fishes, but there is little evidence for ancient WGDs in other metazoan lineages. Here we use recently published transcriptomes and genomes from more than 150 species across the insect phylogeny to investigate whether ancient WGDs occurred during the evolution of Hexapoda, the most diverse clade of animals. Using gene age distributions and phylogenomics, we found evidence for 18 ancient WGDs and six other large-scale bursts of gene duplication during insect evolution. These bursts of gene duplication occurred in the history of lineages such as the Lepidoptera, Trichoptera, and Odonata. To further corroborate the nature of these duplications, we evaluated the pattern of gene retention from putative WGDs observed in the gene age distributions. We found a relatively strong signal of convergent gene retention across many of the putative insect WGDs. Considering the phylogenetic breadth and depth of the insect phylogeny, this observation is consistent with polyploidy as we expect dosage balance to drive the parallel retention of genes. Together with recent research on plant evolution, our hexapod results suggest that genome duplications contributed to the evolution of two of the most diverse lineages of eukaryotes on Earth.

Post-polyploid diploidization and diversification through dysploid changes

Whole-genome duplications are widespread across land plant phylogenies and particularly frequent in ferns and angiosperms. Genome duplications spurred the evolution of key innovations associated with diversification in many angiosperm clades and lineages. Such diversifications are not initiated by genome doubling per se. Rather, differentiation of the primary polyploid populations through a range of processes results in post-polyploid genome diploidization. Structural diploidization gradually reverts the polyploid genome to one functionally diploid-like through chromosomal rearrangements which frequently result in dysploid changes. Dysploidies may lead to reproductive isolation among post-polyploid offspring and significantly contribute to speciation and cladogenetic events.

The long and short of doubling down: polyploidy, epigenetics, and the temporal dynamics of genome fractionation

We consider the rapidly advancing discipline of plant evolutionary genomics, with a focus on the evolution of polyploid genomes. In many lineages, polyploidy is followed by 'biased fractionation', the unequal loss of genes from ancestral progenitor genomes. Mechanistically, it has been proposed that biased fractionation results from changes in the epigenetic landscape near genes, likely mediated by transposable elements. These epigenetic changes result in unequal gene expression between duplicates, establishing differential fitness that leads to biased gene loss with respect to ancestral genomes. We propose a unifying conceptual framework and a set of testable hypotheses based on this model, relating genome size, the proximity of transposable elements to genes, epigenetic reprogramming, chromatin accessibility, and gene expression.

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.

Improved transcriptome sampling pinpoints 26 ancient and more recent polyploidy events in Caryophyllales, including two allopolyploidy events

Studies of the macroevolutionary legacy of polyploidy are limited by an incomplete sampling of these events across the tree of life. To better locate and understand these events, we need comprehensive taxonomic sampling as well as homology inference methods that accurately reconstruct the frequency and location of gene duplications. We assembled a data set of transcriptomes and genomes from 168 species in Caryophyllales, of which 43 transcriptomes were newly generated for this study, representing one of the most densely sampled genomic-scale data sets available. We carried out phylogenomic analyses using a modified phylome strategy to reconstruct the species tree. We mapped the phylogenetic distribution of polyploidy events by both tree-based and distance-based methods, and explicitly tested scenarios for allopolyploidy. We identified 26 ancient and more recent polyploidy events distributed throughout Caryophyllales. Two of these events were inferred to be allopolyploidy. Through dense phylogenomic sampling, we show the propensity of polyploidy throughout the evolutionary history of Caryophyllales. We also provide a framework for utilizing transcriptome data to detect allopolyploidy, which is important as it may have different macroevolutionary implications compared with autopolyploidy.