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

ATAD3 Proteins: Unique Mitochondrial Proteins Essential for Life in Diverse Eukaryotic Lineages

ATPase family AAA domain-containing 3 (ATAD3) proteins are unique mitochondrial proteins that arose deep in the eukaryotic lineage but that are surprisingly absent in Fungi and Amoebozoa. These ∼600-amino acid proteins are anchored in the inner mitochondrial membrane and are essential in metazoans and Arabidopsis thaliana. ATAD3s comprise a C-terminal ATPases Associated with a variety of cellular Activities (AAA+) matrix domain and an ATAD3_N domain, which is located primarily in the inner membrane space but potentially extends to the cytosol to interact with the ER. Sequence and structural alignments indicate that ATAD3 proteins are most similar to classic chaperone unfoldases in the AAA+ family, suggesting that they operate in mitochondrial protein quality control. A. thaliana has four ATAD3 genes in two distinct clades that appear first in the seed plants, and both clades are essential for viability. The four genes are generally coordinately expressed, and transcripts are highest in growing apices and imbibed seeds. Plants with disrupted ATAD3 have reduced growth, aberrant mitochondrial morphology, diffuse nucleoids and reduced oxidative phosphorylation complex I. These and other pleiotropic phenotypes are also observed in ATAD3 mutants in metazoans. Here, we discuss the distribution of ATAD3 proteins as they have evolved in the plant kingdom, their unique structure, what we know about their function in plants and the challenges in determining their essential roles in mitochondria.

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.

Differential expression and evolutionary diversification of RNA helicases in Boechera sexual and apomictic reproduction

In higher plants, sexual reproduction is characterized by meiosis of the first cells of the germlines, and double fertilization of the egg and central cell after gametogenesis. In contrast, in apomicts of the genus Boechera, meiosis is omitted or altered and only the central cell requires fertilization, while the embryo forms parthenogenetically from the egg cell. To deepen the understanding of the transcriptional basis underlying these differences, we applied RNA-seq to compare expression in reproductive tissues of different Boechera accessions. This confirmed previous evidence of an enrichment of RNA helicases in plant germlines. Furthermore, few RNA helicases were differentially expressed in female reproductive ovule tissues harboring mature gametophytes from apomictic and sexual accessions. For some of these genes, we further found evidence for a complex recent evolutionary history. This included a homolog of Arabidopsis thaliana FASCIATED STEM4 (FAS4). In contrast to AtFAS4, which is a single-copy gene, FAS4 is represented by three homologs in Boechera, suggesting a potential for subfunctionalization to modulate reproductive development. To gain first insights into functional roles of FAS4, we studied Arabidopsis lines carrying mutant alleles. This identified the crucial importance of AtFAS4 for reproduction, as we observed developmental defects and arrest during male and female gametogenesis.

Halophytes and heavy metals: A multi-omics approach to understand the role of gene and genome duplication in the abiotic stress tolerance of Cakile maritima

PREMISE

The origin of diversity is a fundamental biological question. Gene duplications are one mechanism that provides raw material for the emergence of novel traits, but evolutionary outcomes depend on which genes are retained and how they become functionalized. Yet, following different duplication types (polyploidy and tandem duplication), the events driving gene retention and functionalization remain poorly understood. Here we used Cakile maritima, a species that is tolerant to salt and heavy metals and shares an ancient whole-genome triplication with closely related salt-sensitive mustard crops (Brassica), as a model to explore the evolution of abiotic stress tolerance following polyploidy.

METHODS

Using a combination of ionomics, free amino acid profiling, and comparative genomics, we characterize aspects of salt stress response in C. maritima and identify retained duplicate genes that have likely enabled adaptation to salt and mild levels of cadmium.

RESULTS

Cakile maritima is tolerant to both cadmium and salt treatments through uptake of cadmium in the roots. Proline constitutes greater than 30% of the free amino acid pool in C. maritima and likely contributes to abiotic stress tolerance. We find duplicated gene families are enriched in metabolic and transport processes and identify key transport genes that may be involved in C. maritima abiotic stress tolerance.

CONCLUSIONS

These findings identify pathways and genes that could be used to enhance plant resilience and provide a putative understanding of the roles of duplication types and retention on the evolution of abiotic stress response.

Genomes of Meniocus linifolius and Tetracme quadricornis reveal the ancestral karyotype and genomic features of core Brassicaceae

Brassicaceae represents an important plant family from both a scientific and economic perspective. However, genomic features related to the early diversification of this family have not been fully characterized, especially upon the uplift of the Tibetan Plateau, which was followed by increasing aridity in the Asian interior, intensifying monsoons in Eastern Asia, and significantly fluctuating daily temperatures. Here, we reveal the genomic architecture that accompanied early Brassicaceae diversification by analyzing two high-quality chromosome-level genomes for Meniocus linifolius (Arabodae; clade D) and Tetracme quadricornis (Hesperodae; clade E), together with genomes representing all major Brassicaceae clades and the basal Aethionemeae. We reconstructed an ancestral core Brassicaceae karyotype (CBK) containing 9 pseudochromosomes with 65 conserved syntenic genomic blocks and identified 9702 conserved genes in Brassicaceae. We detected pervasive conflicting phylogenomic signals accompanied by widespread ancient hybridization events, which correlate well with the early divergence of core Brassicaceae. We identified a successive Brassicaceae-specific expansion of the class I TREHALOSE-6-PHOSPHATE SYNTHASE 1 (TPS1) gene family, which encodes enzymes with essential regulatory roles in flowering time and embryo development. The TPS1s were mainly randomly amplified, followed by expression divergence. Our results provide fresh insights into historical genomic features coupled with Brassicaceae evolution and offer a potential model for broad-scale studies of adaptive radiation under an ever-changing environment.

Transposable elements contribute to the establishment of the glycine shuttle in Brassicaceae species

C3 -C4 intermediate photosynthesis has evolved at least five times convergently in the Brassicaceae, despite this family lacking bona fide C4 species. The establishment of this carbon concentrating mechanism is known to require a complex suite of ultrastructural modifications, as well as changes in spatial expression patterns, which are both thought to be underpinned by a reconfiguration of existing gene-regulatory networks. However, to date, the mechanisms which underpin the reconfiguration of these gene networks are largely unknown. In this study, we used a pan-genomic association approach to identify genomic features that could confer differential gene expression towards the C3 -C4 intermediate state by analysing eight C3 species and seven C3 -C4 species from five independent origins in the Brassicaceae. We found a strong correlation between transposable element (TE) insertions in cis-regulatory regions and C3 -C4 intermediacy. Specifically, our study revealed 113 gene models in which the presence of a TE within a gene correlates with C3 -C4 intermediate photosynthesis. In this set, genes involved in the photorespiratory glycine shuttle are enriched, including the glycine decarboxylase P-protein whose expression domain undergoes a spatial shift during the transition to C3 -C4 photosynthesis. When further interrogating this gene, we discovered independent TE insertions in its upstream region which we conclude to be responsible for causing the spatial shift in GLDP1 gene expression. Our findings hint at a pivotal role of TEs in the evolution of C3 -C4 intermediacy, especially in mediating differential spatial gene expression.

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.

Angiosperm-wide analysis of fruit and ovary evolution aided by a new nuclear phylogeny supports association of the same ovary type with both dry and fleshy fruits

Fruit functions in seed protection and dispersal and belongs to many dry and fleshy types, yet their evolutionary pattern remains unclear in part due to uncertainties in the phylogenetic relationships among several orders and families. Thus we used nuclear genes of 502 angiosperm species representing 231 families to reconstruct a well supported phylogeny, with resolved relationships for orders and families with previously uncertain placements. Using this phylogeny as a framework, molecular dating supports a Triassic origin of the crown angiosperms, followed by the emergence of most orders in the Jurassic and Cretaceous and their rise to ecological dominance during the Cretaceous Terrestrial Revolution. The robust phylogeny allowed an examination of the evolutionary pattern of fruit and ovary types, revealing a trend of parallel carpel fusions during early diversifications in eudicots, monocots, and magnoliids. Moreover, taxa in the same order or family with the same ovary type can develop either dry or fleshy fruits with strong correlations between specific types of dry and fleshy fruits; such associations of ovary, dry and fleshy fruits define several ovary-fruit "modules" each found in multiple families. One of the frequent modules has an ovary containing multiple ovules, capsules and berries, and another with an ovary having one or two ovules, achenes (or other single-seeded dry fruits) and drupes. This new perspective of relationships among fruit types highlights the closeness of specific dry and fleshy fruit types, such as capsule and berry, that develop from the same ovary type and belong to the same module relative to dry and fleshy fruits of other modules (such as achenes and drupes). Further analyses of gene families containing known genes for ovary and fruit development identified phylogenetic nodes with multiple gene duplications, supporting a possible role of whole-genome duplications, in combination with climate changes and animal behaviors, in angiosperm fruit and ovary diversification.

Genome-wide identification, evolution, and expression analysis of the NAC gene family in chestnut (Castanea mollissima)

The NAC gene family is one of the most important transcription factor families specific to plants, responsible for regulating many biological processes, including development, stress response, and signal transduction. However, it has not yet been characterized in chestnut, an important nut tree species. Here, we identified 115 CmNAC genes in the chestnut genome, which were divided into 16 subgroups based on the phylogenetic analysis. Numerous cis-acting elements related to auxin, gibberellin, and abscisic acid were identified in the promoter region of CmNACs, suggesting that they play an important role in the growth and development of chestnut. The results of the collinear analysis indicated that dispersed duplication and whole-genome-duplication were the main drivers of CmNAC gene expansion. RNA-seq data of developmental stages of chestnut nut, bud, and ovule revealed the expression patterns of CmNAC genes. Additionally, qRT-PCR experiments were used to verify the expression levels of some CmNAC genes. The comprehensive analysis of the above results revealed that some CmNAC members may be related to chestnut bud and nut development, as well as ovule fertility. The systematic analysis of this study will help to increase understanding of the potential functions of the CmNAC genes in chestnut growth and development.

Characterization of the complete plastome sequence of Korean endemic, Cardamine glechomifolia H.Lév. (Brassicaceae, Brassicales)

In this study, we report the complete plastome sequence of Cardamine glechomifolia H.Lév. 1913 (NCBI acc. no. OP894664). This plastome shows typical quadripartite structure. The plastome size is 154,307 bp, which consists of 84,015 bp large single-copy (LSC), 17,690 bp small single-copy (SSC), and 26,301 bp inverted repeat (IR) regions. The plastome contains 112 genes, including 78 protein-coding, 30 tRNA, and four rRNA genes. The infA gene is pseudogenized. Sixteen genes contain one intron and two genes (clpP and ycf3) have two introns. The phylogenomic analysis conducted in our study reveals that the genus Cardamine, which encompasses C. glechomifolia, exhibits three distinct clades. In order to elucidate the interrelationship among the three clades, it is imperative to conduct additional investigations by augmenting the number of Cardamine samples.

Remarkable mitochondrial genome heterogeneity in Meniocus linifolius (Brassicaceae)

KEY MESSAGE

Detailed analyses of 16 genomes identified a remarkable acceleration of mutation rate, hence mitochondrial sequence and structural heterogeneity, in Meniocus linifolius (Brassicaceae). The powerhouse, mitochondria, in plants feature high levels of structural variation, while the encoded genes are normally conserved. However, the substitution rates and spectra of mitochondria DNA within the Brassicaceae, a family with substantial scientific and economic importance, have not been adequately deciphered. Here, by analyzing three newly assembled and 13 known mitochondrial genomes (mitogenomes), we report the highly variable genome structure and mutation rates in Brassicaceae. The genome sizes and GC contents are 196,604 bp and 46.83%, 288,122 bp and 44.79%, and 287,054 bp and 44.93%, for Meniocus linifolius (Mli), Crucihimalaya lasiocarpa (Cla), and Lepidium sativum (Lsa), respectively. In total, 29, 33, and 34 protein-coding genes (PCGs) and 14, 18, and 18 tRNAs are annotated for Mli, Cla, and Lsa, respectively, while all mitogenomes contain one complete circular molecule with three rRNAs and abundant RNA editing sites. The Mli mitogenome features four conformations likely mediated by the two pairs of long repeats, while at the same time seems to have an unusual evolutionary history due to higher GC content, loss of more genes and sequences, but having more repeats and plastid DNA insertions. Corroborating with these, an ambiguous phylogenetic position with long branch length and elevated synonymous substitution rate in nearly all PCGs are observed for Mli. Taken together, our results reveal a high level of mitogenome heterogeneity at the family level and provide valuable resources for further understanding the evolutionary pattern of organelle genomes in Brassicaceae.

Intergrative metabolomic and transcriptomic analyses reveal the potential regulatory mechanism of unique dihydroxy fatty acid biosynthesis in the seeds of an industrial oilseed crop Orychophragmus violaceus

BACKGROUND

Orychophragmus violaceus is a potentially important industrial oilseed crop due to the two 24-carbon dihydroxy fatty acids (diOH-FA) that was newly identified from its seed oil via a 'discontinuous elongation' process. Although many research efforts have focused on the diOH-FA biosynthesis mechanism and identified the potential co-expressed diacylglycerol acyltranferase (DGAT) gene associated with triacylglycerol (TAG)-polyestolides biosynthesis, the dynamics of metabolic changes during seed development of O. violaceus as well as its associated regulatory network changes are poorly understood.

RESULTS

In this study, by combining metabolome and transcriptome analysis, we identified that 1,003 metabolites and 22,479 genes were active across four stages of seed development, which were further divided into three main clusters based on the patterns of metabolite accumulation and/or gene expression. Among which, cluster2 was mostly related to diOH-FA biosynthesis pathway. We thus further constructed transcription factor (TF)-structural genes regulatory map for the genes associated with the flavonoids, fatty acids and diOH-FA biosynthesis pathway in this cluster. In particular, several TF families such as bHLH, B3, HD-ZIP, MYB were found to potentially regulate the metabolism associated with the diOH-FA pathway. Among which, multiple candidate TFs with promising potential for increasing the diOH-FA content were identified, and we further traced the evolutionary history of these key genes among species of Brassicaceae.

CONCLUSION

Taken together, our study provides new insight into the gene resources and potential relevant regulatory mechanisms of diOH-FA biosynthesis uniquely in seeds of O. violaceus, which will help to promote the downstream breeding efforts of this potential oilseed crop and advance the bio-lubricant industry.

Phylogenomics insights into gene evolution, rapid species diversification, and morphological innovation of the apple tribe (Maleae, Rosaceae)

Maleae is one of the most widespread tribes of Rosaceae and includes several important fruit crops and ornamental plants. We used nuclear genes from 62 transcriptomes/genomes, including 26 newly generated transcriptomes, to reconstruct a well-supported phylogeny and study the evolution of fruit and leaf morphology and the possible effect of whole genome duplication (WGD). Our phylogeny recovered 11 well-supported clades and supported the monophyly of most genera (except Malus, Sorbus, and Pourthiaea) with at least two sampled species. A WGD was located to the most recent common ancestor (MRCA) of Maleae and dated to c. 54 million years ago (Ma) near the Early Eocene Climatic Optimum, supporting Gillenieae (x = 9) being a parental lineage of Maleae (x = 17) and including duplicate regulatory genes related to the origin of the fleshy pome fruit. Whole genome duplication-derived paralogs that are retained in specific lineages but lost in others are predicted to function in development, metabolism, and other processes. An upshift of diversification and innovations of fruit and leaf morphologies occurred at the MRCA of the Malinae subtribe, coinciding with the Eocene-Oligocene transition (c. 34 Ma), following a lag from the time of the WGD event. Our results provide new insights into the Maleae phylogeny, its rapid diversification, and morphological and molecular evolution.

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.

Evolution of phenotypic disparity in the plant kingdom

The plant kingdom exhibits diverse bodyplans, from single-celled algae to complex multicellular land plants, but it is unclear how this phenotypic disparity was achieved. Here we show that the living divisions comprise discrete clusters within morphospace, separated largely by reproductive innovations, the extinction of evolutionary intermediates and lineage-specific evolution. Phenotypic complexity correlates not with disparity but with ploidy history, reflecting the role of genome duplication in plant macroevolution. Overall, the plant kingdom exhibits a pattern of episodically increasing disparity throughout its evolutionary history that mirrors the evolutionary floras and reflects ecological expansion facilitated by reproductive innovations. This pattern also parallels that seen in the animal and fungal kingdoms, suggesting a general pattern for the evolution of multicellular bodyplans.

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.

The Impact of Whole Genome Duplication on the Evolution of the Arachnids

The proliferation of genomic resources for Chelicerata in the past 10 years has revealed that the evolution of chelicerate genomes is more dynamic than previously thought, with multiple waves of ancient whole genome duplications affecting separate lineages. Such duplication events are fascinating from the perspective of evolutionary history because the burst of new gene copies associated with genome duplications facilitates the acquisition of new gene functions (neofunctionalization), which may in turn lead to morphological novelties and spur net diversification. While neofunctionalization has been invoked in several contexts with respect to the success and diversity of spiders, the overall impact of whole genome duplications on chelicerate evolution and development remains imperfectly understood. The purpose of this review is to examine critically the role of whole genome duplication on the diversification of the extant arachnid orders, as well as assess functional datasets for evidence of subfunctionalization or neofunctionalization in chelicerates. This examination focuses on functional data from two focal model taxa: the spider Parasteatoda tepidariorum, which exhibits evidence for an ancient duplication, and the harvestman Phalangium opilio, which exhibits an unduplicated genome. I show that there is no evidence that taxa with genome duplications are more successful than taxa with unduplicated genomes. I contend that evidence for sub- or neofunctionalization of duplicated developmental patterning genes in spiders is indirect or fragmentary at present, despite the appeal of this postulate for explaining the success of groups like spiders. Available expression data suggest that the condition of duplicated Hox modules may have played a role in promoting body plan disparity in the posterior tagma of some orders, such as spiders and scorpions, but functional data substantiating this postulate are critically missing. Spatiotemporal dynamics of duplicated transcription factors in spiders may represent cases of developmental system drift, rather than neofunctionalization. Developmental system drift may represent an important, but overlooked, null hypothesis for studies of paralogs in chelicerate developmental biology. To distinguish between subfunctionalization, neofunctionalization, and developmental system drift, concomitant establishment of comparative functional datasets from taxa exhibiting the genome duplication, as well as those that lack the paralogy, is sorely needed.

Evolutionary Dynamics of FLC-like MADS-Box Genes in Brassicaceae

MADS-box genes encode transcription factors that play important roles in the development and evolution of plants. There are more than a dozen clades of MADS-box genes in angiosperms, of which those with functions in the specification of floral organ identity are especially well-known. From what has been elucidated in the model plant Arabidopsis thaliana, the clade of FLC-like MADS-box genes, comprising FLC-like genes sensu strictu and MAF-like genes, are somewhat special among the MADS-box genes of plants since FLC-like genes, especially MAF-like genes, show unusual evolutionary dynamics, in that they generate clusters of tandemly duplicated genes. Here, we make use of the latest genomic data of Brassicaceae to study this remarkable feature of the FLC-like genes in a phylogenetic context. We have identified all FLC-like genes in the genomes of 29 species of Brassicaceae and reconstructed the phylogeny of these genes employing a Maximum Likelihood method. In addition, we conducted selection analyses using PAML. Our results reveal that there are three major clades of FLC-like genes in Brassicaceae that all evolve under purifying selection but with remarkably different strengths. We confirm that the tandem arrangement of MAF-like genes in the genomes of Brassicaceae resulted in a high rate of duplications and losses. Interestingly, MAF-like genes also seem to be prone to transposition. Considering the role of FLC-like genes sensu lato (s.l.) in the timing of floral transition, we hypothesize that this rapid evolution of the MAF-like genes was a main contributor to the successful adaptation of Brassicaceae to different environments.

Underestimated diversity in high elevations of a global biodiversity hotspot: two new endemic species of Aethionema (Brassicaceae) from the alpine zone of Iran

Although the mountains in South-West Asia are a global biodiversity hotspot, our understanding of their biodiversity, especially in the commonly remote alpine and subnival zones, is still limited. This is well exemplified here by Aethionema umbellatum (Brassicaceae), a species considered to have a wide yet disjoint distribution in the Zagros and Yazd-Kerman mountains of western and central Iran. Morphological and molecular phylogenetic data (based on plastid trnL-trnF and nuclear ITS sequences) show that A. umbellatum is restricted to a single mountain range in southwestern Iran (Dena Mts., southern Zagros), whereas populations from central Iran (Yazd-Kerman and central Zagros) and from western Iran (central Zagros) belong to species new to science, A. alpinum and A. zagricum, respectively. Both new species are phylogenetically and morphologically close to A. umbellatum, with which they share unilocular fruits and one-seeded locules. However, they are easily distinguishable by leaf shape, petal size, and fruit characters. This study confirms that the alpine flora of the Irano-Anatolian region is still poorly known. As the proportion of rare and local endemic species in alpine habitats is high, these habitats are of prime interest for conservation efforts.

Completing the TRB family: newly characterized members show ancient evolutionary origins and distinct localization, yet similar interactions

Telomere repeat binding proteins (TRBs) belong to a family of proteins possessing a Myb-like domain which binds to telomeric repeats. Three members of this family (TRB1, TRB2, TRB3) from Arabidopsis thaliana have already been described as associated with terminal telomeric repeats (telomeres) or short interstitial telomeric repeats in gene promoters (telo-boxes). They are also known to interact with several protein complexes: telomerase, Polycomb repressive complex 2 (PRC2) E(z) subunits and the PEAT complex (PWOs-EPCRs-ARIDs-TRBs). Here we characterize two novel members of the TRB family (TRB4 and TRB5). Our wide phylogenetic analyses have shown that TRB proteins evolved in the plant kingdom after the transition to a terrestrial habitat in Streptophyta, and consequently TRBs diversified in seed plants. TRB4-5 share common TRB motifs while differing in several others and seem to have an earlier phylogenetic origin than TRB1-3. Their common Myb-like domains bind long arrays of telomeric repeats in vitro, and we have determined the minimal recognition motif of all TRBs as one telo-box. Our data indicate that despite the distinct localization patterns of TRB1-3 and TRB4-5 in situ, all members of TRB family mutually interact and also bind to telomerase/PRC2/PEAT complexes. Additionally, we have detected novel interactions between TRB4-5 and EMF2 and VRN2, which are Su(z)12 subunits of PRC2.

The Gynandropsis gynandra genome provides insights into whole-genome duplications and the evolution of C4 photosynthesis in Cleomaceae

Gynandropsis gynandra (Cleomaceae) is a cosmopolitan leafy vegetable and medicinal plant, which has also been used as a model to study C4 photosynthesis due to its evolutionary proximity to C3 Arabidopsis (Arabidopsis thaliana). Here, we present the genome sequence of G. gynandra, anchored onto 17 main pseudomolecules with a total length of 740 Mb, an N50 of 42 Mb and 30,933 well-supported gene models. The G. gynandra genome and previously released genomes of C3 relatives in the Cleomaceae and Brassicaceae make an excellent model for studying the role of genome evolution in the transition from C3 to C4 photosynthesis. Our analyses revealed that G. gynandra and its C3 relative Tarenaya hassleriana shared a whole-genome duplication event (Gg-α), then an addition of a third genome (Th-α, +1×) took place in T. hassleriana but not in G. gynandra. Analysis of syntenic copy number of C4 photosynthesis-related gene families indicates that G. gynandra generally retained more duplicated copies of these genes than C3T. hassleriana, and also that the G. gynandra C4 genes might have been under positive selection pressure. Both whole-genome and single-gene duplication were found to contribute to the expansion of the aforementioned gene families in G. gynandra. Collectively, this study enhances our understanding of the polyploidy history, gene duplication and retention, as well as their impact on the evolution of C4 photosynthesis in Cleomaceae.

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.

An updated classification of the Brassicaceae (Cruciferae)

Based on recent achievements in phylogenetic studies of the Brassicaceae, a novel infrafamilial classification is proposed that includes major improvements at the subfamilial and supertribal levels. Herein, the family is subdivided into two subfamilies, Aethionemoideae (subfam. nov.) and Brassicoideae. The Brassicoideae, with 57 of the 58 tribes of Brassicaceae, are further partitioned into five supertribes, including the previously recognized Brassicodae and the newly established Arabodae, Camelinodae, Heliophilodae, and Hesperodae. Additional tribus-level contributions include descriptions of the newly recognized Arabidopsideae, Asperuginoideae, Hemilophieae, Schrenkielleae, and resurrection of the Chamireae and Subularieae. Further detailed comments on 17 tribes in need of clarifications are provided.

Synteny Identifies Reliable Orthologs for Phylogenomics and Comparative Genomics of the Brassicaceae

Large genomic data sets are becoming the new normal in phylogenetic research, but the identification of true orthologous genes and the exclusion of problematic paralogs is still challenging when applying commonly used sequencing methods such as target enrichment. Here, we compared conventional ortholog detection using OrthoFinder with ortholog detection through genomic synteny in a data set of 11 representative diploid Brassicaceae whole-genome sequences spanning the entire phylogenetic space. Then, we evaluated the resulting gene sets regarding gene number, functional annotation, and gene and species tree resolution. Finally, we used the syntenic gene sets for comparative genomics and ancestral genome analysis. The use of synteny resulted in considerably more orthologs and also allowed us to reliably identify paralogs. Surprisingly, we did not detect notable differences between species trees reconstructed from syntenic orthologs when compared with other gene sets, including the Angiosperms353 set and a Brassicaceae-specific target enrichment gene set. However, the synteny data set comprised a multitude of gene functions, strongly suggesting that this method of marker selection for phylogenomics is suitable for studies that value downstream gene function analysis, gene interaction, and network studies. Finally, we present the first ancestral genome reconstruction for the Core Brassicaceae which predating the Brassicaceae lineage diversification ∼25 million years ago.

Comparative Analysis of the Chloroplast Genome of Cardamine hupingshanensis and Phylogenetic Study of Cardamine

Cardamine hupingshanensis (K. M. Liu, L. B. Chen, H. F. Bai and L. H. Liu) is a perennial herbal species endemic to China with narrow distribution. It is known as an important plant for investigating the metabolism of selenium in plants because of its ability to accumulate selenium. However, the phylogenetic position of this particular species in Cardamine remains unclear. In this study, we reported the chloroplast genome (cp genome) for the species C. hupingshanensis and analyzed its position within Cardamine. The cp genome of C. hupingshanensis is 155,226 bp in length and exhibits a typical quadripartite structure: one large single copy region (LSC, 84,287 bp), one small single copy region (17,943 bp) and a pair of inverted repeat regions (IRs, 26,498 bp). Guanine-Cytosine (GC) content makes up 36.3% of the total content. The cp genome contains 111 unique genes, including 78 protein-coding genes, 29 tRNA genes and 4 rRNA genes. A total of 115 simple sequences repeats (SSRs) and 49 long repeats were identified in the genome. Comparative analyses among 17 Cardamine species identified the five most variable regions (trnH-GUG-psbA, ndhK-ndhC, trnW-CCA-trnP-UGG, rps11-rpl36 and rpl32-trnL-UAG), which could be used as molecular markers for the classification and phylogenetic analyses of various Cardamine species. Phylogenetic analyses based on 79 protein coding genes revealed that the species C. hupingshanensis is more closely related to the species C. circaeoides. This relationship is supported by their shared morphological characteristics.

A high-quality chromosome-level genome assembly of Pelteobagrus vachelli provides insights into its environmental adaptation and population history

Pelteobagrus vachelli is a freshwater fish with high economic value, but the lack of genome resources has severely restricted its industrial development and population conservation. Here, we constructed the first chromosome-level genome assembly of P. vachelli with a total length of approximately 662.13 Mb and a contig N50 was 14.02 Mb, and scaffolds covering 99.79% of the assembly were anchored to 26 chromosomes. Combining the comparative genome results and transcriptome data under environmental stress (high temperature, hypoxia and Edwardsiella. ictaluri infection), the MAPK signaling pathway, PI3K-Akt signaling pathway and apelin signaling pathway play an important role in environmental adaptation of P. vachelli, and these pathways were interconnected by the ErbB family and involved in cell proliferation, differentiation and apoptosis. Population evolution analysis showed that artificial interventions have affected wild populations of P. vachelli. This study provides a useful genomic information for the genetic breeding of P. vachelli, as well as references for further studies on fish biology and evolution.

Divergent evolutionary trajectories of bryophytes and tracheophytes from a complex common ancestor of land plants

The origin of plants and their colonization of land fundamentally transformed the terrestrial environment. Here we elucidate the basis of this formative episode in Earth history through patterns of lineage, gene and genome evolution. We use new fossil calibrations, a relative clade age calibration (informed by horizontal gene transfer) and new phylogenomic methods for mapping gene family origins. Distinct rooting strategies resolve tracheophytes (vascular plants) and bryophytes (non-vascular plants) as monophyletic sister groups that diverged during the Cambrian, 515-494 million years ago. The embryophyte stem is characterized by a burst of gene innovation, while bryophytes subsequently experienced an equally dramatic episode of reductive genome evolution in which they lost genes associated with the elaboration of vasculature and the stomatal complex. Overall, our analyses reveal that extant tracheophytes and bryophytes are both highly derived from a more complex ancestral land plant. Understanding the origin of land plants requires tracing character evolution across a diversity of modern lineages.

The evolutionary history of Cardamine bulbifera shows a successful rapid postglacial Eurasian range expansion in the absence of sexual reproduction

BACKGROUND AND AIMS

Sexual reproduction is known to drive plant diversification and adaptation. Here we investigate the evolutionary history and spatiotemporal origin of a dodecaploid (2n = 12x = 96) Eurasian deciduous woodland species, Cardamine bulbifera, which reproduces and spreads via vegetative bulb-like structures only. The species has been among the most successful range-expanding understorey woodland plants in Europe, which raises the question of the genetic architecture of its gene pool, since its hexaploid (2n = 6x = 48) but putatively outcrossing closest relative, C. quinquefolia, displays a smaller distribution range in Eastern Europe towards the Caucasus region. Cardamine bulbifera belongs to a small monophyletic clade of four species comprising also C. abchasica (2n = 2x = 16) and C. bipinnata (unknown ploidy) from the Caucasus region.

METHODS

We sequenced the genomes of the two polyploids and their two putative ancestors using Illumina short-read sequencing technology (×7-8 coverage). Covering the entire distribution range, genomic data were generated for 67 samples of the two polyploids (51 samples of C. bulbifera, 16 samples of C. quinquefolia) and 6 samples of the putative diploid taxa (4 samples of C. abchasica, 2 samples of C. bipinnata) to unravel the evolutionary origin of the polyploid taxa using phylogenetic reconstructions of biparentally and maternally inherited genetic sequence data. Ploidy levels of C. bulbifera and C. quinquefolia were analysed by comparative chromosome painting. We used genetic assignment analysis (STRUCTURE) and approximate Bayesian computation (ABC) modelling to test whether C. bulbifera represents genetically differentiated lineages and addressed the hypothesis of its hybrid origin. Comparative ecological modelling was applied to unravel possible niche differentiation among the two polyploid species.

KEY RESULTS

Cardamine bulbifera was shown to be a non-hybridogenous, auto-dodecaploid taxon of early Pleistocene origin, but with a history of past gene flow with its hexaploid sister species C. quinquefolia, likely during the last glacial maximum in shared refuge areas in Eastern Europe towards Western Turkey and the Crimean Peninsula region. The diploid Caucasian endemic C. abchasica is considered an ancestral species, which also provides evidence for the origin of the species complex in the Caucasus region. Cardamine bulbifera successfully expanded its distribution range postglacially towards Central and Western Europe accompanied by a transition to exclusively vegetative propagation.

CONCLUSIONS

A transition to vegetative propagation in C. bulbifera is hypothesized as the major innovation to rapidly expand its distribution range following postglacially progressing woodland vegetation throughout Europe. Preceding and introgressive gene flow from its sister species C. quinquefolia in the joint refuge area is documented. This transition and ecological differentiation may have been triggered by preceding introgressive gene flow from its sister species in the joint East European refuge areas.

The Brassicaceae genome resource (TBGR): A comprehensive genome platform for Brassicaceae plants

The Brassicaceae is an important plant family. We built a user-friendly, web-based, comparative, and functional genomic database, The Brassicaceae Genome Resource (TBGR, http://www.tbgr.org.cn), based on 82 released genomes from 27 Brassicaceae species. The TBGR database contains a large number of important functional genes, including 4,096 glucosinolate genes, 6,625 auxin genes, 13,805 flowering genes, 36,632 resistance genes, 1,939 anthocyanin genes, and 1,231 m6A genes. A total of 1,174,049 specific guide sequences for clustered regularly interspaced short palindromic repeats and 5,856,479 transposable elements were detected in Brassicaceae. TBGR also provides information on synteny, duplication, and orthologs for 27 Brassicaceae species. The TBGR database contains 1,183,851 gene annotations obtained using the TrEMBL, Swiss-Prot, Nr, GO, and Pfam databases. The BLAST, Synteny, Primer Design, Seq_fetch, and JBrowse tools are provided to help users perform comparative genomic analyses. All the genome assemblies, gene models, annotations, and bioinformatics results can be easily downloaded from the TBGR database. We plan to improve and continuously update the database with newly assembled genomes and comparative genomic studies. We expect the TBGR database to become a key resource for the study of the Brassicaceae.

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.

A Reappraisal of Polyploidy Events in Grasses (Poaceae) in a Rapidly Changing World

Around 80% of megaflora species became extinct at the Cretaceous–Paleogene (K–Pg) boundary. Subsequent polyploidy events drove the survival of thousands of plant species and played a significant historical role in the development of the most successful modern cereal crops. However, current and rapid global temperature change poses an urgent threat to food crops worldwide, including the world’s big three cereals: rice, wheat, and maize, which are members of the grass family, Poaceae. Some minor cereals from the same family (such as teff) have grown in popularity in recent years, but there are important knowledge gaps regarding the similarities and differences between major and minor crops, including how polyploidy affects their biological processes under natural and (a)biotic stress conditions and thus the potential to harness polyploidization attributes for improving crop climate resilience. This review focuses on the impact of polyploidy events on the Poaceae family, which includes the world’s most important food sources, and discusses the past, present, and future of polyploidy research for major and minor crops. The increasing accessibility to genomes of grasses and their wild progenitors together with new tools and interdisciplinary research on polyploidy can support crop improvement for global food security in the face of climate change.

The Phytotoxin Myrigalone A Triggers a Phased Detoxification Programme and Inhibits Lepidium sativum Seed Germination via Multiple Mechanisms including Interference with Auxin Homeostasis

Molecular responses of plants to natural phytotoxins comprise more general and compound-specific mechanisms. How phytotoxic chalcones and other flavonoids inhibit seedling growth was widely studied, but how they interfere with seed germination is largely unknown. The dihydrochalcone and putative allelochemical myrigalone A (MyA) inhibits seed germination and seedling growth. Transcriptome (RNAseq) and hormone analyses of Lepidium sativum seed responses to MyA were compared to other bioactive and inactive compounds. MyA treatment of imbibed seeds triggered the phased induction of a detoxification programme, altered gibberellin, cis-(+)-12-oxophytodienoic acid and jasmonate metabolism, and affected the expression of hormone transporter genes. The MyA-mediated inhibition involved interference with the antioxidant system, oxidative signalling, aquaporins and water uptake, but not uncoupling of oxidative phosphorylation or p-hydroxyphenylpyruvate dioxygenase expression/activity. MyA specifically affected the expression of auxin-related signalling genes, and various transporter genes, including for auxin transport (PIN7, ABCG37, ABCG4, WAT1). Responses to auxin-specific inhibitors further supported the conclusion that MyA interferes with auxin homeostasis during seed germination. Comparative analysis of MyA and other phytotoxins revealed differences in the specific regulatory mechanisms and auxin transporter genes targeted to interfere with auxin homestasis. We conclude that MyA exerts its phytotoxic activity by multiple auxin-dependent and independent molecular mechanisms.

Hemilophiacardiocarpa (Brassicaceae), a new species from Yunnan, southwest China

Hemilophiacardiocarpa (Brassicaceae), the sixth species of the Chinese endemic genus Hemilophia, is described and illustrated. This plant is found in the Jiaozishan Mountains in Dongchuan District, northern Yunnan, southwest China. Morphologically, it shows close relationships with H.rockii and H.pulchella, but differs from it in the leaf shape and size, inflorescence, flower size, shape of fruit and length of its pedicel. An updated key to the taxa of Hemilophia is provided.

Phenotypic plasticity guides Moricandia arvensis divergence and convergence across the Brassicaceae floral morphospace

Many flowers exhibit phenotypic plasticity. By inducing the production of several phenotypes, plasticity may favour the rapid exploration of different regions of the floral morphospace. We investigated how plasticity drives Moricandia arvensis, a species displaying within-individual floral polyphenism, across the floral morphospace of the entire Brassicaceae family. We compiled the multidimensional floral phenotype, the phylogenetic relationships, and the pollination niche of over 3000 species to construct a family-wide floral morphospace. We assessed the disparity between the two M. arvensis floral morphs (as the distance between the phenotypic spaces occupied by each morph) and compared it with the family-wide disparity. We measured floral divergence by comparing disparity with the most common ancestor, and estimated the convergence of each floral morph with other species belonging to the same pollination niches. Moricandia arvensis exhibits a plasticity-mediated floral disparity greater than that found between species, genera and tribes. The novel phenotype of M. arvensis moves outside the region occupied by its ancestors and relatives, crosses into a new region where it encounters a different pollination niche, and converges with distant Brassicaceae lineages. Our study suggests that phenotypic plasticity favours floral divergence and rapid appearance of convergent flowers, a process which facilitates the evolution of generalist pollination systems.

Pulvinatusia (Brassicaceae), a new cushion genus from China and its systematic position

The new genus and species Pulvinatusiaxuegulaensis (Brassicaceae) are described and illustrated. The species is a cushion plant collected from Xuegu La, Xizang, China. Its vegetative parts are most similar to those of Arenariabryophylla (Caryophyllaceae) co-occurring in the same region, while its leaves and fruits closely resemble those of Xerodrabapatagonica (Brassicaceae) from Patagonian Argentina and Chile. Family-level phylogenetic analyses based on both nuclear ITS and plastome revealed that it is a member of the tribe Crucihimalayeae, but the infra-/intergeneric relationships within the tribe are yet to be resolved.

Autopolyploidization affects transcript patterns and gene targeting frequencies in Physcomitrella

In Physcomitrella, whole-genome duplications affected the expression of about 3.7% of the protein-encoding genes, some of them relevant for DNA repair, resulting in a massively reduced gene-targeting frequency. Qualitative changes in gene expression after an autopolyploidization event, a pure duplication of the whole genome (WGD), might be relevant for a different regulation of molecular mechanisms between angiosperms growing in a life cycle with a dominant diploid sporophytic stage and the haploid-dominant mosses. Whereas angiosperms repair DNA double-strand breaks (DSB) preferentially via non-homologous end joining (NHEJ), in the moss Physcomitrella homologous recombination (HR) is the main DNA-DSB repair pathway. HR facilitates the precise integration of foreign DNA into the genome via gene targeting (GT). Here, we studied the influence of ploidy on gene expression patterns and GT efficiency in Physcomitrella using haploid plants and autodiploid plants, generated via an artificial WGD. Single cells (protoplasts) were transfected with a GT construct and material from different time-points after transfection was analysed by microarrays and SuperSAGE sequencing. In the SuperSAGE data, we detected 3.7% of the Physcomitrella genes as differentially expressed in response to the WGD event. Among the differentially expressed genes involved in DNA-DSB repair was an upregulated gene encoding the X-ray repair cross-complementing protein 4 (XRCC4), a key player in NHEJ. Analysing the GT efficiency, we observed that autodiploid plants were significantly GT suppressed (p < 0.001) attaining only one third of the expected GT rates. Hence, an alteration of global transcript patterns, including genes related to DNA repair, in autodiploid Physcomitrella plants correlated with a drastic suppression of HR.

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.

Evolutionary Analysis of the YABBY Gene Family in Brassicaceae

The YABBY gene family is one of the plant transcription factors present in all seed plants. The family members were extensively studied in various plants and shown to play important roles in plant growth and development, such as the polarity establishment in lateral organs, the formation and development of leaves and flowers, and the response to internal plant hormone and external environmental stress signals. In this study, a total of 364 YABBY genes were identified from 37 Brassicaceae genomes, of which 15 were incomplete due to sequence gaps, and nine were imperfect (missing C2C2 zinc-finger or YABBY domain) due to sequence mutations. Phylogenetic analyses resolved these YABBY genes into six compact clades except for a YAB3-like gene identified in Aethionema arabicum. Seventeen Brassicaceae species each contained a complete set of six basic YABBY genes (i.e., 1 FIL, 1 YAB2, 1 YAB3, 1 YAB5, 1 INO and 1 CRC), while 20 others each contained a variable number of YABBY genes (5-25) caused mainly by whole-genome duplication/triplication followed by gene losses, and occasionally by tandem duplications. The fate of duplicate YABBY genes changed considerably according to plant species, as well as to YABBY gene type. These YABBY genes were shown to be syntenically conserved across most of the Brassicaceae species, but their functions might be considerably diverged between species, as well as between paralogous copies, as demonstrated by the promoter and expression analysis of YABBY genes in two Brassica species (B. rapa and B. oleracea). Our study provides valuable insights for understanding the evolutionary story of YABBY genes in Brassicaceae and for further functional characterization of each YABBY gene across the Brassicaceae species.

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.

Genome evolution of the psammophyte Pugionium for desert adaptation and further speciation

Plants’ adaptations to and divergence in arid deserts have long fascinated scientists and the general public. Here, we present a genomic analysis of two congeneric desert plant species that clarifies their evolutionary history and shows that their common ancestor arose from a hybrid polyploidization, which provided genomic foundations for their survival in deserts. The whole-genome duplication was followed by translocation-based rearrangements of the ancestral chromosomes. Rapid evolution of genes in these reshuffled chromosomes contributed greatly to the divergences of the two species in desert microhabitats during which gene flow was continuous. Our results provide insights into plant adaptation in the arid deserts and highlight the significance of polyploidy-driven chromosomal structural variations in species divergence.

Deserts exert strong selection pressures on plants, but the underlying genomic drivers of ecological adaptation and subsequent speciation remain largely unknown. Here, we generated de novo genome assemblies and conducted population genomic analyses of the psammophytic genus Pugionium (Brassicaceae). Our results indicated that this bispecific genus had undergone an allopolyploid event, and the two parental genomes were derived from two ancestral lineages with different chromosome numbers and structures. The postpolyploid expansion of gene families related to abiotic stress responses and lignin biosynthesis facilitated environmental adaptations of the genus to desert habitats. Population genomic analyses of both species further revealed their recent divergence with continuous gene flow, and the most divergent regions were found to be centered on three highly structurally reshuffled chromosomes. Genes under selection in these regions, which were mainly located in one of the two subgenomes, contributed greatly to the interspecific divergence in microhabitat adaptation.

Pleistocene dynamics of the Eurasian steppe as a driving force of evolution: Phylogenetic history of the genus Capsella (Brassicaceae)

Capsella is a model plant genus of the Brassicaceae closely related to Arabidopsis. To disentangle its biogeographical history and intrageneric phylogenetic relationships, 282 individuals of all five currently recognized Capsella species were genotyped using a restriction digest-based next-generation sequencing method. Our analysis retrieved two main lineages within Capsella that split c. one million years ago, with western C. grandiflora and C. rubella forming a sister lineage to the eastern lineage consisting of C. orientalis. The split was attributed to continuous latitudinal displacements of the Eurasian steppe belt to the south during Early Pleistocene glacial cycles. During the interglacial cycles of the Late Pleistocene, hybridization of the two lineages took place in the southwestern East European Plain, leading to the allotetraploid C. bursa-pastoris. Extant genetic variation within C. orientalis postdated any extensive glacial events. Ecological niche modeling showed that suitable habitat for C. orientalis existed during the Last Glacial Maximum around the north coast of the Black Sea and in southern Kazakhstan. Such a scenario is also supported by population genomic data that uncovered the highest genetic diversity in the south Kazakhstan cluster, suggesting that C. orientalis originated in continental Asia and migrated north- and possibly eastwards after the last ice age. Post-glacial hybridization events between C. bursa-pastoris and C. grandiflora/rubella in the southwestern East European Plain and the Mediterranean gave rise to C. thracica. Introgression of C. grandiflora/rubella into C. bursa-pastoris resulted in a new Mediterranean cluster within the already existing Eurasian C. bursa-pastoris cluster. This study shows that the continuous displacement and disruption of the Eurasian steppe belt during the Pleistocene was the driving force in the evolution of Capsella.

Arabis alpina: A perennial model plant for ecological genomics and life-history evolution

Many model organisms were chosen and achieved prominence because of an advantageous combination of their life‐history characteristics, genetic properties and also practical considerations. Discoveries made in Arabidopsis thaliana, the most renowned noncrop plant model species, have markedly stimulated studies in other species with different biology. Within the family Brassicaceae, the arctic–alpine Arabis alpina has become a model complementary to Arabidopsis thaliana to study the evolution of life‐history traits, such as perenniality, and ecological genomics in harsh environments. In this review, we provide an overview of the properties that facilitated the rapid emergence of A. alpina as a plant model. We summarize the evolutionary history of A. alpina, including genomic aspects, the diversification of its mating system and demographic properties, and we discuss recent progress in the molecular dissection of developmental traits that are related to its perennial life history and environmental adaptation. From this published knowledge, we derive open questions that might inspire future research in A. alpina, other Brassicaceae species or more distantly related plant families.

Gene duplications and phylogenomic conflict underlie major pulses of phenotypic evolution in gymnosperms

Inferring the intrinsic and extrinsic drivers of species diversification and phenotypic disparity across the tree of life is a major challenge in evolutionary biology. In green plants, polyploidy (or whole-genome duplication, WGD) is known to play a major role in microevolution and speciation, but the extent to which WGD has shaped macroevolutionary patterns of diversification and phenotypic innovation across plant phylogeny remains an open question. Here, we examine the relationship of various facets of genomic evolution-including gene and genome duplication, genome size, and chromosome number-with macroevolutionary patterns of phenotypic innovation, species diversification, and climatic occupancy in gymnosperms. We show that genomic changes, such as WGD and genome-size shifts, underlie the origins of most major extant gymnosperm clades, and notably, our results support an ancestral WGD in the gymnosperm lineage. Spikes of gene duplication typically coincide with major spikes of phenotypic innovation, while increased rates of phenotypic evolution are typically found at nodes with high gene-tree conflict, representing historic population-level dynamics during speciation. Most shifts in gymnosperm diversification since the rise of angiosperms are decoupled from putative WGDs and instead are associated with increased rates of climatic occupancy evolution, particularly in cooler and/or more arid climatic conditions, suggesting that ecological opportunity, especially in the later Cenozoic, and environmental heterogeneity have driven a resurgence of gymnosperm diversification. Our study provides critical insight on the processes underlying diversification and phenotypic evolution in gymnosperms, with important broader implications for the major drivers of both micro- and macroevolution in plants.

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.

Comparison of glucosinolate diversity in the crucifer tribe Cardamineae and the remaining order Brassicales highlights repetitive evolutionary loss and gain of biosynthetic steps

We review glucosinolate (GSL) diversity and analyze phylogeny in the crucifer tribe Cardamineae as well as selected species from Brassicaceae (tribe Brassiceae) and Resedaceae. Some GSLs occur widely, while there is a scattered distribution of many less common GSLs, tentatively sorted into three classes: ancient, intermediate and more recently evolved. The number of conclusively identified GSLs in the tribe (53 GSLs) constitute 60% of all GSLs known with certainty from any plant (89 GSLs) and apparently unique GSLs in the tribe constitute 10 of those GSLs conclusively identified (19%). Intraspecific, qualitative GSL polymorphism is known from at least four species in the tribe. The most ancient GSL biosynthesis in Brassicales probably involved biosynthesis from Phe, Val, Leu, Ile and possibly Trp, and hydroxylation at the β-position. From a broad comparison of families in Brassicales and tribes in Brassicaceae, we estimate that a common ancestor of the tribe Cardamineae and the family Brassicaceae exhibited GSL biosynthesis from Phe, Val, Ile, Leu, possibly Tyr, Trp and homoPhe (ancient GSLs), as well as homologs of Met and possibly homoIle (intermediate age GSLs). From the comparison of phylogeny and GSL diversity, we also suggest that hydroxylation and subsequent methylation of indole GSLs and usual modifications of Met-derived GSLs (formation of sulfinyls, sulfonyls and alkenyls) occur due to conserved biochemical mechanisms and was present in a common ancestor of the family. Apparent loss of homologs of Met as biosynthetic precursors was deduced in the entire genus Barbarea and was frequent in Cardamine (e.g. C. pratensis, C. diphylla, C. concatenata, possibly C. amara). The loss was often associated with appearance of significant levels of unique or rare GSLs as well as recapitulation of ancient types of GSLs. Biosynthetic traits interpreted as de novo evolution included hydroxylation at rare positions, acylation at the thioglucose and use of dihomoIle and possibly homoIle as biosynthetic precursors. Biochemical aspects of the deduced evolution are discussed and testable hypotheses proposed. Biosyntheses from Val, Leu, Ile, Phe, Trp, homoPhe and homologs of Met are increasingly well understood, while GSL biosynthesis from mono- and dihomoIle is poorly understood. Overall, interpretation of known diversity suggests that evolution of GSL biosynthesis often seems to recapitulate ancient biosynthesis. In contrast, unprecedented GSL biosynthetic innovation seems to be rare.

Charting the genomic landscape of seed-free plants

During the past few years several high-quality genomes has been published from Charophyte algae, bryophytes, lycophytes and ferns. These genomes have not only elucidated the origin and evolution of early land plants, but have also provided important insights into the biology of the seed-free lineages. However, critical gaps across the phylogeny remain and many new questions have been raised through comparing seed-free and seed plant genomes. Here, we review the reference genomes available and identify those that are missing in the seed-free lineages. We compare patterns of various levels of genome and epigenomic organization found in seed-free plants to those of seed plants. Some genomic features appear to be fundamentally different. For instance, hornworts, Selaginella and most liverworts are devoid of whole-genome duplication, in stark contrast to other land plants. In addition, the distribution of genes and repeats appear to be less structured in seed-free genomes than in other plants, and the levels of gene body methylation appear to be much lower. Finally, we highlight the currently available (or needed) model systems, which are crucial to further our understanding about how changes in genes translate into evolutionary novelties.

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.

Genome-wide macroevolutionary signatures of key innovations in butterflies colonizing new host plants

The mega-diversity of herbivorous insects is attributed to their co-evolutionary associations with plants. Despite abundant studies on insect-plant interactions, we do not know whether host-plant shifts have impacted both genomic adaptation and species diversification over geological times. We show that the antagonistic insect-plant interaction between swallowtail butterflies and the highly toxic birthworts began 55 million years ago in Beringia, followed by several major ancient host-plant shifts. This evolutionary framework provides a valuable opportunity for repeated tests of genomic signatures of macroevolutionary changes and estimation of diversification rates across their phylogeny. We find that host-plant shifts in butterflies are associated with both genome-wide adaptive molecular evolution (more genes under positive selection) and repeated bursts of speciation rates, contributing to an increase in global diversification through time. Our study links ecological changes, genome-wide adaptations and macroevolutionary consequences, lending support to the importance of ecological interactions as evolutionary drivers over long time periods.